Biomedical subjects

Gregory A Voth

Publications and source records attributed to Gregory A Voth.

At least 19 recordsLinked to original sources

PubMed2006

An accurate and simple quantum model for liquid water.

The path-integral molecular dynamics and centroid molecular dynamics methods have been applied to investigate the behavior of liquid water at ambient conditions starting from a recently developed simple point charge/flexible (SPC/Fw) model. Several quantum structural, thermodynamic, and dynamical properties have been computed and compared to the corresponding classical values, as well as to the available experimental data. The path-integral molecular dynamics simulations show that the inclusion of quantum effects results in a less structured liquid with a reduced amount of hydrogen bonding in comparison to its classical analog. The nuclear quantization also leads to a smaller dielectric constant and a larger diffusion coefficient relative to the corresponding classical values. Collective and single molecule time correlation functions show a faster decay than their classical counterparts. Good agreement with the experimental measurements in the low-frequency region is obtained for the quantum infrared spectrum, which also shows a higher intensity and a redshift relative to its classical analog. A modification of the original parametrization of the SPC/Fw model is suggested and tested in order to construct an accurate quantum model, called q-SPC/Fw, for liquid water. The quantum results for several thermodynamic and dynamical properties computed with the new model are shown to be in a significantly better agreement with the experimental data. Finally, a force-matching approach was applied to the q-SPC/Fw model to derive an effective quantum force field for liquid water in which the effects due to the nuclear quantization are explicitly distinguished from those due to the underlying molecular interactions. Thermodynamic and dynamical properties computed using standard classical simulations with this effective quantum potential are found in excellent agreement with those obtained from significantly more computationally demanding full centroid molecular dynamics simulations. The present results suggest that the inclusion of nuclear quantum effects into an empirical model for water enhances the ability of such model to faithfully represent experimental data, presumably through an increased ability of the model itself to capture realistic physical effects.

Computer Simulation↗

PubMed2006

Multiscale simulation of transmembrane proteins.

Multiscale simulation is employed to examine changes in atomistic-level protein structure due to long wavelength membrane undulations and plane stress fields. An ensemble of atomistic-level simulations of a model of a transmembrane influenza A virus M2 proton channel in a dimyristoylphosphatidylcholine (DMPC) bilayer is coupled to a corresponding mesoscopic model of a DMPC bilayer in an explicit mesoscopic solvent. Structural variations in the key proton gating His37 residues of the M2 channel are examined. Small, but distinct variations in the structure of the His37 residues are observed in both the open and closed states of the channel as a result of the coupling to mesoscopic-level membrane motions.

Dimyristoylphosphatidylcholine↗

PubMed2006

Mechanisms of passive ion permeation through lipid bilayers: insights from simulations.

Multistate empirical valence bond and classical molecular dynamics simulations were used to explore mechanisms for passive ion leakage through a dimyristoyl phosphatidylcholine lipid bilayer. In accordance with a previous study on proton leakage (Biophys. J. 2005, 88, 3095), it was found that the permeation mechanism must be a highly concerted one, in which ion, solvent, and membrane coordinates are coupled. The presence of the ion itself significantly alters the response of those coordinates, suggesting that simulations of transmembrane water structures without explicit inclusion of the ionic solute are insufficient for elucidating transition mechanisms. The properties of H(+), Na(+), OH(-), and bare water molecules in the membrane interior were compared, both by biased sampling techniques and by constructing complete and unbiased transition paths. It was found that the anomalous difference in leakage rates between protons and other cations can be largely explained by charge delocalization effects rather than the usual kinetic picture (Grotthuss hopping of the proton). Permeability differences between anions and cations through phosphatidylcholine bilayers are correlated with suppression of favorable membrane breathing modes by cations.

Computer Simulation↗

PubMed2006

Modeling real dynamics in the coarse-grained representation of condensed phase systems.

This work presents a systematic multiscale methodology to provide a more faithful representation of real dynamics in coarse-grained molecular simulation models. The theoretical formalism is based on the recently developed multiscale coarse-graining (MS-CG) method [S. Izvekov and G. A. Voth, J. Phys. Chem. B. 109, 2469 (2005); J. Chem. Phys. 123, 134105 (2005)] and relies on the generalized Langevin equation approach and its simpler Langevin equation limit. The friction coefficients are determined in multiscale fashion from the underlying all-atom molecular dynamics simulations using force-velocity and velocity-velocity correlation functions for the coarse-grained sites. The diffusion properties in the resulting CG Brownian dynamics simulations are shown to be quite accurate. The time dependence of the velocity autocorrelation function is also well-reproduced relative to the all-atom model if sufficient resolution of the CG sites is implemented.

Computer Simulation↗

PubMed2006

Charge delocalization in proton channels, I: the aquaporin channels and proton blockage.

The explicit contribution to the free energy barrier and proton conductance from the delocalized nature of the excess proton is examined in aquaporin channels using an accurate all-atom molecular dynamics computer simulation model. In particular, the channel permeation free energy profiles are calculated and compared for both a delocalized (fully Grotthuss shuttling) proton and a classical (nonshuttling) hydronium ion along two aquaporin channels, Aqp1 and GlpF. To elucidate the effects of the bipolar field thought to arise from two alpha-helical macrodipoles on proton blockage, free energy profiles were also calculated for computational mutants of the two channels where the bipolar field was eliminated by artificially discharging the backbone atoms. Comparison of the free energy profiles between the proton and hydronium cases indicates that the magnitude of the free energy barrier and position of the barrier peak for the fully delocalized and shuttling proton are somewhat different from the case of the (localized) classical hydronium. The proton conductance through the two aquaporin channels is also estimated using Poisson-Nernst-Planck theory for both the Grotthuss shuttling excess proton and the classical hydronium cation.

Aquaporin 1↗

PubMed2006

Charge delocalization in proton channels, II: the synthetic LS2 channel and proton selectivity.

In this study, the minimalist synthetic LS2 channel is used as a prototype to examine the selectivity of protons over other cations. The free-energy profiles along the transport pathway of LS2 are calculated for three cation species: a realistic delocalized proton (including Grotthuss shuttling)--H(+), a classical (nonshuttling) hydronium--H(3)O(+), and a potassium cation--K(+). The overall barrier for K(+) is approximately twice as large as that for H(+), explaining the >100 times larger maximal ion conductance for the latter, in qualitative agreement with the experimental result. The profile for the classical hydronium is quantitatively intermediate between those of H(+) and K(+) and qualitatively more similar to that of H(+), for which the locations of the peaks are well correlated with the troughs of the pore radius profile. There is a strong correlation between the free-energy profiles and the very different characteristic hydration structures of the three cation species. This work suggests that the passage of various cations through ion channels cannot always be explained by simple electrostatic desolvation considerations.

Biophysics↗

PubMed2006

Atomistic and coarse-grained analysis of double spectrin repeat units: the molecular origins of flexibility.

Spectrin is an ubiquitous protein in metazoan cells, and its flexibility is one of the keys to maintaining cellular structure and organization. Both alpha-spectrin and beta-spectrin polypeptides consist primarily of triple coiled-coil modular repeat units, and two important factors that determine spectrin flexibility are the bending flexibility between two consecutive repeat units and the conformational flexibility of individual repeat units. Atomistic molecular dynamics (MD) simulations are used here to study double spectrin repeat units (DSRUs) from the human erythrocyte beta-spectrin (HEbeta89) and the chicken brain alpha-spectrin (CBalpha1617). From the results of MD simulations, a highly conserved Trp residue in the A-helix of most repeat units that has been suggested to be important in conferring stability to the coiled-coil structures is found not to have a significant effect on the conformational flexibility of individual repeat units. Characterization of the bending flexibility for two consecutive repeats of spectrin via atomistic simulations and coarse-grained (CG) modeling indicate that the bending flexibility is governed by the interactions between the AB-loop of the first repeat unit, the BC-loop of the second repeat unit and the linker region. Specifically, interactions between residues in these regions can lead to a strong directionality in the bending behavior of two repeat units. The biological implications of these finding are discussed.

Amino Acid Sequence↗

PubMed2006

Evaluation of nonlinear quantum time correlation functions within the centroid dynamics formulation.

A method to evaluate nonlinear centroid correlation functions is presented that is amenable to simple numerical computation. It can be implemented with the centroid molecular dynamics method for approximate quantum dynamics with no additional assumptions. Two nonlinear correlation functions are evaluated for a model potential using this scheme and compared with results from exact quantum calculations.

Models, Molecular↗

PubMed2006

Direct observation of Bin/amphiphysin/Rvs (BAR) domain-induced membrane curvature by means of molecular dynamics simulations.

The process of membrane curvature generation by BAR (Bin/amphiphysin/Rvs) domains is thought to involve the plastering of the negatively charged cell membrane to the positively charged concave surface of the BAR domain. Recent work [Peter, B. J., et al. (2004) Science, 303,495-499; Masuda, M., et al. (2006) EMBO J. 25, 2889-2897; and Gallop, J. L., et al. (2006) EMBO J. 25, 2898-2910] has demonstrated the importance of the charged, crescent-shaped surface and the N-terminal amphipathic helices (present in N-BAR domains) for generating membrane curvature. These experiments suggest that curvature is generated by the synergistic action of the N-terminal helices embedding in the lipid bilayer and the charged crescent-shaped dimer acting to "scaffold" membrane curvature. Here, we present atomistic molecular dynamics simulations that directly show membrane binding to the concave face of N-BAR domains, resulting in the generation of local membrane curvature that matches the curvature presented by the BAR domain. These simulations provide direct molecular-scale evidence that BAR domains create curvature by acting as a scaffold, forcing the membrane to locally adopt the intrinsic shape of the BAR domain. We find that BAR domains bind strongly through the maximum curvature surface and, additionally, at an orientation that presents a lesser degree of curvature, thus enabling N-BAR domains to induce a range of local curvatures. Finally, we find that the N-terminal region may play a role in biasing the orientations of N-BAR domains on the membrane surface to those that favor binding to the concave face and subsequent membrane bending.

Animals↗

PubMed2006

Characterization of the solvation and transport of the hydrated proton in the perfluorosulfonic acid membrane nafion.

The solvation and transport properties of the sulfonate-hydronium ion pair have been studied in hydrated Nafion through molecular dynamics simulation. Explicit proton and charge delocalization of the excess proton transport, via the Grotthuss hopping mechanism, were treated using the self-consistent multistate empirical valence bond (SCI-MS-EVB) method. The nature of the sulfonate-hydronium ion pair was characterized through analysis of free-energy profiles. It was found that, in general, the excess proton is solvated between two water molecules of a Zundel moiety while in the contact ion pair position, but then it transitions to an Eigen-like configuration in the solvent-separated pair position. Furthermore, the positive charge associated with the excess proton passes between the contact and solvent-separated ion pair positions through the Grotthuss mechanism rather than simple vehicular diffusion. The total proton diffusion was decomposed into vehicular and Grotthuss components and were found to be of the same relative magnitude, but with a strong negative correlation resulting in a smaller overall diffusion. Correlated motions between the ion pair were examined through the distinct portion of the van Hove correlation function, and a characteristic time scale of approximately 425 ps was observed. Additionally, the association of the hydrated proton with the hydrophobic polymer backbone suggests its amphiphile-like behavior (see Acc. Chem. Res. 2006, 39, 143; Phys. Rev. 1954, 95, 249; J. Chem. Phys. 2005, 123, 084309).

Chemistry, Physical↗

PubMed2006

Tail aggregation and domain diffusion in ionic liquids.

An extended multiscale coarse-graining model for ionic liquids is used to investigate the liquid crystal-like phase in certain ionic liquids. The tail groups of the cations with a sufficient side-chain length are found to aggregate, forming spatially heterogeneous domains, due to the competition between the electrostatic interactions between the charged head groups and the anions and the collective short-range interactions between the neutral tail groups. With a sufficiently long alkyl chain at a low enough temperature, the tail domains remain relatively stable, despite the diffusion of individual ions in the liquid phase. With increasing temperature, the average tail domains begin to diffuse, while beyond a transition temperature, their average density has an almost uniform distribution, although the tail groups still form instantaneous domains.

Carbon↗

PubMed2006

Mesoscopic modeling of bacterial flagellar microhydrodynamics.

A particle-based hybrid method of elastic network model and smooth-particle hydrodynamics has been employed to describe the propulsion of bacterial flagella in a viscous hydrodynamic environment. The method explicitly models the two aspects of bacterial propulsion that involve flagellar flexibility and long-range hydrodynamic interaction of low-Reynolds-number flow. The model further incorporates the molecular organization of the flagellar filament at a coarse-grained level in terms of the 11 protofilaments. Each of these protofilaments is represented by a collection of material points that represent the flagellin proteins. A computational model of a single flexible helical segment representing the filament of a bacterial flagellum is presented. The propulsive dynamics and the flow fields generated by the motion of the model filament are examined. The nature of flagellar deformation and the influence of hydrodynamics in determining the shape of deformations are examined based on the helical filament.

Bacterial Physiological Phenomena↗

PubMed2006

Mixed atomistic and coarse-grained molecular dynamics: simulation of a membrane-bound ion channel.

The recently developed multiscale coarse-graining (MS-CG) method (Izvekov, S.; Voth, G. A. J. Phys. Chem. B 2005, 109, 2469; J. Chem. Phys. 2005, 123, 134105) is used to build a mixed all-atom and coarse-grained (AA-CG) model of the gramicidin A (gA) ion channel embedded in a dimyristoylphosphatidylcholine (DMPC) lipid bilayer and water environment. In this model, the gA peptide was described in full atomistic detail, while the lipid and water molecules were described using coarse-grained representations. The atom-CG and CG-CG interactions in the mixed AA-CG model were determined using the MS-CG method. Molecular dynamics (MD) simulations were performed using the resulting AA-CG model. The results from simulations of the AA-CG model compare very favorably to those from all-atom MD simulations of the entire system. Since the MS-CG method employs a general and systematic approach to obtain effective interactions from the underlying all-atom models, the present approach to rigorously develop mixed AA-CG models has the potential to be extended to many other systems.

Computer Simulation↗

PubMed2006

Examining the influence of linkers and tertiary structure in the forced unfolding of multiple-repeat spectrin molecules.

The unfolding pathways of multiple-repeat spectrin molecules were examined using steered molecular dynamics (SMD) simulations to forcibly unfold double- and triple-repeat spectrin molecules. Although SMD has previously been used to study other repeating-domain proteins, spectrin offers a unique challenge in that the linker connecting repeat units has a definite secondary structure, that of an alpha-helix. Therefore, the boundary conditions imposed on a double- or triple-repeat spectrin must be carefully considered if any relationship to the real system is to be deduced. This was accomplished by imposing additional forces on the system which ensure that the terminal alpha-helices behave as if there were no free noncontiguous helical ends. The results of the SMD simulations highlight the importance of the rupture of the alpha-helical linker on the subsequent unfolding events. Rupture of the linker propagates unfolding in the adjacent repeat units by destabilizing the tertiary structure, ultimately resulting in complete unfolding of the affected repeat unit. Two dominant classes of unfolding pathways are observed after the initial rupture of a linker which involve either rupture of another linker (possibly adjacent) or rupture of the basic tertiary structure of a repeat unit. The relationship between the force response observed on simulation timescales and those of experiment or physiological conditions is also discussed.

Binding Sites↗

PubMed2006

Free energy profiles for H+ conduction in the D-pathway of Cytochrome c Oxidase: a study of the wild type and N98D mutant enzymes.

The molecular mechanism for proton conduction in the D-pathway of Cytochrome c Oxidase (CcO) is investigated through the free energy profile, i.e., potential of mean force (PMF) calculations of both the native enzyme and the N98D mutant. The multistate empirical valence bond (MS-EVB) model was applied to simulate the interaction of an excess proton with the channel environment. In the study of the wild type enzyme, the PMF reveals the previously proposed proton trap inside the channel; it also shows a high free energy barrier against the passage of proton at the entry of the channel, where two conserved asparagines (ASN80/98) may be essential for the gating of proton uptake. We also present data from an investigation of the N98D mutant, which has been previously shown to completely eliminate proton pumping but significantly enhance the oxidase activity in Rhodobacter sphaeroides. These results suggest that mutating Asn98 to negatively charged aspartate will create an unfavorable energy barrier sufficiently high to prevent the overall proton uptake through the D-pathway, whereas with a protonated aspartic acid the proton conduction was found to be accelerated. Plausible explanations for the origin of the uncoupling of proton pumping from the oxidase activity will be discussed.

Animals↗

PubMed2006

A comparative study of imaginary time path integral based methods for quantum dynamics.

The recently introduced approximate many-body quantum simulation method, ring polymer molecular dynamics (RPMD), is compared to the centroid molecular dynamics method (CMD). Comparisons of simulation results for liquid para-hydrogen at two state points and liquid ortho-deuterium at one state point are presented. The calculated quantum correlation functions for the two methods are shown to be in good agreement with one another for a large portion of the time spectrum. However, as the quantum mechanical nature of the system increases, RPMD is less accurate in predicting the kinetic energy of the system than is CMD. A simplified and highly efficient algorithm is proposed which largely corrects this deficiency.

Journal Article↗

PubMed2006

Amphiphilic character of the hydrated proton in methanol-water solutions.

The hydrated proton was studied in methanol-water solutions of varying methanol concentrations using the multistate empirical valence bond simulation method. Amphiphile-like behavior of the hydrated proton was noted from its anisotropic association with the methanol methyl groups. Molecular length immiscibility was also characterized through the enumeration of water and protonated water clusters. Excess proton diffusion was calculated across the varying methanol concentrations and found to be in good agreement with experiment after correcting for nuclear quantum effects.

Anisotropy↗

PubMed2006

Origins of proton transport behavior from selectivity domain mutations of the aquaporin-1 channel.

The permeation free-energy profile and maximum ion conductance of proton transport along the channel of three aquaporin-1 (AQP1) mutants (H180A/R195V, H180A, and R195V) are calculated via molecular dynamics simulations and Poisson-Nernst-Planck theory. The proton dynamics was described by the multistate empirical valence bond (MS-EVB) model. The results reveal three major contributions to the overall free-energy barrier for proton transport in AQP1: 1), the bipolar field, 2), the electrostatic repulsion due to the Arg-195 residue, and 3), the dehydration penalty due to the narrow channel pore. The double mutation (H180A/R195V) drastically drops the overall free-energy barrier by roughly 20 kcal/mol via simultaneously relaxing the direct electrostatic interaction (by R195V) and dehydration effect (by H180A).

Amino Acid Substitution↗