Dataset Information

Exploring ion permeation energetics in gramicidin A using polarizable charge equilibration force fields.

ABSTRACT: All-atom molecular dynamics simulations have been applied in the recent past to explore the free energetics underlying ion transport processes in biological ion channels. Roux and co-workers, Kuyucak and co-workers, Busath and co-workers, and others have performed rather elegant and extended time scale molecular dynamics simulations using current state-of-the-art fixed-charge (nonpolarizable) force fields to calculate the potential of mean force defining the equilibrium flux of ions through prototypical channels such as gramicidin A. An inescapable conclusion of such studies has been the gross overestimation of the equilibrium free energy barrier, generally predicted to be from 10 to 20 kcal/mol depending on the force field and simulation protocol used in the calculation; this translates to an underestimation of experimentally measurable single channel conductances by several orders of magnitude. Next-generation polarizable force fields have been suggested as possible alternatives for more quantitative predictions of the underlying free energy surface in such systems. (1) Presently, we consider ion permeation energetics in the gramicidin A channel using a novel polarizable force field. Our results predict a peak barrier height of 6 kcal/mol relative to the channel entrance; this is significantly lower than the uncorrected value of 12 kcal/mol for nonpolarizable force fields such as GROMOS and CHARMM27 which do not account for electronic polarization. These results provide promising initial indications substantiating the long-conjectured importance of polarization effects in describing ion-protein interactions in narrow biological channels.

SUBMITTER: Patel S

PROVIDER: S-EPMC2818143 | biostudies-literature | 2009 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Exploring ion permeation energetics in gramicidin A using polarizable charge equilibration force fields.

Patel Sandeep S Davis Joseph E JE Bauer Brad A BA

Journal of the American Chemical Society 20091001 39

All-atom molecular dynamics simulations have been applied in the recent past to explore the free energetics underlying ion transport processes in biological ion channels. Roux and co-workers, Kuyucak and co-workers, Busath and co-workers, and others have performed rather elegant and extended time scale molecular dynamics simulations using current state-of-the-art fixed-charge (nonpolarizable) force fields to calculate the potential of mean force defining the equilibrium flux of ions through prot ...[more]

PMID: 19788320

Similar Datasets

Project description:Polarizable force fields for lipid and solvent environments are used for molecular dynamics simulations of a fully hydrated dipalmitoylphosphatidylcholine (DPPC) bilayer and gramicidin A (gA) dimer embedded in a dimyristoylphosphatidylcholine (DMPC) bilayer. The lipid bilayer is modeled using the CHARMM charge equilibration (CHEQ) polarizable force field for lipids and the TIP4P-FQ force field to represent solvent. For the DPPC bilayer system, results are compared to the same system simulated using the nonpolarizable CHARMM27r (C27r) force field and TIP3P water. Calculated atomic and electron density profiles, head group orientations as measured by the phosphorus-nitrogen vector orientation, and deuterium order parameters are found to be consistent with previous simulations and with experiment. The CHEQ model exhibits greater water penetration into the bilayer interior, as demonstrated by the potential of mean force calculated from the water density profile. This is a result of the variation of the water molecular dipole from 2.55 D in the bulk to 1.88 D in the interior. We discuss this finding in the context of previous studies (both simulation and experiment) that have investigated the extent of penetration of water into DPPC bilayers. We also discuss the effects of including explicit polarization on the water dipole moment variation as a function of distance from the bilayer. We show distributions of atomic charges over the course of the simulation since the CHEQ model allows the charges to fluctuate. We have calculated the interfacial dipole potential, which the CHEQ model predicts to be 0.95 V compared to 0.86 V as predicted by the C27r model. We also discuss dielectric permittivity profiles and the differences arising between the two models. We obtain bulk values of 72.77 for the CHEQ model (TIP4P-FQ water) and 91.22 for C27r (TIP3P), and values approaching unity in the membrane interior. Finally, we present results of simulations of gA embedded in a DMPC bilayer using the CHEQ model and discuss structural properties.

Project description:Ion descriptors in molecular mechanics models are calibrated against reference data on ion-water interactions. It is then typically assumed that these descriptors will also satisfactorily describe interactions of ions with other functional groups, such as those present in biomolecules. However, several studies now demonstrate that this transferability assumption produces, in many different cases, large errors. Here we address this issue in a representative polarizable model and focus on transferability of cationic interactions from water to a series of alcohols. Both water and alcohols use hydroxyls for ion-coordination, and, therefore, this set of molecules constitutes the simplest possible case of transferability. We obtain gas phase reference data systematically from "gold-standard" quantum Monte Carlo and CCSD(T) methods, followed by benchmarked vdW-corrected DFT. We learn that the original polarizable model yields large gas phase water ? alcohol transferability errors - the RMS and maximum errors are 2.3 and 5.1 kcal/mol, respectively. These errors are, nevertheless, systematic in that ion-alcohol interactions are overstabilized, and systematic errors typically imply that some essential physics is either missing or misrepresented. A comprehensive analysis shows that when both low- and high-field responses of ligand dipole polarization are described accurately, then transferability improves significantly - the RMS and maximum errors in the gas phase reduce, respectively, to 0.9 and 2.5 kcal/mol. Additionally, predictions of condensed phase transfer free energies also improve. Nevertheless, within the limits of the extrathermodynamic assumptions necessary to separate experimental estimates of salt dissolution into constituent cationic and anionic contributions, we note that the error in the condensed phase is systematic, which we attribute, at least, partially to the parametrization in long-range electrostatics. Overall, this work demonstrates a rational approach to boosting transferability of ionic interactions that will be applicable broadly to improving other polarizable and nonpolarizable models.

Dataset Information

Exploring ion permeation energetics in gramicidin A using polarizable charge equilibration force fields.

Publications

Exploring ion permeation energetics in gramicidin A using polarizable charge equilibration force fields.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure