Analysis of solvent structure in proteins using neutron D2O-H2O solvent maps: pattern of primary and secondary hydration of trypsin.

A method of determining the water structure in protein crystals is described using neutron solvent difference maps. These maps are obtained by comparing the changes in diffracted intensities between two data sets, one in which H2O is the major solvent constituent, and a second in which D2O is the solvent medium. To a good first approximation, the protein atom contributions to the scattering intensities in both data sets are equal and cancel, but since H2O and D2O have very different neutron-scattering properties, their differences are accentuated to reveal an accurate representation of the solvent structure. The method also employs a series of density modification steps that impose known physical constraints on the density distribution function in the unit cell by making real space modifications directly to the density maps. Important attributes of the method are that (1) it is less subjective in the assignment of water positions than X-ray analysis; (2) there is threefold improvement in the signal-to-noise ratio for the solvent density; and (3) the iterative density modification produces a low-biased representation of the solvent density. Tests showed that water molecules with as low as 10% occupancy could be confidently assigned. About 300 water sites were assigned for trypsin from the refined solvent density; 140 of these sites were defined in the maps as discrete peaks, while the remaining were found within less-ordered channels of density. There is a very good correspondence between the sites in the primary hydration layer and waters found in the X-ray structure. Most water sites are clustered into H-bonding networks, many of which are found along intermolecular contact zones. The bound water is equally distributed between contacting apolar and polar atoms at the protein interface. A common occurrence at hydrophobic surfaces is that apolar atoms are circumvented by one or more waters that are part of a larger water network. When the effects on surface accessibility by neighboring molecules in the crystal lattice are taken into consideration, only about 29% of the surface does not interface ordered water. About 25% of the ordered water is found in the second hydration sphere. In many instances these waters bridge larger clusters of primary layer waters. It is apparent that, in certain regions of the crystal, the organization of ordered water reflects the characteristics of the crystal environment more than those of trypsin's surface alone.     

Monte carlo simulation of the solvent structure in crystals of a hydrated cyclic peptide

The Monte Carlo simulation of the structure of the 16 ordered and disordered waters in the unit cell of crystals of the cyclic peptide cyclo(-L -Ala-L -Pro-D -Phe)₂ is reported. The water structure has been characterized in terms of the statistically averaged positions of the individual molecules, their root-mean-square movements about these positions, the probability of finding a water in a given spatial position in the crystal (probability maps), and examination of individual configurations of the system. In this way a picture is obtained of the water structure, including water orientations (hydrogen positions), the hydrogen-bonding network, and fluctuations in these structural features, to a degree hitherto unavailable either from experimental or theoretical studies. In addition, the variation in water structure in various peptide environments was studied and correlated with the energetics of the individual water molecules. Variations in the crystalline environment of different water molecules lead to energy differences of as much as 4–5 kcal/mol in their average energies. Similarly, differences are observed in the water–peptide and water–water components of the energy. Two different water potentials were tested. The results were compared with experimental data in terms of mean positions, root-mean-square movements, and the Fourier transform of the simulated water structure. The agreement factor (R factor) calculated from the theoretical water probability distribution was 18.8% compared to the x-ray value of 14.5%, and the value of 28% when the water is omitted.     

Water structure in vitamin B12 coenzyme crystals. I. Analysis of the neutron and x-ray solvent densities.

The disordered solvent distribution in crystals of vitamin B12 coenzyme was examined using the methods of high-resolution neutron and x-ray diffraction. One set of neutron (0.95 A) and two sets of x-ray (0.94 and 1.1 A) data were collected and the resulting models were extensively refined using least-squares and Fourier syntheses. The solvent regions were analyzed in two stages: first, main sites were assigned to the well defined regions of solvent density and refined using least squares; second, continuous sites were assigned representing the more disordered diffuse and elongated regions of solvent density. During the analysis an acetone molecule was also located. Water networks were formulated from the assigned sites in the above models and also from those assigned in the original structure determination (Lenhert, 1968), using criteria that included hydrogen bonding (derived from small crystal hydrates), van der Waals contact distances, side-chain disorder, water molecule orientations, and the presence or absence of foreign solvent. The well established networks extend throughout all the solvent regions of the crystal with interesting orientational arrangements of the individual waters around both polar and apolar groups of the coenzyme molecule. The networks were seen to be consistent among each of the four models in terms of occupying relatively similar positions. However, the occupancy values of the individual networks varied between the models; some networks were clearly visible in one but attenuated in another. The specific details of the water structure (bonding geometries, short-range nonbonded contacts, orientations of the waters, polar and apolar interactions, etc.) are described in the following paper.     

Continuum electrostatic analysis of preferred solvation sites around proteins in solution

To understand water-protein interactions in solution, the electrostatic field is calculated by solving the Poisson-Boltzmann equation, and the free energy surface of water is mapped by translating and rotating an explicit water molecule around the protein. The calculation is applied to T4 lysozyme with data available on the conservation of solvent binding sites in 18 crystallographically independent molecules. The free energy maps around the ordered water sites provide information on the relationship between water positions in crystal structure and in solution. Results show that almost all conserved sites and the majority of nonconserved sites are within 1.3 A of local free energy minima. This finding is in sharp contrast to the behavior of randomly placed water molecules in the boundary layer, which, on the average, must travel more than 3 A to the nearest free energy minimum. Thus, the solvation sites are at least partially determined by protein-water interactions rather than by crystal packing alone. The characteristic water residence times, obtained from the free energies at the local minima, are in good agreement with nuclear magnetic resonance experiments. Only about half of the potential sites show up as ordered water in the 1.7 A resolution X-ray structure. Crystal packing interactions can stabilize weak or mobile potential sites (in fact, some ordered water positions are not close to free energy minima) or can prevent water from occupying certain sites. Apart from a few buried water molecules that are strong binders, the free energies are not very different for conserved and nonconserved sites. We show that conservation of a water site between two crystals occurs if the positions of protein atoms, primarily contributing to the free energy at the local minimum, do not substantially change from one structure to the other. This requirement can be correlated with the nature of the side chain contacting the water molecule in the site.     

Conserved water molecules in a large family of microbial ribonucleases.

We systematically analyzed the crystallographically determined water molecules of all known structures of RNase T1 and compared them to the ordered solvent in a large number of related microbial nucleases. To assess the crystallographers' impact on the interpretation of the solvent structure, we independently refined five validation structures from diffraction data derived from five isomorphous crystals of RNase T1. We also compared the positions of water molecules found in 11 published isomorphous RNase T1 inhibitor complexes. These data suggest that the positions of most of the waters located on the surface of a protein and that are well-determined in the experimental electron density maps are determined primarily by crystal packing forces. Water molecules with less well-defined electron density are in general unique to one or a small number of crystal structures. Only a small number of the well-defined waters are found to be independent of the crystal environment. These waters have a low accessible surface area and B-factor, and tend to be conserved in the crystal structures of a number of evolutionary related ribonucleases as well. A single water molecule is found conserved in all known microbial ribonucleases.     

The structure of 2Zn pig insulin crystals at 1.5 A resolution

The paper describes the arrangement of the atoms within rhombohedral crystals of 2Zn pig insulin as seen in electron density maps calculated from X-ray data extending to 1.5 A (1 A = 10(-10) m = 10(-1) nm) at room temperature and refined to R = 0.153. The unit cell contains 2 zinc ions, 6 insulin molecules and about 3 x 283 water molecules. The atoms in the protein molecules appear well defined, 7 of the 102 side chains in the asymmetric unit have been assigned alternative disordered positions. The electron density over the water molecules has been interpreted in terms of 349 sites, 217 weighted 1.0, 126 weighted 0.5, 5 at 0.33 and 1 at 0.25 giving ca. 282 molecules. The positions and contacts of all the residues belonging to the two A and B chains of the asymmetric unit are shown first and then details of their arrangement in the two insulin molecules, 1 and 2, which are different. The formation from these molecules of a compact dimer and the further aggregation of three dimers to form a hexamer around two zinc ions, follows. It appears that in the packing of the hexamers in the crystal there are conflicting influences; too-close contacts between histidine B5 residues in neighbouring hexamers are probably responsible for movements of atoms at the beginning of the A chain of one of the two molecules of the dimer that initiate movements in other parts, particularly near the end of the B chain. At every stage of the building of the protein structure, residues to chains of definite conformation, molecules, dimers, hexamers and crystals, we can trace the effect of the packing of like groups to like, aliphatic groups together, aromatic groups together, hydrogen-bonded structures, positive and negative ions. Between the protein molecules, the water is distributed in cavities and channels that are continuous throughout the crystals. More than half the water molecules appear directly hydrogen bonded to protein atoms. These are generally in contact with other water molecules in chains and rings of increasing disorder, corresponding with their movement through the crystals. Within the established crystal structure we survey next the distribution of hydrogen bonds within the protein molecules and between water and protein and water and water; all but eight of the active atoms in the protein form at least one hydrogen bond.(ABSTRACT TRUNCATED AT 400 WORDS)     

Solvent interactions in nucleic acid crystal hydrates

A detailed knowledge of structural and energetic aspects of water-nucleic acid interactions is essential for understanding the role of solvent in stabilizing the various helical forms of nucleic acids. In this study, computer simulation techniques have been used to predict structural properties of solvent networks in small nucleic acid crystal hydrates. A detailed comparison of predicted and experimental results on the structure of the solvent networks is presented and includes an analysis of both the local environment and hydrogen bond pattern of each water molecule. A correlation between the environment of each unique water molecule and its energetic properties (such a dipole moment and binding energy) is seen. As in the previous studies on small amino acid hydrate crystals, non-pair additive (cooperative) effects are found to be non-negligible. It is concluded that the potential functions used in this initial study lead to simulated solvent networks in reasonable agreement with experimental data. Thus, it is now feasible to use them in studies of hydration of larger helical fragments of nucleic acids of more direct biological interest.     

Monte Carlo studies on water in the dCpG/proflavin crystal hydrate

The extensive water network identified in the crystallographic studies of the dCpG/Proflavin hydrate by Neidle, Berman and Shieh (Nature 288, 129, 1980) forms an ideal test case for a) assessing the accuracy of theoretical calculations on nucleic acid--water systems based on statistical thermodynamic computer simulation, and b) the possible use of computer simulation in predicting the water positions in crystal hydrates for use in the further refinement and interpretation of diffraction data. Monte Carlo studies have been carried out on water molecules in the unit cell of dCpG/proflavin, with the nucleic acid complex fixed and the condensed phase environment of the system treated by means of periodic boundary conditions. Intermolecular interactions are described by potential functions representative of quantum mechanical calculations developed by Clementi and coworkers, and widely used in recent studies of the aqueous hydration of various forms of DNA fragments. The results are analyzed in terms of hydrogen bond topology, hydrogen bond distances and energies, mean water positions, and water crystal probability density maps. Detailed comparison of calculated and experimentally observed results are given, and the sensitivity of results to choice of potential is determined by comparison with simulation results based on a set of empirical potentials.     

Monte Carlo Studies on Water in the dCpG/Proflavin Crystal Hydrate

Abstract

The extensive water network identified in the crystallographic studies of the dCpG/Proflavin hydrate by Neidle, Berman and Shieh (Nature 288, 129, 1980) forms an ideal test case for a) assessing the accuracy of theoretical calculations on nucleic acid—water systems based on statistical thermodynamic computer simulation, and b) the possible use of computer simulation in predicting the water positions in crystal hydrates for use in the further refinement and interpretation of diffraction data. Monte Carlo studies have been carried out on water molecules in the unit cell of dCpG/proflavin, with the nucleic acid complex fixed and the condensed phase environment of the system treated by means of periodic boundary conditions. Intermolecular interactions are described by potential functions representative of quantum mechanical calculations developed by Clementi and coworkers, and widely used in recent studies of the aqueous hydration of various forms of DNA fragments. The results are analyzed in terms of hydrogen bond topology, hydrogen bond distances and energies, mean water positions, and water crystal probability density maps. Detailed comparison of calculated and experimentally observed results are given, and the sensitivity of results to choice of potential is determined by comparison with simulation results based on a set of empirical potentials.     

Display and interpretation of solvent electron density distributions in insulin crystals

In macromolecular crystallography, three-dimensional contour surfaces are useful for interactive computer graphics displays of the protein electron density but are less effective for presenting static images of large volumes of solvent density. A raster-based computer graphics program which displays depth-cued projections of continuous density distributions has been developed to analyze the distribution of solvent atoms in macromolecular crystals. Maps of the water distribution in the cubic insulin crystal show some well-ordered waters, which are bound to surrounding protein atoms by multiple hydrogen bonds, and an ill-defined solvent structure at a greater distance from the protein surface. Molecular dynamics calculations were used to assist in the interpretation of the time-varying solvent structure within two enclosed cavities in the crystal. Two water molecules that ligate a sodium ion were almost immobile during the simulation but the majority of water molecules were found to move rapidly between the density maxima identified from the crystallographic refinement.     

Crystal structures of the signal transducing protein GlnK from Thermus thermophilus HB8

The Thermus thermophilus HB8 genome encodes a signal transducing PII protein, GlnK. The crystal structures of GlnK have been determined in two different space groups, P2(1)2(1)2(1) and P3(1)21. The PII protein has the T-loop, which is essential for interactions with receptor proteins. In both crystal forms, three GlnK molecules form a trimer in the asymmetric unit. In one P2(1)2(1)2(1) crystal form, the three T-loops in the trimer are disordered, while in another P2(1)2(1)2(1) crystal form, the T-loop from one molecule in the trimer is ordered. In the P3(1)21 crystal, one T-loop is ordered while the other two T-loops are disordered. The conformations of the ordered T-loops significantly differ between the two crystal forms; one makes the alpha-helix in the middle of the T-loop, while the other has an extension of the beta-hairpin. Two different conformations are captured by the crystal contacts. The observation of multiple T-loop conformations suggests that the T-loop could potentially exhibit "polysterism," which would be important for interactions with receptor proteins. The crystal structures of the nucleotide-bound forms, GlnK.ATP and GlnK.ADP, have also been determined. ATP/ADP binding within a cleft at the interface of two adjacent T. thermophilus GlnK monomers might affect the conformation of the T-loop.     

The contribution of electrostatics to hydrogen exchange in the unfolded protein state

Although electrostatics have long been recognized to play an important role in hydrogen exchange (HX) with solvent, the quantitative assessment of its magnitude in the unfolded state has hitherto been lacking. This limits the utility of HX as a quantitative method to study protein stability, folding, and dynamics. Using the intrinsically disordered human protein α-synuclein as a proxy for the unfolded state, we show that a hybrid mean-field approach can effectively compute the electrostatic potential at all backbone amide positions along the chain. From the electrochemical potential, a fourfold reduction in hydroxide concentration near the protein backbone is predicted for the C-terminal domain, a prognosis that is in direct agreement with experimentally derived protection factors from NMR spectroscopy. Thus, impeded HX for the C-terminal region of α-synuclein is not the result of intramolecular hydrogen bonding and/or structure formation.     

Ordered water structure in an A-DNA octamer at 1.7 A resolution

The crystal structure of the deoxyoctamer d(G-G-Br U-A-BrU-A-C-C) was refined to a resolution of 1.7 A using combined diffractometer and synchrotron data. The analysis was carried out independently in two laboratories using different procedures. Although the final results are identical the comparison of the two approaches highlights potential problems in the refinement of oligonucleotides when only limited data are available. As part of the analysis the positions of 84 solvent molecules in the asymmetric unit were established. The DNA molecule is highly solvated, particularly the phosphate-sugar back-bone and the functional groups of the bases. The major groove contains, in the central BrU-A-BrU-A region, a ribbon of water molecules forming closed pentagons with shared edges. These water molecules are linked to the base O and N atoms and to the solvent chains connecting the O-1 phosphate oxygen atoms on each strand. The minor groove is also extensively hydrated with a continuous network in the central region and other networks at each end. The pattern of hydration is briefly compared with that observed in the structure of a B-dodecamer.     

Atomic crystal and molecular dynamics simulation structures of human carbonic anhydrase II: insights into the proton transfer mechanism

Human carbonic anhydrase II (HCA II) is a zinc-metalloenzyme that catalyzes the reversible interconversion of CO2 and HCO3-. The rate-limiting step of this catalysis is the transfer of a proton between the Zn-bound solvent molecule and residue His64. In order to fully characterize the active site structural features implicated in the proton transfer mechanism, the refined X-ray crystal structure of uncomplexed wild type HCA II to 1.05 A resolution with an Rcryst value of 12.0% and an Rfree value of 15.1% has been elucidated. This structure provides strong clues as to the pathway of the intramolecular proton transfer between the Zn-bound solvent and His64. The structure emphasizes the role of the solvent network, the unique positioning of solvent molecule W2, and the significance of the dual conformation of His64 in the active site. The structure is compared with molecular dynamics (MD) simulation calculations of the Zn-bound hydroxyl/His64+ (charged) and the Zn-bound water/His64 (uncharged) HCA II states. A comparison of the crystallographic anisotropic atomic thermal parameters and MD simulation root-mean-square fluctuation values show excellent agreement in the atomic motion observed between the two methods. It is also interesting that the observed active site solvent positions in the crystal structure are also the most probable positions of the solvent during the MD simulations. On the basis of the comparative study of the MD simulation results, the HCA II crystal structure observed is most likely in the Zn-bound water/His64 state. This conclusion is based on the following observations: His64 is mainly (80%) orientated in an inward conformation; electron density omit maps infer that His64 is not charged in an either inward or outward conformation; and the Zn-bound solvent is most likely a water molecule.     

Refined atomic model of the four-layer aggregate of the tobacco mosaic virus coat protein at 2.4-A resolution.

Previous x-ray studies (2.8-A resolution) on crystals of tobacco mosaic virus coat protein grown from solutions containing high salt have characterized the structure of the protein aggregate as a dimer of a bilayered cylindrical disk formed by 34 chemically identical subunits. We have determined the crystal structure of the disk aggregate at 2.4-A resolution using x-ray diffraction from crystals maintained at cryogenic temperatures. Two regions of interest have been extensively refined. First, residues of the low-radius loop region, which were not modeled previously, have been traced completely in our electron density maps. Similar to the structure observed in the virus, the right radial helix in each protomer ends around residue 87, after which the protein chain forms an extended chain that extends to the left radial helix. The left radial helix appears as a long alpha-helix with high temperature factors for the main-chain atoms in the inner portion. The side-chain atoms in this region (residues 90-110) are not visible in the electron density maps and are assumed to be disordered. Second, interactions between subunits in the symmetry-related central A pair have been determined. No direct protein-protein interactions are observed in the major overlap region between these subunits; all interactions are mediated by two layers of ordered solvent molecules. The current structure emphasizes the importance of water in biological macromolecular assemblies.     

Alternation of DNA and solvent layers in the A form of d(GGCGCC) obtained by ethanol crystallization

We have determined by X-ray crystallography the structure of the hexamer duplex d(GGCGCC)2 in the A-form using ethanol as a precipitant. The same sequence had previously been crystallized in the B-form, but with 2-methyl-2,4-pentanediol as a precipitant. It appears that ethanol precipitation is a useful method to induce the formation of A-form crystals of DNA. Packing of the molecules in the crystal has unique features: the known interaction of A-DNA duplexes between terminal base-pairs and the minor groove of neighbor molecules is combined with a superstructure consisting in an alternation of DNA layers and solvent layers (water/ions). This organization in layers has been observed before, also with hexamers in the A conformation which crystallize in the same space group (C2221). The solvent layer has a precise thickness, although very few ordered water molecules can be detected. Another feature of this crystal is its large unit cell, which gives rise to an asymmetric unit with three hexamer duplexes. One of the three duplexes is quite different from the other two in several aspects: the number of base pairs per turn, the twist pattern, the mean value of the twist angle and the fact that one terminal base-pair is not stacked as part of the duplex and appears to be disordered. So the variability in conformation of this sequence is remarkable.     

Networks of water molecules in a proflavine deoxydinucleoside phosphate complex

Using previously reported ab-initio atom-atom potentials for the interactions of a water molecule with phosphates, sugars and bases and newly computed ab-initio atom-atom potentials for the interaction between a proflavine ion and water, we have analyzed with the Monte-Carlo Metropolis method networks of water molecules hydrating a 2:2 complex of proflavine and deoxycytidylyl-3',5'-guanosine, recently studied with X-ray crystallography. From our simulations we have i) verified the quality of our atom-atom potentials by obtaining patterns of oxygen atoms in very good agreement with the X-ray patterns for the minor groove and in reasonable agreement in the major groove, ii) predicted the water's hydrogen atoms positions and iii) preliminarily predicted the number of water molecules not reported in the X-ray study but present in the major groove. The above data, even if preliminary, and the analyses on the energetics of the water-water, water-proflavine and water-dCpG interactions indicate that very detailed accounts on the water filaments in the above crystal can be obtained optimally by merging computer and X-ray experiments.     

Protonation status of metal-bound ligands can be determined by quantum refinement

The protonation status of key residues and bound ligands are often important for the function of a protein. Unfortunately, protons are not discerned in normal protein crystal structures, so their positions have to be determined by more indirect methods. We show that the recently developed quantum refinement method can be used to determine the position of protons in crystal structures. By replacing the molecular-mechanics potential, normally used in crystallographic refinement, by more accurate quantum chemical calculations, we get information about the ideal structure of a certain protonation state. By comparing the refined structures of different protonation states, the one that fits the crystallographic raw data best can be decided using four criteria: the R factors, electron density maps, strain energy, and divergence from the unrestrained quantum chemical structure. We test this method on alcohol dehydrogenase, for which the pK(a) of the zinc-bound solvent molecule is experimentally known. We show that we can predict the correct protonation state for both a deprotonated alcohol and a neutral water molecule.     

Life: Self-Directing with Unlimited Variability on Self-Speeding

Abstract

The crystal structure of the deoxyoctamer d(G-G-^Br U-A-^BrU-A-C-C) was refined to a resolution of 1.7Å using combined diffractometer and synchrotron data. The analysis was carried out independently in two laboratories using different procedures. Although the final results are identical the comparison of the two approaches highlights potential problems in the refinement of oligonucleotides when only limited data are available.

As part of the analysis the positions of 84 solvent molecules in the asymmetric unit were established. The DNA molecule is highly solvated, particularly the phosphate-sugar backbone and the functional groups of the bases. The major groove contains, in the central ^BrU-A-^BrU-A region, a ribbon of water molecules forming closed pentagons with shared edges. These water molecules are linked to the base O and N atoms and to the solvent chains connecting the O-1 phosphate oxygen atoms on each strand. The minor groove is also extensively hydrated with a continuous network in the central region and other networks at each end. The pattern of hydration is briefly compared with that observed in the crystal structure of a B-dodecamer.     

Networks of Water Molecules in a Proflavine Deoxydinucleoside Phosphate Complex

Abstract

Using previously reported ab-initio atom-atom potentials for the interactions of a water molecule with phosphates, sugars and bases and newly computed ab-initio atom-atom potentials for the interaction between a proflavine ion and water, we have analyzed with the Monte-Carlo Metropolis method networks of water molecules hydrating a 2:2 complex of proflavine and deoxycytidylyl-3′,5′-guanosine, recently studied with X-ray crystallography. From our simulations we have i) verified the quality of our atom-atom potentials by obtaining patterns of oxygen atoms in very good agreement with the X-ray patterns for the minor groove and in reasonable agreement in the major groove, ii) predicted the water's hydrogen atoms positions and iii) preliminarily predicted the number of water molecules not reported in the X-ray study but present in the major groove. The above data, even if preliminary, and the analyses on the energetics of the water-water, water-proflavine and water-dCpG interactions indicate that very detailed accounts on the water filaments in the above crystal can be obtained optimally by merging computer and X-ray experiments.