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

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

Continuum Approaches to Understanding Ion and Peptide Interactions with the Membrane

Experimental and computational studies have shown that cellular membranes deform to stabilize the inclusion of transmembrane (TM) proteins harboring charge. Recent analysis suggests that membrane bending helps to expose charged and polar residues to the aqueous environment and polar head groups. We previously used elasticity theory to identify membrane distortions that minimize the insertion of charged TM peptides into the membrane. Here, we extend our work by showing that it also provides a novel, computationally efficient method for exploring the energetics of ion and small peptide penetration into membranes. First, we show that the continuum method accurately reproduces energy profiles and membrane shapes generated from molecular simulations of bare ion permeation at a fraction of the computational cost. Next, we demonstrate that the dependence of the ion insertion energy on the membrane thickness arises primarily from the elastic properties of the membrane. Moreover, the continuum model readily provides a free energy decomposition into components not easily determined from molecular dynamics. Finally, we show that the energetics of membrane deformation strongly depend on membrane patch size both for ions and peptides. This dependence is particularly strong for peptides based on simulations of a known amphipathic, membrane binding peptide from the human pathogen Toxoplasma gondii. In total, we address shortcomings and advantages that arise from using a variety of computational methods in distinct biological contexts