A detailed representation of electrostatic energy in prediction of sequence and pH dependence of protein stability

A molecular mechanics model, previously validated in applications to structure prediction, is shown to reproduce experiment in predictions of protein ionization state, and in predictions of sequence and pH dependence of protein stability. Over a large dataset, 1876 values of ΔΔG of folding, the RMSD is 1.34 kcal/mol. Using an alternative measure of accuracy, either the sign of the calculated ΔΔG agrees with experiment or the absolute value of the deviation is less than 1.0 kcal/mol, 1660 of 1876 data points (88.5%) pass the condition. Relative to models used previously in computer‐aided protein design, the concept, we propose, most responsible for the performance of our model, and for the extensibility to non‐neutral values of pH, is the treatment of electrostatic energy. The electronic structure of the protein is modeled using distributed atomic multipoles. The structured liquid state of the solvent is modeled using a dielectric continuum. A modification to the energetics of the reaction field, induced by the protein in the dielectric continuum, attempts to account for preformed multipoles of solvent water molecules and ions. An adjustable weight (with optimal value.141) applied to the total vacuum energy accounts implicitly for electronic polarization. A threshold distance, beyond which pairwise atomic interactions are neglected, is not used. In searches through subspaces of sequences and conformations, efficiency remains acceptable for useful applications. Proteins 2014; 82:2497–2511. © 2014 Wiley Periodicals, Inc.     

Implicit solvent simulations of peptide interactions with anionic lipid membranes

A recently developed implicit membrane model (IMM1) is supplemented with a Gouy-Chapman term describing counterion-screened electrostatic interactions of a solute with negatively charged membrane lipids. The new model is tested on peptides that bind to anionic membranes. Pentalysine binds just outside the plane of negative charge, whereas Lys-Phe peptides insert their aromatic rings into the hydrophobic core. Melittin and magainin 2 bind more strongly to anionic than to neutral membranes and in both cases insert their hydrophobic residues into the hydrocarbon core. The third domain of Antennapedia homeodomain (penetratin) binds as an alpha-helix in the headgroup region. Cardiotoxin II binds strongly to anionic membranes but marginally to neutral ones. In all cases, the location and configuration of the peptides are consistent with experimental data, and the effective energy changes upon binding compare favorably with experimental binding free energies. The model opens the way to exploring the effect of membrane charge on the location, conformation, and dynamics of a large variety of biologically active peptides on membranes.     

Energetics of protein folding

The energetics of protein folding determine the 3D structure of a folded protein. Knowledge of the energetics is needed to predict the 3D structure from the amino acid sequence or to modify the structure by protein engineering. Recent developments are discussed: major factors are reviewed and auxiliary factors are discussed briefly. Major factors include the hydrophobic factor (burial of non-polar surface area) and van der Waals interactions together with peptide hydrogen bonds and peptide solvation. The long-standing model for the hydrophobic factor (free energy change proportional to buried non-polar surface area) is contrasted with the packing-desolvation model and the approximate nature of the proportionality between free energy and apolar surface area is discussed. Recent energetic studies of forming peptide hydrogen bonds (gas phase) are reviewed together with studies of peptide solvation in solution. Closer agreement is achieved between the 1995 values for protein unfolding enthalpies in vacuum given by Lazaridis-Archontis-Karplus and Makhatadze-Privalov when the solvation enthalpy of the peptide group is taken from electrostatic calculations. Auxiliary factors in folding energetics include salt bridges and side-chain hydrogen bonds, disulfide bridges, and propensities to form alpha-helices and beta-structure. Backbone conformational entropy is a major energetic factor which is discussed only briefly for lack of knowledge.     

Monte Carlo folding of trans-membrane helical peptides in an implicit generalized Born membrane

An efficient Monte Carlo (MC) algorithm using concerted backbone rotations is combined with a recently developed implicit membrane model to simulate the folding of the hydrophobic transmembrane domain M2TM of the M2 protein from influenza A virus and Sarcolipin at atomic resolution. The implicit membrane environment is based on generalized Born theory and has been calibrated against experimental data. The MC sampling has previously been used to fold several small polypeptides and been shown to be equivalent to molecular dynamics (MD). In combination with a replica exchange algorithm, M2TM is found to form continuous membrane spanning helical conformations for low temperature replicas. Sarcolipin is only partially helical, in agreement with the experimental NMR structures in lipid bilayers and detergent micelles. Higher temperature replicas exhibit a rapidly decreasing helicity, in agreement with expected thermodynamic behavior. To exclude the possibility of an erroneous helical bias in the simulations, the model is tested by sampling a synthetic Alanine-rich polypeptide of known helicity. The results demonstrate there is no overstabilization of helical conformations, indicating that the implicit model captures the essential components of the native membrane environment for M2TM and Sarcolipin.     

Simple solvation potential for coarse-grained models of proteins

We formulate a simple solvation potential based on a coarsed-grained representation of amino acids with two spheres modeling the C(alpha) atom and an effective side-chain centroid. The potential relies on a new method for estimating the buried area of residues, based on counting the effective number of burying neighbors in a suitable way. This latter quantity shows a good correlation with the buried area of residues computed from all atom crystallographic structures. We check the discriminatory power of the solvation potential alone to identify the native fold of a protein from a set of decoys and show the potential to be considerably selective.     

Nonpolar Alkanes Modify Lithium‐Ion Solvation for Improved Lithium Deposition and Stripping

Lithium metal batteries have been plagued by the high reactivity of lithium. Reactive additives that can passivate the lithium metal surface and limit electrolyte accessibility to a fresh lithium surface have been widely explored, but can have limited utility with continuous consumption of the additive. In this work, an alternative strategy is explored. The use of nonreactive cosolvents such as nonpolar alkanes is studied and its is shown that hexane and cyclohexane addition to ether solvents (1,3‐dioxolane and 1,2‐dimethoxyethane) halves the nucleation and growth overpotentials for lithium deposition, increases the cell coulombic efficiency, improves the lithium deposition morphology, increases the electrolyte oxidative stability (>0.2 V), and doubles the cycle life—even when compared to a widely used fluorinated ether. The nonpolar alkanes modify the lithium‐ion solvation environment and reduce the solvation free energy; hence reducing the reaction barrier for lithium deposition. Exploration of nonpolar alkanes as part of the electrolyte mixture is a promising strategy for controlling metal deposition.     

Analysis on long‐range residue–residue communication using molecular dynamics

We investigated the possibility of inter‐residue communication of side chains in barstar, an 89 residue protein, using mutual information theory. The normalized mutual information (NMI) of the dihedral angles of the side chains was obtained from all‐atom molecular dynamics simulations. The accumulated NMI from an explicit solvent equilibrated trajectory (600 ns) with free backbone exhibits a parabola‐shaped distribution over the inter‐residue distances (0–36 Å): smaller at the end regimes but larger in the middle regime. This analysis, plus several other measures, does not find unusual long‐range communication for free backbone in explicit solvent simulations. Proteins 2014; 82:2896–2901. © 2014 Wiley Periodicals, Inc.     

Water penetration in the low and high pressure native states of ubiquitin

Theoretical studies on the solvation of methane molecules in water have shown that the effect of increased pressure is to stabilize solvent separated contacts relative to direct contacts. This suggests that high pressure stabilizes waters that have penetrated into a protein's core, indicating a mechanism for the high pressure denaturation of proteins. We test this theory on a folded protein by studying the penetration of water into the native state of ubiquitin at low and high pressures, using molecular dynamics. An ensemble of conformations sampled in the folded state of ubiquitin has been determined by NMR at two pressures below the protein's denaturation pressure, 30 atm and 3000 atm. We find that 1-5 more waters penetrate the high pressure conformations than the low pressure conformations. Low volume configurations of the system are favored at high pressures, but different components of the system may experience increases or decreases in their specific volumes. We find that penetrating waters have a higher volume per water than bulk waters, but that the volume per protein residue may be lowered by solvation. Furthermore, we find that penetration of the protein by water at high pressures is driven by the difference in the pressure dependence of the probability of cavity opening in the protein and pressure dependence of the probability of cavity opening in the bulk solvent. The volume changes associated with cavity opening and closing indicate that each penetrating water reduces the volume of the system by about 12 mL/mol. The experimental volume change going from the low pressure to the high pressure native state of ubiquitin is 24 mL/mol. Our results indicate that this volume change can be explained by penetration of the protein by two water molecules.     

α‐chymotrypsin in water–acetone and water–dimethyl sulfoxide mixtures: Effect of preferential solvation and hydration

We investigated water/organic solvent sorption and residual enzyme activity to simultaneously monitor preferential solvation/hydration of protein macromolecules in the entire range of water content at 25°C. We applied this approach to estimate protein destabilization/stabilization due to the preferential interactions of bovine pancreatic α‐chymotrypsin with water‐acetone (moderate‐strength H‐bond acceptor) and water‐DMSO (strong H‐bond acceptor) mixtures. There are three concentration regimes for the dried α‐chymotrypsin. α‐Chymotrypsin is preferentially hydrated at high water content. The residual enzyme activity values are close to 100%. At intermediate water content, the dehydrated α‐chymotrypsin has a higher affinity for acetone/DMSO than for water. Residual enzyme activity is minimal in this concentration range. The acetone/DMSO molecules are preferentially excluded from the protein surface at the lowest water content, resulting in preferential hydration. The residual catalytic activity in the water‐poor acetone is ～80%, compared with that observed after incubation in pure water. This effect is very small for the water‐poor DMSO. Two different schemes are operative for the hydrated enzyme. At high and intermediate water content, α‐chymotrypsin exhibits preferential hydration. However, at intermediate water content, in contrast to the dried enzyme, the initially hydrated α‐chymotrypsin possesses increased preferential hydration parameters. At low water content, no residual enzyme activity was observed. Preferential binding of DMSO/acetone to α‐chymotrypsin was detected. Our data clearly demonstrate that the hydrogen bond accepting ability of organic solvents and the protein hydration level constitute key factors in determining the stability of protein–water–organic solvent systems.     

A theoretical study on the optical absorption of green fluorescent protein chromophore in solutions

Ethyl 4-(4-hydroxyphenyl) methylidene- 2-methyl-5-oxoimidazolacetate (HBMIA) is a model chromophore of green fluorescent protein. The electronic structure of HBMIA in aqueous solution phase is studied using a hybrid method of quantum chemistry and statistical mechanics, RISM-SCF-SEDD. The solvatochromic shift is correctly reproduced by the present computations.