Microinjected ribonuclease A as a probe for lysosomal pathways of intracellular protein degradation

There are multiple pathways of intracellular protein degradation, and molecular determinants within proteins appear to target them for particular pathways of breakdown. We use red cell-mediated microinjection to introduce radiolabeled proteins into cultured human fibroblasts in order to follow their catabolism. A well-characterized protein, bovine pancreatic ribonuclease A (RNase A), is localized initially in the cytosol of cells after microinjection, but it is subsequently taken up and degraded by lysosomes. This lysosomal pathway of proteolysis is subject to regulation in that RNase A is taken up and degraded by lysosomes at twice the rate when serum is omitted from the culture medium. Subtilisin cleaves RNase A between residues 20 and 21, and the separated fragments are termed RNase S-peptide (residues 1–20) and RNase S-protein (residues 21–124). Microinjected RNase S-protein is degraded in a serum-independent manner, while RNase S-peptide microinjected alone shows a twofold increase in degradation in response to serum withdrawal. Furthermore, covalent linkage of S-peptide to other proteins prior to microinjection causes degradation of the conjugate to become serum responsive. These results show that recognition of RNase A and certain other proteins for enhanced lysosomal degradation during serum withdrawal is based on some feature of the amino-terminal 20 amino acids. The entire S-peptide is not required for enhanced lysosomal degradation during serum withdrawal because degradation of certain fragments is also responsive to serum. We have identified the essential region to be within residues 7–11 of RNase S-peptide (Lys-Phe-Glu-Arg-Gln; KFERQ). To determine whether related peptides exist in cellular proteins, we raised antibodies to the pentapeptide. Affinity-purified antibodies to KFERQ specifically precipitate 25–35% of cellular proteins, and these proteins are preferentially degraded in response to serum withdrawal. Computer analyses of known protein sequences indicate that proteins degraded by lysosomes at an enhanced rate in response to serum withdrawal contain peptide regions related, but not identical, to KFERQ. We suggest two possible peptide motifs related to KFERQ and speculate about possible mechanisms of selective delivery of proteins to lysosomes based on such peptide regions.     

Limited proteolysis of ribonuclease A with thermolysin in trifluoroethanol.

We have examined the proteolysis of bovine pancreatic ribonuclease A (RNase) by thermolysin when dissolved in aqueous buffer, pH 7.0, in the presence of 50% (v/v) trifluoroethanol (TFE). Under these solvent conditions, RNase acquires a conformational state characterized by an enhanced content of secondary structure (helix) and reduced tertiary structure, as given by CD measurements. It was found that the TFE-resistant thermolysin, despite its broad substrate specificity, selectively cleaves the 124-residue chain of RNase in its TFE state (20-42 degrees C, 6-24 h) at peptide bond Asn 34-Leu 35, followed by a slower cleavage at peptide bond Thr 45-Phe 46. In the absence of TFE, native RNase is resistant to proteolysis by thermolysin. Two nicked RNase species, resulting from cleavages at one or two peptide bonds and thus constituted by two (1-34 and 35-124) (RNase Th1) or three (1-34, 35-45 and 46-124) (RNase Th2) fragments linked covalently by the four disulfide bonds of the protein, were isolated to homogeneity by chromatography and characterized. CD measurements provided evidence that RNase Th1 maintains the overall conformational features of the native protein, but shows a reduced thermal stability with respect to that of the intact species (-delta Tm 16 degrees C); RNase Th2 instead is fully unfolded at room temperature. That the structure of RNase Th1 is closely similar to that of the intact protein was confirmed unambiguously by two-dimensional NMR measurements. Structural differences between the two protein species are located only at the level of the chain segment 30-41, i.e., at residues nearby the cleaved Asn 34-Leu 35 peptide bond. RNase Th1 retained about 20% of the catalytic activity of the native enzyme, whereas RNase Th2 was inactive. The 31-39 segment of the polypeptide chain in native RNase forms an exposed and highly flexible loop, whereas the 41-48 region forms a beta-strand secondary structure containing active site residues. Thus, the conformational, stability, and functional properties of nicked RNase Th1 and Th2 are in line with the concept that proteins appear to tolerate extensive structural variations only at their flexible or loose parts exposed to solvent. We discuss the conformational features of RNase in its TFE-state that likely dictate the selective proteolysis phenomenon by thermolysin.     

A semisynthetic bovine pancreatic ribonuclease containing a unique nitrotyrosine residue

A fully active semisynthetic ribonuclease, RNase 1-118:111-124, may be prepared by enzymatically removing six residues from the COOH terminus of the protein (positions 119-124) and then complementing the inactive RNase 1-118 with a chemically synthesized peptide containing the COOH-terminal 14 residues of the molecule (RNase 111-124) [M. C. Lin, B. Gutte, S. Moore, and R. B. Merrifield (1970) J. Biol. Chem. 245, 5169-5170]. Nitration of tyrosine-115 in the peptide followed by complex formation with RNase 1-118 affords a fully active enzyme containing a unique nitrotyrosine residue in a position which is known and which is very likely to be completely exterior to the active site region. The binding constant between the tetradecapeptide and RNase 1-118 (5 X 10(6) M-1 at pH 6.0) is not changed by the nitration. Crystals of the nitrated complex are isomorphous with those of RNase 1-118:111-124, for which a refined 1.8-A structure has recently been obtained.     

Protein stabilization through phage display

RNase S consists of two proteolytic fragments of RNase A, residues 1-20 (S20) and residues 21-124 (S pro). A 15-mer peptide (S15p) with high affinity for S pro was selected from a phage display library. Peptide residues that are buried in the structure of the wild type complex are conserved in S15p though there are several changes at other positions. Isothermal titration calorimetry studies show that the affinity of S15p is comparable to that of the wild type peptide at 25 degrees C. However, the magnitudes of DeltaH(o) and DeltaC(p) are lower for S15p, suggesting that the thermal stability of the complex is enhanced. In agreement with this prediction, at pH 6, the T(m) of the S15p complex was found to be 10 degrees C higher than that of the wild type complex. This suggests that for proteins where fragment complementation systems exist, phage display can be used to find mutations that increase protein thermal stability.     

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.     

Osmolytes stabilize ribonuclease S by stabilizing its fragments S protein and S peptide to compact folding-competent states

Osmolytes stabilize proteins to thermal and chemical denaturation. We have studied the effects of the osmolytes sarcosine, betaine, trimethylamine-N-oxide, and taurine on the structure and stability of the protein.peptide complex RNase S using x-ray crystallography and titration calorimetry, respectively. The largest degree of stabilization is achieved with 6 m sarcosine, which increases the denaturation temperatures of RNase S and S pro by 24.6 and 17.4 degrees C, respectively, at pH 5 and protects both proteins against tryptic cleavage. Four crystal structures of RNase S in the presence of different osmolytes do not offer any evidence for osmolyte binding to the folded state of the protein or any perturbation in the water structure surrounding the protein. The degree of stabilization in 6 m sarcosine increases with temperature, ranging from -0.52 kcal mol(-1) at 20 degrees C to -5.4 kcal mol(-1) at 60 degrees C. The data support the thesis that osmolytes that stabilize proteins, do so by perturbing unfolded states, which change conformation to a compact, folding competent state in the presence of osmolyte. The increased stabilization thus results from a decrease in conformational entropy of the unfolded state.     

Studies on the interaction of ribonuclease inhibitor with pancreatic ribonuclease involving differential labeling of cysteinyl residues.

Ribonuclease inhibitor (RI) is a protein that forms a very tight complex with ribonucleases (RNases) of the pancreatic type. RI contains 30 thiol groups, some of which are important for the enzyme-inhibitor interaction. To examine which thiols are affected by the binding of RNase, differential labeling experiments were performed. Reaction of porcine RI with the cysteine-specific labeling reagent 4-N,N-dimethylaminoazobenzene-4'-iodoacetamido-2'-sulfonic acid resulted in labeling of an average of 7.4 of the 30 cysteinyl residues. Binding of bovine pancreatic RNase A caused a 3.2-fold reduction in the extent of modification. Peptide mapping showed that in free RI, Cys-57, -371, and -404 were labeled to the greatest extent (yield, 0.4-0.6 mol/mol). RNase A did not protect Cys-57 against modification, whereas the labeling of Cys-371 and -404 was reduced by more than 90%. A second group of residues was labeled to a lesser extent in free RI (yield, 0.04-0.2 mol/mol). Within this group 11 residues were protected by RNase A by more than 90%, 2 were not affected at all, and 7 were protected between 10 and 90%. Seven cysteinyl residues in RI that were protected in the RI.RNase A complex were no longer protected in the RI.S-protein complex. These residues were mainly present in the N-terminal region of RI. However, when the S-peptide was included to yield the RI.RNase S complex, the same pattern of labeling was obtained as with the RI.RNase A complex. Addition of the S-peptide alone had no effect on the labeling. The implications of these observations with respect to RNase binding areas of RI are discussed in relation to the results obtained from the analysis of active RI molecules that contain deletions.     

Molecular modeling of an antigenic complex between a viral peptide and a class I major histocompatibility glycoprotein.

Computer simulation of the conformations of short antigenic peptides (5-10 residues) either free or bound to their receptor, the major histocompatibility complex (MHC)-encoded glycoprotein H-2 Ld, was employed to explain experimentally determined differences in the antigenic activities within a set of related peptides. Starting for each sequence from the most probable conformations disclosed by a pattern-recognition technique, several energy-minimized structures were subjected to molecular dynamics simulations (MD) either in vacuo or solvated by water molecules. Notably, antigenic potencies were found to correlate to the peptides propensity to form and maintain an overall alpha-helical conformation through regular i,i + 4 hydrogen bonds. Accordingly, less active or inactive peptides showed a strong tendency to form i,i + 3 hydrogen bonds at their N-terminal end. Experimental data documented that the C-terminal residue is critical for interaction of the peptide with H-2 Ld. This finding could be satisfactorily explained by a 3-D Q.S.A.R. analysis postulating interactions between ligand and receptor by hydrophobic forces. A 3-D model is proposed for the complex between a high-affinity nonapeptide and the H-2 Ld receptor. First, the H-2 Ld molecule was built from X-ray coordinates of two homologous proteins: HLA-A2 and HLA-Aw68, energy-minimized and studied by MD simulations. With HLA-A2 as template, the only realistic simulation was achieved for a solvated model with minor deviations of the MD mean structure from the X-ray conformation. Water simulation of the H-2 Ld protein in complex with the antigenic nonapeptide was then achieved with the template-derived optimal parameters. The bound peptide retains mainly its alpha-helical conformation and binds to hydrophobic residues of H-2 Ld that correspond to highly polymorphic positions of MHC proteins. The orientation of the nonapeptide in the binding cleft is in accordance with the experimentally determined distribution of its MHC receptor-binding residues (agretope residues). Thus, computer simulation was successfully employed to explain functional data and predicts alpha-helical conformation for the bound peptide.     

Thermodynamics of protein-peptide interactions in the ribonuclease S system studied by titration calorimetry

Two fragments of pancreatic ribonuclease A, a truncated version of S-peptide (residues 1-15) and S-protein (residues 21-124), combine to give a catalytically active complex designated ribonuclease S. Residue 13 in the peptide is methionine. According to the X-ray structure of the complex of S-protein and S-peptide (1-20), this residue is almost fully buried. We have substituted Met-13 with seven other hydrophobic residues ranging in size from glycine to phenylalanine and have determined the thermodynamic parameters associated with the binding of these analogues to S-protein by titration calorimetry at 25 degrees C. These data should provide useful quantitative information for evaluating the contribution of hydrophobic interactions in the stabilization of protein structures.     

A developmentally regulated Streptomyces endoribonuclease resembles ribonuclease E of Escherichia coli

We report that the Streptomyces species S. lividans and S. coelicolor , morphologically complex Gram-positive soil bacteria, contain a developmentally regulated endoribonuclease activity (here named RNase ES) that functionally and immunologically resembles Escherichia coli RNase E. In Streptomyces cells, RNA I — the antisense repressor of replication of ColE1-type plasmids — is cleaved at sites attacked by RNase E. A Mg²⁺-dependent endonuclease that produces RNase E-like cleavages in RNA I and 9S ribosomal RNA was identified in S. lividans cell extracts. A Streptomyces peptide migrating at 70 kDa in SDS/polyacrylamide gels binds to RNase E substrates and reacts with three separate anti-RNase E monoclonal antibodies; the endonucleolytic cleavage activity co-purified with the immunoreactive 70 kDa peptide. We show that RNase ES activity is regulated during the Streptomyces life cycle: activity increased as cells progressed from exponential growth to stationary phase in liquid culture, or from mycelial growth to sporulation on solid media. While mutations that interfere with S. coelicolor development late in its life cycle did not prevent this developmentally associated increase in RNase ES activity, the increase was blocked by a mutation ( bldA ) that interferes early with both morphological and physiological differentiation.     

Subtle functional collective motions in pancreatic-like ribonucleases: from ribonuclease A to angiogenin

The analysis of the dynamic behavior of enzymes is fundamental to structural biology. A direct relationship between protein flexibility and biological function has been shown for bovine pancreatic ribonuclease (RNase A) (Rasmussen et al., Nature 1992;357:423-424). More recently, crystallographic studies have shown that functional motions in RNase A involve the enzyme beta-sheet regions that move concertedly on substrate binding and release (Vitagliano et al., Proteins 2002;46:97-104). These motions have been shown to correspond to intrinsic dynamic properties of the native enzyme by molecular dynamics (MD) simulations. To unveil the occurrence of these collective motions in other members of pancreatic-like superfamily, we carried out MD simulations on human angiogenin (Ang). Essential dynamics (ED) analyses performed on the trajectories reveal that Ang exhibits collective motions similar to RNase A, despite the limited sequence identity (33%) of the two proteins. Furthermore, we show that these collective motions are also present in ensembles of experimentally determined structures of both Ang and RNase A. Finally, these subtle concerted beta-sheet motions were also observed for other two members of the pancreatic-like superfamily by comparing the ligand-bound and ligand-free structures of these enzymes. Taken together, these findings suggest that pancreatic-like ribonucleases share an evolutionary conserved dynamic behavior consisting of subtle beta-sheet motions, which are essential for substrate binding and release.     

Dynamics of the native and the ligand-bound structures of eosinophil cationic protein: network of hydrogen bonds at the catalytic site

Eosinophil Cationic Protein (ECP) is sequentially and structurally similar to ribonuclease A (RNase A). It belongs to the RNase A family of proteins and the RNA catalysis is essential to its biological function. In the present study, we have generated the dinucleotide-bound structures of ECP by docking the dinucleotides to a number of molecular dynamics (MD) generated ECP structures. The stability of the docked enzyme-ligand complexes was ascertained by extensive MD simulations. The modes of ligand binding are explored by essential dynamics studies. The role of water molecules in the stability of the complex and in the catalysis was investigated. The active site residues form a complex network of connections with the ligand and with a water molecule. The catalytic mechanism of the RNA cleavage is examined on the basis of the active site geometry obtained by the simulations.     

Pretransitional structural changes in the thermal denaturation of ribonuclease S and S protein

Two mechanisms have been proposed for the thermal unfolding of ribonuclease S (RNase S). The first is a sequential partial unfolding of the S peptide/S protein complex followed by dissociation, whereas the second is a concerted denaturation/dissociation. The thermal denaturation of ribonuclease S and its fragment, the S protein, were followed with circular dichroism and infrared spectra. These spectra were analyzed by the principal component method of factor analysis. The use of multiple spectral techniques and of factor analysis monitored different aspects of the denaturation simultaneously. The unfolding pathway was compared with that of the parent enzyme ribonuclease A (RNase A), and a model was devised to assess the importance of the dissociation in the unfolding. The unfolding patterns obtained from the melting curves of each protein imply the existence of multiple intermediate states and/or processes. Our data provide evidence that the pretransition in the unfolding of ribonuclease S is due to partial unfolding of the S protein/S peptide complex and that the dissociation occurs at higher temperature. Our observations are consistent with a sequential denaturation mechanism in which at least one partial unfolding step comes before the main conformational transition, which is instead a concerted, final unfolding/dissociation step.     

Primary structure of a base non-specific and adenylic acid preferential ribonuclease from Aspergillus saitoi

The complete primary structure of a base non-specific and adenylic acid preferential RNase (RNase M) from Aspergillus saitoi was determined. The sequence was determined by analysis of the peptides generated by digestion of heat-denatured RNase M with lysylendopeptidase, and the peptides generated from RCM RNase M by digestion with staphylococcal V8 protease or chemical cleavage with BrCN. It consisted of 238 amino acid residues and carbohydrate moiety attached to the 74th asparagine residue. The molecular weight of the protein moiety deduced from the sequence was 26,596. The locations of 10 half cystine residues are almost superimposable on those of RNase Rh from Rhizopus niveus and RNase T2 from Aspergillus oryzae which have similar base specificity. The homology between RNase M and RNase Rh and RNase T2 amounted to 97 and 160 amino acid residues, respectively. The amino acid sequences conserved in the three RNases are concentrated around the three histidine residues, which are supposed to form part of the active sites of these RNases.     

Characterization of the histidine proton nuclear magnetic resonances of a semisynthetic ribonuclease

The proton magnetic resonance spectrum at 300 MHz of the histidine residues in a semisynthetic derivative of bovine pancreatic ribonuclease (RNase A) has been determined. The derivative RNase 1-118 . 111-124 was prepared by enzymically removing six residues from the COOH terminus of the protein (positions 119-124) and then complementing the inactive RNase 1-118 with a chemically synthesized peptide containing the COOH-terminal 14 residues of ribonuclease (RNase 111-124) [Lin, M.C., Gutte, B., Moore, S., & Merrifield, R.B. (1970) J. Biol. Chem. 245, 5169-5170]. Comparison of the line positions of the C(2)-1H resonances of these residues and of their pH dependence with those reported by other workers has allowed assignment of the resonances to individual residues, as well as the determination of individual pK values for histidine-12, histidine-105, and histidine-119. The assignment of histidine-119 was confirmed by the use of a selectively deuterated derivative. The titration behavior of all four histidine residues is indistinguishable from that observed by others for bovine pancreatic ribonuclease A. Partial dissociation of the noncovalent semisynthetic complex was evident at 30 degrees C, pH 4.0, 0.3 M NaCl; pertinent spectra were analyzed to provide an estimate of the association constant between the component chains under these conditions of 1.9 X 10(3) M-1.     

A covalent angiogenin/ribonuclease hybrid with a fourth disulfide bond generated by regional mutagenesis

Human angiogenin is a blood vessel inducing protein whose primary structure displays 33% identity to that of bovine pancreatic ribonuclease A (RNase A). Angiogenin catalyzes limited cleavage of 18S and 28S ribosomal RNA and is several orders of magnitude less potent than RNase A toward conventional substrates. A striking structural difference between angiogenin and RNase is the virtual absence of sequence similarity within the region of RNase that contains the Cys-65--Cys-72 disulfide bond. Indeed, angiogenin lacks this disulfide linkage. The present report describes the use of regional mutagenesis to generate a covalent angiogenin/RNase hybrid protein, ARH-I, where residues 58-70 of angiogenin have been replaced by the corresponding segment of RNase A (residues 59-73). The protein expressed in Escherichia coli readily folds at pH 8.5 to form the four expected disulfide bonds. The in vivo angiogenic potency of ARH-I is markedly diminished compared with that of angiogenin when examined using the chick chorioallantoic membrane assay. In contrast, its enzymatic activity is dramatically increased. With high molecular weight wheat germ RNA and tRNA, ARH-I is 660- and 300-fold more active than angiogenin, respectively, while with poly(uridylic acid), poly(cytidylic acid), cytidylyl(3'----5')adenosine (CpA), and uridylyl(3'----5')adenosine (UpA) activity is enhanced by about 200-fold. In addition, the specificity of ARH-I toward dinucleoside 3',5'-phosphates is qualitatively similar to RNase A; while angiogenin prefers cytidylyl(3'----5')guanosine (CpG) to UpA, both RNase and the hybrid prefer UpA to CpG. ARH-I also displays greater than 10-fold enhanced activity toward rRNA in intact ribosomes, while abolishing the capacity of the ribosome to support cell-free protein synthesis. The enhanced enzymatic properties of ARH-I parallel a 2-fold increase in chemical reactivity of active-site lysine and histidine residues based on rates of chemical modification. The data indicate that introduction of a region of RNase A containing the Cys-65--Cys-72 disulfide bond into angiogenin dramatically increases RNase-like enzymatic activity while reducing its angiogenicity.     

Inhibition of human pancreatic ribonuclease by the human ribonuclease inhibitor protein

The ribonuclease inhibitor protein (RI) binds to members of the bovine pancreatic ribonuclease (RNase A) superfamily with an affinity in the femtomolar range. Here, we report on structural and energetic aspects of the interaction between human RI (hRI) and human pancreatic ribonuclease (RNase 1). The structure of the crystalline hRI x RNase 1 complex was determined at a resolution of 1.95 A, revealing the formation of 19 intermolecular hydrogen bonds involving 13 residues of RNase 1. In contrast, only nine such hydrogen bonds are apparent in the structure of the complex between porcine RI and RNase A. hRI, which is anionic, also appears to use its horseshoe-shaped structure to engender long-range Coulombic interactions with RNase 1, which is cationic. In accordance with the structural data, the hRI.RNase 1 complex was found to be extremely stable (t(1/2)=81 days; K(d)=2.9 x 10(-16) M). Site-directed mutagenesis experiments enabled the identification of two cationic residues in RNase 1, Arg39 and Arg91, that are especially important for both the formation and stability of the complex, and are thus termed "electrostatic targeting residues". Disturbing the electrostatic attraction between hRI and RNase 1 yielded a variant of RNase 1 that maintained ribonucleolytic activity and conformational stability but had a 2.8 x 10(3)-fold lower association rate for complex formation and 5.9 x 10(9)-fold lower affinity for hRI. This variant of RNase 1, which exhibits the largest decrease in RI affinity of any engineered ribonuclease, is also toxic to human erythroleukemia cells. Together, these results provide new insight into an unusual and important protein-protein interaction, and could expedite the development of human ribonucleases as chemotherapeutic agents.     

Dissecting the effect of trifluoroethanol on ribonuclease A. Subtle structural changes detected by nonspecific proteases.

With the aim to distinguish between local and global conformational changes induced by trifluoroethanol in RNase A, spectroscopic and activity measurements in combination with proteolysis by unspecific proteases have been exploited for probing structural transitions of RNase A as a function of trifluoroethanol concentration. At > 30% (v/v) trifluoroethanol (pH 8.0; 25 degrees C), circular dichroism and fluorescence spectroscopy indicate a cooperative collapse of the tertiary structure of RNase A coinciding with the loss of its enzymatic activity. In contrast to the denaturation by guanidine hydrochloride, urea or temperature, the breakdown of the tertiary structure in trifluoroethanol is accompanied by an induction of secondary structure as detected by far-UV circular dichroism spectroscopy. Proteolysis with the nonspecific proteases subtilisin Carlsberg or proteinase K, both of which attack native RNase A at the Ala20-Ser21 peptide bond, yields refined information on conformational changes, particularly in the pretransition region. While trifluoroethanol at concentrations > 40% results in a strong increase of the rate of proteolysis and new primary cleavage sites (Tyr76-Ser77, Met79-Ser80) were identified, the rate of proteolysis at trifluoroethanol concentrations < 40% (v/v) is much smaller (up to two orders of magnitude) than that of the native RNase A. The proteolysis data point to a decreased flexibility in the surrounding of the Ala20-Ser21 peptide bond, which we attribute to subtle conformational changes of the ribonuclease A molecule. These changes, however, are too marginal to alter the overall catalytic and spectroscopic properties of ribonuclease A.     

Recognition of ribonuclease A by 3'-5'-pyrophosphate-linked dinucleotide inhibitors: a molecular dynamics/continuum electrostatics analysis

The proteins of the pancreatic ribonuclease A (RNase A) family catalyze the cleavage of the RNA polymer chain. The development of RNase inhibitors is of significant interest, as some of these compounds may have a therapeutic effect in pathological conditions associated with these proteins. The most potent low molecular weight inhibitor of RNase reported to date is the compound 5′-phospho-2′-deoxyuridine-3-pyrophosphate (P→5)-adenosine-3-phosphate (pdUppA-3′-p). The 3′,5′-pyrophosphate group of this compound increases its affinity and introduces structural features which seem to be unique in pyrophosphate-containing ligands bound to RNase A, such as the adoption of a syn conformation by the adenosine base at RNase subsite B₂ and the placement of the 5′-β-phosphate of the adenylate (instead of the α-phosphate) at subsite P₁ where the phosphodiester bond cleavage occurs. In this work, we study by multi-ns molecular dynamics simulations the structural properties of RNase A complexes with the ligand pdUppA-3′-p and the related weaker inhibitor dUppA, which lacks the 3′ and 5′ terminal phosphate groups of pdUppA-3′-p. The simulations show that the adenylate 5′-β-phosphate binding position and the adenosine syn orientation constitute robust structural features in both complexes, stabilized by persistent interactions with specific active-site residues of subsites P₁ and B₂. The simulation structures are used in conjunction with a continuum-electrostatics (Poisson-Boltzmann) model, to evaluate the relative binding affinity of the two complexes. The computed relative affinity of pdUppA-3′-p varies between −7.9 kcal/mol and −2.8 kcal/mol for a range of protein/ligand dielectric constants (ε_p) 2–20, in good agreement with the experimental value (−3.6 kcal/mol); the agreement becomes exact with ε_p = 8. The success of the continuum-electrostatics model suggests that the differences in affinity of the two ligands originate mainly from electrostatic interactions. A residue decomposition of the electrostatic free energies shows that the terminal phosphate groups of pdUppA-3′-p make increased interactions with residues Lys⁷ and Lys⁶⁶ of the more remote sites P₂ and P₀, and His¹¹⁹ of site P₁.     

Structural insight into the mechanism of double-stranded RNA processing by ribonuclease III

Members of the ribonuclease III (RNase III) family are double-stranded RNA (dsRNA) specific endoribonucleases characterized by a signature motif in their active centers and a two-base 3' overhang in their products. While Dicer, which produces small interfering RNAs, is currently the focus of intense interest, the structurally simpler bacterial RNase III serves as a paradigm for the entire family. Here, we present the crystal structure of an RNase III-product complex, the first catalytic complex observed for the family. A 7 residue linker within the protein facilitates induced fit in protein-RNA recognition. A pattern of protein-RNA interactions, defined by four RNA binding motifs in RNase III and three protein-interacting boxes in dsRNA, is responsible for substrate specificity, while conserved amino acid residues and divalent cations are responsible for scissile-bond cleavage. The structure reveals a wealth of information about the mechanism of RNA hydrolysis that can be extrapolated to other RNase III family members.