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

Identification and Characterization of Novel Mutations in the Human Gene Encoding the Catalytic Subunit Calpha of Protein Kinase A (PKA)

The genes PRKACA and PRKACB encode the principal catalytic (C) subunits of protein kinase A (PKA) Cα and Cβ, respectively. Cα is expressed in all eukaryotic tissues examined and studies of Cα knockout mice demonstrate a crucial role for Cα in normal physiology. We have sequenced exon 2 through 10 of PRKACA from the genome of 498 Norwegian donors and extracted information about PRKACA mutations from public databases. We identified four interesting nonsynonymous point mutations, Arg45Gln, Ser109Pro, Gly186Val, and Ser263Cys, in the Cα1 splice variant of the kinase. Cα variants harboring the different amino acid mutations were analyzed for kinase activity and regulatory (R) subunit binding. Whereas mutation of residues 45 and 263 did not alter catalytic activity or R subunit binding, mutation of Ser109 significantly reduced kinase activity while R subunit binding was unaltered. Mutation of Cα Gly186 completely abrogated kinase activity and PKA type I but not type II holoenzyme formation. Gly186 is located in the highly conserved DFG motif of Cα and mutation of this residue to Val was predicted to result in loss of binding of ATP and Mg2+, which may explain the kinetic inactivity. We hypothesize that individuals born with mutations of Ser109 or Gly186 may be faced with abnormal development and possibly severe disease

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

Tools to map target genes of bacterial two‐component system response regulators

Studies on bacterial physiology are incomplete without knowledge of the signalling and regulatory systems that a bacterium uses to sense and respond to its environment. Two-component systems (TCSs) are among the most prevalent bacterial signalling systems, and they control essential and secondary physiological processes; however, even in model organisms, we lack a complete understanding of the signals sensed, the phosphotransfer partners and the functions regulated by these systems. In this review, we discuss several tools to map the genes targeted by transcriptionally acting TCSs. Many of these tools have been used for studying individual TCSs across diverse species, but systematic approaches to delineate entire signalling networks have been very few. Since genome sequences and high-throughput technologies are now readily available, the methods presented here can be applied to characterize the entire DNA-binding TCS signalling network in any bacterial species and are especially useful for non-model environmental bacteria

Quantitative Analysis of Repertoire-Scale Immunoglobulin Properties in Vaccine-Induced B-Cell Responses

Recent advances in the next-generation sequencing of B-cell receptors (BCRs) enable the characterization of humoral responses at a repertoire-wide scale and provide the capability for identifying unique features of immune repertoires in response to disease, vaccination, or infection. Immunosequencing now readily generates 103–105 sequences per sample; however, statistical analysis of these repertoires is challenging because of the high genetic diversity of BCRs and the elaborate clonal relationships among them. To date, most immunosequencing analyses have focused on reporting qualitative trends in immunoglobulin (Ig) properties, such as usage or somatic hypermutation (SHM) percentage of the Ig heavy chain variable (IGHV) gene segment family, and on reducing complex Ig property distributions to simple summary statistics. However, because Ig properties are typically not normally distributed, any approach that fails to assess the distribution as a whole may be inadequate in (1) properly assessing the statistical significance of repertoire differences, (2) identifying how two repertoires differ, and (3) determining appropriate confidence intervals for assessing the size of the differences and their potential biological relevance. To address these issues, we have developed a technique that uses Wilcox’ robust statistics toolbox to identify statistically significant vaccine-specific differences between Ig repertoire properties. The advantage of this technique is that it can determine not only whether but also where the distributions differ, even when the Ig repertoire properties are non-normally distributed. We used this technique to characterize murine germinal center (GC) B-cell repertoires in response to a complex Ebola virus-like particle (eVLP) vaccine candidate with known protective efficacy. The eVLP-mediated GC B-cell responses were highly diverse, consisting of thousands of clonotypes. Despite this staggering diversity, we identified statistically significant differences between non-immunized, vaccine only, and vaccine-plus-adjuvant groups in terms of Ig properties, including IGHV-family usage, SHM percentage, and characteristics of the BCR complementarity-determining region. Most notably, our analyses identified a robust eVLP-specific feature—enhanced IGHV8-family usage in B-cell repertoires. These findings demonstrate the utility of our technique in identifying statistically significant BCR repertoire differences following vaccination. More generally, our approach is potentially applicable to a wide range of studies in infection, vaccination, auto-immunity, and cancer

Probing the Donor and Acceptor Substrate Specificity of the γ-Glutamyl Transpeptidase

γ-Glutamyl transpeptidase (GGT) is a two-substrate enzyme that plays a central role in glutathione metabolism and is a potential target for drug design. GGT catalyzes the cleavage of γ-glutamyl donor substrates and the transfer of the γ-glutamyl moiety to an amine of an acceptor substrate or water. Although structures of bacterial GGT have revealed details of the protein–ligand interactions at the donor site, the acceptor substrate site is relatively undefined. The recent identification of a species-specific acceptor site inhibitor, OU749, suggests that these inhibitors may be less toxic than glutamine analogues. Here we investigated the donor and acceptor substrate preferences of <i>Bacillus anthracis</i> GGT (CapD) and applied computational approaches in combination with kinetics to probe the structural basis of the enzyme’s substrate and inhibitor binding specificities and compare them with human GGT. Site-directed mutagenesis studies showed that the R432A and R520S variants exhibited 6- and 95-fold decreases in hydrolase activity, respectively, and that their activity was not stimulated by the addition of the l-Cys acceptor substrate, suggesting an additional role in acceptor binding and/or catalysis of transpeptidation. Rat GGT (and presumably HuGGT) has strict stereospecificity for l-amino acid acceptor substrates, while CapD can utilize both l- and d-acceptor substrates comparably. Modeling and kinetic analysis suggest that R520 and R432 allow two alternate acceptor substrate binding modes for l- and d-acceptors. R432 is conserved in <i>Francisella tularensis</i>, <i>Yersinia pestis</i>, <i>Burkholderia mallei</i>, <i>Helicobacter pylori</i> and <i>Escherichia coli</i>, but not in human GGT. Docking and MD simulations point toward key residues that contribute to inhibitor and acceptor substrate binding, providing a guide to designing novel and specific GGT inhibitors

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