Folding and stability of the three-stranded beta-sheet peptide Betanova: insights from molecular dynamics simulations

Development and in vivo evaluation of intranasal formulations of parathyroid hormone (1-34)

For efficient intranasal transport of parathyroid hormone (1-34) [PTH(1-34)], there is a great medical need to investigate permeation enhancers for intranasal formulations. In this study, the development of PTH(1-34) intranasal formulations was conducted. Based on conformation and chemical stability studies, the most preferable aqueous environment was determined to be 0.008 M acetate buffer solution (ABS). Subsequently, citric acid and Kolliphor^® HS·15 were compared as permeation enhancers. The mechanisms of action of citric acid and Kolliphor^® HS·15 were investigated using an in vitro model of nasal mucosa, and Kolliphor^® HS·15 led to higher permeability of fluorescein isothiocyanate-labeled PTH(1-34) (FITC-PTH) by enhancing both the transcellular and paracellular routes. Moreover, citric acid showed severe mucosal toxicity resulting in cilia shedding, while Kolliphor^® HS·15 did not cause obvious mucosa damage. Finally, Kolliphor^® HS·15 was studied as a permeation enhancer using a liquid chromatography tandem mass spectrometry (LC-MS/MS) method. The results showed that 5% and 10% Kolliphor^® HS·15 increased the bioavailability of PTH(1-34) to 14.76% and 30.87%, respectively. In conclusion, an effective and biosafe PTH(1-34) intranasal formulation was developed by using 10% Kolliphor^® HS·15 as a permeation enhancer. Intranasal formulations with higher concentrations of Kolliphor^® HS·15 for higher bioavailability of PTH(1-34) could be further researched.

Detection of protein conformational changes with multilayer graphene nanopore sensors

Detecting conformational change in protein or peptide is imperative in understanding their dynamic function and diagnosing diseases. Existing techniques either rely on ensemble average that lacks the necessary sensitivity or require florescence labeling. Here we propose to discriminate between different protein conformations with multiple layers of graphene nanopore sensors by measuring the effect of protein-produced electrostatic potential (EP) on electric transport. Using conformations of the octapeptide Angiotensin II obtained through molecular dynamics simulations, we show that the EP critically depends on the geometries of constituent atoms and each conformation carries a unique EP signature. We then, using quantum transport simulations, reveal that these characteristic EP profiles cause distinctive modulation to electric charge densities of the graphene nanopores, leading to distinguishable changes in conductivity. Our results open the potential of label-free, single-molecule, and real-time detection of protein conformational changes.

Molecular dynamics simulations to investigate the aggregation behaviors of the Abeta(17-42) oligomers

The amyloid beta-peptides (Abetas) are the main protein components of amyloid deposits in Alzheimer's disease (AD). Detailed knowledge of the structure and assembly dynamics of Abeta is important for the development of properly targeted AD therapeutics. So far, the process of the monomeric Abeta assembling into oligomeric fibrils and the mechanism underlying the aggregation process remain unclear. In this study, several molecular dynamics simulations were conducted to investigate the aggregation behaviors of the Abeta(17-42) oligomers associated with various numbers of monomers (dimer, trimer, tetramer, and pentamer). Our results showed that the structural stability of the Abeta(17-42) oligomers increases with increasing the number of monomer. We further demonstrated that the native hydrophobic contacts are positive correlated with the beta-sheet contents, indicating that hydrophobic interaction plays an important role in maintaining the structural stability of the Abeta(17-42) oligomers, particularly for those associated with more monomers. Our results also showed that the stability of the C-terminal hydrophobic segment 2 (residues 30-42) is higher than that of the N-terminal hydrophobic segment 1 (residues 17-21), suggesting that hydrophobic segment 2 may act as the nucleation site for aggregation. We further identified that Met35 residue initiates the hydrophobic interactions and that the intermolecular contact pairs, Gly33-Gly33 and Gly37-Gly37, form a stable "molecular notch", which may mediate the packing of the beta-sheet involving many other hydrophobic residues during the early stage of amyloid-like fibril formation.

Structural and dynamical analysis of an engineered FhuA channel protein embedded into a lipid bilayer or a detergent belt

Engineered channel proteins are promising nano-components with applications in nanodelivery and nanoreactors technology. Because few of the engineered channel proteins have been crystallized, solution studies based on Neutron Scattering, Circular Dichroism and NMR play a major role. Consequently, the understanding of membrane proteins dynamics in water/detergent solutions or when embedded in a lipid membrane, can clarify how the environment affects protein behavior. In this study, molecular dynamics simulations of the FhuA Escherichia coli outer membrane channel protein and its engineered FhuA Δ1-159 variant have been performed in two different environments: a DNPC (1,2-dinervonyl-sn-glycero-3-phosphocholine) lipid bilayer and a water/OES (N-octyl-2-hydroxyethyl sulfoxide) detergent solution. Furthermore the FhuA Δ1-159 variant has been simulated in the open and closed states, the last induced by the presence of six 3-(2-pyridyldithio)-propionic-acid in the channel inner core. Differences in protein structural and dynamical behavior between the two environments have been found. Considering the FhuA protein characterized by an elliptical-cylindrical symmetry: (a) neither variations on the secondary structure nor axial deformation have been observed in any of the systems; (b) the ellipticity of the channel section (open state) and its fluctuations are enhanced in presence of water/OES, while diminished or suppressed in the DNPC bilayer; (c) the insertion of hydrophobic pyridyl groups into the FhuA Δ1-159 channel (closed state) induces a higher ellipticity in water/OES solution, while shifting to a circular section in the DNPC membrane; (d) the cork domain represented by the first 159 amino acids does not play a major role for protein stability.

Modelling proteins: conformational sampling and reconstruction of folding kinetics

In the last decades biomolecular simulation has made tremendous inroads to help elucidate biomolecular processes in-silico. Despite enormous advances in molecular dynamics techniques and the available computational power, many problems involve long time scales and large-scale molecular rearrangements that are still difficult to sample adequately. In this review we therefore summarise recent efforts to fundamentally improve this situation by decoupling the sampling of the energy landscape from the description of the kinetics of the process. Recent years have seen the emergence of many advanced sampling techniques, which permit efficient characterisation of the relevant family of molecular conformations by dispensing with the details of the short-term kinetics of the process. Because these methods generate thermodynamic information at best, they must be complemented by techniques to reconstruct the kinetics of the process using the ensemble of relevant conformations. Here we review recent advances for both types of methods and discuss their perspectives to permit efficient and accurate modelling of large-scale conformational changes in biomolecules. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

LpxA: a natural nanotube

UDP-N-acetylglucosamine 3-O-acyltransferase is a protein with a left-handed parallel beta-helix, which is a natural nanotube. They are associated with unusual high stability. To identify the reason behind the structural stability of beta-helical nanotubular structure, we have performed a total of 4 mus molecular dynamics simulations of the protein in implicit solvent at four different temperatures and monitored the unfolding pathway. The correlation in movement between different regions of the nanotubular structure has been identified from the dynamical cross-correlation map and contribution of some specific residues towards unfolding transition has been identified by principal component analysis. Difference in stability of the three loop regions has also been characterized. Construction of the unfolding conformational energy landscape identifies the probable intermediates that can appear in the unfolding pathway of the protein.

Amyloid-beta42 oligomer structures from fibrils: a systematic molecular dynamics study

Recent experimental data demonstrate that small, soluble amyloid-beta42 oligomers play an important role in Alzheimer's disease because they exhibit neurotoxic properties and also act as seed for fibril growth. We performed all-atom molecular dynamics simulations in explicit solvent of 0.7 micros in total on five Abeta9-42 oligomers (monomer through pentamer) starting from the fibril conformation. The initial conformation proves to be stable in the trimer to pentamer, and the two parallel in-register beta-sheets as well as the connecting turn are preserved. The dimer undergoes larger conformational changes in its C-terminus, and the predominant conformation detected exhibits an additional antiparallel beta-sheet in one of the subunits. This conformational rearrangement allows efficient shielding of hydrophobic residues from the solvent, which is not possible for a dimer in the fibril conformation. In addition to the presence of the hydrogen bonds in the beta-sheets, the larger oligomers are stabilized by interchain D23-K28 salt bridges, whereas a D23-N27 interaction is found in the dimer. The degree of structural similarity to the fibril conformation detected for the oligomers in the simulation may also offer a structural explanation for the experimental finding that trimers and tetramers act as more potent seeds in fibril formation than dimers because only small conformational changes will be required for fibril growth. The fact that the dimer predominantly exists in conformations distinct from the larger oligomers and the fibril is also interesting for the design of anti-Alzheimer drugs, because it suggests that multiple drugs might be required to target the structurally different neurotoxic oligomers.

PEP-FOLD: an online resource for de novo peptide structure prediction

Rational peptide design and large-scale prediction of peptide structure from sequence remain a challenge for chemical biologists. We present PEP-FOLD, an online service, aimed at de novo modelling of 3D conformations for peptides between 9 and 25 amino acids in aqueous solution. Using a hidden Markov model-derived structural alphabet (SA) of 27 four-residue letters, PEP-FOLD first predicts the SA letter profiles from the amino acid sequence and then assembles the predicted fragments by a greedy procedure driven by a modified version of the OPEP coarse-grained force field. Starting from an amino acid sequence, PEP-FOLD performs series of 50 simulations and returns the most representative conformations identified in terms of energy and population. Using a benchmark of 25 peptides with 9–23 amino acids, and considering the reproducibility of the runs, we find that, on average, PEP-FOLD locates lowest energy conformations differing by 2.6 Å Cα root mean square deviation from the full NMR structures. PEP-FOLD can be accessed at http://bioserv.rpbs.univ-paris-diderot.fr/PEP-FOLD

Probing the nanosecond dynamics of a designed three-stranded beta-sheet with a massively parallel molecular dynamics simulation

Recently a temperature-jump FTIR study of a designed three-stranded sheet showing a fast relaxation time of approximately 140 +/- 20 ns was published. We performed massively parallel molecular dynamics simulations in explicit solvent to probe the structural events involved in this relaxation. While our simulations produce similar relaxation rates, the structural ensemble is broad. We observe the formation of turn structure, but only very weak interaction in the strand regions, which is consistent with the lack of strong backbone-backbone NOEs in previous structural NMR studies. These results suggest that either (D)P(D)P-II folds at time scales longer than 240 ns, or that (D)P(D)P-II is not a well-defined three-stranded beta-sheet. This work also provides an opportunity to compare the performance of several popular forcefield models against one another.

Molecular-dynamics simulations of C- and N-terminal peptide derivatives of GCN4-p1 in aqueous solution

We report the investigation of two 16-residue peptides in aqueous solution by means of molecular-dynamics simulations. The peptides constitute the C- and N-terminal halves of the 33-residue monomer whose dimer constitutes the leucine zipper of the yeast transcriptional activator, denoted GCN4-p1. To examine a hypothesis about coiled-coil formation, in which the C-terminal half contains a helix-formation trigger site absent in the N-terminal half, experimental studies of the two peptides have determined their helix propensities under several conditions of temperature, pH, and salt concentration with circular dichroism. An NMR experiment provides additional evidence. At temperatures of 278 and 325 K and pH 7.5, mixtures of alpha- and pi-helical secondary structure constitute the most probable conformations in both C- and N-terminal halves. A bifurcated salt bridge between Arg25 and Glu22/20 correlates with the structural fluctuations of the C-terminal half. It also exhibits a persistent loop at the N-terminal end involving the side chains of His18 and Glu22, which is reminiscent of helix-capping boxes. Nonreversible unfolding appears to occur abruptly in the Arg25 mutant, suggesting a cooperative event. Analysis does not indicate that the N-terminal half is less stable than the C-terminal half, indicating that 100 ns is too short a period to observe complete unfolding.

Characterization of the conformational space of a triple-stranded beta-sheet forming peptide with molecular dynamics simulations

Molecular dynamics (MD) simulations have been performed on a series of mutants of the 20 amino acid peptide Betanova in order to critically assess the ability of MD simulations to reproduce the folding and stability of small beta-sheet-forming peptides on currently accessible timescales. Simulations were performed in both water and in 40% methanol solution, using an explicit solvent model. The simulations suggest that all mutants adopt a wide range of conformations in solution, that the structures are highly flexible, and that stabilization of compact structures is due to a delicate balance of hydrophobic and polar side-chain interactions. Simulations longer than 100 ns, although not sufficient for a complete thermodynamic and kinetic description of the system, sample an ensemble of compact conformations characterized by the loss of ordered beta-sheet secondary structure. This suggests that no significant free energy barrier separates the different conformations available.

Folding cooperativity in a three-stranded beta-sheet model

The thermodynamic behavior of a previously designed three-stranded beta-sheet was studied via several microseconds of standard and replica exchange molecular dynamics simulations. The system is shown to populate at least four thermodynamic minima, including two partially folded states in which only a single hairpin is formed. Simulated melting curves show different profiles for the C and N-terminal hairpins, consistent with differences in secondary structure content in published NMR and CD/FTIR measurements, which probed different regions of the chain. Individual beta-hairpins that comprise the three-stranded beta-sheet are observed to form cooperatively. Partial folding cooperativity between the component hairpins is observed, and good agreement between calculated and experimental values quantifying this cooperativity is obtained when similar analysis techniques are used. However, the structural detail in the ensemble of conformations sampled in the simulations permits a more direct analysis of this cooperativity than has been performed on the basis of experimental data. The results indicate the actual folding cooperativity perpendicular to strand direction is significantly larger than the lower bound obtained previously.

Interpreting NMR Data for β-Peptides Using Molecular Dynamics Simulations

NMR is one of the most used techniques to resolve structure of proteins and peptides in solution. However, inconsistencies may occur due to the fact that a polypeptide may adopt more than one conformation. Since the NOE distance bounds and (3)J-values used in such structure determination represent a nonlinear average over the total ensemble of conformers, imposition of NOE or (3)J-value restraints to obtain one unique conformation is not an appropriate procedure in such cases. Here, we show that unrestrained MD simulation of a solute in solution using a high-quality force field yields a conformational ensemble that is largely compatible with the experimental NMR data on the solute. Four 100 ns MD simulations of two forms of a nine-residue beta-peptide in methanol at two temperatures produced conformational ensembles that were used to interpret the NMR data on this molecule and resolve inconsistencies between the experimental NOEs. The protected and unprotected forms of the beta-peptide adopt predominantly a 12/10-helix in agreement with the qualitative interpretation of the NMR data. However, a particular NOE was not compatible with this helix indicating the presence of other conformations. The simulations showed that 3(14)()-helical structures were present in the ensemble of the unprotected form and that their presence correlates with the fulfillment of the particular NOE. Additionally, all inter-hydrogen distances were calculated to compare NOEs predicted by the simulations to the ones observed experimentally. The MD conformational ensembles allowed for a detailed and consistent interpretation of the experimental data and showed the small but specific conformational differences between the protected and unprotected forms of the peptide.

Effect of hexafluoroisopropanol alcohol on the structure of melittin: a molecular dynamics simulation study

The molecular mechanism by which HFIP stabilizes the alpha-helical structure of peptides is not well understood. In the present study, we use melittin as a model to gain insight into the details of the atomistic interactions of HFIP with the peptide. We have performed extensive comparative molecular dynamics simulations (up to 100 nsec) in the absence and in the presence of HFIP. In agreement with recent NMR experiments, the simulations show rapid loss of tertiary structure in water at pH 2 but much higher helicity in 35% HFIP. The MD simulations also indicate that melittin adopts a highly dynamic global structure in 35% HFIP solution with two alpha-helical segments sampling a wide range of angular orientations. The analysis of the HFIP distribution shows the tendency of HFIP to aggregate around the peptide, increasing the local cosolvent concentration to more than two times that in the bulk concentration. The correlation of local peptide structure with HFIP coating suggests that displacement of water at the peptide surface is the main contribution of HFIP in stabilizing the secondary structure of melittin. Finally, a stabilizing effect promoted by the presence of counter-ions was also observed in the simulations.

The role of calcium in the conformational dynamics and thermal stability of the D-galactose/D-glucose-binding protein from Escherichia coli

We have characterized stability and conformational dynamics of the calcium depleted D-galactose/D-glucose-binding protein (GGBP) from Escherichia coli. The structural stability of the protein was investigated by steady state and time resolved fluorescence, and far-UV circular dichroism in the temperature range from 20 degrees C to 70 degrees C. We have found that the absence of the Ca(2+) ion results in a significant destabilization of the C-terminal domain of the protein. In particular, the melting temperature decreases by about 10 degrees C with the simultaneous loss of the melting cooperativity. Time resolved fluorescence quenching revealed significant loosening of the protein when highly shielded Trp residue(s) became accessible to acrylamide at higher temperatures. We have documented a significant stabilizing effect of glucose that mostly reverts the effect of calcium, that is, the thermal stability of the protein increases by about 10 degrees C and the melting cooperativity is restored. Moreover, the protein structure remains compact with low amplitude of the segmental mobility up to high temperatures. We have used molecular dynamics to identify the structural feature responsible for changes in the temperature stability. Disintegration of the Ca(2+)-binding loop seems to be responsible for the loss of the stability in the absence of calcium. The new insights on the structural properties and temperature stability of the calcium depleted GGBP contribute to better understanding of the protein function and constitute important information for the development of new biotechnological applications of this class of proteins.

Principles of carbopeptoid folding: a molecular dynamics simulation study

The conformational spaces of five oligomers of tetrahydrofuran-based carbopeptoids in chloroform and dimethyl sulfoxide were investigated through nine molecular dynamics simulations. Prompted by nuclear magnetic resonance experiments that indicated various stable folds for some but not all of these carbopeptoids, their folding behaviour was investigated as a function of stereochemistry, chain length and solvent. The conformational distributions of these molecules were analysed in terms of occurrence of hydrogen bonds, backbone torsional-angle distributions, conformational clustering and solute configurational entropy. While a cis-linkage across the tetrahydrofuran ring favours right-handed helical structures, a trans-linkage results in a larger conformational variability. Intra-solute hydrogen bonding is reduced with increasing chain length and with increasing solvent polarity. Solute configurational entropies confirm the picture obtained: they are smaller for cis- than for trans-linked peptides, for chloroform than for dimethyl sulfoxide as solvent and for shorter peptide chains. The simulations provide an atomic picture of molecular conformational variability that is consistent with the available experimental data.

Sequence Dependence of Amyloid Fibril Formation: Insights from Molecular Dynamics Simulations

The clarification of the physico-chemical determinants underlying amyloid deposition is critical for our understanding of misfolding diseases. With this purpose we have performed a systematic all-atom molecular dynamics (MD) study of a series of single point mutants of the de novo designed amyloidogenic peptide STVIIE. Sixteen different 50ns long simulations using explicit solvent have been carried out starting from four different conformations of a polymeric six-stranded beta-sheet. The simulations have provided evidence for the influence of a small number of site-specific hydrophobic interactions on the packing and stabilization of nascent aggregates, as well as the interplay between side-chain interactions and the net charge of the molecule on the strand arrangement of polymeric beta-sheets. This MD analysis has also shed light into the origin of the position dependence on mutation of beta-sheet polymerization that was found experimentally for this model system. Our results suggest that MD can be applied to detect critical positions for beta-sheet aggregation within a given amyloidogenic stretch. Studies similar to the one presented here can guide site-directed mutations or the design of drugs that specifically disrupt the key stabilizing interactions of beta-sheet aggregates.

Folding thermodynamics of peptides

A simplified interaction potential for protein folding studies at the atomic level is discussed and tested on a set of peptides with approximately 20 residues each. The test set contains both alpha-helical (Trp cage, F(s)) and beta-sheet (GB1p, GB1m2, GB1m3, Betanova, LLM) peptides. The model, which is entirely sequence-based, is able to fold these different peptides for one and the same choice of model parameters. Furthermore, the melting behavior of the peptides is in good quantitative agreement with experimental data. Apparent folded populations obtained using different observables are compared, and are found to be very different for some of the peptides (e.g., Betanova). In other cases (in particular, GB1m2 and GB1m3), the different estimates agree reasonably well, indicating a more two-state-like melting behavior.

Folding thermodynamics of three beta-sheet peptides: a model study

We study the folding thermodynamics of a beta-hairpin and two three-stranded beta-sheet peptides using a simplified sequence-based all-atom model, in which folding is driven mainly by backbone hydrogen bonding and effective hydrophobic attraction. The native populations obtained for these three sequences are in good agreement with experimental data. We also show that the apparent native population depends on which observable is studied; the hydrophobicity energy and the number of native hydrogen bonds give different results. The magnitude of this dependence matches well with the results obtained in two different experiments on the beta-hairpin.

Molecular dynamics simulations of peptides containing an unnatural amino acid: dimerization, folding, and protein binding

We have performed molecular dynamics (MD) simulations to study the dimerization, folding, and binding to a protein of peptides containing an unnatural amino acid. NMR studies have shown that the substitution of one residue in a tripeptide beta-strand by the unnatural amino acid Hao (5-HO2CCONH-2-MeO-C6H3-CO-NHNH2) modifies the conformational flexibility of the beta-strand and the hydrogen-bonding properties of its two edges: The number of hydrogen-bond donors and acceptors increases at one edge, whereas at the other, they are sterically hindered. In simulations in chloroform, the Hao-containing peptide 9 (i-PrCO-Phe-Hao-Val-NHBu) forms a beta-sheet-like hydrogen-bonded dimer, in good agreement with the available experimental data. Addition of methanol to the solution induces instability of this beta-sheet, as confirmed by the experiments. MD simulations also reproduce the folding of the synthetic peptide 1a (i-PrCO-Hao-Ut-Phe-Ile-Leu-NHMe) into a beta-hairpin-like structure in chloroform. Finally, the Hao-containing peptide, Ac-Ala-Hao-Ala-NHMe, is shown to form a stable complex with the Ras analogue, Rap1 A, in water at room temperature. Together with the available experimental data, these simulation studies indicate that Hao-containing peptides may serve as inhibitors of beta-sheet interactions between proteins.

Orientation restraints in molecular dynamics simulations using time and ensemble averaging

In this article we present methodology for simulating protein dynamics while imposing restraints derived from NMR measurements on partially ordered molecules. Such measurements may include residual dipolar couplings and chemical-shift anisotropies. We define a restraint potential for use in molecular dynamics and energy minimization. The presented potential is consistent with the simultaneously optimized molecular order tensor. Restraining can be performed with time and ensemble averaging. We performed a large number of molecular dynamics simulations of the histidine containing phosphocarrier protein with restraints on backbone N-H vector orientations derived from residual dipolar couplings. From these simulations it is evident that the use of time- or ensemble-averaged restraints is essential to leave the fluctuations of the restrained vectors unaffected. Without averaging the fluctuations of the restrained vectors are reduced significantly. This also has the effect of decreasing the apparent molecular order-parameter tensor.

Possible locally driven folding pathways of TC5b, a 20-residue protein

A novel computational procedure for modeling possible locally driven folding pathways by stepwise elongations of the peptide chain was successfully applied to TC5b, a 20-residue miniprotein. Systematic exploration of the possible locally driven pathways showed that the Trp-cage structure of TC5b could be obtained by stepwise elongation starting from the noncentral local nucleation centers preexisting in the unfolded state of TC5b. The probable locally driven folding pathway starts with folding of alpha-helical fragment 4-9, followed by formation of the proper three-dimensional structure of fragment 4-12, and then 4-18. Accordingly, the Trp-cage-forming interactions emerge successively, first Trp(6)-Pro(12), then Trp(6)-Pro(18), and then Trp(6)-Tyr(3). The Trp-cage-like structures of TC5b found in this study by independent energy calculations are in excellent agreement with the NMR experimental data. The same procedure rationalizes the incomplete Trp-cage formation observed for two analogs of TC5b. Generally, the success of this novel approach is encouraging and provides some justification for the use of computational simulations of locally driven protein folding.

The structural basis for biphasic kinetics in the folding of the WW domain from a formin-binding protein: lessons for protein design?

The mechanism of formation of beta-sheets is of great importance because of the significant role of such structures in the initiation and propagation of amyloid diseases. In this study we examine the folding of a series of three-stranded antiparallel beta-sheets known as WW domains. Whereas other WW domains have been shown to fold with single-exponential kinetics, the WW domain from murine formin-binding protein 28 has recently been shown to fold with biphasic kinetics. By using a combination of kinetics and thermodynamics to characterize a simple model for this protein, the origins of the biphasic kinetics is found to lie in the fact that most of the protein is able to fold without requiring one of the beta-hairpins to be correctly registered. The correct register of this hairpin is enforced by a surface-exposed hydrophobic contact, which is not present in other WW domains. This finding suggests the use of judiciously chosen surface-exposed hydrophobic pairs as a protein design strategy for enforcing the desired strand registry.

Interplay between hydrophobic cluster and loop propensity in beta-hairpin formation: a mechanistic study

We investigated the structural determinants of the stability of a designed beta-hairpin containing a natural hydrophobic cluster from the protein GB1 and a D-Pro-Gly turn forming sequence. The results of our simulations shed light on the factors leading to an ordered secondary structure in a model peptide: in particular, the importance of the so-called diagonal interactions in forming a stable hydrophobic nucleus in the beta-hairpin, together with the more obvious lateral interactions, is examined. With the use of long timescale MD simulations in explicit water, we show the role of diagonal interactions in driving the peptide to the correct folded structure (formation of the hydrophobic core with Trp 2, Tyr 4, and Phe 9 in the first stages of refolding) and in keeping it in the ensemble of folded conformations. The combination of the stabilizing effects of the D-Pro-Gly turn sequence and of the hydrophobic nucleus formation thus favors the attainment of an ordered secondary structure compatible with the one determined experimentally. Moreover, our data underline the importance of the juxtapositions of the side chains of amino acids not directly facing each other in the three-dimensional structure. The combination of these interactions forces the peptide to sample a nonrandom portion of the conformational space, as can be seen in the rapid collapse to an ordered structure in the refolding simulation, and shows that the unfolded state can be closely correlated to the folded ensemble of structures, at least in the case of small model peptides.

Response of SCP-2L domain of human MFE-2 to ligand removal: binding site closure and burial of peroxisomal targeting signal

In the study of the structure and function relationship of human MFE-2, we have investigated the dynamics of human MFE-2SCP-2L (hSCP-2L) and its response to ligand removal. A comparison was made with homologous rabbit SCP-2. Breathing and a closing motion are found, identifiable with an adjustment in size and a closing off of the binding pocket. Crucial residues for structural integrity have been identified. Particularly mobile areas of the protein are loop 1 that is connecting helices A and C in space, and helix D, next to the entrance of the pocket. In hSCP-2L, the binding pocket gets occupied by Phe93, which is making a tight hydrophobic contact with Trp36. In addition, it is found that the C-terminal peroxisomal targeting signal (PTS1) that is solvent exposed in the complexed structure becomes buried when no ligand is present. Moreover, an anti-correlation exists between burial of PTS1 and the size of the binding pocket. The results are in accordance with plant nsLTPs, where a similar accommodation of binding pocket size was found after ligand binding/removal. Furthermore, the calculations support the suggestion of a ligand-assisted targeting mechanism.

Mechanism by which 2,2,2-trifluoroethanol/water mixtures stabilize secondary-structure formation in peptides: a molecular dynamics study

Molecular dynamics simulation techniques have been used to investigate the effect of 2,2,2-trifluoroethanol (TFE) as a cosolvent on the stability of three different secondary structure-forming peptides: the alpha-helix from Melittin, the three-stranded beta-sheet peptide Betanova, and the beta-hairpin 41-56 from the B1 domain of protein G. The peptides were studied in pure water and 30% (vol/vol) TFE/water mixtures at 300 K. The simulations suggest that the stabilizing effect of TFE is induced by the preferential aggregation of TFE molecules around the peptides. This coating displaces water, thereby removing alternative hydrogen-bonding partners and providing a low dielectric environment that favors the formation of intrapeptide hydrogen bonds. Because TFE interacts only weakly with nonpolar residues, hydrophobic interactions within the peptides are not disrupted. As a consequence, TFE promotes stability rather than inducing denaturation.