Water-protein interactions from high-resolution protein crystallography

CryoEM-sampling of metastable conformations appearing in cofactor-ligand association and catalysis of glutamate dehydrogenase

Kinetic aspects of enzymatic reactions are described by equations based on the Michaelis-Menten theory for the initial stage. However, the kinetic parameters provide little information on the atomic mechanism of the reaction. In this study, we analyzed structures of glutamate dehydrogenase in the initial and steady stages of the reaction using cryoEM at near-atomic resolution. In the initial stage, four metastable conformations displayed different domain motions and cofactor/ligand association modes. The most striking finding was that the enzyme-cofactor-substrate complex, treated as a single state in the enzyme kinetic theory, comprised at least three different metastable conformations. In the steady stage, seven conformations, including derivatives from the four conformations in the initial stage, made the reaction pathway complicated. Based on the visualized conformations, we discussed stage-dependent pathways to illustrate the dynamics of the enzyme in action.

Coenzyme-binding pathway on glutamate dehydrogenase suggested from multiple-binding sites visualized by cryo-electron microscopy

The structure of hexameric glutamate dehydrogenase (GDH) in the presence of the coenzyme nicotinamide adenine dinucleotide phosphate (NADP) was visualized using cryogenic transmission electron microscopy to investigate the ligand-binding pathways to the active site of the enzyme. Each subunit of GDH comprises one hexamer-forming core domain and one nucleotide-binding domain (NAD domain), which spontaneously opens and closes the active-site cleft situated between the two domains. In the presence of NADP, the potential map of GDH hexamer, assuming D3 symmetry, was determined at a resolution of 2.4 Å, but the NAD domain was blurred due to the conformational variety. After focused classification with respect to the NAD domain, the potential maps interpreted as NADP molecules appeared at five different sites in the active-site cleft. The subunits associated with NADP molecules were close to one of the four metastable conformations in the unliganded state. Three of the five binding sites suggested a pathway of NADP molecules to approach the active-site cleft for initiating the enzymatic reaction. The other two binding modes may rarely appear in the presence of glutamate, as demonstrated by the reaction kinetics. Based on the visualized structures and the results from the enzymatic kinetics, we discussed the binding modes of NADP to GDH in the absence and presence of glutamate.

Structural analysis of wild-type and Val120Thr mutant Candida boidinii formate dehydrogenase by X-ray crystallography

Candida boidinii NAD+-dependent formate dehydrogenase (CbFDH) has gained significant attention for its potential application in the production of biofuels and various industrial chemicals from inorganic carbon dioxide. The present study reports the atomic X-ray crystal structures of wild-type CbFDH at cryogenic and ambient temperatures, as well as that of the Val120Thr mutant at cryogenic temperature, determined at the Turkish Light Source `Turkish DeLight'. The structures reveal new hydrogen bonds between Thr120 and water molecules in the active site of the mutant CbFDH, suggesting increased stability of the active site and more efficient electron transfer during the reaction. Further experimental data is needed to test these hypotheses. Collectively, these findings provide invaluable insights into future protein-engineering efforts that could potentially enhance the efficiency and effectiveness of CbFDH.

De novo prediction of explicit water molecule positions by a novel algorithm within the protein design software MUMBO

By mediating interatomic interactions, water molecules play a major role in protein-protein, protein-DNA and protein-ligand interfaces, significantly affecting affinity and specificity. This notwithstanding, explicit water molecules are usually not considered in protein design software because of high computational costs. To challenge this situation, we analyzed the binding characteristics of 60,000 waters from high resolution crystal structures and used the observed parameters to implement the prediction of water molecules in the protein design and side chain-packing software MUMBO. To reduce the complexity of the problem, we incorporated water molecules through the solvation of rotamer pairs instead of relying on solvated rotamer libraries. Our validation demonstrates the potential of our algorithm by achieving recovery rates of 67% for bridging water molecules and up to 86% for fully coordinated waters. The efficacy of our algorithm is highlighted further by the prediction of 3 different proteinligand complexes. Here, 91% of water-mediated interactions between protein and ligand are correctly predicted. These results suggest that the new algorithm could prove highly beneficial for structure-based protein design, particularly for the optimization of ligand-binding pockets or protein-protein interfaces.

Heterogeneous and Allosteric Role of Surface Hydration for Protein-Ligand Binding

Atomistic-level understanding of surface hydration mediating protein-protein interactions and ligand binding has been a challenge due to the dynamic nature of water molecules near the surface. We develop a computational method to evaluate the solvation free energy based on the density map of the first hydration shell constructed from all-atom molecular dynamics simulation and use it to examine the binding of two intrinsically disordered ligands to their target protein domain. One ligand is from the human protein, and the other is from the 1918 Spanish flu virus. We find that the viral ligand incurs a 6.9 kcal/mol lower desolvation penalty upon binding to the target, which is consistent with its stronger binding affinity. The difference arises from the spatially fragmented and nonuniform water density profiles of the first hydration shell. In particular, residues that are distal from the ligand-binding site contribute to a varying extent to the desolvation penalty, among which the "entropy hotspot" residues contribute significantly. Thus, ligand binding alters hydration on remote sites in addition to affecting the binding interface. The nonlocal effect disappears when the conformational motion of the protein is suppressed. The present results elucidate the interplay between protein conformational dynamics and surface hydration. Our approach of measuring the solvation free energy based on the water density of the first hydration shell is tolerant of the conformational fluctuation of protein, and we expect it to be applicable to investigating a broad range of biomolecular interfaces.

Prediction of hydrophilic and hydrophobic hydration structure of protein by neural network optimized using experimental data

The hydration structures of proteins, which are necessary for their folding, stability, and functions, were visualized using X-ray and neutron crystallography and transmission electron microscopy. However, complete visualization of hydration structures over the entire protein surface remains difficult. To compensate for this incompleteness, we developed a three-dimensional convolutional neural network to predict the probability distribution of hydration water molecules on the hydrophilic and hydrophobic surfaces, and in the cavities of proteins. The neural network was optimized using the distribution patterns of protein atoms around the hydration water molecules identified in the high-resolution X-ray crystal structures. We examined the feasibility of the neural network using water sites in the protein crystal structures that were not included in the datasets. The predicted distribution covered most of the experimentally identified hydration sites, with local maxima appearing in their vicinity. This computational approach will help to highlight the relevance of hydration structures to the biological functions and dynamics of proteins.

QM/MM study of N501 involved intermolecular interaction between SARS-CoV-2 receptor binding domain and antibody of human origin

Intermolecular interaction between key residue N501 of the epitope on SARS-CoV-2 RBD and screening antibody B38 was studied using the QM/MM and QM approach. The QM/MM optimized geometry shows that angle X-H---Y is 165° for O-H---O between mAb light chain S30 and RBD N501. High level MP2 calculations indicated the interaction between RBD N501 and S30 of B38 Fab light chain provide a relatively strong attractive force of - 3.32 kcal/mol, whereas the hydrogen bond between RBD Q498 and S30 was quantified as 0.10 kcal/mol. The decrease in ESP partial charge on hydrogen atom of hydroxyl group on S30 drops from 0.38 a.u. to 0.31 a.u., exhibiting the sharing of 0.07 a.u. from the lone pair electron oxygen of N501 due to hydrogen bond formation. The NBO occupancy of hydrogen atom also decreases from 25.79 % to 22.93 % in the hydroxyl H-O NBO bond of S30. However, the minor change of NBO hybridization of hydroxyl oxygen of S30 from sp3.00 to sp3.05 implies the rigidity of hydrogen bond tetrahedral geometry in the relative dynamic protein complex. The O-H---O angle is 165° which is close but not exactly linear. The structural requirement for sp3 hybridization of oxygen for hydroxyl group on S30 and dimension of protein likely prevent O-H---O from adopting linear geometry. The hydrogen bond strengths were also calculated using a variety of DFT methods, and the result of - 3.33 kcal/mol from the M06L method is the closest to that of the MP2 calculation. Results of this work may aid in the COVID-19 vaccine and drug screening.

General anesthetic binding mode via hydration with weak affinity and molecular discrimination: General anesthetic dissolution in interfacial water of the common binding site of GABAA receptor

The GABAA receptor (GABAAR) is a target channel for the loss of awareness of general anesthesia. General anesthetic (GA) spans a wide range of chemical structures, such as monatomic molecules, barbital acids, phenols, ethers, and alkanes. GA has a weak binding affinity, and the affinity has a characteristic that correlates with the solubility in olive oil rather than the molecular shape. The GA binding site of GABAAR is common to GAs and exists in the transmembrane domain of the GABAAR intersubunit. In this study, the mechanism of GA binding, which allows binding of various GAs with intersubunit selectivity, was elucidated from the hydration analysis of the binding site. Regardless of the diverse GA chemical structures, a strong correlation was observed between the binding free energy and total dehydration number of the binding process. The GA binding free energy was more involved in the binding dehydration and showed molecular recognition that allowed for the binding of various GA structures via binding site hydration. We regarded the GA substitution for the interfacial water molecule of the binding site as a dissolution into the interfacial hydration layer. The elucidation of the GA binding mechanism mediated by hydration at the GABAAR common binding site provides a rationale for the combined use of anesthetics in medical practice and its combination adjustments via drug interactions.

Diffusion Coefficient of Intracrystalline Water in Intrinsic Hen Egg-White Lysozyme Crystals Determined by Confocal Raman Spectroscopy

Protein crystals composed of protein molecules are expected as a novel porous material. They have high porosity, and the knowledge of the diffusion of intracrystalline water is important. In this study, the diffusion coefficient of intracrystalline water in intrinsic hen egg-white lysozyme (HEWL) crystals was determined by a method that combines confocal Raman spectroscopy and air convection with controlled relative humidity. Similar to common porous materials, the drying process of the protein crystals includes three periods: constant-rate drying, falling-rate drying, and equilibrium state. During the falling-rate drying period, the drying rate depends on the diffusion of intracrystalline water in the protein crystal. The gradient of the water content was measured using confocal Raman spectroscopy. The diffusion coefficient of the intrinsic HEWL crystals was determined as 3.1 × 10^-7 cm²/s with a water content of 36.3 vol %. The estimated diffusion coefficients of the intrinsic HEWL crystals without cross-linking were in close agreement with those of the cross-linked protein crystals. This study is timely as the knowledge of the intrinsic diffusion coefficient is useful not only for understanding the mechanism of hydration of proteins but also in practical applications such as porous materials, drug binding, and cryoprotectant soaks.

Knowledge-based prediction of DNA hydration using hydrated dinucleotides as building blocks

Database-derived water probability densities around structurally and sequentially distinct DNA dinucleotide fragments reproduce the known hydration motifs, which thus can be used as building blocks to predict DNA hydration.

Water plays an important role in stabilizing the structure of DNA and mediating its interactions. Here, the hydration of DNA was analyzed in terms of dinucleotide fragments from an ensemble of 2727 nonredundant DNA chains containing 41 853 dinucleotides and 316 265 associated first-shell water molecules. The dinucleotides were classified into categories based on their 16 sequences and the previously determined structural classes known as nucleotide conformers (NtCs). The construction of hydrated dinucleotide building blocks allowed dinucleotide hydration to be calculated as the probability of water density distributions. Peaks in the water densities, known as hydration sites (HSs), uncovered the interplay between base and sugar-phosphate hydration in the context of sequence and structure. To demonstrate the predictive power of hydrated DNA building blocks, they were then used to predict hydration in an independent set of crystal and NMR structures. In ten tested crystal structures, the positions of predicted HSs and experimental waters were in good agreement (more than 40% were within 0.5 Å) and correctly reproduced the known features of DNA hydration, for example the ‘spine of hydration’ in B-DNA. Therefore, it is proposed that hydrated building blocks can be used to predict DNA hydration in structures solved by NMR and cryo-EM, thus providing a guide to the interpretation of experimental data and computer models. The data for the hydrated building blocks and the predictions are available for browsing and visualization at the website https://watlas.datmos.org/watna/.

Wide-Line NMR Melting Diagrams, Their Thermodynamic Interpretation, and Secondary Structure Predictions for A30P and E46K α-Synuclein

Parkinson's disease is thought to be caused by aggregation of the intrinsically disordered protein, α-synuclein. Two amyloidogenic variants, A30P, and E46K familial mutants were investigated by wide-line ¹H NMR spectrometry as a completion of our earlier work on wild-type and A53T α-synuclein (Bokor M. et al. WT and A53T α-synuclein systems: melting diagram and its new interpretation. Int. J. Mol. Sci.2020, 21, 3997.). A monolayer of mobile water molecules hydrates A30P α-synuclein at the lowest potential barriers (temperatures), while E46K α-synuclein has here third as much mobile hydration, insufficient for functionality. According to wide-line ¹H NMR results and secondary structure predictions, E46K α-synuclein is more compact than the A30P variant and they are more compact than the wild type (WT) and A53T variants. Linear hydration vs potential barrier sections of A30P α-synuclein shows one and E46K shows two slopes. The different slopes of the latter between potential barriers E _a,1 and E _a,2 reflect a change in water-protein interactions. The 31-32% homogeneous potential barrier distribution of the protein-water bonds refers to a non-negligible amount of secondary structures in all four α-synuclein variants. The secondary structures detected by wide-line ¹H NMR are solvent-exposed α-helices, which are predicted by secondary structure models. β-sheets are only minor components of the protein structures as three- and eight-state predicted secondary structures are dominated by α-helices and coils.

Surface-tethered planar membranes containing the β-barrel assembly machinery: a platform for investigating bacterial outer membrane protein folding

The outer membrane of Gram-negative bacteria presents a robust physicochemical barrier protecting the cell from both the natural environment and acting as the first line of defense against antimicrobial materials. The proteins situated within the outer membrane are responsible for a range of biological functions including controlling influx and efflux. These outer membrane proteins (OMPs) are ultimately inserted and folded within the membrane by the β-barrel assembly machine (Bam) complex. The precise mechanism by which the Bam complex folds and inserts OMPs remains unclear. Here, we have developed a platform for investigating Bam-mediated OMP insertion. By derivatizing a gold surface with a copper-chelating self-assembled monolayer, we were able to assemble a planar system containing the complete Bam complex reconstituted within a phospholipid bilayer. Structural characterization of this interfacial protein-tethered bilayer by polarized neutron reflectometry revealed distinct regions consistent with known high-resolution models of the Bam complex. Additionally, by monitoring changes of mass associated with OMP insertion by quartz crystal microbalance with dissipation monitoring, we were able to demonstrate the functionality of this system by inserting two diverse OMPs within the membrane, pertactin, and OmpT. This platform has promising application in investigating the mechanism of Bam-mediated OMP insertion, in addition to OMP function and activity within a phospholipid bilayer environment.

Hydration and its Hydrogen Bonding State on a Protein Surface in the Crystalline State as Revealed by Molecular Dynamics Simulation

Protein hydration is crucial for the stability and molecular recognition of a protein. Water molecules form a hydration water network on a protein surface via hydrogen bonds. This study examined the hydration structure and hydrogen bonding state of a protein, staphylococcal nuclease, at various hydration levels in its crystalline state by all-atom molecular dynamics (MD) simulation. Hydrophilic residues were more hydrated than hydrophobic residues. As the water content increases, both types of residues were uniformly more hydrated. The number of hydrogen bonds per single water asymptotically approaches 4, the same as bulk water. The distances and angles of hydrogen bonds in hydration water in the protein crystal were almost the same as those in the tetrahedral structure of bulk water regardless of the hydration level. The hydrogen bond structure of hydration water observed by MD simulations of the protein crystalline state was compared to the Hydrogen and Hydration Database for Biomolecule from experimental protein crystals.

High-resolution structure of the amino acid transporter AdiC reveals insights into the role of water molecules and networks in oligomerization and substrate binding

Background

The L-arginine/agmatine transporter AdiC is part of the arginine-dependent extreme acid resistance system of the bacterium Escherichia coli and its pathogenic varieties such as strain E. coli O157:H7. At the present time, there is a lack of knowledge concerning the role of water molecules and networks for the structure and function of AdiC, and solute transporters in general.

Results

The structure of the L-arginine/agmatine transporter AdiC was determined at 1.7 Å resolution by X-ray crystallography. This high resolution allowed for the identification of numerous water molecules buried in the structure. In combination with molecular dynamics (MD) simulations, we demonstrate that water molecules play an important role for stabilizing the protein and key residues, and act as placeholders for atoms of the AdiC substrates L-arginine and agmatine. MD simulations unveiled flexibility and restrained mobility of gating residues W202 and W293, respectively. Furthermore, a water-filled cavity was identified at the dimer interface of AdiC. The two monomers formed bridging interactions through water-mediated hydrogen bonds. The accessibility and presence of water molecules in this cavity was confirmed with MD simulations. Point mutations disrupting the interfacial water network validated the importance of water molecules for dimer stabilization.

Conclusions

This work gives new insights into the role and importance of water molecules in the L-arginine/agmatine transporter AdiC for protein stabilization and substrate-binding site shaping and as placeholders of substrate atoms. Furthermore, and based on the observed flexibility and restrained mobility of gating residues, a mechanistic role of the gate flexibility in the transport cycle was proposed. Finally, we identified a water-filled cavity at the dimeric interface that contributes to the stability of the amino acid transporter oligomer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01102-4.

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.

Structural basis of the protochromic green/red photocycle of the chromatic acclimation sensor RcaE

Cyanobacteriochromes (CBCRs) are bilin-binding photosensors of the phytochrome superfamily that show remarkable spectral diversity. The green/red CBCR subfamily is important for regulating chromatic acclimation of photosynthetic antenna in cyanobacteria and is applied for optogenetic control of gene expression in synthetic biology. It is suggested that the absorption change of this subfamily is caused by the bilin C15-Z/C15-E photoisomerization and a subsequent change in the bilin protonation state. However, structural information and direct evidence of the bilin protonation state are lacking. Here, we report a high-resolution (1.63Å) crystal structure of the bilin-binding domain of the chromatic acclimation sensor RcaE in the red-absorbing photoproduct state. The bilin is buried within a "bucket" consisting of hydrophobic residues, in which the bilin configuration/conformation is C5-Z,syn/C10-Z,syn/C15-E,syn with the A- through C-rings coplanar and the D-ring tilted. Three pyrrole nitrogens of the A- through C-rings are covered in the α-face with a hydrophobic lid of Leu249 influencing the bilin pK _a, whereas they are directly hydrogen bonded in the β-face with the carboxyl group of Glu217. Glu217 is further connected to a cluster of waters forming a hole in the bucket, which are in exchange with solvent waters in molecular dynamics simulation. We propose that the "leaky bucket" structure functions as a proton exit/influx pathway upon photoconversion. NMR analysis demonstrated that the four pyrrole nitrogen atoms are indeed fully protonated in the red-absorbing state, but one of them, most likely the B-ring nitrogen, is deprotonated in the green-absorbing state. These findings deepen our understanding of the diverse spectral tuning mechanisms present in CBCRs.

Molecular Recognition and Self-Organization in Life Phenomena Studied by a Statistical Mechanics of Molecular Liquids, the RISM/3D-RISM Theory

There are two molecular processes that are essential for living bodies to maintain their life: the molecular recognition, and the self-organization or self-assembly. Binding of a substrate by an enzyme is an example of the molecular recognition, while the protein folding is a good example of the self-organization process. The two processes are further governed by the other two physicochemical processes: solvation and the structural fluctuation. In the present article, the studies concerning the two molecular processes carried out by Hirata and his coworkers, based on the statistical mechanics of molecular liquids or the RISM/3D-RISM theory, are reviewed.

Cryo-EM and MD infer water-mediated proton transport and autoinhibition mechanisms of Vo complex

Rotary vacuolar adenosine triphosphatases (V-ATPases) drive transmembrane proton transport through a V_o proton channel subcomplex. Despite recent high-resolution structures of several rotary ATPases, the dynamic mechanism of proton pumping remains elusive. Here, we determined a 2.7-Å cryo-electron microscopy (cryo-EM) structure of yeast V_o proton channel in nanodisc that reveals the location of ordered water molecules along the proton path, details of specific protein-lipid interactions, and the architecture of the membrane scaffold protein. Moreover, we uncover a state of V_o that shows the c-ring rotated by ~14°. Molecular dynamics simulations demonstrate that the two rotary states are in thermal equilibrium and depict how the protonation state of essential glutamic acid residues couples water-mediated proton transfer with c-ring rotation. Our cryo-EM models and simulations also rationalize a mechanism for inhibition of passive proton transport as observed for free V_o that is generated as a result of V-ATPase regulation by reversible disassembly in vivo.

Molecular Motions and Interactions in Aqueous Solutions of Thymosin-β4 , Stabilin CTD and Their 1 : 1 Complex, Studied by 1 H-NMR Spectroscopy

Wide-line ¹ H NMR measurements were extended and all results were interpreted in a thermodynamics-based new approach on aqueous solutions of thymosin-β₄ (Tβ₄ ), stabilin cytoplasmic domain (CTD), and their 1 : 1 complex. Energy distributions of potential barriers controlling the motion of protein-bound water molecules were determined. Heterogeneous and homogeneous regions were found in the protein-water interface. The measure of heterogeneity of this interface gives quantitative value for the portion of disordered parts in the protein. Ordered structural elements were found extending up to ∼20 % of the individual whole proteins. About 40 % of the binding sites of free Tβ₄ get involved in bonds holding the complex together. The complex has the most heterogeneous solvent accessible surface (SAS) in terms of protein-water interactions. The complex is more disordered than Tβ₄ or stabilin CTD. The greater SAS area of the complex is interpreted as a clear sign of its open structure.

Hydration of Carboxyl Groups: A Route toward Molecular Recognition?

On Earth, water plays an active role in cellular life, over several scales of distance and time. At a nanoscale, water drives macromolecular conformation through hydrophobic forces and at short times acts as a proton donor/acceptor providing charge carriers for signal transmission. At longer times and larger distances, water controls osmosis, transport, and protein mobility. Neutron diffraction experiments augmented by computer simulation, show that the three-dimensional shape of the hydration shell of carboxyl and carboxylate groups belonging to different molecules is characteristic of each molecule. Different hydration shells identify and distinguish specific sites with the same chemical structure. This experimental evidence suggests an active role of water also in controlling, modulating, and mediating chemical reactions involving carboxyl and carboxylate groups.

Accurate prediction of chemical shifts for aqueous protein structure on "Real World" data

Here we report a new machine learning algorithm for protein chemical shift prediction that outperforms existing chemical shift calculators on realistic data that is not heavily curated, nor eliminates test predictions ad hoc. Our UCBShift predictor implements two modules: a transfer prediction module that employs both sequence and structural alignment to select reference candidates for experimental chemical shift replication, and a redesigned machine learning module based on random forest regression which utilizes more, and more carefully curated, feature extracted data. When combined together, this new predictor achieves state-of-the-art accuracy for predicting chemical shifts on a randomly selected dataset without careful curation, with root-mean-square errors of 0.31 ppm for amide hydrogens, 0.19 ppm for Hα, 0.84 ppm for C', 0.81 ppm for Cα, 1.00 ppm for Cβ, and 1.81 ppm for N. When similar sequences or structurally related proteins are available, UCBShift shows superior native state selection from misfolded decoy sets compared to SPARTA+ and SHIFTX2, and even without homology we exceed current prediction accuracy of all other popular chemical shift predictors.

Applications of water molecules for analysis of macromolecule properties

Water molecules maintain proteins' structures, functions, stabilities and dynamics. They can occupy certain positions or pass quickly via a protein's interior. Regardless of their behaviour, water molecules can be used for the analysis of proteins' structural features and biochemical properties. Here, we present a list of several software programs that use the information provided by water molecules to: i) analyse protein structures and provide the optimal positions of water molecules for protein hydration, ii) identify high-occupancy water sites in order to analyse ligand binding modes, and iii) detect and describe tunnels and cavities. The analysis of water molecules' distribution and trajectories sheds a light on proteins' interactions with small molecules, on the dynamics of tunnels and cavities, on protein composition and also on the functionality, transportation network and location of functionally relevant residues. Finally, the correct placement of water molecules in protein crystal structures can significantly improve the reliability of molecular dynamics simulations.

Mobility of water and of protein atoms at the protein-water interface, monitored by anisotropic atomic displacement parameters, are largely uncorrelated

A non-redundant set of 231 protein crystal structures refined at a resolution better than (or equal to) 1 Å was extracted from the Protein Data Bank and the degree of conformational rigidity at the protein-water interface was examined by means of the Hirshfeld test and by comparing the orientations of the anisotropic Us for contacting protein and water atoms. Contacts between protein and water atoms are more rigid that contacts between water atoms and the degree of rigidity increases for shorter contacts and for more hydrogen-bonded atoms. Nevertheless, water and protein atoms are not rigidly held together. On the contrary, they seem to have little influence on their mobility to such an extent that hydration water, different from the protein atoms, cannot be considered to be properly in the solid state.

Efficient asymmetric synthesis of chiral alcohols using high 2-propanol tolerance alcohol dehydrogenase SmADH2 via an environmentally friendly TBCR system

Based on substrate-coupled cofactor regeneration system, a high 2-propanol tolerance SmADH2 together with TBCR system can synthesise structurally diverse chiral alcohols at a high substrate loading with only 1.25 equivalents of 2-propanol.

Dipole-dipole interactions between tryptophan side chains and hydration water molecules dominate the observed dynamic stokes shift of lysozyme

The fluorescence intensity of tryptophan residues in hen egg-white lysozyme was measured up to 500 ps after the excitation by irradiation pulses at 290 nm. From the time-dependent variation of fluorescence intensity in a wavelength range of 320-370 nm, the energy relaxation in the dynamic Stokes shift was reconstructed as the temporal variation in wavenumber of the estimated fluorescence maximum. The relaxation was approximated by two exponential curves with decay constants of 1.2 and 26.7 ps. To interpret the relaxation, a molecular dynamics simulation of 75 ns was conducted for lysozyme immersed in a water box. From the simulation, the energy relaxation in the electrostatic interactions of each tryptophan residue was evaluated by using a scheme derived from the linear response theory. Dipole-dipole interactions between each of the Trp62 and Trp123 residues and hydration water molecules displayed an energy relaxation similar to that experimentally observed regarding time constants and magnitudes. The side chains of these residues were partly or fully exposed to the solvent. In addition, by inspecting the variation in dipole moments of the hydration water molecules around lysozyme, it was suggested that the observed relaxation could be attributed to the orientational relaxation of hydration water molecules participating in the hydrogen-bond network formed around each of the two tryptophan residues.

Unit-cell response of tetragonal hen egg white lysozyme upon controlled relative humidity variation

Variation of relative humidity (rH) greatly affects the internal order of solvent-based protein crystals, and the rearrangement of molecules can be efficiently recorded in distinct diffraction patterns. This study focuses on this topic, reporting the effect of rH variation experiments on hen egg white lysozyme (HEWL) polycrystalline precipitates of tetragonal symmetry using X-ray powder diffraction (XRPD). In situ XRPD data were collected on HEWL specimens during dehydration and rehydration processes using laboratory instrumentation. A known polymorph [space group P43212, a = 79.07181 (1), c = 38.0776 (1) Å] was identified during gradual dehydration from 95 to 63% rH and vice versa. Pawley analysis of collected data sets and accurate extraction of unit-cell parameters indicated a characteristic evolution of the tetragonal axes with rH. In addition, there is a low humidity level below which samples do not retain their crystallinity. This work illustrates the accuracy of laboratory XRPD as a probe for time-resolved studies of proteins and in situ investigations of gradual structural modifications upon rH variation. These experiments provide essential information for improving production and post-production practices of microcrystalline protein-based pharmaceuticals.

Ice formation and solvent nanoconfinement in protein crystals

Ice formation within protein crystals is a major obstacle to the cryocrystallographic study of protein structure, and has limited studies of how the structural ensemble of a protein evolves with temperature in the biophysically interesting range from ∼260 K to the protein-solvent glass transition near 200 K. Using protein crystals with solvent cavities as large as ∼70 Å, time-resolved X-ray diffraction was used to study the response of protein and internal solvent during rapid cooling. Solvent nanoconfinement suppresses freezing temperatures and ice-nucleation rates so that ice-free, low-mosaicity diffraction data can be reliably collected down to 200 K without the use of cryoprotectants. Hexagonal ice (Ih) forms in external solvent, but internal crystal solvent forms stacking-disordered ice (Isd) with a near-random stacking of cubic and hexagonal planes. Analysis of powder diffraction from internal ice and single-crystal diffraction from the host protein structure shows that the maximum crystallizable solvent fraction decreases with decreasing crystal solvent-cavity size, and that an ∼6 Å thick layer of solvent adjacent to the protein surface cannot crystallize. These results establish protein crystals as excellent model systems for the study of nanoconfined solvent. By combining fast cooling, intense X-ray beams and fast X-ray detectors, complete structural data sets for high-value targets, including membrane proteins and large complexes, may be collected at ∼220-240 K that have much lower mosaicities and comparable B factors, and that may allow more confident identification of ligand binding than in current cryocrystallographic practice.

Functional Hydration Behavior: Interrelation between Hydration and Molecular Properties at Lipid Membrane Interfaces

Water is an abundant commodity and has various important functions. It stabilizes the structure of biological macromolecules, controls biochemical activities, and regulates interfacial/intermolecular interactions. Common aspects of interfacial water can be obtained by overviewing fundamental functions and properties at different temporal and spatial scales. It is important to understand the hydrogen bonding and structural properties of water and to evaluate the individual molecular species having different hydration properties. Water molecules form hydrogen bonds with biomolecules and contribute to the adjustment of their properties, such as surface charge, hydrophilicity, and structural flexibility. In this review, the fundamental properties of water molecules and the methods used for the analyses of water dynamics are summarized. In particular, the interrelation between the hydration properties, determined by molecules, and the properties of molecules, determined by their hydration properties, are discussed using the lipid membrane as an example. Accordingly, interesting water functions are introduced that provide beneficial information in the fields of biochemistry, medicine, and food chemistry.

The Role of Hydrogen Bonding in the Folding/Unfolding Process of Hydrated Lysozyme: A Review of Recent NMR and FTIR Results

The biological activity of proteins depends on their three-dimensional structure, known as the native state. The main force driving the correct folding mechanism is the hydrophobic effect and when this folding kinetics is altered, aggregation phenomena intervene causing the occurrence of illnesses such as Alzheimer and Parkinson’s diseases. The other important effect is performed by water molecules and by their ability to form a complex network of hydrogen bonds whose dynamics influence the mobility of protein amino acids. In this work, we review the recent results obtained by means of spectroscopic techniques, such as Fourier Transform Infrared (FTIR) and Nuclear Magnetic Resonance (NMR) spectroscopies, on hydrated lysozyme. In particular, we explore the Energy Landscape from the thermal region of configurational stability up to that of the irreversible denaturation. The importance of the coupling between the solute and the solvent will be highlighted as well as the different behaviors of hydrophilic and hydrophobic moieties of protein amino acid residues.

Understanding the function of water during the gelation of globular proteins by temperature-dependent near infrared spectroscopy

Water plays an indispensable role in the gelation of proteins, but its function still remains unclear. In this work, the variation of water species with the structural changes of globular proteins was investigated using temperature-dependent near infrared (NIR) spectroscopy. Ovalbumin (OVA) was used as a model protein, which forms a gel-like structure as the temperature increases through three phases, i.e., phase I (native), phase II (molten globule state), and phase III (gel state). The structural change and the content variation of different water species in the three phases of gelation were analyzed by two-dimensional correlation NIR spectroscopy and Gaussian fitting. A decrease in the water species with two hydrogen bonds (S2) was found and the change follows the same phases as OVA. In the first two phases, the change occurs after those of other water species but in the third phase, the change is faster than that of free water species. The result indicates that in the native and molten globule states, S2 is located in the hydration shell of OVA to maintain the stability of the protein structure, and then in the gel state, high temperature weakens the hydrogen bonding of S2 and leads to the destruction of the hydration shell, making OVA clusters form a gel structure.

Toward comprehensive measurement of protein hydration dynamics: Facilitation of NMR-based methods by reverse micelle encapsulation

Protein-water interactions are a fundamental determinant of protein structure and function. Despite their importance, the molecular details of water orientations and dynamics near protein surfaces remain poorly understood, largely due to the difficulty of measuring local water mobility near the protein in a site-resolved fashion. Solution NMR-based measurement of water mobility via the nuclear Overhauser effect was presented as a method for performing comprehensive, site-resolved measurements of water dynamics many years ago. Though this approach yielded extensive insight on the dynamics and locations of waters buried within proteins, its promise for measuring surface hydration dynamics was impeded by various technical barriers. Over the past several years, however, this approach has been pursued anew with the aid of reverse micelle encapsulation of proteins of interest. The confined environment of the reverse micelle resolves many of these barriers and permits site-resolved measurement of relative water dynamics across much of the protein surface. Here, the development of this strategy for measuring hydration dynamics is reviewed with particular focus on the important remaining challenges to its widespread application.

Effects of protein-crystal hydration and temperature on side-chain conformational heterogeneity in monoclinic lysozyme crystals

The modulation of main-chain and side-chain conformational heterogeneity and solvent structure in monoclinic lysozyme crystals by dehydration (related to water activity) and temperature is examined. Decreasing the relative humidity (from 99 to 11%) and decreasing the temperature both lead to contraction of the unit cell, to an increased area of crystal contacts and to remodeling of primarily contact and solvent-exposed residues. Both lead to the depopulation of some minor side-chain conformers and to the generation of new conformations. Side-chain modifications and main-chain r.m.s.d.s associated with cooling from 298 to 100 K depend on relative humidity and are minimized at 85% relative humidity (r.h.). Dehydration from 99 to 93% r.h. and cooling from 298 to 100 K result in a comparable number of remodeled residues, with dehydration-induced remodeling somewhat more likely to arise from contact interactions. When scaled to equivalent temperatures based on unit-cell contraction, the evolution of side-chain order parameters with dehydration shows generally similar features to those observed on cooling to T = 100 K. These results illuminate the qualitative and quantitative similarities between structural perturbations induced by modest dehydration, which routinely occurs in samples prepared for 298 and 100 K data collection, and cryocooling. Differences between these perturbations in terms of energy landscapes and occupancies, and implications for variable-temperature crystallography between 180 and 298 K, are discussed. It is also noted that remodeling of a key lysozyme active-site residue by dehydration, which is associated with a radical decrease in the enzymatic activity of lysozyme powder, arises due to a steric clash with the residue of a symmetry mate.

Molecular Motions and Interactions in Aqueous Solutions of Thymosin-β4 , Stabilin C-Terminal Domain (CTD) and Their 1:1 Complex Studied by 1 H NMR Spectroscopy

Wide-line ¹ H NMR measurements were extended and all results were reinterpreted in a new thermodynamics-based approach to study aqueous solutions of thymosin-β₄ (Tβ₄ ), stabilin C-terminal domain (CTD) and their 1:1 complex. The energy distributions of the potential barriers, which control motion of protein-bound water molecules, were determined. Heterogeneous and homogeneous regions were found at the protein-water interface. The measure of heterogeneity gives a quantitative value for the portion of disordered parts in the protein. Ordered structural elements were found extending up to 20 % of the whole proteins. About 40 % of the binding sites of free Tβ₄ become involved in bonds holding the complex together. The complex has the most heterogeneous solvent accessible surface (SAS) in terms of protein-water interactions. The complex is more disordered than Tβ₄ or stabilin CTD. The greater SAS area of the complex is interpreted as a clear sign of its open structure.

Slow Dynamics of Tryptophan-Water Networks in Proteins

Water has a profound effect on the dynamics of biomolecules and governs many biological processes, leading to the concept that function is slaved to solvent dynamics within and surrounding the biomolecule. Protein conformational changes on μs-ms time scales are frequently associated with protein function, but little is known about the behavior of protein-bound water on these time scales. Here we have used NMR relaxation dispersion measurements to probe the tryptophan indoles in the enzyme dihydrofolate reductase (DHFR). We find that during structural changes on the μs-ms time scale, large chemical shift changes are often observed for the NH proton on the indole ring, while relatively smaller chemical shift changes are observed for the ring nitrogen atom. Comparison with experimental chemical shifts and density functional theory-based chemical shift predictions show that during the structural change the tryptophan indole NHs remain bound to water, but the geometry of the protein-bound water networks changes. These results establish that relaxation dispersion measurements can indirectly probe water dynamics and indicate that water can influence, or be influenced by, protein conformational changes on the μs-ms time scale. Our data show that structurally conserved bound water molecules can play a critical role in transmitting information between functionally important regions of the protein and provide evidence that internal protein motions can be coupled through the mediation of hydrogen-bonded water bound in the protein structure.

Biosynthesis of Neokestose Laurate Catalyzed by Candida antarctica Lipase B and Its Antimicrobial Activity against Food Pathogenic and Spoilage Bacteria

To increase the functionality and broaden the potential application of neokestose, neokestose laurate was biosynthesized using Candida antarctica lipase B as biocatalyst, for which a mixture of 20% DMSO in 2-methyl-2-butanol (v/v) was chosen as the reaction medium. The optimum conditions for biosynthesis were as follows: a molar ratio of vinyl laurate to neokestose of 12, a temperature of 50 °C, molecular sieves of 100 g/L, and enzyme loading of 10 g/L. Under the optimal conditions, the conversion rate was achieved over 80%. The synthesized chemical 6'-O-lauroylneokestose confirmed by nuclear magnetic resonance (NMR) exhibited good emulsification with critical micelle concentration (CMC) of 352 μM and broad antibacterial activity against Gram-positive bacteria such as Staphylococcus aureus, Listeria monocytogenes, Streptococcus mutans, Bacillus subtilis, and Bacillus cereus. Conclusively, 6'-O-lauroylneokestose was evidenced to be a dual-functional agent with emulsification and antibacterial activity, showing promising application potential in the food industry.

Influences of lone-pair electrons on directionality of hydrogen bonds formed by hydrophilic amino acid side chains in molecular dynamics simulation

The influence of lone-pair electrons on the directionality of hydrogen bonds that are formed by oxygen and nitrogen atoms in the side chains of nine hydrophilic was investigated using molecular dynamics simulations. The simulations were conducted using two types of force fields; one incorporated lone-pair electrons placed at off-atom sites and the other did not. The density distributions of the hydration water molecules around the oxygen and nitrogen atoms were calculated from the simulation trajectories, and were compared with the empirical hydration distribution functions, which were constructed from a large number of hydration water molecules found in the crystal structures of proteins. Only simulations using the force field explicitly incorporating lone-pair electrons reproduced the directionality of hydrogen bonds that is observed in the empirical distribution functions for the deprotonated oxygen and nitrogen atoms in the sp 2-hybridization. The amino acids that include such atoms are functionally important glutamate, aspartate, and histidine. Therefore, a set of force field that incorporates lone-pair electrons as off-atom charge sites would be effective for considering hydrogen bond formation by these amino acids in molecular dynamics simulation studies.

Water Distribution, Dynamics, and Interactions with Alzheimer's β-Amyloid Fibrils Investigated by Solid-State NMR

Water is essential for protein folding and assembly of amyloid fibrils. Internal water cavities have been proposed for several amyloid fibrils, but no direct structural and dynamical data have been reported on the water dynamics and site-specific interactions of water with the fibrils. Here we use solid-state NMR spectroscopy to investigate the water interactions of several Aβ40 fibrils. 1H spectral lineshapes, T2 relaxation times, and two-dimensional (2D) 1H-13C correlation spectra show that there are five distinct water pools: three are peptide-bound water, while two are highly dynamic water that can be assigned to interfibrillar water and bulk-like matrix water. All these water pools are associated with the fibrils on the nanometer scale. Water-transferred 2D correlation spectra allow us to map out residue-specific hydration and give evidence for the presence of a water pore in the center of the three-fold symmetric wild-type Aβ40 fibril. In comparison, the loop residues and the intramolecular strand-strand interface have low hydration, excluding the presence of significant water cavities in these regions. The Osaka Aβ40 mutant shows lower hydration and more immobilized water than wild-type Aβ40, indicating the influence of peptide structure on the dynamics and distribution of hydration water. Finally, the highly mobile interfibrillar and matrix water exchange with each other on the time scale of seconds, suggesting that fibril bundling separates these two water pools, and water molecules must diffuse along the fibril axis before exchanging between these two environments. These results provide insights and experimental constraints on the spatial distribution and dynamics of water pools in these amyloid fibrils.

How many water molecules are detected in X-ray protein crystal structures?

The positions of several water molecules can be determined in protein crystallography, either buried in internal cavities or at the protein surface. It is important to be able to estimate the expected number of these water molecules to facilitate crystal structure determination. Here, a multiple Poisson regression model implemented on nearly 10 000 protein crystal structures shows that the number of detectable water molecules depends on eight variables: crystallographic resolution, R factor, percentage of solvent in the crystal, average B factor of the protein atoms, percentage of amino acid residues in loops, average solvent-accessible surface area of the amino acid residues, grand average of hydropathy of the protein(s) in the asymmetric unit and normalized number of heteroatoms that are not water molecules. Furthermore, a secondary analysis tested the effect of different software packages. Given the values of these eight variables, it is possible to compute the expected number of water molecules detectable in electron-density maps with reasonable accuracy (as suggested by an external validation of the model).

The "Sticky Patch" Model of Crystallization and Modification of Proteins for Enhanced Crystallizability

Crystallization of macromolecules has long been perceived as a stochastic process, which cannot be predicted or controlled. This is consistent with another popular notion that the interactions of molecules within the crystal, i.e., crystal contacts, are essentially random and devoid of specific physicochemical features. In contrast, functionally relevant surfaces, such as oligomerization interfaces and specific protein-protein interaction sites, are under evolutionary pressures so their amino acid composition, structure, and topology are distinct. However, current theoretical and experimental studies are significantly changing our understanding of the nature of crystallization. The increasingly popular "sticky patch" model, derived from soft matter physics, describes crystallization as a process driven by interactions between select, specific surface patches, with properties thermodynamically favorable for cohesive interactions. Independent support for this model comes from various sources including structural studies and bioinformatics. Proteins that are recalcitrant to crystallization can be modified for enhanced crystallizability through chemical or mutational modification of their surface to effectively engineer "sticky patches" which would drive crystallization. Here, we discuss the current state of knowledge of the relationship between the microscopic properties of the target macromolecule and its crystallizability, focusing on the "sticky patch" model. We discuss state-of-the-art in silico methods that evaluate the propensity of a given target protein to form crystals based on these relationships, with the objective to design variants with modified molecular surface properties and enhanced crystallization propensity. We illustrate this discussion with specific cases where these approaches allowed to generate crystals suitable for structural analysis.

Unifying the microscopic picture of His-containing turns: from gas phase model peptides to crystallized proteins

The presence in crystallized proteins of a local anchoring between the side chain of a His residue, located in the central position of a γ- or β-turn, and its local main chain environment, is assessed by the comparison of protein structures with relevant isolated model peptides.

Waters in room temperature and cryo protein crystal structures

Since it has been observed that low temperature protein crystal structures may differ from room temperature structures, it is necessary to compare systematically the protein hydration structure in low and room protein crystal structures. High quality data sets of protein structures were built in an extremely rigorous manner and crystal symmetry was included in the identification of four types of water molecules (buried in the protein core, deeply inserted into crevices at the protein surface, first and second hydration layers). More water molecules are observed at low temperature only if the resolution is better than 2.1–2.3 Å. At worse resolution, temperature does not play any role. The numerous water molecules that become detectable at low temperature and at higher resolution are more mobile, relative to the protein average flexibility. Despite that, the occupancy does not depend on temperature. It can be hypothesized that water structure and around proteins and hydrogen bond network do not depend on the temperature, at least in the temperature range examined here. At low temperature more water molecules are detected because the average flexibility of all the atoms decreases, so that also water molecules that are considerably more mobile than the average atoms become observable in the electron density maps.

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.

SAXS/SANS on Supercharged Proteins Reveals Residue-Specific Modifications of the Hydration Shell

Water molecules in the immediate vicinity of biomacromolecules, including proteins, constitute a hydration layer characterized by physicochemical properties different from those of bulk water and play a vital role in the activity and stability of these structures, as well as in intermolecular interactions. Previous studies using solution scattering, crystallography, and molecular dynamics simulations have provided valuable information about the properties of these hydration shells, including modifications in density and ionic concentration. Small-angle scattering of x-rays (SAXS) and neutrons (SANS) are particularly useful and complementary techniques to study biomacromolecular hydration shells due to their sensitivity to electronic and nuclear scattering-length density fluctuations, respectively. Although several sophisticated SAXS/SANS programs have been developed recently, the impact of physicochemical surface properties on the hydration layer remains controversial, and systematic experimental data from individual biomacromolecular systems are scarce. Here, we address the impact of physicochemical surface properties on the hydration shell by a systematic SAXS/SANS study using three mutants of a single protein, green fluorescent protein (GFP), with highly variable net charge (+36, -6, and -29). The combined analysis of our data shows that the hydration shell is locally denser in the vicinity of acidic surface residues, whereas basic and hydrophilic/hydrophobic residues only mildly modify its density. Moreover, the data demonstrate that the density modifications result from the combined effect of residue-specific recruitment of ions from the bulk in combination with water structural rearrangements in their vicinity. Finally, we find that the specific surface-charge distributions of the different GFP mutants modulate the conformational space of flexible parts of the protein.

Changes in hydration structure are necessary for collective motions of a multi-domain protein

Conformational motions of proteins are necessary for their functions. To date, experimental studies measuring conformational fluctuations of a whole protein structure have revealed that water molecules hydrating proteins are necessary to induce protein functional motions. However, the underlying microscopic mechanism behind such regulation remains unsolved. To clarify the mechanism, multi-domain proteins are good targets because it is obvious that water molecules between domains play an important role in domain motions. Here, we show how changes in hydration structure microscopically correlate with large-amplitude motions of a multi-domain protein, through molecular dynamics simulation supported by structural analyses and biochemical experiments. We first identified collective domain motions of the protein, which open/close an active-site cleft between domains. The analyses on changes in hydration structure revealed that changes in local hydration in the depth of the cleft are necessary for the domain motion and vice versa. In particular, ‘wetting’/‘drying’ at a hydrophobic pocket and ‘adsorption’/‘dissociation’ of a few water molecules at a hydrophilic crevice in the cleft were induced by dynamic rearrangements of hydrogen-bond networks, and worked as a switch for the domain motions. Our results microscopically demonstrated the importance of hydrogen-bond networks of water molecules in understanding energy landscapes of protein motions.