Polarizable Water Model for the Coarse-Grained MARTINI Force Field

Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models

Free Energy Simulations of a GTPase: GTP and GDP Binding to Archaeal Initiation Factor 2

International audienceArchaeal initiation factor 2 (aIF2) is a protein involved in the initiation of protein biosynthesis. In its GTP-bound, "ON" conformation, aIF2 binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To understand the specific binding of GTP and GDP and its dependence on the ON or OFF conformational state of aIF2, molecular dynamics free energy simulations (MDFE) are a tool of choice. However, the validity of the computed free energies depends on the simulation model, including the force field and the boundary conditions, and on the extent of conformational sampling in the simulations. aIF2 and other GTPases present specific difficulties; in particular, the nucleotide ligand coordinates a divalent Mg(2+) ion, which can polarize the electronic distribution of its environment. Thus, a force field with an explicit treatment of electronic polarizability could be necessary, rather than a simpler, fixed charge force field. Here, we begin by comparing a fixed charge force field to quantum chemical calculations and experiment for Mg(2+):phosphate binding in solution, with the force field giving large errors. Next, we consider GTP and GDP bound to aIF2 and we compare two fixed charge force fields to the recent, polarizable, AMOEBA force field, extended here in a simple, approximate manner to include GTP. We focus on a quantity that approximates the free energy to change GTP into GDP. Despite the errors seen for Mg(2+):phosphate binding in solution, we observe a substantial cancellation of errors when we compare the free energy change in the protein to that in solution, or when we compare the protein ON and OFF states. Finally, we have used the fixed charge force field to perform MDFE simulations and alchemically transform GTP into GDP in the protein and in solution. With a total of about 200 ns of molecular dynamics, we obtain good convergence and a reasonable statistical uncertainty, comparable to the force field uncertainty, and somewhat lower than the predicted GTP/GDP binding free energy differences. The sign and magnitudes of the differences can thus be interpreted at a semiquantitative level, and are found to be consistent with the experimental binding preferences of ON- and OFF-aIF2

Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations

Alkali (Li+, Na+, K+, Rb+, and Cs+) and halide (F−, Cl−, Br−, and I−) ions play an important role in many biological phenomena, roles that range from stabilization of biomolecular structure, to influence on biomolecular dynamics, to key physiological influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomolecular structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in solution and the interactions of ions with other molecules. At present, the best force fields for biomolecules employ a simple additive, nonpolarizable, and pairwise potential for atomic interaction. In this work, we describe our efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and solution properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mechanical treatment, our goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodology is general and can be extended to other ions and to polarizable force-field models. Our starting point centered on observations from long simulations of biomolecules in salt solution with the AMBER force fields where salt crystals formed well below their solubility limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Åqvist cation parameters. To provide a more appropriate balance, we reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, we calculated hydration free energies of the solvated ions and also lattice energies (LE) and lattice constants (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4PEW, and SPC/E. In addition to well reproducing the solution and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells

Separation of Betti Reaction Product Enantiomers: Absolute Configuration and Inhibition of Botulinum Neurotoxin A

The racemic product of the Betti reaction of 5-chloro-8-hydroxyquinoline, benzaldehyde, and 2-aminopyridine was separated by chiral HPLC to determine which enantiomer inhibited botulinum neurotoxin serotype A. When the enantiomers unexpectedly proved to have comparable activity, the absolute structures of (+)-(R)-1 and (-)-(S)-1 were determined by comparison of calculated and observed circular dichroism spectra. Molecular modeling studies were undertaken in an effort to understand the observed bioactivity and revealed different ensembles of binding modes, with roughly equal binding energies, for the two enantiomers

Separation of betti reaction product enantiomers: Absolute configuration and inhibition of botulinum neurotoxin A

The racemic product of the Betti reaction of 5-chloro-8-hydroxyquinoline, benzaldehyde, and 2-aminopyridine was separated by chiral HPLC to determine which enantiomer inhibited botulinum neurotoxin serotype A. When the enantiomers unexpectedly proved to have comparable activity, the absolute structures of (+)-(R)-1 and (-)-(S)-1 were determined by comparison of calculated and observed circular dichroism spectra. Molecular modeling studies were undertaken in an effort to understand the observed bioactivity and revealed different ensembles of binding modes, with roughly equal binding energies, for the two enantiomers. © 2011 American Chemical Society

Conformational selection by the aIF2 GTPase: a molecular dynamics study of functional pathways.

International audienceArchaeal initiation factor 2 (aIF2) is a GTPase involved in protein biosynthesis. In its GTP-bound, "ON" conformation, it binds an initiator tRNA and carries it to the ribosome. In its GDP-bound, "OFF" conformation, it dissociates from tRNA. To improve our understanding of the role of each conformational state in the aIF2 "life cycle", we start from the state immediately after GTP hydrolysis, ON:GDP:P(i) (where P(i) is inorganic phosphate), and consider the possible next steps on the pathway to the OFF:GDP product. The first possibility is P(i) dissociation, leading to ON:GDP, which could then relax into OFF:GDP. We use molecular dynamics simulations to compute the P(i) dissociation free energy and show that dissociation is highly favorable. The second possibility is conformational relaxation into the OFF state before P(i) dissociation, to form OFF:GDP:P(i). We estimate the corresponding free energy approximately, 2 ± 3.5 kcal/mol, so that this is an uphill or weakly downhill process. A third possibility is relaxation into another conformation, neither ON nor OFF. Indeed, a third, "MIXED" conformation was seen recently in a crystal structure of the aIF2:GDP:P(i) complex. For this conformational state, P(i) dissociation is weakly unfavorable, in contrast to the ON and OFF states. From this, we will deduce that if the MIXED:GDP complex is not too unstable, the ON:GDP:P(i) → MIXED:GDP:P(i) transformation is a downhill process, which can occur spontaneously. This suggests that the MIXED state could be a functional intermediate