Interaction of substrate uridyl 3'',5''-adenosine with ribonuclease A: a molecular dynamics study.

A wealth of information available from x-ray crystallographic structures of enzyme-ligand complexes makes it possible to study interactions at the molecular level. However, further investigation is needed when i) the binding of the natural substrate must be characterized, because ligands in the stable enzyme-ligand complexes are generally inhibitors or the analogs of substrate and transition state, and when ii) ligand binding is in part poorly characterized. We have investigated these aspects in the binding of substrate uridyl 3',5'-adenosine (UpA) to ribonuclease A (RNase A). Based on the systematically docked RNase A-UpA complex resulting from our previous study, we have undertaken a molecular dynamics simulation of the complex with solvent molecules. The molecular dynamics trajectories of this complex are analyzed to provide structural explanations for varied experimental observations on the ligand binding at the B2 subsite of ribonuclease A. The present study suggests that B2 subsite stabilization can be effected by different active site groups, depending on the substrate conformation. Thus when adenosine ribose pucker is O4'-endo, Gln69 and Glu111 form hydrogen-bonding contacts with adenine base, and when it is C2'-endo, Asn71 is the only amino acid residue in direct contact with this base. The latter observation is in support of previous mutagenesis and kinetics studies. Possible roles for the solvent molecules in the binding subsites are described. Furthermore, the substrate conformation is also examined along the simulation pathway to see if any conformer has the properties of a transition state. This study has also helped us to recognize that small but concerted changes in the conformation of the substrate can result in substrate geometry favorable for 2',3' cyclization. The identified geometry is suitable for intraligand proton transfer between 2'-hydroxyl and phosphate oxygen atom. The possibility of intraligand proton transfer as suggested previously and the mode of transfer before the formation of cyclic intermediate during transphosphorylation are discussed.     

Limits to Catalysis by Ribonuclease A

Bovine pancreatic ribonuclease A (RNase A) catalyzes the cleavage of the P-O(5') bond in RNA. Although this enzyme has been the object of much landmark work in bioorganic chemistry, the nature of its rate-limiting transition state and its catalytic rate enhancement had been unknown. Here, the value of k(cat)/K(m) for the cleavage of UpA by wild-type RNase A was found to be inversely related to the concentration of added glycerol. In contrast, the values of k(cat)/K(m) for the cleavage of UpA by a sluggish mutant of RNase A and the cleavage of the poor substrate UpOC(6)H(4)-p-NO(2) by wild-type RNase A were found to be independent of glycerol concentration. Yet, UpA cleavage by the wild-type and mutant enzymes was found to have the same dependence on sucrose concentration, indicating that catalysis of UpA cleavage by RNase A is limited by desolvation. The rate of UpA cleavage by RNase A is maximal at pH 6.0, where k(cat) = 1.4 × 10(3) s(-1) and k(cat)/K(m) = 2.3 × 10(6) M(-1)s(-1) at 25°C. At pH 6.0 and 25°C, the uncatalyzed rate of [5,6-(3)H]Up[3,5,8-(3)H]A cleavage was found to be k(uncat) = 5 × 10(-9) s(-1) (t(1/2) = 4 years). Thus, RNase A enhances the rate of UpA cleavage by 3 × 10(11)-fold by binding to the transition state for P-O(5') bond cleavage with a dissociation constant of <2 × 10(-15) M.     

Probing the substrate specificity of Escherichia coli RNase E using a novel oligonucleotide-based assay

Endoribonuclease RNase E has a central role in both processing and decay of RNA in Escherichia coli, and apparently in many other organisms, where RNase E homologs were identified or their existence has been predicted from genomic data. Although the biochemical properties of this enzyme have been already studied for many years, the substrate specificity of RNase E is still poorly characterized. Here, I have described a novel oligonucleotide-based assay to identify specific sequence determinants that either facilitate or impede the recognition and cleavage of RNA by the catalytic domain of the enzyme. The knowledge of these determinants is crucial for understanding the nature of RNA–protein interactions that control the specificity and efficiency of RNase E cleavage and opens new perspectives for further studies of this multi-domain protein. Moreover, the simplicity and efficiency of the proposed assay suggest that it can be a valuable tool not only for the characterization of RNase E homologs but also for the analysis of other site-specific nucleases.     

Computer modeling of human angiogenin-dinucleotide substrate interaction

Structures of substrate bound human angiogenin complexes have been obtained for the first time by computer modeling. The dinucleotides CpA and UpA have been docked onto human angiogenin using a systematic grid search procedure in torsion and Eulerian angle space. The docking was guided throughout by the similarity of angiogenin-substrate interactions with interactions of RNase A and its substrate. The models were subjected to 1 nanosecond of molecular dynamics to access their stability. Structures extracted from MD simulations were refined by simulated annealing. Stable hydrogen bonds that bridged protein and ligand residues during the MD simulations were taken as restraints for simulated annealing. Our analysis on the MD structures and annealed models explains the substrate specificity of human angiogenin and is in agreement with experimental results. This study also predicts the B2 binding site residues of angiogenin, for which no experimental information is available so far. In the case of one of the substrates, CpA, we have also identified the presence of a water molecule that invariantly bridges the B2 base with the protein. We have compared our results to the RNase A-substrate complex and highlight the similarities and differences.     

Plant mitochondrial RNase P

Molecular investigations in mitochondria of higher plants have to take in account the complicated genomic structure of these organelles and their complex mode of gene expression. Recently tRNA processing activities and particulary RNase P-like activities have been described for mitochondria of mono- and dicot plants. The determined biochemical characteristics of these plant mitochondrial tRNA processing enzymes now allow a comparison to the bacterial prototype from which they evolved. The substrate specifity of the plant mitochondrial RNase P in particular has unique selection parameters distinct from theE. coli RNase P.     

Binding of phosphorothioate oligonucleotides with RNase H1 can cause conformational changes in the protein and alter the interactions of RNase H1 with other proteins

We recently found that toxic PS-ASOs can cause P54nrb and PSF nucleolar mislocalization in an RNase H1-dependent manner. To better understand the underlying mechanisms of these observations, here we utilize different biochemical approaches to demonstrate that PS-ASO binding can alter the conformations of the bound proteins, as illustrated using recombinant RNase H1, P54nrb, PSF proteins and various isolated domains. While, in general, binding of PS-ASOs or ASO/RNA duplexes stabilizes the conformations of these proteins, PS-ASO binding may also cause the unfolding of RNase H1, including both the hybrid binding domain and the catalytic domain. The extent of conformational change correlates with the binding affinity of PS-ASOs to the proteins. Consequently, PS-ASO binding to RNase H1 induces the interaction of RNase H1 with P54nrb or PSF in a 2′-modification and sequence dependent manner, and toxic PS-ASOs tend to induce more interactions than non-toxic PS-ASOs. PS-ASO binding also enhances the interaction between P54nrb and PSF. However, the interaction between RNase H1 and P32 protein can be disrupted upon binding of PS-ASOs. Together, these results suggest that stronger binding of PS-ASOs can cause greater conformational changes of the bound proteins, subsequently affecting protein–protein interactions. These observations thus provide deeper understanding of the molecular basis of PS-ASO-induced protein mislocalization or degradation observed in cells and advance our understanding of why some PS-ASOs are cytotoxic.     

Replacement of residues 8-22 of angiogenin with 7-21 of RNase A selectively affects protein synthesis inhibition and angiogenesis.

The region of human angiogenin containing residues 8-21 is highly conserved in angiogenins from four mammalian species but differs substantially from the corresponding region of the homologous protein ribonuclease A (RNase A). Regional mutagenesis has been employed to replace this segment of angiogenin with the corresponding RNase A sequence, and the activities of the resulting covalent angiogenin/RNase hybrid, designated ARH-III, have been examined. The ribonucleolytic activity of ARH-III is unchanged toward most substrates, including tRNA, naked 18S and 28S rRNA, CpA, CpG, UpA, and UpG. In contrast, the capacity of ARH-III to inhibit cell-free protein synthesis is decreased 20-30-fold compared to that of angiogenin. The angiogenic activity of ARH-III is also different; it is actually more potent. It induces a maximal response in the chick chorioallantoic membrane assay at 0.1 ng per egg, a 10-fold lower dose than required for angiogenin. In addition, binding of ARH-III to the placental ribonuclease inhibitor is increased by at least 1 order of magnitude (Ki less than or equal to 7 x 10(-17) M) compared to angiogenin. Thus, mutation of a highly conserved region of angiogenin markedly affects those properties likely involved in its biological function(s); it does not, however, alter ribonucleolytic activity toward most substrates.     

Zebrafish RNase T2 genes and the evolution of secretory ribonucleases in animals

Background  

Members of the Ribonuclease (RNase) T2 family are common models for enzymological studies, and their evolution has been well characterized in plants. This family of acidic RNases is widespread, with members in almost all organisms including plants, animals, fungi, bacteria and even some viruses. While several biological functions have been proposed for these enzymes in plants, their role in animals is unknown. Interestingly, in vertebrates most of the biological roles of plant RNase T2 proteins are carried out by members of a different family, RNase A. Still, RNase T2 proteins are conserved in these animals     

A computational analysis of ordering in ABCHN at the nematic–isotropic transition

A computational analysis of ordering in the nematogenic compound 4-alkenyl bicyclohexylnitrile has been carried out based on quantum mechanics and intermolecular forces. The evaluation of atomic charge and dipole moment at each atomic center has been carried out using the complete neglect differential overlap (CNDO/2) method. Modified Rayleigh–Schrodinger perturbation theory along with a multicentered- multipole expansion method has been employed to evaluate long-range intermolecular interactions, while a ‘6-exp‘ potential function has been assumed for short-range interactions. The total interaction energy values obtained through these computations were used to calculate the probability of each configuration at room temperature (300 K), the nematic–isotropic transition temperature (364.7 K) and above transition temperature (450 K) using the Maxwell–Boltzmann formula. The various possible configurations during the different modes (i.e., stacking, in-plane and terminal) of interactions have been studied in terms of variation of probability due to small departures from the most probable configurations. An attempt has been made to analyze the characteristic features of liquid crystallinity in terms of their relative order with molecular parameters introduced in this paper.     

RNase 8, a novel RNase A superfamily ribonuclease expressed uniquely in placenta

We report the identification and characterization of the gene encoding the eighth and final human ribonuclease (RNase) of the highly diversified RNase A superfamily. The RNase 8 gene is linked to seven other RNase A superfamily genes on chromosome 14. It is expressed prominently in the placenta, but is not detected in any other tissues examined. Phylogenetic analysis suggests that RNase 7 is the closest relative of RNase 8 and that the pair likely resulted from a recent gene duplication event in primates. Further analysis reveals that the RNase 8 gene has incorporated non-silent mutations at an elevated rate (1.3 × 10^–9 substitutions/site/year) and that orthologous RNase 8 genes from 6 of 10 primate species examined have been deactivated by frameshifting deletions or point mutations at crucial structural or catalytic residues. The ribonucleolytic activity of recombinant human RNase 8 is among the lowest of members of this superfamily and it exhibits neither antiviral nor antibacterial activities characteristic of some other RNase A ribonucleases. The rapid evolution, species-limited deactivation and tissue-specific expression of RNase 8 suggest a unique physiological function and reiterates the evolutionary plasticity of the RNase A superfamily.     

Solution structure of RNase P RNA

The ribonucleoprotein enzyme ribonuclease P (RNase P) processes tRNAs by cleavage of precursor-tRNAs. RNase P is a ribozyme: The RNA component catalyzes tRNA maturation in vitro without proteins. Remarkable features of RNase P include multiple turnovers in vivo and ability to process diverse substrates. Structures of the bacterial RNase P, including full-length RNAs and a ternary complex with substrate, have been determined by X-ray crystallography. However, crystal structures of free RNA are significantly different from the ternary complex, and the solution structure of the RNA is unknown. Here, we report solution structures of three phylogenetically distinct bacterial RNase P RNAs from Escherichia coli, Agrobacterium tumefaciens, and Bacillus stearothermophilus, determined using small angle X-ray scattering (SAXS) and selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) analysis. A combination of homology modeling, normal mode analysis, and molecular dynamics was used to refine the structural models against the empirical data of these RNAs in solution under the high ionic strength required for catalytic activity.     

Molecular modeling of the three-dimensional structure of the bacterial RNase P holoenzyme

Bacterial ribonuclease P (RNase P), an enzyme involved in tRNA maturation, consists of a catalytic RNA subunit and a protein cofactor. Comparative phylogenetic analysis and molecular modeling have been employed to derive secondary and tertiary structure models of the RNA subunits from Escherichia coli (type A) and Bacillus subtilis (type B) RNase P. The tertiary structure of the protein subunit of B.subtilis and Staphylococcus aureus RNase P has recently been determined. However, an understanding of the structure of the RNase P holoenzyme (i.e. the ribonucleoprotein complex) is lacking. We have now used an EDTA-Fe-based footprinting approach to generate information about RNA-protein contact sites in E.coli RNase P. The footprinting data, together with results from other biochemical and biophysical studies, have furnished distance constraints, which in turn have enabled us to build three-dimensional models of both type A and B versions of the bacterial RNase P holoenzyme in the absence and presence of its precursor tRNA substrate. These models are consistent with results from previous studies and provide both structural and mechanistic insights into the functioning of this unique catalytic RNP complex.     

Interactions of the cytotoxic RNase A dimers with the cytosolic ribonuclease inhibitor

Ribonuclease A (RNase A) dimers have been recently found to be endowed with some of the special, i.e., non-catalytic biological activities of RNases, such as antitumor and aspermatogenic activities. These activities have been so far attributed to RNases which can escape the neutralizing action of the cytosolic RNase inhibitor (cRI). However, when the interactions of the two cytotoxic RNase A dimers with cRI were investigated in a quantitative fashion and at the molecular level, the dimers were found to bind cRI with high affinity and to form tight complexes.     

A molecular dynamics study based post facto free energy analysis of the binding of bovine angiogenin with UMP and CMP ligands

Angiogenin is a protein belonging to the superfamily of RNase A. The RNase activity of this protein is essential for its angiogenic activity. Although members of the RNase A family carry out RNase activity, they differ markedly in their strength and specificity. In this paper, we address the problem of higher specificity of angiogenin towards cytosine against uracil in the first base binding position. We have carried out extensive nano-second level molecular dynamics(MD) computer simulations on the native bovine angiogenin and on the CMP and UMP complexes of this protein in aqueous medium with explicit molecular solvent. The structures thus generated were subjected to a rigorous free energy component analysis to arrive at a plausible molecular thermodynamic explanation for the substrate specificity of angiogenin.     

Mutual interactions between subunits of the human RNase MRP ribonucleoprotein complex

The eukaryotic ribonuclease for mitochondrial RNA processing (RNase MRP) is mainly located in the nucleoli and belongs to the small nucleolar ribonucleoprotein (snoRNP) particles. RNase MRP is involved in the processing of pre-rRNA and the generation of RNA primers for mitochondrial DNA replication. A closely related snoRNP, which shares protein subunits with RNase MRP and contains a structurally related RNA subunit, is the pre-tRNA processing factor RNase P. Up to now, 10 protein subunits of these complexes have been described, designated hPop1, hPop4, hPop5, Rpp14, Rpp20, Rpp21, Rpp25, Rpp30, Rpp38 and Rpp40. To get more insight into the assembly of the human RNase MRP complex we studied protein–protein and protein–RNA interactions by means of GST pull-down experiments. A total of 19 direct protein–protein and six direct protein–RNA interactions were observed. The analysis of mutant RNase MRP RNAs showed that distinct regions are involved in the direct interaction with protein subunits. The results provide insight into the way the protein and RNA subunits assemble into a ribonucleoprotein particle. Based upon these data a new model for the architecture of the human RNase MRP complex was generated.     

A positive zymogram for distinguishing among RNase and phosphodiesterases I and II

The positive zymogram, which depends upon indirect production of a formazan from the adenosine released by action of RNase upon UpA, has been modified so that the phosphodiesterases I and II may also be detected. After electrophoretic separation of protein, each of three strips of supporting medium is overlayed with one of three agarose gels containing the enzyme train, adjuncts and (a) adenylyl (3′ → 5′)-uridine, (b) adenylyl (3′ → 5′)-adenosine, or (c), either uridylyl (3′ → 5′)-adenosine or cytidylyl (3′ → 5′)-adenosine or both. The location of purple spots is indicative of the various enzymes as follows: On both (a) and (b) phosphodiesterase I; on both (b) and (c), phosphodiesterase II; on (c) only, RNase (pancreatic type). Positive reactions on all three overlays suggest a combination of enzymes or “nothing dehydrogenase.” Presence of the latter is proved when formazan appears in a fourth overlay devoid of dinucleoside monophosphate.     

Modes of binding of 2'-AMP to RNase T1. A computer modeling study.

The modes of binding of adenosine 2'-monophosphate (2'-AMP) to the enzyme ribonuclease (RNase) T1 were determined by computer modelling studies. The phosphate moiety of 2'-AMP binds at the primary phosphate binding site. However, adenine can occupy two distinct sites--(1) The primary base binding site where the guanine of 2'-GMP binds and (2) The subsite close to the N1 subsite for the base on the 3'-side of guanine in a guanyl dinucleotide. The minimum energy conformers corresponding to the two modes of binding of 2'-AMP to RNase T1 were found to be of nearly the same energy implying that in solution 2'-AMP binds to the enzyme in both modes. The conformation of the inhibitor and the predicted hydrogen bonding scheme for the RNase T1-2'-AMP complex in the second binding mode (S) agrees well with the reported x-ray crystallographic study. The existence of the first mode of binding explains the experimental observations that RNase T1 catalyses the hydrolysis of phosphodiester bonds adjacent to adenosine at high enzyme concentrations. A comparison of the interactions of 2'-AMP and 2'-GMP with RNase T1 reveals that Glu58 and Asn98 at the phosphate binding site and Glu46 at the base binding site preferentially stabilise the enzyme-2'-GMP complex.