Biomedical subjects

Gregory A Voth

Publications and source records attributed to Gregory A Voth.

At least 55 records · Page 3Linked to original sources

PubMed2005

Multi-scale modeling of phase separation in mixed lipid bilayers.

An approach to bridging the phenomenological field theory description of phase separation in binary mixed lipid bilayers with coarse-grained (CG) molecular dynamics (MD) simulation is presented. CG MD simulation is carried out for a 1:1 dipalmitoylphosphatidylcholine/dipalmitoylphosphatidylethanolamine lipid mixture at the liquid-gel phase coexistence condition. The liquid-gel phase separation can be characterized by the bilayer thickness, area per lipid molecule, and orientation parameter of the lipid tails. After a local order parameter is defined using the lipid tail bond orientation parameter, the CG MD data are bridged to a mesoscopic model based on the phenomenological Landau-Ginzberg free-energy functional. All parameters in this mesoscopic model are defined from the information of the phase boundary structure and the distributions of the order parameter in the liquid and gel phases. It is found that the mesoscopic model reproduces the equilibrium properties of the system very well, including collective fluctuations in both phases, spatial correlation functions of the order parameter, and the line tension. The possibility of using a time-dependent Landau-Ginzberg model to mimic the phase-separation dynamics is also investigated, using the relaxation time constant obtained by fitting the time-dependent correlation functions of the order parameter.

Cell Membrane↗

PubMed2005

Ab initio molecular-dynamics simulation of aqueous proton solvation and transport revisited.

The solvation and transport of the hydrated excess proton is studied using the Car-Parrinello molecular-dynamics (CPMD) simulation method. The simulations were performed using BLYP and HCTH gradient-corrected exchange-correlation energy functionals. The fictitious electronic mass was chosen to be small enough so that the underlying water structural and dynamical properties were converged with respect to this important CPMD simulation parameter. An unphysical overstructuring of liquid water in the CPMD simulations using the BLYP functional resulted in the formation of long-lived hydrogen-bonding structures involving the excess proton and a particular (special) water oxygen. The excess proton was observed to be attracted to the special oxygen through the entire length of the BLYP CPMD simulations. Consequently, the excess proton diffusion was limited by the mobility of the special oxygen in the slowly diffusing water network and, in turn, the excess proton self-diffusion coefficient was found to be significantly below the experimental value. On the other hand, the structural properties of liquid water in the HCTH CPMD simulation were seen to be in better agreement with experiment, although the water and excess proton diffusions were still well below the experimental value.

Journal Article↗

PubMed2005

A computational study of the closed and open states of the influenza a M2 proton channel.

In this study, four possible conformations of the His-37 and Trp-41 residues for the closed state of the influenza M2 ion channel were identified by a conformation scan based on a solid-state NMR restraint. In the four conformations, the His-37 residue can be of either the t-160 or t60 rotamer, whereas Trp-41 can be of either the t-105 or t90 rotamer. These conformations were further analyzed by density functional theory calculations and molecular dynamics simulations, and the data indicate that the His-37 residue most likely adopts the t60 rotamer and should be monoprotonated at the delta-nitrogen site, whereas Trp-41 adopts the t90 rotamer. This result is consistent with published experimental data and points to a simple gating mechanism: in the closed state, the His-37 and Trp-41 residues adopt the (t60, t90) conformation, which nearly occludes the pore, preventing nonproton ions from passing through due to the steric and desolvation effects. Moreover, the His-37 tetrad interrupts the otherwise continuous hydrogen-bonding network of the pore water by forcing the water molecules above and below it to adopt opposite orientations, thus adding to the blockage of proton shuttling. The channel can be easily opened by rotating the His-37 chi2 angle from 60 to 0 degrees . This open structure allows pore water to penetrate the constrictive region and to form a continuous water wire for protons to shuttle through, while being still narrow enough to exclude other ions.

Computer Simulation↗

PubMed2005

Multiscale coupling of mesoscopic- and atomistic-level lipid bilayer simulations.

A multiscale method is presented to bridge between the atomistic and mesoscopic membrane systems. The atomistic model in this case is the united atom dimyristoylphosphatidylcholine membrane system, although the method is completely general. Atomistic molecular dynamics provides the expansion modulus which is used to parametrize a mesoscopic elastic membrane model. The resulting elastic membrane model, including explicit mesoscopic solvent, shows appropriate static and dynamic undulation behaviors. Large membranes of approximately 100 nm in length can then be easily simulated using the mesoscopic membrane system. The critical feedback from the mesoscopic system back down to the atomistic-scale system is accomplished by bridging the stress (or surface tension) of a small region in the mesoscopic membrane to the corresponding atomistic membrane system. Because of long length-scale modes of membranes such as undulation and buckling, the local tension responds differently from the frame tension, when subjected to external perturbations. The effect of these membrane modes is shown for the stress response of a local membrane region and therefore the atomistic membrane system. In addition, certain equilibrium static and dynamic properties of stand-alone and multiscale coupled systems are presented for several different membrane sizes. Although static properties such as two-dimensional pair-correlation function and order parameters show no noticeable discrepancy for the different systems, lipid self-diffusion and the rotational relaxation of lipid dipoles have a strong dependence on the membrane size (or long-wavelength membrane motions), which is properly modeled by the present multiscale method.

Algorithms↗

PubMed2005

Computer simulation of explicit proton translocation in cytochrome c oxidase: the D-pathway.

Proton translocation in the D-pathway of cytochrome c oxidase has been studied by a combination of classical molecular dynamics and the multistate empirical valence bond methodology. This approach allows for explicit Grotthuss proton hopping between water molecules. According to mutagenesis experiments, the role of proton donor/acceptor along the D-pathway is carried by the highly conserved residue Glu-242. The present multistate empirical valence bond simulations indicate that the protonation/deprotonation state of Glu-242 is strongly coupled to the distance of proton propagation in the D-pathway. The proton was seen to travel the full length of the D-pathway when Glu-242 was deprotonated; however, it was trapped halfway along the path when Glu-242 was protonated. Further investigation in terms of both proton dynamical properties and free energy calculations for the pathway of proton transport provides evidence for a two-step proton transport mechanism in the D-pathway.

Animals↗

PubMed2005

Effective force field for liquid hydrogen fluoride from ab initio molecular dynamics simulation using the force-matching method.

A recently developed force-matching method for obtaining effective force fields for condensed matter systems from ab initio molecular dynamics (MD) simulations has been applied to fit a simple nonpolarizable two-site pairwise force field for liquid hydrogen fluoride. The ab initio MD in this case was a Car-Parrinello (CP) MD simulation of 64 HF molecules at nearly ambient conditions within the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure included a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The performance of the force-match potential was examined for the gas-phase dimer and for the liquid phase at various temperatures. The model was able to reproduce correctly the bent structure and energetics of the gas-phase dimer, while the results for the structural properties, self-diffusion, vibrational spectra, density, and thermodynamic properties of liquid HF were compared to both experiment and the CP MD simulation. The force-matching model performs well in reproducing nearly all of the liquid properties as well as the aggregation behavior at different temperatures. The model is computationally cheap and compares favorably to many more computationally expensive potential energy functions for liquid HF.

Chemistry, Physical↗

PubMed2005

A linear-scaling self-consistent generalization of the multistate empirical valence bond method for multiple excess protons in aqueous systems.

An extension to the multistate empirical valence bond (MS-EVB) method is presented in this paper that is capable of treating multiple excess protons within the context of molecular-dynamics simulation. The computational cost of the method scales linearly with respect to the number of excess protons. Calculations for a 0.44 M HCl systems are carried out to illustrate the multiproton extension of the MS-EVB method. A significant decrease in the Eigen-type H(9)O(4)(+) cation is observed in the contact ion-pair configuration formed between Cl(-) and hydronium ions.

Journal Article↗

PubMed2005

Coupling field theory with continuum mechanics: a simulation of domain formation in giant unilamellar vesicles.

Domain formation is modeled on the surface of giant unilamellar vesicles using a Landau field theory model for phase coexistence coupled to elastic deformation mechanics (e.g., membrane curvature). Smooth particle applied mechanics, a form of smoothed particle continuum mechanics, is used to solve either the time-dependent Landau-Ginzburg or Cahn-Hilliard free-energy models for the composition dynamics. At the same time, the underlying elastic membrane is modeled using smooth particle applied mechanics, resulting in a unified computational scheme capable of treating the response of the composition fields to arbitrary deformations of the vesicle and vice versa. The results indicate that curvature coupling, along with the field theory model for composition free energy, gives domain formations that are correlated with surface defects on the vesicle. In the case that external deformations are included, the domain structures are seen to respond to such deformations. The present simulation capability provides a significant step forward toward the simulation of realistic cellular membrane processes.

Biomechanical Phenomena↗

PubMed2005

A coarse-grained model for double-helix molecules in solution: spontaneous helix formation and equilibrium properties.

A new reductionist coarse-grained model is presented for double-helix molecules in solution. As with such models for lipid bilayers and micelles, the level of description is both particulate and mesoscopic. The particulate (bead-and-spring) nature of the model makes for a simple implementation in standard molecular dynamics simulation codes and allows for investigation of thermomechanic properties without preimposing any (form of) response function. The mesoscopic level of description--where groups of atoms are condensed into coarse-grained beads--causes long-range interactions to be effectively screened, which greatly enhances the efficiency and scalability of simulations. Without imposing local or global order parameters, a linear initial configuration of the model molecule spontaneously assembles into a double helix due to the interplay between three contributions: hydrophobic/hydrophilic interactions between base pairs, backbone, and solvent; phosphate-phosphate repulsion along the backbone; and favorable base-pair stacking energy. We present results for the process of helix formation as well as for the equilibrium properties of the final state, and investigate how both depend on the input parameters. The current model holds promise for two routes of investigation: First, within a limited set of generic parameters, the effect of local (atomic-scale) perturbations on overall helical properties can be systematically studied. Second, since the efficiency allows for a direct simulation of both small and large (>100 base pairs) systems, the model presents a testground for systematic coarse-graining methods.

Base Pairing↗

PubMed2005

Excess proton solvation and delocalization in a hydrophilic pocket of the proton conducting polymer membrane nafion.

Solvation properties of the hydrated excess proton are studied in a hydrophilic pocket of Nafion 117 through a series of molecular dynamics simulations. The multistate empirical valence bond (MS-EVB) methodology, which enables the delocalization of the excess proton through the Grotthuss hopping mechanism, was employed for one of the excess protons in the simulation cell. Simulations were performed such that "classical" nondissociable hydronium cations and a single excess proton treated with the MS-EVB methodology were at a concentration ratio of 39:1. Two degrees of hydration of the Nafion polymer electrolyte membrane were simulated, each displaying the same marked difference between the solvation structures of the classical versus MS-EVB treated (Grotthuss shuttling) excess proton species. These differences are attributed to the solvent dynamics needed to transfer the cation between the solvent separated and contact pair positions about the sulfonic acid counterion. The results demonstrate that it is generally impossible to describe the low pH conditions in the hydrophilic domains of Nafion without the explicit treatment of Grotthuss delocalization in the underlying molecular dynamics model for the excess protons.

Fluorocarbon Polymers↗

PubMed2005

A multiscale coarse-graining method for biomolecular systems.

A new approach is presented for obtaining coarse-grained (CG) force fields from fully atomistic molecular dynamics (MD) trajectories. The method is demonstrated by applying it to derive a CG model for the dimyristoylphosphatidylcholine (DMPC) lipid bilayer. The coarse-graining of the interparticle force field is accomplished by an application of a force-matching procedure to the force data obtained from an explicit atomistic MD simulation of the biomolecular system of interest. Hence, the method is termed a "multiscale" CG (MS-CG) approach in which explicit atomistic-level forces are propagated upward in scale to the coarse-grained level. The CG sites in the lipid bilayer application were associated with the centers-of-mass of atomic groups because of the simplicity in the evaluation of the forces acting on them from the atomistic data. The resulting CG lipid bilayer model is shown to accurately reproduce the structural properties of the phospholipid bilayer.

Chemistry, Physical↗

PubMed2005

Protons may leak through pure lipid bilayers via a concerted mechanism.

Protons are known to permeate pure lipid bilayers at a rate that is anomalous compared to those of other small monovalent cations. The prevailing mechanism via which they cross the membrane is still unclear, and it is unknown how to probe the mechanism directly by experiment. One of the more popular theories assumes the formation of membrane-spanning single-file water wires providing a matrix along which the protons can "hop" over the barrier. However, free energy calculations on such structures (without the presence of an excess proton) suggest that this mechanism alone cannot account for the observed permeation rates. We use the multistate empirical valence bond method to directly study water structures surrounding a (delocalized) excess proton on its way through the membrane. We find that membrane-spanning networks, rather than single-file chains, are formed around the proton. We also find that such structures are considerably stabilized in the presence of the proton, with lifetimes of several hundreds of picoseconds. The observed structures are suggestive of a new, concerted, mechanism and provide some direction for further investigation.

Biophysics↗

PubMed2005

Fast centroid molecular dynamics: a force-matching approach for the predetermination of the effective centroid forces.

A fast centroid molecular dynamics (CMD) methodology is proposed in which the effective centroid forces are predetermined through a force-matching algorithm applied to a standard path integral molecular dynamics simulation. The resulting method greatly reduces the computational cost of generating centroid trajectories, thus extending the applicability of CMD. The method is applied to the study of liquid para-hydrogen at two state points and liquid ortho-deuterium at one state point. The static and dynamical results are compared to those obtained from full adiabatic CMD simulations and found to be in excellent agreement for all three systems; the transport properties are also compared to experiment and found to have a similar level of agreement.

Journal Article↗

PubMed2005

A bond-order analysis of the mechanism for hydrated proton mobility in liquid water.

Bond-order analysis is introduced to facilitate the study of cooperative many-molecule effects on proton mobility in liquid water, as simulated using the multistate empirical valence-bond methodology. We calculate the temperature dependence for proton mobility and the total effective bond orders in the first two solvation shells surrounding the H(5)O(2) (+) proton-transferring complex. We find that proton-hopping between adjacent water molecules proceeds via this intermediate, but couples to hydrogen-bond dynamics in larger water clusters than previously anticipated. A two-color classification of these hydrogen bonds leads to an extended mechanism for proton mobility.

Journal Article↗

PubMed2004

The dynamic stress responses to area change in planar lipid bilayer membranes.

The viscoelastic properties of planar phospholipid (dimyristoylphosphatidylcholine) bilayer membranes at 308 K are studied, many of them for the first time, using the nonequilibrium molecular dynamics simulation (NEMD) method for membrane area change. First, we present a unified formulation of the intrinsic three-dimensional (3D) and apparent in-plane viscoelastic moduli associated with area change based on the constitutive relations for a uniaxial system. The NEMD simulations of oscillatory area change process are then used to obtain the frequency-domain moduli. In the 4-250 GHz range, the intrinsic 3D elastic moduli of 20-27 kbar and viscous moduli of 0.2-9 kbar are found with anisotropy and monotonic frequency dispersion. In contrast, the apparent in-plane elastic moduli (1-9 kbar) are much smaller than, and the viscous moduli (2-6 kbar) comparable to, their 3D counterparts, due to the interplay between the lateral and normal relaxations. The time-domain relaxation functions, separately obtained by applying stepwise strains, can be fit by 4-6 exponential decay modes spanning subpicosecond to nanosecond timescale and are consistent with the frequency-domain results. From NEMD with varying strain amplitude, the linear constitutive model is shown to be valid up to 6 and 20% area change for the intrinsic 3D elastic and viscous responses, respectively, and up to 20% area change for the apparent in-plane viscoelasticity. Inclusion of a gramicidin A dimer (approximately 1 mol %) yields similar response properties with possibly smaller (<10%) viscous moduli. Our results agree well with available data from ultrasonic experiments, and demonstrate that the third dimension (thickness) of the planar lipid bilayer is integral to the in-plane viscoelasticity.

Dimyristoylphosphatidylcholine↗

PubMed2004

A centroid molecular dynamics study of liquid para-hydrogen and ortho-deuterium.

Centroid molecular dynamics (CMD) is applied to the study of collective and single-particle dynamics in liquid para-hydrogen at two state points and liquid ortho-deuterium at one state point. The CMD results are compared with the results of classical molecular dynamics, quantum mode coupling theory, a maximum entropy analytic continuation approach, pair-product forward- backward semiclassical dynamics, and available experimental results. The self-diffusion constants are in excellent agreement with the experimental measurements for all systems studied. Furthermore, it is shown that the method is able to adequately describe both the single-particle and collective dynamics of quantum liquids.

Journal Article↗

PubMed2004

Coupling field theory with mesoscopic dynamical simulations of multicomponent lipid bilayers.

A method for simulating a two-component lipid bilayer membrane in the mesoscopic regime is presented. The membrane is modeled as an elastic network of bonded points; the spring constants of these bonds are parameterized by the microscopic bulk modulus estimated from earlier atomistic nonequilibrium molecular dynamics simulations for several bilayer mixtures of DMPC and cholesterol. The modulus depends on the composition of a point in the elastic membrane model. The dynamics of the composition field is governed by the Cahn-Hilliard equation where a free energy functional models the coupling between the composition and curvature fields. The strength of the bonds in the elastic network are then modulated noting local changes in the composition and using a fit to the nonequilibrium molecular dynamics simulation data. Estimates for the magnitude and sign of the coupling parameter in the free energy model are made treating the bending modulus as a function of composition. A procedure for assigning the remaining parameters in the free energy model is also outlined. It is found that the square of the mean curvature averaged over the entire simulation box is enhanced if the strength of the bonds in the elastic network are modulated in response to local changes in the composition field. We suggest that this simulation method could also be used to determine if phase coexistence affects the stress response of the membrane to uniform dilations in area. This response, measured in the mesoscopic regime, is already known to be conditioned or renormalized by thermal undulations.

Computer Simulation↗

PubMed2004

Mesoscopic lateral diffusion in lipid bilayers.

The lateral diffusion in bilayers is modeled with a multiscale mesoscopic simulation. The methodology consists of two simulations, where the first employs atomistic models to obtain exact results for the mesoscopic model. The second simulation takes the results obtained from the first to parameterize an effective force field that is employed in a new coarse-grained model. The multiscale aspect of this scheme occurs at the point where the microscopic time-averaged results of the first simulation are employed to parameterize the second simulation that operates in a higher spatial and temporal domain. The results of both simulation schemes give quantitative information on the details associated with lipid lateral diffusion.

Cholesterol↗