Dissection of the structural and functional role of a conserved hydration site in RNase T1

A secreted ribonuclease effector from Verticillium dahliae localizes in the plant nucleus to modulate host immunity

The arms race between fungal pathogens and plant hosts involves recognition of fungal effectors to induce host immunity. Although various fungal effectors have been identified, the effector functions of ribonucleases are largely unknown. Herein, we identified a ribonuclease secreted by Verticillium dahliae (VdRTX1) that translocates into the plant nucleus to modulate immunity. The activity of VdRTX1 causes hypersensitive response (HR)-related cell death in Nicotiana benthamiana and cotton. VdRTX1 possesses a signal peptide but is unlikely to be an apoplastic effector because its nuclear localization in the plant is necessary for cell death induction. Knockout of VdRTX1 significantly enhanced V. dahliae virulence on tobacco while V. dahliae employs the known suppressor VdCBM1 to escape the immunity induced by VdRTX1. VdRTX1 homologs are widely distributed in fungi but transient expression of 24 homologs from other fungi did not yield cell death induction, suggesting that this function is specific to the VdRTX1 in V. dahliae. Expression of site-directed mutants of VdRTX1 in N. benthamiana leaves revealed conserved ligand-binding sites that are important for VdRTX1 function in inducing cell death. Thus, VdRTX1 functions as a unique HR-inducing effector in V. dahliae that contributes to the activation of plant immunity.

Decoupling between the translation and rotation of water in the proximity of a protein molecule

The interaction between water and biomacromolecules is of fundamental interest in biophysics, biochemistry and physical chemistry. By combining neutron scattering and molecular dynamics simulations on a perdeuterated protein at a series of hydration levels, we demonstrated that the translational motion of water is slowed down more significantly than its rotation, when water molecules approach the protein molecule. Further analysis of the simulation trajectories reveals that the observed decoupling results from the fact that the translational motion of water is more correlated over space and more retarded by the charged/polar residues and spatial confinement on the protein surface, than the rotation. Moreover, around the stable protein residues (with smaller atomic fluctuations), water exhibits more decoupled dynamics, indicating a connection between the observed translation-rotation decoupling in hydration water and the local stability of the protein molecule.

The fungal ribonuclease-like effector protein CSEP0064/BEC1054 represses plant immunity and interferes with degradation of host ribosomal RNA

The biotrophic fungal pathogen Blumeria graminis causes the powdery mildew disease of cereals and grasses. We present the first crystal structure of a B. graminis effector of pathogenicity (CSEP0064/BEC1054), demonstrating it has a ribonuclease (RNase)-like fold. This effector is part of a group of RNase-like proteins (termed RALPHs) which comprise the largest set of secreted effector candidates within the B. graminis genomes. Their exceptional abundance suggests they play crucial functions during pathogenesis. We show that transgenic expression of RALPH CSEP0064/BEC1054 increases susceptibility to infection in both monocotyledonous and dicotyledonous plants. CSEP0064/BEC1054 interacts in planta with the pathogenesis-related protein PR10. The effector protein associates with total RNA and weakly with DNA. Methyl jasmonate (MeJA) levels modulate susceptibility to aniline-induced host RNA fragmentation. In planta expression of CSEP0064/BEC1054 reduces the formation of this RNA fragment. We propose CSEP0064/BEC1054 is a pseudoenzyme that binds to host ribosomes, thereby inhibiting the action of plant ribosome-inactivating proteins (RIPs) that would otherwise lead to host cell death, an unviable interaction and demise of the fungus.

Powdery mildews are common plant diseases which affect important crop plants including cereals such as wheat and barley. The fungi that cause this disease are obligate biotrophs: they have an absolute requirement for living host cells which they penetrate with feeding structures called haustoria. These fungi must be highly effective at avoiding immune recognition which would lead to death of the host cell and the pathogen. We assume they do this by delivering effector proteins to the host. While several hundred secreted effectors have been described in cereal powdery mildews, it is unknown how they work. Here, we use X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy to determine the structure and interactions of the effector CSEP0064/BEC1054, representative of the largest class of effectors resembling fungal RNases. We find that this effector binds nucleic acids. Expression of the effector in plants increases susceptibility to infection. Moreover, transgenic CSEP0064/BEC1054 expression in wheat inhibits the degradation of host ribosomal RNA induced by ribosome-inactivating proteins (RIPs). We propose a novel mechanism of action for the RNase-like effectors in powdery mildews: they may act as pseudoenzymes to inhibit the host RIPs, known components of plant immune responses that lead to host cell death.

How proteins modify water dynamics

Much of biology happens at the protein-water interface, so all dynamical processes in this region are of fundamental importance. Local structural fluctuations in the hydration layer can be probed by ¹⁷O magnetic relaxation dispersion (MRD), which, at high frequencies, measures the integral of a biaxial rotational time correlation function (TCF)-the integral rotational correlation time. Numerous ¹⁷O MRD studies have demonstrated that this correlation time, when averaged over the first hydration shell, is longer than in bulk water by a factor 3-5. This rotational perturbation factor (RPF) has been corroborated by molecular dynamics simulations, which can also reveal the underlying molecular mechanisms. Here, we address several outstanding problems in this area by analyzing an extensive set of molecular dynamics data, including four globular proteins and three water models. The vexed issue of polarity versus topography as the primary determinant of hydration water dynamics is resolved by establishing a protein-invariant exponential dependence of the RPF on a simple confinement index. We conclude that the previously observed correlation of the RPF with surface polarity is a secondary effect of the correlation between polarity and confinement. Water rotation interpolates between a perturbed but bulk-like collective mechanism at low confinement and an exchange-mediated orientational randomization (EMOR) mechanism at high confinement. The EMOR process, which accounts for about half of the RPF, was not recognized in previous simulation studies, where only the early part of the TCF was examined. Based on the analysis of the experimentally relevant TCF over its full time course, we compare simulated and measured RPFs, finding a 30% discrepancy attributable to force field imperfections. We also compute the full ¹⁷O MRD profile, including the low-frequency dispersion produced by buried water molecules. Computing a local RPF for each hydration shell, we find that the perturbation decays exponentially with a decay "length" of 0.3 shells and that the second and higher shells account for a mere 3% of the total perturbation measured by ¹⁷O MRD. The only long-range effect is a weak water alignment in the electric field produced by an electroneutral protein (not screened by counterions), but this effect is negligibly small for ¹⁷O MRD. By contrast, we find that the ¹⁷O TCF is significantly more sensitive to the important short-range perturbations than the other two TCFs examined here.

Extensive Assessment of Various Computational Methods for Aspartate's pKa Shift

A series of computational methods for pK_a shift prediction are extensively tested on a set of benchmark protein systems, aiming at identifying pitfalls and evaluating their performance on high variants. Including 19 ASP residues in 10 protein systems, the benchmark set consists of both residues with highly shifted pK_a values as well as those varying little from the reference value, with an experimental RMS free energy differences of 2.49 kcal/mol with respect to blocked amino acid, namely the RMS pK_a shift being 1.82 pK_a units. The constant pH molecular dynamics (MD), alchemical methods, PROPKA3.1, and multiconformation continuum electrostatics give RMSDs of 1.52, 2.58, 1.37, and 3.52 pK_a units, respectively, on the benchmark set. The empirical scoring method is the most accurate one with extremely low computational cost, and the pH-dependent model is also able to provide accurate results, while the accuracy of MD sampling incorporating alchemical free energy simulation is prohibited by convergence achievement and the performance of conformational search incorporating multiconformation continuum electrostatics is bad. Former research works did not define statistical uncertainty with care and yielded the questionable conclusion that alchemical methods perform well in most benchmarks. In this work the traditional alchemical methods are thoroughly tested for high variants. We also performed the first application of nonequilibrium alchemical methods to the pK_a cases.

Mechanism of endonuclease cleavage by the HigB toxin

Bacteria encode multiple type II toxin-antitoxin modules that cleave ribosome-bound mRNAs in response to stress. All ribosome-dependent toxin family members structurally characterized to date adopt similar microbial RNase architectures despite possessing low sequence identities. Therefore, determining which residues are catalytically important in this specialized RNase family has been a challenge in the field. Structural studies of RelE and YoeB toxins bound to the ribosome provided significant insights but biochemical experiments with RelE were required to clearly demonstrate which residues are critical for acid-base catalysis of mRNA cleavage. Here, we solved an X-ray crystal structure of the wild-type, ribosome-dependent toxin HigB bound to the ribosome revealing potential catalytic residues proximal to the mRNA substrate. Using cell-based and biochemical assays, we further determined that HigB residues His54, Asp90, Tyr91 and His92 are critical for activity in vivo, while HigB H54A and Y91A variants have the largest effect on mRNA cleavage in vitro Comparison of X-ray crystal structures of two catalytically inactive HigB variants with 70S-HigB bound structures reveal that HigB active site residues undergo conformational rearrangements likely required for recognition of its mRNA substrate. These data support the emerging concept that ribosome-dependent toxins have diverse modes of mRNA recognition.

Molecular basis of ribosome recognition and mRNA hydrolysis by the E. coli YafQ toxin

Bacterial type II toxin-antitoxin modules are protein–protein complexes whose functions are finely tuned by rapidly changing environmental conditions. E. coli toxin YafQ is suppressed under steady state growth conditions by virtue of its interaction with its cognate antitoxin, DinJ. During stress, DinJ is proteolytically degraded and free YafQ halts translation by degrading ribosome-bound mRNA to slow growth until the stress has passed. Although structures of the ribosome with toxins RelE and YoeB have been solved, it is unclear what residues among ribosome-dependent toxins are essential for mediating both recognition of the ribosome and the mRNA substrate given their low sequence identities. Here we show that YafQ coordinates binding to the 70S ribosome via three surface-exposed patches of basic residues that we propose directly interact with 16S rRNA. We demonstrate that YafQ residues H50, H63, D67 and H87 participate in acid-base catalysis during mRNA hydrolysis and further show that H50 and H63 functionally complement as general bases to initiate the phosphodiester cleavage reaction. Moreover YafQ residue F91 likely plays an important role in mRNA positioning. In summary, our findings demonstrate the plasticity of ribosome-dependent toxin active site residues and further our understanding of which toxin residues are important for function.

Structural coupling between FKBP12 and buried water

Globular proteins often contain structurally well-resolved internal water molecules. Previously, we reported results from a molecular dynamics study that suggested that buried water (Wat3) may play a role in modulating the structure of the FK506 binding protein-12 (FKBP12) (Park and Saven, Proteins 2005; 60:450-463). In particular, simulations suggested that disrupting a hydrogen bond to Wat3 by mutating E60 to either A or Q would cause a structural perturbation involving the distant W59 side chain, which rotates to a new conformation in response to the mutation. This effectively remodels the ligand-binding pocket, as the side chain in the new conformation is likely to clash with bound FK506. To test whether the protein structure is in effect modulated by the binding of a buried water in the distance, we determined high-resolution (0.92-1.29 A) structures of wild-type FKBP12 and its two mutants (E60A, E60Q) by X-ray crystallography. The structures of mutant FKBP12 show that the ligand-binding pocket is indeed remodeled as predicted by the substitution at position 60, even though the water molecule does not directly interact with any of the amino acids of the binding pocket. Thus, these structures support the view that buried water molecules constitute an integral, noncovalent component of the protein structure. Additionally, this study provides an example in which predictions from molecular dynamics simulations are experimentally validated with atomic precision, thus showing that the structural features of protein-water interactions can be reliably modeled at a molecular level.

Temperature dependence of protein-hydration hydrodynamics by molecular dynamics simulations

The dynamics of water molecules near the protein surface are different from those of bulk water and influence the structure and dynamics of the protein itself. To elucidate the temperature dependence hydration dynamics of water molecules, we present results from the molecular dynamic simulation of the water molecules surrounding two proteins (Carboxypeptidase inhibitor and Ovomucoid) at seven different temperatures (T=273 to 303 K, in increments of 5 K). Translational diffusion coefficients of the surface water and bulk water molecules were estimated from 2 ns molecular dynamics simulation trajectories. Temperature dependence of the estimated bulk water diffusion closely reflects the experimental values, while hydration water diffusion is retarded significantly due to the protein. Protein surface induced scaling of translational dynamics of the hydration waters is uniform over the temperature range studied, suggesting the importance protein-water interactions.

Molecular basis of the recognition process: hydrogen-bonding patterns in the guanine primary recognition site of ribonuclease T1

Investigation of the intrinsic H-bonding pattern of the guanine complex with a sizable segment (from Asn43 to Glu46) of the primary recognition site (PRS) in RNase T1 at the B3LYP/6-311G(d,p) level of theory enables the electronic density characteristics of the H-bonding patterns of the guanine-PRS complexes to be identified. The perfect H-bonding pattern in the guanine recognition site is achieved through the guanine complex interactions with the large segment of the PRS. Two significant short H-bonds, O epsilon 1...HN1 and O epsilon 2...HN2, have been identified. The similar short H-bond distances found in the anionic GC- base pair and in this study suggest that the short hydrogen-bond distances may be characteristic of the multiple H-bonded anionic nucleobases. The H-bonding energy distribution, the geometric analysis of the H-bonding pattern, and the electron structure characteristics of the H-bonds in the guanine PRS of RNase T1 all suggest that the O epsilon 1...HN1 and O epsilon 2...HN2 side-chain H-bonds dominate the binding at the guanine recognition site of RNase T1. Also, the geometry evidence, the electron structure characteristics, and the properties of the bond critical points of the H-bonds reveal that the side-chain H-bonding and the main-chain H-bonding are mutually intensifying. Thus the positive cooperativity between Asn43 to Tyr45 and Glu46 is proposed.

Minimizing frustration by folding in an aqueous environment

Although life as we know it evolved in an aqueous medium, the properties of water are not completely understood. In this review, we focus on the role of water in guiding protein folding and stability. Specifically, we discuss the mechanisms of protein folding in an aqueous environment, the effects of water on the folding energy landscape as well as the transition state ensemble, and interactions of water with the folded state. We show that water cannot be viewed as a passive solvent, but rather, plays a very active role in the life of a protein.

Role of flexibility and polarity as determinants of the hydration of internal cavities and pockets in proteins

Molecular dynamics simulations of Staphylococcal nuclease and of 10 variants with internal polar or ionizable groups were performed to investigate systematically the molecular determinants of hydration of internal cavities and pockets in proteins. In contrast to apolar cavities in rigid carbon structures, such as nanotubes or buckeyballs, internal cavities in proteins that are large enough to house a few water molecules will most likely be dehydrated unless they contain a source of polarity. The water content in the protein interior can be modulated by the flexibility of protein elements that interact with water, which can impart positional disorder to water molecules, or bias the pattern of internal hydration that is stabilized. This might explain differences in the patterns of hydration observed in crystal structures obtained at cryogenic and room temperature conditions. The ability of molecular dynamics simulations to determine the most likely sites of water binding in internal pockets and cavities depends on its efficiency in sampling the hydration of internal sites and alternative protein and water conformations. This can be enhanced significantly by performing multiple molecular dynamics simulations as well as simulations started from different initial hydration states.

First structural investigation of the restriction ribonuclease RegB: NMR spectroscopic conditions, 13C/15N double-isotopic labelling and two-dimensional heteronuclear spectra

The bacteriophage T4 genome-encoded ribonuclease RegB is the unique well-defined restriction endoribonuclease. This protein cleaves with an almost absolute specificity its RNA substrate in the middle of the GGAG tetranucleotide mainly found in the Shine-Dalgarno sequence (required for the prokaryotic initiation of the translation). This protein has no significant homology to any known ribonuclease and its structure has never been investigated. The extreme toxicity of this ribonuclease prevents the expression of large quantities for structural studies. Here, we show that the toxicity of RegB can be bypassed by using the RegB H48A point mutant and explain why resolving the structure of this mutant is relevant. For nuclear magnetic resonance (NMR) purposes, we report the preparation of highly pure (13)C/(15)N double-labelled 1.2mM samples of RegB H48A using a high yield expression procedure in minimal medium (30 mg/L). We also present a set of solution conditions that maintain the concentrated samples of this protein stable for long periods at the NMR-required temperature. Finally, we present the first (1)H/(15)N and (1)H/(13)C two-dimensional NMR spectra of RegB H48A. These spectra show that the protein is folded and that the full structural analysis of RegB by NMR is feasible.

Dynamics of water molecules buried in cavities of apolipoprotein E studied by molecular dynamics simulations and continuum electrostatic calculations

Molecular dynamics (MD) simulations of several nanoseconds each were used to monitor the dynamic behavior of the five crystal water molecules buried in the interior of the N-terminal domain of apolipoprotein E. These crystal water molecules are fairly well conserved in several apolipoprotein E structures, suggesting that they are not an artifact of the crystal and that they may have a structural and/or functional role for the protein. All five buried crystal water molecules leave the protein interior in the course of the longest simulations and exchange with water molecules from the bulk. The free energies of binding evaluated from the electrostatic binding free energy computed using a continuum model and estimates of the binding entropy changes represent shallow minima. The corresponding calculated residence times of the buried water molecules range from tens of picoseconds to hundreds of nanoseconds, which denote rather short times as for buried water molecules. Several water exchanges monitored in the simulations show that water molecules along the exit/entrance pathway use a relay of H bonds primarily formed with charged residues which helps either the exit or the entrance from or into the buried site. The exit/entrance of water molecules from/into the sites is permitted essentially by local motions of, at most, two side chains, indicating that, in these cases, complex correlated atomic motions are not needed to open the buried site toward the surface of the protein. This provides a possible explanation for the short residence times.

Charge-charge interactions are the primary determinants of the pK values of the ionizable groups in Ribonuclease T1

Coulomb's law and a finite difference Poisson-Boltzmann based analysis are used to predict the pK values for 15 ionizable side chains (6 Asp, 6 Glu and 3 His) in ribonuclease T1. These predicted values are compared to the measured pK values to gain insight into the most important factors that influence the pK values of the ionizable groups in proteins. Charge-charge interactions are clearly the most important factor that determines the pK values of most ionizable groups in ribonuclease T1. However, pK values can be shifted by several pK units by the Born self energy associated with burying ionizable groups and by favorable intramolecular hydrogen bonding.

Protein-water interactions in a dynamic world

Protein-water interactions are key to biological function. They have an underlying dynamic component that pervades the functional roles associated both with particular systems and with the properties of proteins in general. This article focuses on the specific ways in which the dynamics of water are important to protein structure, motion and adaptability to changes in the protein environment.

Magnetic relaxation dispersion studies of biomolecular solutions

Although the MRD method has a long record in biomolecular systems, it has undergone a renaissance in the past few years as methodological developments have provided access to new types of information. In particular, MRD studies of quadrupolar nuclei such as 17O and 23Na have yielded valuable insights about the interactions of proteins and oligonucleotides with their solvent environment. The biomolecular MRD literature is still dominated by hydration studies, but the method has also been used to study the interaction of organic cosolvents and inorganic counterions with biomolecules. The MRD method can potentially make important contributions to the understanding of the mechanisms whereby protein conformational stability is affected by nonaqueous solvent components, such as denaturants, stabilizers, and helix promoters. Residence times of water molecules and other low molecular weight species in association with biomolecules can be determined by MRD. Such residence times are of general interest for understanding the kinetics of biomolecule-ligand interactions and, when exchange is gated by the biomolecule, can be used to characterize large-scale conformational fluctuations on nanosecond-millisecond time scales. By monitoring the integrity and specific internal hydration sites as well as the global solvent exposure, the MRD method can also shed light on the structure and dynamics of biomolecules in fluctuating nonnative states. Because it does not rely on high resolution, the MRD method is also applicable to very large biomolecules and complexes and has even been used to investigate protein crystals, gels, and biological tissues. In fact, dynamic studies of solids and liquid crystals were among the earliest applications of the MRD method. In many of its diverse applications, the MRD method provides unique information, complementing that available from high-resolution NMR.

The contribution of metal ions to the conformational stability of ribonuclease T1: crystal versus solution

In the crystalline state, ribonuclease T1 binds calcium ions at different lattice-dependent positions. In solution, its conformational stability is also remarkably increased in the presence of divalent metal ions. Combining urea unfolding studies and X-ray crystallography, we compared the presence of several metal ions at specific sites in the protein to their contribution to the overall stabilizing effect in solution. We constructed thermodynamic cycles involving particular metal ions and specific carboxylate functions. The resulting coupling energies indicate that some (but not all) metal ions found at lattice contacts in crystal structures may indeed significantly contribute to stability enhancement in the presence of metal ions in solution.

Structural and thermodynamic analysis of the binding of solvent at internal sites in T4 lysozyme

To investigate the structural and thermodynamic basis of the binding of solvent at internal sites within proteins a number of mutations were constructed in T4 lysozyme. Some of these were designed to introduce new solvent-binding sites. Others were intended to displace solvent from preexisting sites. In one case Val-149 was replaced with alanine, serine, cysteine, threonine, isoleucine, and glycine. Crystallographic analysis shows that, with the exception of isoleucine, each of these substitutions results in the binding of solvent at a polar site that is sterically blocked in the wild-type enzyme. Mutations designed to perturb or displace a solvent molecule present in the native enzyme included the replacement of Thr-152 with alanine, serine, cysteine, valine, and isoleucine. Although the solvent molecule was moved in some cases by up to 1.7 A, in no case was it completely removed from the folded protein. The results suggest that hydrogen bonds from the protein to bound solvent are energy neutral. The binding of solvent to internal sites within proteins also appears to be energy neutral except insofar as the bound solvent may prevent a loss of energy due to potential hydrogen bonding groups that would otherwise be unsatisfied. The introduction of a solvent-binding site appears to require not only a cavity to accommodate the water molecule but also the presence of polar groups to help satisfy its hydrogen-bonding potential. It may be easier to design a site to accommodate two or more water molecules rather than one as the solvent molecules can then hydrogen-bond to each other. For similar reasons it is often difficult to design a point mutation that will displace a single solvent molecule from the core of a protein.

Protein-bound water molecule counting by resolution of (1)H spin-lattice relaxation mechanisms

Water proton spin-lattice relaxation is studied in dilute solutions of bovine serum albumin as a function of magnetic field strength, oxygen concentration, and solvent deuteration. In contrast to previous studies conducted at high protein concentrations, the observed relaxation dispersion is accurately Lorentzian with an effective correlation time of 41 +/- 3 ns when measured at low proton and low protein concentrations to minimize protein aggregation. Elimination of oxygen flattens the relaxation dispersion profile above the rotational inflection frequency, nearly eliminating the high field tail previously attributed to a distribution of exchange times for either whole water molecules or individual protons at the protein-water interface. The small high-field dispersion that remains is attributed to motion of the bound water molecules on the protein or to internal protein motions on a time scale of order one ns. Measurements as a function of isotope composition permit separation of intramolecular and intermolecular relaxation contributions. The magnitude of the intramolecular proton-proton relaxation rate constant is interpreted in terms of 25 +/- 4 water molecules that are bound rigidly to the protein for a time long compared with the rotational correlation time of 42 ns. This number of bound water molecules neglects the possibility of local motions of the water in the binding site; inclusion of these effects may increase the number of bound water molecules by 50%.

A "structural" water molecule in the family of fatty acid binding proteins

A single water molecule (w135), buried within the structure of rat intestinal fatty acid binding protein (I-FABP), is investigated by NMR, molecular dynamics simulations, and analysis of known crystal structures. An ordered water molecule was found in structurally analogous position in 24 crystal structures of nine different members of the family of fatty acid binding proteins. There is a remarkable conservation of the local structure near the w135 binding site among different proteins from this family. NMR cross-relaxation measurements imply that w135 is present in the I-FABP:ANS (1-sulfonato-8-(1')anilinonaphthalene) complex in solution with the residence time of >300 ps. Mean-square positional fluctuations of w135 oxygen observed in MD simulations (0.18 and 0.13 A2) are comparable in magnitude to fluctuations exhibited by the backbone atoms and result from highly constrained binding pocket as revealed by Voronoi volumes (averages of 27.0 +/- 1.8 A3 and 24.7 +/- 2.2 A3 for the two simulations). Escape of w135 from its binding pocket was observed only in one MD simulation. The escape process was initiated by interactions with external water molecules and was accompanied by large deformations in beta-strands D and E. Immediately before the release, w135 assumed three distinct states that differ in hydrogen bonding topology and persisted for about 15 ps each. Computer simulations suggest that escape of w135 from the I-FABP matrix is primarily determined by conformational fluctuations of the protein backbone and interactions with external water molecules.

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.