Energy landscape of a peptide consisting of alpha-helix, 3(10)-helix, beta-turn, beta-hairpin, and other disordered conformations

The energy landscape of a peptide [Ace-Lys-Gln-Cys-Arg-Glu-Arg-Ala-Nme] in explicit water was studied with a multicanonical molecular dynamics simulation, and the AMBER parm96 force field was used for the energy calculation. The peptide was taken from the recognition helix of the DNA-binding protein, c-MYB: A rugged energy landscape was obtained, in which the random-coil conformations were dominant at room temperature. The CD spectra of the synthesized peptide revealed that it is in the random state at room temperature. However, the 300 K canonical ensemble, Q(300K), contained alpha-helix, 3(10)-helix, beta-turn, and beta-hairpin structures with small but notable probabilities of existence. The complete alpha-helix, imperfect alpha-helix, and random-coil conformations were separated from one another in the conformational space. This means that the peptide must overcome energy barriers to form the alpha-helix. The overcoming process may correspond to the hydrogen-bond rearrangements from peptide-water to peptide-peptide interactions. The beta-turn, imperfect 3(10)-helix, and beta-hairpin structures, among which there are no energy barriers at 300 K, were embedded in the ensemble of the random-coil conformations. Two types of beta-hairpin with different beta-turn regions were observed in Q(300K). The two beta-hairpin structures may have different mechanisms for the beta-hairpin formation. The current study proposes a scheme that the random state of this peptide consists of both ordered and disordered conformations. In contrast, the energy landscape obtained from the parm94 force field was funnel like, in which the peptide formed the helical conformation at room temperature and random coil at high temperature.     

Beta-hairpin folding mechanism of a nine-residue peptide revealed from molecular dynamics simulations in explicit water

The beta-hairpin fold mechanism of a nine-residue peptide, which is modified from the beta-hairpin of alpha-amylase inhibitor tendamistat (residues 15-23), is studied through direct folding simulations in explicit water at native folding conditions. Three 300-nanosecond self-guided molecular dynamics (SGMD) simulations have revealed a series of beta-hairpin folding events. During these simulations, the peptide folds repeatedly into a major cluster of beta-hairpin structures, which agree well with nuclear magnetic resonance experimental observations. This major cluster is found to have the minimum conformational free energy among all sampled conformations. This peptide also folds into many other beta-hairpin structures, which represent some local free energy minimum states. In the unfolded state, the N-terminal residues of the peptide, Tyr-1, Gln-2, and Asn-3, have a confined conformational distribution. This confinement makes beta-hairpin the only energetically favored structure to fold. The unfolded state of this peptide is populated with conformations with non-native intrapeptide interactions. This peptide goes through fully hydrated conformations to eliminate non-native interactions before folding into a beta-hairpin. The folding of a beta-hairpin starts with side-chain interactions, which bring two strands together to form interstrand hydrogen bonds. The unfolding of the beta-hairpin is not simply the reverse of the folding process. Comparing unfolding simulations using MD and SGMD methods demonstrate that SGMD simulations can qualitatively reproduce the kinetics of the peptide system.     

Free energy landscape and folding mechanism of a beta-hairpin in explicit water: a replica exchange molecular dynamics study

The free energy landscape and the folding mechanism of the C-terminal beta-hairpin of protein G is studied by extensive replica exchange molecular dynamics simulations (40 replicas and 340 ns total simulation time), using the GROMOS96 force field and the SPC explicit water solvent. The study reveals that the system preferentially adopts a beta-hairpin structure at biologically important temperatures, and that the helix content is low at all temperatures studied. Representing the free energy landscape as a function of several types of reaction coordinates, four local minima corresponding to the folded, partially folded, molten globule, and unfolded states are identified. The findings suggest that the folding of the beta-hairpin occurs as the sequence: collapse of hydrophobic core --> formation of H-bond --> formation of the turn. Identifying the folded and molten globule states as the main conformations, the free energy landscape of the beta-hairpin is consistent with a two-state behavior with a broad transition state. The temperature dependence of the folding-unfolding transition is investigated in some detail. The enthalpy and entropy jumps at the folding transition temperature are found to be about three times lower than the experimental estimates, indicating that the folding-unfolding transition in silico is less cooperative than its in vitro counterpart.     

beta-hairpin stability and folding: molecular dynamics studies of the first beta-hairpin of tendamistat

The stability and (un)folding of the 19-residue peptide, SCVTLYQSWRYSQADNGCA, corresponding to the first beta-hairpin (residues 10 to 28) of the alpha-amylase inhibitor tendamistat (PDB entry 3AIT) has been studied by molecular dynamics simulations in explicit water under periodic boundary conditions at several temperatures (300 K, 360 K and 400 K), starting from various conformations for simulation lengths, ranging from 10 to 30 ns. Comparison of trajectories of the reduced and oxidized native peptides reveals the importance of the disulphide bridge closing the beta-hairpin in maintaining a proper turn conformation, thereby insuring a proper side-chain arrangement of the conserved turn residues. This allows rationalization of the conservation of those cysteine residues among the family of alpha-amylase inhibitors. High temperature simulations starting from widely different initial configurations (native beta-hairpin, alpha and left-handed helical and extended conformations) begin sampling similar regions of the conformational space within tens of nanoseconds, and both native and non-native beta-hairpin conformations are recovered. Transitions between conformational clusters are accompanied by an increase in energy fluctuations, which is consistent with the increase in heat capacity measured experimentally upon protein folding. The folding events observed in the various simulations support a model for beta-hairpin formation in which the turn is formed first, followed by hydrogen bond formation closing the hairpin, and subsequent stabilization by side-chain hydrophobic interactions.     

Folding-unfolding thermodynamics of a beta-heptapeptide from equilibrium simulations

The thermodynamics of folding and unfolding of a beta-heptapeptide in methanol solution has been studied at four different temperatures, 298 K, 340 K, 350 K, and 360 K, by molecular dynamics simulation. At each of these temperatures, the 50-ns simulations were sufficient to generate an equilibrium distribution between a relatively small number of conformations (approximately 10(2)), showing that, even above the melting temperature (approximately 340 K), the peptide does not randomly sample conformational space. The free energy of folding and the free energy difference between pairs of conformations have been calculated from their relative populations. The experimentally determined folded conformation at 298 K, a left-handed 3(1)-helix, is at each of the four temperatures the predominant conformation, with its probability and average lifetime decreasing with increasing temperature. The most common intermediates of folding and unfolding are also the same at the four temperatures. Paths and rates of interconversion between different conformations have been determined. It has been found that folding can occur through multiple pathways, not necessarily downhill in free energy, although the final step involves a reduced number of intermediates.     

Thermodynamic and kinetic characterization of a beta-hairpin peptide in solution: an extended phase space sampling by molecular dynamics simulations in explicit water

The folding of the amyloidogenic H1 peptide MKHMAGAAAAGAVV taken from the syrian hamster prion protein is explored in explicit aqueous solution at 300 K using long time scale all-atom molecular dynamics simulations for a total simulation time of 1.1 mus. The system, initially modeled as an alpha-helix, preferentially adopts a beta-hairpin structure and several unfolding/refolding events are observed, yielding a very short average beta-hairpin folding time of approximately 200 ns. The long time scale accessed by our simulations and the reversibility of the folding allow to properly explore the configurational space of the peptide in solution. The free energy profile, as a function of the principal components (essential eigenvectors) of motion, describing the main conformational transitions, shows the characteristic features of a funneled landscape, with a downhill surface toward the beta-hairpin folded basin. However, the analysis of the peptide thermodynamic stability, reveals that the beta-hairpin in solution is rather unstable. These results are in good agreement with several experimental evidences, according to which the isolated H1 peptide adopts very rapidly in water beta-sheet structure, leading to amyloid fibril precipitates [Nguyen et al., Biochemistry 1995;34:4186-4192; Inouye et al., J Struct Biol 1998;122:247-255]. Moreover, in this article we also characterize the diffusion behavior in conformational space, investigating its relations with folding/unfolding conditions.     

Essential dynamics of reversible peptide folding: memory-free conformational dynamics governed by internal hydrogen bonds

A principal component analysis has been applied on equilibrium simulations of a beta-heptapeptide that shows reversible folding in a methanol solution. The analysis shows that the configurational space contains only three dense sub-states. These states of relatively low free energy correspond to the "native" left-handed helix, a partly helical intermediate, and a hairpin-like structure. The collection of unfolded conformations form a relatively diffuse cloud with little substructure. Internal hydrogen-bonding energies were found to correlate well with the degree of folding. The native helical structure folds from the N terminus; the transition from the major folding intermediate to the native helical structure involves the formation of the two most C-terminal backbone hydrogen bonds. A four-state Markov model was found to describe transition frequencies between the conformational states within error limits, indicating that memory-effects are negligible beyond the nanosecond time-scale. The dominant native state fluctuations were found to be very similar to unfolding motions, suggesting that unfolding pathways can be inferred from fluctuations in the native state. The low-dimensional essential subspace, describing 69% of the collective atomic fluctuations, was found to converge at time-scales of the order of one nanosecond at all temperatures investigated, whereas folding/unfolding takes place at significantly longer time-scales, even above the melting temperature.     

Understanding the roles of amino acid residues in tertiary structure formation of chignolin by using molecular dynamics simulation

Chignolin is a 10-residue peptide (GYDPETGTWG) that forms a stable beta-hairpin structure in water. However, its design template, GPM12 (GYDDATKTFG), does not have a specific structure. To clarify which amino acids give it the ability to form the beta-hairpin structure, we calculated the folding free-energy landscapes of chignolin, GPM12, and their chimeric peptides using multicanonical molecular dynamics (MD) simulation. Cluster analysis of the conformational ensembles revealed that the native structure of chignolin was the lowest in terms of free energy while shallow local minima were widely distributed in the free energy landscape of GPM12, in agreement with experimental observations. Among the chimeric peptides, GPM12(D4P/K7G) stably formed the same beta-hairpin structure as that of chignolin in the MD simulation. This was confirmed by nuclear magnetic resonance (NMR) spectroscopy. A comparison of the free-energy landscapes showed that the conformational distribution of the Asp3-Pro4 sequence was inherently biased in a way that is advantageous both to forming hydrogen bonds with another beta-strand and to initiating loop structure. In addition, Gly7 helps stabilize the loop structure by having a left-handed alpha-helical conformation. Such a conformation is necessary to complete the loop structure, although it is not preferred by other amino acids. Our results suggest that the consistency between the short-range interactions that determine the local geometries and the long-range interactions that determine the global structure is important for stable tertiary structure formation.     

Free energy landscapes of peptides by enhanced conformational sampling

The free energy landscapes of peptide conformations in water have been observed by the enhanced conformational sampling method, applying the selectively enhanced multicanonical molecular dynamics simulations. The conformations of the peptide dimers, -Gly-Gly-, -Gly-Ala-, -Gly-Ser-, -Ala-Gly-, -Asn-Gly-, -Pro-Gly-, -Pro-Ala-, and -Ala-Ala-, which were all blocked with N-terminal acetyl and C-terminal N-methyl groups, were individually sampled with the explicit TIP3P water molecules. From each simulation trajectory, we obtained the canonical ensemble at 300 K, from which the individual three-dimensional landscape was drawn by the potential of mean force using the three reaction coordinates: the backbone dihedral angle, psi, of the first amino acid, the backbone dihedral angle, phi, of the second amino acid, and the distance between the carbonyl oxygen of the N-terminal acetyl group and the C-terminal amide proton. The most stable state and several meta-stable states correspond to extended conformations and typical beta-turn conformations, and their free energy values were accounted for from the potentials of mean force at the states. In addition, the contributions from the intra-molecular energies of peptides and those from the hydration effects were analyzed. Consequently, the stable beta-turn conformations in the free energy landscape were consistent with the empirically preferred beta-turn types for each amino acid sequence. The thermodynamic values for the hydration effect were decomposed and they correlated well with the empirical values estimated from the solvent accessible surface area of each molecular conformation during the trajectories. The origin of the architecture of protein local fragments was analyzed from the viewpoint of the free energy and its decomposed factors.     

Ligand binding to a high‐energy partially unfolded protein

The conformational energy landscape of a protein determines populations of all possible conformations of the protein and also determines the kinetics of the conversion between the conformations. Interaction with ligands influences the conformational energy landscapes of proteins and shifts populations of proteins in different conformational states. To investigate the effect of ligand binding on partial unfolding of a protein, we use Escherichia coli dihydrofolate reductase (DHFR) and its functional ligand NADP⁺ as a model system. We previously identified a partially unfolded form of DHFR that is populated under native conditions. In this report, we determined the free energy for partial unfolding of DHFR at varying concentrations of NADP⁺ and found that NADP⁺ binds to the partially unfolded form as well as the native form. DHFR unfolds partially without releasing the ligand, though the binding affinity for NADP⁺ is diminished upon partial unfolding. Based on known crystallographic structures of NADP⁺‐bound DHFR and the model of the partially unfolded protein we previously determined, we propose that the adenosine‐binding domain of DHFR remains folded in the partially unfolded form and interacts with the adenosine moiety of NADP⁺. Our result demonstrates that ligand binding may affect the conformational free energy of not only native forms but also high‐energy non‐native forms.     

Comparative analysis of thermal unfolding simulations of RNA recognition motifs (RRMs) of TAR DNA-binding protein 43 (TDP-43)

TAR DNA-binding protein 43 (TDP-43) inclusions have been found in Amyotrophic lateral sclerosis (ALS) and several other neurodegenerative diseases. Many studies suggest the involvement of RNA recognition motifs (RRMs) in TDP-43 proteinopathy. To elucidate the structural stability and the unfolding dynamics of RRMs, we have carried out atomistic molecular dynamics simulations at two different temperatures (300 and 500 K). The simulations results indicate that there are distinct structural differences in the unfolding pathway between the two domains and RRM1 unfolds faster than RRM2 in accordance with the lower thermal stability found experimentally. The unfolding behaviors of secondary structures showed that the α-helix was more stable than β-sheet and structural rearrangements of β-sheets results in formation of additional α-helices. At higher temperature, RRM1 exhibit increased overall flexibility and unfolding than RRM2. The temperature-dependent free energy landscapes consist of multiple metastable states stabilized by non-native contacts and hydrogen bonds in RRM2, thus rendering the RRM2 more prone to misfolding. The structural rearrangements of RRM2 could lead to aberrant protein–protein interactions that may account for enhanced aggregation and toxicity of TDP-43. Our analysis, thus identify the structural and thermodynamic characteristics of the RRMs of TDP-43, which will serve to uncover molecular mechanisms and driving forces in TDP-43 misfolding and aggregation.     

Phi-analysis at the experimental limits: mechanism of beta-hairpin formation

The 37-residue Formin-binding protein, FBP28, is a canonical three-stranded beta-sheet WW domain. Because of its small size, it is so insensitive to chemical denaturation that it is barely possible to determine accurately a denaturation curve, as the transition spans 0-7 M guanidinium hydrochloride (GdmCl). It is also only marginally stable, with a free energy of denaturation of just 2.3 kcal/mol at 10 degrees Celsius so only small changes in energy upon mutation can be tolerated. But these properties and relaxation times for folding of 25 micros-400 micros conspire to allow the rapid acquisition of accurate and reproducible kinetic data for Phi-analysis using classical temperature-jump methods. The transition state for folding is highly polarized with some regions having Phi-values of 0 and others 1, as readily seen in chevron plots, with Phi-values of 0 having the refolding arms overlaying and those of 1 the unfolding arms superimposable. Good agreement is seen with transition state structures identified from independent molecular dynamics (MD) simulations at 60, 75, and 100 degrees Celsius, which allows us to explore further the details of the folding and unfolding pathway of FBP28. The first beta-turn is near native-like in the transition state for folding (experimental) and unfolding (MD and experiment). The simulations show that there are transient contacts between the aromatic side-chains of the beta-strands in the denatured state and that these interactions provide the driving force for folding of the first beta-hairpin of this three-stranded sheet. Only after the backbone hydrogen bonds are formed between beta1 and beta2 does a hydrogen bond form to stabilize the intervening turn, or the first beta-turn.     

Molecular dynamics simulations of beta-hairpin folding

Molecular dynamics simulations of beta-hairpin folding have been carried out with a solvent-referenced potential at 274 K. The model peptide V4DPGV4 formed stable beta-hairpin conformations and the beta-hairpin ratio calculated by the DSSP algorithm was about 56% in the 50-ns simulation. Folding into beta-hairpin conformations is independent of the initial conformations. The simulations provided insights into the folding mechanism. The hydrogen bond often formed in a beta-turn first, and then propagated by forming more hydrogen bonds along the strands. Unfolding and refolding occurred repeatedly during the simulations. Both the hydrogen bonding and the hydrophobic interaction played important roles in forming the ordered structure. Without the hydrophobic effect, stable beta-hairpin conformations did not form in the simulations. With the same energy functions, the alanine-based peptide (AAQAA)3Y folded into helical conformations, in agreement with experiments. Folding into an alpha-helix or a beta-hairpin is amino acid sequence-dependent.     

Structure of the 21-30 fragment of amyloid beta-protein

Folding and self-assembly of the 42-residue amyloid beta-protein (Abeta) are linked to Alzheimer's disease (AD). The 21-30 region of Abeta, Abeta(21-30), is resistant to proteolysis and is believed to nucleate the folding of full-length Abeta. The conformational space accessible to the Abeta(21-30) peptide is investigated by using replica exchange molecular dynamics simulations in explicit solvent. Conformations belonging to the global free energy minimum (the "native" state) from simulation are in good agreement with reported NMR structures. These conformations possess a bend motif spanning the central residues V24-K28. This bend is stabilized by a network of hydrogen bonds involving the side chain of residue D23 and the amide hydrogens of adjacent residues G25, S26, N27, and K28, as well as by a salt bridge formed between side chains of K28 and E22. The non-native states of this peptide are compact and retain a native-like bend topology. The persistence of structure in the denatured state may account for the resistance of this peptide to protease degradation and aggregation, even at elevated temperatures.     

Point mutations in membrane proteins reshape energy landscape and populate different unfolding pathways

Using single-molecule force spectroscopy, we investigated the effect of single point mutations on the energy landscape and unfolding pathways of the transmembrane protein bacteriorhodopsin. We show that the unfolding energy barriers in the energy landscape of the membrane protein followed a simple two-state behavior and represent a manifestation of many converging unfolding pathways. Although the unfolding pathways of wild-type and mutant bacteriorhodopsin did not change, indicating the presence of same ensemble of structural unfolding intermediates, the free energies of the rate-limiting transition states of the bacteriorhodopsin mutants decreased as the distance of those transition states to the folded intermediate states decreased. Thus, all mutants exhibited Hammond behavior and a change in the free energies of the intermediates along the unfolding reaction coordinate and, consequently, their relative occupancies. This is the first experimental proof showing that point mutations can reshape the free energy landscape of a membrane protein and force single proteins to populate certain unfolding pathways over others.     

Energy landscape and dynamics of the beta-hairpin G peptide and its isomers: Topology and sequences

We have investigated free energy landscape [MM/PBSA + normal modes entropy] of permutations in the G peptide (41-56) from the protein G B1 domain by studying six isomers corresponding to moving the hydrophobic cluster along the beta-strands (toward the turn: T1, AGEWTYDDKTFTVTET; T2, GEDTWDYATFTVTKTE; T3, GEDDWTYATFTVTKTE; toward the end: E1, WTYDDAGETKTFTVT; E2, WEYTGDDATKTETFTV; E3, WTYEGDDATKTETFTV). The free energy terms include molecular mechanics energy, Poisson-Boltzmann electrostatic solvation energy, surface area solvation energy, and conformational entropy estimated by using normal mode analysis. From the wild type to T1, then T3, and finally T2, we see a progressively changing energy landscape, toward a less stable beta-hairpin structure. Moving the hydrophobic cluster outside toward the end region causes a greater change in the energy landscape. alpha-Helical instead of a beta-hairpin structure was the most stable form for the E2 isomer. However, no matter how much the sequence changes, for all variants studied, ideal "native" beta-hairpin topologies remain as minima (regardless of whether global or local) in the energy landscape. In general, we find that the energy landscape is dependent on the hydrophobic cluster topology and on the sequence. Our present study indicates that the key is the relative conformational energies of the different conformations. Changes in the sequence strongly modulate the relative stabilities of topologically similar regions in the energy landscape, rather than redefine the topology space. This finding is consistent with a population redistribution in the process of protein folding. The limited variation of topological space, compared with the number of possible sequence changes, may relate to the observation that the number of known protein folds are far less than the sequential allowance.     

Structural and thermodynamics characters of isolated α-syn12 peptide: long-time temperature replica-exchange molecular dynamics in aqueous solution

The structural and thermodynamics characters of α-syn12 (residues 1-12 of the human α-synuclein protein) peptide in aqueous solution were investigated through temperature replica-exchange molecular dynamics (T-REMD) simulations with the GROMOS 43A1 force field. The two independent T-REMD simulations were completed starting from an initial conformational α-helix and an irregular structure, respectively. Each replica was run for 300 ns. The structural and thermodynamics characters were studied based on parameters such as distributions of backbone dihedral angles, free energy surface, stability of folded β-hairpin structure, and favorite conformations. The results showed that the isolated α-syn12 peptide in water adopted four different conformational states: the first state was a β-hairpin ensemble with Turn(9-6) and four hydrogen bonds, the second state was a β-hairpin ensemble with two turns (Turn(9-6) and Turn(5-2)) and three hydrogen bonds, the third state was a disordered structure with both Turn(8-5) and Turn(5-2), and the last state was a π-helix ensemble. Meanwhile, we studied the free energy change of α-syn12 peptide from the unfolded state to the β-hairpin state, which was in good agreement with the experiments and molecular dynamics simulations for some other peptides. We also analyzed the driving force of the peptide transition. The results indicated that the driving forces were high solvent exposure of hydrophobic Leu8 and hydrophobic residues in secondary structure. To our knowledge, this was the first report to study the isolated α-syn12 peptide in water by T-REMD.     

Exploring the influence of hyperthermophilic protein Ssh10b on the stability and conformation of RNA by molecular dynamics simulation

The hyperthermophilic Ssh10b from Sulfolobus shibatae is a member of the Sac10b family, which binds RNA in vivo as a physiological substrate, and it has been postulated to play a key role in chromosomal organization in Archaea. Even though the crystal structure of Ssh10b‐RNA was resolved successively by X‐ray diffraction (Protein Data Bank [PDB] code: 3WBM), the detailed dynamic characteristics of Ssh10b‐RNA are still unclear. In this study, molecular dynamics (MDs) simulations at 6 temperatures (300, 350, 375, 400, 450, and 500 K) and molecular mechanics Generalized‐Born surface area (MM‐GB/SA) free energy calculations were performed to investigate the mechanism of how Ssh10b protects and stabilizes RNA. The simulation results indicate that RNA is stabilized by Ssh10b when the temperature rises up to 375 K. RNA is found to undergo conformational transition between A‐RNA and A′‐RNA when Ssh10b binds to RNA at 3 different temperatures (300, 350, and 375 K). Salt bridges, hydrogen bonds and hydrophobic interactions are observed, and some residues have significant impact on the structural stability of the complex. This study increases our understanding of the dynamics and interaction mechanism of hyperthermophilic proteins and RNA at the atomic level, and offers a model for studying the structural biology of hyperthermophilic proteins and RNA.