How adding a single methylene to dihydrofolate reductase can change its conformational dynamics.

Studies of the effects of pressure on proteins from piezophilic (pressure-loving) microbes compared with homologous proteins from mesophilic microbes have been relatively rare. Interestingly, such studies of dihydrofolate reductase show that a single-site mutation from an aspartic acid to a glutamic acid can reverse the pressure-dependent monotonic decrease in activity to that in a monotonic pressure-dependent activation. This residue is near the active site but is not thought to directly participate in the catalytic mechanism. Here, the ways that addition of one carbon to the entire protein could lead to such a profound difference in pressure effects are explored using molecular dynamics simulations. The results indicate that the glutamate changes the coupling between a helix and the ß-sheet due to the extra flexibility of the side chain, which further changes correlated motions of other regions of the protein.

Enzymes from piezophiles.

The discovery of microbial communities in extreme conditions that would seem hostile to life leads to the question of how the molecules making up these microbes can maintain their structure and function. While microbes that live under extremes of temperature have been heavily studied, those that live under extremes of pressure, or "piezophiles", are now increasingly being studied because of advances in sample collection and high-pressure cells for biochemical and biophysical measurements. Here, adaptations of enzymes in piezophiles against the effects of pressure are discussed in light of recent experimental and computational studies. However, while concepts from studies of enzymes from temperature extremophiles can provide frameworks for understanding adaptations by piezophile enzymes, the effects of temperature and pressure on proteins differ in significant ways. Thus, the state of the knowledge of adaptation in piezophile enzymes is still in its infancy and many more experiments and computational studies on different enzymes from a variety of piezophiles are needed.

Understanding how water models affect the anomalous pressure dependence of their diffusion coefficients.

The self-diffusion coefficient of water shows an anomalous increase with increasing hydrostatic pressure up to a broad maximum (PmD) near 1 kbar at 298 K, which has been attributed to pressure effects on the tetrahedral hydrogen bond network of water. Moreover, the ability of a water model to reproduce anomalous properties of water is a signature that it is reproducing the network. Here, water was simulated between 1 bar and 5 kbar using three water models, two four-site (with all charges in the molecular plane) and one single-site multipole (which accounts for out-of-molecular plane charge), that have reasonable pressure-temperature properties. For these three models, the diffusion coefficients display a maximum in the pressure dependence and the radial distribution functions show good agreement with the limited experimental structural data at high pressure that are available. In addition, a variety of properties associated with the network are examined, including hydrogen bond lifetimes and occupancies, three-body angle distributions, and tetrahedral order parameters. Results suggest that the initial increasing diffusion with pressure is because hydrogen bonds are distorted and thus weakened by pressure, but above PmD, the hydrogen bonds are weakened to the point it behaves more like a normal liquid. In other words, the PmD may be a measure of the angular strength of hydrogen bonds. In addition, since the four-site models over-predict the values of PmD while the multipole model under-predicts it, out-of-plane charge may improve four-site models.

Reduction potential calculations of the Fe-S clusters in Thermus thermophilus respiratory complex I.

Respiratory complex I facilitates electron transfer from NADH to quinone over ~95 Å through a chain of seven iron-sulfur (Fe-S) clusters in the respiratory chain. In this study, the reduction potentials of the Fe-S clusters in Thermus thermophilus complex I are calculated using a Density Functional Theory + Poisson-Boltzmann method. Our results indicate that the reduction potentials are influenced by a variety of factors including the clusters being deeply buried in the complex and the protonation state of buried ionizable residues. In addition, as several of the ionizable side chains have predicted pKa values near pH 7, relatively small structural fluctuations could lead to significant (0.2 V) shifts in the reduction potential of several of the Fe-S clusters, suggesting a dynamic mechanism for electron transfer. Moreover, the method used here is a useful computational tool to study other questions about complex I. © 2019 Wiley Periodicals, Inc.

Effects of Pressure and Temperature on the Atomic Fluctuations of Dihydrofolate Reductase from a Psychropiezophile and a Mesophile.

Determining the effects of extreme conditions on proteins from "extremophilic" and mesophilic microbes is important for understanding how life adapts to living at extremes as well as how extreme conditions can be used for sterilization and food preservation. Previous molecular dynamics simulations of dihydrofolate reductase (DHFR) from a psychropiezophile (cold- and pressure-loving), Moritella profunda (Mp), and a mesophile, Escherichia coli (Ec), at various pressures and temperatures indicate that atomic fluctuations, which are important for enzyme function, increase with both temperature and pressure. Here, the factors that cause increases in atomic fluctuations in the simulations are examined. The fluctuations increase with temperature not only because of greater thermal energy and thermal expansion of the protein but also because hydrogen bonds between protein atoms are weakened. However, the increase in fluctuations with pressure cannot be due to thermal energy, which remains constant, nor the compressive effects of pressure, but instead, the hydrogen bonds are also weakened. In addition, increased temperature causes larger increases in fluctuations of the loop regions of MpDHFR than EcDHFR, and increased pressure causes both increases and decreases in fluctuations of the loops, which differ between the two.

Diffusion of aqueous solutions of ionic, zwitterionic, and polar solutes.

The properties of aqueous solutions of ionic, zwitterionic, and polar solutes are of interest to many fields. For instance, one of the many anomalous properties of aqueous solutions is the behavior of water diffusion in different monovalent salt solutions. In addition, solutes can affect the stabilities of macromolecules such as proteins in aqueous solution. Here, the diffusivities of aqueous solutions of sodium chloride, potassium chloride, tri-methylamine oxide (TMAO), urea, and TMAO-urea are examined in molecular dynamics simulations. The decrease in the diffusivity of water with the concentration of simple ions and urea can be described by a simple model in which the water molecules hydrogen bonded to the solutes are considered to diffuse at the same rate as the solutes, while the remainder of the water molecules are considered to be bulk and diffuse at almost the same rate as pure water. On the other hand, the decrease in the diffusivity of water with the concentration of TMAO is apparently affected by a decrease in the diffusion rate of the bulk water molecules in addition to the decrease due to the water molecules hydrogen bonded to TMAO. In other words, TMAO enhances the viscosity of water, while urea barely affects it. Overall, this separation of water molecules into those that are hydrogen bonded to solute and those that are bulk can provide a useful means of understanding the short- and long-range effects of solutes on water.

Extreme biophysics: Enzymes under pressure.

A critical question about piezophilic (pressure-loving) microbes is how their constituent molecules maintain function under high pressure. Here, factors are examined that may lead to the increased activity under pressure in dihydrofolate reductase from the piezophilic Moritella profunda compared to the homologous enzyme from the mesophilic Escherichia coli. Molecular dynamics simulations are performed at various temperatures and pressures to examine how pressure affects the flexibility of the enzymes from these two microbes, since both stability and flexibility are necessary for enzyme activity. The results suggest that collective motions on the 10-ns timescale are responsible for the flexibility necessary for "corresponding states" activity at the growth conditions of the parent organism. In addition, the results suggest that while the lower stability of many enzymes from deep-sea microbes may be an adaptation for greater flexibility at low temperatures, high pressure may enhance their adaptation to low temperatures. © 2017 Wiley Periodicals, Inc.

Building better water models using the shape of the charge distribution of a water molecule.

The unique properties of liquid water apparently arise from more than just the tetrahedral bond angle between the nuclei of a water molecule since simple three-site models of water are poor at mimicking these properties in computer simulations. Four- and five-site models add partial charges on dummy sites and are better at modeling these properties, which suggests that the shape of charge distribution is important. Since a multipole expansion of the electrostatic potential describes a charge distribution in an orthogonal basis set that is exact in the limit of infinite order, multipoles may be an even better way to model the charge distribution. In particular, molecular multipoles up to the octupole centered on the oxygen appear to describe the electrostatic potential from electronic structure calculations better than four- and five-site models, and molecular multipole models give better agreement with the temperature and pressure dependence of many liquid state properties of water while retaining the computational efficiency of three-site models. Here, the influence of the shape of the molecular charge distribution on liquid state properties is examined by correlating multipoles of non-polarizable water models with their liquid state properties in computer simulations. This will aid in the development of accurate water models for classical simulations as well as in determining the accuracy needed in quantum mechanical/molecular mechanical studies and ab initio molecular dynamics simulations of water. More fundamentally, this will lead to a greater understanding of how the charge distribution of a water molecule leads to the unique properties of liquid water. In particular, these studies indicate that p-orbital charge out of the molecular plane is important.

Quasiharmonic analysis of protein energy landscapes from pressure-temperature molecular dynamics simulations.

Positional fluctuations of an atom in a protein can be described as motion in an effective local energy minimum created by the surrounding protein atoms. The dependence of atomic fluctuations on both temperature (T) and pressure (P) has been used to probe the nature of these minima, which are generally described as harmonic in experiments such as x-ray crystallography and neutron scattering. Here, a quasiharmonic analysis method is presented in which the P-T dependence of atomic fluctuations is in terms of an intrinsic isobaric thermal expansivity αP and an intrinsic isothermal compressibility κT. The method is tested on previously reported mean-square displacements from P-T molecular dynamics simulations of lysozyme, which were interpreted to have a pressure-independent dynamical transition Tg near 200 K and a change in the pressure dependence near 480 MPa. Our quasiharmonic analysis of the same data shows that the P-T dependence can be described in terms of αP and κT where below Tg, the temperature dependence is frozen at the Tg value. In addition, the purported transition at 480 MPa is reinterpreted as a consequence of the pressure dependence of Tg and the quasiharmonic frequencies. The former also indicates that barrier heights between substates are pressure dependent in these data. Furthermore, the insights gained from this quasiharmonic analysis, which was of the energy landscape near the native state of a protein, suggest that similar analyses of other simulations may be useful in understanding such phenomena as pressure-induced protein unfolding.

What makes proteins work: exploring life in P-T-X.

Although considerable progress has been made in the molecular biophysics of proteins, it is still not possible to reliably design an enzyme for a given function. The current understanding of enzyme function is that both structure and flexibility are important. Much attention has been focused recently on protein folding and thus structure, spurred on by insights from the folding funnel concept. For experimental studies of protein folding, variations in temperature (T) and chemical composition (X) of the solution have been traditionally exploited, although more recent studies using variations in pressure (P) made possible through new instrumentation have led to a deeper understanding of the energy landscape of protein folding. Other work has shown that flexibility is also essential for enzymes, although it is still not clear what type is important. Another avenue has been to take advantage of 'Nature's laboratory' by exploring homologous proteins from organisms that live in extreme conditions, or 'extremophiles'. While the most studied extremophiles live at extremes of T and X, recent exploration of deep-sea environments has led to the discovery of organisms living under high P, or 'piezophiles'. An exploration of targeted enzymes from organisms with various P-T-X growth conditions coupled with advances in biophysical instrumentation and computer simulations that allow studies of these enzymes at different P-T-X conditions may lead to a better understanding of 'flexibility' and to general design criteria for active enzymes. Preface. Kamal Shukla's great contribution to science has been his vision that physical sciences could bring new insights to biological sciences, and that the marriage of methodologies, particularly theoretical/computational with experimental, was needed to tackle the complexities of biology. Furthermore, his openness to new methods and different ideas outside the current fad has helped make his vision a reality. In my remarks below, I have not tried to limit myself to projects that I know Kamal had sponsored, nor have I tried to highlight all that he has sponsored. Instead, everything I mention has been influenced directly or indirectly by his efforts. Perhaps the indirect influences are most telling, because they would not have happened without Kamal.

A single-site multipole model for liquid water.

Accurate and efficient empirical potential energy models that describe the atomistic interactions between water molecules in the liquid phase are essential for computer simulations of many problems in physics, chemistry, and biology, especially when long length or time scales are important. However, while models with non-polarizable partial charges at four or five sites in a water molecule give remarkably good values for certain properties, deficiencies have been noted in other properties and increasing the number of sites decreases computational efficiency. An alternate approach is to utilize a multipole expansion of the electrostatic potential due to the molecular charge distribution, which is exact outside the charge distribution in the limits of infinite distances or infinite orders of multipoles while partial charges are a qualitative representation of electron density as point charges. Here, a single-site multipole model of water is presented, which is as fast computationally as three-site models but is also more accurate than four- and five-site models. The dipole, quadrupole, and octupole moments are from quantum mechanical-molecular mechanical calculations so that they account for the average polarization in the liquid phase, and represent both the in-plane and out-of-plane electrostatic potentials of a water molecule in the liquid phase. This model gives accurate thermodynamic, dynamic, and dielectric properties at 298 K and 1 atm, as well as good temperature and pressure dependence of these properties.

Web-based computational chemistry education with CHARMMing III: Reduction potentials of electron transfer proteins.

A module for fast determination of reduction potentials, E°, of redox-active proteins has been implemented in the CHARMM INterface and Graphics (CHARMMing) web portal (www.charmming.org). The free energy of reduction, which is proportional to E°, is composed of an intrinsic contribution due to the redox site and an environmental contribution due to the protein and solvent. Here, the intrinsic contribution is selected from a library of pre-calculated density functional theory values for each type of redox site and redox couple, while the environmental contribution is calculated from a crystal structure of the protein using Poisson-Boltzmann continuum electrostatics. An accompanying lesson demonstrates a calculation of E°. In this lesson, an ionizable residue in a [4Fe-4S]-protein that causes a pH-dependent E° is identified, and the E° of a mutant that would test the identification is predicted. This demonstration is valuable to both computational chemistry students and researchers interested in predicting sequence determinants of E° for mutagenesis.

Web-based computational chemistry education with CHARMMing I: Lessons and tutorial.

This article describes the development, implementation, and use of web-based "lessons" to introduce students and other newcomers to computer simulations of biological macromolecules. These lessons, i.e., interactive step-by-step instructions for performing common molecular simulation tasks, are integrated into the collaboratively developed CHARMM INterface and Graphics (CHARMMing) web user interface (http://www.charmming.org). Several lessons have already been developed with new ones easily added via a provided Python script. In addition to CHARMMing's new lessons functionality, web-based graphical capabilities have been overhauled and are fully compatible with modern mobile web browsers (e.g., phones and tablets), allowing easy integration of these advanced simulation techniques into coursework. Finally, one of the primary objections to web-based systems like CHARMMing has been that "point and click" simulation set-up does little to teach the user about the underlying physics, biology, and computational methods being applied. In response to this criticism, we have developed a freely available tutorial to bridge the gap between graphical simulation setup and the technical knowledge necessary to perform simulations without user interface assistance.

Hydrophobic hydration and the anomalous partial molar volumes in ethanol-water mixtures.

The anomalous behavior in the partial molar volumes of ethanol-water mixtures at low concentrations of ethanol is studied using molecular dynamics simulations. Previous work indicates that the striking minimum in the partial molar volume of ethanol VE as a function of ethanol mole fraction XE is determined mainly by water-water interactions. These results were based on simulations that used one water model for the solute-water interactions but two different water models for the water-water interactions. This is confirmed here by using two more water models for the water-water interactions. Furthermore, the previous work indicates that the initial decrease is caused by association of the hydration shells of the hydrocarbon tails, and the minimum occurs at the concentration where all of the hydration shells are touching each other. Thus, the characteristics of the hydration of the tail that cause the decrease and the features of the water models that reproduce this type of hydration are also examined here. The results show that a single-site multipole water model with a charge distribution that mimics the large quadrupole and the p-orbital type electron density out of the molecular plane has "brittle" hydration with hydrogen bonds that break as the tails touch, which reproduces the deep minimum. However, water models with more typical site representations with partial charges lead to flexible hydration that tends to stay intact, which produces a shallow minimum. Thus, brittle hydration may play an essential role in hydrophobic association in water.

The molecular charge distribution, the hydration shell, and the unique properties of liquid water.

The most essential features of a water molecule that give rise to its unique properties are examined using computer simulations of different water models. The charge distribution of a water molecule characterized by molecular multipoles is quantitatively linked to the liquid properties of water via order parameters for the degree (S(2)) and symmetry (ΔS(2)) of the tetrahedral arrangement of the nearest neighbors, or "hydration shell." ΔS(2) also appears to determine the long-range tetrahedral network and interfacial structure. From the correlations, some models are shown to be unable to reproduce certain properties due to the limitations of the model itself rather than the parameterization, which indicates that they are lacking essential molecular features. Moreover, since these properties depend not only on S(2) but also on ΔS(2), the long-range structure in these models may be incorrect. Based on the molecular features found in the models that are best able to reproduce liquid properties, the most essential features of a water molecule in liquid water appear to be a charge distribution with a large dipole, a large quadrupole, and negative charge out of the molecular plane, as well as a symmetrically ordered tetrahedral hydration shell that results from this charge distribution. The implications for modeling water are also discussed.

Identifying residues that cause pH-dependent reduction potentials.

The pH dependence of the reduction potential E° for a metalloprotein indicates that the protonation state of at least one residue near the redox site changes and may be important for its activity. The responsible residue is usually identified by site-specific mutagenesis, which may be time-consuming. Here, the titration of E° for Chromatium vinosum high-potential iron-sulfur protein is predicted to be in good agreement with experiment using density functional theory and Poisson-Boltzmann calculations if only the sole histidine undergoes changes in protonation. The implementation of this approach into CHARMMing, a user-friendly web-based portal, allows users to identify residues in other proteins causing similar pH dependence.

Effects of microcomplexity on hydrophobic hydration in amphiphiles.

Hydrophobic hydration is critical in biology as well as many industrial processes. Here, computer simulations of ethanol/water mixtures show that a three-stage mechanism of dehydration of ethanol explains the anomalous concentration dependence of the thermodynamic partial molar volumes, as well as recent data from neutron diffraction and Raman scattering. Moreover, the simulations show that the breakdown of hydrophobic hydration shells, whose structure is determined by the unique molecular properties of water, is caused by the microcomplexity of the environment and may be representative of early events in protein folding and structure stabilization in aqueous solutions.

Calculating standard reduction potentials of [4Fe-4S] proteins.

The oxidation-reduction potentials of electron transfer proteins determine the driving forces for their electron transfer reactions. Although the type of redox site determines the intrinsic energy required to add or remove an electron, the electrostatic interaction energy between the redox site and its surrounding environment can greatly shift the redox potentials. Here, a method for calculating the reduction potential versus the standard hydrogen electrode, E°, of a metalloprotein using a combination of density functional theory and continuum electrostatics is presented. This work focuses on the methodology for the continuum electrostatics calculations, including various factors that may affect the accuracy. The calculations are demonstrated using crystal structures of six homologous HiPIPs, which give E° that are in excellent agreement with experimental results.

Characterizing the effects of the protein environment on the reduction potentials of metalloproteins.

The reduction potentials of electron transfer proteins are critically determined by the degree of burial of the redox site within the protein and the degree of permanent polarization of the polypeptide around the redox site. Although continuum electrostatics calculations of protein structures can predict the net effect of these factors, quantifying each individual contribution is a difficult task. Here, the burial of the redox site is characterized by a dielectric radius R(p) (a Born-type radius for the protein), the polarization of the polypeptide is characterized by an electret potential ϕ(p) (the average electrostatic potential at the metal atoms), and an electret-dielectric spheres (EDS) model of the entire protein is then defined in terms of R(p) and ϕ(p). The EDS model shows that for a protein with a redox site of charge Q, the dielectric response free energy is a function of Q(2), while the electret energy is a function of Q. In addition, R(p) and ϕ(p) are shown to be characteristics of the fold of a protein and are predictive of the most likely redox couple for redox sites that undergo different redox couples.

Identifying sequence determinants of reduction potentials of metalloproteins.

The reduction potential of an electron transfer protein is one of its most important functional characteristics. Although the type of redox site and the protein fold are the major determinants of the reduction potential of a redox-active protein, its amino acid sequence may tune the reduction potential as well. Thus, homologous proteins can often be divided into different classes, with each class characterized by a biological function and a reduction potential. Site-specific mutagenesis of the sequence determinants of the differences in the reduction potential between classes should change the reduction potential of a protein in one class to that of the other class. Here, a procedure is presented that combines energetic and bioinformatic analysis of homologous proteins to identify sequence determinants that are also good candidates for site-specific mutations, using the [4Fe-4S] ferredoxins and the [4Fe-4S] high-potential iron-sulfur proteins as examples. This procedure is designed to guide site-specific mutations or more computationally expensive studies, such as molecular dynamics simulations. To make the procedure more accessible to the general scientific community, it is being implemented into CHARMMing, a Web-based portal, with a library of density functional theory results for the redox site that are used in the setting up of Poisson-Boltzmann continuum electrostatics calculations for the protein energetics.