The structures of RNase A complexed with 3'-CMP and d(CpA): active site conformation and conserved water molecules

The three-dimensional structure of human RNase 4, unliganded and complexed with d(Up), reveals the basis for its uridine selectivity

The RNase 4 family is unique among RNase enzymes, displaying the highest level of sequence similarity and encompassing the shortest polypeptide chain. It is the only one showing high specificity. The human representative is an intracellular and plasma enzyme, first isolated from colon adenocarcinoma cell line HT-29. The crystal structures of human recombinant RNase 4, unliganded and in complex with d(Up), have been determined, revealing in the unique active site an explanation for the uridine specificity. Arg101, at a position not involved in catalysis in the other RNase enzymes, penetrates the enzyme moiety shaping the recognition pocket, a flip that is mediated by the interaction with the (shorter chain) C-terminal carboxylate group, providing an anchoring point for the O4 atom of the substrate uridine. The bulky Phe42 side-chain forces Asp80 to be in the chi1=-72.49 degrees rotamer, accepting a hydrogen bond from Thr44, further converting the latter into a hydrogen bond acceptor. This favours an interaction with the -NH-donor group of uridine at position 3 over that with the =N-acceptor of cytidine. The two chemical groups that distinguish uracyl from cytosine are used by the enzyme to discriminate between these two bases.

A new remote subsite in ribonuclease A

The interaction between bovine pancreatic ribonuclease A (RNase A) and its RNA substrate extends beyond the scissile bond. Enzymic subsites interact with the bases and the phosphoryl groups of a bound substrate. We evaluated the four cationic residues closest to known subsites for their abilities to interact with a bound nucleic acid. Lys-37, Arg-39, Arg-85, and Lys-104 were replaced individually by an alanine residue, and the resulting enzymes were assayed as catalysts of poly(cytidylic acid) (poly(C)) cleavage. The values of Km and kcat/Km for poly(C) cleavage were affected only by replacing Arg-85. Moreover, the contribution of Arg-85 to the binding of the ground state and the transition state was uniform---Km increased by 15-fold and kcat/Km decreased by 10-fold. The contribution of Arg-85 to binding was also apparent in the values of Kd for complexes with oligonucleotides of different length. This contribution was dependent on salt concentration, as expected from a coulombic interaction between a cationic side chain and an anionic phosphoryl group. Together, these data indicate that Arg-85 interacts with a particular phosphoryl group of a bound nucleic acid. We propose that Arg-85 comprises a new distal subsite in RNase A---the P(-1) subsite.

Modeling of angiogenin - 3-NMP complex

Angiogenin belongs to the Ribonuclease superfamily and has a weak enzymatic activity that is crucial for its biological function of stimulating blood vessel growth. Structural studies on ligand bound Angiogenin will go a long way in understanding the mechanism of the protein as well as help in designing drugs against it. In this study we present the first available structure of nucleotide ligand bound Angiogenin obtained by computer modeling. The importance of this study in itself notwithstanding, is a precursor to modeling a full dinucleotide substrate onto Angiogenin. Bovine Angiogenin, the structure of which has been solved at a high resolution, was earlier subjected to Molecular Dynamics simulations for a nanosecond. The MD structures offer better starting points for docking as they offer lesser obstruction than the crystal structure to ligand binding. The MD structure with the least serious short contacts was modeled to obtain a steric free Angiogenin - 3' mononucleotide complex structure. The structures were energetically minimized and subjected to a brief spell of Molecular Dynamics. The results of the simulation show that all the ligand-Angiogenin interactions and hydrogen bonds are retained, redeeming the structure and docking procedure. Further, following ligand - protein interactions in the case of the ligands 3'-CMP and 3'-UMP we were able to speculate on how Angiogenin, a predominantly prymidine specific ribonuclease prefers Cytosine to Uracil in the first base position.

The solution structure of a cytotoxic ribonuclease from the oocytes of Rana catesbeiana (bullfrog)

RC-RNase is a pyrimidine-guanine sequence-specific ribonuclease and a lectin possessing potent cell cytotoxicity. It was isolated from the oocytes of Rana catesbeiana (bull frog). From analysis of an extensive set of 1H homonuclear 2D NMR spectra we have completed the resonance assignments. Determination of the three-dimensional structure was carried out with the program X-PLOR using a total of 951 restraints including 814 NMR-derived distances, 61 torsion angles, and 76 hydrogen bond restraints. In the resultant family of 15 best structures, selected from a total of 150 calculated structures, the root-mean-square deviation from the average structure for the backbone heavy-atoms involved in well-defined secondary structure is 0.48 A, while that for all backbone heavy-atoms is 0.91 A. The structure of RC-RNase consists of three alpha-helices and two triple-stranded anti-parallel beta-sheets and folds in a kidney-shape, very similar to the X-ray crystal structure of a homolo gous protein, onconase isolated from Rana pipiens. We have also investigated the interaction between RC-RNase and two inhibitors, cytidylyl(2'-->5')guanosine (2',5'-CpG) and 2'-deoxycytidylyl(3'-->5')-2'-deoxyguanosine (3',5'-dCpdG). Based on the ligand-induced chemical shift changes in RC-RNase and the NOE cross-peaks between RC-RNase and the inhibitors, the key residues involved in protein-inhibitor interaction have been identified. The inhibitors were found to bind in a "retro-binding" mode, with the guanine base bonded to the B1 subsite. The His103 residue was found to occupy the B state with the imidazole ring pointing away from the active site. The structure coordinates and the NMR restraints have been deposited in the Brookhaven Protein Data Bank (1bc4 and 1bc4mr, respectively).

Binding of a substrate analog to a domain swapping protein: X-ray structure of the complex of bovine seminal ribonuclease with uridylyl(2',5')adenosine

Bovine seminal ribonuclease (BS-RNase) is a unique member of the pancreatic-like ribonuclease superfamily. The native enzyme is a mixture of two dimeric forms with distinct structural features. The most abundant form is characterized by the swapping of N-terminal fragments. In this paper, the crystal structure of the complex between the swapping dimer and uridylyl(2',5')adenosine is reported at 2.06 A resolution. The refined model has a crystallographic R-factor of 0.184 and good stereochemistry. The quality of the electron density maps enables the structure of both the inhibitor and active site residues to be unambiguously determined. The overall architecture of the active site is similar to that of RNase A. The dinucleotide adopts an extended conformation with the pyrimidine and purine base interacting with Thr45 and Asn71, respectively. Several residues (Gln11, His12, Lys41, His119, and Phe120) bind the oxygens of the phosphate group. The structural similarity of the active sites of BS-RNase and RNase A includes some specific water molecules believed to be relevant to catalytic activity. Upon binding of the dinucleotide, small but significant modifications of the tertiary and quaternary structure of the protein are observed. The ensuing correlation of these modifications with the catalytic activity of the enzyme is discussed.

Structure and stability of the P93G variant of ribonuclease A

The peptide bonds preceding Pro 93 and Pro 114 of bovine pancreatic ribonuclease A (RNase A) are in the cis conformation. The trans-to-cis isomerization of these bonds had been indicted as the slow step during protein folding. Here, site-directed mutagenesis was used to replace Pro 93 or Pro 114 with a glycine residue, and the crystalline structure of the P93G variant was determined by X-ray diffraction analysis to a resolution of 1.7 A. This structure is essentially identical to that of the wild-type protein, except for the 91-94 beta-turn containing the substitution. In the wild-type protein, the beta-turn is of type VIa. In the P93G variant, this turn is of type II with the peptide bond preceding Gly 93 being trans. The thermal stabilities of the P93G and P114G variants were assessed by differential scanning calorimetry and thermal denaturation experiments monitored by ultraviolet spectroscopy. The value of delta deltaGm which reports on the stability lost in the variants, is 1.5-fold greater for the P114G variant than for the P93G variant. The greater stability of the P93G variant is likely due to the relatively facile accommodation of residues 91-94 in a type II turn, which has a preference for a glycine residue in its i + 2 position.

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

The side chains of histidine and aspartate residues form a hydrogen bond in the active sites of many enzymes. In serine proteases, the His...Asp hydrogen bond of the catalytic triad is known to contribute greatly to catalysis, perhaps via the formation of a low-barrier hydrogen bond. In bovine pancreatic ribonuclease A (RNase A), the His...Asp dyad is composed of His119 and Asp121. Previously, site-directed mutagenesis was used to show that His119 has a fundamental role, to act as an acid during catalysis of RNA cleavage [Thompson, J. E., and Raines, R. T. (1994) J. Am. Chem. Soc. 116, 5467-5468]. Here, Asp121 was replaced with an asparagine or alanine residue. The crystalline structures of the two variants were determined by X-ray diffraction analysis to a resolution of 1.6 A with an R-factor of 0.18. Replacing Asp121 with an asparagine or alanine residue does not perturb the overall conformation of the enzyme. In the structure of D121N RNase A, Ndelta rather than Odelta of Asn121 faces His119. This alignment in the crystalline state is unlikely to exist in solution because catalysis by the D121N variant is not compromised severely. The steady-state kinetic parameters for catalysis by the wild-type and variant enzymes were determined for the cleavage of uridylyl(3'-->5')adenosine and poly(cytidylic acid), and for the hydrolysis of uridine 2',3'-cyclic phosphate. Replacing Asp121 decreases the values of kcat/Km and kcat for cleavage by 10-fold (D121N) and 10(2)-fold (D121A). Replacing Asp121 also decreases the values of kcat/Km and kcat for hydrolysis by 10(0. 5)-fold (D121N) and 10-fold (D121A) but has no other effect on the pH-rate profiles for hydrolysis. There is no evidence for the formation of a low-barrier hydrogen bond between His119 and either an aspartate or an asparagine residue at position 121. Apparently, the major role of Asp121 is to orient the proper tautomer of His119 for catalysis. Thus, the mere presence of a His...Asp dyad in an enzymic active site is not a mandate for its being crucial in effecting catalysis.

Conserved water molecules contribute to the extensive network of interactions at the active site of protein kinase A

Protein kinases constitute a large family of regulatory enzymes, each with a distinct specificity to restrict its action to its physiological target(s) only. The catalytic (C) subunit of protein kinase A, regarded as a structural prototype for this family, is composed of a conserved core flanked by two nonconserved segments at the amino and carboxyl termini. Here we summarize evidence to show that (i) the active site consists of an extended network of interactions that weave together both domains of the core as well as both segments that flank the core; (ii) the opening and closing of the active site cleft, including the dynamic and coordinated movement of the carboxyl terminal tail, contributes directly to substrate recognition and catalysis; and (iii) in addition to peptide and ATP, the active site contains six structured water molecules that constitute a conserved structural element of the active site. One of these active-site conserved water molecules is locked into place by its interactions with the nucleotide, the peptide substrate/inhibitor, the small and large domains of the conserved core, and Tyr-330 from the carboxyl-terminal "tail."

Three-dimensional solution structure of human angiogenin determined by 1H,15N-NMR spectroscopy--characterization of histidine protonation states and pKa values

Human angiogenin is a member of the pancreatic ribonuclease superfamily that induces blood vessel formation. Its three-dimensional solution structure has been determined to high resolution by heteronuclear NMR spectroscopy. 30 structures were calculated, based on a total of 1441 assigned NOE correlations, 64 coupling constants and 50 hydrogen bonds. The backbone atomic rms difference from the mean coordinates is 0.067 +/- 0.012 nm and 0.13 nm from the previously determined crystal structure. The side-chain of Gln117 was found to obstruct the active site as observed in the crystal state. There was no evidence of an alternative open form of angiogenin, although two sets of chemical shifts were observed for some residues, mainly around the active site and in the C-terminal segment. The topology of the ribonucleolytic active site is described with a particular emphasis on the conformation and protonation of active-site His residues. The side-chain of His114 adopts two main conformations in solution. In contrast to pancreatic ribonuclease A, His13 was shown to be more basic than His114, with pKa values of 6.65 and 6.05 respectively. The His47 residue is located in an environment very resistant to protonation with a pKa lower than 4.

Molecular mechanics calculations on Rous sarcoma virus protease with peptide substrates

Molecular models of Rous sarcoma virus (RSV) protease and 20 peptide substrates with single amino acid substitutions at positions from P4 to P3', where the scissile bond is between P1 and P1'. were built and compared with kinetic measurements. The unsubstituted peptide substrate. Pro-Ala-Val-Ser-Leu-Ala-Met-Thr, represents the NC-PR cleavage site of RSV protease. Models were built of two intermediates in the catalytic reaction, RSV protease with peptide substrate and with the tetrahedral intermediate. The energy minimization used an algorithm that increased the speed and eliminated a cutoff for nonbonded interactions. The calculated protease-substrate interaction energies showed correlation with the relative catalytic efficiency of peptide hydrolysis. The calculated interaction energies for the 8 RSV protease-substrate models with changes in P1 to P1' next to the scissile bond gave the highest correlation coefficient of 0.79 with the kinetic measurements, whereas all 20 substrates showed the lower, but still significant correlation of 0.46. Models of the tetrahedral reaction intermediates gave a correlation of 0.72 for the 8 substrates with changes next to the scissile bond, whereas a correlation coefficient of only 0.34 was observed for all 20 substrates. The differences between the energies calculated for the tetrahedral intermediate and the bound peptide gave the most significant correlation coefficients of 0.90 for models with changes in P1 and P1', and 0.56 for all substrates. These results are compared to those from similar calculations on HIV-1 protease and discussed in relation to the rate-limiting steps in the catalytic mechanism and the entropic contributions.

Crystallographic analysis of human immunodeficiency virus 1 protease with an analog of the conserved CA-p2 substrate -- interactions with frequently occurring glutamic acid residue at P2' position of substrates

Human immunodeficiency virus type 1 (HIV-1) protease hydrolysis of the Gag CA-p2 cleavage site is crucial for virion maturation and is optimal at acidic pH. To understand the processing of the CA-p2 site, we have determined the structure of HIV-1 protease complexed with an analog of the CA-p2 site, the reduced peptide inhibitor Arg-Val-Leu-r-Phe-Glu-Ala-Ahx-NH2 [r denotes the reduced peptide bond and Ahx 2-aminohexanoic acid (norleucine), respectively]. The crystal structure was refined to an R-factor of 0.17 at 0.21-nm resolution. The crystals have nearly the same lattice as related complexes in P2(1)2(1)2(1) which have twofold disordered inhibitor, but are in space group P2(1). and the asymmetric unit contains two dimers of HIV-1 protease related by 180 degrees rotation. An approximate non-crystallographic symmetry has replaced the exact crystal symmetry resulting in well-ordered inhibitor structure. Each protease dimer binds one ordered inhibitor molecule, but in opposite orientations. The interactions of the inhibitor with the two dimers are very similar for the central P2 Val to P2' Glu residues, but show more variation for the distal P3 Arg and P4' Ahx residues. Importantly, the carboxylate oxygens of Glu at P2' in the inhibitor are within hydrogen-bonding distance of a carboxylate oxygen of Asp30 of the protease suggesting that the two side chains share a proton. This interaction suggests that the enzyme-substrate complex is additionally stabilized at lower pH. The importance of this interaction is emphasized by the absence of polymorphisms of Asp30 in the protease and variants of P2' Glu in the critical CA-p2 cleavage site.

Crystal structures of the copper and nickel complexes of RNase A: metal-induced interprotein interactions and identification of a novel copper binding motif

We report the crystal structures of the copper and nickel complexes of RNase A. The overall topology of these two complexes is similar to that of other RNase A structures. However, there are significant differences in the mode of binding of copper and nickel. There are two copper ions per molecule of the protein, but there is only one nickel ion per molecule of the protein. Significant changes occur in the interprotein interactions as a result of differences in the coordinating groups at the common binding site around His-105. Consequently, the copper- and nickel-ion-bound dimers of RNase A act as nucleation sites for generating different crystal lattices for the two complexes. A second copper ion is present at an active site residue His-119 for which all the ligands are from one molecule of the protein. At this second site, His-119 adopts an inactive conformation (B) induced by the copper. We have identified a novel copper binding motif involving the alpha-amino group and the N-terminal residues.

A decade of protein engineering on ribonuclease T1--atomic dissection of the enzyme-substrate interactions

During the last decade, protein engineering has been used to identify the residues that contribute to the ribonuclease-T1-catalyzed transesterification. His40, Glu58 and His92 accelerate the associative nucleophilic displacement at the phosphate atom by the entering 2'-oxygen downstream guanosines in a highly cooperative manner. Glu58, assisted by the protonated His40 imidazole, abstracts a proton from the 2'-oxygen, while His92 protonates the leaving group. Tyr38, Arg77 and Phe100 further stabilize the transition state of the reaction. A functionally independent subsite, including Asn36 and Asn98, contributes to chemical turnover by aligning the substrate relative to the catalytic side chains upon binding of the leaving group. An invariant structural motive, involving residues 42-46, renders ribonuclease T1 guanine specific through a series of intermolar hydrogen bonds. Tyr42 contributes significantly to guanine binding through a parallel face-to-face stacking interaction. Tyr45, often referred to as the lid of the guanine-binding site, does not contribute to the binding of the base.

5'-Diphosphoadenosine 3'-phosphate is a potent inhibitor of bovine pancreatic ribonuclease A

As a first step toward the development of stable, selective, and potent inhibitors of those members of the pancreatic RNase superfamily that induce biological responses, we have focussed on low molecular weight compounds and studied their interactions with the active-site of bovine pancreatic ribonuclease A (RNase A). A new inhibitor is described, 5'-diphosphoadenosine 3'-phosphate, which binds to RNase A more tightly than any previous low molecular weight compound: its Ki value of 1.3 microM at pH 7 is 8-fold lower than that for uridine-vanadate, a transition-state analog, and 110-fold lower than that for 2'-CMP, one of the best-characterized RNase A ligands. The new inhibitor is found to contact RNase A residues that are conserved in several homologous mammalian RNases and hence should be able to serve as a basis for the design of even tighter-binding inhibitors of these enzymes.

All in one: a highly detailed rotamer library improves both accuracy and speed in the modelling of sidechains by dead-end elimination

BACKGROUND

About a decade ago, the concept of rotamer libraries was introduced to model sidechains given known mainchain coordinates. Since then, several groups have developed methods to handle the challenging combinatorial problem that is faced when searching rotamer libraries. To avoid a combinatorial explosion, the dead-end elimination method detects and eliminates rotamers that cannot be members of the global minimum energy conformation (GMEC). Several groups have applied and further developed this method in the fields of homology modelling and protein design.

RESULTS

This work addresses at the same time increased prediction accuracy and calculation speed improvements. The proposed enhancements allow the elimination of more than one-third of the possible rotameric states before applying the dead-end elimination method. This is achieved by using a highly detailed rotamer library allowing the safe application of an energy-based rejection criterion without risking the elimination of a GMEC rotamer. As a result, we gain both in modelling accuracy and in computational speed. Being completely automated, the current implementation of the dead-end elimination prediction of protein sidechains can be applied to the modelling of sidechains of proteins of any size on the high-end computer systems currently used in molecular modelling. The improved accuracy is highlighted in a comparative study on a collection of proteins of varying size for which score results have previously been published by multiple groups. Furthermore, we propose a new validation method for the scoring of the modelled structure versus the experimental data based upon the volume overlap of the predicted and observed sidechains. This overlap criterion is discussed in relation to the classic RMSD and the frequently used +/- 40 degrees window in comparing chi 1 and chi 2 angles.

CONCLUSIONS

We have shown that a very detailed library allows the introduction of a safe energy threshold rejection criterion, thereby increasing both the execution speed and the accuracy of the modelling program. We speculate that the current method will allow the sidechain prediction of medium-sized proteins and complex protein interfaces involving up to 150 residues on low-end desktop computers.

Mechanism of ribonuclease inhibition by ribonuclease inhibitor protein based on the crystal structure of its complex with ribonuclease A

We describe the mechanism of ribonuclease inhibition by ribonuclease inhibitor, a protein built of leucine-rich repeats, based on the crystal structure of the complex between the inhibitor and ribonuclease A. The structure was determined by molecular replacement and refined to an Rcryst of 19.4% at 2.5 A resolution. Ribonuclease A binds to the concave region of the inhibitor protein comprising its parallel beta-sheet and loops. The inhibitor covers the ribonuclease active site and directly contacts several active-site residues. The inhibitor only partially mimics the RNase-nucleotide interaction and does not utilize the p1 phosphate-binding pocket of ribonuclease A, where a sulfate ion remains bound. The 2550 A2 of accessible surface area buried upon complex formation may be one of the major contributors to the extremely tight association (Ki = 5.9 x 10(-14) M). The interaction is predominantly electrostatic; there is a high chemical complementarity with 18 putative hydrogen bonds and salt links, but the shape complementarity is lower than in most other protein-protein complexes. Ribonuclease inhibitor changes its conformation upon complex formation; the conformational change is unusual in that it is a plastic reorganization of the entire structure without any obvious hinge and reflects the conformational flexibility of the structure of the inhibitor. There is a good agreement between the crystal structure and other biochemical studies of the interaction. The structure suggests that the conformational flexibility of RI and an unusually large contact area that compensates for a lower degree of complementarity may be the principal reasons for the ability of RI to potently inhibit diverse ribonucleases. However, the inhibition is lost with amphibian ribonucleases that have substituted most residues corresponding to inhibitor-binding residues in RNase A, and with bovine seminal ribonuclease that prevents inhibitor binding by forming a dimer.

Ribonuclease A mutant His119 Asn: the role of histidine in catalysis

Bovine pancreatic ribonuclease A (RNase A) has been widely used as a convenient model for structural and functional studies. The enzyme catalyzes cleavage of phosphodiester bonds in RNA and related substrates. Three amino acid residues located at the active site of RNase A (His12, His119, and Lys41) are known to be involved in catalysis. Mutation of His119 to asparagine was generated to study the role of His119 in RNase A catalysis. The mutant enzyme has been isolated and characterized. The mutation significantly decreases the rate of the transesterification reaction and has no effect on substrate affinity of the enzyme. An analysis of the enzymatic properties of H119N RNase A suggests that the imidazole ring of His119 of the wild-type enzyme must be protonated in an enzyme-substrate productive complex. Thus our results indicate that a contribution of protonated His119 into the catalysis is not restricted to protonation of oxygen atom of the substrate leaving group and that His119 participates directly in a transition state stabilization via hydrogen bonding.

Packing at the protein-water interface

We have determined the packing efficiency at the protein-water interface by calculating the volumes of atoms on the protein surface and nearby water molecules in 22 crystal structures. We find that an atom on the protein surface occupies, on average, a volume approximately 7% larger than an atom of equivalent chemical type in the protein core. In these calculations, larger volumes result from voids between atoms and thus imply a looser or less efficient packing. We further find that the volumes of individual atoms are not related to their chemical type but rather to their structural location. More exposed atoms have larger volumes. Moreover, the packing around atoms in locally concave, grooved regions of protein surfaces is looser than that around atoms in locally convex, ridge regions. This as a direct manifestation of surface curvature-dependent hydration. The net volume increase for atoms on the protein surface is compensated by volume decreases in water molecules near the surface. These waters occupy volumes smaller than those in the bulk solvent by up to 20%; the precise amount of this decrease is directly related to the extent of contact with the protein.

A catalytic function for the structurally conserved residue Phe 100 of ribonuclease T1

The function of the conserved Phe 100 residue of RNase T1 (EC 3.1.27.3) has been investigated by site-directed mutagenesis and X-ray crystallography. Replacement of Phe 100 by alanine results in a mutant enzyme with kcat reduced 75-fold and a small increase in Km for the dinucleoside phosphate substrate GpC. The Phe 100 Ala substitution has similar effects on the turnover rates of GpC and its minimal analogue GpOMe, in which the leaving cytidine is replaced by methanol. The contribution to catalysis is independent of the nature of the leaving group, indicating that Phe 100 belongs to the primary site. The contribution of Phe 100 to catalysis may result from a direct van der Waals contact between its aromatic ring and the phosphate moiety of the substrate. Phe 100 may also contribute to the positioning of the pentacovalent phosphorus of the transition state, relative to other catalytic residues. If compared to the corresponding wild-type data, the structural implications of the mutation in the present crystal structure of Phe 100 Ala RNase T1 complexed with the specific inhibitor 2'-GMP are restricted to the active site. Repositioning of 2'-GMP, caused by the Phe 100 Ala mutation, generates new or improved contacts of the phosphate moiety with Arg 77 and His 92. In contrast, interactions with the Glu 58 carboxylate appear to be weakened. The effects of the His 92 Gln and Phe 100 Ala mutations on GpC turnover are additive in the corresponding double mutant, indicating that the contribution of Phe 100 to catalysis is independent of the catalytic acid His 92. The present results lead to the conclusion that apolar residues may contribute considerably to catalyze conversions of charged molecules to charged products, involving even more polar transition states.

Three-dimensional structure of the complexes of ribonuclease A with 2',5'-CpA and 3',5'-d(CpA) in aqueous solution, as obtained by NMR and restrained molecular dynamics

The three-dimensional structure of the complexes of ribonuclease A with cytidyl-2',5'-adenosine (2',5'-CpA) and deoxycytidyl-3',5'-deoxyadenosine [3',5'-d(CpA)] in aqueous solution has been determined by 1H NMR methods in combination with restrained molecular dynamics calculations. Twenty-three intermolecular NOE cross-corrections for the 3',5'-d(CpA) complex and 19 for the 2',5'-CpA, together with about 1,000 intramolecular NOEs assigned for each complex, were translated into distance constraints and used in the calculation. No significant changes in the global structure of the enzyme occur upon complex formation. The side chains of His 12, Thr 45, His 119, and the amide backbone group of Phe 120 are involved directly in the binding of the ligands at the active site. The conformation of the two bases is anti in the two complexes, but differs from the crystal structure in the conformation of the two sugar rings in 3',5'-d(CpA), shown to be in the S-type region, as deduced from an analysis of couplings between the ribose protons. His 119 is found in the two complexes in only one conformation, corresponding to position A in the free protein. Side chains of Asn 67, Gln 69, Asn 71, and Glu 111 from transient hydrogen bonds with the adenine base, showing the existence of a pronounced flexibility of these enzyme side chains at the binding site of the downstream adenine. All other general features on the structures coincide clearly with those observed in the crystal state.

Analysis of the mechanism of the Serratia nuclease using site-directed mutagenesis

Based on crystal structure analysis of the Serratia nuclease and a sequence alignment of six related nucleases, conserved amino acid residues that are located in proximity to the previously identified catalytic site residue His89 were selected for a mutagenesis study. Five out of 12 amino acid residues analyzed turned out to be of particular importance for the catalytic activity of the enzyme: Arg57, Arg87, His89, Asn119 and Glu127. Their replacement by alanine, for example, resulted in mutant proteins of very low activity, < 1% of the activity of the wild-type enzyme. Steady-state kinetic analysis of the mutant proteins demonstrates that some of these mutants are predominantly affected in their kcat, others in their Km. These results and the determination of the pH and metal ion dependence of selected mutant proteins were used for a tentative assignment for the function of these amino acid residues in the mechanism of phosphodiester bond cleavage by the Serratia nuclease.

Two crystal structures of the leupeptin-trypsin complex

Three-dimensional structures of trypsin with the reversible inhibitor leupeptin have been determined in two different crystal forms. The first structure was determined at 1.7 A resolution with R-factor = 17.7% in the trigonal crystal space group P3(1)21, with unit cell dimensions of a = b = 55.62 A, c = 110.51 A. The second structure was determined at a resolution of 1.8 A with R-factor = 17.5% in the orthorhombic space group P2(1)2(1)2(1), with unit cell dimensions of a = 63.69 A, b = 69.37 A, c = 63.01 A. The overall protein structure is very similar in both crystal forms, with RMS difference for main-chain atoms of 0.27 A. The leupeptin backbone forms four hydrogen bonds with trypsin and a fifth hydrogen bond interaction is mediated by a water molecule. The aldehyde carbonyl of leupeptin forms a covalent bond of 1.42 A length with side-chain oxygen of Ser-195 in the active site. The reaction of trypsin with leupeptin proceeds through the formation of stable tetrahedral complex in which the hemiacetal oxygen atom is pointing out of the oxyanion hole and forming a hydrogen bond with His-57.

Human immunodeficiency virus, type 1 protease substrate specificity is limited by interactions between substrate amino acids bound in adjacent enzyme subsites

The specificity of the retroviral protease is determined by the ability of substrate amino acid side chains to bind into eight individual subsites within the enzyme. Although the subsites are able to act somewhat independently in selection of amino acid side chains that fit into each pocket, significant interactions exist between individual subsites that substantially limit the number of cleavable amino acid sequences. The substrate peptide binds within the enzyme in an extended anti-parallel beta sheet conformation with substrate amino acid side chains adjacent in the linear sequence extending in opposite directions in the enzyme-substrate complex. From this geometry, we have defined both cis and trans steric interactions, which have been characterized by a steady state kinetic analysis of human immunodeficiency virus, type-1 protease using a series of peptide substrates that are derivatives of the avian leukosis/sarcoma virus nucleocapsid-protease cleavage site. These peptides contain both single and double amino acid substitutions in seven positions of the minimum length substrate required by the retroviral protease for specific and efficient cleavage. Steady state kinetic data from the single amino acid substituted peptides were used to predict effects on protease-catalyzed cleavage of corresponding double substituted peptide substrates. The calculated Gibbs' free energy changes were compared with actual experimental values in order to determine how the fit of a substrate amino acid in one subsite influences the fit of amino acids in adjacent subsites. Analysis of these data shows that substrate specificity is limited by steric interactions between pairs of enzyme subsites. Moreover, certain enzyme subsites are relatively tolerant of substitutions in the substrate and exert little effect on adjacent subsites, whereas others are more restrictive and have marked influence on adjacent cis and trans subsites.

Analysis of six protein structures predicted by comparative modeling techniques

The protein structures of six comparative modeling targets were predicted in a procedure that relied on improved energy minimization, without empirical rules, to position all new atoms. The structures of human nucleoside diphosphate kinase NM23-H2, HPr from Mycoplasma capricolum, 2Fe-2S ferredoxin from Haloarcula marismortui, eosinophil-derived neurotoxin (EDN), mouse cellular retinoic acid protein I (CRABP1), and P450eryf were predicted with root mean square deviations on C alpha atoms of 0.69, 0.73, 1.11, 1.48, 1.69, and 1.73 A, respectively, compared to the target crystal structures. These differences increased as the sequence similarity between the target and parent proteins decreased from about 60 to 20% identity. More residues were predicted than form the common region shared by the two crystal structures. In most cases insertions or deletions between the target and the related protein of known structure were not correctly positioned. One two residue insertion in CRABP1 was predicted in the correct conformation, while a nine residue insertion in EDN was predicted in the correct spatial region, although not in the correct conformation. The positions of common cofactors and their binding sites were predicted correctly, even when overall sequence similarity was low.

Bovine pancreatic ribonuclease A as a model of an enzyme with multiple substrate binding sites

Bovine pancreatic ribonuclease A is an enzyme that catalyses the depolymerization of RNA. This process involves the interaction of the enzyme with the polymeric substrate in the active site and its correct alignment on the surface of the enzyme through multiple binding subsites that essentially recognize the negatively charged phosphate groups of the substrate. The enzyme shows a strong specificity for pyrimidine bases at the 3'-position of the phosphodiester bond that is cleaved and a preference for purine bases at the 5'-position and, probably, for guanine at the next position. On the other hand, the enzyme shows a clear preference for polynucleotide substrates over oligonucleotides. In this review the contributions to the catalytic mechanism of some amino-acid residues that are located at non catalytic binding subsites are analysed.

A critical assessment of comparative molecular modeling of tertiary structures of proteins

In spite of the tremendous increase in the rate at which protein structures are being determined, there is still an enormous gap between the numbers of known DNA-derived sequences and the numbers of three-dimensional structures. In order to shed light on the biological functions of the molecules, researchers often resort to comparative molecular modeling. Earlier work has shown that when the sequence alignment is in error, then the comparative model is guaranteed to be wrong. In addition, loops, the sites of insertions and deletions in families of homologous proteins, are exceedingly difficult to model. Thus, many of the current problems in comparative molecular modeling are minor versions of the global protein folding problem. In order to assess objectively the current state of comparative molecular modeling, 13 groups submitted blind predictions of seven different proteins of undisclosed tertiary structure. This assessment shows that where sequence identity between the target and the template structure is high (> 70%), comparative molecular modeling is highly successful. On the other hand, automated modeling techniques and sophisticated energy minimization methods fail to improve upon the starting structures when the sequence identity is low (approximately 30%). Based on these results it appears that insertions and deletions are still major problems. Successfully deducing the correct sequence alignment when the local similarity is low is still difficult. We suggest some minimal testing of submitted coordinates that should be required of authors before papers on comparative molecular modeling are accepted for publication in journals.