Assisted peptide folding by surface pattern recognition

Natively disordered proteins belong to a unique class of biomolecules whose function is related to their flexibility and their ability to adopt desired conformations upon binding to substrates. In some cases these proteins can bind multiple partners, which can lead to distinct structures and promiscuity in functions. In other words, the capacity to recognize molecular patterns on the substrate is often essential for the folding and function of intrinsically disordered proteins. Biomolecular pattern recognition is extremely relevant both in vivo (e.g., for oligomerization, immune response, induced folding, substrate binding, and molecular switches) and in vitro (e.g., for biosensing, catalysis, chromatography, and implantation). Here, we use a minimalist computational model system to investigate how polar/nonpolar patterns on a surface can induce the folding of an otherwise unstructured peptide. We show that a model peptide that exists in the bulk as a molten globular state consisting of many interconverting structures can fold into either a helix-coil-helix or an extended helix structure in the presence of a complementary designed patterned surface at low hydrophobicity (3.7%) or a uniform surface at high hydrophobicity (50%). However, we find that a carefully chosen surface pattern can bind to and catalyze the folding of a natively unfolded protein much more readily or effectively than a surface with a noncomplementary or uniform distribution of hydrophobic residues.     

The In Vivo Association of BiP with Newly Synthesized Proteins Is Dependent on the Rate and Stability of Folding and Not Simply on the Presence of Sequences That Can Bind to BiP

Immunoglobulin heavy chain-binding protein (BiP) is a member of the hsp70 family of chaperones and one of the most abundant proteins in the ER lumen. It is known to interact transiently with many nascent proteins as they enter the ER and more stably with protein subunits produced in stoichiometric excess or with mutant proteins. However, there also exists a large number of secretory pathway proteins that do not apparently interact with BiP. To begin to understand what controls the likelihood that a nascent protein entering the ER will associate with BiP, we have examined the in vivo folding of a murine λI immunoglobulin (Ig) light chain (LC). This LC is composed of two Ig domains that can fold independent of the other and that each possess multiple potential BiP-binding sequences. To detect BiP binding to the LC during folding, we used BiP ATPase mutants, which bind irreversibly to proteins, as “kinetic traps.” Although both the wild-type and mutant BiP clearly associated with the unoxidized variable region domain, we were unable to detect binding of either BiP protein to the constant region domain. A combination of in vivo and in vitro folding studies revealed that the constant domain folds rapidly and stably even in the absence of an intradomain disulfide bond. Thus, the simple presence of a BiP-binding site on a nascent chain does not ensure that BiP will bind and play a role in its folding. Instead, it appears that the rate and stability of protein folding determines whether or not a particular site is recognized, with BiP preferentially binding to proteins that fold slowly or somewhat unstably.     

Insect juvenile hormone binding protein shows ancestral fold present in human lipid-binding proteins

Low molecular weight juvenile hormone binding proteins (JHBPs) are specific carriers of juvenile hormone (JH) in the hemolymph of butterflies and moths. As hormonal signal transmitters, these proteins exert a profound effect on insect development. The crystal structure of JHBP from Galleria mellonella shows an unusual fold consisting of a long α-helix wrapped in a highly curved antiparallel β-sheet. JHBP structurally resembles the folding pattern found in tandem repeats in some mammalian lipid-binding proteins, with similar organization of one cavity and a disulfide bond between the long helix and the β-sheet. JHBP reveals, therefore, an archetypal fold used by nature for hydrophobic ligand binding. The JHBP molecule possesses two hydrophobic cavities. Several lines of experimental evidence conclusively indicate that JHBP binds JH in only one cavity, close to the N- and C-termini, and that this binding induces a structural change. The second cavity, located at the opposite end of the molecule, could bind another ligand.     

Mechanism of induced folding: Both folding before binding and binding before folding can be realized in staphylococcal nuclease mutants

A considerable number of functional proteins are unstructured under physiological condition. These "intrinsically disordered" proteins exhibit induced folding when they bind their targets. The induced folding comprises two elementary processes: folding and binding. Two mechanisms are possible for the induced folding: either folding before binding or binding before folding. We found that these two mechanisms can be distinguished by the target-concentration dependence of folding kinetics. We also created two types of mutants of staphylococcal nuclease showing the different inhibitor-concentration dependence of induced folding kinetics. One mutant obeys the scheme of binding before folding, while the other the folding before binding. This is the first experimental evidence demonstrating that both mechanisms are realized for a single protein. Binding before folding is possible, when the protein lacks essential nonlocal interaction to stabilize the native conformation. The results cast light on the protein folding mechanism involved in the intrinsically disordered proteins.     

GroEL Can Unfold Late Intermediates Populated on the Folding Pathways of Monellin

The modulation of the folding mechanism of the small protein single-chain monellin (MNEI) by the Escherichia coli chaperone GroEL has been studied. In the absence of the chaperone, the folding of monellin occurs via three parallel routes. When folding is initiated in the presence of a saturating concentration of GroEL, only 50-60% of monellin molecules fold completely. The remaining 40-50% of the monellin molecules remain bound to the GroEL and are released only upon addition of ATP. It is shown that the basic folding mechanism of monellin is not altered by the presence of GroEL, but that it occurs via only one of the three available routes when folding is initiated in the presence of saturating concentrations of GroEL. Two pathways become nonoperational because GroEL binds very tightly to early intermediates that populate these pathways in a manner that makes the GroEL-bound intermediates incompetent to fold. This accounts for the monellin molecules that remain GroEL-bound at the end of the folding reaction. The third pathway remains operational because the GroEL-bound early intermediate on this pathway is folding-competent, suggesting that this early intermediate binds to GroEL in a manner that is different from that of the binding of the early intermediates on the other two pathways. It appears, therefore, that the same protein can bind GroEL in more than one way. The modulation of the folding energy landscape of monellin by GroEL occurs because GroEL binds folding intermediates on parallel folding pathways, in different ways, and with different affinities. Moreover, when GroEL is added to refolding monellin at different times after commencement of refolding, the unfolding of two late kinetic intermediates on two of the three folding pathways can be observed. It appears that the unfolding of late folding intermediates is enabled by a thermodynamic coupling mechanism, wherein GroEL binds more tightly to an early intermediate than to a late intermediate on a folding pathway, with preferential binding energy being larger than the stability of the late intermediate. Hence, it is shown that GroEL can inadvertently and passively cause, through its ability to bind different folding intermediates differentially, the unfolding of late productive intermediates on folding pathways, and that its unfolding action is not restricted solely to misfolded or kinetically trapped intermediates.     

Dynamics of protein folding and cofactor binding monitored by single-molecule force spectroscopy

Many proteins in living cells require cofactors to carry out their biological functions. To reach their functional states, these proteins need to fold into their unique three-dimensional structures in the presence of their cofactors. Two processes, folding of the protein and binding of cofactors, intermingle with each other, making the direct elucidation of the folding mechanism of proteins in the presence of cofactors challenging. Here we use single-molecule atomic force microscopy to directly monitor the folding and cofactor binding dynamics of an engineered metal-binding protein G6-53 at the single-molecule level. Using the mechanical stability of different conformers of G6-53 as sensitive probes, we directly identified different G6-53 conformers (unfolded, apo- and Ni²⁺-bound) populated along the folding pathway of G6-53 in the presence of its cofactor Ni²⁺. By carrying out single-molecule atomic force microscopy refolding experiments, we monitored kinetic evolution processes of these different conformers. Our results suggested that the majority of G6-53 folds through a binding-after-folding mechanism, whereas a small fraction follows a binding-before-folding pathway. Our study opens an avenue to utilizing force spectroscopy techniques to probe the folding dynamics of proteins in the presence of cofactors at the single-molecule level, and we anticipated that this method can be used to study a wide variety of proteins requiring cofactors for their function.     

Intrinsically unstructured proteins by design—electrostatic interactions can control binding,folding, and function of a helix‐loop‐helix heterodimer

Intrinsically disordered proteins that exist as unordered monomeric structures in aqueous solution at pH 7 but fold into four‐helix bundles upon binding to recognized polypeptide targets have been designed. NMR and CD spectra of the monomeric polypeptides show the hallmarks of unordered structures, whereas in the bound state they are highly helical. Analytical ultracentrifugation data shows that the polypeptides bind to their targets to form exclusively heterodimers at neutral pH. To demonstrate the relationship between binding, folding, and function, a catalytic site for ester hydrolysis was introduced into an unordered and largely inactive monomer, but that was structured and catalytically active in the presence of a specific polypeptide target. Electrostatic interactions between surface‐exposed residues inhibited the binding and folding of the monomers at pH 7. Charge–charge repulsion between ionizable amino acids was thus found to be sufficient to disrupt binding between polypeptide chains despite their inherent propensities for structure formation and may be involved in the folding and function of inherently disordered proteins in biology. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.     

Mechanism of Membrane Permeation Induced by Synthetic β-Hairpin Peptides

We have investigated the membrane destabilizing properties of synthetic amphiphilic cationic peptides, MAX1 and MAX35, which have the propensity to form β-hairpin structures under certain conditions, and a control non-β-hairpin-forming peptide MAX8V16E. All three peptides bind to liposomes containing a mixture of zwitterionic POPC and negatively charged POPS lipids as determined by Zeta potential measurements. Circular dichroism measurements indicated folding of MAX1 and MAX35 in the presence of the POPC/POPS liposomes, whereas no such folding was observed with MAX8V16E. There was no binding or folding of these peptides to liposomes containing only POPC. MAX1 and MAX35 induced release of contents from negatively charged liposomes, whereas MAX8V16E failed to promote solute release under identical conditions. Thus, MAX1 and MAX35 bind to, and fold at the surface of negatively charged liposomes adopting a lytic conformation. We ruled out leaky fusion as a mechanism of release by including 2 mol % PEG-PE in the liposomes, which inhibits aggregation/fusion but not folding of MAX or MAX-induced leakage. Using a concentration-dependent quenching probe (calcein), we determined that MAX-induced leakage of liposome contents was an all-or-none process. At MAX1 concentrations, which cause release of ∼50% of the liposomes that contain small (R_h <1.5 nm) markers, only ∼15% of those liposomes release a fluorescent dextran of 40 kDa. A multimeric model of the pore is presented based on these results. Atomistic molecular dynamics simulations show that barrels consisting of 10 β-hairpin MAX1 and MAX35 peptides are relatively more stable than MAX8V16E barrels in the bilayer, suggesting that barrels of this size are responsible for the peptides lytic action.     

Tackling Force-Field Bias in Protein Folding Simulations: Folding of Villin HP35 and Pin WW Domains in Explicit Water

The ability to fold proteins on a computer has highlighted the fact that existing force fields tend to be biased toward a particular type of secondary structure. Consequently, force fields for folding simulations are often chosen according to the native structure, implying that they are not truly “transferable.” Here we show that, while the AMBER ff03 potential is known to favor helical structures, a simple correction to the backbone potential (ff03^∗) results in an unbiased energy function. We take as examples the 35-residue α-helical Villin HP35 and 37 residue β-sheet Pin WW domains, which had not previously been folded with the same force field. Starting from unfolded configurations, simulations of both proteins in Amber ff03^∗ in explicit solvent fold to within 2.0 Å RMSD of the experimental structures. This demonstrates that a simple backbone correction results in a more transferable force field, an important requirement if simulations are to be used to interpret folding mechanism.     

A pairwise residue contact area-based mean force potential for discrimination of native protein structure

Background  

Considering energy function to detect a correct protein fold from incorrect ones is very important for protein structure prediction and protein folding. Knowledge-based mean force potentials are certainly the most popular type of interaction function for protein threading. They are derived from statistical analyses of interacting groups in experimentally determined protein structures. These potentials are developed at the atom or the amino acid level. Based on orientation dependent contact area, a new type of knowledge-based mean force potential has been developed.     

Iterative annealing mechanism explains the functions of the GroEL and RNA chaperones

Molecular chaperones are ATP‐consuming machines, which facilitate the folding of proteins and RNA molecules that are kinetically trapped in misfolded states. Unassisted folding occurs by the kinetic partitioning mechanism according to which folding to the native state, with low probability as well as misfolding to one of the many metastable states, with high probability, occur rapidly. GroEL is an all‐purpose stochastic machine that assists misfolded substrate proteins to fold. The RNA chaperones such as CYT‐19, which are ATP‐consuming enzymes, help the folding of ribozymes that get trapped in metastable states for long times. GroEL does not interact with the folded proteins but CYT‐19 disrupts both the folded and misfolded ribozymes. The structures of GroEL and RNA chaperones are strikingly different. Despite these differences, the iterative annealing mechanism (IAM) quantitatively explains all the available experimental data for assisted folding of proteins and ribozymes. Driven by ATP binding and hydrolysis and GroES binding, GroEL undergoes a catalytic cycle during which it samples three allosteric states, T (apo), R (ATP bound), and R^″ (ADP bound). Analyses of the experimental data show that the efficiency of the GroEL–GroES machinery and mutants is determined by the resetting rate k _{R ″ → T} , which is largest for the wild‐type (WT) GroEL. Generalized IAM accurately predicts the folding kinetics of Tetrahymena ribozyme and its variants. Chaperones maximize the product of the folding rate and the steady‐state native state fold by driving the substrates out of equilibrium. Neither the absolute yield nor the folding rate is optimized.     

A New Computational Approach for Mechanical Folding Kinetics of RNA Hairpins

Based on an ensemble of kinetically accessible conformations, we propose a new analytical model for RNA folding kinetics. The model gives populational kinetics, kinetic rates, transition states, and pathways from the rate matrix. Applications of the new kinetic model to mechanical folding of RNA hairpins such as trans-activation-responsive RNA reveal distinct kinetic behaviors in different force regimes, from zero force to forces much stronger than the critical force for the folding-unfolding transition. In the absence of force or a low force, folding can be initiated (nucleated) at any position by forming the first base stack and there exist many pathways for the folding process. In contrast, for a higher force, the folding/unfolding would predominantly proceed along a single zipping/unzipping pathway. Studies for different hairpin-forming sequences indicate that depending on the nucleotide sequence, a kinetic intermediate can emerge in the low force regime but disappear in high force regime, and a new kinetic intermediate, which is absent in the low and high force regimes, can emerge in the medium force range. Variations of the force lead to changes in folding cooperativity and rate-limiting steps. The predicted network of pathways for trans-activation-responsive RNA suggests two parallel dominant pathways. The rate-limiting folding steps (at f = 8 pN) are the formation of specific basepairs that are 2-4 basepairs away from the loop. At a higher force (f = 11 pN), the folding rate is controlled by the formation of the bulge loop. The predicted rates and transition states are in good agreement with the experimental data for a broad force regime.     

Model for intracellular folding of the human immunodeficiency virus type 1 gp120.

The intracellular folding of the human immunodeficiency virus type 1 gp120 has been assessed by analyzing the ability of the glycoprotein to bind to the viral receptor CD4. Pulse-chase experiments revealed that the glycoprotein was initially produced in a conformation that was unable to bind to CD4 and that the protein attained the appropriate tertiary structure for binding with a half-life of approximately 30 min. The protein appears to fold within the rough endoplasmic reticulum, since blocking of transport to the Golgi apparatus by the oxidative phosphorylation inhibitor carbonyl cyanide m-chlorophenylhydrazone did not appear to perturb the folding kinetics of the molecule. The relatively lengthy folding time was not due to modification of the large number of N-linked glycosylation sites on gp120, since inhibition of the first steps in oligosaccharide modification by the inhibitors deoxynojirimycin or deoxymannojirimycin did not impair the CD4-binding activity of the glycoprotein. However, production of the glycoprotein in the presence of tunicamycin and removal of the N-linked sugars by endoglycosidase H treatment both resulted in deglycosylated proteins that were unable to bind to CD4, suggesting in agreement with previous results, that glycosylation contributes to the ability of gp120 to bind to CD4. Interestingly, incomplete endoglycosidase H treatment revealed that a partially glycosylated glycoprotein could bind to the receptor, implying that a subset of glycosylation sites, perhaps some of those conserved in different isolates of human immunodeficiency virus type 1, might be important for binding of the viral glycoprotein to the CD4 receptor.     

Single-molecule dynamics reveals cooperative binding-folding in protein recognition

The study of associations between two biomolecules is the key to understanding molecular function and recognition. Molecular function is often thought to be determined by underlying structures. Here, combining a single-molecule study of protein binding with an energy-landscape–inspired microscopic model, we found strong evidence that biomolecular recognition is determined by flexibilities in addition to structures. Our model is based on coarse-grained molecular dynamics on the residue level with the energy function biased toward the native binding structure (the Go model). With our model, the underlying free-energy landscape of the binding can be explored. There are two distinct conformational states at the free-energy minimum, one with partial folding of CBD itself and significant interface binding of CBD to Cdc42, and the other with native folding of CBD itself and native interface binding of CBD to Cdc42. This shows that the binding process proceeds with a significant interface binding of CBD with Cdc42 first, without a complete folding of CBD itself, and that binding and folding are then coupled to reach the native binding state. The single-molecule experimental finding of dynamic fluctuations among the loosely and closely bound conformational states can be identified with the theoretical, calculated free-energy minimum and explained quantitatively in the model as a result of binding associated with large conformational changes. The theoretical predictions identified certain key residues for binding that were consistent with mutational experiments. The combined study identified fundamental mechanisms and provided insights about designing and further exploring biomolecular recognition with large conformational changes.     

Kinetic partitioning mechanism governs the folding of the third FnIII domain of tenascin-C: evidence at the single-molecule level

Statistical mechanics and molecular dynamics simulations proposed that the folding of proteins can follow multiple parallel pathways on a rugged energy landscape from unfolded state en route to their folded native states. Kinetic partitioning mechanism is one of the possible mechanisms underlying such complex folding dynamics. Here, we use single-molecule atomic force microscopy technique to directly probe the multiplicity of the folding pathways of the third fibronectin type III domain from the extracellular matrix protein tenascin-C (TNfn3). By stretching individual (TNfn3)₈ molecules, we forced TNfn3 domains to undergo mechanical unfolding and refolding cycles, allowing us to directly observe the folding pathways of TNfn3. We found that, after being mechanically unraveled and then relaxed to zero force, TNfn3 follows multiple parallel pathways to fold into their native states. The majority of TNfn3 fold into the native state in a simple two-state fashion, while a small percentage of TNfn3 were found to be trapped into kinetically stable folding intermediate states with well-defined three-dimensional structures. Furthermore, the folding of TNfn3 was also influenced by its neighboring TNfn3 domains. Complex misfolded states of TNfn3 were observed, possibly due to the formation of domain-swapped dimeric structures. Our studies revealed the ruggedness of the folding energy landscape of TNfn3 and provided direct experimental evidence that the folding dynamics of TNfn3 are governed by the kinetic partitioning mechanism. Our results demonstrated the unique capability of single-molecule AFM to probe the folding dynamics of proteins at the single-molecule level.     

Evaluation of models for the evolution of protein sequences and functions under structural constraint

In the field of evolutionary structural genomics, methods are needed to evaluate why genomes evolved to contain the fold distributions that are observed. In order to study the effects of population dynamics in the evolved genomes we need fast and accurate evolutionary models which can analyze the effects of selection, drift and fixation of a protein sequence in a population that are grounded by physical parameters governing the folding and binding properties of the sequence. In this study, various knowledge-based, force field, and statistical methods for protein folding have been evaluated with four different folds: SH2 domains, SH3 domains, Globin-like, and Flavodoxin-like, to evaluate the speed and accuracy of the energy functions. Similarly, knowledge-based and force field methods have been used to predict ligand binding specificity in SH2 domain. To demonstrate the applicability of these methods, the dynamics of evolution of new binding capabilities by an SH2 domain is demonstrated.     

Transition Path Times Measured by Single-Molecule Spectroscopy

The transition path is a tiny fraction of a molecular trajectory during which the free-energy barrier is crossed. It is a single-molecule property and contains all mechanistic information of folding processes of biomolecules such as proteins and nucleic acids. However, the transition path has been difficult to probe because it is short and rarely visited when transitions actually occur. Recent technical advances in single-molecule spectroscopy have made it possible to directly probe transition paths, which has opened up new theoretical and experimental approaches to investigating folding mechanisms. This article reviews recent single-molecule fluorescence and force spectroscopic measurements of transition path times and their connection to both theory and simulations.     

The final stages of folding of the membrane protein bacteriorhodopsin occur by kinetically indistinguishable parallel folding paths that are mediated by pH

The folding of the transmembrane protein bacteriorhodopsin that occurs during the binding of its retinal cofactor is investigated in a membrane-like environment. Changes in the retinal absorption band reveal two transient retinal-protein intermediate states, with apparent absorption maxima at 380 nm and 440 nm, respectively. Studies on a bacteriorhodopsin mutant of Lys216, which cannot bind retinal covalently, add to evidence that retinal is non-covalently bound in these intermediate states. The two retinal-protein intermediates are genuine intermediate states that form in parallel, each with an observed rate constant of 1.1 s-1. Meanwhile no formation of the folded state is detected. Folded bacteriorhodopsin, with all trans retinal covalently bound, forms from both retinal-bound intermediates with the same apparent rate constant of 0.0070 s-1 that is independent of retinal concentration. Retinal isomerisation then occurs with a rate constant of 0.00033 s-1 to give bacteriorhodopsin containing all trans and 13 cis-retinal. These results provide experimental evidence for multiple folding routes for a membrane protein that are pH dependent, with pH conditions determining the apparent folding route. These observed parallel folding paths are kinetically indistinguishable, which contrasts with most other observations of parallel folding pathways where only pathways with different kinetics have been reported. Furthermore, together with previous work, this study shows that bacteriorhodopsin has to populate at least two folding intermediates, during folding in the mixed lipid micelles investigated here, before the final fold is attained.     

Strategies for folding of affinity tagged proteins using GroEL and osmolytes

Obtaining a proper fold of affinity tagged chimera proteins can be difficult. Frequently, the protein of interest aggregates after the chimeric affinity tag is cleaved off, even when the entire chimeric construct is initially soluble. If the attached protein is incorrectly folded, chaperone proteins such as GroEL bind to the misfolded construct and complicate both folding and affinity purification. Since chaperonin/osmolyte mixtures facilitate correct folding from the chaperonin, we explored the possibility that we could use this intrinsic binding reaction to advantage to refold two difficult-to-fold chimeric constructs. In one instance, we were able to recover activity from a properly folded construct after the construct was released from the chaperonin in the presence of osmolytes. As an added advantage, we have also found that this method involving chaperonins can enable researchers to decide (1) if further stabilization of the folded product is required and (2) if the protein construct in question will ever be competent to fold with osmolytes.