Article Information

• PDF (1309 kb)

PubMed

• Articles by P. Auffinger
• Articles by E. Westhof

Related Articles

• Dynamical Properties of a Hydrated Lipid Bilayer from a Mult...

  Dynamical Properties of a Hydrated Lipid Bilayer from a Multinanosecond Molecular Dynamics Simulation Biophysical Journal, Volume 81, Issue 5, 1 November 2001, Pages 2484-2494 Preston B. Moore, Carlos F. Lopez and Michael L. Klein Abstract A fully hydrated dimiristoylphosphatidylcholine (DMPC) bilayer has been studied by a molecular dynamics simulation. The system, which consisted of 64 DMPC molecules and 1792 water molecules, was run in the NVE ensemble at a temperature of 333K for a total of 10ns. The resulting trajectory was used to analyze structural and dynamical quantities. The electron density, bilayer spacing, and order parameters (), based on the AMBER forcefield and SPCE water model are in good agreement with previous calculations and experimental data. The simulation reveals evidence for two types of lateral diffusive behavior: cage hopping and that of a two-dimensional liquid. The lateral diffusion coefficient is 8×10 cm/s. We characterize the rotational motion, and find that the lipid tail rotation (=−0.04rad/ns) is slower then the head group rotation (=2.2rad/ns), which is slower than the overall in plane (=3.2rad/ns) for the lipid molecule. Abstract | Full Text | PDF (213 kb)

• Rotation of Lipids in Membranes: Molecular Dynamics Simulati...

  Rotation of Lipids in Membranes: Molecular Dynamics Simulation, P Spin-Lattice Relaxation, and Rigid-Body Dynamics Biophysical Journal, Volume 94, Issue 8, 15 April 2008, Pages 3074-3083 Jeffery B. Klauda, Mary F. Roberts, Alfred G. Redfield, Bernard R. Brooks and Richard W. Pastor Abstract Molecular dynamics simulations and P-NMR spin-lattice () relaxation rates from 0.022 to 21.1 T of fluid phase dipalmitoylphosphatidylcholine bilayers are compared. Agreement between experiment and direct prediction from simulation indicates that the dominant slow relaxation (correlation) times of the dipolar and chemical shift anisotropy spin-lattice relaxation are ∼10ns and 3ns, respectively. Overall reorientation of the lipid body, consisting of the phosphorus, glycerol, and acyl chains, is well described within a rigid-body model. Wobble, with =1–2×10s, is the primary component of the 10ns relaxation; this timescale is consistent with the tumbling of a lipid-sized cylinder in a medium with the viscosity of liquid hexadecane. The value for the diffusion constant for rotation about the long axis of the lipid body, is difficult to determine precisely because of averaging by fast motions and wobble; it is tentatively estimated to be 1×10s. The resulting /≈0.1 implies that axial rotation is strongly modulated by interactions at the lipid/water interface. Rigid-body modeling and potential of mean force evaluations show that the choline group is relatively uncoupled from the rest of the lipid. This is consistent with the ratio of chemical shift anisotropy and dipolar correlation times reported here and the previous observations that P-NMR lineshapes are axially symmetric even in the gel phase of dipalmitoylphosphatidylcholine. Abstract | Full Text | PDF (322 kb)

• The Mg Binding Sites of the 5S rRNA Loop E Motif as Investig...

  The Mg Binding Sites of the 5S rRNA Loop E Motif as Investigated by Molecular Dynamics Simulations Chemistry & Biology, Volume 10, Issue 6, 1 June 2003, Pages 551-561 Pascal Auffinger, Lukasz Bielecki and Eric Westhof Summary Molecular dynamics simulations have been used to investigate the binding of Mg ions to the deep groove of the eubacterial 5S rRNA loop E. The simulations suggest that long-lived and specific water-mediated interactions established between the hydrated ions and the RNA atoms lining up the binding sites contribute to the stabilization of this motif. The Mg binding specificity is modulated by two factors: (i) a required electrostatic complementarity and (ii) a structural correspondence between the hydrated ion and its binding pocket that can be estimated by its degree of dehydration and the resulting number and lifetime of the intervening water-mediated contacts. Two distinct binding modes for pentahydrated Mg ions that result in a significant freezing of the tumbling motions of the ions are described, and mechanistic details related to the stabilization of nucleic acids by divalent ions are provided. Summary | Full Text | PDF (849 kb)

• …more

Abstract

Multiple molecular dynamics trajectories of the solvated and neutralized 17-residue tRNA(Asp) anticodon hairpin were generated for a total of 3 ns. Explicit treatment of all long-ranged electrostatic interactions by the particle mesh Ewald algorithm, as implemented in the AMBER MD software package, effected a degree of structural stabilization not previously achieved by use of a long 16-A solvent interaction truncation scheme. The increased stability of this multiple molecular dynamics set was appropriate for an in-depth analysis of the six 500-ps-long trajectories and allowed the characterization of a number of key structural interactions. The dynamical behavior of the standard Watson-Crick base pairs, the noncanonical G30-U40 "wobble" base pair, and the psi 32-C38 pseudo-base pair is presented as well as that of two C--H... O hydrogen bonds found to contribute to the array of tertiary interactions that stabilize the seven-nucleotide native loop conformation. The least mobile residue in the loop is U33, which forms the U-turn motif and which participates in several hydrogen-bonding interactions, whereas the most mobile residue is the apical residue G34 at the wobble position, a factor undoubtedly important in its biological function. The set of multiple molecular dynamics trajectories obtained does not converge on a 500-ps time scale to a unique dynamical model but instead describes an ensemble of structural microstates accessible to the system under the present simulation protocol, which is the result of local structural heterogeneity rather than of global conformational changes.