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

Effects of altered backbone composition on the folding kinetics and mechanism of an ultrafast-folding protein

Sequence-encoded protein folding is a ubiquitous biological process that has been successfully engineered in a range of oligomeric molecules with artificial backbone chemical connectivity. A remarkable aspect of protein folding is the contrast between the rapid rates at which most sequences in nature fold and the vast number of conformational states possible in an unfolded chain with hundreds of rotatable bonds. Research efforts spanning several decades have sought to elucidate the fundamental chemical principles that dictate the speed and mechanism of natural protein folding. In contrast, little is known about how protein mimetic entities transition between an unfolded and folded state. Here, we report effects of altered backbone connectivity on the folding kinetics and mechanism of the B domain of Staphylococcal protein A (BdpA), an ultrafast-folding sequence. A combination of experimental biophysical analysis and atomistic molecular dynamics simulations performed on the prototype protein and several heterogeneous-backbone variants reveal the interplay among backbone flexibility, folding rates, and structural details of the transition state ensemble. Collectively, these findings suggest a significant degree of plasticity in the mechanisms that can give rise to ultrafast folding in the BdpA sequence and provide atomic level insights into how protein mimetic chains adopt an ordered folded state.

Implications of the unfolded state in the folding energetics of heterogeneous-backbone protein mimetics

Sequence-encoded folding is the foundation of protein structure and is also possible in synthetic chains of artificial chemical composition. In natural proteins, the characteristics of the unfolded state are as important as those of the folded state in determining folding energetics. While much is known about folded structures adopted by artificial protein-like chains, corresponding information about the unfolded states of these molecules is lacking. Here, we report the consequences of altered backbone composition on the structure, stability, and dynamics of the folded and unfolded states of a compact helix-rich protein. Characterization through a combination of biophysical experiments and atomistic simulation reveals effects of backbone modification that depend on both the type of artificial monomers employed and where they are applied in sequence. In general, introducing artificial connectivity in a way that reinforces characteristics of the unfolded state ensemble of the prototype natural protein minimizes the impact of chemical changes on folded stability. These findings have implications in the design of protein mimetics and provide an atomically detailed picture of the unfolded state of a natural protein and artificial analogues under non-denaturing conditions.

Biophysical experiments and atomistic simulation reveal impacts of protein backbone alteration on the ensemble that defines the unfolded state. These effects have implications on folded stability of protein mimetics.

Low folding cooperativity of HP35 revealed by single-molecule force spectroscopy and molecular dynamics simulation

Some small proteins, such as HP35, fold at submicrosecond timescale with low folding cooperativity. Although these proteins have been extensively investigated, still relatively little is known about their folding mechanism. Here, using single-molecule force spectroscopy and steered molecule dynamics simulation, we study the unfolding of HP35 under external force. Our results show that HP35 unfolds at extremely low forces without a well-defined unfolding transition state. Subsequently, we probe the structure of unfolded HP35 using the persistence length obtained in the force spectroscopy. We found that the persistence length of unfolded HP35 is around 0.72 nm, >40% longer than typical unstructured proteins, suggesting that there are a significant amount of residual secondary structures in the unfolded HP35. Molecular dynamics simulation further confirmed this finding and revealed that many native contacts are preserved in HP35, even its two ends have been extended up to 8 nm. Our results therefore suggest that retaining a significant amount of secondary structures in the unfolded state of HP35 may be an efficient way to reduce the entropic cost for the formation of tertiary structure and increase the folding speed, although the folding cooperativity is compromised. Moreover, we anticipate that the methods we used in this work can be extended to the study of other proteins with complex folding behaviors and even intrinsically disordered ones.

Secondary structures of native and pathogenic huntingtin N-terminal fragments

Huntington's disease is a neurodegenerative disorder caused by a polyglutamine (polyQ) expansion in the N-terminal fragment of the Huntingtin (Htt) protein. Structural properties of Htt N-terminal regions and the molecular mechanism leading to protein aggregation have not been fully explained yet. We performed all-atom replica exchange molecular dynamics to investigate the structures of Htt N-terminal parts with polyQ tracts of nonpathogenic and pathogenic lengths. The monomers were composed of the headpiece (17 N-terminal residues), a polyQ tract (polyQ(17) for native and polyQ(55) for pathogenic sequence), and a polyP(11) region, followed by 17 amino acids of mixed sequence. We found that corresponding regions in both fragments fold to similar secondary structures; the headpiece and polyQ stretch adopt mainly α-helical conformations, and polyP(11) forms the PP II-type helix. The native N-terminal fragment is more compact and stabilized by hydrophobic interactions between the surface of polyP(11) and the amphipathic helix of the headpiece. In the pathogenic fragment the headpiece is solvent exposed and does not interact with polyP(11). The predicted structure of the native N-terminal fragment agrees with the X-ray structure of the Htt first exon containing polyQ(17). The structure of the pathogenic fragment adheres to an aggregation model that is mediated by the Htt headpiece.

Folding of a helix at room temperature is critically aided by electrostatic polarization of intraprotein hydrogen bonds

We report direct folding of a 17-residue helix protein (pdb:2I9M) by standard molecular dynamics simulation (single trajectory) at room temperature with implicit solvent. Starting from a fully extended structure, 2I9M successfully folds into the native conformation within 16 ns using adaptive hydrogen bond-specific charges to take into account the electrostatic polarization effect. Cluster analysis shows that conformations in the native state cluster have the highest population (78.4%) among all sampled conformations. Folding snapshots and the secondary structure analysis demonstrate that the folding of 2I9M begins at terminals and progresses toward the center. A plot of the free energy landscape indicates that there is no significant free energy barrier during folding, which explains the observed fast folding speed. For comparison, exactly the same molecular dynamics simulation but carried out under existing AMBER charges failed to fold 2I9M into native-like structures. The current study demonstrates that electrostatic polarization of intraprotein hydrogen bonding, which stabilizes the helix, is critical to the successful folding of 2I9m.

Native like structure in the unfolded state of the villin headpiece helical subdomain, an ultrafast folding protein

The villin headpiece subdomain, HP36, is the smallest naturally occurring protein that folds cooperatively. Its small size, rapid folding, and simple three-helix topology have made it an extremely popular system for computational studies of protein folding. The role of unfolded state structure in rapid folding is an area of active investigation, but relatively little is known about the properties of unfolded states under native conditions. A peptide fragment, HP21, which contains the first and second helices of HP36 has been shown to be a good model for structure in the unfolded state of the intact domain but a detailed description of the conformational propensities of HP21 is lacking and the balance between native and nonnative interactions is not known. A series of three-dimensional NMR experiments were performed on (13)C, (15)N-labeled HP21 to investigate in detail its conformational propensities. Analysis of (13)C(alpha), (13)C(beta), (13)CO chemical shifts, Deltadelta(13)C(alpha) - Deltadelta(13)C(beta) secondary shifts, the secondary structure propensity scores, NOEs, (15)N R(2) values and comparison of experimental chemical shifts with those of HP36 and with chemical shifts calculated using the SHIFTS and SHIFTX programs all indicate that there is significant native like structure in the HP21 ensemble, and thus by implication in the unfolded state of HP36.

Folding processes of the B domain of protein A to the native state observed in all-atom ab initio folding simulations

Reaching the native states of small proteins, a necessary step towards a comprehensive understanding of the folding mechanisms, has remained a tremendous challenge to ab initio protein folding simulations despite the extensive effort. In this work, the folding process of the B domain of protein A (BdpA) has been simulated by both conventional and replica exchange molecular dynamics using AMBER FF03 all-atom force field. Started from an extended chain, a total of 40 conventional (each to 1.0 micros) and two sets of replica exchange (each to 200.0 ns per replica) molecular dynamics simulations were performed with different generalized-Born solvation models and temperature control schemes. The improvements in both the force field and solvent model allowed successful simulations of the folding process to the native state as demonstrated by the 0.80 A C(alpha) root mean square deviation (RMSD) of the best folded structure. The most populated conformation was the native folded structure with a high population. This was a significant improvement over the 2.8 A C(alpha) RMSD of the best nativelike structures from previous ab initio folding studies on BdpA. To the best of our knowledge, our results demonstrate, for the first time, that ab initio simulations can reach the native state of BdpA. Consistent with experimental observations, including Phi-value analyses, formation of helix II/III hairpin was a crucial step that provides a template upon which helix I could form and the folding process could complete. Early formation of helix III was observed which is consistent with the experimental results of higher residual helical content of isolated helix III among the three helices. The calculated temperature-dependent profile and the melting temperature were in close agreement with the experimental results. The simulations further revealed that phenylalanine 31 may play critical to achieve the correct packing of the three helices which is consistent with the experimental observation. In addition to the mechanistic studies, an ab initio structure prediction was also conducted based on both the physical energy and a statistical potential. Based on the lowest physical energy, the predicted structure was 2.0 A C(alpha) RMSD away from the experimentally determined structure.

Distal charge transport in peptides

Biological systems often transport charges and reactive processes over substantial distances. Traditional models of chemical kinetics generally do not describe such extreme distal processes. In this Review, an atomistic model for a distal transport of information, which was specifically developed for peptides, is considered. Chemical reactivity is taken as the result of distal effects based on two-step bifunctional kinetics involving unique, very rapid motional properties of peptides in the subpicosecond regime. The bifunctional model suggests highly efficient transport of charge and reactivity in an isolated peptide over a substantial distance; conversely, a very low efficiency in a water environment was found. The model suggests ultrafast transport of charge and reactivity over substantial molecular distances in a peptide environment. Many such domains can be active in a protein.

Two-stage folding of HP-35 from ab initio simulations

Accurate ab initio simulation of protein folding is a critical step toward elucidation of protein-folding mechanisms. Here, we demonstrate highly accurate folding of the 35 residue villin headpiece subdomain (HP35) by all-atom molecular dynamics simulations using AMBER FF03 and the generalized-Born solvation model. In a set of 20 micros long simulations, the protein folded to the native state in multiple trajectories, with the lowest C(alpha) RMSD being 0.39 A for residues 2-34 (excluding residues 1 and 35). The native state had the highest population among all sampled conformations, and the center of most populated cluster had a C(alpha) RMSD of 1.63 A. Folding of this protein can be described as a two-stage process that followed a well-defined pathway. In the first stage, formation of helices II and III as a folding intermediate constituted the rate-limiting step and was initiated at a folding nucleus around residues Phe17 and Pro21. The folding intermediate further acted as a template that facilitated the folding and docking of helix I in the second stage. Detailed descriptions of the folding kinetics and the roles of key residues are presented.

Examination of the folding of a short alanine-based helical peptide with salt bridges using molecular dynamics simulation

A molecular dynamics simulation of the folding of a short alanine-based helical peptide of 17 residues with three Glu...Lys (i, i + 4) salt bridge pairs, referred to as the AEK17 peptide, was carried out. The simulation gave an estimated simulation folding time of 2.5 ns, shorter than 12 ns for an alanine-based peptide of 16 residues with three Lys residues only, referred to as the AK16 peptide, simulated previously. After folded, the AEK17 peptide had a helical content of 77%, in excellent agreement with the experimentally determined value of 80%. An examination of the folding pathways of AEK17 indicated that the peptide proceeded via three-turn helix conformations more than the helix-turn-helix conformation in the folding pathways. An analysis of interactions indicated that the formation of hydrogen bonds between Lys residue side chains and backbone carbonyls is a major factor in the abundant conformation of the three-turn helix intermediate. The substitution of three Ala with Glu residues reduces the extent of hydrophobic interaction in alanine-based AK peptides with the result that the breaking of the interactions of Lys epsilon-NH3+(side chain)...C=O(backbone) is a major activation action for the AEK17 to achieve a complete fold, in contrast to the AK16 peptide, in which breaking non-native hydrophobic interaction is the rate-determining step.

Folding free-energy landscape of villin headpiece subdomain from molecular dynamics simulations

High-accuracy ab initio folding has remained an elusive objective despite decades of effort. To explore the folding landscape of villin headpiece subdomain HP35, we conducted two sets of replica exchange molecular dynamics for 200 ns each and three sets of conventional microsecond-long molecular dynamics simulations, using AMBER FF03 force field and a generalized-Born solvation model. The protein folded consistently to the native state; the lowest C(alpha)-rmsd from the x-ray structure was 0.46 A, and the C(alpha)- rmsd of the center of the most populated cluster was 1.78 A at 300 K. ab initio simulations have previously not reached this level. The folding landscape of HP35 can be partitioned into the native, denatured, and two intermediate-state regions. The native state is separated from the major folding intermediate state by a small barrier, whereas a large barrier exists between the major folding intermediate and the denatured states. The melting temperature T(m) = 339 K extracted from the heat-capacity profile was in close agreement with the experimentally derived T(m) = 342 K. A comprehensive picture of the kinetics and thermodynamics of HP35 folding emerges when the results from replica exchange and conventional molecular dynamics simulations are combined.

Influence of temperature, friction, and random forces on folding of the B-domain of staphylococcal protein A: All-atom molecular dynamics in implicit solvent

The influences of temperature, friction, and random forces on the folding of protein A have been analyzed. A series of all-atom molecular dynamics folding simulations with the Amber ff99 potential and Generalized Born solvation, starting from the fully extended chain, were carried out for temperatures from 300 to 500 K, using (a) the Berendsen thermostat (with no explicit friction or random forces) and (b) Langevin dynamics (with friction and stochastic forces explicitly present in the system). The simulation temperature influences the relative time scale of the major events on the folding pathways of protein A. At lower temperatures, helix 2 folds significantly later than helices 1 and 3. However, with increasing temperature, the folding time of helix 2 approaches the folding times of helices 1 and 3. At lower temperatures, the complete formation of secondary and tertiary structure is significantly separated in time whereas, at higher temperatures, they occur simultaneously. These results suggest that some earlier experimental and theoretical observations of folding events, e.g., the order of helix formation, could depend on the temperature used in those studies. Therefore, the differences in temperature used could be one of the reasons for the discrepancies among published experimental and computational studies of the folding of protein A. Friction and random forces do not change the folding pathway that was observed in the simulations with the Berendsen thermostat, but their explicit presence in the system extends the folding time of protein A.

NMR Characterization of a Peptide Model Provides Evidence for Significant Structure in the Unfolded State of the Villin Headpiece Helical Subdomain

The villin headpiece subdomain (HP36) is the smallest naturally occurring protein that folds cooperatively. The protein folds on a microsecond time scale. Its small size and very rapid folding have made it a popular target for biophysical studies of protein folding. Temperature-dependent one-dimensional (1D) NMR studies of the full-length protein together with CD and 1D NMR studies of the 21-residue peptide fragment (HP21) derived from HP36 have shown that there is significant structure in the unfolded state of HP36 and have demonstrated that HP21 is a good model of these interactions. Here, we characterized the model peptide HP21 in detail by two-dimensional NMR. Strongly upfield shifted C(alpha) protons, the magnitude of the 3J(NH,alpha) coupling constants, and the pattern of backbone-backbone and backbone-side chain NOEs indicate that the ensemble of structures populated by HP21 contains alpha-helical structure and native as well as non-native hydrophobic contacts. The hydrogen-bonded secondary structure inferred from the NOEs is, however, not sufficient to confer significant protection against amide H-D exchange. These studies indicate that there is significant secondary structure and hydrophobic clustering in the unfolded state of HP36. The implications for the folding of HP36 are discussed.