The role of lysine-41 in RNase A catalysis - A quantum chemical study on the active site ligand complex

Several mechanisms have been proposed to explain the action of enzymes at the atomic level. Among them, the recent proposals involving short hydrogen bonds as a step in catalysis by Gerlt and Gassman [1] and proton transfer through low barrier hydrogen bonds (LBHBs) [2, 3] have attracted attention. There are several limitations to experimentally testing such hypotheses, Recent developments in computational methods facilitate the study of active site-ligand complexes to high levels of accuracy, Our previous studies, which involved the docking of the dinucleotide substrate UpA to the active site of RNase A [4, 5], enabled us to obtain a realistic model of the ligand-bound active site of RNase A. From these studies, based on empirical potential functions, we were able to obtain the molecular dynamics averaged coordinates of RNase A, bound to the ligand UpA. A quantum mechanical study is required to investigate the catalytic process which involves the cleavage and formation of covalent bonds. In the present study, we have investigated the strengths of some of the hydrogen bonds between the active site residues of RNase A and UpA at the ab initio quantum chemical level using the molecular dynamics averaged coordinates as the starting point. The 49 atom system and other model systems were optimized at the 3-21G level and the energies of the optimized systems were obtained at the 6-31G* level. The results clearly indicate the strengthening of hydrogen bonds between neutral residues due to the presence of charged species at appropriate positions. Such a strengthening manifests itself in the form of short hydrogen bonds and a low barrier for proton transfer. In the present study, the proton transfer between the 2'-OH of ribose (from the substrate) and the imidazole group from the H12 of RNase A is influenced by K41, which plays a crucial role in strengthening the neutral hydrogen bond, reducing the barrier for proton transfer.

Characterization of substrate UpA binding to RNase A--computer modelling and energetics approach.

In the past two decades RNase A has been the focus of diverse investigations in order to understand the nature of substrate binding and to know the mechanism of enzyme action. Although this system is reasonably well characterized from the view point of some of the binding sites, the details of interactions in the second base binding (B2) site is insufficient. Further, the nature of ligand-protein interaction is elucidated generally by studies on RNase A-substrate analog complexes (mainly with the help of X-ray crystallography). Hence, the details of interactions at atomic level arising due to substrates are inferred indirectly. In the present paper, the dinucleotide substrate UpA is fitted into the active site of RNase A Several possible substrate conformations are investigated and the binding modes have been selected based on Contact Criteria. Thus identified RNase A-UpA complexes are energy minimized in coordinate space and are analysed in terms of conformations, energetics and interactions. The best possible ligand conformations for binding to RNase A are identified by experimentally known interactions and by the energetics. Upon binding of UpA to RNase A the changes associated,with protein back bone, Side chains in general and at the binding sites in particular are described. Further, the detailed interactions between UpA and RNase A are characterized in terms of hydrogen bonds and energetics. An extensive study has helped in interpreting the diverse results obtained from a number of experiments and also in evaluating the extent of changes the protein and the substrate undergo in order to maximize their interactions.

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.

MS_RHII-RSD, a Dual-Function RNase HII-(p)ppGpp Synthetase from Mycobacterium smegmatis

In the noninfectious soil saprophyte Mycobacterium smegmatis, intracellular levels of the stress alarmones guanosine tetraphosphate and guanosine pentaphosphate, together termed (p)ppGpp, are regulated by the enzyme Rel(Msm). This enzyme consists of a single, bifunctional polypeptide chain that is capable of both synthesizing and hydrolyzing (p)ppGpp. The rel(Msm), knockout strain of M. smegmatis (Delta rel(Msm)) is expected to show a (p)ppGpp null (p)ppGpp(0)] phenotype. Contrary to this expectation, the strain is capable of synthesizing (p)ppGpp in vivo. In this study, we identify and functionally characterize the open reading frame (ORF), MSMEG_5849, that encodes a second functional (p)ppGpp synthetase in M. smegmatis. In addition to (p)ppGpp synthesis, the 567-amino-acid-long protein encoded by this gene is capable of hydrolyzing RNA(.)DNA hybrids and bears similarity to the conventional RNase HII enzymes. We have classified this protein as actRel(Msm) in accordance with the recent nomenclature proposed and have named it MS_RHII-RSD, indicating the two enzymatic activities present RHII, RNase HII domain, originally identified as (d) under bar omain of (u) under bar nknown (f) under bar unction 429 (DUF429), and RSD, RelA_SpoT nucleotidyl transferase domain, the SYNTH domain responsible for (p)ppGpp synthesis activity]. MS_RHII-RSD is expressed and is constitutively active in vivo and behaves like a monofunctional (p)ppGpp synthetase in vitro. The occurrence of the RNase HII and (p)ppGpp synthetase domains together on the same polypeptide chain is suggestive of an in vivo role for this novel protein as a link connecting the essential life processes of DNA replication, repair, and transcription to the highly conserved stress survival pathway, the stringent response.

Dynamics of ribonuclease A and ribonuclease S: Computational and experimental studies

RNase S is a complex consisting of two proteolytic fragments of RNase A: the S peptide (residues 1-20) and S protein (residues 21-124). RNase S and RNase A have very similar X-ray structures and enzymatic activities. previous experiments have shown increased rates of hydrogen exchange and greater sensitivity to tryptic cleavage for RNase S relative to RNase A. It has therefore been asserted that the RNase S complex is considerably more dynamically flexible than RNase A. In the present study we examine the differences in the dynamics of RNaseS and RNase A computationally, by MD simulations, and experimentally, using trypsin cleavage as a probe of dynamics. The fluctuations around the average solution structure during the simulation were analyzed by measuring the RMS deviation in coordinates. No significant differences between RNase S and RNase A dynamics were observed in the simulations. We were able to account for the apparent discrepancy between simulation and experiment by a simple model, According to this model, the experimentally observed differences in dynamics can be quantitatively explained by the small amounts of free S peptide and S protein that are present in equilibrium with the RNase S complex. Thus, folded RNase A and the RNase S complex have identical dynamic behavior, despite the presence of a break in polypeptide chain between residues 20 and 21 in the latter molecule. This is in contrast to what has been widely believed for over 30 years about this important fragment complementation system.

Ammonium ions prevent methylation of uridine to ribothymidine in Azotobacter vinelandii tRNA

It was shown that tRNA from Azotobacter vinelandii grown in the presence of ammonium chloride lacks ribothymidine while that grown in the absence of the ammonium salt contains this modified nucleoside. [32P]-Labelled tRNA from this organism grown in a medium containing the ammonium salt was digested with RNase T1 and the pseudouridinecontaining tetranucleotide, common to all tRNAs was isolated and analysed for the nucleoside replacing the ribothymidine. It was found to be uridine. Cells previously labelled with [32P]- phosphate in the ammonium salt medium were washed and incubated in the ammonium saltfree medium to test whether ribothymidine would be formed upon removal of the ammoniumions. Methylation of the uridine did not take place.

Influence of deamidation(s) in the 67-74 region of ribonuclease on its refolding

The influence of chemical mutation featuring the selective conversion of asparagine or glutamine to aspartic or glutamic acid, respectively, on the kinetics of refolding of reduced RNase has been studied. The monodeamidated derivatives of RNase A, viz. RNase Aa1a, Aa1b, and Aa1c having their deamidations in the region 67-74, were found to regain nearly their original enzymatic activity. However, a marked difference in the kinetics of refolding is seen, the order of regain of enzymic activity being RNase A greater than Aa1c congruent to Aa1a greater than Aa1b. The similarities in the distinct elution positions on Amberlite XE-64, gel electrophoretic mobilities, and u.v. spectra of reoxidized and native derivatives indicated that the native structures are formed. The slower rate of reappearance of enzymic activity in the case of the monodeamidated derivatives appears to result from altered interactions in the early stages of refolding. The roles of some amino acid residues of the 67-74 region in the pathway of refolding of RNase A are discussed.

solation and characterization of monodeamidated derivatives of bovine pancreatic Ribonuclease A

The isolation and characterization of the initial intermediates formed during the irreversible acid denaturation of enzyme Ribonuclease A are described. The products obtained when RNase A is maintained in 0.5 M HCl at 30° for periods up to 20 h have been analyzed by ion-exchange chromatography on Amberlite XE-64. Four distinct components were found to elute earlier to RNase A; these have been designated RNase Aa2, Aa1c, Aa1b, and Aa1a in order of their elution. With the exception of RNase Aa2, the other components are nearly as active as RNase A. Polyacrylamide gel electrophoresis at near-neutral pH indicated that RNase Aa1a, Aa1b, and Aa1c are monodeamidated derivatives of RNase A; RNase Aa2 contains, in addition, a small amount of a dideamidated component. RNase Aa2, which has 75% enzymic activity as compared to RNase A, consists of dideamidated and higher deamidated derivatives of RNase A. Except for differences in the proteolytic susceptibilities at an elevated temperature or acidic pH, the monodeamidated derivatives were found to have very nearly the same enzymic activity and the compact folded structure as the native enzyme. Fingerprint analyses of the tryptic peptides of monodeamidated derivatives have shown that the deamidations are restricted to an amide cluster in the region 67–74 of the polypeptide chain. The initial acid-catalyzed deamidation occurs in and around the 65–72 disulfide loop giving rise to at least three distinct monodeamidated derivatives of RNase A without an appreciable change in the catalytic activity and conformation of the ribonuclease molecule. Significance of this specific deamidation occurring in highly acidic conditions, and the biological implications of the physiological deamidation reactions of proteins are discussed.

Effect of cupric-isonicotinic acid hydrazide complex on avian myeloblastosis virus reverse transcriptase and its associated enzyme activities

Cupric complex of isonicotinic acid hydrazide inhibits DNA synthesis by avian myloblastosis virus reverse transcriptase. This inhibition occurs in the presence of either ribonucleotide or deoxyribonucleotide templates. The inhibition of reverse transcriptase by cupric-INH complex is considerably reduced when stored or proteolytically cleaved enzyme was used in the reaction. The complex also inhibits the reverse transciptase-associated RNase H activity. The cupric-isonicotinic acid hydrazide complex cleaves pBR 322 from I DNA into smaller molecules in the presence or absence of reverse transcriptase-associated endonuclease. However, in the presence of the enzyme the DNA is cleaved to a greater extent.

RNA polymerase activity in isolated nuclei of Nicotiana sanderae callus: Characteristics and modulation during differentiation

Isolated nuclei from differentiating cultures of Nicotiana sanderae showed increased levels of RNA polymerase activity as compared to the nuclei from callus cultures. The RNA synthetic activity was dependent on nucleotide triphosphates and Mg2+ and was destroyed by RNase. Maximum activity was obtained in the presence of 50 mM (NH4)2 SO4 and α-amanitin inhibited 40% and 55% of the activity in the nuclei from callus and differentiating tissue respectively. The nuclei from differentiating tissue elicited a 3-fold increase in RNA polymerase I and a 4-fold augmentation in RNA polymerase II activities.

The influence of esterification of carboxyl groups of ribonuclease-a on its structure and immunological activity

The effect of modification of carboxyl groups of Ribonuclease-Aa on the enzymatic activity and the antigenic structure of the protein has been studied. Modification of four of the eleven free carboxyl groups of the protein by esterification in anhydrous methanol/0.1 M hydrochloric acid resulted in nearly 80% loss in enzymatic activity but had very little influence on the antigenic structure of the protein. Further increases in the modification of the carboxyl groups caused a progressive loss in immunological activity, and the fully methylated RNase-A exhibited nearly 30% immunological activity. Concomitant with this change in the antigenic structure of the protein, the ability of the molecule to complement with RNase-S-protein increased, clearly indicating the unfolding of the peptide "tail" from the remainder of the molecule. The susceptibility to proteolysis, accessibility of methionine residues for orthobenzoquinone reaction and the loss in immunological activity of the more extensively esterified derivatives of RNase-A are suggestive of the more flexible conformation of these derivatives as compared with the compact native conformation. The fact that even the fully methylated RNase-A retains nearly 30% of its immunological activity suggested that the modified protein contained antibody recognizable residual native structure, which presumably accommodates some antigenic determinants.

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 i? 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.

Thermodynamic and structural studies of cavity formation in proteins suggest that loss of packing interactions rather than the hydrophobic effect dominates the observed energetics

The hydrophobic effect is widely believed to be an important determinant of protein stability. However, it is difficult to obtain unambiguous experimental estimates of the contribution of the hydrophobic driving force to the overall free energy of folding. Thermodynamic and structural studies of large to small substitutions in proteins are the most direct method of measuring this contribution. We have substituted the buried residue Phe8 in RNase S with alanine, methionine, and norleucine, Binding thermodynamics and structures were characterized by titration calorimetry and crystallography, respectively. The crystal structures of the RNase S F8A, F8M, and F8Nle mutants indicate that the protein tolerates the changes without any main chain adjustments, The correlation of structural and thermodynamic parameters associated with large to small substitutions was analyzed for nine mutants of RNase S as well as 32 additional cavity-containing mutants of T4 lysozyme, human lysozyme, and barnase. Such substitutions were typically found to result in negligible changes in Delta C-p and positive values of both Delta Delta H degrees and aas of folding. Enthalpic effects were dominant, and the sign of Delta Delta S is the opposite of that expected from the hydrophobic effect. Values of Delta Delta G degrees and Delta Delta H degrees correlated better with changes in packing parameters such as residue depth or occluded surface than with the change in accessible surface area upon folding. These results suggest that the loss of packing interactions rather than the hydrophobic effect is a dominant contributor to the observed energetics for large to small substitutions. Hence, estimates of the magnitude of the hydrophobic driving force derived from earlier mutational studies are likely to be significantly in excess of the actual value.

Computer modelling studies on the mechanism of action of ribonuclease T1

The mechanism of action of ribonuclease (RNase) T1 is still a matter of considerable debate as the results of x-ray, 2-D nmr and site-directed mutagenesis studies disagree regarding the role of the catalytically important residues. Hence computer modelling studies were carried out by energy minimisation of the complexes of RNase T1 and some of its mutants (His40Ala, His40Lys, and Glu58Ala) with the substrate guanyl cytosine (GpC), and of native RNase T1 with the reaction intermediate guanosine 2',3'-cyclic phosphate (G greater than p). The puckering of the guanosine ribose moiety in the minimum energy conformer of the RNase T1-GpC (substrate) complex was found to be O4'-endo and not C3'-endo as in the RNase T1-3'-guanylic acid (inhibitor/product) complex. A possible scheme for the mechanism of action of RNase T1 has been proposed on the basis of the arrangement of the catalytically important amino acid residues His40, Glu58, Arg77, and His92 around the guanosine ribose and the phosphate moiety in the RNase T1-GpC and RNase T1-G greater than p complexes. In this scheme, Glu58 serves as the general base group and His92 as the general acid group in the transphosphorylation step. His40 may be essential for stabilising the negatively charged phosphate moiety in the enzyme-transition state complex.

Modes of binding of guanosine monophosphates to ribonuclease T1 - A computer modelling study

Computer-modelling studies on the modes of binding of the three guanosine monophosphate inhibitors 2'-GMP, 3'-GMP, and 5'-GMP to ribonuclease (RNase) T1 have been carried out by energy minimization in Cartesian-coordinate space. The inhibitory power was found to decrease in the order 2'-GMP > 3'-GMP > 5'-GMP in agreement with the experimental observations. The ribose moiety was found to form hydrogen bonds with the protein in all the enzyme-inhibitor complexes, indicating that it contributes to the binding energy and does not merely act as a spacer between the base and the phosphate moieties as suggested earlier. 2'-GMP and 5'-GMP bind to RNase T1 in either of the two ribose puckered forms (with C3'-endo more favoured over the C2'-endo) and 3'-GMP binds to RNase T1 predominantly in C3'-endo form. The catalytically important residue His-92 was found to form hydrogen bond with the phosphate moiety in all the enzyme-inhibitor complexes, indicating that this residue may serve as a general acid group during catalysis. Such an interaction was not found in either X-ray or two-dimensional NMR studies.