Biological nanopores having tunable pore diameters and uses thereof as analytical tools

Abstract

The invention relates to the field of nanopores, in particular to engineered Fragaceatoxin C (FraC) nanopores and their application in analyzing biopolymers and other (biological) compounds, such as single-molecule (protein) sequencing. Provided is a system comprising oligomeric FraC nanopores comprised in a lipid bilayer, wherein the sum of the nanopore fraction in the heptameric (Type II) state and the nanopore fraction in the hexameric (Type III) state represents at least 60% of the total number of FraC nanopores.

Description

CROSS-REFERENCE

This application is a National Stage Entry of International Patent Application No. PCT/NL2019/050588, filed Sep. 11, 2019, which claims the benefit of European Application No. 18193722.8, filed Sep. 11, 2018, each of which is entirely incorporated herein by reference.

SEQUENCE LISTING

The instant application contains a sequence listing which has been submitted electronically in ASCII format and is hereby incorporated by reference in its entirety. Said ASCII copy, created on Feb. 19, 2021, is named Sequence-Listing-as-filed-19feb.2021-20885-4501.txt, and is 9,181 bytes in size.

BACKGROUND

The invention relates generally to the field of nanopores and the use thereof in analyzing biopolymers and other (biological) compounds. In particular, it relates to engineered Fragaceatoxin C (FraC) nanopores and their application in single molecule analysis, such as single molecule peptide sequencing.

DETAILED DESCRIPTION

Biological nanopores are proteins that open nanoscale water conduits on biological or synthetic membranes. Under an external potential, the ionic current across single nanopore is used to recognize analytes traversing the nanopore. Most notably, nanopores are now used to sequence nucleic acids at the single-molecule level. In nanopore DNA sequencing, individual nucleic acid trains are threaded base-by-base through nanopores, while the ionic current is used to identify individual nucleobases^1,2.

Watanabe et al. (Analytical Chemistry 2017 89 (21), 11269-11277) describe the analysis of pore formation and detection of a single protein molecule using a large nanopore among five different pore-forming proteins, including FraC. It is demonstrated that the identification of appropriate pores for nanopore sensing can be achieved by classifying the channel current signals and performing noise analysis.

However, the sequencing of proteins with nanopores presents a new set of challenges. Amino acids have a larger chemically variability compared to nucleobases, and they cannot be uniformly captured or stretched by the electrical potential inside the nanopore. Furthermore, enzymes that process proteins or polypeptides amino acid-by-amino acid are not yet known.

Alternatively, proteins might be first fragmented and then the mass of individual peptides identified by nanopore currents. This approach would be similar to conventional protein sequencing approaches using tandem mass spectrometry. A nanopore peptide mass identifier, however, would have the advantage of being low-cost and portable and single-molecule. The latter is important because it would allow the analysis of the chemical heterogeneity in post-translational modifications and, especially when coupled to high-throughput devices, permit the detection of low-abundance proteins. Previous work with PEG molecules^3,4,5,6,7,8, neutral peptides⁹ or oligosaccharides¹⁰, uniformly charged peptides^{11,12,13,14,15} and other peptides¹⁶ revealed that there might be a direct correlation between the depth of the current blockade and the molecular weight of polymers, when the composition of the analyte is uniform. On the other hand, a wealth of other studies, including work with DNA^17,18 and amino acid enantiomers¹⁹ revealed that the chemical identity of molecules and especially the charge inside the nanopore²⁰ have a strong and unpredictable effect on the ionic current, suggesting that the identification of the mass of complex biopolymers such as peptides might not be possible. An additional complication is that peptides of opposite charge are not efficiently captured and analysed at a fixed potential^{16,21,22,23,24}. Finally, the diameter and geometry of biological nanopores cannot be easily adapted to study the array of sizes, shapes and chemical composition of polypeptides in solution.

Recently we have shown that octameric Fragaceatoxin C (FraC) nanopores²⁵ from the sea anemone Actinia fragacea can be used to study DNA²⁶, proteins and peptides²⁷. See also WO2018/012963 in the name of the present applicant. The transmembrane region of FraC is unique compared to other nanopores used in biopolymer analysis as it is formed by α-helices that describe a sharp and narrow constriction at the trans exit of the nanopore. Crucially, we showed that an electro-osmotic flow across the nanopore can be engineered to capture polypeptides at a fixed potential despite their charge composition²⁷. However, peptides smaller than 1.6 kDa in size translocated too fast across the nanopore to be sampled.

Based on these studies, the present inventors realized and recognized that nanopores with a smaller diameter are required to detect peptides with lower molecular weight. Therefore, they aimed at providing a strategy that allows for tuning the diameter of FraC nanopores, such that a larger range of peptides sizes can be identified.

It was surprisingly found that the FraC nanopore can be engineered to induce the formation of different nanopore types (herein referred to Type II and/or Type III) when comprised in the context of a lipid bilayer, thereby creating a biological nanopore with sub-nm constriction. Importantly, these novel, narrow types of nanopores allow for distinguishing (small) peptides differing by the substitution of one amino acid with a ˜40 Da resolution, while previous nanopore studies only reported differences of about 200 Da. Moreover, at selected pH conditions the FraC nanopore signal directly correlated to the mass of the peptide. The invention herewith provides a new and unique approach for the single-molecule identification of proteins based on nanopores.

In one embodiment, the invention relates to a system comprising oligomeric Fragaceatoxin C (FraC) nanopores comprised in a lipid bilayer, wherein the sum of the nanopore fraction in the Type II state and the nanopore fraction in the Type III state represents at least 60% of the total number of FraC nanopores.

For example, the sum of the Type II and Type III state nanopores represents at least 65%, preferably at least 70%, of the total number of FraC nanopores.

As used herein, the term “Type II” state refers to nanopores having an apparent heptameric stoichiometry, and/or a conductance of about 1.22-1.26 nS when assayed at pH 7.5 in a 1M NaCl solution or about 0.99-1.08 nS when assayed at pH 4.5 in a 1 M KCl solution. Conductance values are suitably determined by collecting single channels under −50 mV applied potential using 1 M NaCl, 15 mM Tris pH 7.5, or 1 M KCl, 0.1 M citric acid, 180 mM Tris base pH 4.5.

Type II FraC nanopores are furthermore characterized by an apparent pore size (at the narrowest constriction) of about 1.1 nm as calculated from homology modeling.

As used herein, the term “Type III” state refers to nanopores having an apparent hexameric stoichiometry, and/or a conductance of about 0.37-0.43 nS when assayed at pH 4.5 in a 1M KCl solution.

Type III FraC nanopores are furthermore characterized by a pore size (at the narrowest constriction) of about 0.8 nm as shown by homology modeling.

Accordingly, in one embodiment, the invention provides a system comprising oligomeric FraC nanopores comprised in a lipid bilayer, wherein the sum of the nanopore fraction in the heptameric (Type II) state and the nanopore fraction in the hexameric (Type III) state represents at least 60% of the total number of FraC nanopores.

In another embodiment, the invention provides a system comprising oligomeric FraC nanopores comprised in a lipid bilayer, wherein the sum of (i) the nanopore fraction showing a conductance of about 0.99-1.08 nS (Type II) when assayed at pH 4.5 in a 1 M KCl solution and (ii) the nanopore fraction showing a conductance of about 0.37-0.43 nS (Type III) when assayed at pH 4.5 in a 1 M KCl solution represents at least 60% of the total number of FraC nanopores.

Still further, the invention provides a system comprising oligomeric FraC nanopores comprised in a lipid bilayer, wherein the sum of the nanopore fraction having an apparent pore size of about 1.1 nm (Type II) and the nanopore fraction having an apparent pore size of about 0.8 nM (Type III) represents at least 60% of the total number of FraC nanopores.

The relative amounts of Type II and Type III nanopores can vary according to needs. In one aspect, at least 40%, preferably at least 50%, of the FraC nanopores is in the Type II state. Alternatively, or additionally, at least 15%, preferably at least 20%, of the FraC nanopores is in the Type III state.

In one embodiment, at least 60%, preferably at least 70%, of the FraC nanopores is in the Type II state. In another embodiment, at least 60%, preferably at least 70%, of the FraC nanopores is in the Type III state.

Also encompassed are systems comprising essentially one oligomeric form/Type of FraC. For example, in one embodiment, at least 90%, preferably at least 95%, of the FraC nanopores is present in the Type II state. In a specific aspect, all of the FraC nanopores are in the Type II state. In another embodiment, at least 90%, preferably at least 95%, of the FraC nanopores is present in the Type III state. In a specific aspect, all of the FraC nanopores are in the Type III state. The different oligomeric forms of FraC can be readily isolated using liquid chromatographic techniques, including size-exclusion, affinity, reverse-phase or ion exchange chromatography.

In a specific aspect, the FraC nanopores comprise or consist of mutant FraC monomers comprising one or more mutations that weaken the interaction between the nanopore and the lipid bilayer, i.e. the lipid interface.

Very good results are obtained when FraC is mutated at position W112 and/or W116. For example, in one embodiment, FraC is mutated at position W112, preferably while W116 is not mutated, or at position W116, preferably while W112 is not mutated. In a further embodiment, the FraC mutant comprises a mutation at both positions W112 and W116. According to the present invention, the W residues are substituted with either S, T, A, N, Q or G, preferably with S or T. FraC contains 179 amino acids with relative molecular weight of 20 kDa. The cDNA for FraC is available under the accession number FM958450 in GenBank. The polypeptide sequence of FraC is available under the accession number B9W5G6 in UniProt (SEQ ID NO: 1).

Moreover, the crystal structure of FraC was resolved in 2010 and deposited in the RCSB PDB under the accession number 3LIM. In a mutant according to the invention, the residue numbering corresponds to the residue numbering as in the FraC sequence according to UniProtKB accession number B9W5G6 (SEQ ID NO: 1).

The importance of (conserved) tryptophan residues for the functioning of pore-forming toxins has been previously studied. Tanaka et al.²⁵ revealed structures of FraC corresponding to four different stages of its activation route, namely the water-soluble form, the lipid-bound form, an assembly intermediate and the transmembrane pore. Mutational analysis revealed that mutant W112R/W116F lacks the ability to bind to lipid membrane, thus becoming completely inactive. Garcia-Linares et al., 2016 studied the role of the tryptophan residues in the specific interaction of the sea anemone Stichodactyla helianthus's Actinoporin Sticholysin II (StnI) with biological membranes.

It was found that residues W110 and W114 (corresponding to W112 and W116 of FraC) sustain the hydrophobic effect, which is one of the major driving forces for membrane binding in the presence of cholesterol. Notably, while the authors state that “the results obtained support actinoporins' Trp residues playing a major role in membrane recognition and binding”, they also conclude that “their residues have an only minor influence on the diffusion and oligomerization steps needed to assemble a functional pore”. Herewith, the findings of the present inventors that W112 and/or W116 are suitably engineered to tune the oligomeric state of FraC could not have reasonably been predicted in view of the prior art.

Hence, in one aspect the invention provides a system comprising oligomeric FraC nanopores comprised in a lipid bilayer, wherein the FraC nanopores comprise mutant FraC monomers comprising a mutation at position W112 and/or W116, preferably wherein the W residue(s) is/are independently substituted with either S, T, A, N, Q or G, preferably with S or T. In one embodiment, it comprises or consists of FraC mutant W112S or W112T. In another embodiment, it comprises or consists of mutant W116S or W116T. In a still further embodiment, the system comprises or consists of FraC mutant W112S/W116S, W112T/W116S or W112S/W116T.

The inventors noticed that type II FraC nanopores inserted in the lipid bilayer more efficiently at low pH. Therefore, to increase the production of type II nanopores at physiological pH, aspartic acid at position 109, which is located at the lipid interface, was exchanged for an uncharged residue. Satisfactorily, the fraction of type II nanopores at pH 7.5 increased from 23.0±4.9% to 48±3.6%, and a small fraction of type III nanopores appeared.

The invention also provides a system comprising oligomeric FraC nanopores comprised in a lipid bilayer, wherein the FraC nanopores comprise mutant FraC monomers comprising a mutation at position D109, preferably herein said mutation comprises the substitution of D with an uncharged residue, such as S or T, more preferably with S. The concomitant substitution of tryptophan at position 116 with serine showed a further small increased in the fraction of type I and type II nanopores at pH 7.5.

Therefore, the invention further relates to a system comprising oligomeric FraC nanopores comprised in a lipid bilayer, wherein the FraC nanopores comprise mutant FraC monomers comprising mutation D109S and one or both of W112S and W116S.

A system with the nanopores of the invention can accommodate peptides ranging from 22 to 4 amino acids in length. Even smaller peptides can be detected using further fine-tuning of the transmembrane region of the nanopore, for example by introducing amino acids with bulky side-chains. We also showed that the nanopores can discriminate differences between an alanine and glutamate (˜40 Da) in mixture of peptides. Furthermore, the inventors found that at exactly pH 3.8 the ionic signal of the peptides depended on the mass the analyte, while at higher pH values the current signal of negatively charged peptides was higher than expected from their mass alone.

Without wishing to be bound by any theory, the inventors' explanation is that the peptides analyzed lost their charge, while the constriction of the nanopore still retained enough negative charge to recognize the peptide charge. Most likely, a negatively charged constriction is important for creating an electrophoretic environment for peptide-mass recognition. At the same time, the electrostatic interaction of the constriction with negatively charged analytes might prevent the correct position of the analyte within the reading frame of the nanopore.

Presumably, peptides need to be uniformly charged which can be achieved by lowering the pH of the solution. At the same time, however, the constriction of the nanopore should be negatively charged in order to obtain optimal mass recognition. Obtaining both effects may be challenging, because by lowering the pH the charge of the constriction also becomes less charged.

Therefore, in addition to the amino acid substitution(s) disclosed herein above, mutant FraC monomers may comprise one or more unnatural amino acids comprising a moiety that holds a negative charge at low pH, preferably wherein said moiety is a sulfate or phosphate group. In one embodiment, a reside (e.g. at position 10) is mutated to cysteine and then oxidized, hence introducing a sulfonic (or sulfinic) group at that position. The charge of such group remains negative over the all pH range. With this approach, the recognition of peptides could be improved. Alternatively, peptides might be chemically modified (e.g. by esterification) to neutralize the negative charge.

In a system according to the invention, the FraC nanopores are comprised in a lipid bilayer. The reconstitution of FraC nanopores in lipid bilayers has been described in the art. Typically, the lipid bilayer comprises phosphatidylcholine (PC), preferably 1,2-diphytanoyl-sn-glycero-3-phosphocholine (DPhPC), optionally in combination with sphingomyelin (SM). Very good results are obtained when DPhPC and SM are present in about equal amounts by mass.

When a system according to the invention is in use, the nanopore is typically positioned between a first liquid medium and a second liquid medium, wherein at least one liquid medium comprises an analyte, and wherein the system is operative to detect a property of the analyte. In one embodiment, the system is operative to detect a property of the analyte comprises subjecting the nanopore to an electric field such that the analyte electrophoretically and/or electroosmotically translocates through the nanopore. As exemplified herein below, a system provided herein is particularly suitable for the analysis of a proteinaceous substance, preferably a peptide, more preferably a peptide up to about 30 amino acids in length. More in particular, a system of the invention provides for capture of peptides with different charge, recognition of the mass of the peptide and a resolution up to only 40 Da.

However, this is in no way to be understood that the invention is limited to applications relating to peptide analysis. For example, other analytes that can be detected using a system of the invention include (non-proteinaceous) biomarkers, antibiotics or other drugs, DNA, metabolites and small biological molecules.

The invention further relates to a mutant Fragaceatoxin C (FraC) polypeptide comprising one or more of the above mutations. These polypeptides are advantageously used in an (analytical) system herein disclosed. Also provided is an isolated nucleic acid molecule encoding a mutant FraC polypeptide according to the invention, and an expression vector comprising the isolated nucleic acid molecule. Still further, the invention provides a host cell comprising said expression vector.

In one embodiment, the mutant FraC polypeptide comprises a mutation at position D109, W112 and/or W116. For example, it comprises a mutation at W112 (optionally while W116 is not mutated) or it comprises a mutation at W116 (optionally while W112 is not mutated). In one embodiment, it comprises a mutation at both W112 and W116. As indicated herein above, the mutation(s) may comprise the substitution of D or W with S, T, A, N, Q or G, preferably with S or T. In a specific aspect, the invention provides mutant FraC W112S, FraC W116S, or FraC W112S/W116S, or its encoding nucleic acid molecule, or vector comprising the same. Still further, it provides a polypeptide comprising mutation D109S, preferably wherein the mutant is D109S/W116S, or its encoding nucleic acid molecule, or vector comprising the same.

Any one of these mutations may be supplemented with one or more unnatural amino acids comprising a moiety that holds a negative charge at low pH, for example wherein said moiety is a sulfate, sulfonic acid or phosphate group. Preferred positions for introducing such negative charge residue(s) include one or more of positions 10, 17 and 24.

In a specific aspect, the mutation W116 is supplemented with mutation D10C. The thiol group of the cysteine is then oxidized to sulfonic acid e.g. by incubation of FraC double mutant monomers with 10% (v/v) hydrogen peroxide. It was found that the introduction of a sulfonic acid moiety at position 10 of FraC gives rise to oligomerised pores that show a quiet signal in electrophysiology recordings as compared to a more noisy signal observed for nanopores that had not been subjected to oxidation. Accordingly, in one embodiment the invention provides a mutant comprising the D10C substitution, preferably in combination with one or more of W112S, W116S and D109S, more preferably in combination with at least W116S.

A further embodiment relates to a method for providing a system according to the invention, comprising the steps of

• • providing recombinant FraC monomers;
  • contacting said monomers with liposomes to assemble them into oligomers;
  • recovering the oligomers from the liposomes; and
  • contacting the oligomers with a lipid bilayer, which may contain sphingomyelin, to allow the formation of FraC nanopores.

In one embodiment, the contacting with a lipid bilayer is performed at a pH below 4.5, preferably below 4.0.

A method of the invention may furthermore comprise the step of isolating a fraction comprising FraC nanopores in the Type II state, and/or a fraction comprising FraC nanopores in the Type III state. In one aspect, it comprises isolating different oligomeric forms of FraC using a liquid chromatographic technique, including size-exclusion, affinity, reverse-phase or ion exchange chromatography.

A peptide mass-detecting FraC nanopore system of the present invention is advantageously integrated in real-time protein sequencing system. To that end, the system preferably comprises one or more further modifications.

In one embodiment, a protease-unfoldase pair is attached directly above (i.e. on the cis side of) the FraC nanopore. Then, cleaved peptides will be sequentially recognized and translocated across the nanopore. For example, the barrel-shaped ATP-dependent ClpXP protease is an ideal candidate because it can encase the digested peptides preventing its release in solution. Another approach is based on a protein complex that constitutes the proteasome, or any other protease. For example, the complex includes the 20S alpha/beta subunits of the proteasome, and the 19S regulatory particle, of which the ATPase is the minimal required unit. However, other proteases could also be used. The protease will cleave the polypeptide specifically (for example it could cut after a positively charged residue or a negatively charged residue or a aromatic residue etcetera), or it will be engineered to cut specifically, or it will cut at nonspecific locations within the polypeptide chain. The protease will ideally encase the substrate and will allow the docking of other components (e.g. unfoldases) to feed the polypeptide to the unfoldase active site.

The attachment may be of a covalent or non-covalent nature. For example, it can achieved by chemical attachment, by genetic fusion, or by introducing a binding loop into the FraC nanopore that can interact non-covalently with the peptidase.

We demonstrated that the peptides entering the cis side of the nanopore have a high probability of exiting the nanopore to the trans chamber, which will prevent duplicate detection events. Furthermore, we showed that at low pH peptides are likely to be captured and their mass recognized by the nanopore at a fixed applied potential irrespectively of their chemical composition. If such low pH values will not be compatible with enzymatic activity, asymmetric solutions on both side of the nanopore can be used^34,35. In such system, conditions in the cis side can be tuned to optimize the ATPase activity of the unfoldase-peptidase, while the pH and ionic strength of the trans side can be optimized to capture and recognize individual peptides.

A FraC nanopore or a mutant FraC polypeptide as provided herein is advantageously used in all sorts of analyte analysis, including peptide or DNA analysis, preferably wherein peptide analysis comprises peptide mass detection and/or peptide sequencing. However, whereas the advantageous properties of a system provided herein are demonstrated in the context of peptide analysis, a person skilled in the art will appreciate that it can be used for various applications. Other possible applications of the invention include the following:

• • Peptide post-modification detection (glycosylation etc), proteomics;
  • DNA sequencing with higher accuracy, or DNA post-modification detection (methylation etc);
  • Other small analytes detection or polymers analysis with higher resolution.
  • Single molecule protein and DNA sequencing;
  • Directly peptides, biomarkers, antibiotics and small molecules detection in human samples;
  • Trapping of different size proteins for binding ligands analysis like glucose etc.

Mass spectrometry is the workhorse of the proteomics field. At present, the nanopore system falls short from the resolution of commercial mass spectrometers. A peptide mass-analyzer device based on FraC nanopores as herein disclosed has distinctive advantages compared to conventional mass spectrometers, which are expensive, extremely complex and unwieldy, and are not single-molecule. By contrast, nanopores can be integrated in portable and low-cost devices containing hundreds of thousands of individual sensors. Hence, in one embodiment the FraC nanopore system is integrated in a portable device comprising a plurality of individual FraC nanopore systems as herein described.

In addition, the electrical nature of the signal allows sampling biological samples in real-time. Furthermore, since the nanopore reads individual molecules, the signal contains additional information not available for ensemble techniques. In particular, single-molecule detection, especially when coupled to high throughput analysis, will allow detecting low abundance peptides and to unravel the chemical heterogeneity in post translational modifications, challenges that are hard to address with conventional mass spectrometry. For example, the invention also provides for the use of a FraC nanopore system or a mutant FraC polypeptide in single molecule detection, preferably in combination with high throughput analysis.

LEGEND TO THE FIGURES

FIGS. 1A-1G. Preparation and characterization of type I, type II and type III FraC nanopores. FIG. 1A, Cut through of a surface representation of WT-FraC oligomer (PDB: 4TSY)²⁵ colored according to the vacuum electrostatic potential as calculated by Pymol. One protomer is shown as a carton presentation with tryptophans 112 and 116 displayed as spheres. FIG. 1B, Percentage of the distribution of type I, type II and type III for WT-FraC, W112S-FraC, W116S-FraC and W112S-W116S-FraC at pH 7.5 and 4.5. FIG. 1C, IV curves of type II nanopores formed by WT-FraC, W116S-FraC and W112S-W116S-FraC at pH 7.5 (15 mM Tris·HCl, 1 M KCl). FIG. 1D, Single nanopore conductance of W116S-FraC in 1 M KCl (0.1 M citric acid and 180 mM Tris base) at pH 4.5. FIG. 1E, Typical current traces for the three nanopore types of W116S-FraC in 1 M KCl at pH 4.5 under −50 mV applied potential. FIG. 1F, Reversal potentials measured under asymmetric condition of KCl (1960 mM cis, 467 mM trans) at pH 4.5 for the three W116S-FraC nanopore types. The ion selectivity was given by using the Goldman-Hodgkin-Katz equation (equation 1). FIG. 1G, Molecular models of the three type FraC nanopores constructed from the FraC crystals structure using the symmetrical docking function of Rosetta. The electrophysiology recordings were performed with 10 kHz sampling and 2 kHz filter. The error bars and color shadow in the I-V curves are standard deviations from three repeats at least.

FIGS. 2A-2F. Single channel conductance distributions of FraC nanopores at pH 7.5 and 4.5. FIG. 2A, The table reports the average conductance values which were obtained by fitting Gaussian functions to conductance histograms. S.D. represents the standard deviation of all single channels (number given as n). FIGS. 2B-2F, Each panel represents a different batch of FraC nanopores as indicated. Single channels were collected under −50 mV applied potential. pH 7.5 and 4.5 were obtained using 1 M NaCl, 15 mM Tris, or 1 M KCl, 0.1 M citric acid, 180 mM Tris base respectively.

FIGS. 3A-3B. Discrimination of angiotensin peptides in mixture with type II W116S-FraC nanopores. FIG. 3A, (i) Sequences of angiotensin I, II, III and IV with corresponding Ires % measured at −30 mV. (ii) Typical blockades provoked by the four angiotensin peptides. (iii) Color density plot of the Ires % versus the standard deviation of the current amplitude for angiotensin I added to the cis compartment, and after the further addition of angiotensin II, angiotensin III, and angiotensin IV to the cis chamber (iv). FIG. 3B, Discrimination of angiotensin II and angiotensin A. (i) Table showing the sequences, the molecular weights and the Ires % of the peptides. The peptides differ by one amino acid as shown in red. (ii) Representative traces of the peptide blockades. Color density plot of the Ires % versus the standard deviation of the current amplitude for angiotensin II blockades prior (iii) and after (iv) the further addition of angiotensin A to the cis chamber. All measurements and recordings were performed in pH 4.5 buffer containing 1 M KCl, 0.1 M citric acid, 180 mM Tris base with a 50 kHz sampling and 10 kHz filter. Standard deviations were calculated from minimum three repeats. Color density plot was created with Origin. FIG. 3A discloses SEQ ID NOS 14, 17, 20, and 22, respectively, in order of appearance. FIG. 3B discloses SEQ ID NOS 17 and 19, respectively, in order of appearance.

FIGS. 4A-4C. Evaluation of biological peptides having different chemical compositions. Relation between the molecular weight and Ires % of peptide using: (FIG. 4A) type I WT-FraC nanopores, (FIG. 4B) type II W116S-FraC nanopores and (FIG. 4C) type III W112S-W116S-FraC nanopores at pH 4.5. The solid line represents a second order polynomial fitting. Current blockades were measured at −30 mV for type I and II pore, and at −50 mV for type III pore. Error bars are standard deviations obtained from at least three measurements.

FIGS. 5A-5D. A nanopore peptide mass spectrometer at pH 3.8. FIG. 5A, Amino acid sequences of four different peptides and their overall charge at different pH. The chargeable amino acids are underlined. FIG. 5B, pH dependence of the Ires % for the four peptides (cis) shown in a using type II W116S-FraC nanopores under −30 mV applied potential. FIG. 5C, Comparison of the Ires % versus the mass of peptides at pH 4.5 and 3.8. FIG. 5D, Voltage dependence of c-Myc dwell times at different pHs. All electrophysiology measurements were carried out in 1 M KCl, 0.1 M citric acid, and pH was adjusted with 1 M Tris base to desired values. 50 kHz sampling rate and 10 kHz filter was used for collecting the current events. Error bars are standard deviations obtained from at least three measurements. The charges of the peptides were calculated according to the pKa for individual amino acids³⁶. FIG. 5A discloses SEQ ID NOS 14-15, 17 and 20, respectively, in order of appearance.

FIGS. 6A-6C: Discrimination of short peptide mixture with type III FraC nanopores comprising mutant W112S-W116S-FraC. FIG. 6A, Sequence, Ires % (˜50 mV) and MW of angiotensin IV, angiotensin 4-8, endomorphin I and leucine enkephalin. FIG. 6B, Typical blockades provoked by the different peptides. FIG. 6C, Color density plot showing the Ires % versus the standard deviation of the current blockade for the mixture of angiotensin IV, angiotensin 4-8, endomorphin I and leucine-enkephalin. All measurements and recordings were performed in pH 4.5 buffer containing 1 M KCl, 0.1 M citric acid, 180 mM Tris base with a 50 kHz sampling and 10 kHz filter. Standard deviations were calculated from three repeats at least. FIG. 6A discloses SEQ ID NOS 22-24 and 26, respectively, in order of appearance.

FIGS. 7A-7B. Characterization of type II FraC nanopores comprising an oxidized cysteine at position 10. Difference between the D10C/W116S type II pore (FIG. 7A) and the oxidized D10C/W116S type II pore (FIG. 7B). Recordings were performed in a buffer containing 1 M NaCl, pH 7.5, ±50 mV.

EXPERIMENTAL SECTION

Materials and Methods

Chemicals

Endothelin 1 (≥97%, CAS #117399-94-7), endothelin 2 (≥97%, CAS #123562-20-9), dynorphin A porcine (≥95%, CAS #80448-90-4), angiotensin I (90%, CAS #70937-97-2), angiotensin II (≥93%, CAS #4474-91-3), c-Myc 410-419 (≥97%, #M2435), Asn1-Val5-Angiotensin II (≥97%, CAS #20071-00-5), Ile7 Angiotensin III (≥95%, #A0911), leucine enkephalin (≥95%, #L9133), 5-methionine enkephalin (≥95%, CAS #82362-17-2), endomorphin I (≥95%, CAS #189388-22-5), pentane (≥99%, CAS #109-66-0), hexadecane (99%, CAS #544-76-3), Trizma®hydrochloride (≥99%, CAS #1185-53-1), Trizma®base (≥99%, CAS #77-86-1), Potassium chloride (≥99%, CAS #7447-40-7), N,N-Dimethyldodecylamine N-oxide (LADO, ≥99%, CAS #1643-20-5) were obtained from Sigma-Aldrich. Pre angiotensin 1-14 (≥97%, #002-45), angiotensin 1-9 (≥95%, #002-02), angiotensin A (≥95%, #002-36), angiotensin III (≥95%, #002-31), angiotensin IV (≥95%, #002-28) were purchased from Pheonix Pharmaceuticals. Angiotensin 4-8 (≥95%) was synthesized by BIOMATIK. 1,2-diphytanoyl-sn-glycero-3-phosphocholine (DPhPC, #850356P) and sphingomyelin (Porcine brain, #860062) were purchased from Avanti Polar Lipids. Citric acid (99.6%, CAS #77-92-9) was obtained from ACROS. n-Dodecyl (3-D-maltoside (DDM, ≥99.5%, CAS #69227-93-6) was bought from Glycon Biochemical EmbH. DNA primers were synthesized from Integrated DNA Technologies (IDT), enzymes from Thermo scientific. All peptides were dissolved with Milli-Q water without further purification and stored in −20° C. freezer. pH 7.5 buffer containing 15 mM Tris in this study was prepared by dissolving 1.902 g Trizma® HCl and 0.354 g Trizma® base in 1 litre Milli-Q water (Millipore, Inc).

FraC Monomer Expression and Purification

FraC gene containing NcoI and HindIII restriction site at the 5′ and 3′ ends, respectively, and a sequence encoding for a poly histidine tag at the 3′ terminus was cloned to a pT7-SC1 plasmid. Plasmids were transformed into BL21(DE3) E.cloni® competent cell by electroporation. Cells were grown on LB agar plate containing 100 μg/mL ampicillin overnight at 37° C. The entire plate was then harvested and inoculated into 200 mL fresh 2YT media and the culture was grown with 220 rpm shaking at 37° C. until the optical density at 600 nm of the cell culture reached 0.8. Then 0.5 mM IPTG was added to the media and the culture was transferred to 25° C. for overnight growth with 220 rpm shaking. The next day the cells were centrifuged (2000×g, 30 minutes) and the pellet stored at −80° C. Cell pellets harvested from 100 mL culture media were used to purify FraC monomer. 30 mL of cell lysis buffer (150 mM NaCl, 15 mM Tris, 1 mM MgCl₂, 4 M urea, 0.2 mg/mL lysozyme and 0.05 unit/mL DNase) was added to resuspend the pellet and vigorously mixed for 1 hour. Cell lysate was then sonicated with Branson Sonifier 450 for 2 minutes (duty cycle 10%, output control 3). Then the crude lysate was centrifuged down at 4° C. for 30 minutes (5400×g), and the supernatant incubated with 100 μL Ni-NTA beads (Qiagen) for 1 hour with gentle shaking. Beads were spun down and loaded to a Micro Bio-spin column (Bio-rad). 10 mL of SDEX buffer (150 mM NaCl, 15 mM Tris, pH 7.5) containing 20 mM imidazole was used to wash the beads, and proteins were eluded with 150 μL elution buffer (SDEX buffer, 300 mM imidazole). The concentration of protein was measured by measuring the absorption at 280 nm with Nano-drop 2000 (Thermo scientific) using the elution buffer as blank. To further confirm the purity of monomer, monomeric FraC was diluted to 0.5 mg/mL using the elution buffer and 9 μL of the diluted sample was mixed with 3 μL of 4× loading buffer (250 mM Tris HCl, pH 6.8. 8% SDS, 0.01% bromophenol blue and 40% glycerol) and then loaded to 12% SDS-PAGE gel. Gels were run for 30 min with 35 mA constant applied current, and stained with coomassie dye (InstantBlue™, Expdedeon) for more than 1 hour before viewing using a gel imager (Gel Doc™, Bio-rad).

≥FraC|B9W5G6 (amino acid sequence)
SADVAGAVIDGAGLGFDVLKTVLEALGNVKRKIAVGIDNESGKTWTAMNTY
FRSGTSDIVLPHKVAHGKALLYNGQKNRGPVATGVVGVIAYSMSDGNTLAV
LFSVPYDYNWYSNWWNVRVYKGQKRADQRMYEELYYHRSPFRGDNGWHSRG
LGYGLKSRGFMNSSGHAILEIHVTKA

• • 6×His-WtFraC (amino acid sequence) as used in the present invention. Bold residues indicate residues of the N- and C-terminal end that were added to the original sequence.

MA              SADVAGAVIDGAGLGFDVLKTVLEALGNVKRKIAVGIDNESGKTWTAMN
TYFRSGTSDIVLPHKVAHGKALLYNGQKNRGPVATGVVGVIAYSMSDGNTL
AVLFSVPYDYNWYSNWWNVRVYKGQKRADQRMYEELYYHRSPFRGDNGWHS
RGLGYGLKSRGFMNSSGHAILEIHVTKAGSAHHHHHH
>6xHis-WtFraC (DNA sequence)
ATGGCGAGCGCCGATGTCGCGGGTGCGGTAATCGACGGTGCGGGTCTGGGC
TTTGACGTACTGAAAACCGTGCTGGAGGCCCTGGGCAACGTTAAACGCAAA
ATTGCGGTAGGGATTGATAACGAATCGGGCAAGACCTGGACAGCGATGAAT
ACCTATTTCCGTTCTGGTACGAGTGATATTGTGCTCCCACATAAGGTGGCG
CATGGTAAGGCGCTGCTGTATAACGGTCAAAAAAATCGCGGTCCTGTCGCG
ACCGGCGTAGTGGGTGTGATTGCCTATAGTATGTCTGATGGGAACACACTG
GCGGTACTGTTCTCCGTGCCGTACGATTATAATTGGTATAGCAATTGGTGG
AACGTGCGTGTCTACAAAGGCCAGAAGCGTGCCGATCAGCGCATGTACGAG
GAGCTGTACTATCATCGCTCGCCGTTTCGCGGCGACAACGGTTGGCATTCC
CGGGGCTTAGGTTATGGACTCAAAAGTCGCGGCTTTATGAATAGTTCGGGC
CACGCAATCCTGGAGATTCACGTTACCAAAGCAGGCTCTGCGCATCATCAC
CACCATCACTGATAAGCTT

FraC Mutation Preparation

FraC mutants were prepared according to MEGAWHOP method. 25 μL REDTaq® ReadyMix™ was mixed with 4 μM primer (see Table 1) containing the desired mutation with 50 ng plasmid (pT7-SC1 with wild type FraC gene) as template and the final volume was brought to 50 μL with MilliQ water.

TABLE 1
Primer sequences used in this study for
preparing FraC mutants.
Primer name   DNA sequences
T7 promoter   5′ TAATACGACTCACTATAGGG 3′
T7 terminator 5′ GCTAGTTATTGCTCAGCGG 3′
W112S Fw      5′ ACGATTATAATAGCTATAGCAATTGGTGG 3′
W116S Fw      5′ ATTGGTATAGCAATAGCTGGAACGTG 3′
W112116S Fw   5′ GTACGATTATAATAGCTATAGCAATAGCTGGA
              ACGTGC 3′
D109S ReV     5′ TGCTATACCAATTATAGCTGTACGGCA 3′

The PCR protocol was initiated by 150 seconds denature step at 95° C., followed by 30 cycles of denaturing (95° C., 15 s), annealing (55° C., 15 s), and extension (72° C., 60 s). The PCR products (MEGA primer) were combined and purified using a QIAquick PCR purification kit with final DNA concentration around 200 ng/μL. The second PCR was performed for whole plasmid amplification. 2 μL of MEGA primer, 1 μL Phire II enzyme, 10 μL 5× Phire buffer, 1 μL 10 mM dNTPs, were mixed with PCR water to 50 μL final volume. PCR started with pre-incubated at 98° C. (30 s) and then 25 cycles of denaturing (98° C., 5 s), annealing (72° C., 90 s), extension (72° C., 150 s). When the PCR was completed, 1 μL Dpn I enzyme was added and the mixture kept at 37° C. for 1 hour. Then the temperature was raised to 65° C. for 1 minute to inactivate the enzyme. Products were then transformed into E. cloni® 10G cells (Lucigen) competent cell by electroporation. Cells were plated on LB agar plates containing 100 μg/mL ampicillin and grew at 37° C. overnight. Single clones were enriched and sent for sequencing.

Sphingomyelin-DPhPC Liposome Preparation

20 mg sphingomyelin and 20 mg DPhPC (1,2-diphytanoyl-sn-glycero-3-phosphocholine) were dissolved in 4 mL pentane with 0.5% v/v ethanol and brought to a 50 mL round flask. The solvent was then evaporated by rotation and using a hair-dryer to warm-up the flask. After evaporation, the flask was kept at ambient temperature for an additional 30 minutes. The lipid film was resuspended with 4 mL SDEX buffer (150 mM NaCl, 15 mM Tris, pH 7.5) and the solution immersed in a sonication bath for 5 minutes. Liposome suspensions were stored at −20° C.

FraC Oligomerization

FraC oligomerization was triggered by incubation of FraC monomers with sphingomelyin-DPhPC liposomes. Frozen liposome were thawed and sonicated in a water bath for one minute. FraC monomers were diluted to one mg/mL using SDEX buffer, and then 50 μL of FraC monomers were added to 50 μl of a 10 mg/mL liposome solution to obtain a mass ratio of 10:1 (liposome:protein). The lipoprotein solution was incubated at 37° C. for 30 min to allow oligomerization. Then 10 μl of 5% (w/v, 0.5% final) LADO was added to the lipoprotein solution to solubilize the liposomes. After clarification (typically 1 minute) the solution was transferred to a 50 mL Falcon tube. Then 10 mL of SDEX buffer containing 0.02% DDM and 100 μL of pre-washed Ni-NTA beads were added to the Falcon tube and mixed gently in shaker for 1 hour at room temperature. The beads were then spun down and loaded to a Micro Bio-spin column. 10 mL wash buffer (150 mM NaCl, 15 mM Tris, 20 mM imidazole, 0.02% DDM, pH 7.5) was used to wash the beads and oligomers eluded with 100 μL elution buffer (typically 200 mM EDTA, 75 mM NaCl, 7.5 mM Tris pH 7.5, 0.02% DDM). The FraC oligomers were stored at 4° C. Under these conditions the nanopores are stable for several months.

W112S-W116S-FraC Oligomer Separation with His-Trap Chromatography

200 μL of W112S-W116S-FraC monomers (3 mg/mL) were incubated with 300 μL of Sphingomyelin-DPhPC liposome (10 mg/mL) and kept at 4° C. for 48 hours after which 0.5% LADO (final concentration) was added to solubilize the lipoprotein. Then the buffer was exchanged to the 500 mM NaCl, 15 mM Tris, 0.01% DDM, 30 mM imidazole, pH 7.5 (binding buffer) using a PD SpinTrap G-25 column. W112S-W116S-FraC oligomers were then loaded to Histrap HP 1 mL column (General Electric) using an AKTA pure FPLC system (General Electric). The loaded oligomers were washed with 10 column volumes of 500 mM NaCl, 15 mM Tris, 0.01% DDM, 30 mM imidazole, pH 7.5, prior applying an imidazole gradient (from 30 mM to 1 M imidazole, 500 mM NaCl, 15 mM Tris, 0.01% DDM, pH 7.5) over 30 column volumes. The signal was monitored with the absorbance at 280 nm and fractions were collected when the absorbance was higher than 5 mAu.

Electrophysiology Measurement and Data Analysis

Electrical recordings were performed as explained in details previously^27,37. I_O, referring to open pore current, were measured by fitting Gaussian functions to event amplitude histograms. Residual current values (Ires %) were calculated by dividing the blockade current (I_B) by open pore current (I_B/I_O×100%). Dwell times and inter-event times were measured by fitting single exponentials to histograms of cumulative distribution.

Ion Permeability Measurement

In order to measure reversal potentials, a single channel was obtained under symmetric conditions (840 mM KCl, 500 μL in each electrophysiology chamber) and the electrodes were balanced. The 400 μL of a buffered stock solution of 3.36 M KCl was slowly added to cis chamber, while 400 μL of salt free buffered solution was added to the trans chamber to obtain a total volume of 900 μL (trans:cis, 467 mM KCl:1960 mM KCl). After the equilibrium was reached, IV curves were collected from −30 to +30 mV. The resulting voltage at zero current is the reversal potential (Vr). The ion selectivity (P_K⁺/P_Cl⁻) was then calculated using the Goldman-Hodgkin-Katz equation (equation 1) where [a_K+/Cl−]_cis/trans is the activity of the K⁺ or Cl⁻ in the cis or trans compartment, R the gas constant, T the temperature and F the Faraday's constant.

$\begin{array}{cc}\frac{P_{K^{+}}}{P_{C{l}^{-}}}=\frac{{\left[a_{{\mathrm{Cl}}^{-}}\right]}_{\mathrm{trans}}-{\left[a_{{\mathrm{Cl}}^{-}}\right]}_{cis}{e}^{V_{r}F/RT}}{{\left[a_{K^{+}}\right]}_{\mathrm{trans}}{e}^{V_{r}F/RT}-{\left[a_{K^{+}}\right]}_{\mathrm{cis}}}& \left(1\right)\end{array}$

The activity of ions was calculated by multiplying the molar concentration of the ion for the mean ion mobility (0.649 for 500 mM KCl, and 0.573 for 2000 mM). Ag/AgCl electrodes were surrounded by 2.5% agarose bridge in 2.5 M NaCl.

Molecular Models of Type I, II and III FraC Nanopores

The 3D models with different multimeric order, ranging from five to nine monomers, were constructed with the symmetrical docking function of Rosetta³⁸. A monomer without lipids was extracted from the crystal structure of FraC with lipids (PDB_ID 4tsy). Symmetrical docking arranged this monomer around a central rotational axis ranging in order form 5 to 9. In total Rosetta generated and scored 10 000 copies for each symmetry. In all cases, a multimeric organization with a symmetry similar to the crystal structure could be identified as a top scoring solution. However, in the pentameric assembly the multimer interface was not fully satisfied as compared to the crystal structure, with large portions left exposed. The 9-fold symmetric model however exhibited a significant drop in Rosetta score compared to the 6- 7- and 8-fold symmetric models indicating an unfavored assembly of the nonameric assembly with the 6-7- and 8-fold assemblies as the most plausible. To create lipid bound models, the crystal structure with lipids was superimposed on each monomer of the generated models, allowing the lipid coordinates to be transferred. The residues within 4.5 angstrom of the lipids were minimized with the Amber10 force field.

Example 1: Engineering the Size of FraC Nanopores

One of the main challenges in biological nanopores analysis is to obtain nanopores with different size and shape. Most of biological nanopores are formed by multiple repeats of individual monomers. Hence, different nanopore sizes might be obtained by engineering the protein oligomeric composition²⁸. We noticed that at pH 7.5 a small fraction of Wild Type FraC (WT-FraC) nanopores showed a lower conductance (1.26±0.08 nS, −50 mV, type II WT-FraC, FIG. 1B) compared to the dominant fraction (2.26±0.08 nS, −50 mV, type I WT-FraC), suggesting that FraC might be able spontaneously to assemble into nanopores with smaller size. At pH 4.5 yet a smaller nanopore conductance was observed (0.42±0.03 nS, type III WT-FraC, −50 mV, FIG. 1B). We noticed that the reconstitution of lower conductance nanopores depended to several purification conditions. In particular, we observed that the occurrence of type II and type III nanopores increased when the oligomers were stored in solution for several weeks or when the concentration of monomeric WT-FraC was reduced during oligomerisation.

In an attempt to enrich for type II and type III FraC nanopores, the interaction between the nanopore and the lipid interface was weakened by substituting W112 and W116 at the lipid interface of FraC (FIG. 1A) with serine. The inventors reasoned that a lower concentration of monomers during liposome-triggered oligomerisation would increase the population of lower molecular mass oligomers. Surprisingly, it was found that at pH 4.5 using W116S-FraC and W112S-W116S-FraC oligomers, type II and type III FraC nanopores were the dominant species, respectively (FIG. 1B, FIG. 2). Conveniently, the different nanopore types could be separated by Ni-NTA affinity chromatography using an imidazole gradient. Furthermore, it was found that enrichment of type II and type III FraC nanopores could also be obtained at pH 7.5 by replacing aspartate 109 at the lipid interface with serine (see FIG. 2E, Table 2).

TABLE 2
relative amounts of Type I, Type II and Type III for each of the
FraC nanopores investigated at neutral and acidic pH.
		Type I(%)	SD	Type II(%)	SD	Type III(%)	SD
pH 7.5	Wild type	85.7	3.8	14.3	3.8	0.0	0.0
	W112S	61.9	4.3	38.1	4.3	0.0	0.0
	W116S	61.1	5.7	38.9	5.7	0.0	0.0
	W112116S	27.1	3.9	72.9	3.9	0.0	0.0
	D109S	50.3	3.8	48.0	3.6	1.7	1.5
	D109SW116S	29.3	9.5	66.7	8.4	4.0	4.0
pH 4.5	Wild type	42.5	10.6	51.9	7.7	5.6	7.9
	W116S	29.0	4.0	47.0	4.6	24.0	2.3
	W112SW116S	21.7	4.7	38.0	8.5	40.3	9.3
	D109S	35.7	2.1	56.3	9.1	8.0	7.0
	D109SW116S	19.3	8.4	64.3	6.7	16.3	3.2

Among FraC nanopores of the same type, the lipid interface modifications caused by W112S and W116S substitutions did not alter the conductance and ion selectivity as compared to that of wild type (FIG. 1C, FIG. 2, Table 3) suggesting that the modifications did not altered the overall fold of the nanopores. When characterised in lipid bilayers, type I, type II and type III nanopores showed a well-defined single conductance distribution and a steady open pore current (FIGS. 1D-1E).

Interestingly, Type I, Type II and Type III nanopores showed increasing cation selectivity (from 2.0 for type I to 4.2 for type III W116S-FraC nanopores at pH 4.5 (FIG. 1F, Table 3), most likely reflecting a larger overlap of the electrical double layer in the nanopores with a narrower constriction.

TABLE 3
Ion selectivity of different FraC pores at pH 7.5 and 4.5.
		pH 7.5	pH 4.5
		Reversal		Reversal
		potential (mV)	PK+/PCl−	potential (mV)	PK+/PCl−
WT-FraC	Type I	17.2 ± 1.2	3.6 ± 0.4	10.5 ± 1.4	2.1 ± 0.2
	Type II	20.8 ± 1.6	5.2 ± 0.9	12.3 ± 1.2	2.4 ± 0.2
	Type III	/	/	20.6 ± 1.1	5.0 ± 0.6
W116S-FraC	Type I	/	/	10.1 ± 0.9	2.0 ± 0.1
	Type II	/	/	12.8 ± 0.7	2.5 ± 0.2
	Type III	/	/	18.8 ± 0.5	4.2 ± 0.2
W112S-W116S-FraC	Type I	/	/	8.8 ± 1.2	1.9 ± 0.2
	Type II	/	/	14.0 ± 0.1	2.8 ± 0.1
	Type III	/	/	20.1 ± 0.6	4.8 ± 0.3

The ion selectivity (P_K+/P_Cl−) was calculated from the reversal potential according to the Goldman-Hodgkin-Katz equation:

$\frac{P_{K^{+}}}{P_{{\mathrm{Cl}}^{-}}}=\frac{{\left[a_{{\mathrm{Cl}}^{-}}\right]}_{\mathrm{trans}}-{\left[a_{{\mathrm{Cl}}^{-}}\right]}_{\mathrm{cis}}{e}^{V_{r}F/\mathrm{RT}}}{{\left[a_{K^{+}}\right]}_{\mathrm{trans}}{e^{V_{r}F/\mathrm{RT\_}}\left[a_{K^{+}}\right]}_{\mathrm{cis}}},$ where Vr is the reversal potential, P_K+/P_Cl− the ion selectivity, a the activity of ions and F the Farady constant. Electrophysiology recordings were carried out with 1960 mM KCl in the cis solution and 467 mM KCl in the trans solution. The activity of ions was calculated by multiplying the molar concentration of the ion for the mean ion mobility (0.649 for 500 mM KCl, and 0.573 for 2000 mM)³. Errors are given as standard deviations calculated from 3 experiments at least.

These findings strongly suggest that the three types of FraC nanopores represent nanopores with different protomeric compositions. Molecular modelling allowed predicting the diameter of type II (1.1 nm) and type III (0.8 nm) nanopores (FIG. 1G); and revealed that type III FraC is most likely the biological nanopore with the smallest constriction known to date.

Example 2: Identification of Peptides Containing Single Amino Acid Substitutions Using Type II or Type III FraC Nanopores as Sensor

Type II FraC nanopores were used to sample a series of angiotensin peptides (which in blood regulate blood pressure and fluid balance. The peptides were added to the cis side of type II W116S-FraC nanopores and the induced ionic current blockades (I_B) was measured. Residual currents percent (Ires %, defined as I_B/I_O×100) were used instead of current blockades because they provided more reliable values when comparing different nanopores. Results are shown in FIG. 3A and Table 4.

Angiotensin I (DRVYIHPFHL (SEQ ID NO: 14), 1296.5 Da), showed the deepest blockade (Ires %=8.8±0.2) and angiotensin IV (VYIHPF (SEQ ID NO: 22), 774.9 Da) the shallowest blockade (Ires %=38.9±4.0). The residual current of angiotensin II (DRVYIIPF (SEQ ID NO: 17), 1046.2 Da, Ires %=17.9±1.3) and angiotensin III (RVYIIPF (SEQ ID NO: 20), 931.1 Da, Ires %=22.1±0.5) fell at intermediate values. When the four peptides were tested simultaneously, individual peptides could be readily discriminated (FIG. 3A).

TABLE 4
Peptide analysis using different types of FraC nanopores at pH 4.5
		Molecular			Ires %
		weight	pH		(IB/IO) %	Dwell time
Peptide	Sequence	(g/mol)	7.5	pH 4.5	pH4.5	(ms)
WT-FraC type I pore, −30 mV
Endothelin 2	CSCSSWLDKECVYFCHLDIIW	2546.9	−2.15	0.36	6.1 ± 1.8	104.0 ± 29.9
Endothelin 1	CSCSSLMDKECVYFCHLDIIW	2491.9	−2.15	0.36	7.5 ± 0.5	19.73 ± 1.95
Dynorphin A	YGGFLRRIRPKLKWDNQ	2147.5	3.76	4.48	15.1 ± 2.6	3.68 ± 0.76
Pre angiotensinogen	DRVYIHPFHLVIHN	1758.9	0.03	3.45	24.6 ± 2.3	0.29 ± 0.04
Angiotensin I	DRVYIHPFHL	1296.5	−0.06	2.46	43.4 ± 0.9	0.15 ± 0.04
W116S-FraC type II pore, -30 mV
Angiotensin I	DRVYIHPFHL	1296.5	−0.06	2.46	8.8 ± 0.2	0.54 ± 0.01
c-Myc 410-419	EQKLISEEDL	1203.3	−3.24	−1.19	30.0 ± 3.4	0.12 ± 0.01
Angiotensin 1-9	DRVYIHPFH	1183.3	−0.06	2.46	14.0 ± 0.2	0.37 ± 0.04
Angiotensin II	DRVYIHPF	1046.2	−0.15	1.47	17.9 ± 1.3	0.37 ± 0.04
AsnlVal5
Angiotensin II	NRVYVHPF	1031.2	0.85	2.03	19.6 ± 0.2	0.34 ± 0.06
Angiotensin A	ARVYIHPF	1002.2	0.85	2.03	21.0 ± 0.6	0.34 ± 0.02
Angiotensin III	RVYIHPF	931.1	0.85	2.03	22.1 ± 0.5	0.35 ± 0.04
Ile7 Angiotensin III	RVYIHPI	897.1	0.85	2.03	24.3 ± 0.4	0.19 ± 0.05
Angiotensin IV	VYIHPF	774.9	−0.15	1.02	38.9 ± 4.0	0.15 ± 0.06
W112S-W116S-FraC type III pore, −50 mV
Angiotensin IV	VYIHPF	774.9	−0.15	1.02	1.1 ± 0.8	0.61 ± 0.07
Angiotensin 4-8	YIHPF	675.8	−0.15	1.02	8.2 ± 0.4	0.40 ± 0.04
Endomorphin I	YPWF	610.7	−0.24	0.04	19.2 ± 0.5	0.32 ± 0.04
Met5 Enkephalin	YGGFM	573.7	−0.24	0.04	33.5 ± 0.7	0.16 ± 0.02
Leucine Enkephalin	YGGFL	555.6	−0.24	0.04	34.5 ± 2.4	0.20 ± 0.05
The electrophysiology solution contained 1M KCl, 0.1M citric acid, 180 mM Tris base at pH 4.5.
Recordings were performed using a 50 kHz sampling and applying 10 kHz Bassel filter.
Standard deviations were obtained for at least three measurements.
The charges of the peptides were calculated according to the pKa for individual amino acid36.

The resolution limit of the nanopore sensor was challenged by sampling a mixture of peptides. Remarkably, angiotensin II and angiotensin A, having an identical composition except for the initial amino acid (aspartate in angiotensin II vs. alanine in angiotensin A), appeared as distinctive peaks in Ires % plots (FIG. 3B). Smaller differences in peptide mass, e.g. the 34 Da difference between phenylalanine and leucine in angiotensin III and Ile7 angiotensin III, were observed but not easily detected (data not shown), indicating the resolution of our system at ˜40 Da. Smaller peptides such as angiotensin II 4-8 (YIHPF (SEQ ID NO: 23), 675.8 Da), endomorphin I (YPWF (SEQ ID NO: 24), 610.7 Da) or leucine enkephalin (YGGFL (SEQ ID NO: 26), 555.6 Da) translocated too quickly across type II W116S-FraC nanopores to be sampled. However, they could be readily measured using type III W112S-W116S-FraC nanopores (Table 4;

FIGS. 6A-6C).

Example 3: A Nanopore Mass Spectrometer for Peptides

Although the ability of biological nanopores to distinguish between known analytes is useful, a more powerful application would be the identification of peptide masses directly from ionic current blockades without holding prior knowledge of the analyte identity. In an effort to assess FraC nanopores as peptide mass analyzer, additional peptides were tested at pH 4.5 and 1 M KCl using type I, type II and type III FraC nanopores (FIGS. 4A-4C, Table 4). Crucially, analytes with largely different charge compositions were included.

It was found that for most of peptides there was a direct correlation between the size and the residual current (FIGS. 4A-4C). A notable exception was c-Myc 410-419 (1203.3 Da), an intentionally selected peptide because it includes a long stretch of negatively charged residues (FIG. 5A). The overall negative charge of the peptide at pH 4.5 (see Table 4) was expected to have an effect on both peptide capture and recognition. c-Myc 410-419 could be readily captured at negative applied potentials (trans), indicating that the cis to trans electroosmotic flow across the nanopore can overcome the electrostatic energy barrier opposing peptide capture. However, the Ires % of c-Myc 410-419 (30.0±3.4) was significantly higher than the expected value (FIG. 4b).

We reasoned that such anomaly might be due to the interaction between the acidic amino acids of the peptide and the negatively charged constriction of FraC nanopores. Thus, we lowered the pH solution to values where the aspartate and glutamate side chains in the peptides are expected to be protonated, hence become neutral (FIG. 5A). Rewardingly, at pH 3.8, the signal corresponding to c-Myc 410-419 (1203.3 Da) fell between the signal of angiotensin I (1296.5 Da), and angiotensin II (1046.2 Da, FIGS. 5B-5C). This indicates that, after losing its negative charges, the peptide blockades scaled with the expected mass of the peptides.

It has been assumed^1,30,31 and experimentally proven³² that the voltage dependence of the average dwell time (off) can report on the translocation of a molecule across a nanopore. Under a negative bias (trans) for positively charged peptides (added in cis) both electrophoretic and electroosmotic forces (from cis to trans) promote the entry and translocation²⁷ across the nanopore. For negatively charged peptides, such as c-Myc 410-419 at pH 4.5 (FIG. 5A), the electroosmotic driving force must be stronger than the opposing electrophoretic force. The voltage dependence of τ_off was examined for c-Myc 410-419 at different pH values (FIG. 5D). At pH 4.5 the peptide exhibited a maximum in τ_off at −50 mV, suggesting that at low potentials c-Myc 410-419 returns to the cis chamber (<50 mV), and at higher potentials (≥50 mV) c-Myc 410-419 exits to the trans chamber. At pH 3.8 and lower, we observed a decrease in τ_off at higher potentials, indicating that c-Myc 410-419 crosses the membrane to the trans chamber.

As shown in FIGS. 6A-6C, type III FraC nanopore can detect differences in peptide length down to 4 amino acids (mass around 500 Dalton) in a peptide mixture. It was also found that the residual current signal correlated well with the mass of peptides, suggesting that Type III can be used as a detector for a peptide having a mass down to −500 Da.

Example 4: Peptide Mass Identifier at pH 3

This example shows that mutation D10C can be used as additional mutation to obtain a FraC pore showing a quiet signal in electrophysiology recordings.

Using mutant W116S as exemplary mutant, the aspartic acid at position 10 of FraC was converted to cysteine by site-directed mutagenesis. The thiol group of cysteine was then oxidized to sulfonic acid by incubation of FraC monomers with 10% hydrogen peroxide (v/v), which was dissolved in regular buffer (e.g. 10 mM Tris buffer pH 7.5, 150 mM NaCl). As a control, the double mutant was left without oxidation.

D10C/W116S FraC was oligomerized, and the oligomers tested in electrical recordings. FIGS. 7A-7B shows the trace comparison between the D10C/W116S pore and oxidized D10C/W116S pore to demonstrate the difference after oxidization. Oligomerised pores from oxidized D10C/W116S FraC monomers showed a quiet signal in electrophysiology recordings, as compared to a more noisy signal observed for nanopores that had not been subjected to oxidation.

REFERENCES

• 1. Clarke, J. et al. Continuous base identification for single-molecule nanopore DNA sequencing. Nat. Nanotechnol. 4, 265-270 (2009).
• 2. Derrington, I. M. et al. Nanopore DNA sequencing with MspA. Proc. Natl. Acad. Sci. 107, 16060-16065 (2010).
• 3. Bezrukov, S. M., Vodyanoy, I., Brutyan, R. A. & Kasianowicz, J. J. Dynamics and free energy of polymers partitioning into a nanoscale pore. Macromolecules 29, 8517-8522 (1996).
• 4. Robertson, J. W. F. et al. Single-molecule mass spectrometry in solution using a solitary nanopore. Proc. Natl. Acad. Sci. 104, 8207-8211 (2007).
• 5. Baaken, G. et al. High-Resolution Size-Discrimination of Single Nonionic Synthetic Polymers with a Highly Charged Biological Nanopore. ACS Nano 9, 6443-6449 (2015).
• 6. Aksoyoglu, M. A. et al. Size-dependent forced PEG partitioning into channels: VDAC, OmpC, and α-hemolysin. Proc. Natl. Acad. Sci. 113, 9003-9008 (2016).
• 7. Krasilnikov, O. V., Rodrigues, C. G. & Bezrukov, S. M. Single polymer molecules in a protein nanopore in the limit of a strong polymer-pore attraction. Phys. Rev. Lett. 97, 1-4 (2006).
• 8. Oukhaled, A. G., Biance, A. L., Pelta, J., Auvray, L. & Bacri, L. Transport of long neutral polymers in the semidilute regime through a protein nanopore. Phys. Rev. Lett. 108, 1-4 (2012).
• 9. Zhao, Q., Jayawardhana, D. A., Wang, D. & Guan, X. Study of peptide transport through engineered protein channels. J. Phys. Chem. B 113, 3572-3578 (2009).
• 10. Bacri, L. et al. Discrimination of neutral oligosaccharides through a nanopore. Biochem. Biophys. Res. Commun. 412, 561-564 (2011).
• 11. Movileanu, L., Schmittschmitt, J. P., Scholtz, J. M. & Bayley, H. Interactions of peptides with a protein pore. Biophys. J. 89, 1030-1045 (2005).
• 12. Mohammad, M. M., Prakash, S., Matouschek, A. & Movileanu, L. Controlling a single protein in a nanopore through electrostatic traps. J. Am. Chem. Soc. 130, 4081-4088 (2008).
• 13. Piguet, F. et al. Identification of single amino acid differences in uniformly charged homopolymeric peptides with aerolysin nanopore. Nat. Commun. 9, (2018).
• 14. Lamichhane, U. et al. Peptide translocation through the mesoscopic channel: Binding kinetics at the single molecule level. Eur. Biophys. J. 42, 363-369 (2013).
• 15. Stefureac, R., Long, Y. T., Kraatz, H. B., Howard, P. & Lee, J. S. Transport of α-helical peptides through α-hemolysin and aerolysin pores. Biochemistry 45, 9172-9179 (2006).
• 16. Chavis, A. E. et al. Single Molecule Nanopore Spectrometry for Peptide Detection. ACS Sensors 2, 1319-1328 (2017).
• 17. Maglia, G., Restrepo, M. R., Mikhailova, E. & Bayley, H. Enhanced translocation of single DNA molecules through-hemolysin nanopores by manipulation of internal charge. Proc. Natl. Acad. Sci. 105, 19720-19725 (2008).
• 18. Stoddart, D., Heron, A. J., Mikhailova, E., Maglia, G. & Bayley, H. Single-nucleotide discrimination in immobilized DNA oligonucleotides with a biological nanopore. Proc. Natl. Acad. Sci. 106, 7702-7707 (2009).
• 19. Boersma, A. J. & Bayley, H. Continuous stochastic detection of amino acid enantiomers with a protein nanopore. Angew. Chemie—Int. Ed. 51, 9606-9609 (2012).
• 20. Stoddart, D. et al. Nucleobase recognition in ssDNA at the central constriction of the hemolysin pore. Nano Lett. 10, 3633-3637 (2010).
• 21. Kennedy, E., Dong, Z., Tennant, C. & Timp, G. Reading the primary structure of a protein with 0.07 nm 3 resolution using a subnanometre-diameter pore. Nat. Nanotechnol. 11, 968-976 (2016).
• 22. Li, S., Cao, C., Yang, J. & Long, Y. Detection of Peptides with Different Charges and Lengths by Using the Aerolysin Nanopore. 4, 1-5 (2018).
• 23. Asandei, A. et al. Electroosmotic Trap Against the Electrophoretic Force Near a Protein Nanopore Reveals Peptide Dynamics during Capture and Translocation. ACS Appl. Mater. Interfaces 8, 13166-13179 (2016).
• 24. Chinappi, M. & Cecconi, F. Protein sequencing via nanopore based devices: a nanofluidics perspective. J. Phys. Condens. Matter in press, (2018).
• 25. Tanaka, K., Caaveiro, J. M. M., Morante, K., Gonzilez-Manãs, J. M. & Tsumoto, K. Structural basis for self-assembly of a cytolytic pore lined by protein and lipid. Nat. Commun. 6, 4-6 (2015).
• 26. Wloka, C., Mutter, N. L., Soskine, M. & Maglia, G. Alpha-Helical Fragaceatoxin C Nanopore Engineered for Double-Stranded and Single-Stranded Nucleic Acid Analysis. Angew. Chemie—Int. Ed. 55, 12494-12498 (2016).
• 27. Huang, G., Willems, K., Soskine, M., Wloka, C. & Maglia, G. Electro-osmotic capture and ionic discrimination of peptide and protein biomarkers with FraC nanopores. Nat. Commun. 8, 1-13 (2017).
• 28. Soskine, M., Biesemans, A., De Maeyer, M. & Maglia, G. Tuning the size and properties of ClyA nanopores assisted by directed evolution. J. Am. Chem. Soc. 135, 13456-13463 (2013).
• 29. Aqvist, J. et al. Dipoles Localized at Helix Termini of Proteins Stabilize Charges. Proc. Natl. Acad. Sci. 88, 2026-2030 (1991).
• 30. Rincon-Restrepo, M., Mikhailova, E., Bayley, H. & Maglia, G. Controlled translocation of individual DNA molecules through protein nanopores with engineered molecular brakes. Nano Lett. 11, 746-750 (2011).
• 31. Wanunu, M., Sutin, J., McNally, B., Chow, A. & Meller, A. DNA translocation governed by interactions with solid-state nanopores. Biophys. J. 95, 4716-4725 (2008).
• 32. Biesemans, A., Soskine, M. & Maglia, G. A Protein Rotaxane Controls the Translocation of Proteins Across a ClyA Nanopore. Nano Lett. 15, 6076-6081 (2015).
• 33. Ho, C. W. et al. Engineering a nanopore with co-chaperonin function. Sci. Adv. 1, 1-9 (2015).
• 34. Wanunu, M., Morrison, W., Rabin, Y., Grosberg, A. Y. & Meller, A. Electrostatic focusing of unlabelled DNA into nanoscale pores using a salt gradient. Nat. Nanotechnol. 5, 160-165 (2010).
• 35. Stoddart, D., Franceschini, L., Heron, A., Bayley, H. & Maglia, G. DNA stretching and optimization of nucleobase recognition in enzymatic nanopore sequencing. Nanotechnology 26, 10-16 (2015).
• 36. Stryer, L. Biochemistry. Biochemistry (4th ed.) (1995).
• 37. Soskine, M., Biesemans, A. & Maglia, G. Single-molecule analyte recognition with ClyA nanopores equipped with internal protein adaptors. J. Am. Chem. Soc. 137, 5793-5797 (2015).
• 38. Andre, I., Bradley, P., Wang, C. & Baker, D. Prediction of the structure of symmetrical protein assemblies. Proc. Natl. Acad. Sci. 104, 17656-17661 (2007).

Claims

1. A system comprising oligomeric Fragaceatoxin C (FraC) nanopores disposed in a lipid bilayer, wherein a sum of a first number of nanopores in a heptameric (Type II) state and a second number of nanopores in a hexameric (Type III) state is at least 60% of a total number of the oligomeric FraC nanopores.

2. The system of claim 1, wherein the sum of the first number of nanopores in the heptameric (Type II) state and the second number of nanopores in the hexameric (Type III) state is at least 65% of the total number of the oligomeric FraC nanopores.

3. The system of claim 1, wherein the first number of nanopores in the heptameric (Type II) state is at least 60% of the total number of the oligomeric FraC.

4. The system of claim 1, wherein the second number of nanopores in the hexameric (Type III) state is at least 60% of the total number of the oligomeric FraC nanopores.

5. The system of claim 1, wherein the oligomeric FraC nanopores comprise mutant FraC monomers comprising a mutation at position W112 or W116.

6. The system of claim 5, wherein said mutation comprises a substitution of a tryptophan (W) with a serine (S), a threonine (T), an alanine (A), an asparagine (N), a glutamine (Q) or a glycine (G).

7. The system of claim 5, wherein said mutant FraC monomers comprise mutation W112S or W116S.

8. The system of claim 7, wherein the mutant FraC monomers comprise mutations (a) W112S and W116S, (b) D109S and W116S, or (c) D10C and W116S.

9. The system of claim 1, wherein the oligomeric FraC nanopores comprise mutant FraC monomers comprising a mutation at position D109.

10. The system of claim 9, wherein the mutant FraC monomers comprise mutation D109S or D109T.

11. The system of claim 1, wherein the system is operative to detect a property of an analyte by subjecting the oligomeric FraC nanopores to an electric field such that the analyte electrophoretically or electro-osmotically translocates through a nanopore of the oligomeric FraC nanopores.

12. The system of claim 11, wherein the analyte is a peptide.

13. A method for providing the system of claim 1, comprising: (i) providing recombinant FraC monomers;(ii) contacting said recombinant FraC monomers with liposomes to assemble the recombinant FraC monomers into oligomers;(iii) recovering the oligomers from the liposomes; and(iv) contacting the oligomers with the lipid bilayer to allow formation of the oligomeric FraC nanopores.

14. The method of claim 13, further comprising isolating the first number of nanopores in the heptameric (Type II) state and the second number of nanopores in the hexameric (Type III) state.

15. A method of performing peptide analysis, the method comprising (a) providing the system of claim 1 and a solution comprising a peptide for analysis, wherein the solution is disposed on one side of the lipid bilayer; and(b) subjecting a nanopore of the oligomeric FraC nanopores to an electric field such that the peptide electrophoretically or electro-osmotically translocates through the nanopore.

16. The method of claim 15, wherein the system is integrated in a portable device comprising a plurality of individual systems.

17. A method of performing single molecule detection, the method comprising: (a) providing the system of claim 1 and a solution comprising a plurality of single molecules, wherein the solution is disposed on one side of the lipid bilayer;(b) subjecting a nanopore of the oligomeric FraC nanopores to an electric field such that a single molecule of the plurality of single molecules electrophoretically or electro-osmotically translocates through the nanopore; and(c) detecting the single molecule translocated through the nanopore.

18. The system of claim 1, wherein a Type II FraC nanopore in the first number of nanopores has (i) a conductance from about 50% to about 68% of a conductance of a Type I FraC nanopore when assayed at pH 7.5 in a 1M NaCl solution, or (ii) a conductance from about 41% to about 62% of the conductance of the Type I FraC nanopore when assayed at pH 4.5 in a 1M KCl solution.

19. The system of claim 18, wherein a Type III FraC nanopore in the second number of nanopores has a conductance from about 16% to about 24% of the conductance of the Type I FraC nanopore when assayed at pH 4.5 in a 1M KCl solution.

20. The system of claim 1, wherein a Type II FraC nanopore in the first number of nanopores or a Type III FraC nanopore in the second number of nanopores comprises a monomer comprising an amino acid sequence as set forth in SEQ ID NO: 1.

21. The system of claim 1, wherein a Type II FraC nanopore in the first number of nanopores or a Type III FraC nanopore in the second number of nanopores comprises a monomer comprising an engineered variation of an amino acid sequence as set forth in SEQ ID NO: 1.

22. The system of claim 21, wherein the engineered variation comprises a mutant sequence comprising a mutation at positions 109, 112, or 116 when aligned to the amino acid sequence as set forth in SEQ ID NO: 1.

23. The system of claim 22, wherein the mutation comprises a substitution of aspartic acid (D) with an uncharged amino acid residue.

24. The system of claim 1, wherein the system further comprises a first electrode on a first side of the lipid bilayer and a second electrode on a second side of the lipid bilayer.