Effect of Sequence and Stereochemistry Reversal on p53 Peptide Mimicry

Peptidomimetics effective in modulating protein-protein interactions and resistant to proteolysis have potential in therapeutic applications. An appealing yet underperforming peptidomimetic strategy is to employ D-amino acids and reversed sequences to mimic a lead peptide conformation, either separately or as the combined retro-inverso peptide. In this work, we examine the conformations of inverse, reverse and retro-inverso peptides of p53(15–29) using implicit solvent molecular dynamics simulation and circular dichroism spectroscopy. In order to obtain converged ensembles for the peptides, we find enhanced sampling is required via the replica exchange molecular dynamics method. From these replica exchange simulations, the D-peptide analogues of p53(15–29) result in a predominantly left-handed helical conformation. When the parent sequence is reversed sequence as either the L-peptide and D-peptide, these peptides display a greater helical propensity, feature reflected by NMR and CD studies in TFE/water solvent. The simulations also indicate that, while approximately similar orientations of the side-chains are possible by the peptide analogues, their ability to mimic the parent peptide is severely compromised by backbone orientation (for D-amino acids) and side-chain orientation (for reversed sequences). A retro-inverso peptide is disadvantaged as a mimic in both aspects, and further chemical modification is required to enable this concept to be used fruitfully in peptidomimetic design. The replica exchange molecular simulation approach adopted here, with its ability to provide detailed conformational insights into modified peptides, has potential as a tool to guide structure-based design of new improved peptidomimetics.     

Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential

During replica exchange molecular dynamics (RexMD) simulations, several replicas of a system are simulated at different temperatures in parallel allowing for exchange between replicas at frequent intervals. This technique allows significantly improved sampling of conformational space and is increasingly being used for structure prediction of peptides and proteins. A drawback of the standard temperature RexMD is the rapid increase of the replica number with increasing system size to cover a desired temperature range. In an effort to limit the number of replicas, a new Hamiltonian-RexMD method has been developed that is specifically designed to enhance the sampling of peptide and protein conformations by applying various levels of a backbone biasing potential for each replica run. The biasing potential lowers the barrier for backbone dihedral transitions and promotes enhanced peptide backbone transitions along the replica coordinate. The application on several peptide cases including in all cases explicit solvent indicates significantly improved conformational sampling when compared with standard MD simulations. This was achieved with a very modest number of 5-7 replicas for each simulation system making it ideally suited for peptide and protein folding simulations as well as refinement of protein model structures in the presence of explicit solvent.     

Peptide folding simulations

Developments in the design of small peptides that mimic proteins in complexity, recent advances in nanosecond time-resolved spectroscopy methods to study peptides and the development of modern, highly parallel simulation algorithms have come together to give us a detailed picture of peptide folding dynamics. Two newly implemented simulation techniques, parallel replica dynamics and replica exchange molecular dynamics, can now describe directly from simulations the kinetics and thermodynamics of peptide formation, respectively. Given these developments, the simulation community now has the tools to verify and validate simulation protocols and models (forcefields).     

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.     

Role of backbone hydration and salt-bridge formation in stability of alpha-helix in solution

We test molecular level hypotheses for the high thermal stability of alpha-helical conformations of alanine-based peptides by performing detailed atomistic simulations of a 20-amino-acid peptide with explicit treatment of water. To assess the contribution of large side chains to alpha-helix stability through backbone desolvation and salt-bridge formation, we simulate the alanine-rich peptide, Ac-YAEAAKAAEAAKAAEAAKAF-Nme, referred to as the EK peptide, that has three pairs of "i, i + 3" glutamic acid(-) and lysine(+) substitutions. Efficient configurational sampling of the EK peptide over a wide temperature range enabled by the replica exchange molecular dynamics technique allows characterization of the stability of alpha-helix with respect to heat-induced unfolding. We find that near ambient temperatures, the EK peptide predominately samples alpha-helical configurations with 80% fractional helicity at 300 K. The helix melts over a broad range of temperatures with melting temperature, T(m), equal to 350 K, that is significantly higher than the T(m) of a 21-residue polyalanine peptide, A(21). Salt-bridges between oppositely charged Glu(-) and Lys(+) side chains can, in principle, provide thermal stability to alpha-helical conformers. For the specific EK peptide sequence, we observe infrequent formation of Glu-Lys salt-bridges (with approximately 10-20% probability) and therefore we conclude that salt-bridge formation does not contribute significantly to the EK peptide's helical stability. However, lysine side chains are found to shield specific "i, i + 4" backbone hydrogen bonds from water, indicating that large side-chain substituents can play an important role in stabilizing alpha-helical configurations of short peptides in aqueous solution through mediation of water access to backbone hydrogen bonds. These observations have implications on molecular engineering of peptides and biomolecules in the design of their thermostable variants where the shielding mechanism can act in concert with other factors such as salt-bridge formation, thereby increasing thermal stability considerably.     

Self-assembly of rationally designed peptides under two-dimensional confinement

The rational design of interfacially confined biomolecules offers a unique opportunity to explore the cooperative relationship among self-assembly, nucleation, and growth processes. This article highlights the role of electrostatics in the self-assembly of β-sheet-forming peptides at the air-water interface. We characterize the phase behavior of a periodically sequenced sheet-forming peptide by using Langmuir techniques, Brewster angle microscopy, attenuated total reflection Fourier transform infrared spectroscopy, and circular dichroism spectroscopy. We find that peptides with an alternating binary sequence transition at high pressures from discrete circular domains to fibrous domains. The qualitative behavior is independent of surface pressure but dependent on molecular areas. In addition, thermodynamic models are employed to specifically quantify differences in electrostatics by obtaining parameters for the critical aggregation area, the limiting molecular area, and the dimensionless ratio of line tension/dipole density. Using these parameters, we are able to relate localized charge distribution to phase transitions, which will allow us to apply these molecules to examine how the dynamics of self-assembly can be directly coupled to the formation of composite nanostructures in biology.     

Mapping the conformational dynamics and pathways of spontaneous steric zipper Peptide oligomerization

The process of protein misfolding and self-assembly into various, polymorphic aggregates is associated with a number of important neurodegenerative diseases. Only recently, crystal structures of several short peptides have provided detailed structural insights into -sheet rich aggregates, known as amyloid fibrils. Knowledge about early events of the formation and interconversion of small oligomeric states, an inevitable step in the cascade of peptide self-assembly, however, remains still limited. We employ molecular dynamics simulations in explicit solvent to study the spontaneous aggregation process of steric zipper peptide segments from the tau protein and insulin in atomistic detail. Starting from separated chains with random conformations, we find a rapid formation of structurally heterogeneous, -sheet rich oligomers, emerging from multiple bimolecular association steps and diverse assembly pathways. Furthermore, our study provides evidence that aggregate intermediates as small as dimers can be kinetically trapped and thus affect the structural evolution of larger oligomers. Alternative aggregate structures are found for both peptide sequences in the different independent simulations, some of which feature characteristics of the known steric zipper conformation (e.g., -sheet bilayers with a dry interface). The final aggregates interconvert with topologically distinct oligomeric states exclusively via internal rearrangements. The peptide oligomerization was analyzed through the perspective of a minimal oligomer, i.e., the dimer. Thereby all observed multimeric aggregates can be consistently mapped onto a space of reduced dimensionality. This novel method of conformational mapping reveals heterogeneous association and reorganization dynamics that are governed by the characteristics of peptide sequence and oligomer size.     

The investigation of the secondary structure propensities and free-energy landscapes of peptide ligands by replica exchange molecular dynamics simulations

The conformational states of two peptide sequences that bind to staphylococcal enterotoxin B are sampled by replica exchange molecular dynamic (REMD) simulations in explicit water. REMD simulations were treated with 52 replicas in the range of 280–501 K for both peptides. The conformational ensembles of both peptides are dominated by random coil, bend and turn structures with a small amount of helical structures for each temperature. In addition, while an insignificant presence of β-bridge structures were observed for both peptides, the β-sheet structure was observed only for peptide 3. The results obtained from simulations at 300 K are consistent with the experimental results obtained from circular dichroism spectroscopy. From the analysis of REMD results, we also calculated hydrophobic and hydrophilic solvent accessible surface areas for both peptides, and it was observed that the hydrophobic segments of the peptides tend to form bend or turn structures. Moreover, the free-energy landscapes of both peptides were obtained by principal component analysis to understand how the secondary structural properties change according to their complex space. From the free-energy analysis, we have found several minima for both peptides at decreased temperature. For these obvious minima of both peptides, it was observed that the random coil, bend and turn structures are still dominant and the helix, β-bridge or β-sheet structures can appear or disappear with respect to minima. On the other hand, when we compare the results of REMD with conventional MD simulations for these peptides, the configurations of peptide 3 might be trapped in energy minima during the conventional MD simulations. Hence, it can be said that the REMD simulations have provided a sufficiently high sampling efficiency.     

The Alzheimer's peptides Abeta40 and 42 adopt distinct conformations in water: a combined MD / NMR study

The role of peptides Abeta40 and Abeta42 in the early pathogenesis of Alzheimer's disease (AD) is frequently emphasized in the literature. It is known that Abeta42 is more prone to aggregation than Abeta40, even though they differ in only two (IA) amino acid residues at the C-terminal end. A direct comparison of the ensembles of conformations adopted by the monomers in solution has been limited by the inherent flexibility of the unfolded peptides. Here, we characterize the conformations of Abeta40 and Abeta42 in water by using a combination of molecular dynamics (MD) and measured scalar (3)J(HNHalpha) data from NMR experiments. We perform replica exchange MD (REMD) simulations and find that classical forcefields reproduce the NMR data quantitatively when the sampling is extended to the microseconds time-scale. Using the quantitative agreement of the NMR data as a validation of the model, we proceed to compare the conformational ensembles of the Abeta40 and Abeta42 peptide monomers. Our analysis confirms the existence of structured regions within the otherwise flexible Abeta peptides. We find that the C terminus of Abeta42 is more structured than that of Abeta40. The formation of a beta-hairpin in the sequence (31)IIGLMVGGVVIA involving short strands at residues 31-34 and 38-41 (in bold) reduces the C-terminal flexibility of the Abeta42 peptide and may be responsible for the higher propensity of this peptide to form amyloids.     

Coiled coils at the edge of configurational heterogeneity. Structural analyses of parallel and antiparallel homotetrameric coiled coils reveal configurational sensitivity to a single solvent-exposed amino acid substitution

A detailed understanding of the mechanisms by which particular amino acid sequences can give rise to more than one folded structure, such as for proteins that undergo large conformational changes or misfolding, is a long-standing objective of protein chemistry. Here, we describe the crystal structures of a single coiled-coil peptide in distinct parallel and antiparallel tetrameric configurations and further describe the parallel or antiparallel crystal structures of several related peptide sequences; the antiparallel tetrameric assemblies represent the first crystal structures of GCN4-derived peptides exhibiting such a configuration. Intriguingly, substitution of a single solvent-exposed residue enabled the parallel coiled-coil tetramer GCN4-pLI to populate the antiparallel configuration, suggesting that the two configurations are close enough in energy for subtle sequence changes to have important structural consequences. We present a structural analysis of the small changes to the helix register and side-chain conformations that accommodate the two configurations and have supplemented these results using solution studies and a molecular dynamics energetic analysis using a replica exchange methodology. Considering the previous examples of structural nonspecificity in coiled-coil peptides, the findings reported here not only emphasize the predisposition of the coiled-coil motif to adopt multiple configurations but also call attention to the associated risk that observed crytstal structures may not represent the only (or even the major) species present in solution.     

Comparison of molecular dynamics simulation methods for amyloid β(1-42) monomers containing D-aspartic acid residues for predicting retention times in chromatography

Molecular dynamics simulations of amyloid β(1-42) containing D-aspartic acid residues were performed using several continuous solvent models to investigate the usefulness of simulation methods for D-amino acid-containing proteins and peptides. Normal molecular dynamics simulations and replica exchange molecular dynamics simulations, which are one of the generalized-ensemble algorithms, were performed. Because the β-structure contents of amyloid β(1-42) peptides obtained by replica exchange molecular dynamics simulations with Onufriev-Bashford-Case generalized Born implicit solvent were qualitatively consistent with experimental data, replica exchange molecular dynamics rather than other methods appeared to be more reasonable for calculations of amyloid β(1-42) containing D-aspartic acid residues. Computational results revealed that peptides with stereoinversion of Asp23 tend to form β-sheet structures by themselves, in contrast to the wild-type peptides that form β-sheet structures only after aggregation. These results are expected to be useful for computational investigations of proteins and peptides such as prediction of retention time of peptides and proteins containing D-aspartic acid residues.     

An inhibitory peptide derived from the α-subunit of the epithelial sodium channel (ENaC) shows a helical conformation

Proteolytic activation of the heteromeric epithelial sodium channel (ENaC) is thought to involve the release of inhibitory peptides from the extracellular domains of its α- and γ-subunit. Recently, we demonstrated that an α-13-mer peptide, corresponding to a putative inhibitory region within the extracellular domain of human αENaC, inhibits human αβγENaC. The aim of the present study was to investigate the structural basis of the inhibitory effect of this α-13-mer peptide. Analysis of the peptide by replica exchange molecular dynamics method, circular dichroism spectroscopy, nuclear magnetic resonance spectroscopy, and molecular dynamics simulations suggested that a helical turn at the carboxy-terminus is the preferred conformational state of the α-13-mer peptide. From this we predicted that a specific mutation (leucine 188 to alanine) should have a strong effect on the conformational preferences of the peptide. To functionally test this, we compared the effect of the wild-type α-13-mer with that of a mutant α-L188A-13-mer on ENaC currents in Xenopus laevis oocytes heterologously expressing human αβγENaC. We demonstrated that replacing the leucine 188 by alanine abolished the inhibitory effect of the α-13-mer peptide on ENaC. These findings suggest that a helical conformation in its carboxyterminal part is functionally important to mediate ENaC inhibition by the α-13-mer peptide. However, high resolution structural information on the complex of the inhibitory αENaC peptide and the channel are needed to confirm this conclusion.     

Insights into nucleic acid conformational dynamics from massively parallel stochastic simulations

The helical hairpin is one of the most ubiquitous and elementary secondary structural motifs in nucleic acids, capable of serving functional roles and participating in long-range tertiary contacts. Yet the self-assembly of these structures has not been well-characterized at the atomic level. With this in mind, the dynamics of nucleic acid hairpin formation and disruption have been studied using a novel computational tool: large-scale, parallel, atomistic molecular dynamics simulation employing an inhomogeneous distributed computer consisting of more than 40,000 processors. Using multiple methodologies, over 500 micro s of atomistic simulation time has been collected for a large ensemble of hairpins (sequence 5'-GGGC[GCAA]GCCU-3'), allowing characterization of rare events not previously observable in simulation. From uncoupled ensemble dynamics simulations in unperturbed folding conditions, we report on 1), competing pathways between the folded and unfolded regions of the conformational space; 2), observed nonnative stacking and basepairing traps; and 3), a helix unwinding-rewinding mode that is differentiated from the unfolding and folding dynamics. A heterogeneous transition state ensemble is characterized structurally through calculations of conformer-specific folding probabilities and a multiplexed replica exchange stochastic dynamics algorithm is used to derive an approximate folding landscape. A comparison between the observed folding mechanism and that of a peptide beta-hairpin analog suggests that although native topology defines the character of the folding landscape, the statistical weighting of potential folding pathways is determined by the chemical nature of the polymer.     

Membrane assembly of simple helix homo-oligomers studied via molecular dynamics simulations

The assembly of simple transmembrane helix homo-oligomers is studied by combining a generalized Born implicit membrane model with replica exchange molecular dynamics simulations to sample the conformational space of various oligomerization states and the native oligomeric conformation. Our approach is applied to predict the structures of transmembrane helices of three proteins--glycophorin A, the M2 proton channel, and phospholamban--using only peptide sequence and the native oligomerization state information. In every case, the methodology reproduces native conformations that are in good agreement with available experimental structural data. Thus, our method should be useful in the prediction of native structures of transmembrane domains of other peptides. When we ignore the experimental constraint on the native oligomerization state and attempt de novo prediction of the structure and oligomerization state based only on sequence and simple energetic considerations, we identify the pentamer as the most stable oligomer for phospholamban. However, for the glycophorin A and the M2 proton channels, we tend to predict higher oligomers as more stable. Our studies demonstrate that reliable predictions of the structure of transmembrane helical oligomers can be achieved when the observed oligomerization state is imposed as a constraint, but that further efforts are needed for the de novo prediction of both structure and oligomeric state.     

NMR and molecular modeling characterization of RGD containing peptides.

The tripeptide sequence arginine-glycine-aspartic acid (RGD) has been shown to be the key recognition segment in numerous cell adhesion proteins. The solution conformation and dynamics in DMSO-d6 of the cyclic pentapeptides, [formula: see text], a potent fibrinogen receptor antagonist, and [formula: see text], a weak fibrinogen receptor antagonist, have been characterized by nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. 1H-1H distance constraints derived from two-dimensional NOE spectroscopy and torsional angle constraints obtained from 3JNH-H alpha coupling constants, combined with computer-assisted modeling using conformational searching algorithms and energy minimization have allowed several low energy conformations of the peptides to be determined. Low temperature studies in combination with molecular dynamics simulations suggest that each peptide does not exist in a single, well-defined conformation, but as an equilibrating mixture of conformers in fast exchange on the NMR timescale. The experimental results can be fit by considering pairs of low energy conformers. Despite this inherent flexibility, distinct conformational preferences were found which may be related to the biological activity of the peptides.     

Structure and energy landscape of a photoswitchable peptide: a replica exchange molecular dynamics study

A replica exchange molecular dynamics (REMD) simulation of a bicyclic azobenzene peptide in explicit dimethyl sulfoxide solution is presented in order to characterize the conformational structures and energy landscape of a photoswitchable peptide. It is shown that an enhanced-sampling technique such as the REMD method is essential to obtain a converged conformational sampling of the peptide at room temperature. This is because conventional MD simulations of less than approximately 100-ns length are either trapped in local minima (at 295 K) or-if run at high temperature-do not resemble the room-temperature REMD results. Calculating various nuclear Overhauser effects (NOEs) and (3)J-couplings, a good overall agreement between the REMD simulations and the NMR experiments of Renner et al. (Biopolymers 2000;54:501-514) is found. In particular, the REMD study confirms the general picture drawn by Renner et al. that the trans-isomer of the azobenzene peptide exhibits a well-defined structure, while the cis-isomer is a conformational heterogeneous system; that is, the trans-isomer occurs in 2 well-defined conformers, while the cis-isomer represents an energetically frustrated system that leads to an ensemble of conformational structures. Employing a principal component analysis of the REMD data, the free energy landscape of the systems is studied at various temperatures. The implications for the folding and unfolding pathways of the system are discussed.     

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.     

Study of the aggregation mechanism of polyglutamine peptides using replica exchange molecular dynamics simulations

Polyglutamine (polyQ, a peptide) with an abnormal repeat length is the causative agent of polyQ diseases, such as Huntington’s disease. Although glutamine is a polar residue, polyQ peptides form insoluble aggregates in water, and the mechanism for this aggregation is still unclear. To elucidate the detailed mechanism for the nucleation and aggregation of polyQ peptides, replica exchange molecular dynamics simulations were performed for monomers and dimers of polyQ peptides with several chain lengths. Furthermore, to determine how the aggregation mechanism of polyQ differs from those of other peptides, we compared the results for polyQ with those of polyasparagine and polyleucine. The energy barrier between the monomeric and dimeric states of polyQ was found to be relatively low, and it was observed that polyQ dimers strongly favor the formation of antiparallel β-sheet structures. We also found a characteristic behavior of the monomeric polyQ peptide: a turn at the eighth residue is always present, even when the chain length is varied. We previously showed that a structure including more than two sets of β-turns is stable, so a long monomeric polyQ chain can act as an aggregation nucleus by forming several pairs of antiparallel β-sheet structures within a single chain. Since the aggregation of polyQ peptides has some features in common with an amyloid fibril, our results shed light on the mechanism for the aggregation of polyQ peptides as well as the mechanism for the formation of general amyloid fibrils, which cause the onset of amyloid diseases.     

Mapping Conformational Ensembles of Aβ Oligomers in Molecular Dynamics Simulations

Although the oligomers formed by Aβ peptides appear to be the primary cytotoxic species in Alzheimer's disease, detailed information about their structures appears to be lacking. In this article, we use exhaustive replica exchange molecular dynamics and an implicit solvent united-atom model to study the structural properties of Aβ monomers, dimers, and tetramers. Our analysis suggests that the conformational ensembles of Aβ dimers and tetramers are very similar, but sharply distinct from those sampled by the monomers. The key conformational difference between monomers and oligomers is the formation of β-structure in the oligomers occurring together with the loss of intrapeptide interactions and helix structure. Our simulations indicate that, independent of oligomer order, the Aβ aggregation interface is largely confined to the sequence region 10-23, which forms the bulk of interpeptide interactions. We show that the fractions of β structure computed in our simulations and measured experimentally are in good agreement.