Molecular dynamics simulations of individual alpha-helices of bacteriorhodopsin in dimyristoylphosphatidylcholine. II. Interaction energy analysis.
Ontology highlight
ABSTRACT:
SUBMITTER: Woolf TB
PROVIDER: S-EPMC1299368 | biostudies-other | 1998 Jan
REPOSITORIES: biostudies-other
Publications
Molecular dynamics simulations of individual alpha-helices of bacteriorhodopsin in dimyristoylphosphatidylcholine. II. Interaction energy analysis.
Biophysical journal 19980101 1
The concepts of hydrophobicity and hydrophobic moments have been applied in attempts to predict membrane protein secondary and tertiary structure. The current paper uses molecular dynamics computer calculations of individual bacteriorhodopsin helices in explicit dimyristoylphosphatidylcholine bilayers to examine the atomic basis of these approaches. The results suggest that the types of interactions between a particular amino acid and the surrounding bilayer depend on the position and type of th ...[more]
Similar Datasets
Project description:BackgroundAnalysis of known protein structures reveals that identical sequence fragments in proteins can adopt different secondary structure conformations. The extent of this conformational diversity is influenced by various factors like the intrinsic sequence propensity, sequence context and other environmental factors such as pH, site directed mutations or alteration of the binding ligands. Understanding the mechanism by which the environment affects the structural ambivalence of these peptides has potential implications for protein design and reliable local structure prediction algorithms. Identification of the structurally ambivalent sequence fragments and determining the rules which dictate their conformational preferences play an important role in understanding the conformational changes observed in misfolding diseases. However, a systematic classification of their intrinsic sequence patterns or a statistical analysis of their properties and sequence context in relation to the origin of their structural diversity have largely remained unexplored.ResultsIn this work, the conformational variability of α-helices is studied by mapping sequences from the non-redundant database to identical sequences across all classes of the SCOP (Structural Classification of Proteins) database. Some helices retain their conformations when mapped in the SCOP database while others exhibit a complete/partial switch to non-helical conformations. The results clearly depict the differences in the propensities of amino acids for the variable and conserved helices. Sequences flanking these ambivalent sequence fragments have anisotropic propensities at the N- and C-termini. This structural variability is depicted by molecular dynamics simulations in explicit solvent, which show that the short conserved helices retain their conformations while their longer counterparts fray into two or more shorter helices. Variable helices in the non-redundant database exhibit a trend of retaining helical conformations while their corresponding non-helical conformations in SCOP database show large deviations from their respective initial structures by adopting partial or full helical conformations. Partially ambivalent helices are also found to retain their respective conformations.ConclusionsAll sequence fragments which show structural diversity in different proteins of the non-redundant database are investigated. The final conformation of these ambivalent sequences are dictated by a fine tuning of their intrinsic sequence propensity and the anisotropic amino acid propensity of the flanking sequences. This analysis may unravel the connection between diverse secondary structures, which conserve the overall structural fold of the protein thus determining its function.
| S-EPMC2973962 | biostudies-literature
Project description:The interactions of transmembrane (TM) α-helices with the phospholipid membrane and with one another are central to understanding the structure and stability of integral membrane proteins. These interactions may be analyzed via coarse grained molecular dynamics (CGMD) simulations. To obtain statistically meaningful analysis of TM helix interactions, large (N ca. 100) ensembles of CGMD simulations are needed. To facilitate the running and analysis of such ensembles of simulations, we have developed Sidekick, an automated pipeline software for performing high throughput CGMD simulations of α-helical peptides in lipid bilayer membranes. Through an end-to-end approach, which takes as input a helix sequence and outputs analytical metrics derived from CGMD simulations, we are able to predict the orientation and likelihood of insertion into a lipid bilayer of a given helix of a family of helix sequences. We illustrate this software via analyses of insertion into a membrane of short hydrophobic TM helices containing a single cationic arginine residue positioned at different positions along the length of the helix. From analyses of these ensembles of simulations, we estimate apparent energy barriers to insertion which are comparable to experimentally determined values. In a second application, we use CGMD simulations to examine the self-assembly of dimers of TM helices from the ErbB1 receptor tyrosine kinase and analyze the numbers of simulation repeats necessary to obtain convergence of simple descriptors of the mode of packing of the two helices within a dimer. Our approach offers a proof-of-principle platform for the further employment of automation in large ensemble CGMD simulations of membrane proteins.
| S-EPMC4871227 | biostudies-literature
Project description:Extensive all-atom molecular dynamics (?24 ?s total) allowed exploration of configurational space and calculation of lateral diffusion coefficients of the components of a protein-embedded, cholesterol-containing model bilayer. The three model membranes are composed of an ?50/50 (by mole) dipalmitoylphosphatidylcholine (DPPC)/cholesterol bilayer and contained an ?-helical transmembrane protein (HIV-1 gp41 TM). Despite the high concentration of cholesterol, normal Brownian motion was observed and the calculated diffusion coefficients (on the order of 10(-9) cm(2)/s) are consistent with experiments. Diffusion is sensitive to a variety of parameters, and a temperature difference of ?4 K from thermostat artifacts resulted in 2-10-fold differences in diffusion coefficients and significant differences in lipid order, membrane thickness, and unit cell area. Also, the specific peptide sequence likely underlies the consistently observed faster diffusion in one leaflet. Although the simulations here present molecular dynamics (MD) an order of magnitude longer than those from previous studies, the three systems did not approach ergodicity. The distributions of cholesterol and DPPC around the peptides changed on the microsecond time scale, but not significantly enough to thoroughly explore configurational space. These simulations support conclusions of other recent microsecond MD in that even longer time scales are needed for equilibration of model membranes and simulations of more realistic cellular or viral bilayers.
| S-EPMC4254001 | biostudies-literature
Project description:Multiple molecular dynamics simulations of bacterioopsin pulling from its C-terminus show that its alpha-helices unfold individually. In the first metastable state observed in the simulations, helix G is unfolded at its C-terminal segment while the rest of helix G (residues 200-216) is folded and opposes resistance because of a salt-bridge network consisting of Asp-212 and Lys-216 on helix G and Arg-82 and Asp-85 on helix C. Helix G unfolds inside the bundle because the external force is applied to its C-terminal end in a direction perpendicular to the surface of the membrane. Inversely, helix F has to flip by 180 degrees to exit from the membrane because the applied force and the helical N-C axis point in opposite directions. At the highest peak of the force, which cannot be interpreted in single-molecule force spectroscopy experiments, helix F has a pronounced kink at Pro-186. Mutation of Pro-186 and/or the charged side chains mentioned above, which are involved in very favorable electrostatic interactions in the low-dielectric region of the membrane, are expected to reduce the highest peak of the force. Helices E and D unfold in a similar way to helices G and F, respectively. Hence, the force-distance profile and sequence of events during forced unfolding of bacterioopsin are influenced by the up-and-down topology of the seven-helix bundle. The sequential extraction of individual helices from the membrane suggests that the spontaneous (un)folding of bacterioopsin proceeds through metastable bundles of fewer than seven helices. The metastable states observed in the simulations provide atomic level evidence that corroborates the interpretation of very recent force spectroscopy experiments of bacteriorhodopsin refolding.
| S-EPMC1614499 | biostudies-literature
Project description:The significant work that has been invested toward understanding protein-protein interaction has not translated into significant advances in structure-based predictions. In particular redesigning protein surfaces to bind to unrelated receptors remains a challenge, partly due to receptor flexibility, which is often neglected in these efforts. In this work, we computationally graft the binding epitope of various small proteins obtained from the RCSB database to bind to barnase, lysozyme, and trypsin using a previously derived and validated algorithm. In an effort to probe the protein complexes in a realistic environment, all native and designer complexes were subjected to a total of nearly 400 ns of explicit-solvent molecular dynamics (MD) simulation. The MD data led to an unexpected observation: some of the designer complexes were highly unstable and decomposed during the trajectories. In contrast, the native and a number of designer complexes remained consistently stable. The unstable conformers provided us with a unique opportunity to define the structural and energetic factors that lead to unproductive protein-protein complexes. To that end we used free energy calculations following the MM-PBSA approach to determine the role of nonpolar effects, electrostatics and entropy in binding. Remarkably, we found that a majority of unstable complexes exhibited more favorable electrostatics than native or stable designer complexes, suggesting that favorable electrostatic interactions are not prerequisite for complex formation between proteins. However, nonpolar effects remained consistently more favorable in native and stable designer complexes reinforcing the importance of hydrophobic effects in protein-protein binding. While entropy systematically opposed binding in all cases, there was no observed trend in the entropy difference between native and designer complexes. A series of alanine scanning mutations of hot-spot residues at the interface of native and designer complexes showed less than optimal contacts of hot-spot residues with their surroundings in the unstable conformers, resulting in more favorable entropy for these complexes. Finally, disorder predictions revealed that secondary structures at the interface of unstable complexes exhibited greater disorder than the stable complexes.
| S-EPMC2754190 | biostudies-literature
Project description:Gating of mechanosensitive (MS) channels is driven by a hierarchical cascade of movements and deformations of transmembrane helices in response to bilayer tension. Determining the intrinsic mechanical properties of the individual transmembrane helices is therefore central to understanding the intricacies of the gating mechanism of MS channels. We used a constant-force steered molecular dynamics (SMD) approach to perform unidirectional pulling tests on all the helices of MscL in M. tuberculosis and E. coli homologs. Using this method, we could overcome the issues encountered with the commonly used constant-velocity SMD simulations, such as low mechanical stability of the helix during stretching and high dependency of the elastic properties on the pulling rate. We estimated Young's moduli of the α-helices of MscL to vary between 0.2 and 12.5 GPa with TM2 helix being the stiffest. We also studied the effect of water on the properties of the pore-lining TM1 helix. In the absence of water, this helix exhibited a much stiffer response. By monitoring the number of hydrogen bonds, it appears that water acts like a 'lubricant' (softener) during TM1 helix elongation. These data shed light on another physical aspect underlying hydrophobic gating of MS channels, in particular MscL.
| S-EPMC5463901 | biostudies-literature
Project description:The influence of the solvent on the main-chain conformation (phi and Psi dihedral angles) of alpha-helices has been studied by complementary approaches. A first approach consisted in surveying crystal structures of both soluble and membrane proteins. The residues of analysis were further classified as exposed to either the water (polar solvent) or the lipid (apolar solvent) environment or buried to the core of the protein (intermediate polarity). The statistical results show that the more polar the environment, the lower the value of phi(i) and the higher the value of Psi(i) are. The intrahelical hydrogen bond distance increases in water-exposed residues due to the additional hydrogen bond between the peptide carbonyl oxygen and the aqueous environment. A second approach involved nanosecond molecular dynamics simulations of poly-Ala alpha-helices in environments of different polarity: water to mimic hydrophilic environments that can form hydrogen bonds with the peptide carbonyl oxygen and methane to mimic hydrophobic environments without this hydrogen bond capabilities. These simulations reproduce similar effects in phi and Psi angles and intrahelical hydrogen bond distance and angle as observed in the protein survey analysis. The magnitude of the intrahelical hydrogen bond in the methane environment is stronger than in the water environment, suggesting that alpha-helices in membrane-embedded proteins are less flexible than in soluble proteins. There is a remarkable coincidence between the phi and Psi angles obtained in the analysis of residues exposed to the lipid in membrane proteins and the results from computer simulations in methane, which suggests that this simulation protocol properly mimic the lipidic cell membrane and reproduce several structural characteristics of membrane-embedded proteins. Finally, we have compared the phi and Psi torsional angles of Pro kinks in membrane protein crystal structures and in computer simulations.
| S-EPMC1302110 | biostudies-other
Project description:Protein folding occurs as a set of transitions between structural states within an energy landscape. An oversimplified view of the folding process emerges when transiently populated states are undetected because of limited instrumental resolution. Using force spectroscopy optimized for 1-microsecond resolution, we reexamined the unfolding of individual bacteriorhodopsin molecules in native lipid bilayers. The experimental data reveal the unfolding pathway in unprecedented detail. Numerous newly detected intermediates-many separated by as few as two or three amino acids-exhibited complex dynamics, including frequent refolding and state occupancies of <10 μs. Equilibrium measurements between such states enabled the folding free-energy landscape to be deduced. These results sharpen the picture of the mechanical unfolding of membrane proteins and, more broadly, enable experimental access to previously obscured protein dynamics.
| S-EPMC5436802 | biostudies-literature
Project description:We describe application of the implicit solvation model (see the first paper of this series), to Monte Carlo simulations of several peptides in bilayer- and water-mimetic environments, and in vacuum. The membrane-bound peptides chosen were transmembrane segments A and B of bacteriorhodopsin, the hydrophobic segment of surfactant lipoprotein, and magainin2. Their conformations in membrane-like media are known from the experiments. Also, molecular dynamics study of surfactant lipoprotein with different explicit solvents has been reported (Kovacs, H., A. E. Mark, J. Johansson, and W. F. van Gunsteren. 1995. J. Mol. Biol. 247:808-822). The principal goal of this work is to compare the results obtained in the framework of our solvation model with available experimental and computational data. The findings could be summarized as follows: 1) structural and energetic properties of studied molecules strongly depend on the solvent; membrane-mimetic media significantly promote formation of alpha-helices capable of traversing the bilayer, whereas a polar environment destabilizes alpha-helical conformation via reduction of solvent-exposed surface area and packing; 2) the structures calculated in a membrane-like environment agree with the experimental ones; 3) noticeable differences in conformation of surfactant lipoprotein assessed via Monte Carlo simulation with implicit solvent (this work) and molecular dynamics in explicit solvent were observed; 4) in vacuo simulations do not correctly reproduce protein-membrane interactions, and hence should be avoided in modeling membrane proteins.
| S-EPMC1300218 | biostudies-other
Project description:The energetics and hydrogen bonding profiles of the helix-to-coil transition were found to be an additive property and to increase linearly with chain length, respectively, in alanine-rich α-helical peptides. A model system of polyalanine repeats was used to establish this hypothesis for the energetic trends and hydrogen bonding profiles. Numerical measurements of a synthesized polypeptide Ac-Y(AEAAKA)kF-NH2 and a natural α-helical peptide a2N (1-17) provide evidence of the hypothesis's generality. Adaptive steered molecular dynamics was employed to investigate the mechanical unfolding of all of these alanine-rich polypeptides. We found that the helix-to-coil transition is primarily dependent on the breaking of the intramolecular backbone hydrogen bonds and independent of specific side-chain interactions and chain length. The mechanical unfolding of the α-helical peptides results in a turnover mechanism in which a 310-helical structure forms during the unfolding, remaining at a near constant population and thereby maintaining additivity in the free energy. The intermediate partially unfolded structures exhibited polyproline II helical structure as previously seen by others. In summary, we found that the average force required to pull alanine-rich α-helical peptides in between the endpoints-namely the native structure and free coil-is nearly independent of the length or the specific primary structure.
| S-EPMC8204395 | biostudies-literature