The role of plastic beta-hairpin and weak hydrophobic core in the stability and unfolding of a full sequence design protein

In this study, the thermal stability of a designed alpha/beta protein FSD (full sequence design) was studied by explicit solvent simulations at three moderate temperatures, 273 K, 300 K, and 330 K. The average properties of the ten trajectories at each temperature were analyzed. The thermal unfolding, as judged by backbone root-mean-square deviation and percentage of native contacts, was displayed with increased sampling outside of the native basin as the temperature was raised. The positional fluctuation of the hairpin residues was significantly higher than that of the helix residues at all three temperatures. The hairpin segment displayed certain plasticity even at 273 K. Apart from the terminal residues, the highest fluctuation was shown in the turn residues 7-9. Secondary structure analysis manifested the structural heterogeneity of the hairpin segment. It was also revealed by the simulation that the hydrophobic core was vulnerable to thermal denaturation. Consistent with the experiment, the I7Y mutation in the double mutant FSD-EY (FSD with mutations Q1E and I7Y) dramatically increased the protein stability in the simulation, suggesting that the plasticity of the hairpin can be partially compensated by a stronger hydrophobic core. As for the unfolding pathway, the breathing of the hydrophobic core and the separation of the two secondary structure elements (alpha helix and beta hairpin) was the initiation step of the unfolding. The loss of global contacts from the separation further destabilized the hairpin structure and also led to the unwinding of the helix.     

Stability and Unfolding Mechanism of the N‐terminal β‐Hairpin from [2Fe‐2S] Ferredoxin I by Molecular Dynamics Simulations

The stability and unfolding mechanism of the N‐terminal β‐hairpin of the [2Fe‐2S] ferredoxin I from the blue‐green alga Aphanothece sacrum in pure methanol, 40% (v/v) methanol‐water, and pure water systems were investigated by 10 ns molecular dynamics simulations under periodic boundary conditions. The β‐hairpin was mostly in its native‐like state in pure methanol, whereas it unfolds dramatically following the ‘zip‐up’ mechanism when it was placed in pure water. Both interstrand and inside‐turn hydrogen bonds account for the stability of the β‐hairpin in its native‐like conformation, whereas hydrophobic interactions among nonpolar side chains are responsible for maintaining its stable loop‐like intermediate structures in 40% (v/v) methanol‐water. Reducing solvent polarity seems to increase the stability of the β‐hairpin in its native‐like structure. Methanol is likely to mimic the partially hydrophobic environment around the N‐terminal β‐hairpin by the subsequent α‐helix.     

Folding of a miniprotein with mixed fold

Using the 28 residue betabetaalpha protein FSD-EY as a target system, we examine correction terms for the ECEPP/3 force field. We find an increased probability of formation of the native state at low temperatures resulting from a reduced propensity to form alpha helices and increased formation of beta sheets. Our analysis of the observed folding events suggests that the C-terminal helix of FSD-EY is much more stable than the N-terminal beta hairpin and forms first. The hydrophobic groups of the helix provide a template which promotes the formation of the beta hairpin that is never observed to form without the helix.     

非天然相互作用对均聚多肽链螺旋形成动力学过程的影响的计算机模拟研究

以螺旋结构的形成过程为研究对象,基于粗粒化的格子模型和动态蒙特卡罗模拟方法,初步探讨了非天然氢键相互作用对均聚多肽链螺旋折叠动力学过程的影响.研究发现,非天然氢键的引入虽然延缓了其热力学转变的发生,但也从整体上降低了折叠动力学过程的能垒,在某一特定温度之下,反而可以提高折叠速率.对其折叠路径分布的分析表明,非天然氢键可以减少慢速折叠路径的发生,而后者是导致折叠时间增加的主要因素.另一方面,比较特定温度下多肽链链构象及螺旋片断随时间的演化进程,发现非天然氢键在一定程度上影响了天然氢键的形成以及天然态构象的稳定存在,同时也加快了其部分解折叠过程.这说明,非天然相互作用的存在有利于蛋白质构象的快速动态调整,从而行使其相应的生物功能.     

Denatured-state ensemble and the early-stage folding of the G29A mutant of the B-domain of protein A

The folding mechanism of the G29A mutant of the B-domain of protein A (BdpA) has been studied by all-atom molecular dynamics simulation using AMBER force field (ff03) and generalized Born continuum solvent model. Started from the extended chain conformation, a total of 16 simulations (400 ns each) at 300 K captured some early folding events of the G29A mutant of BdpA. In one of the 16 trajectories, the G29A mutant folded within 2.8 A (root mean square) of the wild-type NMR structure. We observed that the fast burial of hydrophobic residues was the driving force to bring the distant residues into close proximity. The initiation of the helix I and III occurred during the stage of hydrophobic collapse. The initiation and growth of the helix II was slow. Both the secondary structure formation and the development of the native tertiary contacts suggested a multistage folding process. Clustering analysis indicated that two helix species (helices I and III) could be intermediates. Further analysis revealed that the hydrophobic residues of partially folded helix II formed nativelike hydrophobic contacts with helices I and III that stabilized a nativelike state and delayed the completion of folding of the entire protein. The details of the early folding process were compared with other theoretical and experimental studies. It was found that a nativelike hydrophobic cluster was formed by residues including F(30), I(31), L(34), L(44), L(45), and A(48) that prevented further development of the native structures, and breaking the hydrophobic cluster like this one contributed to the rate-limiting step. This was in complete agreement with the recent kinetic measurements in which mutations of these residues to Gly and Ala substantially increased the folding rates by as much as 60 times. Apparently, destabilization of nonnative states dramatically enhanced the folding rates.     

Dynamic folding pathway models of the villin headpiece subdomain (HP‐36) structure

We have investigated the folding pathway of the 36‐residue villin headpiece subdomain (HP‐36) by action‐derived molecular dynamics simulations. The folding is initiated by hydrophobic collapse, after which the concurrent formation of full tertiary structure and α‐helical secondary structure is observed. The collapse is observed to be associated with a couple of specific native contacts contrary to the conventional nonspecific hydrophobic collapse model. Stable secondary structure formation after the collapse suggests that the folding of HP‐36 follows neither the framework model nor the diffusion‐collision model. The C‐terminal helix forms first, followed by the N‐terminal helix positioned in its native orientation. The short middle helix is shown to form last. Signs for multiple folding pathways are also observed. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

Determination of an ensemble of structures representing the denatured state of the bovine acyl-coenzyme a binding protein

The denatured state of a protein contains important information about the determinants of the folding process. By combining site-directed spin-labeling NMR experiments and restrained computer simulations, we have determined ensembles of conformations that represent the denatured state of the bovine acyl-coenzyme A binding protein (ACBP) at three different concentrations of guanidine hydrochloride. As the experimentally determined distance information corresponds to weighted averages over a broad ensemble of structures, we applied the experimental restraints to a system of noninteracting replicas of the protein by using a Monte Carlo sampling scheme. This procedure permits us to sample ensembles of conformations that are compatible with the experimental data and thus to obtain information regarding the distribution of structures in the denatured state. Our results show that the denatured state of ACBP is highly heterogeneous. The high sensitivity of the computational method that we present, however, enabled us to identify long-range interactions between two regions, located near the N- and C-termini, that include both native and non-native elements. The preferential formation of these contacts suggests that the sequence-dependent patterns of helical propensity and hydrophobicity are important determinants of the structure in the denatured state of ACBP.     

A free-energy landscape for coupled folding and binding of an intrinsically disordered protein in explicit solvent from detailed all-atom computations

The N-terminal repressor domain of neural restrictive silencer factor (NRSF) is an intrinsically disordered protein (IDP) that binds to the paired amphipathic helix (PAH) domain of mSin3. An NMR experiment revealed that the minimal binding unit of NRSF is a 15-residue segment that adopts a helical structure upon binding to a cleft of mSin3. We computed a free-energy landscape of this system by an enhanced conformational sampling method, all-atom multicanonical molecular dynamics. The simulation started from a configuration where the NRSF segment was fully disordered and distant from mSin3 in explicit solvent. In the absence of mSin3, the disordered NRSF segment thermally fluctuated between hairpins, helices, and bent structures. In the presence of mSin3, the segment bound to mSin3 by adopting the structures involved in the isolated state, and non-native and native complexes were formed. The free-energy landscape comprised three superclusters, and free-energy barriers separated the superclusters. The native complex was located at the center of the lowest free-energy cluster. When NRSF landed in the largest supercluster, the generated non-native complex moved on the landscape to fold into the native complex, by increasing the interfacial hydrophobic contacts and the helix content. When NRSF landed in other superclusters, the non-native complex overcame the free-energy barriers between the various segment orientations in the binding cleft of mSin3. Both population-shift and induced-fit (or induced-folding) mechanisms work cooperatively in the coupled folding and binding. The diverse structural adaptability of NRSF may be related to the hub properties of the IDP.     

The interplay of turn formation and hydrophobic interactions on the early kinetic events in protein folding

While both turn formation and hydrophobic interactions play dominant roles in the initiation of protein folding, their individual contributions to the folding kinetics and to the structural stability of the protein still remain poorly understood. Here, we applied a photolabile linker to "cage" some important structural motifs, including both α-helices and β-sheets, into their non-native states. These "caged" structural motifs are then relaxed by laser-flash photolysis and their refolding events followed by photoacoustic calorimetry (PAC) and photothermal beam deflection (PBD). These experiments, combined with our previous results, revealed that spontaneous α-helix formation can occur extremely rapidly (10(8)-10(9) s(-1)) if the process is driven solely by turn formation followed by helix propagation. However, if sequestering of the side chains of hydrophobic amino acid residues participates in the refolding process, which may provide additional driving force beyond that afforded by turn formation alone, the refolding rate will be retarded, often by many orders of magnitude. This is usually the case in the formation of three-stranded β-sheets (10(7)-10(8) s(-1)) and β-hairpins (10(5)-10(6) s(-1)). Thus, we propose that proteins take advantage of the hierarchy of timescales associated with either turn formation, hydrophobic interactions, or global collapse of tertiary structure to accomplish the folding process in an orderly fashion, as these events are sufficiently separated in time and do not interfere with one another.     

Folding Dynamics of 10‐Residue β‐Hairpin Peptide Chignolin

Short peptides that fold into β‐hairpins are ideal model systems for investigating the mechanism of protein folding because their folding process shows dynamics typical of proteins. We performed folding, unfolding, and refolding molecular dynamics simulations (total of 2.7 μs) of the 10‐residue β‐hairpin peptide chignolin, which is the smallest β‐hairpin structure known to be stable in solution. Our results revealed the folding mechanism of chignolin, which comprises three steps. First, the folding begins with hydrophobic assembly. It brings the main chain together; subsequently, a nascent turn structure is formed. The second step is the conversion of the nascent turn into a tight turn structure along with interconversion of the hydrophobic packing and interstrand hydrogen bonds. Finally, the formation of the hydrogen‐bond network and the complete hydrophobic core as well as the arrangement of side‐chain–side‐chain interactions occur at approximately the same time. This three‐step mechanism appropriately interprets the folding process as involving a combination of previous inconsistent explanations of the folding mechanism of the β‐hairpin, that the first event of the folding is formation of hydrogen bonds and the second is that of the hydrophobic core, or vice versa.     

Probing the formation of stable tertiary structure in a model miniprotein at atomic resolution: determinants of stability of a helical hairpin

The minimal model system to study the basic principles of protein folding is the hairpin. The formation of beta-hairpins, which are the basic components of antiparallel beta-sheets, has been studied extensively in the past decade, but much less is known about helical hairpins. Here, we probe hairpin formation between a polyproline type-II helix and an alpha-helix as present in the natural miniprotein peptide YY (PYY). Both turn sequence and interactions of aromatic side chains from the C-terminal alpha-helix with the pockets formed by N-terminal Pro residues are shown by site-directed mutagenesis and solution NMR spectroscopy in different solvent systems to be important determinants of backbone dynamics and hairpin stability, suggesting a close analogy with some beta-hairpin structures. It is shown that multiple relatively weak contacts between the helices are necessary for the formation of the helical hairpin studied here, whereas the type-I beta-turn acts like a hinge, which through certain single amino acid substitutions is destabilized such that hairpin formation is completely abolished. Denaturation and renaturation of tertiary structure by temperature or cosolvents were probed by measuring changes of chemical shifts. Folding of PYY is both reversible and cooperative as inferred from the sigmoidal denaturation curves displayed by residues at the interface of the helical hairpin. Such miniproteins thus feature an important hallmark of globular proteins and should provide a convenient system to study basic aspects of helical hairpin folding that are complementary to those derived from studies of beta-hairpins.     

Micelle‐Triggered β‐Hairpin to α‐Helix Transition in a 14‐Residue Peptide from a Choline‐Binding Repeat of the Pneumococcal Autolysin LytA

Choline‐binding modules (CBMs) have a ββ‐solenoid structure composed of choline‐binding repeats (CBR), which consist of a β‐hairpin followed by a short linker. To find minimal peptides that are able to maintain the CBR native structure and to evaluate their remaining choline‐binding ability, we have analysed the third β‐hairpin of the CBM from the pneumococcal LytA autolysin. Circular dichroism and NMR data reveal that this peptide forms a highly stable native‐like β‐hairpin both in aqueous solution and in the presence of trifluoroethanol, but, strikingly, the peptide structure is a stable amphipathic α‐helix in both zwitterionic (dodecylphosphocholine) and anionic (sodium dodecylsulfate) detergent micelles, as well as in small unilamellar vesicles. This β‐hairpin to α‐helix conversion is reversible. Given that the β‐hairpin and α‐helix differ greatly in the distribution of hydrophobic and hydrophilic side chains, we propose that the amphipathicity is a requirement for a peptide structure to interact and to be stable in micelles or lipid vesicles. To our knowledge, this “chameleonic” behaviour is the only described case of a micelle‐induced structural transition between two ordered peptide structures.     

Structure change of β-hairpin induced by turn optimization: an enhanced sampling molecular dynamics simulation study

Our previous study showed that for the tested polypeptides which have similar β-hairpin structures but different sequences, their folding free energy pathways are dominantly determined by the turn conformational propensity. In this study, we study how the turn conformational propensity affects the structure of hairpins. The folding of two mutants of GB1p peptide (GB1m2 and GB1m3), which have the optimized turn sequence ((6)DDATK(11)T?→?(6)NPATG(11)K) with native structures unsolved, were simulated using integrated tempering sampling molecular dynamics simulations and the predicted stable structures were compared to wild-type GB1p. It was observed that the turn optimization of GB1p generates a more favored 5-residue type I(') turn in addition to the 6-residue type I turn in wild-type GB1p. As a result two distinctly different hairpin structures are formed corresponding to the "misfolded" (M) and the "folded" (F) states. M state is a one-residue-shifted asymmetric β-hairpin structure whereas F state has the similar symmetric hairpin structure as wild-type GB1p. The formation of the favored type I(') turn has a small free energy barrier and leads to the shifted β-hairpin structure, following the modified "zipping" model. The presence of disfavored type I turn structure makes the folding of a β-hairpin consistent with the "hydrophobic-core-centric" model. On the other hand, the folding simulations on other two GB1p mutants (GB1r1 and GBr2), which have the position of the hydrophobic core cluster further away from the turn compared to wild-type GB1p, showed that moving the hydrophobic core cluster away from the turn region destabilizes but does not change the hairpin structure. Therefore, the present study showed that the turn conformational propensity is a key factor in affecting not only the folding pathways but also the stable structure of β-hairpins, and the turn conformational change induced by the turn optimization leads to significant changes of β-hairpin structure.     

Folding mechanisms of individual beta-hairpins in a Go model of Pin1 WW domain by all-atom molecular dynamics simulations

This paper examines the folding mechanism of an individual beta-hairpin in the presence of other hairpins by using an off-lattice model of a small triple-stranded antiparallel beta-sheet protein, Pin1 WW domain. The turn zipper model and the hydrophobic collapse model originally developed for a single beta-hairpin in literature is confirmed to be useful in describing beta-hairpins in model Pin1 WW domain. We find that the mechanism for folding a specific hairpin is independent of whether it folds first or second, but the formation process are significantly dependent on temperature. More specifically, beta1-beta2 hairpin folds via the turn zipper model at a low temperature and the hydrophobic collapse model at a high temperature, while the folding of beta2-beta3 hairpin follows the turn zipper model at both temperatures. The change in folding mechanisms is interpreted by the interplay between contact stability (enthalpy) and loop lengths (entropy), the effect of which is temperature dependent.     

Contact pair dynamics during folding of two small proteins: chicken villin head piece and the Alzheimer protein beta-amyloid

The folding of an extended protein to its unique native state requires establishment of specific, predetermined, often distant, contacts between amino acid residue pairs. The dynamics of contact pair formation between various hydrophobic residues during folding of two different small proteins, the chicken villin head piece (HP-36) and the Alzheimer protein beta-amyloid (betaA-40), are investigated by Brownian dynamics (BD) simulations. These two proteins represent two very different classes-HP-36 being globular while betaA-40 is nonglobular, stringlike. Hydropathy scale and nonlocal helix propensity of amino acids are used to model the complex interaction potential among the various amino acid residues. The minimalistic model we use here employs a connected backbone chain of atoms of equal size while an amino acid is attached to each backbone atom as an additional atom of differing sizes and interaction parameters, determined by the characteristics of each amino acid. Even for such simple models, we find that the low-energy structures obtained by BD simulations of both the model proteins mimic the native state of the real protein rather well, with a best root-mean-square deviation of 4.5 A for HP-36. For betaA-40 (where a single well-defined structure is not available), the simulated structures resemble the reported ensemble rather well, with the well-known beta-bend correctly reproduced. We introduce and calculate a contact pair distance time correlation function, C(P) (ij)(t), to quantify the dynamical evolution of the pair contact formation between the amino acid residue pairs i and j. The contact pair time correlation function exhibits multistage dynamics, including a two stage fast collapse, followed by a slow (microsecond long) late stage dynamics for several specific pairs. The slow late stage dynamics is in accordance with the findings of Sali et al. Analysis of the individual trajectories shows that the slow decay is due to the attempt of the protein to form energetically more favorable pair contacts to replace the less favorable ones. This late stage contact formation is a highly cooperative process, involving participation of several pairs and thus entropically unfavorable and expected to face a large free energy barrier. This is because any new pair contact formation among hydrophobic pairs will require breaking of several contacts, before the favorable ones can be formed. This aspect of protein folding dynamics is similar to relaxation in glassy liquids, where also alpha relaxation requires highly cooperative process of hopping. The present analysis suggests that waiting time for the necessary pair contact formation may obey the Poissonian distribution. We also study the dynamics of Forster energy transfer during folding between two tagged amino acid pairs. This dynamics can be studied by fluorescence resonance energy transfer (FRET). It is found that suitably placed donor-acceptor pairs can capture the slow dynamics during folding. The dynamics probed by FRET is predicted to be nonexponential.     

Structural and thermodynamic features of complexes formed by DNA and synthetic polynucleotides with dodecylamine and dodecyltrimethylammonium bromide

Complex formation of native and denatured DNA, single-stranded polyribonucleotides poly(A) and poly(U), as well as double-stranded poly(A).poly(U) with dodecylamine (DDA) and dodecyltrimethylammonium bromide (DTAB) has been studied by UV-, CD-, IR-spectroscopy and fluorescence analysis of hydrophobic probe pyrene. DDA and DTAB were shown to bind cooperatively with DNA and polyribonucleotides, resulting in the formation of complexes containing hydrophobic micelle-like clusters. Critical aggregation concentration (CAC) of DDA and DTAB shifts sharply to lower values (30-50 times) in the presence of DNA and polynucleotides as compared to critical micelle concentration (CMC) of free DDA and DTAB in solution. The analysis of binding isotherms within the frame of the model of cooperative binding of low-molecular ligands to linear polymers allowed us to determine the thermodynamic parameters of complex formation and estimate the contribution of electrostatic interaction of positively charged heads of amphiphiles with negatively charged phosphate groups of DNA and polyribonucleotides, and hydrophobic interaction of aliphatic chains to complex stability. Electrostatic interaction was shown to make the main contribution to the stability of DNA complexes with DDA, while preferential contribution of hydrophobic interactions is characteristic of DTAB complexes with DNA. The opposite effect of DDA and DTAB on the thermal stability of DNA double helix was demonstrated from UV-melting of DNA-while DTAB stabilizes the DNA helix, DDA, to the contrary, destabilizes it. The destabilizing effect of DDA seems to originate from the displacement of intramolecular hydrogen bonds in complementary Watson-Crick A.T and G.C base pairs with intermolecular H-bonds between unsubstituted DDA amino groups and proton-accepting sites of nucleic bases.