His...Asp catalytic dyad of ribonuclease A: structure and function of the wild-type, D121N, and D121A enzymes

An active site mutation induces oxygen reactivity in D-arginine dehydrogenase: A case of superoxide diverting protons

Enzymes are potent catalysts that increase biochemical reaction rates by several orders of magnitude. Flavoproteins are a class of enzymes whose classification relies on their ability to react with molecular oxygen (O2) during catalysis using ionizable active site residues. Pseudomonas aeruginosa D-arginine dehydrogenase (PaDADH) is a flavoprotein that oxidizes D-arginine for P. aeruginosa survival and biofilm formation. The crystal structure of PaDADH reveals the interaction of the glutamate 246 (E246) side chain with the substrate and at least three other active site residues, establishing a hydrogen bond network in the active site. Additionally, E246 likely ionizes to facilitate substrate binding during PaDADH catalysis. This study aimed to investigate how replacing the E246 residue with leucine affects PaDADH catalysis and its ability to react with O2 using steady-state kinetics coupled with pH profile studies. The data reveal a gain of O2 reactivity in the E246L variant, resulting in a reduced flavin semiquinone species and superoxide (O2•-) during substrate oxidation. The O2•- reacts with active site protons, resulting in an observed nonstoichiometric slope of 1.5 in the enzyme's log (kcat/Km) pH profile with D-arginine. Adding superoxide dismutase results in an observed correction of the slope to 1.0. This study demonstrates how O2•- can alter the slopes of limbs in the pH profiles of flavin-dependent enzymes and serves as a model for correcting nonstoichiometric slopes in elucidating reaction mechanisms of flavoproteins.

Crystal structure and activity of a de novo enzyme, ferric enterobactin esterase Syn-F4

Producing novel enzymes that are catalytically active in vitro and biologically functional in vivo is a key goal of synthetic biology. Previously, we reported Syn-F4, the first de novo protein that meets both criteria. Syn-F4 hydrolyzed the siderophore ferric enterobactin, and expression of Syn-F4 allowed an inviable strain of Escherichia coli (Δfes) to grow in iron-limited medium. Here, we describe the crystal structure of Syn-F4. Syn-F4 forms a dimeric 4-helix bundle. Each monomer comprises two long α-helices, and the loops of the Syn-F4 dimer are on the same end of the bundle (syn topology). Interestingly, there is a penetrated hole in the central region of the Syn-F4 structure. Extensive mutagenesis experiments in a previous study showed that five residues (Glu26, His74, Arg77, Lys78, and Arg85) were essential for enzymatic activity in vivo. All these residues are located around the hole in the central region of the Syn-F4 structure, suggesting a putative active site with a catalytic dyad (Glu26-His74). The complete inactivity of purified proteins with mutations at the five residues supports the putative active site and reaction mechanism. Molecular dynamics and docking simulations of the ferric enterobactin siderophore binding to the Syn-F4 structure demonstrate the dynamic property of the putative active site. The structure and active site of Syn-F4 are completely different from native enterobactin esterase enzymes, thereby demonstrating that proteins designed de novo can provide life-sustaining catalytic activities using structures and mechanisms dramatically different from those that arose in nature.

A method for quantifying how the activity of an enzyme is affected by the net charge of its nearest crowded neighbor

The electrostatic effects of protein crowding have not been systematically explored. Rather, protein crowding is generally studied with co‐solvents or crowders that are electrostatically neutral, with no methods to measure how the net charge (Z) of a crowder affects protein function. For example, can the activity of an enzyme be affected electrostatically by the net charge of its neighbor in crowded milieu? This paper reports a method for crowding proteins of different net charge to an enzyme via semi‐random chemical crosslinking. As a proof of concept, RNase A was crowded (at distances ≤ the Debye length) via crosslinking to different heme proteins with Z = +8.50 ± 0.04, Z = +6.39 ± 0.12, or Z = −10.30 ± 1.32. Crosslinking did not disrupt the structure of proteins, according to amide H/D exchange, and did not inhibit RNase A activity. For RNase A, we found that the electrostatic environment of each crowded neighbor had significant effects on rates of RNA hydrolysis. Crowding with cationic cytochrome c led to increases in activity, while crowding with anionic “supercharged” cytochrome c or myoglobin diminished activity. Surprisingly, electrostatic crowding effects were amplified at high ionic strength (I = 0.201 M) and attenuated at low ionic strength (I = 0.011 M). This salt dependence might be caused by a unique set of electric double layers at the dimer interspace (maximum distance of 8 Å, which cannot accommodate four layers). This new method of crowding via crosslinking can be used to search for electrostatic effects in protein crowding.

Not making the cut: Techniques to prevent RNA cleavage in structural studies of RNase-RNA complexes

RNases are varied in the RNA structures and sequences they target for cleavage and are an important type of enzyme in cells. Despite the numerous examples of RNases known, and of those with determined three-dimensional structures, relatively few examples exist with the RNase bound to intact cognate RNA substrate prior to cleavage. To better understand RNase structure and sequence specificity for RNA targets, in vitro methods used to assemble these enzyme complexes trapped in a pre-cleaved state have been developed for a number of different RNases. We have surveyed the Protein Data Bank for such structures and in this review detail methodologies that have successfully been used and relate them to the corresponding structures. We also offer ideas and suggestions for future method development. Many strategies within this review can be used in combination with X-ray crystallography, as well as cryo-EM, and other structure-solving techniques. Our hope is that this review will be used as a guide to resolve future yet-to-be-determined RNase-substrate complex structures.

Theoretical Studies on the Effect of Isomerized Aspartic Acid Residues on the Three-Dimensional Structures of Bovine Pancreatic Ribonucleases A

Isomerized aspartic acid (Asp) residues have previously been identified in various aging tissues, and are suspected to contribute to age-related diseases. Asp-residue isomerization occurs nonenzymatically under physiological conditions, resulting in the formation of three types of isomerized Asp (i.e., L-isoAsp, D-Asp, and D-isoAsp) residues. Asp-residue isomerization often accelerates protein aggregation and insolubilization, making structural biology analyses difficult. Recently, Sakaue et al. reported the synthesis of a ribonuclease A (RNase A) in which Asp121 was artificially replaced with different isomerized Asp residues, and experimentally demonstrated that the enzymatic activities of these artificial mutants were completely lost. However, their structural features have not yet been elucidated. In the present study, the three-dimensional (3D) structures of these artificial-mutant RNases A were predicted using molecular dynamics (MD) simulations. The 3D structures of wild-type and artificial-mutant RNases A were converged by 3000-ns MD simulations. Our computational data show that the structures of the active site and the formation frequencies of the appropriate catalytic dyad structures in the artificial-mutant RNases A were quite different from wild-type RNase A. These computational findings may provide an explanation for the experimental data which show that artificial-mutant RNases A lack enzymatic activity. Herein, MD simulations have been used to evaluate the influences of isomerized Asp residues on the 3D structures of proteins.

Crystal structures of non-oxidative decarboxylases reveal a new mechanism of action with a catalytic dyad and structural twists

Hydroxybenzoic acids, like gallic acid and protocatechuic acid, are highly abundant natural compounds. In biotechnology, they serve as critical precursors for various molecules in heterologous production pathways, but a major bottleneck is these acids’ non-oxidative decarboxylation to hydroxybenzenes. Optimizing this step by pathway and enzyme engineering is tedious, partly because of the complicating cofactor dependencies of the commonly used prFMN-dependent decarboxylases. Here, we report the crystal structures (1.5–1.9 Å) of two homologous fungal decarboxylases, AGDC1 from Arxula adenivorans, and PPP2 from Madurella mycetomatis. Remarkably, both decarboxylases are cofactor independent and are superior to prFMN-dependent decarboxylases when heterologously expressed in Saccharomyces cerevisiae. The organization of their active site, together with mutational studies, suggests a novel decarboxylation mechanism that combines acid–base catalysis and transition state stabilization. Both enzymes are trimers, with a central potassium binding site. In each monomer, potassium introduces a local twist in a β-sheet close to the active site, which primes the critical H86-D40 dyad for catalysis. A conserved pair of tryptophans, W35 and W61, acts like a clamp that destabilizes the substrate by twisting its carboxyl group relative to the phenol moiety. These findings reveal AGDC1 and PPP2 as founding members of a so far overlooked group of cofactor independent decarboxylases and suggest strategies to engineer their unique chemistry for a wide variety of biotechnological applications.

Covalent binding of uracil DNA glycosylase UdgX to abasic DNA upon uracil excision

Uracil DNA glycosylases (UDGs) are important DNA repair enzymes that excise uracil from DNA, yielding an abasic site. Recently, UdgX, an unconventional UDG with extremely tight binding to DNA containing uracil, was discovered. The structure of UdgX from Mycobacterium smegmatis in complex with DNA shows an overall similarity to that of family 4 UDGs except for a protruding loop at the entrance of the uracil-binding pocket. Surprisingly, H109 in the loop was found to make a covalent bond to the abasic site to form a stable intermediate, while the excised uracil remained in the pocket of the active site. H109 functions as a nucleophile to attack the oxocarbenium ion, substituting for the catalytic water molecule found in other UDGs. To our knowledge, this change from a catalytic water attack to a direct nucleophilic attack by the histidine residue is unprecedented. UdgX utilizes a unique mechanism of protecting cytotoxic abasic sites from exposure to the cellular environment.

Enhanced Azidolysis by the Formation of Stable Ser-His Catalytic Dyad in a Glycoside Hydrolase Family 10 Xylanase Mutant

Glycoside hydrolases require carboxyl groups as catalysts for their activity. A retaining xylanase from Streptomyces olivaceoviridis E-86 belonging to glycoside hydrolase family 10 possesses Glu128 and Glu236 that respectively function as acid/base and nucleophile. We previously developed a unique mutant of the retaining xylanase, N127S/E128H, whose deglycosylation is triggered by azide. A crystallographic study reported that the transient formation of a Ser–His catalytic dyad in the reaction cycle possibly reduced the azidolysis reaction. In the present study, we engineered a catalytic dyad with enhanced stability by site-directed mutagenesis and crystallographic study of N127S/E128H. Comparison of the Michaelis complexes of N127S/E128H with pNP-X₂ and with xylopentaose showed that Ser127 could form an alternative hydrogen bond with Thr82, which disrupts the formation of the Ser–His catalytic dyad. The introduction of T82A mutation in N127S/E128H produces an enhanced first-order rate constant (6 times that of N127S/E128H). We confirmed the presence of a stable Ser–His hydrogen bond in the Michaelis complex of the triple mutant, which forms the productive tautomer of His128 that acts as an acid catalyst. Because the glycosyl azide is applicable in the bioconjugation of glycans by using click chemistry, the enzyme-assisted production of the glycosyl azide may contribute to the field of glycobiology.

Tautomeric stabilities of 4-fluorohistidine shed new light on mechanistic experiments with labeled ribonuclease A

Ribonuclease A is the oldest model for studying enzymatic mechanisms, yet questions remain about proton transfer within the active site. Seminal work by Jackson et al. (Science, 1994) labeled Ribonuclease A with 4-fluorohistidine, concluding that active-site histidines act as general acids and bases. Calculations of 4-fluorohistidine indicate that the π-tautomer is predominant in all simulated environments (by ~17 kJ/mol), strongly suggesting that fluoro-labeled ribonuclease A functions with His119 in π-tautomer. The tautomeric form of His119 during proton transfer and tautomerism as a putative mechanistic step in wild-type RNase A remain open questions and should be considered in future mechanistic studies.

The first crystal structure of human RNase 6 reveals a novel substrate-binding and cleavage site arrangement

We describe the first human RNase 6 crystal structure in complex with sulfate anions. Kinetic analysis, site-directed mutagenesis and molecular dynamics simulations identified novel substrate recognition and cleavage sites.

Human RNase 6 is a cationic secreted protein that belongs to the RNase A superfamily. Its expression is induced in neutrophils and monocytes upon bacterial infection, suggesting a role in host defence. We present here the crystal structure of RNase 6 obtained at 1.72 Å (1 Å=0.1 nm) resolution, which is the first report for the protein 3D structure and thereby setting the basis for functional studies. The structure shows an overall kidney-shaped globular fold shared with the other known family members. Three sulfate anions bound to RNase 6 were found, interacting with residues at the main active site (His¹⁵, His¹²² and Gln¹⁴) and cationic surface-exposed residues (His³⁶, His³⁹, Arg⁶⁶ and His⁶⁷). Kinetic characterization, together with prediction of protein–nucleotide complexes by molecular dynamics, was applied to analyse the RNase 6 substrate nitrogenous base and phosphate selectivity. Our results reveal that, although RNase 6 is a moderate catalyst in comparison with the pancreatic RNase type, its structure includes lineage-specific features that facilitate its activity towards polymeric nucleotide substrates. In particular, enzyme interactions at the substrate 5′ end can provide an endonuclease-type cleavage pattern. Interestingly, the RNase 6 crystal structure revealed a novel secondary active site conformed by the His³⁶–His³⁹ dyad that facilitates the polynucleotide substrate catalysis.

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.

Structures of the Apo and FAD-bound forms of 2-hydroxybiphenyl 3-monooxygenase (HbpA) locate activity hotspots identified by using directed evolution

The FAD-dependent monooxygenase HbpA from Pseudomonas azelaica HBP1 catalyses the hydroxylation of 2-hydroxybiphenyl (2HBP) to 2,3-dihydroxybiphenyl (23DHBP). HbpA has been used extensively as a model for studying flavoprotein hydroxylases under process conditions, and has also been subjected to directed-evolution experiments that altered its catalytic properties. The structure of HbpA has been determined in its apo and FAD-complex forms to resolutions of 2.76 and 2.03 Å, respectively. Comparisons of the HbpA structure with those of homologues, in conjunction with a model of the reaction product in the active site, reveal His48 as the most likely acid/base residue to be involved in the hydroxylation mechanism. Mutation of His48 to Ala resulted in an inactive enzyme. The structures of HbpA also provide evidence that mutants achieved by directed evolution that altered activity are comparatively remote from the substrate-binding site.

A catalyst selection protocol that identifies biomimetic motifs from β-hairpin libraries

Assaying a solid-phase library of histidine-containing β-hairpin peptides by a reactive tagging scheme in organic solvents selects for catalysts that reproduce the strategies used by His-based enzyme active sites to accelerate acyl- and phosphonyl-transfer reactions. Rate accelerations (k(rel)) in organic solvents of up to 2.4 × 10(8) are observed.

What's in your buffer? Solute altered millisecond motions detected by solution NMR

To date, little work has been conducted on the relationship between solute and buffer molecules and conformational exchange motion in enzymes. This study uses solution NMR to examine the effects of phosphate, sulfate, and acetate in comparison to MES- and HEPES-buffered references on the chemical shift perturbation and millisecond, chemical, or conformational exchange motions in the enzyme ribonuclease A (RNase A), triosephosphate isomerase (TIM) and HisF. The results indicate that addition of these solutes has a small effect on (1)H and (15)N chemical shifts for RNase A and TIM but a significant effect for HisF. For RNase A and TIM, Carr-Purcell-Meiboom-Gill relaxation dispersion experiments, however, show significant solute-dependent changes in conformational exchange motions. Some residues show loss of millisecond motions relative to the reference sample upon addition of solute, whereas others experience an enhancement. Comparison of exchange parameters obtained from fits of dispersion data indicates changes in either or both equilibrium populations and chemical shifts between conformations. Furthermore, the exchange kinetics are altered in many cases. The results demonstrate that common solute molecules can alter observed enzyme millisecond motions and play a more active role than what is routinely believed.

Structural and mechanistic insights into guanylylation of RNA-splicing ligase RtcB joining RNA between 3'-terminal phosphate and 5'-OH

The RtcB protein has recently been identified as a 3'-phosphate RNA ligase that directly joins an RNA strand ending with a 2',3'-cyclic phosphate to the 5'-hydroxyl group of another RNA strand in a GTP/Mn(2+)-dependent reaction. Here, we report two crystal structures of Pyrococcus horikoshii RNA-splicing ligase RtcB in complex with Mn(2+) alone (RtcB/ Mn(2+)) and together with a covalently bound GMP (RtcB-GMP/Mn(2+)). The RtcB/ Mn(2+) structure (at 1.6 Å resolution) shows two Mn(2+) ions at the active site, and an array of sulfate ions nearby that indicate the binding sites of the RNA phosphate backbone. The structure of the RtcB-GMP/Mn(2+) complex (at 2.3 Å resolution) reveals the detailed geometry of guanylylation of histidine 404. The critical roles of the key residues involved in the binding of the two Mn(2+) ions, the four sulfates, and GMP are validated in extensive mutagenesis and biochemical experiments, which also provide a thorough characterization for the three steps of the RtcB ligation pathway: (i) guanylylation of the enzyme, (ii) guanylyl-transfer to the RNA substrate, and (iii) overall ligation. These results demonstrate that the enzyme's substrate-induced GTP binding site and the putative reactive RNA ends are in the vicinity of the binuclear Mn(2+) active center, which provides detailed insight into how the enzyme-bound GMP is tansferred to the 3'-phosphate of the RNA substrate for activation and subsequent nucleophilic attack by the 5'-hydroxyl of the second RNA substrate, resulting in the ligated product and release of GMP.

Boronate-mediated biologic delivery

Inefficient cellular delivery limits the landscape of macromolecular drugs. Boronic acids readily form boronate esters with the 1,2- and 1,3-diols of saccharides, such as those that coat the surface of mammalian cells. Here pendant boronic acids are shown to enhance the cytosolic delivery of a protein toxin. Thus, boronates are a noncationic carrier that can deliver a polar macromolecule into mammalian cells.

On the information expressed in enzyme structure: more lessons from ribonuclease A

Brownian computations were directed at Ribonuclease A (RNase A) and variants in folded states so as to quantify information of the statistical type at the atom/covalent bond level. This advanced the research reported in this journal last year on the information properties of enzyme primary structure. Brownian computation data are illustrated for a sixteen-member library. The results identify signature traits that distinguish the folded wild type (WT) molecule from variants. The distinctions are explainable in terms of correlated information and dispersion energy. The Brownian tools used for this study can be directed at other protein families (e.g., kinases, isomerases, etc.) in rapid screening, QSAR, and design applications.

A new RNase sheds light on the RNase/angiogenin subfamily from zebrafish

Recently, extracellular RNases of the RNase A superfamily, with the characteristic CKxxNTF sequence signature, have been identified in fish. This has led to the recognition that these RNases are present in the whole vertebrate subphylum. In fact, they comprise the only enzyme family unique to vertebrates. Four RNases from zebrafish (Danio rerio) have been previously reported and have a very low RNase activity; some of these are endowed, like human angiogenin, with powerful angiogenic and bactericidal activities. In the present paper, we report the three-dimensional structure, the thermodynamic behaviour and the biological properties of a novel zebrafish RNase, ZF-RNase-5. The investigation of its structural and functional properties, extended to all other subfamily members, provides an inclusive description of the whole zebrafish RNase subfamily.

Structural and kinetic analysis of Bacillus subtilis N-acetylglucosaminidase reveals a unique Asp-His dyad mechanism

Three-dimensional structures of NagZ of Bacillus subtilis, the first structures of a two-domain β-N-acetylglucosaminidase of family 3 of glycosidases, were determined with and without the transition state mimicking inhibitor PUGNAc bound to the active site, at 1.84- and 1.40-Å resolution, respectively. The structures together with kinetic analyses of mutants revealed an Asp-His dyad involved in catalysis: His(234) of BsNagZ acts as general acid/base catalyst and is hydrogen bonded by Asp(232) for proper function. Replacement of both His(234) and Asp(232) with glycine reduced the rate of hydrolysis of the fluorogenic substrate 4'-methylumbelliferyl N-acetyl-β-D-glucosaminide 1900- and 4500-fold, respectively, and rendered activity pH-independent in the alkaline range consistent with a role of these residues in acid/base catalysis. N-Acetylglucosaminyl enzyme intermediate accumulated in the H234G mutant and β-azide product was formed in the presence of sodium azide in both mutants. The Asp-His dyad is conserved within β-N-acetylglucosaminidases but otherwise absent in β-glycosidases of family 3, which instead carry a "classical" glutamate acid/base catalyst. The acid/base glutamate of Hordeum vulgare exoglucanase (Exo1) superimposes with His(234) of the dyad of BsNagZ and, in contrast to the latter, protrudes from a second domain of the enzyme into the active site. This is the first report of an Asp-His catalytic dyad involved in hydrolysis of glycosides resembling in function the Asp-His-Ser triad of serine proteases. Our findings will facilitate the development of mechanism-based inhibitors that selectively target family 3 β-N-acetylglucosaminidases, which are involved in bacterial cell wall turnover, spore germination, and induction of β-lactamase.

Investigation of the catalytic mechanism of a synthetic DNAzyme with protein-like functionality: an RNaseA mimic?

The protein enzyme ribonuclease A (RNaseA) cleaves RNA with catalytic perfection, although with little sequence specificity, by a divalent metal ion (M(2+))-independent mechanism in which a pair of imidazoles provides general acid and base catalysis, while a cationic amine provides electrostatic stabilization of the transition state. Synthetic imitation of this remarkable organo-catalyst ("RNaseA mimicry") has been a longstanding goal in biomimetic chemistry. The 9(25)-11 DNAzyme contains synthetically modified nucleotides presenting both imidazole and cationic amine side chains, and catalyzes RNA cleavage with turnover in the absence of M(2+) similarly to RNaseA. Nevertheless, the catalytic roles, if any, of the "protein-like" functional groups have not been defined, and hence the question remains whether 9(25)-11 engages any of these functionalities to mimic aspects of the mechanism of RNaseA. To address this question, we report a mechanistic investigation of 9(25)-11 catalysis wherein we have employed a variety of experiments, such as DNAzyme functional group deletion, mechanism-based affinity labeling, and bridging and nonbridging phosphorothioate substitution of the scissile phosphate. Several striking parallels exist between the results presented here for 9(25)-11 and the results of analogous experiments applied previously to RNaseA. Specifically, our results implicate two particular imidazoles in general acid and base catalysis and suggest that a specific cationic amine stabilizes the transition state via diastereoselective interaction with the scissile phosphate. Overall, 9(25)-11 appears to meet the minimal criteria of an RNaseA mimic; this demonstrates how added synthetic functionality can expand the mechanistic repertoire available to a synthetic DNA-based catalyst.

Mutation analysis in UGT1A9 suggests a relationship between substrate and catalytic residues in UDP-glucuronosyltransferases

UDP-glucuronosyltransferases (UGTs) catalyze the transfer of glucuronic acid from UDP-glucuronic acid to endo- and xenobiotics in our body. UGTs belong to the GT1 family of glycosyltransferases and many GT1s use a serine protease-like catalytic mechanism in which an Asp-His pair deprotonates a hydroxyl on the aglycone for nucleophilic attack on the sugar donor. The pair in human UGTs could be H37 and either D143 or D148 (UGT1A9 numbering). However, H37 is not totally conserved, being replaced by either Pro or Leu in UGT1A4 and UGT2B10. We therefore investigated the role of H37, D143 and D148 in UGT1A9 by site-directed mutagenesis, activity and kinetic measurements with several substrates. The results suggest that H37 is not critical in N-glucuronidation, but is so in O-glucuronidation. The V(max) of the H37A mutant was much less affected in N- than O-glucuronidation, while the reverse was true for the Asp mutations, particularly D143A. We suggest that this is due to the opposing properties of O- and N- nucleophiles. O-nucleophiles require the histidine to deprotonate them so that they become effective nucleophiles, while N-nucleophiles develop a formal positive charge during the reaction (RNH(2)(+)-GlcA), and thus require a negatively charged residue to stabilize the transition state.

Ribonuclease A homologues of the zebrafish: polymorphism, crystal structures of two representatives and their evolutionary implications

The widespread and functionally varied members of the ribonuclease A (RNase A) superfamily provide an excellent opportunity to study evolutionary forces at work on a conserved protein scaffold. Representatives from the zebrafish are of particular interest as the evolutionary distance from non-ichthyic homologues is large. We conducted an exhaustive survey of available zebrafish DNA sequences and found significant polymorphism among its four known homologues. In an extension of previous nomenclature, the variants have been named RNases ZF-1a–c,-2a–d,-3a–e and-4. We present the first X-ray crystal structures of zebrafish ribonucleases, RNases ZF-1a and-3e at 1.35-and 1.85 Å resolution, respectively. Structure-based clustering with ten other ribonuclease structures indicates greatest similarity to mammalian angiogenins and amphibian ribonucleases, and supports the view that all present-day ribonucleases evolved from a progenitor with three disulphide bonds. In their details, the two structures are intriguing melting-pots of features present in ribonucleases from other vertebrate classes. Whereas in RNase ZF-1a the active site is obstructed by the C-terminal segment (as observed in angiogenin), in RNase ZF-3e the same region is open (as observed in more catalytically efficient homologues). The progenitor of present-day ribonucleases is more likely to have had an obstructive C terminus, and the relatively high similarity (late divergence) of RNases ZF-1 and-3 infers that the active site unblocking event has happened independently in different vertebrate lineages.

Enzymatic and Structural Characterisation of Amphinase, a Novel Cytotoxic Ribonuclease from Rana pipiens Oocytes

Besides Onconase (ONC) and its V11/N20/R103-variant, oocytes of the Northern Leopard frog (Rana pipiens) contain another homologue of ribonuclease A, which we named Amphinase (Amph). Four variants (Amph-1-4) were isolated and sequenced, each 114 amino acid residues in length and N-glycosylated at two positions. Sequence identities (a) among the variants and (b) versus ONC are 86.8-99.1% and 38.2-40.0%, respectively. When compared with other amphibian ribonucleases, a typical pattern of cysteine residues is evident but the N-terminal pyroglutamate residue is replaced by a six-residue extension. Amph variants have relatively weak ribonucleolytic activity that is insensitive to human ribonuclease inhibitor protein (RI). Values of k(cat)/K(M) with hypersensitive fluorogenic substrates are 10(4) and 10(2)-fold lower than the maximum values exhibited by ribonuclease A and ONC, respectively, and there is little cytosine/uracil or adenine/guanine discrimination at the B(1) or B(2) subsites, respectively. Amph variants have cytotoxic activity toward A-253 carcinoma cells that requires intact ribonucleolytic activity. The glycan component has little or no influence over single-stranded RNA cleavage, RI evasion or cytotoxicity. The crystal structures of natural and recombinant Amph-2 (determined at 1.8 and 1.9 A resolution, respectively) reveal that the N terminus is unlikely to play a catalytic role (but an unusual alpha2-beta1 loop may do so) and the B(2) subsite is rudimentary. At the active site, structural features that may contribute to the enzyme's low ribonucleolytic activity are the fixture of Lys14 in an obstructive position, the accompanying ejection of Lys42, and a lack of constraints on the conformations of Lys42 and His107.

The mechanism of rate-limiting motions in enzyme function

The ability to use conformational flexibility is a hallmark of enzyme function. Here we show that protein motions and catalytic activity in a RNase are coupled and display identical solvent isotope effects. Solution NMR relaxation experiments identify a cluster of residues, some distant from the active site, that are integral to this motion. These studies implicate a single residue, histidine-48, as the key modulator in coupling protein motion with enzyme function. Mutation of H48 to alanine results in loss of protein motion in the isotope-sensitive region of the enzyme. In addition, k(cat) decreases for this mutant and the kinetic solvent isotope effect on k(cat), which was 2.0 in WT, is near unity in H48A. Despite being located 18 A from the enzyme active site, H48 is essential in coordinating the motions involved in the rate-limiting enzymatic step. These studies have identified, of approximately 160 potential exchangeable protons, a single site that is integral in the rate-limiting step in RNase A enzyme function.

The binding of 3′-N-piperidine-4-carboxyl-3′-deoxy-ara-uridine to ribonuclease A in the crystal

The binding of a moderate inhibitor, 3'-N-piperidine-4-carboxyl-3'-deoxy-ara-uridine, to ribonuclease A has been studied by X-ray crystallography at 1.7A resolution. Two inhibitor molecules are bound in the central RNA binding cavity of RNase A exploiting interactions with residues from peripheral binding sites rather than from the active site of the enzyme. The uracyl moiety of the first inhibitor molecule occupies the purine-preferring site of RNase A, while the rest of the molecule projects to the solvent. The second inhibitor molecule binds with the carboxyl group at the pyrimidine recognition site and the uridine moiety exploits interactions with RNase A residues Lys66, His119 and Asp121. Comparative structural analysis of the 3'-N-piperidine-4-carboxyl-3'-deoxy-ara-uridine complex with other RNase A-ligand complexes provides a structural explanation of its potency. The crystal structure of the RNase A-3'-N-piperidine-4-carboxyl-3'-deoxy-ara-uridine complex provides evidence of a novel ligand-binding pattern in RNase A for 3'-N-aminonucleosides that was not anticipated by modelling studies, while it also suggests ways to improve the efficiency and selectivity of such compounds to develop pharmaceuticals against pathologies associated with RNase A homologues.

Genetic selection for critical residues in ribonucleases

Homologous mammalian proteins were subjected to an exhaustive search for residues that are critical to their structure/function. Error-prone polymerase chain reactions were used to generate random mutations in the genes of bovine pancreatic ribonuclease (RNase A) and human angiogenin, and a genetic selection based on the intrinsic cytotoxicity of ribonucleolytic activity was used to isolate inactive variants. Twenty-three of the 124 residues in RNase A were found to be intolerant to substitution with at least one particular amino acid. Twenty-nine of the 123 residues in angiogenin were likewise intolerant. In both RNase A and angiogenin, only six residues appeared to be wholly intolerant to substitution: two histidine residues involved in general acid/base catalysis and four cysteine residues that form two disulfide bonds. With few exceptions, the remaining critical residues were buried in the hydrophobic core of the proteins. Most of these residues were found to tolerate only conservative substitutions. The importance of a particular residue as revealed by this genetic selection correlated with its sequence conservation, though several non-conserved residues were found to be critical for protein structure/function. Despite voluminous research on RNase A, the importance of many residues identified herein was unknown, and those can now serve as targets for future work. Moreover, a comparison of the critical residues in RNase A and human angiogenin, which share only 35% amino acid sequence identity, provides a unique perspective on the molecular evolution of the RNase A superfamily, as well as an impetus for applying this methodology to other ribonucleases.

A Magnetic Biomimetic Nanocatalyst for Cleaving Phosphoester and Carboxylic Ester Bonds under Mild Conditions

[Structure: see text] As a result of the unique surface structure of nanospheres, Asp and His residues supported on a 12 nm maghemite nanoparticle worked collaboratively as a biomimetic nanocatalyst for hydrolyzing paraoxon (phosphoester) and 4-nitrophenyl acetate (carboxylic ester) in Milli-Q water (pH 7.0) at 37 degrees C, without employing extremes of pH or heavy metals. Our nanocatalyst could be facilely recovered via magnetic concentration. The isolated catalyst exhibited long-term stability, showing no significant loss of its catalytic activity for repeated uses after 3 months.

Role of aspartic acid 121 in human pancreatic ribonuclease catalysis

Bovine pancreatic ribonuclease (RNase A) is one of the most well studied enzymes of the ribonuclease family, unlike its human counterpart, the human pancreatic ribonuclease (HPR), whose physiological role in the body is not clearly understood. Human pancreatic ribonuclease consists of 128 amino acids and the main residues located in the active site of RNase A are also conserved in HPR. In the current study, to investigate the role of Asp-121 in the catalytic activity of human pancreatic ribonuclease, several variants were generated in which Asp-121 was either mutated to an alanine or C-terminal residues beyond Asp-121, and Phe-120 were deleted. The HPR mutants were cloned, expressed in E. coli and purified to homogeneity, and functionally characterized. The mutation D121A in HPR significantly decreased the rate of the enzymatic reaction, however this decrease was not universally observed for all substrates studied. Removal of the seven C-terminal amino acid residues thereby exposing Asp-121 yielded an HPR mutant with enhanced activity, however a further deletion removing Asp-121 resulted in the complete inactivation of HPR. Our results indicate that Asp-121 is crucial for the catalytic activity of HPR and may be involved in the depolymerization activity of the enzyme.

Investigation of protein refolding using a fractional factorial screen: a study of reagent effects and interactions

A recurring obstacle for structural genomics is the expression of insoluble, aggregated proteins. In these cases, the use of alternative salvage strategies, like in vitro refolding, is hindered by the lack of a universal refolding method. To overcome this obstacle, fractional factorial screens have been introduced as a systematic and rapid method to identify refolding conditions. However, methodical analyses of the effectiveness of refolding reagents on large sets of proteins remain limited. In this study, we address this void by designing a fractional factorial screen to rapidly explore the effect of 14 different reagents on the refolding of 33 structurally and functionally diverse proteins. The refolding data was analyzed using statistical methods to determine the effect of each refolding additive. The screen has been miniaturized for automation resulting in reduced protein requirements and increased throughput. Our results show that the choice of pH and reducing agent had the largest impact on protein refolding. Bis-mercaptoacetamide cyclohexane (BMC) and tris (2-carboxyethylphosphine) (TCEP) were superior reductants when compared to others in the screen. BMC was particularly effective in refolding disulfide-containing proteins, while TCEP was better for nondisulfide-containing proteins. From the screen, we successfully identified a positive synergistic interaction between nondetergent sulfobetaine 201 (NDSB 201) and BMC on Cdc25A refolding. The soluble protein resulting from this interaction crystallized and yielded a 2.2 Angstroms structure. Our method, which combines a fractional factorial screen with statistical analysis of the data, provides a powerful approach for the identification of optimal refolding reagents in a general refolding screen.

Structure of 6-Oxo Camphor Hydrolase H122A Mutant Bound to Its Natural Product, (2S,4S)-α-Campholinic Acid

The crotonase homolog, 6-oxo camphor hydrolase (OCH), catalyzes the desymmetrization of bicyclic beta-diketones to optically active keto acids via an enzymatic retro-Claisen reaction, resulting in the cleavage of a carbon-carbon bond. We have previously reported the structure of OCH (Whittingham, J. L., Turkenburg, J. P., Verma, C. S., Walsh, M. A., and Grogan, G. (2003) J. Biol. Chem. 278, 1744-1750), which suggested the involvement of five residues, His-45, His-122, His-145, Asp-154, and Glu-244, in catalysis. Here we report mutation studies on OCH that reveal that H145A and D154N mutants of OCH have greatly reduced values of k(cat)/K(m) derived from a very large increase in K(m) for the native substrate, 6-oxo camphor. In addition, H122A has a greatly reduced value of k(cat), and its K(m) is five times that of the wild-type. The location of the active site is confirmed by the 1.9-A structure of the H122A mutant of OCH complexed with the minor diastereoisomer of (2S,4S)-alpha-campholinic acid, the natural product of the enzyme. This shows the pendant acetate of the product hydrogen bonded to a His-145/Asp-154 dyad and the endocyclic carbonyl of the cyclopentane ring hydrogen bonded to Trp-40. The results are suggestive of a base-catalyzed mechanism of C-C bond cleavage and provide clues to the origin of prochiral selectivity by the enzyme and to the recruitment of the crotonase fold for alternate modes of transition state stabilization to those described for other crotonase superfamily members.

Identification of the catalytic motif of the microbial ribosome inactivating cytotoxin colicin E3

Colicin E3 is a cytotoxic ribonuclease that specifically cleaves 16S rRNA at the ribosomal A-site to abolish protein synthesis in sensitive Escherichia coli cells. We have performed extensive mutagenesis of the 96-residue colicin E3 cytotoxic domain (E3 rRNase), assayed mutant colicins for in vivo cytotoxicity, and tested the corresponding E3 rRNase domains for their ability to inactivate ribosome function in vitro. From 21 alanine mutants, we identified five positions where mutation resulted in a colicin with no measurable cytotoxicity (Y52, D55, H58, E62, and Y64) and four positions (R40, R42, E60, and R90) where mutation caused a significant reduction in cytotoxicity. Mutations that were found to have large in vivo and in vitro effects were tested for structural integrity through circular dichroism and fluorescence spectroscopy using purified rRNase domains. Our data indicate that H58 and E62 likely act as the acid-base pair during catalysis with other residues likely involved in transition state stabilization. Both the Y52 and Y64 mutants were found to be highly destabilized and this is the likely origin of the loss of their cytotoxicity. The identification of important active site residues and sequence alignments of known rRNase homologs has allowed us to identify other proteins containing the putative rRNase active site motif. Proteins that contained this active site motif included three hemagglutinin-type adhesins and we speculate that these have evolved to deliver a cytotoxic rRNase into eukaryotic cells during pathogenesis.

Nickel and extracellular acidification inhibit the water permeability of human aquaporin-3 in lung epithelial cells

Nickel is a common cause of pneumoconiosis. Here, we show that nickel inactivates aquaporin (AQP)-3, the water channel expressed apically in epithelial cells of human terminal airways. Human AQP3 was transiently transfected into human lung cells, and water permeability was measured in transfected and neighboring untransfected cells. Incubation with NiCl2 rapidly, dose-dependently, and reversibly decreased water permeability in AQP3-expressing cells. Acidification of the extracellular medium also caused rapid, dose-dependent, and reversible inhibition of AQP3. Sensitivity of AQP3 to nickel was lower at alkaline pH than at neutral and acidic pH. Cells transfected with human AQP4 and AQP5, which are also expressed in airway epithelia, were insensitive to nickel and extracellular acidification. Zinc and cadmium, other common causes of pneumoconiosis, had no effect on the water permeability of AQP3. Three extracellular residues, Trp128, Ser152, and His241, were responsible for the blocking effect of nickel on human AQP3. Ser152 was identified as a common site for nickel and pH sensitivity. His53, Tyr124, and His154 were also involved in regulation of AQP3 by extracellular pH. In addition, the aromatic side chain of His154 was shown to be important for the water permeability of AQP3. Our results imply that nickel and extracellular pH may modulate lung water clearance and that defective water clearance may be an early component of nickel-induced lung disease.

Global and local motions in ribonuclease A: a molecular dynamics study

The understanding of protein dynamics is one of the major goals of structural biology. A direct link between protein dynamics and function has been provided by x-ray studies performed on ribonuclease A (RNase A) (B. F. Rasmussen et al., Nature, 1992, Vol. 357, pp. 423-424; L. Vitagliano et al., Proteins: Structure, Function, and Genetics, 2002, Vol. 46, pp. 97-104). Here we report a 3 ns molecular dynamics simulation of RNase A in water aimed at characterizing the dynamical behavior of the enzyme. The analysis of local and global motions provides interesting insight on the dynamics/function relationship of RNase A. In agreement with previous crystallographic reports, the present study confirms that the RNase A active site is constituted by rigid (His12, Asn44, Thr45) and flexible (Lys41, Asp83, His119, Asp121) residues. The analysis of the global motions, performed using essential dynamics, shows that the two beta-sheet regions of RNase A move coherently in opposite directions, thus modifying solvent accessibility of the active site, and that the mixed alpha/3(10)-helix (residues 50-60) behaves as a mechanical hinge during the breathing motion of the protein. These data demonstrate that this motion, essential for RNase A substrate binding and release, is an intrinsic dynamical property of the ligand-free enzyme.

Conformational strictness required for maximum activity and stability of bovine pancreatic ribonuclease A as revealed by crystallographic study of three Phe120 mutants at 1.4 A resolution

The replacement of Phe120 with other hydrophobic residues causes a decrease in the activity and thermal stability in ribonuclease A (RNase A). To explain this, the crystal structures of wild-type RNase A and three mutants--F120A, F120G, and F120W--were analyzed up to a 1.4 A resolution. Although the overall backbone structures of all mutant samples were nearly the same as that of wild-type RNase A, except for the C-terminal region of F120G with a high B-factor, two local conformational changes were observed at His119 in the mutants. First, His119 of the wild-type and F120W RNase A adopted an A position, whereas those of F120A and F120G adopted a B position, but the static crystallographic position did not reflect either the efficiency of transphosphorylation or the hydrolysis reaction. Second, His119 imidazole rings of all mutant enzymes were deviated from that of wild-type RNase A, and those of F120W and F120G appeared to be "inside out" compared with that of wild-type RNase A. Only approximately 1 A change in the distance between N(epsilon2) of His12 and N(delta1) of His119 causes a drastic decrease in k(cat), indicating that the active site requires the strict positioning of the catalytic residues. A good correlation between the change in total accessible surface area of the pockets on the surface of the mutant enzymes and enthalpy change in their thermal denaturation also indicates that the effects caused by the replacements are not localized but extend to remote regions of the protein molecule.

Glutamate 170 of human l-3-hydroxyacyl-CoA dehydrogenase is required for proper orientation of the catalytic histidine and structural integrity of the enzyme

l-3-Hydroxyacyl-CoA dehydrogenase (HAD), the penultimate enzyme in the beta-oxidation spiral, reversibly catalyzes the conversion of l-3-hydroxyacyl-CoA to the corresponding 3-ketoacyl-CoA. Similar to other dehydrogenases, HAD contains a general acid/base, His(158), which is within hydrogen bond distance of a carboxylate, Glu(170). To investigate its function in this catalytic dyad, Glu(170) was replaced with glutamine (E170Q), and the mutant enzyme was characterized. Whereas substrate and cofactor binding were unaffected by the mutation, E170Q exhibited diminished catalytic activity. Protonation of the catalytic histidine did not restore wild-type activity, indicating that modulation of the pK(a) of His(158) is not the sole function of Glu(170). The pH profile of charge transfer complex formation, an independent indicator of active site integrity, was unaltered by the amino acid substitution, but the intensity of the charge transfer band was diminished. This observation, coupled with significantly reduced enzymatic stability of the E170Q mutant, implicates Glu(170) in maintenance of active site architecture. Examination of the crystal structure of E170Q in complex with NAD(+) and acetoacetyl-CoA (R = 21.9%, R(free) = 27.6%, 2.2 A) reveals that Gln(170) no longer hydrogen bonds to the side chain of His(158). Instead, the imidazole ring is nearly perpendicular to its placement in the comparable native complex and no longer positioned for efficient catalysis.

On the role of conserved histidine 106 in 10-formyltetrahydrofolate dehydrogenase catalysis: connection between hydrolase and dehydrogenase mechanisms

The enzyme, 10-formyltetrahydrofolate dehydrogenase (FDH), converts 10-formyltetrahydrofolate (10-formyl-THF) to tetrahydrofolate in an NADP(+)-dependent dehydrogenase reaction or an NADP(+)-independent hydrolase reaction. The hydrolase reaction occurs in a 310-amino acid long amino-terminal domain of FDH (N(t)-FDH), whereas the dehydrogenase reaction requires the full-length enzyme. The amino-terminal domain of FDH shares some sequence identity with several other enzymes utilizing 10-formyl-THF as a substrate. These enzymes have two strictly conserved residues, aspartate and histidine, in the putative catalytic center. We have shown recently that the conserved aspartate is involved in FDH catalysis. In the present work we studied the role of the conserved histidine, His(106), in FDH function. Site-directed mutagenesis experiments showed that replacement of the histidine with alanine, asparagine, aspartate, glutamate, glutamine, or arginine in N(t)-FDH resulted in expression of insoluble proteins. Replacement of the histidine with another positively charged residue, lysine, produced a soluble mutant with no hydrolase activity. The insoluble mutants refolded from inclusion bodies adopted a conformation inherent to the wild-type N(t)-FDH, but they did not exhibit any hydrolase activity. Substitution of alanine for three non-conserved histidines located close to the conserved one did not reveal any significant changes in the hydrolase activity of N(t)-FDH. Expressed full-length FDH with the substitution of lysine for the His(106) completely lost both the hydrolase and dehydrogenase activities. Thus, our study showed that His(106), besides being an important structural residue, is also directly involved in both the hydrolase and dehydrogenase mechanisms of FDH. Modeling of the putative hydrolase catalytic center/folate-binding site suggested that the catalytic residues, aspartate and histidine, are unlikely to be adjacent to the catalytic cysteine in the aldehyde dehydrogenase catalytic center. We hypothesize that 10-formyl-THF dehydrogenase reaction is not an independent reaction but is a combination of hydrolase and aldehyde dehydrogenase reactions.

Pentavalent Organo-Vanadates as Transition State Analogues for Phosphoryl Transfer Reactions

Pentavalent organo-vanadates have been put forth as transition state analogues for a variety of phosphoryl transfer reactions. In particular, uridine 2',3'-cyclic vanadate (U>v) has been proposed to resemble the transition state during catalysis by ribonuclease A (RNase A). Here, this hypothesis is tested. Lys41 of RNase A is known to donate a hydrogen bond to a nonbridging phosphoryl oxygen in the transition state during catalysis. Site-directed mutagenesis and semisynthesis were used to create enzymes with natural and nonnatural amino acid residues at position 41. These variants differ by 10(5)-fold in their k(cat)/K(m) values for catalysis, but <40-fold in their K(i) values for inhibition of catalysis by U>v. Plots of logK(i) vs log(K(m)/k(cat)) for three distinct substrates [poly(cytidylic acid), uridine 3'-(p-nitrophenyl phosphate), and cytidine 2',3'-cyclic phosphate] have slopes that range from 0.25 and 0.36. These plots would have a slope of unity if U>v were a perfect transition state analogue. Values of K(i) for U>v correlate weakly with the equilibrium dissociation constant for the enzymic complexes with substrate or product, indicating that U>v bears some resemblance to the substrate and product as well as the transition state. Thus, U>v is a transition state analogue for RNase A, but only a marginal one. This finding indicates that a pentavalent organo-vanadate cannot necessarily be the basis for a rigorous analysis of the transition state for a phosphoryl transfer reaction.

Productive and nonproductive binding to ribonuclease A: X-ray structure of two complexes with uridylyl(2',5')guanosine

Guanine-containing mono- and dinucleotides bind to the active site of ribonuclease A in a nonproductive mode (retro-binding) (Aguilar CF, Thomas PJ, Mills A, Moss DS, Palmer RA. 1992. J Mol Biol 224:265-267). Guanine binds to the highly specific pyrimidine site by forming hydrogen bonds with Thr45 and with the sulfate anion located in the P1 site. To investigate the influence of the anion present in the P1 site on retro-binding, we determined the structure of two new complexes of RNase A with uridylyl(2',5')guanosine obtained by soaking two different forms of pre-grown RNase A crystals. In one case, RNase A was crystallized without removing the sulfate anion strongly bound to the active site; in the other, the protein was first equilibrated with a basic solution to displace the anion from the P1 site. The X-ray structures of the complexes with and without sulfate in P1 were refined using diffraction data up to 1.8 A (R-factor 0.192) and 2.0 A (R-factor 0.178), respectively. The binding mode of the substrate analogue to the protein differs markedly in the two complexes. When the sulfate is located in P1, we observe retro-binding; whereas when the anion is removed from the active site, the uridine is productively bound at the B1 site. In the productive complex, the electron density is very well defined for the uridine moiety, whereas the downstream guanine is disordered. This finding indicates that the interactions of guanine in the B2 site are rather weak and that this site is essentially adenine preferring. In this crystal form, there are two molecules per asymmetric unit, and due to crystal packing, only the active site of one molecule is accessible to the ligand. Thus, in the same crystal we have a ligand-bound and a ligand-free RNase A molecule. The comparison of these two structures furnishes a detailed and reliable picture of the structural alterations induced by the binding of the substrate. These results provide structural information to support the hypotheses on the role of RNase A active site residues that have recently emerged from site-directed mutagenesis studies.

The ultrahigh resolution crystal structure of ribonuclease A containing an isoaspartyl residue: hydration and sterochemical analysis

Crystals of the deamidated form of bovine pancreatic ribonuclease which contains an isoaspartyl residue in position 67 diffract to 0. 87 A at 100 K. We have refined the crystallographic model using anisotropic displacement parameters for all atoms to a conventional crystallographic residual R=0.101 for all observed reflections in the resolution range 61.0-0.87 A. The ratio observations/parameters is 7.2 for the final model. This structure represents one of the highest resolution protein structures to date and interestingly, it is the only example containing more than one molecule in the asymmetric unit with a resolution better than 1.0 A. The non-crystallographic symmetry has been used as a validation check of the geometrical parameters and it has allowed an estimate for an upper limit of errors associated with this high resolution model. In the present structure it was possible to obtain a more accurate picture of the active site whose electron density was not clearly interpretable in the previous 1.9 A resolution structure. In particular, the P1 site is alternatively occupied either by a sulphate anion or by a water molecule network. Most of hydrogen atoms were visible in the electron density maps, including those involved in C(alpha)-H(alpha).O interactions. Analysis of protein-solvent interactions has revealed the occurrence of an extensive cluster of water molecules, predominantly arranged in pentagonal fused rings and surrounding hydrophobic moiety of side-chains. Finally, in spite of the limited sample of residues, we have detected a clear dependence of backbone N-C(alpha)-C angle on residue conformation. This correlation can be fruitfully used as a valuable tool in protein structure validation.