Biomedical subjects

Gregory A Voth

Publications and source records attributed to Gregory A Voth.

At least 37 records · Page 2Linked to original sources

PubMed2006

Multiscale coarse-graining of ionic liquids.

A recently developed multiscale coarse-graining (MS-CG) approach for obtaining coarse-grained force fields from fully atomistic molecular dynamics simulation is applied to the challenging case of the EMIM+NO3- ionic liquid. The force-matching in the MS-CG methodology is accomplished with an explicit separation of bonded and nonbonded forces. While the nonbonded forces are adopted from this force-matching approach, the bonded forces are obtained from fitting the statistical configurational data from the atomistic simulations. The many-body electronic polarizability is also successfully broken into effective pair interactions. With a virial constraint fixing the system pressure, the MS-CG models rebuild satisfactory structural and thermodynamic properties for different temperatures. The MS-CG model developed from a modest atomistic simulation is therefore suitable for simulating much larger systems, because the coarse-grained models show significant time integration efficiency. This approach is expected to be general for coarse-graining other ionic liquids, as well as many other liquid-state systems. The limitations of the present coarse-graining procedure are also discussed.

Algorithms↗

PubMed2006

A second generation mesoscopic lipid bilayer model: connections to field-theory descriptions of membranes and nonlocal hydrodynamics.

A new mesoscopic membrane model is developed in order to examine long-wavelength structural and dynamical membrane phenomena. Two different explicit mesoscopic solvent models are employed. The first mesoscopic solvent is denoted the big liquid oscillating blob system, which is parametrized to model water at a coarse-grained level and is motivated by a Langevin-like approach; the resulting membrane dynamics predict a solvent viscosity dependence consistent with the known viscosity of water. The second mesoscopic solvent is a Weeks-Chandler-Anderson model. Here, it is found that the correct mesoscopic hydrodynamic scaling of the membrane undulation dynamics is still preserved, although accelerated. When the behavior of the two membranes in close proximity to one another is examined, very little correlated motion is observed. However, the theoretically predicted scaling of the entropic undulation energy is confirmed, demonstrating that the entropic interaction between two membranes becomes increasingly repulsive with decreasing separation.

Algorithms↗

PubMed2006

Structure of the liquid-vacuum interface of room-temperature ionic liquids: a molecular dynamics study.

Molecular dynamics simulations for the liquid-vacuum interface of the ionic liquid 1-ethyl-3-methylimidazolium nitrate (EMIM+/NO3-) were performed for both electronically polarizable and nonpolarizable potential energy surfaces. The interfacial structural properties, such as the oscillation in the number density profile, the orientational ordering, and the local clustering in the interfacial region, were calculated. The simulations with both the polarizable and nonpolarizable model demonstrate the existence of an inhomogeneous interfacial structure normal to the surface layer. It was found for both models that the ethyl tail group on EMIM+ is likely to protrude outward from the surface. In the outmost surface layer, the cation is likely to lie on the surface with the imidazolium ring parallel to the interface, while there is a second region with enhanced density from that in the bulk where the cation preferably slants with the imidazolium ring tending to be perpendicular to the surface. The results also reveal that the electronic polarization effect is important for the ionic liquid interface. It is found that the cation is likely to be segregated at the ionic liquid surface for the polarizable model, while for the nonpolarizable model, the anion is found to be more likely to exhibit such behavior. The surface tension of the polarizable model (58.5 +/- 0.5 mN/m) is much smaller than that of the nonpolarizable model (82.7 +/- 0.6 mN/m), in better agreement with extrapolated experimental measurements on similar ionic liquid systems.

Journal Article↗

PubMed2006

A multistate empirical valence bond description of protonatable amino acids.

The multistate empirical valence bond (MS-EVB) model, which was developed for molecular dynamics simulations of proton transport in water and biomolecular systems, is extended for the modeling of protonatable amino acid residues in aqueous environments, specifically histidine and glutamic acid. The parameters of the MS-EVB force field are first determined to reproduce the geometries and energetics of the gas phase amino acid-water clusters. These parameters are then optimized to reproduce experimental pK(a) values. The free energy profiles for acid ionization and the corresponding pK(a) values are calculated by MS-EVB molecular dynamics simulations utilizing the umbrella sampling technique, with the center of excess charge coordinate chosen as the dissociation reaction coordinate. A general procedure for fitting the MS-EVB parameters is formulated, which allows for the parametrization of other amino acid residues with protonatable groups and the subsequent use of the MS-EVB approach for molecular dynamics simulations of proton transfer processes in proteins involving protonation/deprotonation of the protonatable amino acid groups.

Amino Acids↗

PubMed2006

Flexible simple point-charge water model with improved liquid-state properties.

In order to introduce flexibility into the simple point-charge (SPC) water model, the impact of the intramolecular degrees of freedom on liquid properties was systematically studied in this work as a function of many possible parameter sets. It was found that the diffusion constant is extremely sensitive to the equilibrium bond length and that this effect is mainly due to the strength of intermolecular hydrogen bonds. The static dielectric constant was found to be very sensitive to the equilibrium bond angle via the distribution of intermolecular angles in the liquid: A larger bond angle will increase the angle formed by two molecular dipoles, which is particularly significant for the first solvation shell. This result is in agreement with the work of Hochtl et al. [J. Chem. Phys. 109, 4927 (1998)]. A new flexible simple point-charge water model was derived by optimizing bulk diffusion and dielectric constants to the experimental values via the equilibrium bond length and angle. Due to the large sensitivities, the parametrization only slightly perturbs the molecular geometry of the base SPC model. Extensive comparisons of thermodynamic, structural, and kinetic properties indicate that the new model is much improved over the standard SPC model and its overall performance is comparable to or even better than the extended SPC model.

Algorithms↗

PubMed2006

The properties of ion-water clusters. II. Solvation structures of Na+, Cl-, and H+ clusters as a function of temperature.

Ion-water-cluster properties are investigated both through the multistate empirical valence bond potential and a polarizable model. Equilibrium properties of the ion-water clusters H+(H2O)100, Na+(H2O)100, Na+(H2O)20, and Cl-(H2O)17 in the temperature region 100-450 K are explored using a hybrid parallel basin-hopping and tempering algorithm. The effect of the solid-liquid phase transition in both caloric curves and structural distribution functions is investigated. It is found that sodium and chloride ions largely reside on the surface of water clusters below the cluster melting temperature but are solvated into the interior of the cluster above the melting temperature, while the solvated proton was found to have significant propensity to reside on or near the surface in both the liquid- and solid-state clusters.

Algorithms↗

PubMed2006

Computer simulation of proton solvation and transport in aqueous and biomolecular systems.

The excess proton in aqueous media is critical to many aspects of chemistry, biology, and materials science. This species is at the heart of the most elementary of chemical (e.g., acid-base) and biological (e.g., bioenergetics) concepts, yet to this day, it remains mysterious, surprising, and often misunderstood. In this Account, our efforts to describe excess proton solvation and transport through computer modeling and simulation will be described. Results will be summarized for several important systems, as obtained from the multistate empirical valence bond (MS-EVB) approach, which allows for the explicit treatment of (Grotthuss) proton shuttling and charge delocalization.

Biological Transport↗

PubMed2005

Structure of hydrated Na-Nafion polymer membranes.

We use molecular dynamics simulations to investigate the structure of the hydrated Na-Nafion membranes. The membrane is "prepared" by starting with the Nafion chains placed on a cylinder having the water inside it. Minimizing the energy of the system leads to a filamentary hydrophilic domain whose structure depends on the degree of hydration. At 5 wt % water the system does not have enough water molecules to solvate all the ions that could be formed by the dissociation of the -SO3Na groups. As a result, the -SO3Na groups aggregate with the water to form very small droplets that do not join into a continuous phase. The size of the droplets is between 5 and 8 A. As the amount of water present in the membrane is increased, the membrane swells, and SO3Na has an increasing tendency to dissociate into ions. Furthermore, a transition to a percolating hydrophilic network is observed. In the percolating structure, the water forms irregular curvilinear channels branching in all directions. The typical dimension of the cross section of these channels is about 10-20 A. Calculated neutron scattering from the simulated system is in qualitative agreement with experiment. In all simulations, the pendant sulfonated perfluorovinyl side chains of the Nafion hug the walls of the hydrophilic channel, while the sulfonate groups point toward the center of the hydrophilic phase. The expulsion of the side chains from the hydrophilic domain is favored because it allows better interaction between the water molecules. We have also examined the probability of finding water molecules around the Na+ and the -SO3(-) ions as well as the probability of finding other water molecules next to a given water molecule. These probabilities are much broader than those found in bulk water or for one ion in bulk water (calculated with the potentials used in the present simulation). This is due to the highly inhomogeneous nature of the material contained in the small hydrophilic pores.

Journal Article↗

PubMed2005

An efficient and accurate implementation of centroid molecular dynamics using a Gaussian approximation.

An approximate method for Centroid Molecular Dynamics (CMD) is presented which uses a Gaussian approximation. The resulting method, called Gaussian CMD (GCMD), is 100-1000 times faster than CMD because it replaces explicit path-integral sampling, which is the most time-consuming part of CMD, with a Gaussian averaging, which can be done analytically. Several methods for computing the Gaussian width parameter in the GCMD approach are also presented. This new method is shown to give satisfactory results for the position correlation function in simple one-dimensional systems when CMD itself is consistent with the exact result. The GCMD and CMD results are also compared for the case of 1-dimensional systems coupled to harmonic baths, with good success. GCMD is further compared to CMD with good success for liquid para-hydrogen at two different temperatures, 14 K and 25 K, and for ortho-deuterium at 20.7 K.

Journal Article↗

PubMed2005

Coarse-grained modeling of the actin filament derived from atomistic-scale simulations.

A coarse-grained (CG) procedure that incorporates the information obtained from all-atom molecular dynamics (MD) simulations is presented and applied to actin filaments (F-actin). This procedure matches the averaged values and fluctuations of the effective internal coordinates that are used to define a CG model to the values extracted from atomistic MD simulations. The fluctuations of effective internal coordinates in a CG model are computed via normal-mode analysis (NMA), and the computed fluctuations are matched with the atomistic MD results in a self-consistent manner. Each actin monomer (G-actin) is coarse-grained into four sites, and each site corresponds to one of the subdomains of G-actin. The potential energy of a CG G-actin contains three bonds, two angles, and one dihedral angle; effective harmonic bonds are used to describe the intermonomer interactions in a CG F-actin. The persistence length of a CG F-actin was found to be sensitive to the cut-off distance of assigning intermonomer bonds. Effective harmonic bonds for a monomer with its third nearest neighboring monomers are found to be necessary to reproduce the values of persistence length obtained from all-atom MD simulations. Compared to the elastic network model, incorporating the information of internal coordinate fluctuations enhances the accuracy and robustness for a CG model to describe the shapes of low-frequency vibrational modes. Combining the fluctuation-matching CG procedure and NMA, the achievable time- and length scales of modeling actin filaments can be greatly enhanced. In particular, a method is described to compute the force-extension curve using the CG model developed in this work and NMA. It was found that F-actin is easily buckled under compressive deformation, and a writhing mode is developed as a result. In addition to the bending and twisting modes, this novel writhing mode of F-actin could also play important roles in the interactions of F-actin with actin-binding proteins and in the force-generation process via polymerization.

Actins↗

PubMed2005

Extending a spectrin repeat unit. I: linear force-extension response.

Nonequilibrium molecular dynamics simulations were used to calculate the elastic properties of a spectrin repeat unit. A contiguous alpha-helical linker was constructed by employing periodic boundary conditions, allowing a novel scheme for evaluating the thermodynamic force as a function of extension. By measuring the force-extension response under small extensions, spectrin was observed to behave primarily as an elastic material with a spring constant of 1700 +/- 100 pN/nm. The implications of this spring constant, in terms of the properties of the spectrin tetramer, are also discussed.

Algorithms↗

PubMed2005

Extending a spectrin repeat unit. II: rupture behavior.

A spectrin repeat unit was subject to extension using cyclic expansion nonequilibrium molecular dynamics. Periodic boundary conditions were used to examine the effects of the contiguous alpha-helical linker on the force response. The measured force-extension curve shows a linear increase in the force response when the spectrin repeat unit is extended by approximately 0.4 nm. After that point, the force response peaks and subsequently declines. The peak in the force response marks the point where the spectrin repeat unit undergoes a change in its material properties from a strongly elastic material to a mostly viscous one, on the timescales of the simulations. The force peak is also correlated with rupture of the alpha-helical linker, and is likely the event responsible for the peaks in the sawtooth-pattern force-extension curves measured by atomic force microscopy experiments. Rupture of the linker involves simultaneously breaking approximately four hydrogen bonds that maintain the alpha-helical linker. After this initial rupture, the linker undergoes simple helix-to-coil transitions as the spectrin repeat unit continues to be extended. The implications of linker rupture in the interpretation of unfolding and atomic force microscopy experiments are also discussed.

Animals↗

PubMed2005

Probing the molecular-scale lipid bilayer response to shear flow using nonequilibrium molecular dynamics.

A nonequilibrium molecular dynamics simulation of the response of dimyristoylphosphatidylcholine (DMPC) bilayers to a solvent shear flow is presented. Application of shear flow to planar, stationary DMPC bilayers results in a redistribution of the membrane density profile along the bilayer normal due to the alignment of the lipids in the direction of flow and an increase in average lipid chain length. An increase in the intermolecular and intramolecular order of the lipids in response to the shear flow is also observed. This study provides groundwork for understanding the mechanism of the full response of lipid bilayers to externally imposed solvent shear flows, beginning with the response in the absence of collective lipid motions such as undulations and bilayer flow.

Lipid Bilayers↗

PubMed2005

Multiscale coarse graining of liquid-state systems.

A methodology is described to systematically derive coarse-grained (CG) force fields for molecular liquids from the underlying atomistic-scale forces. The coarse graining of an interparticle force field is accomplished by the application of a force-matching method to the trajectories and forces obtained from the atomistic trajectory and force data for the CG sites of the targeted system. The CG sites can be associated with the centers of mass of atomic groups because of the simplicity in the evaluation of forces acting on these sites from the atomistic data. The resulting system is called a multiscale coarse-grained (MS-CG) representation. The MS-CG method for liquids is applied here to water and methanol. For both liquids one-site and two-site CG representations without an explicit treatment of the long-ranged electrostatics have been derived. In addition, for water a two-site model having the explicit long-ranged electrostatics has been developed. To improve the thermodynamic properties (e.g., pressure and density) for the MS-CG models, the constraint for the instantaneous virial was included into the force-match procedure. The performance of the resulting models was evaluated against the underlying atomistic simulations and experiment. In contrast with existing approaches for coarse graining of liquid systems, the MS-CG approach is general, relies only on the interatomic interactions in the reference atomistic system.

Computer Simulation↗

PubMed2005

Systematic coarse-graining of nanoparticle interactions in molecular dynamics simulation.

A recently developed multiscale coarse-graining procedure [Izvekov, S.; Voth, G. A. J. Phys. Chem. B 2005, 109, 2469] is extended to derive coarse-grained models for nanoparticles. The methodology is applied to C(60) and to carbonaceous nanoparticles produced in combustion environments. The coarse-graining of the interparticle force field is accomplished applying a force-matching procedure to data obtained from trajectories and forces from all-atom MD simulations. The CG models are shown to reproduce accurately the structural properties of the nanoparticle systems studied, while allowing for MD simulations of much larger self-assembled nanoparticle systems.

Models, Molecular↗

PubMed2005

Unique spatial heterogeneity in ionic liquids.

A multiscale coarse-graining model for ionic liquids has been extended to investigate the unique aggregation of cations in ionic liquids through computer simulation. It has been found that, with sufficiently long side chains, the tail groups of cations aggregate to form spatially heterogeneous domains, while headgroups of the cations and the anions distribute as uniformly as possible. This is understood as the result of competition between the charged electrostatic interactions between headgroups and anions and the collective short-range interactions between the neutral tail groups. This aggregation can help to explain a number of experimentally observed physical phenomena in ionic liquids.

Journal Article↗

PubMed2005

Allostery of actin filaments: molecular dynamics simulations and coarse-grained analysis.

The structural and mechanical properties of monomeric actin (G-actin), the trimer nucleus, and actin filaments (F-actins) are determined as a function of the conformation of the DNase I-binding loop (DB loop) by using all-atom molecular dynamics simulations and coarse-grained (CG) analysis. Recent x-ray structures of ADP-bound G-actin (G-ADP) by Otterbein et al. [Otterbein, L. R., Graceffa, P. & Dominguez, R. (2001) Science 293, 708-711] and ATP-bound G-actin (G-ATP) by Graceffa and Dominguez [Graceffa, P. & Dominguez, R. (2003) J. Biol. Chem. 278, 34172-34180] indicate that the DB loop of actin does not have a well defined secondary structure in the ATP state but folds into an alpha-helix in the ADP state. MD simulations and CG analysis indicate that such a helical DB loop significantly weakens the intermonomer interactions of actin assemblies and thus leads to a wider, shorter, and more disordered filament. The computed persistence lengths of F-actin composed of G-ATP (16 microm) and of G-ADP (8.5 microm) agree well with the experimental values for the two states. Therefore, the loop-to-helix transition of the DB loop may be one of the factors that lead to the changes in structural and mechanical properties of F-actin after ATP hydrolysis. This result may provide a direct connection between the conformational changes of an actin monomer and the structural and mechanical properties of the cytoskeleton. The information provided by MD simulations also helps to understand the possible origin of the special features of actin dynamics.

Actin Cytoskeleton↗

PubMed2005

The properties of ion-water clusters. I. The protonated 21-water cluster.

The ab initio atom-centered density-matrix propagation approach and the multistate empirical valence bond method have been employed to study the structure, dynamics, and rovibrational spectrum of a hydrated proton in the "magic" 21 water cluster. In addition to the conclusion that the hydrated proton tends to reside on the surface of the cluster, with the lone pair on the protonated oxygen pointing "outwards," it is also found that dynamical effects play an important role in determining the vibrational properties of such clusters. This result is used to analyze and complement recent experimental and theoretical studies.

Adsorption↗