Hydration in proteins observed by high-resolution neutron crystallography

It is well known that water molecules surrounding a protein play important roles in maintaining its structural stability. Water molecules are known to participate in several physiological processes through the formation of hydrogen bonds. However, the hydration structures of most proteins are not known well at an atomic level at present because X-ray protein crystallography has difficulties to localize hydrogen atoms. In contrast, neutron crystallography has no problem in determining the position of hydrogens with high accuracy.1 In this article, the hydration structures of three proteins are described- myoglobin, wild-type rubredoxin, and a mutant rubredoxin-the structures of which were solved at 1.5- or 1.6-A resolution by neutron structure determination. These hydration patterns show fascinating features and the water molecules adopt a variety of shapes in the neutron Fourier maps, revealing details of intermolecular hydrogen bond formation and dynamics of hydration. Our results further show that there are strong relationships between these shapes and the water environments.     

Geometric characteristics of hydrogen bonds involving sulfur atoms in proteins

Sulfur atoms have been known to participate in hydrogen bonds (H‐bonds) and these sulfur‐containing H‐bonds (SCHBs) are suggested to play important roles in certain biological processes. This study aims to comprehensively characterize all the SCHBs in 500 high‐resolution protein structures (≤1.8 Å). We categorized SCHBs into six types according to donor/acceptor behaviors and used explicit hydrogen approach to distinguish SCHBs from those of nonhydrogen bonding interactions. It is revealed that sulfur atom is a very poor H‐bond acceptor, but a moderately good H‐bond donor. In α‐helix, considerable SCHBs were found between the sulphydryl group of cysteine residue i and the carbonyl oxygen of residue i‐4, and these SCHBs exert effects in stabilizing helices. Although for other SCHBs, they possess no specific secondary structural preference, their geometric characteristics in proteins and in free small compounds are significantly distinct, indicating the protein SCHBs are geometrically distorted. Interestingly, sulfur atom in the disulfide bond tends to form bifurcated H‐bond whereas in cysteine‐cysteine pairs prefer to form dual H‐bond. These special H‐bonds remarkably boost the interaction between H‐bond donor and acceptor. By oxidation/reduction manner, the mutual transformation between the dual H‐bonds and disulfide bonds for cysteine‐cysteine pairs can accurately adjust the structural stability and biological function of proteins in different environments. Furthermore, few loose H‐bonds were observed to form between the sulphydryl groups and aromatic rings, and in these cases the donor H is almost over against the rim rather than the center of the aromatic ring. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Hydrogen exchange properties of proteins in native and denatured states monitored by mass spectrometry and NMR.

The extent of deuterium labeling of hen lysozyme, its three-disulfide derivative, and the homologous alpha-lactalbumins, has been measured by both mass spectrometry and NMR. Different conformational states of the proteins were produced by varying the solution conditions. Alternate protein conformers were found to contain different numbers of 2H atoms. Furthermore, measurement in the gas phase of the mass spectrometer or directly in solution by NMR gave consistent results. The unique ability of mass spectrometry to distinguish distributions of 2H atoms in protein molecules is exemplified using samples prepared to contain different populations of 2H-labeled protein. A comparison of the peak widths of bovine alpha-lactalbumin in alternate solution conformations but containing the same average number of 2H atoms showed dramatic differences due to different 2H distributions in the two protein conformers. Measurement of 2H distributions by ESI-MS enabled characterization of conformational averaging and structural heterogeneity. In addition, a time course for hydrogen exchange was examined and the variation in distributions of 2H atom compared with simulations for different hydrogen exchange models. The results clearly show that exchange from the native state of bovine alpha-lactalbumin at 15 degrees C is dominated by local unfolding events.     

A systematic method for analysing the protein hydration structure of T4 lysozyme

A systematic method for the analysis of the hydration structure of proteins is demonstrated on the case study of lysozyme. The method utilises multiple structural data of the same protein deposited in the protein data bank. Clusters of high water occupancy are localised and characterised in terms of their interaction with protein. It is shown that they constitute a network of interconnected hydrogen bonds anchored to the protein molecule. The high occupancy of the clusters does not directly correlate with water–protein interaction energy as was originally hypothesised. The highly occupied clusters rather correspond to the nodes of the hydration network that have the maximum number of hydrogen bonds including both the protein atoms and the surrounding water clusters. Copyright © 2013 John Wiley & Sons, Ltd.     

Empirical solvation models can be used to differentiate native from near-native conformations of bovine pancreatic trypsin inhibitor.

Several hydration models for peptides and proteins based on solvent accessible surface area have been proposed previously. We have evaluated some of these models as well as four new ones in the context of near-native conformations of a protein. In addition, we propose an empirical site-site distance-dependent correction that can be used in conjunction with any of these models. The set of near-native structures consisted of 39 conformations of bovine pancreatic trypsin inhibitor (BPTI) each of which was a local minimum of an empirical energy function (ECEPP) in the absence of solvent. Root-mean-square (rms) deviations from the crystallographically determined structure were in the following ranges: 1.06-1.94 A for all heavy atoms, 0.77-1.36 A for all backbone heavy atoms, 0.68-1.33 A for all alpha-carbon atoms, and 1.41-2.72 A for all side-chain heavy atoms. We have found that there is considerable variation among the solvent models when evaluated in terms of concordance between the solvation free energy and the rms deviations from the crystallographically determined conformation. The solvation model for which the best concordance (0.939) with the rms deviations of the C alpha atoms was found was derived from NMR coupling constants of peptides in water combined with an exponential site-site distance dependence of the potential of mean force. Our results indicate that solvation free energy parameters derived from nonpeptide free energies of hydration may not be transferrable to peptides. Parameters derived from peptide and protein data may be more applicable to conformational analysis of proteins. A general approach to derive parameters for free energy of hydration from ensemble-averaged properties of peptides in solution is described.     

A comparison of successful and failed protein interface designs highlights the challenges of designing buried hydrogen bonds

The accurate design of new protein–protein interactions is a longstanding goal of computational protein design. However, most computationally designed interfaces fail to form experimentally. This investigation compares five previously described successful de novo interface designs with 158 failures. Both sets of proteins were designed with the molecular modeling program Rosetta. Designs were considered a success if a high‐resolution crystal structure of the complex closely matched the design model and the equilibrium dissociation constant for binding was less than 10 μM. The successes and failures represent a wide variety of interface types and design goals including heterodimers, homodimers, peptide‐protein interactions, one‐sided designs (i.e., where only one of the proteins was mutated) and two‐sided designs. The most striking feature of the successful designs is that they have fewer polar atoms at their interfaces than many of the failed designs. Designs that attempted to create extensive sets of interface‐spanning hydrogen bonds resulted in no detectable binding. In contrast, polar atoms make up more than 40% of the interface area of many natural dimers, and native interfaces often contain extensive hydrogen bonding networks. These results suggest that Rosetta may not be accurately balancing hydrogen bonding and electrostatic energies against desolvation penalties and that design processes may not include sufficient sampling to identify side chains in preordered conformations that can fully satisfy the hydrogen bonding potential of the interface.     

Protein-water interactions in ribonuclease A and angiogenin: a molecular dynamics study

It is known that water molecules play an important role in the biological functioning of proteins. The members of the ribonuclease A (RNase A) family of proteins, which are sequentially and structurally similar, are known to carry out the obligatory function of cleaving RNA and individually perform other diverse biological functions. Our focus is on elucidating whether the sequence and structural similarity lead to common hydration patterns, what the common hydration sites are and what the differences are. Extensive molecular dynamics simulations followed by a detailed analysis of protein-water interactions have been carried out on two members of the ribonuclease A superfamily-RNase A and angiogenin. The water residence times are analyzed and their relationship with the characteristic properties of the protein polar atoms, such as their accessible surface area and mean hydration, is studied. The capacity of the polar atoms to form hydrogen bonds with water molecules and participate in protein-water networks are investigated. The locations of such networks are identified for both proteins.     

Large-scale networks of hydration water molecules around proteins investigated by cryogenic X-ray crystallography.

The structures at protein-water interface, i.e. the hydration structure of proteins, have been investigated by cryogenic X-ray crystal structure analyses. Hydration structures appeared far clearer at cryogenic temperature than at ambient temperature, presumably because the motions of hydration water molecules were quenched by cooling. Based on the structural models obtained, the hydration structures were systematically analyzed with respect to the amount of water molecules, the interaction modes between water molecules and proteins, the local and the global distribution of them on the surface of proteins. The standard tetrahedral interaction geometry of water in bulk retained at the interface and enabled the three-dimensional chain connection of hydrogen bonds between hydration water molecules and polar protein atoms. Large-scale networks of hydrogen bonds covering the entire surface of proteins were quite flexible to accommodate to the large-scale conformational changes of proteins and seemed to have great influences on the dynamics and function of proteins. The present observation may provide a new concept for discussing the dynamics of proteins in aqueous solution.     

Water-protein interactions from high-resolution protein crystallography

To understand the role of water in life at molecular and atomic levels, structures and interactions at the protein-water interface have been investigated by cryogenic X-ray crystallography. The method enabled a much clearer visualization of definite hydration sites on the protein surface than at ambient temperature. Using the structural models of proteins, including several hydration water molecules, the characteristics in hydration structures were systematically analysed for the amount, the interaction geometries between water molecules and proteins, and the local and global distribution of water molecules on the surface of proteins. The tetrahedral hydrogen-bond geometry of water molecules in bulk solvent was retained at the interface and enabled the extension of a three-dimensional chain connection of a hydrogen-bond network among hydration water molecules and polar protein atoms over the entire surface of proteins. Networks of hydrogen bonds were quite flexible to accommodate and/or to regulate the conformational changes of proteins such as domain motions. The present experimental results may have profound implications in the understanding of the physico-chemical principles governing the dynamics of proteins in an aqueous environment and a discussion of why water is essential to life at a molecular level.     

Satisfaction of hydrogen-bonding potential influences the conservation of polar sidechains

Although polar amino acids tend to be found on the surface of proteins due to their hydrophilic nature, their important roles within the core of proteins are now becoming better recognized. It has long been understood that a significant number of mainchain functions will not achieve hydrogen bond satisfaction through the formation of secondary structures; in these circumstances, it is generally buried polar residues that provide hydrogen bond satisfaction. Here, we describe an analysis of the hydrogen-bonding of polar amino acids in a set of structurally aligned protein families. This allows us not only to calculate the conservation of each polar residue but also to assess whether conservation is correlated with the hydrogen-bonding potential of polar sidechains. We show that those polar sidechains whose hydrogen-bonding potential is satisfied tend to be more conserved than their unsatisfied or nonhydrogen-bonded counterparts, particularly when buried. Interestingly, these buried and satisfied polar residues are significantly more conserved than buried hydrophobic residues. Forming hydrogen bonds to mainchain amide atoms also influences conservation, with those satisfied buried polar residues that form two hydrogen bonds to mainchain amides being significantly more conserved than those that form only one or none. These results indicate that buried polar residues whose hydrogen-bonding potential is satisfied are likely to have important roles in maintaining protein structure.     

Delineation of a new structural motif involving NHN γ-turn

Macromolecules are characterized by distinctive arrangement of hydrogen bonds. Different patterns of hydrogen bonds give rise to distinct and stable structural motifs. An analysis of 4114 non-redundant protein chains reveals the existence of a three-residue, (i − 1) to (i + 1), structural motif, having two hydrogen-bonded five-membered pseudo rings (the first, an N H···OC involving the first residue, and the second being N H∙∙∙N involving the last two residues), separated by a peptide bond. There could be an additional hydrogen bond between the side-chain at (i-1) and the main-chain NH of (i + 1). The average backbone torsion angles of −76(±21)° and – 12(±17)° at i creates a tight turn in the polypeptide chain, akin to a γ-turn. Indeed, a search of three-residue fragments with restriction on the terminal C^α···C^α distance and the existence of the two pseudo rings on either side revealed the presence 14 846 cases of a variant, termed NHN γ-turn, distinct from the NHO γ-turn (2032 cases) that has traditionally been characterized by the presence of NHO hydrogen bond linking the terminal main-chain atoms. As in the latter, the newly identified γ-turns are also of two types—classical and inverse, occurring in the ratio of 1:6. The propensities of residues to occur in these turns and their secondary structural features have been enumerated. An understanding of these turns would be useful for structure prediction and loop modeling, and may serve as models to represent some of the unfolded state or disordered region in proteins.     

On the satisfaction of backbone‐carbonyl lone pairs of electrons in protein structures

Protein structures are stabilized by a variety of noncovalent interactions (NCIs), including the hydrophobic effect, hydrogen bonds, electrostatic forces and van der Waals’ interactions. Our knowledge of the contributions of NCIs, and the interplay between them remains incomplete. This has implications for computational modeling of NCIs, and our ability to understand and predict protein structure, stability, and function. One consideration is the satisfaction of the full potential for NCIs made by backbone atoms. Most commonly, backbone‐carbonyl oxygen atoms located within α‐helices and β‐sheets are depicted as making a single hydrogen bond. However, there are two lone pairs of electrons to be satisfied for each of these atoms. To explore this, we used operational geometric definitions to generate an inventory of NCIs for backbone‐carbonyl oxygen atoms from a set of high‐resolution protein structures and associated molecular‐dynamics simulations in water. We included more‐recently appreciated, but weaker NCIs in our analysis, such as n→π* interactions, Cα‐H bonds and methyl‐H bonds. The data demonstrate balanced, dynamic systems for all proteins, with most backbone‐carbonyl oxygen atoms being satisfied by two NCIs most of the time. Combinations of NCIs made may correlate with secondary structure type, though in subtly different ways from traditional models of α‐ and β‐structure. In addition, we find examples of under‐ and over‐satisfied carbonyl‐oxygen atoms, and we identify both sequence‐dependent and sequence‐independent secondary‐structural motifs in which these reside. Our analysis provides a more‐detailed understanding of these contributors to protein structure and stability, which will be of use in protein modeling, engineering and design.     

Neutron crystallography for the elucidation of enzyme catalysis

Hydrogen atoms and hydration water molecules in proteins are indispensable for many biochemical processes, especially enzymatic catalysis. The locations of hydrogen atoms in proteins are usually predicted based on X-ray structures, but it is still very difficult to know the ionization states of the catalytic residues, the hydration structure of the protein, and the characteristics of hydrogen-bonding interactions. Neutron crystallography allows the direct observation of hydrogen atoms that play crucial roles in molecular recognition and the catalytic reactions of enzymes. In this review, we present the current status of neutron crystallography in structural biology and recent neutron structural analyses of three enzymes: ascorbate peroxidase, the main protease of severe acute respiratory syndrome coronavirus 2, and copper-containing nitrite reductase.     

Self-association behaviour of atactic polymethacrylic acid in aqueous solution investigated by atomistic molecular dynamics simulations

The self-association behaviour of atactic poly(methacrylic acid) (a-PMA) in water was investigated by atomistic molecular dynamics (MD) simulations. Simulations show that interchain association of a-PMA occurs only in its un-neutralised form, by hydrogen bonding between –COOH groups, which is in agreement with the experimental observation. Chain conformations, dihedral angle distributions, hydration behaviour, scattering structure factor and enthalpy-of-hydration (i.e. aqueous solvation) were analysed as a function of concentration for un-neutralised PMA, across dilute to concentrated regimes. The average 〈R_g〉 of the chain remains unaffected in solution and also for amorphous undissolved a-PMA phase, confirming the occurrence of the approximate theta-solution condition for the first time, as revealed by simulations, in a polar hydrogen-bonding polymer aqueous solution. Chain hydration behaviour and scattering structure factor show significant changes in concentrated regime. Scattering intensity collapse occurs in concentrated PMA solution, due to the existence of the swollen regime captured for the first time by explicit-MD-simulations. The hydration of PMA is driven by H-bonding, specifically between H atoms of the COOH groups and O atoms of water molecules in the closest coordination shell. The enthalpy of hydration of PMA is dominated by PMA–water interactions (charges and H-bonding). The thermodynamic contributions of PMA–PMA and PMA–water interactions towards the electrostatics as well as the dispersion components of the total solvation-enthalpy become more favourable than water–water interactions.     

Wetting of nonconserved residue‐backbones: A feature indicative of aggregation associated regions of proteins

Aggregation is an irreversible form of protein complexation and often toxic to cells. The process entails partial or major unfolding that is largely driven by hydration. We model the role of hydration in aggregation using “Dehydrons.” “Dehydrons” are unsatisfied backbone hydrogen bonds in proteins that seek shielding from water molecules by associating with ligands or proteins. We find that the residues at aggregation interfaces have hydrated backbones, and in contrast to other forms of protein–protein interactions, are under less evolutionary pressure to be conserved. Combining evolutionary conservation of residues and extent of backbone hydration allows us to distinguish regions on proteins associated with aggregation (non‐conserved dehydron‐residues) from other interaction interfaces (conserved dehydron‐residues). This novel feature can complement the existing strategies used to investigate protein aggregation/complexation. Proteins 2016; 84:254–266. © 2015 Wiley Periodicals, Inc.     

A fast and accurate computational approach to protein ionization

We report a very fast and accurate physics-based method to calculate pH-dependent electrostatic effects in protein molecules and to predict the pK values of individual sites of titration. In addition, a CHARMm-based algorithm is included to construct and refine the spatial coordinates of all hydrogen atoms at a given pH. The present method combines electrostatic energy calculations based on the Generalized Born approximation with an iterative mobile clustering approach to calculate the equilibria of proton binding to multiple titration sites in protein molecules. The use of the GBIM (Generalized Born with Implicit Membrane) CHARMm module makes it possible to model not only water-soluble proteins but membrane proteins as well. The method includes a novel algorithm for preliminary refinement of hydrogen coordinates. Another difference from existing approaches is that, instead of monopeptides, a set of relaxed pentapeptide structures are used as model compounds. Tests on a set of 24 proteins demonstrate the high accuracy of the method. On average, the RMSD between predicted and experimental pK values is close to 0.5 pK units on this data set, and the accuracy is achieved at very low computational cost. The pH-dependent assignment of hydrogen atoms also shows very good agreement with protonation states and hydrogen-bond network observed in neutron-diffraction structures. The method is implemented as a computational protocol in Accelrys Discovery Studio and provides a fast and easy way to study the effect of pH on many important mechanisms such as enzyme catalysis, ligand binding, protein-protein interactions, and protein stability.     

A Raman scattering study of the interactions of DNA with its water of hydration

Raman spectroscopy is used to probe the nature of the hydrogen bonds which hold the water of hydration to DNA. The ～ 3450?cm^?1 molecular O–H stretching mode shows that the first six water molecules per base pair of the primary hydration shell are very strongly bound to the DNA. The observed shift in the peak position of this mode permits a determination of the length of the hydrogen bonds for these water molecules. These hydrogen bonds appear to be about 0.3?Å shorter than the hydrogen bonds in bulk water. The linewidth of this mode shows no significant changes above water contents of about 15 water molecules per base pair. This technique of using a vibrational spectroscopy to obtain structural information about the hydration shells of DNA could be used to study the hydration shells of other biomolecules.     

Human recombinant [C22A] FK506-binding protein amide hydrogen exchange rates from mass spectrometry match and extend those from NMR.

Hydrogen/deuterium exchange behavior of human recombinant [C22A] FK506 binding protein (C22A FKBP) has been determined by protein fragmentation, combined with electrospray Fourier transform ion cyclotron resonance mass spectrometry (MS). After a specified period of H/D exchange in solution, C22A FKBP was digested by pepsin under slow exchange conditions (pH 2.4, 0 degree C), and then subjected to on-line HPLC/MS for deuterium analysis of each proteolytic peptide. The hydrogen exchange rate of each individual amide hydrogen was then determined independently by heteronuclear two-dimensional NMR on 15N-enriched C22A FKBP. A maximum entropy method (MEM) algorithm makes it possible to derive the distributions of hydrogen exchange rate constants from the MS-determined deuterium exchange-in curves in either the holoprotein or its proteolytic segments. The MEM-derived rate constant distributions of C22A FKBP and different segments of C22A FKBP are compared to the rate constants determined by NMR for individual amide protons. The rate constant distributions determined by both methods are consistent and complementary, thereby validating protein fragmentation/mass spectrometry as a reliable measure of hydrogen exchange in proteins.     

Kinetic, energetic, and stereochemical factors determining the molecular recognition of proteins and nucleic acids

In the context of stereochemical modeling, it has been shown that damage to the hydration shell of proteins and nucleic acids should be confronted by considerable kinetic barriers caused by the breakage of hydrogen bonds of the shell. Since the structure of the hydration shell is determined by the surface of proteins and nucleic acids, the kinetic barriers arising during the breakage of the shell differ greatly in different regions of the biopolymer surface. In turn, this means that the probability of the participation of different surface regions of proteins and nucleic acids in intermolecular interactions should vary within a wide range; i.e., hydration shells should enhance the selectivity of molecular recognition.