97 research outputs found
The influence of geometry, surface character and flexibility on the permeation of ions and water through biological pores
A hydrophobic constriction site can act as an efficient barrier to ion and
water permeation if its diameter is less than the diameter of an ion's first
hydration shell. This hydrophobic gating mechanism is thought to operate in a
number of ion channels, e.g. the nicotinic receptor, bacterial mechanosensitive
channels (MscL and MscS) and perhaps in some potassium channels (e.g. KcsA,
MthK, and KvAP). Simplified pore models allow one to investigate the primary
characteristics of a conduction pathway, namely its geometry (shape, pore
length, and radius), the chemical character of the pore wall surface, and its
local flexibility and surface roughness. Our extended (ca. 0.1 \mu s) molecular
dynamic simulations show that a short hydrophobic pore is closed to water for
radii smaller than 0.45 nm. By increasing the polarity of the pore wall (and
thus reducing its hydrophobicity) the transition radius can be decreased until
for hydrophilic pores liquid water is stable down to a radius comparable to a
water molecule's radius. Ions behave similarly but the transition from
conducting to non-conducting pores is even steeper and occurs at a radius of
0.65 nm for hydrophobic pores. The presence of water vapour in a constriction
zone indicates a barrier for ion permeation. A thermodynamic model can explain
the behaviour of water in nanopores in terms of the surface tensions, which
leads to a simple measure of "hydrophobicity" in this context. Furthermore,
increased local flexibility decreases the permeability of polar species. An
increase in temperature has the same effect, and we hypothesise that both
effects can be explained by a decrease in the effective solvent-surface
attraction which in turn leads to an increase in the solvent-wall surface free
energy.Comment: Peer reviewed article appeared in Physical Biology
http://www.iop.org/EJ/abstract/1478-3975/1/1/005
Prediction of hydration free energies for aliphatic and aromatic chloro derivatives using molecular dynamics simulations with the OPLS-AA force field
International audienceAll-atom molecular dynamics computer simulations were used to blindly predict the hydration free energies of a range of chloro-organic compounds as part of the SAMPL3 challenge. All compounds were parameterized within the framework of the OPLS-AA force field, using an established protocol to compute the absolute hydration free energy via a windowed free energy perturbation approach and thermodynamic integration. Three different approaches to deriving partial charge parameters were pursued: (1) using existing OPLS-AA atom types and charges with minor adjustments of partial charges on equivalent connecting atoms; (2) calculation of quantum mechanical charges via geometry optimization, followed by electrostatic potential (ESP) fitting, using Jaguar at the LMP2/cc-pVTZ(-F) level; and (3) via geometry optimization and CHelpG charges (Gaussian03 at the HF/6-31G* level), followed by two-stage RESP fitting. Protocol 3 generated the most accurate predictions with a root mean square (RMS) error of 1.2 kcalĀ·mol ā1 for the entire data set. It was found that the deficiency of the standard OPLS-AA parameters, protocol 1 (RMS error 2.4 kcalĀ·mol ā1 overall), was mostly due to compounds with more than three chlorine substituents on an aromatic ring. For this latter subset, the RMS errors were 1.4 kcalĀ·mol ā1 (protocol 3) and 4.3 kcalĀ·mol ā1 (protocol 1), respectively. We propose new OPLS-AA atom types for aromatic carbon and chlorine atoms in rings with ā„ 4 Cl-substituents that perform better than the best QM-based approach, resulting in an RMS error of 1.2 kcalĀ·mol ā1 for these difficult compounds
Path Similarity Analysis: a Method for Quantifying Macromolecular Pathways
Diverse classes of proteins function through large-scale conformational
changes; sophisticated enhanced sampling methods have been proposed to generate
these macromolecular transition paths. As such paths are curves in a
high-dimensional space, they have been difficult to compare quantitatively, a
prerequisite to, for instance, assess the quality of different sampling
algorithms. The Path Similarity Analysis (PSA) approach alleviates these
difficulties by utilizing the full information in 3N-dimensional trajectories
in configuration space. PSA employs the Hausdorff or Fr\'echet path
metrics---adopted from computational geometry---enabling us to quantify path
(dis)similarity, while the new concept of a Hausdorff-pair map permits the
extraction of atomic-scale determinants responsible for path differences.
Combined with clustering techniques, PSA facilitates the comparison of many
paths, including collections of transition ensembles. We use the closed-to-open
transition of the enzyme adenylate kinase (AdK)---a commonly used testbed for
the assessment enhanced sampling algorithms---to examine multiple microsecond
equilibrium molecular dynamics (MD) transitions of AdK in its substrate-free
form alongside transition ensembles from the MD-based dynamic importance
sampling (DIMS-MD) and targeted MD (TMD) methods, and a geometrical targeting
algorithm (FRODA). A Hausdorff pairs analysis of these ensembles revealed, for
instance, that differences in DIMS-MD and FRODA paths were mediated by a set of
conserved salt bridges whose charge-charge interactions are fully modeled in
DIMS-MD but not in FRODA. We also demonstrate how existing trajectory analysis
methods relying on pre-defined collective variables, such as native contacts or
geometric quantities, can be used synergistically with PSA, as well as the
application of PSA to more complex systems such as membrane transporter
proteins.Comment: 9 figures, 3 tables in the main manuscript; supplementary information
includes 7 texts (S1 Text - S7 Text) and 11 figures (S1 Fig - S11 Fig) (also
available from journal site
Prediction of cyclohexane-water distribution coefficients for the SAMPL5 data set using molecular dynamics simulations with the OPLS-AA force field
International audienceAll-atom molecular dynamics (MD) simulations were used to predict water-cyclohexane distribution coefficients D cw of a range of small molecules as part of the SAMPL5 blind prediction challenge. Molecules were parameterized with the trans-ferable all-atom OPLS-AA force field, which required the derivation of new parameters for sulfamides and heterocycles and validation of cyclohexane parameters as a solvent. The distribution coefficient was calculated from the solvation free energies of the compound in water and cyclohexane. Absolute solvation free energies were computed by an established protocol using windowed alchemical free energy perturbation with thermodynamic integration. This protocol resulted in an overall root mean square error (RMSE) in log D cw of almost 4 log units and an overall signed error of ā3 compared to experimental data. There was no substantial overall difference in accuracy between simulating in NV T and NPT ensembles. The signed error suggests a systematic error but the experimental D cw data on their own are insufficient to Manuscript Click here to download Manuscript sampl5-manuscript.pdf Click here to view linked References 2 Ian M. Kenney et al. uncover the source of this error. Preliminary work suggests that the major source of error lies in the hydration free energy calculations
Ligandbook ā an online repository for small and drug-like molecule force field parameters
International audienceLigandbook is a public database and archive for force field parameters of small and drug-like molecules. It is a repository for parameter sets that are part of published work but are not easily available to the community otherwise. Parameter sets can be downloaded and immediately used in molecular dynamics simulations. The sets of parameters are versioned with full histories and carry unique identifiers to facilitate reproducible research. Text-based search on rich metadata and chemical substructure search allow precise identification of desired compounds or functional groups. Ligandbook enables the rapid set up of reproducible molecular dynamics simulations of ligands and protein-ligand complexes. Availability and Implementation: Ligandbook is available online at https://ligandbook.org and supports all modern browsers. Parameters can be searched and downloaded without registration, including access through a programmatic RESTful API. Deposition of files requires free user registration. Ligandbook is implemented in the PHP Symfony2 framework with TCL scripts using the CACTVS toolkit
SAMPL6: calculation of macroscopic pKa values from ab initio quantum mechanical free energies
International audienceMacroscopic pKa values were calculated for all compounds in the SAMPL6 blind prediction challenge, based on quantum chemical calculations with a continuum solvation model and a linear correction derived from a small training set. Microscopic pKa values were derived from the gas-phase free energy difference between protonated and deprotonated forms together with the Conductor-like Polarizable Continuum Solvation Model and the experimental solvation free energy of the proton. pH-dependent microstate free energies were obtained from the microscopic pKas with a maximum likelihood estimator and appropriately summed to yield macroscopic pKa values or microstate populations as function of pH. We assessed the accuracy of three approaches to calculate the microscopic pKas: direct use of the quantum mechanical free energy differences and correction of the direct values for short-comings in the QM solvation model with two different linear models that we independently derived from a small training set of 38 compounds with known pKa. The predictions that were corrected with the linear models had much better accuracy [root-mean-square error (RMSE) 2.04 and 1.95 pKa units] than the direct calculation (RMSE 3.74). Statistical measures indicate that some systematic errors remain, likely due to differences in the SAMPL6 data set and the small training set with respect to their interactions with water. Overall, the current approach provides a viable physics-based route to estimate macroscopic pKa values for novel compounds with reasonable accuracy
- ā¦