Membrane insertion of a lipidated ras peptide studied by FTIR,solid-state NMR,and neutron diffraction spectroscopy

Membrane binding of a doubly lipid modified heptapeptide from the C-terminus of the human N-ras protein was studied by Fourier transform infrared, solid-state NMR, and neutron diffraction spectroscopy. The 16:0 peptide chains insert well into the 1,2-dimyristoyl-sn-glycero-3-phosphocholine phospholipid matrix. This is indicated by a common main phase transition temperature of 21.5 degrees C for both the lipid and peptide chains as revealed by FTIR measurements. Further, (2)H NMR reveals that peptide and lipid chains have approximately the same chain length in the liquid crystalline state. This is achieved by a much lower order parameter of the 16:0 peptide chains compared to the 14:0 phospholipid chains. Finally, proton/deuterium contrast variation of neutron diffraction experiments indicates that peptide chains are localized in the membrane interior analogous to the phospholipid chains. In agreement with this model of peptide chain insertion, the peptide part is localized at the lipid-water interface of the membrane. This is revealed by (1)H nuclear Overhauser enhancement spectra recorded under magic angle spinning conditions. Quantitative cross-peak analysis allows the examination of the average location of the peptide backbone and side chains with respect to the membrane. While the backbone shows the strongest cross-relaxation rates with the phospholipid glycerol, the hydrophobic side chains of the peptide insert deeper into the membrane interior. This is supported by neutron diffraction experiments that reveal a peptide distribution in the lipid-water interface of the membrane. Concurring with these experimental findings, the amide protons of the peptide show strong water exchange as seen in NMR and FTIR measurements. No indications for a hydrogen-bonded secondary structure of the peptide backbone are found. Therefore, membrane binding of the C-terminus of the N-ras protein is mainly due to lipid chain insertion but also supported by interactions between hydrophobic side chains and the lipid membrane. The peptide assumes a mobile and disordered conformation in the membrane. Since the C-terminus of the soluble part of the ras protein is also disordered, we hypothesize that our model for membrane binding of the ras peptide realistically describes the membrane binding of the lipidated C-terminus of the active ras protein.     

Backbone dynamics of proteins studied by two-dimensional heteronuclear NMR spectroscopy and molecular dynamics simulations

To investigate the backbone dynamics of proteins ¹⁵N longitudinal and transverse relaxation experiments combined with {¹H, ¹⁵N{ NOE measurements together with molecular dynamics simulations were carried out using ribonuclease T₁ and the complex of ribonuclease T₁ with 2′GMP as a model protein. The intensity decay of individual amide cross peaks in a series of (¹H, ¹⁵N)HSQC spectra with appropriate relaxation periods was fitted to a single exponential by using a simplex algorithm in order to obtain ¹⁵N T₁ and T₂ relaxation times. The relaxation times were analyzed in terms of the “model-free” approach introduced by Lipari and Szabo. In addition, a nanosecond molecular dynamics (MD ) simulation of ribonuclease T₁ and its 2′GMP complex in water was carried out. The angular reorientations of the backbone amide groups were classified with several coordinate frames following a transformation of NH vector trajectories. In this study, NH librations and backbone dihedral angle fluctuations were distinguished. The NH bond librations were found to be similar for all amides as characterized by correlation times of librational motions in a subpicosecond scale. The angular amplitudes of these motions were found to be about 10°–12° for out-of-plane displacements and 3°–5° for the in-plane displacement. The contributions from the much slower backbone dihedral angle fluctuations strongly depend on the secondary structure. The dependence of the amplitude of local motion on the residue location in the backbone is in good agreement with the results of NMR relaxation measurements and the X-ray data. The protein dynamics is characterized by a highly restricted local motion of those parts of the backbone with defined secondary structure as well as by a high flexibility in loop regions. Comparison of the MD and NMR data of the free liganded enzyme ribonuclease T₁ clearly indicates a restriction of the mobility within certain regions of the backbone upon inhibitor binding. © 1996 John Wiley & Sons, Inc.     

Molecular motions of alpha-cyclodextrin on a dodecyl chain studied by molecular dynamics simulations

Motions of an alpha-cyclodextrin (alpha-CD) molecule on a dodecyl chain adopting the all-trans conformation were investigated in the presence of water by molecular dynamics simulations with CVFF force fields, where the trimethylammonium group of dodecyltrimethylammonium bromide (DTAB) is protruded outside the secondary hydroxyl rim of alpha-CD (the secondary-in structure). The alpha-CD molecule shuttled rapidly on the chain without decomplexation. This rapid motion is consistent with the NMR data. The plane formed by 6 O4 atoms of alpha-CD is most populated between the C6 and C7 atoms of DTAB. This structure is very close to that estimated by NMR. The alpha-CD molecule underwent a restricted rotation in a range of 60 degrees with regard to the plane of the dodecyl chain: this plane at the most population is middle between the two diagonal lines of the normal hexagon formed by 6 O4 atoms of alpha-CD. The published NMR data were reanalyzed in terms of the rotation angle, and a slightly better structure was obtained. The distortion of the alpha-CD cavity from the normal hexagon was decreased upon complex formation with DTAB. The deviation of the center of alpha-CD from the center of the dodecyl chain predicted by molecular dynamics simulations is consistent with the NMR data. The secondary-in structure is energetically more stable than the primary-in structure, as calculated by molecular mechanics with CVFF and Amber force fields. This result is consistent with the NMR data. Molecular dynamics simulations were also carried out for the primary-in structure. Some of the results are close to those of the secondary-in structure.     

Intermolecular polarizability dynamics of aqueous formamide liquid mixtures studied by molecular dynamics simulations

A molecular dynamics simulation study is presented for the relaxation of the polarizability anisotropy in liquid mixtures of formamide and water, using a dipolar induction scheme that involves the intrinsic polarizability and first hyperpolarizability tensors of the molecules, and the dipole-quadrupole polarizability of water species. The long time diffusive decay of the collective polarizability anisotropy correlations exhibits a substantial slowing down as the formamide mole fraction increases in the mixture. The diffusive times for the polarizability relaxation obtained from the authors' simulations are in good agreement with optical Kerr effect experimental data, and they are found to correlate nearly linearly with the estimated mean lifetimes of the hydrogen bonds within the mixture, suggesting that the relaxation of the hydrogen bond network is responsible to some extent for the collective relaxation of the polarizability anisotropy of the mixture. The short time behavior of the polarizability anisotropy relaxation was investigated by computing the nuclear response function, R(t), which is very rapidly dominated by the formamide contribution as it is added to water, due to the much larger polarizability anisotropy of formamide molecules compared to that of water. Several contributions to the Raman spectrum were also analyzed as a function of composition, and the dynamical origin of the different bands was determined.     

Temperature and pressure dependence of alanine dipeptide studied by multibaric-multithermal molecular dynamics simulations

We applied the multibaric-multithermal (MUBATH) molecular dynamics (MD) algorithm to an alanine dipeptide in explicit water. The MUBATH MD simulation covered a wide range of conformational space and sampled the states of PII, C5, alphaR, alphaP, alphaL, and C7(ax). On the other hand, the conventional isobaric-isothermal simulation was trapped in local-minimum free-energy states and sampled only a few of them. We calculated the partial molar enthalpy difference DeltaH and partial molar volume difference DeltaV among these states by the MUBATH simulation using the AMBER parm99 and AMBER parm96 force fields and two sets of initial conditions. We compared these results with those from Raman spectroscopy experiments. The Raman spectroscopy data of DeltaH for the C5 state against the PII state agreed with both MUBATH data with the AMBER parm96 and parm99 force fields. The partial molar enthalpy difference DeltaH for the alphaR state and the partial molar volume difference DeltaV for the C5 state by the Raman spectroscopy agreed with those for the AMBER parm96 force field. On the other hand, DeltaV for the alphaR state by the Raman spectroscopy was consistent with our AMBER-parm99 force-field result. All the experimental results fall between those of simulations using AMBER parm96 and parm99 force fields, suggesting that the ideal force-field parameters lie between those of AMBER parm96 and parm99.     

Effect of acetone accumulation on structure and dynamics of lipid membranes studied by molecular dynamics simulations

The modulation of the properties and function of cell membranes by small volatile substances is important for many biomedical applications. Despite available experimental results, molecular mechanisms of action of inhalants and organic solvents, such as acetone, on lipid membranes remain not well understood. To gain a better understanding of how acetone interacts with membranes, we have performed a series of molecular dynamics (MD) simulations of a POPC bilayer in aqueous solution in the presence of acetone, whose concentration was varied from 2.8 to 11.2 mol%. The MD simulations of passive distribution of acetone between a bulk water phase and a lipid bilayer show that acetone favors partitioning into the water-free region of the bilayer, located near the carbonyl groups of the phospholipids and at the beginning of the hydrocarbon core of the lipid membrane. Using MD umbrella sampling, we found that the permeability barrier of ∼0.5 kcal/mol exists for acetone partitioning into the membrane. In addition, a Gibbs free energy profile of the acetone penetration across a bilayer demonstrates a favorable potential energy well of −3.6 kcal/mol, located at 15–16 Å from the bilayer center. The analysis of the structural and dynamics properties of the model membrane revealed that the POPC bilayer can tolerate the presence of acetone in the concentration range of 2.8–5.6 mol%. The accumulation of the higher acetone concentration of 11.2 mol% results, however, in drastic disordering of phospholipid packing and the increase in the membrane fluidity. The acetone molecules push the lipid heads apart and, hence, act as spacers in the headgroup region. This effect leads to the increase in the average headgroup area per molecule. In addition, the acyl tail region of the membrane also becomes less dense. We suggest, therefore, that the molecular mechanism of acetone action on the phospholipid bilayer has many common features with the effects of short chain alcohols, DMSO, and chloroform.     

Replica exchange molecular dynamics simulations of amyloid peptide aggregation

The replica exchange molecular dynamics (REMD) approach is applied to four oligomeric peptide systems. At physiologically relevant temperature values REMD samples conformation space and aggregation transitions more efficiently than constant temperature molecular dynamics (CTMD). During the aggregation process the energetic and structural properties are essentially the same in REMD and CTMD. A condensation stage toward disordered aggregates precedes the beta-sheet formation. Two order parameters, borrowed from anisotropic fluid analysis, are used to monitor the aggregation process. The order parameters do not depend on the peptide sequence and length and therefore allow to compare the amyloidogenic propensity of different peptides     

The plastic and liquid phases of CCl3Br studied by molecular dynamics simulations

We present a molecular dynamics study of the liquid and plastic crystalline phases of CCl(3)Br. We investigated the short-range orientational order using a recently developed classification method and we found that both phases behave in a very similar way. The only differences occur at very short molecular separations, which are shown to be very rare. The rotational dynamics was explored using time correlation functions of the molecular bonds. We found that the relaxation dynamics corresponds to an isotropic diffusive mode for the liquid phase but departs from this behavior as the temperature is decreased and the system transitions into the plastic phase.     

Drug-dendrimer supramolecular complexation studied from molecular dynamics simulations and NMR spectroscopy

Fully atomistic molecular dynamics (MD) simulations and NMR spectroscopy were employed to get insights about the molecular details of drug-dendrimer supramolecular association phenomena, using piroxicam (PRX) and the third generation poly(amido amine) (PAMAM-G3) dendrimer as model systems. Theoretical results concerning the complex stoichiometry suggest that PRX forms drug-dendrimer complexes of the type 24:1 at pH 7.0. This result was validated with the experimental quantities obtained from aqueous solubility profiles, which led to an empiric stoichiometry of 23:1 for the PRX:PAMAM-G3 system. The predicted binding mode between PRX and PAMAM-G3 accounts for the preferred encapsulation of the drug inside dendrimer cavities, which is mainly driven by van der Waals and hydrogen bonding interactions, and to a lesser extent, for the external association of the guest through electrostatic contacts with the positively charged amino groups of PAMAM periphery. The binding mode obtained from MD simulations was confirmed with 2D-NOESY experiments, which evidence the preferred internal complexation of PRX with PAMAM-G3. The predominance of internal encapsulation over external contacts in the PRX:PAMAM-G3 system differs from the general behaviour expected for acidic anionic guests, for which external electrostatic interactions with the positively charged PAMAM surface have been postulated as the most relevant factor for drug association.     

Water motion in reverse micelles studied by quasielastic neutron scattering and molecular dynamics simulations

Motion of water molecules in Aerosol OT [sodium bis(2-ethylhexyl) sulfosuccinate, AOT] reverse micelles with water content w(0) ranging from 1 to 5 has been explored both experimentally through quasielastic neutron scattering (QENS) and with molecular dynamics (MD) simulations. The experiments were performed at the energy resolution of 85 microeV over the momentum transfer (Q) range of 0.36-2.53 A(-1) on samples in which the nonpolar phase (isooctane) and the AOT alkyl chains were deuterated, thereby suppressing their contribution to the QENS signal. QENS results were analyzed via a jump-diffusion/isotropic rotation model, which fits the results reasonably well despite the fact that confinement effects are not explicitly taken into account. This analysis indicates that in reverse micelles with low-water content (w(0)=1 and 2.5) translational diffusion rate is too slow to be detected, while for w(0)=5 the diffusion coefficient is much smaller than for bulk water. Rotational diffusion coefficients obtained from this analysis increase with w(0) and are smaller than for bulk water, but rotational mobility is less drastically reduced than translational mobility. Using the Faeder/Ladanyi model [J. Phys. Chem. B 104, 1033 (2000)] of reverse micelle interior, MD simulations were performed to calculate the self-intermediate scattering function F(S)(Q,t) for water hydrogens. Comparison of the time Fourier transform of this F(S)(Q,t) with the QENS dynamic structure factor S(Q,omega), shows good agreement between the model and experiment. Separate intermediate scattering functions F(S) (R)(Q,t) and F(S) (CM)(Q,t) were determined for rotational and translational motion. Consistent with the decoupling approximation used in the analysis of QENS data, the product of F(S) (R)(Q,t) and F(S) (CM)(Q,t) is a good approximation to the total F(S)(Q,t). We find that the decay of F(S) (CM)(Q,t) is nonexponential and our analysis of the MD data indicates that this behavior is due to lower water mobility close to the interface and to confinement-induced restrictions on the range of translational displacements. Rotational relaxation also exhibits nonexponential decay. However, rotational mobility of O-H bond vectors in the interfacial region remains fairly high due to the lower density of water-water hydrogen bonds in the vicinity of the interface.     

Diffusion in binary gas mixtures studied by NMR of hyperpolarized gases and molecular dynamics simulations

The dependence of the individual mean square displacement of rare gases in binary mixtures is studied by a combined experimental and theoretical approach. We show that the diffusion constant can be varied in a considerable range by changing the molar fractions of the mixtures. On the experimental side, NMR diffusion measurements are done on hyperpolarized 3He and 129Xe, mixed with several inert buffer gases, in the presence of a magnetic field gradient. The results are compared to diffusion coefficients obtained from atomistic molecular dynamics simulations based on Lennard-Jones type potentials of the corresponding gas mixtures, and to appropriate analytical expressions, yielding very good mutual agreement. This study is the first quantitative validation of the effects of the mutual interactions between gas particles on the individual diffusion properties. It is shown that the dependency of gas phase diffusion properties on the local chemical environment may not be neglected, e.g. in diffusion-controlled chemical reactions.     

Interaction between water molecules and zinc sulfide nanoparticles studied by temperature-programmed desorption and molecular dynamics simulations

We have investigated the bonding of water molecules to the surfaces of ZnS nanoparticles (approximately 2-3 nm sphalerite) using temperature-programmed desorption (TPD). The activation energy for water desorption was derived as a function of the surface coverage through kinetic modeling of the experimental TPD curves. The binding energy of water equals the activation energy of desorption if it is assumed that the activation energy for adsorption is nearly zero. Molecular dynamics (MD) simulations of water adsorption on 3 and 5 nm sphalerite nanoparticles provided insights into the adsorption process and water binding at the atomic level. Water binds with the ZnS nanoparticle surface mainly via formation of Zn-O bonds. As compared with bulk ZnS crystals, ZnS nanoparticles can adsorb more water molecules per unit surface area due to the greatly increased curvature, which increases the distance between adjacent adsorbed molecules. Results from both TPD and MD show that the water binding energy increases with decreasing the water surface coverage. We attribute the increase in binding energy with decreasing surface water coverage to the increasing degree of surface under-coordination as removal of water molecules proceeds. MD also suggests that the water binding energy increases with decreasing particle size due to the further distance and hence lower interaction between adsorbed water molecules on highly curved smaller particle surfaces. Results also show that the binding energy, and thus the strength of interaction of water, is highest in isolated nanoparticles, lower in nanoparticle aggregates, and lowest in bulk crystals. Given that water binding is driven by surface energy reduction, we attribute the decreased binding energy for aggregated as compared to isolated particles to the decrease in surface energy that occurs as the result of inter-particle interactions.     

Ejection of solvated ions from electrosprayed methanol/water nanodroplets studied by molecular dynamics simulations

The ejection of solvated small ions from nanometer-sized droplets plays a central role during electrospray ionization (ESI). Molecular dynamics (MD) simulations can provide insights into the nanodroplet behavior. Earlier MD studies have largely focused on aqueous systems, whereas most practical ESI applications involve the use of organic cosolvents. We conduct simulations on mixed water/methanol droplets that carry excess NH(4)(+) ions. Methanol is found to compromise the H-bonding network, resulting in greatly increased rates of ion ejection and solvent evaporation. Considerable differences in the water and methanol escape rates cause time-dependent changes in droplet composition. Segregation occurs at low methanol concentration, such that layered droplets with a methanol-enriched periphery are formed. This phenomenon will enhance the partitioning of analyte molecules, with possible implications for their ESI efficiencies. Solvated ions are ejected from the tip of surface protrusions. Solvent bridging prior to ion secession is more extensive for methanol/water droplets than for purely aqueous systems. The ejection of solvated NH(4)(+) is visualized as diffusion-mediated escape from a metastable basin. The process involves thermally activated crossing of a ~30 kJ mol(-1) free energy barrier, in close agreement with the predictions of the classical ion evaporation model.     

Nucleation free energy of pore formation in an amphiphilic bilayer studied by molecular dynamics simulations

The formation of a pore in a membrane requires a considerable rearrangement of the amphiphilic molecules about to form the bilayer edge surrounding the pore, and hence is accompanied by a steep increase of the free energy. Recent rupture and conductance experiments suggest that this reshuffling process is also responsible for a small energy barrier that stabilizes "prepores" with diameters of less than 1 nm, rendering both the opening and closing of pores an activated process. We use the potential of mean constraint force method to study this free energy profile, as a function of pore radius, in a coarse grained bilayer model. The calculations show that the free energy rises by (15-20) kT during pore opening, making it an extremely rare nucleation event. Although we do not observe a barrier to pore closure, the results do make the existence of such a barrier plausible. For larger pores we find a smooth transition to Litster's model, from which a line tension coefficient of about 3.7 x 10(-11) J m(-1) is deduced.     

Atomistic molecular dynamics simulations of peptide amphiphile self-assembly into cylindrical nanofibers

Relaxation of a self-assembled structure of 144 peptide amphiphile (PA) molecules into cylindrical nanofibers is studied using atomistic molecular dynamics simulations including explicit water with physiological ion concentration. The PA for these studies includes a hydrophobic alkyl chain that is attached to the N-terminus of the sequence SLSLAAAEIKVAV. The self-assembly is initiated with PA molecules in a roughly cylindrical configuration, as suggested from previous experimental and theoretical investigations, and the cylindrical configuration that results is found to be stable during 40 ns simulations. In the converged structure of the resulting nanofiber, the cylinder radius is ～44 ?, a result that is consistent with experimental results. Water and sodium ions can penetrate into the peptide portion of the fiber but not between the alkyl chains. Even though each PA has an identical sequence, a broad distribution of secondary structure is found in the converged structure of the nanofiber. The β-sheet population for the SLSL and IKV segments of the peptide is ～25%, which is consistent with previous circular dichroism results. We also found that the epitope sequence IKVAV is located on the surface of the nanofiber, as designed for the promotion of the neurite growth. Our findings will be useful for designing new PA fibers that have improved bioactive properties.     

Conformational flexibility of the pentasaccharide LNF-2 deduced from NMR spectroscopy and molecular dynamics simulations

Human milk oligosaccharides (HMOs) are important as prebiotics since they stimulate the growth of beneficial bacteria in the intestine and act as receptor analogues that can inhibit the binding of pathogens. The conformation and dynamics of the HMO Lacto-N-fucopentaose 2 (LNF-2), α-L-Fucp-(1 → 4)[β-D-Galp-(1 → 3)]-β-D-GlcpNAc-(1 → 3)-β-D-Galp-(1 → 4)-D-Glcp, having a Lewis A epitope, has been investigated employing NMR spectroscopy and molecular dynamics (MD) computer simulations. 1D (1)H,(1)H-NOESY experiments were used to obtain proton-proton cross-relaxation rates from which effective distances were deduced and 2D J-HMBC and 1D long-range experiments were utilized to measure trans-glycosidic (3)J(CH) coupling constants. The MD simulations using the PARM22/SU01 force field for carbohydrates were carried out for 600 ns with explicit water as solvent which resulted in excellent sampling for flexible glycosidic torsion angles. In addition, in vacuo MD simulations were performed using an MM3-2000 force field, but the agreement was less satisfactory based on an analysis of heteronuclear trans-glycosidic coupling constants. LNF-2 has a conformationally well-defined region consisting of the terminal branched part of the pentasaccharide, i.e., the Lewis A epitope, and a flexible β-D-GlcpNAc-(1 → 3)-β-D-Galp-linkage towards the lactose unit, which is situated at the reducing end. For this β-(1 → 3)-linkage a negative ψ torsion angle is favored, when experimental NMR data is combined with the MD simulation in the analysis. In addition, flexibility on a similar time scale, i.e., on the order of the global overall molecular reorientation, may also be present for the ? torsion angle of the β-D-Galp-(1 → 4)-D-Glcp-linkage as suggested by the simulation. It was further observed from a temperature variation study that some (1)H NMR chemical shifts of LNF-2 were highly sensitive and this study indicates that Δδ/ΔT may be an additional tool for revealing conformational dynamics of oligosaccharides.     

Kinetics of water filling the hydrophobic channels of narrow carbon nanotubes studied by molecular dynamics simulations

The kinetics of water filling narrow single-walled carbon nanotubes was studied using molecular dynamics simulations. The time required to fully fill a nanotube was linear with respect to the tube length. We observed that water molecules could enter into nanotubes of different lengths, either from one end or from both ends. The probability of having a nanotube filled completely from both ends increased exponentially with the tube length. For short tubes, filling usually proceeded from only one end. For long tubes, filling generally proceeded from both tube ends over three stages, i.e., filling from one end, filling from both ends, and filling from both ends with the dipole reorientation of water molecules to give a concerted ordering within the fully filled tube. The water molecules in the partially filled nanotube were hydrogen bonded similarly to those in the fully filled nanotube. Simulations for the reference Lennard-Jones fluid without hydrogen bonds were also performed and showed that the filling behavior of water molecules can be attributed to strong intermolecular hydrogen bonding.     

The effect of the counterion on water mobility in reverse micelles studied by molecular dynamics simulations

In this study, mobility and structure of water molecules in Aerosol OT (bis(2-ethylhexyl) sulfosuccinate, AOT) reverse micelles with water content w0 = 5 and Na+, K+, Cs+ counterions have been explored with molecular dynamics (MD) simulations. Using the Faeder/Ladanyi model (J. Phys. Chem. B, 2000, 104, 1033) of the reverse micelle interior, MD simulations were performed to calculate the self-intermediate scattering function, FS(Q,t), for water hydrogen atoms that could be measured in a quasielastic neutron scattering experiment. Separate intermediate scattering functions FRS(Q,t) and FCMS(Q,t) were determined for rotational and translational motion. We find that the decay of FCMS(Q,t) is nonexponential and our analysis of the MD data indicates that this behavior arises from decreased water mobility for molecules close to the interface and from confinement-induced restrictions on the range of translational displacements. Rotational relaxation also exhibits nonexponential decay, which is consistent with relatively rapid restricted rotation and slower rotational relaxation over the full angular range. Rotational relaxation is anisotropic, with the O-H bond short-time rotational mobility considerably higher than that of the molecular dipole. This behavior is related to the decreased density of water-water hydrogen bonds in the vicinity of the interface compared to core or bulk water. We find that the interfacial mobility of water molecules is quite different for the three counterion types, but that the core mobility exhibits weak counterion dependence. Differences in interfacial mobility are strongly correlated with structural features, especially ion-water coordination, and the extent of disruption by the counterions of the water hydrogen bond network.