All-Atom Continuous Constant pH Molecular Dynamics With Particle Mesh Ewald and Titratable Water.

Development of a pH stat to properly control solution pH in biomolecular simulations has been a long-standing goal in the community. Toward this goal recent years have witnessed the emergence of the so-called constant pH molecular dynamics methods. However, the accuracy and generality of these methods have been hampered by the use of implicit-solvent models or truncation-based electrostatic schemes.

Introducing titratable water to all-atom molecular dynamics at constant pH.

Recent development of titratable coions has paved the way for realizing all-atom molecular dynamics at constant pH. To further improve physical realism, here we describe a technique in which proton titration of the solute is directly coupled to the interconversion between water and hydroxide or hydronium. We test the new method in replica-exchange continuous constant pH molecular dynamics simulations of three proteins, HP36, BBL, and HEWL.

Charge-leveling and proper treatment of long-range electrostatics in all-atom molecular dynamics at constant pH.

Recent development of constant pH molecular dynamics (CpHMD) methods has offered promise for adding pH-stat in molecular dynamics simulations. However, until now the working pH molecular dynamics (pHMD) implementations are dependent in part or whole on implicit-solvent models. Here we show that proper treatment of long-range electrostatics and maintaining charge neutrality of the system are critical for extending the continuous pHMD framework to the all-atom representation.

Unraveling A Trap-and-Trigger Mechanism in the pH-Sensitive Self-Assembly of Spider Silk Proteins.

When the major ampulate spidroins (MaSp1) are called upon to form spider dragline silk, one of nature's most amazing materials, a small drop in pH must occur. Using a state-of-the-art simulation technique, constant pH molecular dynamics, we discovered a few residues that respond to the pH signal in the dimerization of the N-terminal domain (NTD) of MaSp1 which is an integral step in the fiber assembly. At neutral pH the deprotonation of Glu79 and Glu119 leads to water penetration and structural changes at the monomer-monomer binding interface.

Thermodynamic coupling of protonation and conformational equilibria in proteins: theory and simulation.

Ionization-coupled conformational phenomena are ubiquitous in biology. However, quantitative characterization of the underlying thermodynamic cycle comprised of protonation and conformational equilibria has remained an elusive goal. Here we use theory and continuous constant pH molecular dynamics (CpHMD) simulations to provide a thermodynamic description for the coupling of proton titration and conformational exchange between two distinct states of a protein.

Simulating pH titration of a single surfactant in ionic and nonionic surfactant micelles.

Calculation of surfactant pK(a)'s in micelles is a challenging task using traditional electrostatic methods due to the lack of structural data and information regarding the effective dielectric constant. Here we test the implicit- and hybrid-solvent-based continuous constant pH molecular dynamics (CpHMD) methods for predicting the pK(a) shift of a lauric acid solubilized in three micelles: dodecyl sulfate (DS), dodecyltrimethylammonium (DTA), and dodecyltriethylene glycol ether (DE3). Both types of simulations are able to reproduce the observed positive pK(a) shifts for the anionic DS and nonionic DE3 micelles.

Continuous Constant pH Molecular Dynamics in Explicit Solvent with pH-Based Replica Exchange.

A computational tool that offers accurate pKa values and atomically detailed knowledge of protonation-coupled conformational dynamics is valuable for elucidating mechanisms of energy transduction processes in biology, such as enzyme catalysis and electron transfer as well as proton and drug transport. Toward this goal we present a new technique of embedding continuous constant pH molecular dynamics within an explicit-solvent representation. In this technique we make use of the efficiency of the generalized-Born (GB) implicit-solvent model for estimating the free energy of protein solvation while propagating conformational dynamics using the more accurate explicit-solvent model.

Toward accurate prediction of pKa values for internal protein residues: the importance of conformational relaxation and desolvation energy.

Proton uptake or release controls many important biological processes, such as energy transduction, virus replication, and catalysis. Accurate pK(a) prediction informs about proton pathways, thereby revealing detailed acid-base mechanisms. Physics-based methods in the framework of molecular dynamics simulations not only offer pK(a) predictions but also inform about the physical origins of pK(a) shifts and provide details of ionization-induced conformational relaxation and large-scale transitions.

Molecular dynamics simulations of ionic and nonionic surfactant micelles with a generalized Born implicit-solvent model.

In recent years, all-atom and coarse-grained models have been developed and applied to simulations of micelles and biological membranes. Here, we explore the question of whether a combined all-atom representation of surfactant molecules and continuum description of solvent based on the generalized Born model can be used to study surfactant micelles. Specifically, we report the parameterization of the GBSW model with a surface-area dependent nonpolar solvation energy term for dodecyl sulfate, dodecyl tetramethylammonium, and dodecyl triethyleneglycol ether molecules.

Probing the strand orientation and registry alignment in the propagation of amyloid fibrils.

Detailed knowledge of the structure and growth mechanism of amyloid fibrils is important for understanding the disease process. Recently, solid-state NMR and other spectroscopic data have revealed the equilibrium organization of the tertiary structure of fibrils formed by various segments of beta-amyloid peptides. A three-step "dock-and-lock" mechanism for fibril growth has been proposed on the basis of the kinetic data.

Predicting pKa values with continuous constant pH molecular dynamics.

Knowledge of pK(a) values is important for understanding structure and function relationships in proteins. Over the past two decades, theoretical methods for pK(a) calculations have been mainly based on macroscopic models, in which the protein is considered as a low-dielectric cavity embedded in a high-dielectric continuum. In recent years, constant pH molecular dynamics methods have been developed based on a microscopic description of the protein.