Parallel molecular dynamics simulation in NPT ensemble and its application to general anesthesia

In this paper, we describe a parallel molecular dynamics simulation scheme in NPT ensemble which can be used for studying molecular systems with bond constraints. As an example of application, we choose to investigate general anesthesia whose molecular mechanism is largely unknown. In doing that, we performed two simulations. One is pure lipid bilayer simulation and the other is lipid bilayer simulation with 50 mol% halothane (which is a well-known inhaled anesthetic). From the simulations, we have found that d-spacing increases and area per molecule decreases in the presence of the halothane molecules. We have also found that the halothane molecules are distributed in the hydrophobic tail region, but near headgroup region

[1]  Mitsunori Ikeguchi,et al.  Partial rigid‐body dynamics in NPT, NPAT and NPγT ensembles for proteins and membranes , 2004, J. Comput. Chem..

[2]  David Brown,et al.  A domain decomposition parallel processing algorithm for molecular dynamics simulations of systems of arbitrary connectivity. , 1997 .

[3]  Yan Xu,et al.  Different distribution of fluorinated anesthetics and nonanesthetics in model membrane: a 19F NMR study. , 1997, Biophysical journal.

[4]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[5]  Takao Yamazaki,et al.  A domain decomposition parallel processing algorithm for molecular dynamics simulations of polymers , 1994 .

[6]  Thomas F. Miller,et al.  Symplectic quaternion scheme for biophysical molecular dynamics , 2002 .

[7]  M. Klein,et al.  Effects of anesthetics on the structure of a phospholipid bilayer: molecular dynamics investigation of halothane in the hydrated liquid crystal phase of dipalmitoylphosphatidylcholine. , 1998, Biophysical journal.

[8]  William Smith,et al.  Molecular dynamics on hypercube parallel computers , 1991 .

[9]  Søren Toxvaerd,et al.  Efficient stress relaxation in molecular dynamics simulations of semiflexible n -alkanes , 1998 .

[10]  C. Tanford.  Macromolecules , 1994, Nature.

[11]  K. Laasonen,et al.  Structure, effective pair potential and properties of halothane , 1996 .

[12]  M. Klein,et al.  Distribution of halothane in a dipalmitoylphosphatidylcholine bilayer from molecular dynamics calculations. , 2000, Biophysical journal.

[13]  M Marchi,et al.  Comment on "Efficient stress relaxation in molecular dynamics simulations of semiflexible n-alkanes". , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  R. Eckenhoff.  An inhalational anesthetic binding domain in the nicotinic acetylcholine receptor. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[15]  M. Klein,et al.  Membrane structural perturbations caused by anesthetics and nonimmobilizers: a molecular dynamics investigation. , 2001, Biophysical journal.

[16]  Takao Yamazaki,et al.  A domain decomposition parallel processing algorithm for molecular dynamics simulations of polymers [Comput. Phys. Commun. 83 (1994) 1–13] , 1995 .

[17]  J. Baber,et al.  Distribution of general anesthetics in phospholipid bilayers determined using 2H NMR and 1H-1H NOE spectroscopy. , 1995, Biochemistry.

[18]  Masuhiro Mikami,et al.  Rigid‐body dynamics in the isothermal‐isobaric ensemble: A test on the accuracy and computational efficiency , 2003, J. Comput. Chem..

[19]  Grant S. Heffelfinger,et al.  Parallel atomistic simulations , 2000 .

[20]  Mark E. Tuckerman,et al.  Explicit reversible integrators for extended systems dynamics , 1996 .

[21]  N. Matubayasi,et al.  Reversible molecular dynamics for rigid bodies and hybrid Monte Carlo , 1999 .

[22]  Søren Toxvaerd,et al.  Reply to “Comment on ‘Efficient stress relaxation in molecular dynamics simulations of semiflexible n -alkanes’ ” , 2001 .

[23]  P. Tang,et al.  Amphiphilic sites for general anesthetic action? Evidence from 129Xe-[1H] intermolecular nuclear Overhauser effects. , 1997, Biochimica et biophysica acta.

[24]  Bruce Hendrickson,et al.  A new parallel method for molecular dynamics simulation of macromolecular systems , 1996 .

[25]  M. Tuckerman,et al.  Understanding Modern Molecular Dynamics: Techniques and Applications , 2000 .

[26]  H. C. Andersen.  Rattle: A “velocity” version of the shake algorithm for molecular dynamics calculations , 1983 .

[27]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[28]  W. F. van Gunsteren,et al.  Effect of constraints on the dynamics of macromolecules , 1982 .

[29]  R. Eckenhoff,et al.  Molecular interactions between inhaled anesthetics and proteins. , 1997, Pharmacological reviews.

[30]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[31]  Steve Plimpton,et al.  Fast parallel algorithms for short-range molecular dynamics , 1993 .

[32]  Piero Procacci,et al.  Coordinates scaling and multiple time step algorithms for simulation of solvated proteins in the NPT ensemble , 1998 .

[33]  D. Cafiso,et al.  Contrasting membrane localization and behavior of halogenated cyclobutanes that follow or violate the Meyer-Overton hypothesis of general anesthetic potency. , 1997, Biophysical journal.

[34]  J. Trudell,et al.  Localization of molecular halothane in phospholipid bilayer model nerve membranes. , 1976, Anesthesiology.