Water Determines the Structure and Dynamics of Proteins

Water-Triggered, Irreversible Conformational Change of SARS-CoV-2 Main Protease on Passing from the Solid State to Aqueous Solution

The main protease from SARS-CoV-2 is a homodimer. Yet, a recent 0.1-ms-long molecular dynamics simulation performed by D. E. Shaw's research group shows that it readily undergoes a symmetry-breaking event on passing from the solid state to aqueous solution. As a result, the subunits present distinct conformations of the binding pocket. By analyzing this long simulation, we uncover a previously unrecognized role of water molecules in triggering the transition. Interestingly, each subunit presents a different collection of long-lived water molecules. Enhanced sampling simulations performed here, along with machine learning approaches, further establish that the transition to the asymmetric state is essentially irreversible.

Water molecules at protein-drug interfaces: computational prediction and analysis methods

The fundamental importance of water molecules at drug-protein interfaces is now widely recognised and a significant feature in structure-based drug design. Experimental methods for analysing the role of water in drug binding have many challenges, including the accurate location of bound water molecules in crystal structures, and problems in resolving specific water contributions to binding thermodynamics. Computational analyses of binding site water molecules provide an alternative, and in principle complete, structural and thermodynamic picture, and their use is now commonplace in the pharmaceutical industry. In this review, we describe the computational methodologies that are available and discuss their strengths and weaknesses. Additionally, we provide a critical analysis of the experimental data used to validate the methods, regarding the type and quality of experimental structural data. We also discuss some of the fundamental difficulties of each method and suggest directions for future study.

Hydration of a small protein under carbon nanotube confinement: Adsorbed substates induce selective separation of the dynamical response

The co-involvement of biological molecules and nanomaterials has increasingly come to the fore in modern-day applications. While the "bio-nano" (BN) interface presents physico-chemical characteristics that are manifestly different from those observed in isotropic bulk conditions, the underlying molecular reasons remain little understood; this is especially true of anomalies in interfacial hydration. In this paper, we leverage atomistic simulations to study differential adsorption characteristics of a small protein on the inner (concave) surface of a single-walled carbon nanotube whose diameter exceeds dimensions conducive to single-file water movement. Our findings indicate that the extent of adsorption is decided by the degree of foldedness of the protein conformational substate. Importantly, we find that partially folded substates, but not the natively folded one, induce reorganization of the protein hydration layer into an inner layer water closer to the nanotube axis and an outer layer water in the interstitial space near the nanotube walls. Further analyses reveal sharp dynamical differences between water molecules in the two layers as observed in the onset of increased heterogeneity in rotational relaxation and the enhanced deviation from Fickian behavior. The vibrational density of states reveals that the dynamical distinctions are correlated with differences in crucial bands in the power spectra. The current results set the stage for further systematic studies of various BN interfaces vis-à-vis control of hydration properties.

Structural and dynamical heterogeneity of water trapped inside Na+-pumping KR2 rhodopsin in the dark state

Photoisomerization in the retinal leads to a channel opening in rhodopsins that triggers translocation or pumping of ions/protons. Crystal structures of rhodopsins contain several structurally conserved water molecules. It has been suggested that water plays an active role in facilitating the ion pumping/translocation process by acting as a lubricant in these systems. In this paper, we systematically investigate the localization, structure, dynamics, and energetics of the water molecules along the channel for the resting/dark state of KR2 rhodopsin. By employing several microseconds long atomistic molecular dynamics simulation of this trans-membrane protein system, we demonstrate the presence of five distinct water containing pockets/cavities separated by gateways controlled by protein side-chains. There exists a strong hydrogen bonded network involving these buried water molecules and functionally important key residues. We present evidence of significant structural and dynamical heterogeneity in the water molecules present in these cavities, with very rare exchange between them. The exchange time scale of such buried water with the bulk has an extremely wide range, from tens of nanoseconds to >1.5 µs. The translational and rotational dynamics of buried water are found to be strongly dependent on the protein cavity size and local interactions with a classic signature of trapped diffusion and rotational anisotropy.

Molecular dynamics evidence for nonthermal effects of electric fields on pectin methylesterase activity

Experimental studies relevant to the nonthermal effects of electric fields on biological systems are emerging. However, these effects are poorly understood at the molecular level. The present study investigates pectin methylesterase, a cell wall modifying enzyme in plants, exposed to various electric field strengths. Molecular dynamics (MD) of the enzyme were studied with and without (thermal-only) electric field applications. The measurements were interpreted on the basis of equivalent energy input to gain insights into the effect of electric field treatment time at a constant temperature (50 °C). Results reveal that electric fields exert nonthermal effects on both local and global protein structure. In 1 μs simulations, the results show significant (P ≤ 0.05) shrinkage of the catalytic domain and shortening of enzyme-water hydrogen bond lifetime by a 50 V cm-1 electric field. Unwinding of the helical segments, altered intra- and intermolecular hydrogen bond patterns, and increased hydration are also caused by the 50 V cm-1 electric field. This study serves to understand the electric field influence on the functional role of proteins.

Abnormal Properties of Low-Dimensional Confined Water

Water molecules confined to low-dimensional spaces exhibit unusual properties compared to bulk water. For example, the alternating hydrophilic and hydrophobic nanodomains on flat silicon wafer can induce the abnormal spreading of water (contact angles near 0°) which is caused by the 2D capillary effect. Hence, exploring the physicochemical properties of confined water from the nanoscale is of great value for understanding the challenges in material science and promoting the applications of nanomaterials in the fields of mass transport, nanofluidic designing, and fuel cell. The knowledge framework of confined water can also help to better understand the complex functions of the hydration layer of biomolecules, and even trace the origin of life. In this review, the physical properties, abnormal behaviors, and functions of the confined water are mainly summarized through several common low-dimensional water formats in the fields of solid/air-water interface, nanochannel confinement, and biological hydration layer. These researches indicate that the unusual behaviors of the confined water depend strongly on the confinement size and the interaction between the molecules and confining surface. These diverse properties of confined water open a new door to materials science and may play an important role in the future development of biology.

Age-Dependent Changes in Total and Free Water Content of In Vivo Human Lenses Measured by Magnetic Resonance Imaging

Purpose

To use magnetic resonance imaging (MRI) to measure age-dependent changes in total and free water in human lenses in vivo.

Methods

Sixty-four healthy adults aged 18 to 86 years were recruited, fitted with a 32-channel head receiver coil, and placed in a 3 Tesla clinical MR scanner. Scans of the crystalline lens were obtained using a volumetric interpolated breath-hold examination sequence with dual flip angles, which were corrected for field inhomogeneity post-acquisition using a B1-map obtained using a turbo-FLASH sequence. The spatial distribution and content of corrected total (ρlens) and free (T1) water along the lens optical axis were extracted using custom-written code.

Results

Lens total water distribution and content did not change with age (all P > 0.05). In contrast to total water, a gradient in free water content that was highest in the periphery relative to the center was present in lenses across all ages. However, this initially parabolic free water gradient gradually developed an enhanced central plateau, as indicated by increasing profile shape parameter values (anterior: 0.067/y, P = 0.004; posterior: 0.050/y, P = 0.020) and central free water content (1.932 ms/y, P = 0.022) with age.

Conclusions

MRI can obtain repeatable total and free water measurements of in vivo human lenses. The observation that the lens steady-state free, but not total, water gradient is abolished with age raises the possibility that alterations in protein-water interactions are an underlying cause of the degradation in lens optics and overall vision observed with aging.

Dynamics of aqueous peptide solutions in folded and disordered states examined by dynamic light scattering and dielectric spectroscopy

Characterizing the segmental dynamics of proteins, and intrinsically disordered proteins in particular, is a challenge in biophysics. In this study, by combining data from broadband dielectric spectroscopy (BDS) and both depolarized (DDLS) and polarized (PDLS) dynamic light scattering, we were able to determine the dynamics of a small peptide [ε-poly(lysine)] in water solutions in two different conformations (pure β-sheet at pH = 10 and a more disordered conformation at pH = 7). We found that the segmental (α-) relaxation, as probed by DDLS, is faster in the disordered state than in the folded conformation. The water dynamics, as detected by BDS, is also faster in the disordered state. In addition, the combination of BDS and DDLS results allows us to confirm the molecular origin of water-related processes observed by BDS. Finally, we discuss the origin of two slow processes (A and B processes) detected by DDLS and PDLS in both conformations and usually observed in other types of water solutions. For fully homogeneous ε-PLL solutions at pH = 10, the A-DLS process is assigned to the diffusion of individual β-sheets. The combination of both techniques opens a route for understanding the dynamics of peptides and other biological solutions.

Spectrally Resolved Estimation of Water Entropy in the Active Site of Human Carbonic Anhydrase II

A major challenge in understanding ligand binding to biomacromolecules lies in dissecting the underlying thermodynamic driving forces at the atomic level. Quantifying the contributions of water molecules is often especially demanding, although they can play important roles in biomolecular recognition and binding processes. One example is human carbonic anhydrase II, whose active site harbors a conserved network of structural water molecules that are essential for enzymatic catalysis. Inhibitor binding disrupts this water network and changes the hydrogen-bonding patterns in the active site. Here, we use atomistic molecular dynamics simulations to compute the absolute entropy of the individual water molecules confined in the active site of hCAII using a spectrally resolved estimation (SRE) approach. The entropy decrease of water molecules that remain in the active site upon binding of a dorzolamide inhibitor is caused by changes in hydrogen bonding and stiffening of the hydrogen-bonding network. Overall, this entropy decrease is overcompensated by the gain due to the release of three water molecules from the active site upon inhibitor binding. The spectral density calculations enable the assignment of the changes to certain vibrational modes. In addition, the range of applicability of the SRE approximation is systematically explored by exploiting the gradually changing degree of immobilization of water molecules as a function of the distance to a phospholipid bilayer surface, which defines an "entropy ruler". These results demonstrate the applicability of SRE to biomolecular solvation, and we expect it to become a useful method for entropy calculations in biomolecular systems.

Unfolding of Dynamical Events in the Early Stage of Insulin Dimer Dissociation

The dissociation of an insulin dimer is an important biochemical event that could also serve as a prototype of dissociations in similar biomolecular assemblies. We use a recently developed multidimensional free energy landscape for insulin dimer dissociation to unearth the microscopic and mechanistic aspects of the initial stages of the process that could hold the key to understanding the stability and the rate. The following sequence of events occurs in the initial stages: (i) The backbone hydrogen bonds break partially at the antiparallel β-sheet junction, (ii) the two α-helices (chain B) move away from each other while several residues (chain A) move closer, and (iii) a flow of adjacent water molecules occurs into the junction region. Interestingly, the intermonomeric center-to-center distance does not increase, but the number of native contacts exhibits a sharp decrease. Subsequent steps involve further disengagement of hydrophobic groups. This process is slow because of an entropic bottleneck created by the existence of the large configuration space available in the native state (NS), which is inhabited by low-frequency conformational fluctuations. We carry out a density-of-states analyses in the dimer NS to unearth distinctive features not present in the monomers. These low-frequency modes are also responsible for a large entropic stabilization of the NS. Hydrophobic disengagement in the early stage leads to the formation of a twisted intermediate state which itself is a metastable minimum (IS-1). The subsequent progress leads to another dimeric complex (IS-2), which is on the dissociative pathway and characterized by a further decrease in the native contacts. The dissociation process provides insights into the workings of a biomolecular assembly.

Interfacial Water Many-Body Effects Drive Structural Dynamics and Allosteric Interactions in SARS-CoV-2 Main Protease Dimerization Interface

Following our previous work ( Chem. Sci. 2021, 12, 4889-4907), we study the structural dynamics of the SARS-CoV-2 Main Protease dimerization interface (apo dimer) by means of microsecond adaptive sampling molecular dynamics simulations (50 μs) using the AMOEBA polarizable force field (PFF). This interface is structured by a complex H-bond network that is stable only at physiological pH. Structural correlations analysis between its residues and the catalytic site confirms the presence of a buried allosteric site. However, noticeable differences in allosteric connectivity are observed between PFFs and non-PFFs. Interfacial polarizable water molecules are shown to appear at the heart of this discrepancy because they are connected to the global interface H-bond network and able to adapt their dipole moment (and dynamics) to their diverse local physicochemical microenvironments. The water-interface many-body interactions appear to drive the interface volume fluctuations and to therefore mediate the allosteric interactions with the catalytic cavity.

Increased xerotolerance of Saccharomyces cerevisiae during an osmotic pressure ramp over several generations

Although mechanisms involved in response of Saccharomyces cerevisiae to osmotic challenge are well described for low and sudden stresses, little is known about how cells respond to a gradual increase of the osmotic pressure (reduced water activity; aw ) over several generations as it could encounter during drying in nature or in food processes. Using glycerol as a stressor, we propagated S. cerevisiae through a ramp of the osmotic pressure (up to high molar concentrations to achieve testing-to-destruction) at the rate of 1.5 MPa day-1 from 1.38 to 58.5 MPa (0.990-0.635 aw ). Cultivability (measured at 1.38 MPa and at the harvest osmotic pressure) and glucose consumption compared with the corresponding sudden stress showed that yeasts were able to grow until about 10.5 MPa (0.926 aw ) and to survive until about 58.5 MPa, whereas glucose consumption occurred until 13.5 MPa (about 0.915 aw ). Nevertheless, the ramp conferred an advantage since yeasts harvested at 10.5 and 34.5 MPa (0.778 aw ) showed a greater cultivability than glycerol-shocked cells after a subsequent shock at 200 MPa (0.234 aw ) for 2 days. FTIR analysis revealed structural changes in wall and proteins in the range 1.38-10.5 MPa, which would be likely to be involved in the resistance at extreme osmotic pressure.

Designing cryo-enzymatic reactions in subzero liquid water by lipidic mesophase nanoconfinement

Cryo-enzymology provides the possibility to develop unconventional biological reactions and detect intermediates in ultrafast enzymatic catalysis processes, but also illuminates the understanding of life principles in extremely cold environments. The scarcity of biological or biomimetic host systems that provide liquid water at subzero temperatures inhibits the prosperity of cryo-enzymology. Here we introduce cryo-enzymatic reactions in subzero water nanoconfined within lipid mesophases formed by conventional lipids. We show that the enzymatic reactions that ensue outperform the homologue catalytic processes run at standard temperatures. We use phytantriol-based lipidic mesophases (LMPs), within which water remains in the liquid state down to -120 °C, and combine crystallization and dynamic studies of the confined water to provide a fundamental understanding of the physical status of water at subzero temperatures, which sets the stage for cryo-enzymatic reactions in these environments. In the model horseradish peroxidase oxidization, the cation free-radical product is stabilized in LMPs at -20 °C, in contrast to the fast-consuming reactions at temperatures above 0 °C. Furthermore, the LMP system also supports the cascade reaction and lipase reaction at subzero temperatures, at which enzymatic reactions with both hydrophilic and hydrophobic substrates are successfully carried out. Our designed LMP system opens access to the nature of confined water in the biomimetic environment and provides a platform for low-temperature biomacromolecule reconstitution and the cryogenic control of enzymatic reactions in bionanotechnology.

Hydrogen-bond network distortion of water in the soft confinement of Nafion membrane

A Compton spectroscopy investigation is carried out in hydrated Nafion membranes, enabling identification of distortions in the hydrogen-bond distribution of the polymer hydrating water by means of the subtle changes reflected by the Compton profiles. Indeed, deformations of the Compton profiles are observed when varying hydration, and two different bonding kinds are associated with the water molecules: at low hydration, water surrounds the sulfonic groups, while on increasing hydration, water molecules occupy the interstitial cavities formed upon swelling of the membrane. The analysis is proposed in terms of averaged OH bond length variation. A sizable contraction of the OH distance is observed at low hydration (∼0.09 Å), while at higher hydration levels, the contraction is smaller (∼0.02 Å) and the OH bond length is closer to bulk water. An evaluation of the electron kinetic energy indicates that the spatial changes associated with the water distribution correspond to a consistent binding energy increase. Distinct temperature dependences of each water population are observed, which can be straightly related to water desorption into ice on cooling below the freezing point.

Slowdown of Water Dynamics from the Top to the Bottom of the GroEL Cavity

The GroE molecular chaperone system is a critical protein machine that assists the folding of substrate proteins in its cavity. Water in the cavity is suspected to play a role in substrate protein folding, but the mechanism is currently unknown. Herein, we report measurements of water dynamics in the equatorial and apical domains of the GroEL cavity in the apo and football states, using site-specific tryptophanyl mutagenesis as an intrinsic optical probe with femtosecond resolution combined with molecular dynamics simulations. We observed clearly different water dynamics in the two domains with a slowdown of the cavity water from the apical to equatorial region in the football state. The results suggest that the GroEL cavity provides a unique water environment that may facilitate substrate protein folding.

Residue interaction networks of K-Ras protein with water molecules identifies the potential role of switch II and P-loop

The mutant K-Ras with aberrant signaling is the primary cause of several cancers. The proposed study investigated the influence of water molecules in K-Ras crystal structure, where they have a significant function by understanding their residue interaction networks (RINs). We analyzed the RINs of K-Ras with and without water molecules and determined their interaction properties. RINs were developed with the help of StructureViz2 and RINspector; further, the changes in K-Ras backbone flexibility were predicted with the DynaMine. We found that the residues K42, I142, and L159 are the hotspots from water, including the K-Ras-GTP complex with the highest residue centrality analysis (RCA) Z-score. The DynaMine prediction calculated the NMR S² value for the frequently mutated positions G12, G13, and Q61 showing a minor shift in flexibility, which make up the P-Loop and switch II of the K-Ras protein. This flexibility shift can account for changes in conformational activity and the protein's GTPase activity, making it difficult to recognize by the effectors and exchange factors. Taken together, our study helps in understanding the functional importance of the water molecules in K-Ras protein and the impact of mutation that modulate the conformational state of the protein.

Ab initio molecular dynamics study on energy relaxation path of hydrogen-bonded OH vibration in bulk water

The vibrational energy relaxation paths of hydrogen-bonded (H-bonded) OH excited in pure water and in isotopically diluted (deuterated) water are elucidated via non-equilibrium ab initio molecular dynamics (NE-AIMD) simulations. The present study extends the previous NE-AIMD simulation for the energy relaxation of an excited free OH vibration at an air/water interface [T. Ishiyama, J. Chem. Phys. 154, 104708 (2021)] to the energy relaxation of an excited H-bonded OH vibration in bulk water. The present simulation shows that the excited OH vibration in pure water dissipates its energy on a timescale of 0.1 ps, whereas that in deuterated water relaxes on a timescale of 0.7 ps, consistent with the experimental observations. To decompose these relaxation energies into the components due to intramolecular and intermolecular couplings, constraints are introduced on the vibrational modes except for the target path in the NE-AIMD simulation. In the case of pure water, 80% of the total relaxation is attributed to the pathway due to the resonant intermolecular OH⋯OH stretch coupling, and the remaining 17% and 3% are attributed to intramolecular couplings with the bend overtone and with the conjugate OH stretch, respectively. This result strongly supports a significant role for the Förster transfer mechanism of pure water due to the intermolecular dipole-dipole interactions. In the case of deuterated water, on the other hand, 36% of the total relaxation is due to the intermolecular stretch coupling, and all the remaining 64% arises from coupling with the intramolecular bend overtone.

Characterizing the Interplay between Polymer Solvation and Conformation

Conformational transitions of flexible molecules, especially those driven by hydrophobic effects, tend to be hindered by desolvation barriers. For such transitions, it is thus important to characterize and understand the interplay between solvation and conformation. Using specialized molecular simulations, here we perform such a characterization for a hydrophobic polymer solvated in water. We find that an external potential, which unfavorably perturbs the polymer hydration waters, can trigger a coil-to-globule or collapse transition, and that the relative stabilities of the collapsed and extended states can be quantified by the strength of the requisite potential. Our results also provide mechanistic insights into the collapse transition, highlighting that the bottleneck to polymer collapse is the formation of a sufficiently large cluster, and the collective dewetting of such a cluster. We also study the collapse of the hydrophobic polymer in octane, a nonpolar solvent, and interestingly, we find that the mechanistic details of the transition are qualitatively similar to that in water.

Solubilization of Trans-Resveratrol in Some Mono-Solvents and Various Propylene Glycol + Water Mixtures

This research deals with the determination of solubility, Hansen solubility parameters, dissolution properties, enthalpy-entropy compensation, and computational modeling of a naturally-derived bioactive compound trans-resveratrol (TRV) in water, methanol, ethanol, n-propanol, n-butanol, propylene glycol (PG), and various PG + water mixtures. The solubility of TRV in six different mono-solvents and various PG + water mixtures was determined at 298.2-318.2 K and 0.1 MPa. The measured experimental solubility values of TRV were regressed using six different computational/theoretical models, including van't Hoff, Apelblat, Buchowski-Ksiazczak λh, Yalkowsly-Roseman, Jouyban-Acree, and van't Hoff-Jouyban-Acree models, with average uncertainties of less than 3.0%. The maxima of TRV solubility in mole fraction was obtained in neat PG (2.62 × 10-2) at 318.2 K. However, the minima of TRV solubility in the mole fraction was recorded in neat water (3.12 × 10-6) at 298.2 K. Thermodynamic calculation of TRV dissolution properties suggested an endothermic and entropy-driven dissolution of TRV in all studied mono-solvents and various PG + water mixtures. Solvation behavior evaluation indicated an enthalpy-driven mechanism as the main mechanism for TRV solvation. Based on these data and observations, PG has been chosen as the best mono-solvent for TRV solubilization.

Molecular Tailoring Approach for the Estimation of Intramolecular Hydrogen Bond Energy

Hydrogen bonds (HBs) play a crucial role in many physicochemical and biological processes. Theoretical methods can reliably estimate the intermolecular HB energies. However, the methods for the quantification of intramolecular HB (IHB) energy available in the literature are mostly empirical or indirect and limited only to evaluating the energy of a single HB. During the past decade, the authors have developed a direct procedure for the IHB energy estimation based on the molecular tailoring approach (MTA), a fragmentation method. This MTA-based method can yield a reliable estimate of individual IHB energy in a system containing multiple H-bonds. After explaining and illustrating the methodology of MTA, we present its use for the IHB energy estimation in molecules and clusters. We also discuss the use of this method by other researchers as a standard, state-of-the-art method for estimating IHB energy as well as those of other noncovalent interactions.

MHCBI: a pipeline for calculating peptide-MHC binding energy using semi-empirical quantum mechanical methods with explicit/implicit solvent models

Experimentally estimating peptide-major histocompatibility complex (pMHC) binding affinity has been quite challenging due to the many receptors and the many potential ligands implicated in it. We have thus proposed a straightforward computational methodology considering the different mechanisms involved in pMHC binding to facilitate studying such receptor-ligand interactions. We have developed a pipeline using semi-empirical quantum mechanical methods for calculating pMHC class I and II molecules' binding energy (BE). This pipeline can systematize the methodology for calculating pMHC system BE, enabling the rational design of T-cell epitopes to be used as pharmaceuticals and vaccines.

Water-oriented magnetic anisotropy transition

Water reorientation is essential in a wide range of chemical and biological processes. However, the effects of such reorientation through rotation around the metal–oxygen bond on the chemical and physical properties of the resulting complex are usually ignored. Most studies focus on the donor property of water as a recognized σ donor-type ligand rather than a participant in the π interaction. Although a theoretical approach to study water-rotation effects on the functionality of a complex has recently been conducted, it has not been experimentally demonstrated. In this study, we determine that the magnetic anisotropy of a Co(II) complex can be effectively controlled by the slight rotation of coordinating water ligands, which is achieved by a two-step structural phase transition. When the water molecule is rotated by 21.2 ± 0.2° around the Co–O bond, the directional magnetic susceptibility of the single crystal changes by approximately 30% along the a-axis due to the rotation of the magnetic anisotropy axis through the modification of the π interaction between cobalt(II) and the water ligand. The theoretical calculations further support the hypothesis that the reorientation of water molecules is a key factor contributing to the magnetic anisotropy transition of this complex.

Little is known about how the orientation of coordinated water molecules affects the magnetic properties of single molecule magnets. Here the authors combine experimental data and theoretical calculations to study how the rotation of water molecules alters the magnetic anisotropy of a pyrazine-based cobalt(II) complex.

Resurrection of Ancestral Malate Dehydrogenases Reveals the Evolutionary History of Halobacterial Proteins : Deciphering Gene Trajectories and Changes in Biochemical Properties

Extreme halophilic Archaea thrive in high salt, where, through proteomic adaptation, they cope with the strong osmolarity and extreme ionic conditions of their environment. In spite of wide fundamental interest, however, studies providing insights into this adaptation are scarce, because of practical difficulties inherent to the purification and characterization of halophilic enzymes. In this work, we describe the evolutionary history of malate dehydrogenases (MalDH) within Halobacteria (a class of the Euryarchaeota phylum). We resurrected nine ancestors along the inferred halobacterial MalDH phylogeny, including the Last Common Ancestral MalDH of Halobacteria (LCAHa) and compared their biochemical properties with those of five modern halobacterial MalDHs. We monitored the stability of these various MalDHs, their oligomeric states and enzymatic properties, as a function of concentration for different salts in the solvent. We found that a variety of evolutionary processes such as amino acid replacement, gene duplication, loss of MalDH gene and replacement owing to horizontal transfer resulted in significant differences in solubility, stability and catalytic properties between these enzymes in the three Halobacteriales, Haloferacales and Natrialbales orders since the LCAHa MalDH.We also showed how a stability trade-off might favor the emergence of new properties during adaptation to diverse environmental conditions. Altogether, our results suggest a new view of halophilic protein adaptation in Archaea.

Less Unfavorable Salt Bridges on the Enzyme Surface Result in More Organic Cosolvent Resistance

Biocatalysis for the synthesis of fine chemicals is highly attractive but usually requires organic (co-)solvents (OSs). However, native enzymes often have low activity and resistance in OSs and at elevated temperatures. Herein, we report a smart salt bridge design strategy for simultaneously improving OS resistance and thermostability of the model enzyme, Bacillus subtilits Lipase A (BSLA). We combined comprehensive experimental studies of 3450 BSLA variants and molecular dynamics simulations of 36 systems. Iterative recombination of four beneficial substitutions yielded superior resistant variants with up to 7.6-fold (D64K/D144K) improved resistance toward three OSs while exhibiting significant thermostability (thermal resistance up to 137-fold, and half-life up to 3.3-fold). Molecular dynamics simulations revealed that locally refined flexibility and strengthened hydration jointly govern the highly increased resistance in OSs and at 50-100 °C. The salt bridge redesign provides protein engineers with a powerful and likely general approach to design OSs- and/or thermal-resistant lipases and other α/β-hydrolases.

Hydrophobic dewetting in gating and regulation of transmembrane protein ion channels

Water is at the heart of almost all biological phenomena, without which no life that we know of would have been possible. It is a misleadingly complex liquid that exists in near coexistence with the vapor phase under ambient conditions. Confinement within a hydrophobic cavity can tip this balance enough to drive a cooperative dewetting transition. For a nanometer-scale pore, the dewetting transition leads to a stable dry state that is physically open but impermeable to ions. This phenomenon is often referred to as hydrophobic gating. Numerous transmembrane protein ion channels have now been observed to utilize hydrophobic gating in their activation and regulation. Here, we review recent theoretical, simulation, and experimental studies that together have started to establish the principles of hydrophobic gating and discuss how channels of various sizes, topologies, and biological functions can utilize these principles to control the thermodynamic properties of water within their interior pores for gating and regulation. Exciting opportunities remain in multiple areas, particularly on direct experimental detection of hydrophobic dewetting in biological channels and on understanding how the cell may control the hydrophobic gating in regulation of ion channels.

Embedding Quantum Statistical Excitations in a Classical Force Field

Quantum-mechanically driven charge polarization and charge transfer are ubiquitous in biomolecular systems, controlling reaction rates, allosteric interactions, ligand-protein binding, membrane transport, and dynamically driven structural transformations. Molecular dynamics (MD) simulations of these processes require quantum mechanical (QM) information in order to accurately describe their reactive dynamics. However, current techniques-empirical force fields, subsystem approaches, ab initio MD, and machine learning-vary in their ability to achieve a consistent chemical description across multiple atom types, and at scale. Here we present a physics-based, atomistic force field, the ensemble DFT charge-transfer embedded-atom method, in which QM forces are described at a uniform level of theory across all atoms, avoiding the need for explicit solution of the Schrödinger equation or large, precomputed training data sets. Coupling between the electronic and atomistic length scales is effected through an ensemble density functional theory formulation of the embedded-atom method originally developed for elemental materials. Charge transfer is expressed in terms of ensembles of ionic state basis densities of individual atoms, and charge polarization, in terms of atomic excited-state basis densities. This provides a highly compact yet general representation of the force field, encompassing both local and system-wide effects. Charge rearrangement is realized through the evolution of ensemble weights, adjusted at each dynamical time step via chemical potential equalization.

Nonthermal excitation effects mediated by sub-terahertz radiation on hydrogen exchange in ubiquitin

Water dynamics in the hydration layers of biomolecules play crucial roles in a wide range of biological functions. A hydrated protein contains multiple components of diffusional and vibrational dynamics of water and protein, which may be coupled at ∼0.1-THz frequency (10-ps timescale) at room temperature. However, the microscopic description of biomolecular functions based on various modes of protein-water-coupled motions remains elusive. A novel approach for perturbing the hydration dynamics in the subterahertz frequency range and probing them at the atomic level is therefore warranted. In this study, we investigated the effect of klystron-based, intense 0.1-THz excitation on the slow dynamics of ubiquitin using NMR-based measurements of hydrogen-deuterium exchange. We demonstrated that the subterahertz irradiation accelerated the hydrogen-deuterium exchange of the amides located in the interior of the protein and hydrophobic surfaces while decelerating this exchange in the amides located in the surface loop and short 3₁₀ helix regions. This subterahertz-radiation-induced effect was qualitatively contradictory to the increased-temperature-induced effect. Our results suggest that the heterogeneous water dynamics occurring at the protein-water interface include components that are nonthermally excited by the subterahertz radiation. Such subterahertz-excited components may be linked to the slow function-related dynamics of the protein.

Structure-Function Dataset Reveals Environment Effects within a Fluorescent Protein Model System*

Anisotropic environments can drastically alter the spectroscopy and photochemistry of molecules, leading to complex structure-function relationships. We examined this using fluorescent proteins as easy-to-modify model systems. Starting from a single scaffold, we have developed a range of 27 photochromic fluorescent proteins that cover a broad range of spectroscopic properties, including the determination of 43 crystal structures. Correlation and principal component analysis confirmed the complex relationship between structure and spectroscopy, but also allowed us to identify consistent trends and to relate these to the spatial organization. We find that changes in spectroscopic properties can come about through multiple underlying mechanisms, of which polarity, hydrogen bonding and presence of water molecules are key modulators. We anticipate that our findings and rich structure/spectroscopy dataset can open opportunities for the development and evaluation of new and existing protein engineering methods.

Computing inelastic neutron scattering spectra from molecular dynamics trajectories

Inelastic neutron scattering (INS) provides a weighted density of phonon modes. Currently, INS spectra can only be interpreted for perfectly crystalline materials because of high computational cost for electronic simulations. INS has the potential to provide detailed morphological information if sufficiently large volumes and appropriate structural variety are simulated. Here, we propose a method that allows direct comparison between INS data with molecular dynamics simulations, a simulation method that is frequently used to simulate semicrystalline/amorphous materials. We illustrate the technique by analyzing spectra of a well-studied conjugated polymer, poly(3-hexylthiophene-2,5-diyl) (P3HT) and conclude that our technique provides improved volume and structural variety, but that the classical force field requires improvement before the morphology can be accurately interpreted.

Morphogenic fields: A coming of age

Morphogenesis, the coming-into-being of living organisms, was first described in the 4th century BC by Aristotle, progenitor of biology and embryology. Over the centuries it has been the subject of innumerable commentaries by philosophers, theologians and scientists but no consensus has ever been reached as to its causes. In the late 19th century, along with the emergence of cellular and molecular biology, embryology underwent a renaissance and became a topic of great interest and research. Early on the discipline divided into two opposing factions, those who attempted to explain fetal development on the basis of cellular and molecular mechanisms, and those who invoked the presence of organizing fields. The morphogenic field was first articulated in the early decades of the 20th century by multiple researchers independently of each other. The field became an extremely useful conceptual tool by which to explain a wide range of developmental phenomena. While embryology and genetics originally formed a unified discipline, during the 1930s  40 s geneticists became progressively skeptical of the field notion. The discovery of the DNA structure by Watson and Crick in the early 1950s decisively settled matters and thereafter the two disciplines pursued different lines of inquiry. After World War II embryology and the field concept went into a decades-long decline. By the 1980s an increasing number of scientists began to critically reexamine the morphogenic field concept and it underwent a second renaissance. In this paper I examine the development and evolution of the field concept, both experimentally and conceptually, and highlight the failure of genetic mechanisms to explain morphogenesis. I provide three instances from the medical literature of developmental phenomena which are only explainable on the basis of morphogenic field dynamics and argue that the field concept must be readmitted into mainstream scientific discourse.

Examining the Role of Different Molecular Interactions on Activation Energies and Activation Volumes in Liquid Water

There are a large number of force fields available to model water in molecular dynamics simulations, which each have their own strengths and weaknesses in describing the behavior of the liquid. One particular weakness in many of these models is their description of dynamics away from ambient conditions, where their ability to reproduce measurements is mixed. To investigate this issue, we use the recently developed fluctuation theory for dynamics to directly evaluate measures of the local temperature and pressure dependence: the activation energy and the activation volume. We examine these activation parameters for hydrogen-bond jump exchange times, OH reorientation times, and diffusion coefficients calculated from the SPC/E, SPC/Fw, TIP3P-PME, TIP3P-PME/Fw, OPC3, TIP4P/2005, TIP4P/Ew, E3B2, and E3B3 water models. Activation energy decompositions available through the fluctuation theory approach provide mechanistic insight into the origins of different temperature dependences between the various models, as well as the influence of three-body effects and flexibility.

Direct observation of coherent femtosecond solvent reorganization coupled to intramolecular electron transfer

It is well known that the solvent plays a critical role in ultrafast electron-transfer reactions. However, solvent reorganization occurs on multiple length scales, and selectively measuring short-range solute-solvent interactions at the atomic level with femtosecond time resolution remains a challenge. Here we report femtosecond X-ray scattering and emission measurements following photoinduced charge-transfer excitation in a mixed-valence bimetallic (FeⁱⁱRuⁱⁱⁱ) complex in water, and their interpretation using non-equilibrium molecular dynamics simulations. Combined experimental and computational analysis reveals that the charge-transfer excited state has a lifetime of 62 fs and that coherent translational motions of the first solvation shell are coupled to the back electron transfer. Our molecular dynamics simulations identify that the observed coherent translational motions arise from hydrogen bonding changes between the solute and nearby water molecules upon photoexcitation, and have an amplitude of tenths of ångströms, 120-200 cm^-1 frequency and ~100 fs relaxation time. This study provides an atomistic view of coherent solvent reorganization mediating ultrafast intramolecular electron transfer.

2H NMR study on temperature-dependent water dynamics in amino-acid functionalized silica nanopores

We prepare various amino-acid functionalized silica pores with diameters of ∼6 nm and study the temperature-dependent reorientation dynamics of water in these confinements. Specifically, we link basic Lys, neutral Ala, and acidic Glu to the inner surfaces and combine ²H nuclear magnetic resonance spin-lattice relaxation and line shape analyses to disentangle the rotational motions of the surfaces groups and the crystalline and liquid water fractions coexisting below partial freezing. Unlike the crystalline phase, the liquid phase shows reorientation dynamics, which strongly depends on the chemistry of the inner surfaces. The water reorientation is slowest for the Lys functionalization, followed by Ala and Glu and, finally, the native silica pores. In total, the rotational correlation times of water at the different surfaces vary by about two orders of magnitude, where this span is largely independent of the temperature in the range ∼200-250 K.

Deciphering the evolution of supramolecular nanofibers in solution and solid-state: a combined microscopic and spectroscopic approach

Supramolecular self-assembly of small organic molecules has emerged as a powerful tool to construct well-defined micro- and nanoarchitecture through fine-tuning a range of intermolecular interactions. The size, shape, and optical properties of these nanostructures largely depend on the specific assembly of the molecular building units, temperature and polarity of the medium, and external stimuli. The engineering of supramolecular self-assembled nanostructures with morphology-dependent tunable emission is in high demand due to the promising scope in nanodevices and molecular machines. However, probing the evolution of molecular aggregates from the solution and directing the self-assembly process in a pre-defined fashion are challenging. In the present study, we have deciphered the sequential evolution of supramolecular nanofibers from solution to spherical and oblong-shaped nanoparticles through the variation of solvent polarity, tuning the hydrophobic-hydrophilic interactions. An intriguing case of molecular self-assembly has been elucidated employing a newly designed π-conjugated thiophene derivative (TPAn) through a combination of steady-state absorption, emission measurements, fluorescence correlation spectroscopy (FCS), and electron microscopy. The FCS analysis and microscopy results revealed that the small-sized nanofibers in the dispersion further agglomerated upon solvent evaporation, resulting in a network of nanofibers. Stimuli-responsive reversible interconversion between a network of nanofibers and spherical nanoaggregates was probed both in dispersion and solvent-evaporated state. The evolution of organic nanofibers and a subtle control over the self-assembly process demonstrated in the current investigation provide a general paradigm to correlate the size, shape, and emission properties of fluorescent molecular aggregates in complex heterogeneous media, including a human cell.

Protein flexibility reduces solvent-mediated friction barriers of ligand binding to a hydrophobic surface patch

Solvent fluctuations have been explored in detail for idealized and rigid hydrophobic model systems, but so far it has remained unclear how internal protein motions and their coupling to the surrounding solvent affect the dynamics of ligand binding to biomolecular surfaces. Here, molecular dynamics simulations were used to elucidate the solvent-mediated binding of a model ligand to the hydrophobic surface patch of ubiquitin. The ligand's friction profiles reveal pronounced long-time correlations and enhanced friction in the vicinity of the protein, similar to idealized hydrophobic surfaces. Interestingly, these effects are shaped by internal protein motions. Protein flexibility modulates water density fluctuations near the hydrophobic surface patch and smooths out the friction profile of ligand binding.

Coupling Monte Carlo, Variational Implicit Solvation, and Binary Level-Set for Simulations of Biomolecular Binding

We develop a hybrid approach that combines the Monte Carlo (MC) method, a variational implicit-solvent model (VISM), and a binary level-set method for the simulation of biomolecular binding in an aqueous solvent. The solvation free energy for the biomolecular complex is estimated by minimizing the VISM free-energy functional of all possible solute-solvent interfaces that are used as dielectric boundaries. This functional consists of the solute volumetric, solute-solvent interfacial, solute-solvent van der Waals interaction, and electrostatic free energy. A technique of shifting the dielectric boundary is used to accurately predict the electrostatic part of the solvation free energy. Minimizing such a functional in each MC move is made possible by our new and fast binary level-set method. This method is based on the approximation of surface area by the convolution of an indicator function with a compactly supported kernel and is implemented by simple flips of numerical grid cells locally around the solute-solvent interface. We apply our approach to the p53-MDM2 system for which the two molecules are approximated by rigid bodies. Our efficient approach captures some of the poses before the final bound state. All-atom molecular dynamics simulations with most of such poses quickly reach the final bound state. Our work is a new step toward realistic simulations of biomolecular interactions. With further improvement of coarse graining and MC sampling, and combined with other models, our hybrid approach can be used to study the free-energy landscape and kinetic pathways of ligand binding to proteins.

Hyaluronan orders water molecules in its nanoscale extended hydration shells

Hyaluronan (HA) is an anionic, highly hydrated bio-polyelectrolyte found in the extracellular environment, like the synovial fluid between joints. We explore the extended hydration shell structure of HA in water using femtosecond elastic second-harmonic scattering (fs-ESHS). HA enhances orientational water-water correlations. Angle-resolved fs-ESHS measurements and nonlinear optical modeling show that HA behaves like a flexible chain surrounded by extended shells of orientationally correlated water. We describe several ways to determine the concentration-dependent size and shape of a polyelectrolyte in water, using the amount of water oriented by the polyelectrolyte charges as a contrast agent. The spatial extent of the hydration shell is determined via temperature-dependent measurements and can reach up to 475 nm, corresponding to a length of 1600 water molecules. A strong isotope effect, stemming from nuclear quantum effects, is observed when light water (H2O) is replaced by heavy water (D2O), amounting to a factor of 4.3 in the scattered SH intensity.

Infrared spectroscopy probes ion binding geometries

Infrared (IR) spectroscopy is a well-established technique for probing the structure, behavior, and surroundings of molecules in their native environments. Its characteristics-most specifically high structural sensitivity, ready applicability to aqueous samples, and broad availability-make it a valuable enzymological technique, particularly for the interrogation of ion binding sites. While IR spectroscopy of the "garden variety" (steady state at room temperature with wild-type proteins) is versatile and powerful in its own right, the combination of IR spectroscopy with specialized experimental schemes for leveraging ultrafast time resolution, protein labeling, and other enhancements further extends this utility. This book chapter provides the fundamental physical background and literature context essential for harnessing IR spectroscopy in the general context of enzymology with specific focus on interrogation of ion binding. Studies of lanthanide ions binding to calmodulin are highlighted as illustrative examples of this process. Appropriate sample preparation, data collection, and spectral interpretation are discussed from a detail-oriented and practical perspective with the goal of facilitating the reader's rapid progression from reading words in a book to collecting and analyzing their own data in the lab.

Characterizing the Water Wire in the Gramicidin Channel Found by Monte Carlo Sampling Using Continuum Electrostatics and in Molecular Dynamics Trajectories with Conventional or Polarizable Force Fields

Water molecules play a key role in all biochemical processes. They help define the shape of proteins, and they are reactant or product in many reactions and are released as ligands are bound. They facilitate the transfer of protons through transmembrane proton channel, pump and transporter proteins. Continuum electrostatics (CE) force fields used by program Multiconformation CE (MCCE) capture electrostatic interactions in biomolecules with an implicit solvent, which captures the averaged solvent water equilibrium properties. Hybrid CE methods can use explicit water molecules within the protein surrounded by implicit solvent. These hybrid methods permit the study of explicit hydrogen bond networks within the protein and allow analysis of processes such as proton transfer reactions. Yet hybrid CE methods have not been rigorously tested. Here, we present an explicit treatment of water molecules in the Gramicidin A (gA) channel using MCCE and compare the resulting distributions of water molecules and key hydration features against those obtained with explicit solvent Molecular Dynamics (MD) simulations with the nonpolarizable CHARMM36 and polarizable Drude force fields. CHARMM36 leads to an aligned water wire in the channel characterized by a large absolute net water dipole moment; the MCCE and Drude analysis lead to a small net dipole moment as the water molecules change orientation within the channel. The correct orientation is not as yet known, so these calculations identify an open question.

Investigation of Cooperative Modes for Collective Molecules Using Grid-Based Principal Component Analysis

A grid-based principal component analysis method (GBPCA) has been developed and implemented to investigate the modes of collective molecules present in the cube in the grid system from their trajectories using molecular dynamics (MD) simulation. This method is applied to the simulations of water, methane, and hydrated proteins. In the cases of single molecules, GBPCA demonstrates that while individual molecules interact with other molecules and move randomly, the collective molecules nevertheless produce cooperative principal component (PC) modes. Notably, these PC modes of the collective molecules formed vorticial or symmetrical structures, with some resembling torus-like structures. Such structures were observed not only in simulations with models that reproduced real molecules, such as water or methane, but also in simulations with artificial molecules that were modified to interact weakly or only repulsively. On the contrary, molecules without any interactions did not exhibit the cooperative PC modes that lead to the formation of specific structures. These results imply that when molecules subject to intermolecular forces are located in a space in which they can interact with one another, the PC modes of the collective molecules will form vorticial and torus-like structures. In the cases of hydrated proteins, the results reveal that the collective molecules in the protein and water molecules showed cooperative PC modes, forming vorticial or torus-like structures, as in the cases of the water or methane molecules. The PC modes of the water molecules close to the protein were influenced by the protein motion and indicated cooperative modes with the protein. Moreover, proteins with the same folds had similar PC mode structures. Besides, the PC modes of the proteins composed of multiple domains appeared independently in the protein domain. Hence, it can be understood that the cooperative motions of the collective molecules of the protein and surrounding water molecules contribute strongly to the protein's structure and function. These results are expected to help our understanding of the dynamics of molecules in the liquid state.

Enhancing water sampling of buried binding sites using nonequilibrium candidate Monte Carlo

Water molecules can be found interacting with the surface and within cavities in proteins. However, water exchange between bulk and buried hydration sites can be slow compared to simulation timescales, thus leading to the inefficient sampling of the locations of water. This can pose problems for free energy calculations for computer-aided drug design. Here, we apply a hybrid method that combines nonequilibrium candidate Monte Carlo (NCMC) simulations and molecular dynamics (MD) to enhance sampling of water in specific areas of a system, such as the binding site of a protein. Our approach uses NCMC to gradually remove interactions between a selected water molecule and its environment, then translates the water to a new region, before turning the interactions back on. This approach of gradual removal of interactions, followed by a move and then reintroduction of interactions, allows the environment to relax in response to the proposed water translation, improving acceptance of moves and thereby accelerating water exchange and sampling. We validate this approach on several test systems including the ligand-bound MUP-1 and HSP90 proteins with buried crystallographic waters removed. We show that our BLUES (NCMC/MD) method enhances water sampling relative to normal MD when applied to these systems. Thus, this approach provides a strategy to improve water sampling in molecular simulations which may be useful in practical applications in drug discovery and biomolecular design.

Water-Mediated Interactions Determine Helix Formation of Peptides in Open Nanotubes

Water-mediated interactions (WMIs) play diverse roles in molecular biology. They are particularly relevant in geometrically confined spaces such as the interior of the chaperonin, at the interface between ligands and their binding partners, and in the ribosome tunnel. Inspired in part by the geometry of the ribosome tunnel, we consider confinement effects on the stability of peptides. We describe results from replica exchange molecular dynamics simulations of a system containing a 23-alanine or 23-serine polypeptide confined to nonpolar and polar nanotubes in the gas phase and when open to a water reservoir. We quantify the effect of water in determining the preferred conformational states of these polypeptides by calculating the difference in the solvation free energy for the helix and coil states in the open nanotube in the two phases. Our simulations reveal several possibilities. We find that nanoscopic confinement preferentially stabilizes the helical state of polypeptides with hydrophobic side chains, which is explained by the entropic stabilization mechanism proposed on the basis of polymer physics. Polypeptide chains with hydrophilic side chains can adopt helical structures within nanotubes, but helix formation is sensitive to the nature of the nanotube due to WMIs. We elaborate on the potential implications of our findings to the stability of peptides in the ribosome tunnel.

Hydration and dynamics of L-glutamate ion in aqueous solution

Aqueous solutions of sodium l-glutamate (NaGlu) in the concentration range 0 < c/M ≤ 1.90 at 25 °C were investigated by dielectric relaxation spectroscopy (DRS) and statistical mechanics (1D-RISM and 3D-RISM calculations) to study the hydration and dynamics of the l-glutamate (Glu-) anion. Although at c → 0 water molecules beyond the first hydration shell are dynamically affected, Glu- hydration is rather fragile and for c ⪆ 0.3 M apparently restricted to H2O molecules hydrogen bonding to the carboxylate groups. These hydrating dipoles are roughly parallel to the anion moment, leading to a significantly enhanced effective dipole moment of Glu-. However, l-glutamate dynamics is determined by the rotational diffusion of individual anions under hydrodynamic slip boundary conditions. Thus, the lifetime of the hydrate complexes, as well as of possibly formed [Na+Glu-]0 ionpairs and l-glutamate aggregates, cannot exceed the characteristic timescale for Glu- rotation.