The emerging landscape of single-molecule protein sequencing technologies.

Single-cell profiling methods have had a profound impact on the understanding of cellular heterogeneity. While genomes and transcriptomes can be explored at the single-cell level, single-cell profiling of proteomes is not yet established. Here we describe new single-molecule protein sequencing and identification technologies alongside innovations in mass spectrometry that will eventually enable broad sequence coverage in single-cell profiling. These technologies will in turn facilitate biological discovery and open new avenues for ultrasensitive disease diagnostics.

Role of conservative mutations in protein multi-property adaptation.

Protein physicochemical properties must undergo complex changes during evolution, as a response to modifications in the organism environment, the result of the proteins taking up new roles or because of the need to cope with the evolution of molecular interacting partners. Recent work has emphasized the role of stability and stability-function trade-offs in these protein adaptation processes. In the present study, on the other hand, we report that combinations of a few conservative, high-frequency-of-fixation mutations in the thioredoxin molecule lead to largely independent changes in both stability and the diversity of catalytic mechanisms, as revealed by single-molecule atomic force spectroscopy. Furthermore, the changes found are evolutionarily significant, as they combine typically hyperthermophilic stability enhancements with modulations in function that span the ranges defined by the quite different catalytic patterns of thioredoxins from bacterial and eukaryotic origin. These results suggest that evolutionary protein adaptation may use, in some cases at least, the potential of conservative mutations to originate a multiplicity of evolutionarily allowed mutational paths leading to a variety of protein modulation patterns. In addition the results support the feasibility of using evolutionary information to achieve protein multi-feature optimization, an important biotechnological goal.

Single-aminoacid discrimination in proteins with homogeneous nanopore sensors and neural networks.

A technology capable of sequencing individual protein molecules would revolutionize our understanding of biological processes. Nanopore technology can analyze single heteropolymer molecules such as DNA by measuring the ionic current flowing through a single nanometer hole made in an electrically insulating membrane. This current is sensitive to the monomer sequence. However, proteins are remarkably complex and identifying a single residue change in a protein remains a challenge. In this work, I show that simple neural networks can be trained to recognize protein mutants. Although these networks are quickly and efficiently trained, their ability to generalize in an independent experiment is poor. Using a thermal annealing protocol on the nanopore sample, and examining many mutants with the same nanopore sensor are measures aimed at reducing training data variability which produce an increase in the generalizability of the trained neural network. Using this approach, we obtain a 100% correct assignment among 9 mutants in >50% of the experiments. Interestingly, the neural network performance, compared to a random guess, improves as more mutants are included in the dataset for discrimination. Engineered nanopores prepared with high homogeneity coupled with state-of-the-art analysis of the ionic current signals may enable single-molecule protein sequencing.

Oligonucleotide-Directed Protein Threading Through a Rigid Nanopore.

Nanopore technology enables the detection and analysis of single protein molecules. The technique measures the ionic current passing through a single pore inserted in an electrically insulating membrane. The translocation of the protein molecule through the pore causes a modulation of the ionic current. Analysis of the ionic current reveals the biophysics of co-translocational unfolding and may be used to infer the amino acid sequence and posttranslational modifications of the molecule.

Use of pore-forming toxins to study co-translocational protein folding.

In vivo proteins fold mainly as they emerge from the ribosome or as they emerge from a membrane translocon. Membrane translocation in particular poses technical challenges to the study of the associated protein folding processes. Recently we have developed a single-molecule methodology that allows the capture of a single protein molecule through a membrane translocon with biotinylated oligonucleotides covalently bound at its N- and C- terminus using streptavidin. The resulting rotaxane can be driven forwards and backwards changing the voltage polarity, and carefully planned experiments allow inference of the folding pathway. Here we will discuss the details of a simplified methodological approach.

Proteolytic scanning calorimetry: a novel methodology that probes the fundamental features of protein kinetic stability.

We introduce proteolytic scanning calorimetry, a modification of the differential scanning calorimetry approach to the determination of protein stability in which a proteolytic enzyme (thermolysin) is used to mimic a harsh environment. This methodology allows the straightforward calculation of the rate of irreversible denaturation as a function of temperature and concentration of proteolytic enzyme and, as a result, has the potential to probe efficiently the fundamental biophysical features of protein kinetic stability. In the particular case of Escherichia coli thioredoxin (used as an illustrative example in this article), we find that the rate of irreversible denaturation is determined by 1), the global unfolding mechanism at low thermolysin concentrations, indicating that thermodynamic stability may contribute directly to the kinetic stability of thioredoxin under moderately harsh conditions and 2), the rate of unfolding at high thermolysin concentrations, indicating that the free-energy barrier for unfolding may act as a safety mechanism that ensures significant kinetic stability, even in very harsh environments. This thioredoxin picture, however, is by no means expected to be general and different proteins may show different patterns of kinetic stabilization. Proteolytic scanning calorimetry is particularly well-suited to probe this diversity at a fundamental biophysical level.

Three-way interaction between 14-3-3 proteins, the N-terminal region of tyrosine hydroxylase, and negatively charged membranes.

Tyrosine hydroxylase (TH), the rate-limiting enzyme in the synthesis of catecholamines, is activated by phosphorylation-dependent binding to 14-3-3 proteins. The N-terminal domain of TH is also involved in interaction with lipid membranes. We investigated the binding of the N-terminal domain to its different partners, both in the unphosphorylated (TH-(1-43)) and Ser(19)-phosphorylated (THp-(1-43)) states by surface plasmon resonance. THp-(1-43) showed high affinity for 14-3-3 proteins (K(d) approximately 0.5 microM for 14-3-3gamma and -zeta and 7 microM for 14-3-3eta). The domains also bind to negatively charged membranes with intermediate affinity (concentration at half-maximal binding S(0.5) = 25-58 microM (TH-(1-43)) and S(0.5) = 135-475 microM (THp-(1-43)), depending on phospholipid composition) and concomitant formation of helical structure. 14-3-3gamma showed a preferential binding to membranes, compared with 14-3-3zeta, both in chromaffin granules and with liposomes at neutral pH. The affinity of 14-3-3gamma for negatively charged membranes (S(0.5) = 1-9 microM) is much higher than the affinity of TH for the same membranes, compatible with the formation of a ternary complex between Ser(19)-phosphorylated TH, 14-3-3gamma, and membranes. Our results shed light on interaction mechanisms that might be relevant for the modulation of the distribution of TH in the cytoplasm and membrane fractions and regulation of L-DOPA and dopamine synthesis.

Protein-protein interactions at an enzyme-substrate interface: characterization of transient reaction intermediates throughout a full catalytic cycle of Escherichia coli thioredoxin reductase.

A large collection of structural snapshots along a full catalytic cycle of Escherichia coli thioredoxin reductase (TrxR) has been generated and characterized using a combination of theoretical methods. Molecular models were built starting from the available X-ray crystallographic structures of dimeric wild-type TrxR in the flavin-oxidizing conformation and a C135S TrxR mutant enzyme in a flavin-reducing conformation "trapped" by a cross-link between Cys138 of TrxR and Cys32 of C35S mutant thioredoxin (Trx). The transition between these two extreme states, which is shown to be reproduced in a normal mode analysis, as well as natural cofactor binding and dissociation, were simulated for the wild-type species using unrestrained and targeted molecular dynamics following docking of oxidized Trx to reduced TrxR. The whole set of simulations provides a comprehensive structural framework for understanding the mechanism of disulfide reduction in atomic detail and identifying the most likely intermediates that facilitate entry of NADPH and exit of NADP(+). The crucial role assigned to Arg73 and Lys36 of Trx in substrate binding and complex stabilization was ascertained when R73G, R73D, and K36A site-directed mutants of Trx were shown to be impaired to different extents in their ability to be reduced by TrxR. On the basis of previous findings and the results reported herein, E. coli TrxR appears as a beautifully engineered molecular machine that is capable of synchronizing cofactor capture and ejection with substrate binding and redox activity through an interdomain twisting motion.

Modulation of buried ionizable groups in proteins with engineered surface charge.

Recent work has shown that proteins can tolerate hydrophobic-to-ionizable-residue mutations. Here, we provide experimental evidence that the essential properties (pK value, protonation state, local dynamics) of buried ionizable groups in proteins can be efficiently modulated through the rational design of the surface charge distribution, thus paving the way for the protein engineering exploitation of charge burial.

Free-energy landscapes of membrane co-translocational protein unfolding.

Protein post-translational translocation is found at the plasma membrane of prokaryotes and protein import into organellae. Translocon structures are becoming available, however the dynamics of proteins during membrane translocation remain largely obscure. Here we study, at the single-molecule level, the folding landscape of a model protein while forced to translocate a transmembrane pore. We use a DNA tag to drive the protein into the α-hemolysin pore under a quantifiable force produced by an applied electric potential. Using a voltage-quench approach we find that the protein fluctuates between the native state and an intermediate in the translocation process at estimated forces as low as 1.9 pN. The fluctuation kinetics provide the free energy landscape as a function of force. We show that our stable, ≈15 kBT, substrate can be unfolded and translocated with physiological membrane potentials and that selective divalent cation binding may have a profound effect on the translocation kinetics.

Transmembrane protein rotaxanes reveal kinetic traps in the refolding of translocated substrates.

Understanding protein folding under conditions similar to those found in vivo remains challenging. Folding occurs mainly vectorially as a polypeptide emerges from the ribosome or from a membrane translocon. Protein folding during membrane translocation is particularly difficult to study. Here, we describe a single-molecule method to characterize the folded state of individual proteins after membrane translocation, by monitoring the ionic current passing through the pore. We tag both N and C termini of a model protein, thioredoxin, with biotinylated oligonucleotides. Under an electric potential, one of the oligonucleotides is pulled through a α-hemolysin nanopore driving the unfolding and translocation of the protein. We trap the protein in the nanopore as a rotaxane-like complex using streptavidin stoppers. The protein is subjected to cycles of unfolding-translocation-refolding switching the voltage polarity. We find that the refolding pathway after translocation is slower than in bulk solution due to the existence of kinetic traps.

Beyond Lumry-Eyring: an unexpected pattern of operational reversibility/irreversibility in protein denaturation.

We have found that, contrary to naïve intuition, the degree of operational reversibility in the thermal denaturation of lipase from Thermomyces lanuginosa (an important industrial enzyme) in urea solutions is maximum when the protein is heated several degrees above the end of the temperature-induced denaturation transition. Upon cooling to room temperature, the protein seems to reach a state with enzymatic activity similar to that of the initial native state, but with higher denaturation temperature and radically different behavior in terms of susceptibility to irreversible denaturation. These results show that patterns of operational reversibility/irreversibility in protein denaturation may be more complex than the often-taken-for-granted, two-situation classification (reversible vs. irreversible). Furthermore, they are consistent with the possibility of existence of different native or native-like states separated by high kinetic barriers under native conditions and they suggest experimental procedures to reach and study such "alternative" native states.

Engineering proteins with tunable thermodynamic and kinetic stabilities.

It is widely recognized that enhancement of protein stability is an important biotechnological goal. However, some applications at least, could actually benefit from stability being strongly dependent on a suitable environment variable, in such a way that enhanced stability or decreased stability could be realized as required. In therapeutic applications, for instance, a long shelf-life under storage conditions may be convenient, but a sufficiently fast degradation of the protein after it has performed the planned molecular task in vivo may avoid side effects and toxicity. Undesirable effects associated to high stability are also likely to occur in food-industry applications. Clearly, one fundamental factor involved here is the kinetic stability of the protein, which relates to the time-scale of the irreversible denaturation processes and which is determined to some significant extent by the free-energy barrier for unfolding (the barrier that "separates" the native state from the highly-susceptible-to-irreversible-alterations nonnative states). With an appropriate experimental model, we show that strong environment-dependencies of the thermodynamic and kinetic stabilities can be achieved using robust protein engineering. We use sequence-alignment analysis and simple computational electrostatics to design stabilizing and destabilizing mutations, the latter introducing interactions between like charges which are screened out at high salt. Our design procedures lead naturally to mutating regions which are mostly unstructured in the transition state for unfolding. As a result, the large salt effect on the thermodynamic stability of our consensus plus charge-reversal variant translates into dramatic changes in the time-scale associated to the unfolding barrier: from the order of years at high salt to the order of days at low salt. Certainly, large changes in salt concentration are not expected to occur in biological systems in vivo. Hence, proteins with strong salt-dependencies of the thermodynamic and kinetic stabilities are more likely to be of use in those cases in which high-stability is required only under storage conditions. A plausible scenario is that inclusion of high salt in liquid formulations will contribute to a long protein shelf-life, while the lower salt concentration under the conditions of the application will help prevent the side effects associated with high-stability which may potentially arise in some therapeutic and food-industry applications. From a more general viewpoint, this work shows that consensus engineering and electrostatic engineering can be readily combined and clarifies relevant aspects of the relation between thermodynamic stability and kinetic stability in proteins.

DNA-binding miniproteins based on zinc fingers. Assessment of the interaction using nanopores.

Obtaining artificial proteins that mimic the DNA binding properties of natural transcription factors could open new ways of manipulating gene expression at will. In this context it is particularly interesting to develop simple synthetic systems. Inspired by the modularity of natural transcription factors, we have designed synthetic miniproteins that combine the zinc finger module of the transcription factor GAGA and AT-hook peptide domains. These constructs are capable of binding to composite DNA sequences of up to 14 base pairs with high affinity and good selectivity. In particular, we have synthesized three different chimeras and characterized their DNA binding properties by electrophoresis and fluorescence anisotropy. We have also used, for the first time in the study of peptide-based DNA binders, nanopore force spectroscopy to obtain further data on the DNA interaction.

Role of solvation barriers in protein kinetic stability.

The stability of several protein systems of interest has been shown to have a kinetic basis. Besides the obvious biotechnological implications, the general interest of understanding protein kinetic stability is emphasized by the fact that some emerging molecular approaches to the inhibition of amyloidogenesis focus on the increase of the kinetic stability of protein native states. Lipases are among the most important industrial enzymes. Here, we have studied the thermal denaturation of the wild-type form, four single-mutant variants and two highly stable, multiple-mutant variants of lipase from Thermomyces lanuginosa. In all cases, thermal denaturation was irreversible, kinetically controlled and conformed to the two-state irreversible model. This result supports that the novel molecular-dynamics-focused, directed-evolution approach involved in the preparation of the highly stable variants is successful likely because it addresses kinetic stability and, in particular, because heated molecular dynamics simulations possibly identify regions of disrupted native interactions in the transition state for irreversible denaturation. Furthermore, we find very large mutation effects on activation enthalpy and entropy, which were not accompanied by similarly large changes in kinetic urea m-value. From this we are led to conclude that these mutation effects are associated to some structural feature of the transition state for the irreversible denaturation process that is not linked to large changes in solvent accessibility. Recent computational studies have suggested the existence of solvation/desolvation barriers in at least some protein folding/unfolding processes. We thus propose that a solvation barrier (arising from the asynchrony between breaking of internal contacts and water penetration) may contribute to the kinetic stability of lipase from T. lanuginosa (and, possibly, to the kinetic stability of other proteins as well).

Natural selection for kinetic stability is a likely origin of correlations between mutational effects on protein energetics and frequencies of amino acid occurrences in sequence alignments.

It appears plausible that natural selection constrains, to some extent at least, the stability in many natural proteins. If, during protein evolution, stability fluctuates within a comparatively narrow range, then mutations are expected to be fixed with frequencies that reflect mutational effects on stability. Indeed, we recently reported a robust correlation between the effect of 27 conservative mutations on the thermodynamic stability (unfolding free energy) of Escherichia coli thioredoxin and the frequencies of residues occurrences in sequence alignments. We show here that this correlation likely implies a lower limit to thermodynamic stability of only a few kJ/mol below the unfolding free energy of the wild-type (WT) protein. We suggest, therefore, that the correlation does not reflect natural selection of thermodynamic stability by itself, but of some other factor which is linked to thermodynamic stability for the mutations under study. We propose that this other factor is the kinetic stability of thioredoxin in vivo, since( i) kinetic stability relates to irreversible denaturation, (ii) the rate of irreversible denaturation in a crowded cellular environment (or in a harsh extracellular environment) is probably determined by the rate of unfolding, and (iii) the half-life for unfolding changes in an exponential manner with activation free energy and, consequently, comparatively small free energy effects can have deleterious consequences for kinetic stability. This proposal is supported by the results of a kinetic study of the WT form and the 27 single-mutant variants of E. coli thioredoxin based on the global analyses of chevron plots and equilibrium unfolding profiles determined from double-jump unfolding assays. This kinetic study suggests, furthermore, one of the factors that may contribute to the high activation free energy for unfolding in thioredoxin (required for kinetic stability), namely the energetic optimization of native-state residue environments in regions, which become disrupted in the transition state for unfolding.

Label-Free, Multiplexed, Single-Molecule Analysis of Protein-DNA Complexes with Nanopores.

Protein interactions with specific DNA sequences are crucial in the control of gene expression and the regulation of replication. Single-molecule methods offer excellent capabilities to unravel the mechanism and kinetics of these interactions. Here, we develop a nanopore approach where a target DNA sequence is contained in a hairpin followed by a ssDNA. This system allows DNA-protein complexes to be distinguished from bare DNA molecules as they are pulled through a single nanopore detector, providing both equilibrium and kinetic information. We show that this approach can be used to test the inhibitory effect of small molecules on complex formation and their mechanisms of action. In a proof of concept, we use DNAs with different sequence patterns to probe the ability of the nanopore to distinguish the effects of an inhibitor in a complex mixture of target DNAs and proteins. We anticipate that the use of this technology with arrays of thousands of nanopores will contribute to the development of transcription factor binding inhibitors.

Protein co-translocational unfolding depends on the direction of pulling.

Protein unfolding and translocation through pores occurs during trafficking between organelles, protein degradation and bacterial toxin delivery. In vivo, co-translocational unfolding can be affected by the end of the polypeptide that is threaded into the pore first. Recently, we have shown that co-translocational unfolding can be followed in a model system at the single-molecule level, thereby unravelling molecular steps and their kinetics. Here, we show that the unfolding kinetics of the model substrate thioredoxin, when pulled through an α-haemolysin pore, differ markedly depending on whether the process is initiated from the C terminus or the N terminus. Further, when thioredoxin is pulled from the N terminus, the unfolding pathway bifurcates: some molecules finish unfolding quickly, while others finish ~100 times slower. Our findings have important implications for the understanding of biological unfolding mechanisms and in the application of nanopore technology for the detection of proteins and their modifications.

Single-molecule site-specific detection of protein phosphorylation with a nanopore.

We demonstrate single-molecule, site-specific detection of protein phosphorylation with protein nanopore technology. A model protein, thioredoxin, was phosphorylated at two adjacent sites. Analysis of the ionic current amplitude and noise, as the protein unfolds and moves through an α-hemolysin pore, enables the distinction between unphosphorylated, monophosphorylated and diphosphorylated variants. Our results provide a step toward nanopore proteomics.

Multistep protein unfolding during nanopore translocation.

Cells are divided into compartments and separated from the environment by lipid bilayer membranes. Essential molecules are transported back and forth across the membranes. We have investigated how folded proteins use narrow transmembrane pores to move between compartments. During this process, the proteins must unfold. To examine co-translocational unfolding of individual molecules, we tagged protein substrates with oligonucleotides to enable potential-driven unidirectional movement through a model protein nanopore, a process that differs fundamentally from extension during force spectroscopy measurements. Our findings support a four-step translocation mechanism for model thioredoxin substrates. First, the DNA tag is captured by the pore. Second, the oligonucleotide is pulled through the pore, causing local unfolding of the C terminus of the thioredoxin adjacent to the pore entrance. Third, the remainder of the protein unfolds spontaneously. Finally, the unfolded polypeptide diffuses through the pore into the recipient compartment. The unfolding pathway elucidated here differs from those revealed by denaturation experiments in solution, for which two-state mechanisms have been proposed.