Lessons from computer simulations of Ras proteins in solution and in membrane

Background

A great deal has been learned over the last several decades about the function of Ras proteins in solution and membrane environments. While much of this knowledge has been derived from a plethora of experimental techniques, computer simulations have also played a substantial role.

Scope of review

Our goal here is to summarize the contribution of molecular simulations to our current understanding of normal and aberrant Ras function. We focus on lessons from molecular dynamics simulations in aqueous and membrane environments.

Major conclusions

The central message is that a close interaction between theory and simulation on the one hand and cell-biological, spectroscopic and other experimental approaches on the other has played, and will likely continue to play, a vital role in Ras research.

General significance

Atomistic insights emerging from detailed simulations of Ras in solution and in bilayers may be the key to unlock the secret that to date prevented development of selective anti-Ras inhibitors for cancer therapy.     

Structure and dynamics of interfacial water in an Lalpha phase lipid bilayer from molecular dynamics simulations

Based on molecular dynamics simulations, an analysis of structure and dynamics is performed on interfacial water at a liquid crystalline dipalmitoylphosphatidycholine/water system. Water properties relevant for understanding NMR relaxation are emphasized. The first and second rank orientational order parameters of the water O-H bonds were calculated, where the second rank order parameter is in agreement with experimental determined quadrupolar splittings. Also, two different interfacial water regions (bound water regions) are revealed with respect to different signs of the second rank order parameter. The water reorientation correlation function reveals a mixture of fast and slow decaying parts. The fast (ps) part of the correlation function is due to local anisotropic water reorientation whereas the much slower part is due to more complicated processes including lateral diffusion along the interface and chemical exchange between free and bound water molecules. The 100-ns-long molecular dynamics simulation at constant pressure (1 atm) and at a temperature of 50 degrees C of 64 lipid molecules and 64 x 23 water molecules lack a slow water reorientation correlation component in the ns time scale. The (2)H(2)O powder spectrum of the dipalmitoylphosphatidycholine/water system is narrow and consequently, the NMR relaxation time T(2) is too short compared to experimental results.     

Structure and dynamics of Antarctic fish neuroglobin assessed by computer simulations

Neuroglobin (Ngb) is a heme protein, highly conserved along evolution, predominantly found in the nervous system. It is upregulated by hypoxia and ischemia and may have a neuroprotective role under hypoxic stress. Although many other roles have been proposed, the physiological function is still unclear. Antarctic icefishes lack hemoglobin and some species also lack myoglobin, but all have Ngb and thus may help the elucidation of Ngb function. We present the first theoretically derived structure of fish Ngb and describe its behavior using molecular dynamics simulations. Specifically, we sequenced and analyzed Ngbs from a colorless-blooded Antarctic icefish species Chaenocephalus aceratus and a related red-blooded species (Dissostichus mawsoni). Both fish Ngbs are 6-coordinated but have some peculiarities that differentiate them from mammalian counterparts: they have extensions in the N and C termini that can interact with the EF loop, and a gap in the alignment that changes the CD-region structure/dynamics that has been found to play a key role in human neuroglobin. Our results suggest that a single mutation between both fish Ngbs is responsible for significant difference in the behavior of the proteins. The functional role of these characteristics is discussed.     

An alamethicin channel in a lipid bilayer: molecular dynamics simulations

We present the results of 2-ns molecular dynamics (MD) simulations of a hexameric bundle of Alm helices in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer. These simulations explore the dynamic properties of a model of a helix bundle channel in a complete phospholipid bilayer in an aqueous environment. We explore the stability and conformational dynamics of the bundle in a phospholipid bilayer. We also investigate the effect on bundle stability of the ionization state of the ring of Glu18 side chains. If all of the Glu18 side chains are ionised, the bundle is unstable; if none of the Glu18 side chains are ionized, the bundle is stable. pKA calculations suggest that either zero or one ionized Glu18 is present at neutral pH, correlating with the stable form of the helix bundle. The structural and dynamic properties of water in this model channel were examined. As in earlier in vacuo simulations (Breed et al., 1996 .Biophys. J. 70:1643-1661), the dipole moments of water molecules within the pore were aligned antiparallel to the helix dipoles. This contributes to the stability of the helix bundle.     

Molecular dynamics simulations of the lipid bilayer edge

Phospholipid bilayers have been intensively studied by molecular dynamics (MD) simulation in recent years. The properties of bilayer edges are important in determining the structure and stability of pores formed in vesicles and biomembranes. In this work, we use molecular dynamics simulation to investigate the structure, dynamics, and line tension of the edges of bilayer ribbons composed of pure dimyristoylphosphatidylcholine (DMPC) or palmitoyl-oleoylphosphatidylethanolamine (POPE). As expected, we observe a significant reorganization of lipids at and near the edges. The treatment of electrostatic effects is shown to have a qualitative impact on the structure and stability of the edge, and significant differences are observed in the dynamics and structure of edges formed by DMPC and palmitoyl-oleoylphosphatidylethanolamine. From the pressure anisotropy in the simulation box, we calculate a line tension of approximately 10-30 pN for the DMPC edge, in qualitative agreement with experimental estimates for similar lipids.     

Distribution of amino acids in a lipid bilayer from computer simulations

We have calculated the distribution in a lipid bilayer of small molecules mimicking 17 natural amino acids in atomistic detail by molecular dynamics simulation. We considered both charged and uncharged forms for Lys, Arg, Glu, and Asp. The results give detailed insight in the molecular basis of the preferred location and orientation of each side chain as well the preferred charge state for ionizable residues. Partitioning of charged and polar side chains is accompanied by water defects connecting the side chains to bulk water. These water defects dominate the energetic of partitioning, rather than simple partitioning between water and a hydrophobic phase. Lys, Glu, and Asp become uncharged well before reaching the center of the membrane, but Arg may be either charged or uncharged at the center of the membrane. Phe has a broad distribution in the membrane but Trp and Tyr localize strongly to the interfacial region. The distributions are useful for the development of coarse-grained and implicit membrane potentials for simulation and structure prediction. We discuss the relationship between the distribution in membranes, bulk partitioning to cyclohexane, and several amino acid hydrophobicity scales.     

Structure of crambin in solution, crystal and in the trajectories of molecular dynamics simulations

The mechanisms of the three-dimensional crambin structure alterations in the crystalline environments and in the trajectories of the molecular dynamics simulations in the vacuum and crystal surroundings have been analyzed. In the crystalline state and in the solution the partial regrouping of remote intramolecular packing contacts, involved in the formation and stabilization of the tertiary structure of the crambin molecule, occurs in NMR structures. In the crystalline state it is initiated by the formation of the intermolecular contacts, the conformational influence of its appearance is distributed over the structure. The changes of the conformations and positions of the residues of the loop segments, where the intermolecular contacts of the crystal surroundings are preferably concentrated, are most observable. Under the influence of these contacts the principal change of the regular secondary structure of crambin is taking place: extension of the two-strand β structure to the three-strand structure with the participation of the single last residue N46 of the C-terminal loop. In comparison with the C-terminal loop the more profound changes are observed in the conformation and the atomic positions of the backbone atoms and in the solvent accessibility of the residues of the interhelical loop. In the solution of the ensemble of the 8 NMR structures relative accessibility to the solvent differs more noticeably also in the region of the loop segments and rather markedly in the interhelical loop. In the crambin cryogenic crystal structures the positions of the atoms of the backbone and/or side chain of 14–18 of 46 residues are discretely disordered. The disorganizations of at least 8 of 14 residues occur directly in the regions of the intermolecular contacts and another 5 residues are disordered indirectly through the intramolecular contacts with the residues of the intermolecular contacts. Upon the molecular dynamics simulation in the vacuum surrounding as in the solution of the crystalline structure of crambin the essential changes of the backbone conformation are caused by the intermolecular contacts absence, but partly masked by the structure changes owing to the nonpolar H atoms absence on the simulated structure. The intermolecular contact absence is partly manifested upon the molecular dynamics simulation of the crambin crystal with one protein molecule. Compared to the crystal structure the lengths of the interpeptide hydrogen bonds and other interresidue contacts in an average solution NMR structure are somewhat shorter and accordingly the energy of the interpeptide hydrogen bonds is better. This length shortening can occur at the stage of the refinement of the NMR structures of the crambin and other proteins by its energy minimizations in the vacuum surroundings and not exist in the solution protein structures.     

Interfacing molecular dynamics and macro-scale simulations for lipid bilayer vesicles

A continuum-level model for a giant unilamellar vesicle (GUV) is bridged to a corresponding atomistic model of a dimyristoylphosphatidylcholine (DMPC) bilayer at various cholesterol concentrations via computation of the bulk modulus. The bulk modulus and other microscopically determined parameters are passed to a continuum-level model operating in time- and length-scales orders of magnitude beyond that which is accessible by atomistic-level simulation. The continuum-level simulation method used is the material point method (MPM), and the particular variation used here takes advantage of the spherical nature of many GUVs. An osmotic pressure gradient due to a solvent concentration change is incorporated into the continuum-level simulation, resulting in osmotic swelling of the vesicle. The model is then extended to treat mixtures of DMPC and cholesterol, where small domains of different composition are considered.     

The transmembrane domain of Neu in a lipid bilayer: molecular dynamics simulations

The results of full-atom molecular dynamics simulations of the transmembrane domains (TMDs) of both native, and Glu⁶⁶⁴-mutant (either protonated or unprotonated) Neu in an explicit fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer are presented. For the native TMD peptide, a 10.05 ns trajectory was collected, while for the mutant TMD peptides 5.05 ns trajectories were collected for each. The peptides in all three simulations display stable predominantly -helical hydrogen bonding throughout the trajectories. The only significant exception occurs near the C-terminal end of the native and unprotonated mutant TMDs just outside the level of the lipid headgroups, where -helical hydrogen bonding develops, introducing a kink in the backbone structure. However, there is no indication of the formation of a bulge within the hydrophobic region of either native or mutant peptides. Over the course of the simulation of the mutant peptide, it is found that a significant number of water molecules penetrate the hydrophobic region of the surrounding lipid molecules, effectively hydrating Glu⁶⁶⁴. If the energy cost of such water penetration is significant enough, this may be a factor in the enhanced dimerization affinity of Glu⁶⁶⁴-mutant Neu.     

Molecular dynamics simulations of membrane proteins

Membrane proteins control the traffic across cell membranes and thereby play an essential role in cell function from transport of various solutes to immune response via molecular recognition. Because it is very difficult to determine the structures of membrane proteins experimentally, computational methods have been increasingly used to study their structure and function. Here we focus on two classes of membrane proteins—ion channels and transporters—which are responsible for the generation of action potentials in nerves, muscles, and other excitable cells. We describe how computational methods have been used to construct models for these proteins and to study the transport mechanism. The main computational tool is the molecular dynamics (MD) simulation, which can be used for everything from refinement of protein structures to free energy calculations of transport processes. We illustrate with specific examples from gramicidin and potassium channels and aspartate transporters how the function of these membrane proteins can be investigated using MD simulations.     

Structure and internal dynamics of the bovine pancreatic trypsin inhibitor in aqueous solution from long-time molecular dynamics simulations

Structural and dynamic properties of bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution are investigated using two molecular dynamics (MD) simulations: one of 1.4 ns length and one of 0.8 ns length in which atom-atom distance bounds derived from NMR spectroscopy are included in the potential energy function to make the trajectory satisfy these experimental data more closely. The simulated properties of BPTI are compared with crystal and solution structures of BPTI, and found to be in agreement with the available experimental data. The best agreement with experiment was obtained when atom-atom distance restraints were applied in a time-averaged manner in the simulation. The polypeptide segments found to be most flexible in the MD simulations coincide closely with those showing differences between the crystal and solution structures of BPTI. © 1995 Wiley-Liss, Inc.     

Molecular dynamics simulations of local anesthetic articaine in a lipid bilayer

In order to investigate structural and dynamical properties of local anesthetic articaine in a model lipid bilayer, a series of molecular dynamics simulations have been performed. Simulations were carried out for neutral and charged (protonated) forms of articaine inserted in fully hydrated dimyristoylphosphatidylcholine (DMPC) lipid bilayer. For comparison purpose, a fully hydrated DMPC bilayer without articaine was also simulated. The length of each simulation was 200 ns. Various properties of the lipid bilayer systems in the presence of both charged and uncharged forms of articaine taken at two different concentrations have been examined: membrane area per lipid, mass density distributions, order parameters, radial distribution functions, head group tilt, diffusion coefficients, electrostatic potential, etc, and compared with results of previous simulations of DMPC bilayer in the presence of lidocaine. It was shown that addition of both charged and neutral forms of articaine causes increase of the dipole electrostatic potential in the membrane interior.     

Computer simulations of membrane protein folding: structure and dynamics

A lattice model of membrane proteins with a composite energy function is proposed to study their folding dynamics and native structures using Monte Carlo simulations. This model successfully predicts the seven helix bundle structure of sensory rhodopsin I by practicing a three-stage folding. Folding dynamics of a transmembrane segment into a helix is further investigated by varying the cooperativity in the formation of alpha helices for both random folding and assisted folding. The chain length dependence of the folding time of a hydrophobic segment to a helical state is studied for both free and anchored chains. An unusual length dependence in the folding time of anchored chains is observed.     

Molecular dynamics simulations of lignin peroxidase in solution

The dynamical and structural properties of lignin peroxidase and its Trp171Ala mutant have been investigated in aqueous solution using molecular dynamics (MD) simulations. In both cases, the enzyme retained its overall backbone structure and all its noncovalent interactions in the course of the MD simulations. Very interestingly, the analysis of the MD trajectories showed the presence of large fluctuations in correspondence of the residues forming the heme access channel; these movements enlarge the opening and facilitate the access of substrates to the enzyme active site. Moreover, steered molecular dynamics docking simulations have shown that lignin peroxidase natural substrate (veratryl alcohol) can easily approach the heme edge through the access channel.     

Molecular dynamics simulations of calcium-free calmodulin in solution

A 4-ns molecular dynamics simulation of calcium-free calmodulin in solution has been performed, using Ewald summation to treat electrostatic interactions. Our simulation results were mostly consistent with solution experimental studies, including NMR, fluorescence and x-ray scattering. The secondary structures within the N- and C-terminal domains were conserved in the simulation, with trajectory structures similar to the NMR-derived model structure 1CFD. However, the relative orientations of the domains, for which there are no NMR restraints, differed in details between the simulation and the 1CFD model. The most interesting information provided by the simulations is that the dynamics of calcium-free calmodulin in solution is dominated by slow rigid body reorientations of the domains. The interdomain distance fluctuated between 29 and 39 A, and interdomain orientation angle, defined as the pseudo-dihedral formed by the four calcium binding sites, varied between -2 degrees and 108 degrees. Similarly, the domain linker region also exhibited significant fluctuations, with its length varying in the 34-45 A range and its bend angle in the 10-100 degrees range. The simulations are in accord with fluorescence results suggesting that calcium-free calmodulin is more compact and more flexible than the calcium activated form. Surprisingly, quite similar solvent accessibilities of the hydrophobic patches were seen in the calcium-free trajectory described in this work and previously generated calcium-loaded calmodulin simulations. Thus, our simulations suggest a reexamination of the standard model of the structural change of calmodulin upon calcium binding, involving exposure of the hydrophobic patches to solvent.     

Behaviour of small solutes and large drugs in a lipid bilayer from computer simulations

To reach their biological target, drugs have to cross cell membranes, and understanding passive membrane permeation is therefore crucial for rational drug design. Molecular dynamics simulations offer a powerful way of studying permeation at the single molecule level. Starting from a computer model proven to be able to reproduce the physical properties of a biological membrane, the behaviour of small solutes and large drugs in a lipid bilayer has been studied. Analysis of dihedral angles shows that a few nano seconds are sufficient for the simulations to converge towards common values for those angles, even if the starting structures belong to different conformations. Results clearly show that, despite their difference in size, small solutes and large drugs tend to lie parallel to the bilayer normal and that, when moving from water solution into biomembranes, permeants lose degrees of freedom. This explains the experimental observation that partitioning and permeation are highly affected by entropic effects and are size-dependent. Tilted orientations, however, occur when they make possible the formation of hydrogen bonds. This helps to understand the reason why hydrogen bonding possibilities are an important parameter in cruder approaches which predict drug absorption after administration. Interestingly, hydration is found to occur even in the membrane core, which is usually considered an almost hydrophobic region. Simulations suggest the possibility for highly polar compounds like acetic acid to cross biological membranes while hydrated. These simulations prove useful for drug design in rationalising experimental observations and predicting solute behaviour in biomembranes.     

Structure and dynamics of an amphiphilic peptide in a lipid bilayer: a molecular dynamics study

A molecular dynamics simulation of a simple model membrane system composed of a single amphiphilic helical peptide (ace-K2GL16K2A-amide) in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayer was performed for a total of 1060 ps. The secondary structure of the peptide and its stability were described in terms of average dihedral angles, phi and psi, and the C alpha torsion angles formed by backbone atoms; by the average translation per residue along the helix axis; and by the intramolecular peptide hydrogen bonds. The results indicated that residues 6 through 15 remain in a stable right-handed alpha-helical conformation, whereas both termini exhibit substantial fluctuations. A change in the backbone dihedral angles for residues 16 and 17 is accompanied by the loss of two intramolecular hydrogen bonds, leading to a local but long-lived disruption of the helix. The dynamics of the peptide was characterized in terms of local and global helix motions. The local motions of the N-H bond angles were described in terms of the autocorrelation functions of P2[cos thetaNH(t, t + tau)] and reflected the different degrees of local peptide order as well as a variation in time scale for local motions. The chi1 and chi2 dihedral angles of the leucine side chains underwent frequent transitions between potential minima. No connection between the side-chain positions and their mobility was observed, however. In contrast, the lysine side chains displayed little mobility during the simulation. The global peptide motions were characterized by the tilting and bending motions of the helix. Although the peptide was initially aligned parallel to the bilayer normal, during the simulation it was observed to tilt away from the normal, reaching an angle of approximately 25 degrees by the end of the simulation. In addition, a slight bend of the helix was detected. Finally, the solvation of the peptide backbone and side-chain atoms was also investigated.     

Structure and dynamics of the dilauroylphosphatidylethanolamine lipid bilayer.

A 200-ps molecular dynamics (MD) simulation trajectory of a model dilauroylphosphatidylethanolamine (DLPE) bilayer in water at 315 K has been generated. Segmental order parameters, electron density profiles, and water pair distribution functions have been calculated. Comparison to experiment is made where possible. The dynamics of the system has been studied by analyzing the velocity autocorrelation functions (VAF) of both water and lipid atoms. Furthermore, the diffusive properties of water have been analyzed by computing the mean square displacement (MSD) and orientational correlation function (OCF) of water in two regions around the bilayer. The calculated order parameters show a behavior similar to the liquid crystalline phase of other bilayers, but the region around C1-C3 does not show the expected behavior. The electron density profile shows features that are characteristic of the liquid crystalline phase. The radial distribution functions suggest ordering of water near the charged head groups, which results in about 15 water molecules solvating each lipid molecule. We find from the VAF, MSD, and OCF calculation that the water molecules near the head groups of the lipid bilayer move more slowly than those further away. The VAF of the hydrocarbon chains have features of low-frequency motions that are probably cooperative nature in addition to the high-frequency motions associated with bond angle and torsional motions.