Transition-state stabilization in Escherichia coli ribonuclease P RNA-mediated cleavage of model substrates

We have used model substrates carrying modified nucleotides at the site immediately 5′ of the canonical RNase P cleavage site, the −1 position, to study Escherichia coli RNase P RNA-mediated cleavage. We show that the nucleobase at −1 is not essential but its presence and identity contribute to efficiency, fidelity of cleavage and stabilization of the transition state. When U or C is present at −1, the carbonyl oxygen at C2 on the nucleobase contributes to transition-state stabilization, and thus acts as a positive determinant. For substrates with purines at −1, an exocyclic amine at C2 on the nucleobase promotes cleavage at an alternative site and it has a negative impact on cleavage at the canonical site. We also provide new insights into the interaction between E. coli RNase P RNA and the −1 residue in the substrate. Our findings will be discussed using a model where bacterial RNase P cleavage proceeds through a conformational-assisted mechanism that positions the metal(II)-activated H₂O for an in-line attack on the phosphorous atom that leads to breakage of the phosphodiester bond.     

Catalytic mechanism of Escherichia coli ribonuclease III: kinetic and inhibitor evidence for the involvement of two magnesium ions in RNA phosphodiester hydrolysis

Escherichia coli ribonuclease III (RNase III; EC 3.1.24) is a double-stranded(ds)-RNA-specific endonuclease with key roles in diverse RNA maturation and decay pathways. E.coli RNase III is a member of a structurally distinct superfamily that includes Dicer, a central enzyme in the mechanism of RNA interference. E.coli RNase III requires a divalent metal ion for activity, with Mg²⁺ as the preferred species. However, neither the function(s) nor the number of metal ions involved in catalysis is known. To gain information on metal ion involvement in catalysis, the rate of cleavage of the model substrate R1.1 RNA was determined as a function of Mg²⁺ concentration. Single-turnover conditions were applied, wherein phosphodiester cleavage was the rate-limiting event. The measured Hill coefficient (n_H) is 2.0 ± 0.1, indicative of the involvement of two Mg²⁺ ions in phosphodiester hydrolysis. It is also shown that 2-hydroxy-4H-isoquinoline-1,3-dione—an inhibitor of ribonucleases that employ two divalent metal ions in their catalytic sites—inhibits E.coli RNase III cleavage of R1.1 RNA. The IC₅₀ for the compound is 14 μM for the Mg²⁺-supported reaction, and 8 μM for the Mn²⁺-supported reaction. The compound exhibits noncompetitive inhibitory kinetics, indicating that it does not perturb substrate binding. Neither the O-methylated version of the compound nor the unsubstituted imide inhibit substrate cleavage, which is consistent with a specific interaction of the N-hydroxyimide with two closely positioned divalent metal ions. A preliminary model is presented for functional roles of two divalent metal ions in the RNase III catalytic mechanism.     

Protein primary structure using orthogonal fragmentation techniques in Fourier transform mass spectrometry

Proteomics analysis using tandem mass spectrometry requires informative backbone fragmentation of peptide ions. Collision-activated dissociation (CAD) of cations alone is not sufficiently informative to satisfy all requirements. Thus, there is a need to supplement CAD with a complementary fragmentation technique. Electron capture dissociation (ECD) is complementary to collisional excitation in terms of the cleavage of a different bond (N-C_α versus C-N bond) and other properties. CAD-ECD combination improves protein identification and enables high-throughput de novo sequencing of peptides. ECD and its variants are also useful in mapping labile post-translational modifications in proteins and isomer differentiation; for example, distinguishing Ile from Leu, iso-Asp from Asp and even D- from L-amino acid residues.     

Detection of phosphodiester adducts formed by the reaction of benzo[a]pyrene diol epoxide with 2'-deoxynucleotides using collision-induced dissociation electrospray ionization tandem mass spectrometry

In this study, we investigated the products formed following the reaction of benzo[a]pyrene-7,8-dihydrodiol-9,10-epoxide (B[a]PDE) with 2′-deoxynucleoside 3′-monophosphates. The B[a]PDE plus 2′-deoxynucleotide reaction mixtures were purified using solid phase extraction (SPE) and subjected to HPLC with fluorescence detection. Fractions corresponding to reaction product peaks were collected and desalted using SPE prior to analysis for the presence of molecular ions corresponding to m/z 648, 632, 608 and 623 [M−H]⁻ consistent with B[a]PDE adducted (either on the base or phosphate group) 2′-deoxynucleotides of guanine, adenine, cytosine and thymine, respectively, using LC-ESI-MS/MS collision-induced dissociation (CID). Reaction products were identified having CID product ion spectra containing product ions at m/z 452, 436 and 412 [(B[a]Ptriol+base)−H]⁻, resulting from cleavage of the glycosidic bond between the 2′-deoxyribose and base, corresponding to B[a]PDE adducts of guanine, adenine and cytosine, respectively. Further reaction products were identified having unique CID product ion spectra characteristic of B[a]PDE adduct formation with the phosphate group of the 2′-deoxynucleotide. The presence of product ions at m/z 399 and 497 were observed for all four 2′-deoxynucleotides, corresponding to [(B[a]Ptriol+phosphate)−H]⁻ and [(2′-deoxyribose+phosphate+B[a]Ptriol)−H]⁻, respectively. In conclusion, this investigation provides the first direct evidence for the formation of phosphodiester adducts by B[a]PDE following reaction with 2′-deoxynucleotides.     

Protein primary structure using orthogonal fragmentation techniques in Fourier transform mass spectrometry

Proteomics analysis using tandem mass spectrometry requires informative backbone fragmentation of peptide ions. Collision-activated dissociation (CAD) of cations alone is not sufficiently informative to satisfy all requirements. Thus, there is a need to supplement CAD with a complementary fragmentation technique. Electron capture dissociation (ECD) is complementary to collisional excitation in terms of the cleavage of a different bond (N-Calpha versus C-N bond) and other properties. CAD-ECD combination improves protein identification and enables high-throughput de novo sequencing of peptides. ECD and its variants are also useful in mapping labile post-translational modifications in proteins and isomer differentiation; for example, distinguishing Ile from Leu, iso-Asp from Asp and even D- from L-amino acid residues.     

The bulge region of HIV-1 TAR RNA binds metal ions in solution

Binding of Mg²⁺, Ca²⁺ and Co(NH₃)₆³⁺ ions to the HIV-1 TAR RNA in solution was analysed by ¹⁹F NMR spectroscopy, metal ion-induced RNA cleavages and Brownian dynamics (BD) simulations. Chemically synthesised 29mer oligoribonucleotides of the TAR sequence labelled with 5-fluorouridine (FU) were used for ¹⁹F NMR-monitored metal ion titration. The chemical shift changes of fluorine resonances FU-23, FU-25 and FU-40 upon titration with Mg²⁺ and Ca²⁺ ions indicated specific, although weak, binding at the bulge region with the dissociation constants (K_d) of 0.9 ± 0.6 and 2.7 ± 1.7 mM, respectively. Argininamide, inducing largest ¹⁹F chemical shifts changes at FU-23, was used as a reference ligand (K_d = 0.3 ± 0.1 mM). In the Pb²⁺-induced TAR RNA cleavage experiment, strong and selective cleavage of the C24-U25 phosphodiester bond was observed, while Mg²⁺ and Ca²⁺ induced cuts at all 3-nt residues of the bulge. The inhibition of Pb²⁺-specific TAR cleavage by di- and trivalent metal ions revealed a binding specificity [in the order Co(NH₃)₆³⁺ > Mg²⁺ > Ca²⁺] at the bulge site. A BD simulation search of potential magnesium ion sites within the NMR structure of HIV-1 TAR RNA was conducted on a set of 20 conformers (PDB code 1ANR). For most cases, the bulge region was targeted by magnesium cations.     

Cross talk between the +73/294 interaction and the cleavage site in RNase P RNA mediated cleavage

To monitor functionally important metal ions and possible cross talk in RNase P RNA mediated cleavage we studied cleavage of substrates, where the 2′OH at the RNase P cleavage site (at −1) and/or at position +73 had been replaced with a 2′ amino group (or 2′H). Our data showed that the presence of 2′ modifications at these positions affected cleavage site recognition, ground state binding of substrate and/or rate of cleavage. Cleavage of 2′ amino substituted substrates at different pH showed that substitution of Mg²⁺ by Mn²⁺ (or Ca²⁺), identity of residues at and near the cleavage site, and addition of C5 protein influenced the frequency of miscleavage at −1 (cleavage at the correct site is referred to as +1). From this we infer that these findings point at effects mediated by protonation/deprotonation of the 2′ amino group, i.e. an altered charge distribution, at the site of cleavage. Moreover, our data suggested that the structural architecture of the interaction between the 3′ end of the substrate and RNase P RNA influence the charge distribution at the cleavage site as well as the rate of cleavage under conditions where the chemistry is suggested to be rate limiting. Thus, these data provide evidence for cross talk between the +73/294 interaction and the cleavage site in RNase P RNA mediated cleavage. We discuss the role metal ions might play in this cross talk and the likelihood that at least one functionally important metal ion is positioned in the vicinity of, and use the 2′OH at the cleavage site as an inner or outer sphere ligand.     

Collisionally activated dissociation and infrared multiphoton dissociation of oligonucleotides in a quadrupole ion trap

Infrared multiphoton dissociation (IRMPD) of deprotonated and protonated oligonucleotides ranging from 5 to 40 residues has been performed in a quadrupole ion trap mass spectrometer at normal operating pressure and temperature. Only moderate exposure times and laser powers were required to achieve efficient dissociation. In general, IRMPD and collisionally activated dissociation (CAD) produce comparable sequencing information, indicating that IRMPD is a viable alternative to CAD for oligonucleotide analysis in the quadrupole ion trap. Two major characteristics distinguish CAD and IRMPD spectra for a given parent ion. First, structurally uninformative M-B ions that dominate CAD spectra are generally only low-intensity species in IRMPD spectra because nonresonant activation causes these species to dissociate to backbone cleavage products. Second, phosphate and nucleobase ions can be observed directly in IRMPD experiments because the low-mass cutoff can be set to trap small fragment ions. For this reason IRMPD can sometimes facilitate analysis of sequences containing modified bases.     

Human DNA Polymerase ν Catalyzes Correct and Incorrect DNA Synthesis with High Catalytic Efficiency

DNA polymerase ν (pol ν) is a low fidelity A-family polymerase with a putative role in interstrand cross-link repair and homologous recombination. We carried out pre-steady-state kinetic analysis to elucidate the kinetic mechanism of this enzyme. We found that the mechanism consists of seven steps, similar that of other A-family polymerases. pol ν binds to DNA with a K_d for DNA of 9.2 nm, with an off-rate constant of 0.013 s⁻¹and an on-rate constant of 14 μm⁻¹ s⁻¹. dNTP binding is rapid with K_d values of 20 and 476 μm for the correct and incorrect dNTP, respectively. Pyrophosphorylation occurs with a K_d value for PP_i of 3.7 mm and a maximal rate constant of 11 s⁻¹. Pre-steady-state kinetics, examination of the elemental effect using dNTPαS, and pulse-chase experiments indicate that a rapid phosphodiester bond formation step is flanked by slow conformational changes for both correct and incorrect base pair formation. These experiments in combination with computer simulations indicate that the first conformational change occurs with rate constants of 75 and 20 s⁻¹; rapid phosphodiester bond formation occurs with a K_eq of 2.2 and 1.7, and the second conformational change occurs with rate constants of 2.1 and 0.5 s⁻¹, for correct and incorrect base pair formation, respectively. The presence of a mispair does not induce the polymerase to adopt a low catalytic conformation. pol ν catalyzes both correct and mispair formation with high catalytic efficiency.     

Structural basis for catalysis by onconase

Onconase® (ONC) is a homolog of bovine pancreatic ribonuclease (RNase A) from the frog Rana pipiens. ONC displays antitumoral activity and is in advanced clinical trials for the treatment of cancer. Here, we report the first atomic structures of ONC-nucleic acid complexes: a T89N/E91A ONC-5′-AMP complex at 1.65 Å resolution and a wild-type ONC-d(AUGA) complex at 1.90 Å resolution. The latter structure and site-directed mutagenesis were used to reveal the atomic basis for substrate recognition and turnover by ONC. The residues in ONC that are proximal to the scissile phosphodiester bond (His10, Lys31, and His97) and uracil nucleobase (Thr35, Asp67, and Phe98) are conserved from RNase A and serve to generate a similar bell-shaped pH versus k_cat/K_M profile for RNA cleavage. Glu91 of ONC forms two hydrogen bonds with the guanine nucleobase in d(AUGA), and Thr89 is in close proximity to that nucleobase. Installing a neutral or cationic residue at position 91 or an asparagine residue at position 89 virtually eliminated the 10²-fold guanine:adenine preference of ONC. A variant that combined such substitutions, T89N/E91A ONC, actually preferred adenine over guanine. In contrast, installing an arginine residue at position 91 increased the guanine preference and afforded an ONC variant with the highest known k_cat/K_M value. These data indicate that ONC discriminates between guanine and adenine by using Coulombic interactions and a network of hydrogen bonds. The structure of the ONC-d(AUGA) complex was also used to probe other aspects of catalysis. For example, the T5R substitution, designed to create a favorable Coulombic interaction between ONC and a phosphoryl group in RNA, increased ribonucleolytic activity by twofold. No variant, however, was more toxic to human cancer cells than wild-type ONC. Together, these findings provide a cynosure for understanding catalysis of RNA cleavage in a system of high medicinal relevance.     

Kinetics of tetramolecular quadruplexes

The melting of tetramolecular DNA or RNA quadruplexes is kinetically irreversible. However, rather than being a hindrance, this kinetic inertia allows us to study association and dissociation processes independently. From a kinetic point of view, the association reaction is fourth order in monomer and the dissociation first order in quadruplex. The association rate constant k_on, expressed in M⁻³·s⁻¹ decreases with increasing temperature, reflecting a negative activation energy (E_on) for the sequences presented here. Association is favored by an increase in monocation concentration. The first-order dissociation process is temperature dependent, with a very positive activation energy E_off, but nearly ionic strength independent. General rules may be drawn up for various DNA and RNA sequence motifs, involving 3–6 consecutive guanines and 0–5 protruding bases. RNA quadruplexes are more stable than their DNA counterparts as a result of both faster association and slower dissociation. In most cases, no dissociation is found for G-tracts of 5 guanines or more in sodium, 4 guanines or more in potassium. The data collected here allow us to predict the amount of time required for 50% (or 90%) quadruplex formation as a function of strand sequence and concentration, temperature and ionic strength.     

Active site constraints in the hydrolysis reaction catalyzed by bacterial RNase P: analysis of precursor tRNAs with a single 3'-S-phosphorothiolate internucleotide linkage

Endonucleolytic processing of precursor tRNAs (ptRNAs) by RNase P yields 3′-OH and 5′-phosphate termini, and at least two metal ions are thought to be essential for catalysis. To determine if the hydrolysis reaction catalyzed by bacterial RNase P (RNAs) involves stabilization of the 3′-oxyanion leaving group by direct coordination to one of the catalytic metal ions, ptRNA substrates with single 3′-S-phosphorothiolate linkages at the RNase P cleavage site were synthesized. With a 3′-S-phosphorothiolate-modified ptRNA carrying a 7 nt 5′-flank, a complete shift of the cleavage site to the next unmodified phosphodiester in the 5′-direction was observed. Cleavage at the modified linkage was not restored in the presence of thiophilic metal ions, such as Mn²⁺ or Cd²⁺. To suppress aberrant cleavage, we also constructed a 3′-S-phosphorothiolate-modified ptRNA with a 1 nt 5′-flank. No detectable cleavage of this substrate was seen in reactions catalyzed by RNase P RNAs from Escherichia coli and Bacillus subtilis, independent of the presence of thiophilic metal ions. Ground state binding of modified ptRNAs was not impaired, suggesting that the 3′-S-phosphorothiolate modification specifically prevents formation of the transition state, possibly by excluding catalytic metal ions from the active site.     

Redox and Chemical Activities of the Hemes in the Sulfur Oxidation Pathway Enzyme SoxAX

SoxAX enzymes couple disulfide bond formation to the reduction of cytochrome c in the first step of the phylogenetically widespread Sox microbial sulfur oxidation pathway. Rhodovulum sulfidophilum SoxAX contains three hemes. An electrochemical cell compatible with magnetic circular dichroism at near infrared wavelengths has been developed to resolve redox and chemical properties of the SoxAX hemes. In combination with potentiometric titrations monitored by electronic absorbance and EPR, this method defines midpoint potentials (E_m) at pH 7.0 of approximately +210, −340, and −400 mV for the His/Met, His/Cys⁻, and active site His/CysS⁻-ligated heme, respectively. Exposing SoxAX to S₂O₄²⁻, a substrate analog with E_m ∼−450 mV, but not Eu(II) complexed with diethylene triamine pentaacetic acid (E_m ∼−1140 mV), allows cyanide to displace the cysteine persulfide (CysS⁻) ligand to the active site heme. This provides the first evidence for the dissociation of CysS⁻ that has been proposed as a key event in SoxAX catalysis.     

Hydroxyl radical footprinting in vivo: mapping macromolecular structures with synchrotron radiation

We used a high flux synchrotron X-ray beam to map the structure of 16S rRNA and RNase P in viable bacteria in situ. A 300 ms exposure to the X-ray beam was sufficient for optimal cleavage of the phosphodiester backbone. The in vivo footprints of the 16S rRNA in frozen cells were similar to those obtained in vitro and were consistent with the predicted accessibility of the RNA backbone to hydroxyl radical. Protection or enhanced cleavage of certain nucleotides in vivo can be explained by interactions with tRNA and perturbation of the subunit interface. Thus, short exposures to a synchrotron X-ray beam can footprint the tertiary structure and protein contacts of RNA–protein complexes with nucleotide resolution in living cells.     

Essential Role of the CBD1-CBD2 Linker in Slow Dissociation of Ca2+ from the Regulatory Two-domain Tandem of NCX1

In NCX proteins CBD1 and CBD2 domains are connected through a short linker (3 or 4 amino acids) forming a regulatory tandem (CBD12). Only three of the six CBD12 Ca²⁺-binding sites contribute to NCX regulation. Two of them are located on CBD1 (K_d = ∼0.2 μm), and one is on CBD2 (K_d = ∼5 μm). Here we analyze how the intrinsic properties of individual regulatory sites are affected by linker-dependent interactions in CBD12 (AD splice variant). The three sites of CBD12 and CBD1 + CBD2 have comparable K_d values but differ dramatically in their Ca²⁺ dissociation kinetics. CBD12 exhibits multiphasic kinetics for the dissociation of three Ca²⁺ ions (k_r = 280 s⁻¹, k_f = 7 s⁻¹, and k_s = 0.4 s⁻¹), whereas the dissociation of two Ca²⁺ ions from CBD1 (k_f = 16 s⁻¹) and one Ca²⁺ ion from CBD2 (k_r = 125 s⁻¹) is monophasic. Insertion of seven alanines into the linker (CBD12–7Ala) abolishes slow dissociation of Ca²⁺, whereas the kinetic and equilibrium properties of three Ca²⁺ sites of CBD12–7Ala and CBD1 + CBD2 are similar. Therefore, the linker-dependent interactions in CBD12 decelerate the Ca²⁺ on/off kinetics at a specific CBD1 site by 50–80-fold, thereby representing Ca²⁺ “occlusion” at CBD12. Notably, the kinetic and equilibrium properties of the remaining two sites of CBD12 are “linker-independent,” so their intrinsic properties are preserved in CBD12. In conclusion, the dynamic properties of three sites are specifically modified, conserved, diversified, and integrated by the linker in CBD12, thereby generating a wide range dynamic sensor.     

Mechanistic Studies Reveal Similar Catalytic Strategies for Phosphodiester Bond Hydrolysis by Protein-only and RNA-dependent Ribonuclease P

Ribonuclease P (RNase P) is an endonuclease that catalyzes the essential removal of the 5′ end of tRNA precursors. Until recently, all identified RNase P enzymes were a ribonucleoprotein with a conserved catalytic RNA component. However, the discovery of protein-only RNase P (PRORP) shifted this paradigm, affording a unique opportunity to compare mechanistic strategies used by naturally evolved protein and RNA-based enzymes that catalyze the same reaction. Here we investigate the enzymatic mechanism of pre-tRNA hydrolysis catalyzed by the NYN (Nedd4-BP1, YacP nuclease) metallonuclease of Arabidopsis thaliana, PRORP1. Multiple and single turnover kinetic data support a mechanism where a step at or before chemistry is rate-limiting and provide a kinetic framework to interpret the results of metal alteration, mutations, and pH dependence. Catalytic activity has a cooperative dependence on the magnesium concentration (n_H = 2) under k_cat/K_m conditions, suggesting that PRORP1 catalysis is optimal with at least two active site metal ions, consistent with the crystal structure. Metal rescue of Asp-to-Ala mutations identified two aspartates important for enhancing metal ion affinity. The single turnover pH dependence of pre-tRNA cleavage revealed a single ionization (pK_a ∼ 8.7) important for catalysis, consistent with deprotonation of a metal-bound water nucleophile. The pH and metal dependence mirrors that observed for the RNA-based RNase P, suggesting similar catalytic mechanisms. Thus, despite different macromolecular composition, the RNA and protein-based RNase P act as dynamic scaffolds for the binding and positioning of magnesium ions to catalyze phosphodiester bond hydrolysis.     

Crystal structure of the leadzyme at 1.8 A resolution: metal ion binding and the implications for catalytic mechanism and allo site ion regulation

The leadzyme is a small ribozyme, derived from in vitro selection, which catalyzes site specific, Pb(2+)-dependent RNA cleavage. Pb(2+) is required for activity; Mg(2+) inhibits activity, while many divalent and trivalent ions enhance it. The leadzyme structure consists of an RNA duplex interrupted by a trinucleotide bulge. Here, crystal structures determined to 1.8 A resolution, both with Mg(2+) as the sole divalent counterion and with Mg(2+) and Sr(2+) (which mimics Pb(2+) with respect to binding but not catalysis), reveal the metal ion interactions with both the ground state and precatalytic conformations of the leadzyme. Mg(H(2)O)(6)(2+) ions bridge complementary strands of the duplex at multiple locations by binding tandem purines of one RNA strand in the major groove. At one site, Mg(H(2)O)(6)(2+) ligates the phosphodiester backbone of the trinucleotide bulge in the ground state conformation, but not in the precatalytic conformation, suggesting (a) Mg(2+) may inhibit leadzyme activity by stabilizing the ground state and (b) metal ions which displace Mg(2+) from this site may activate the leadzyme. Binding of Sr(2+) to the presumed catalytic Pb(2+) site in the precatalytic leadzyme induces local structural changes in a manner that would facilitate alignment of the catalytic ribose 2'-hydroxyl with the scissile bond for cleavage. These data support a model wherein binding of a catalytic ion to a precatalytic conformation of the leadzyme, in conjunction with the flexibility of the trinucleotide bulge, may facilitate structural rearrangements around the scissle phosphodiester bond favoring configurations that allow bond cleavage.     

Characterization of a Reduced Form of Plasma Plasminogen as the Precursor for Angiostatin Formation

Plasma plasminogen is the precursor of the tumor angiogenesis inhibitor, angiostatin. Generation of angiostatin in blood involves activation of plasminogen to the serine protease plasmin and facilitated cleavage of two disulfide bonds and up to three peptide bonds in the kringle 5 domain of the protein. The mechanism of reduction of the two allosteric disulfides has been explored in this study. Using thiol-alkylating agents, mass spectrometry, and an assay for angiostatin formation, we show that the Cys⁴⁶²-Cys⁵⁴¹ disulfide bond is already cleaved in a fraction of plasma plasminogen and that this reduced plasminogen is the precursor for angiostatin formation. From the crystal structure of plasminogen, we propose that plasmin ligands such as phosphoglycerate kinase induce a conformational change in reduced kringle 5 that leads to attack by the Cys⁵⁴¹ thiolate anion on the Cys⁵³⁶ sulfur atom of the Cys⁵¹²-Cys⁵³⁶ disulfide bond, resulting in reduction of the bond by thiol/disulfide exchange. Cleavage of the Cys⁵¹²-Cys⁵³⁶ allosteric disulfide allows further conformational change and exposure of the peptide backbone to proteolysis and angiostatin release. The Cys⁴⁶²-Cys⁵⁴¹ and Cys⁵¹²-Cys⁵³⁶ disulfides have −/+RHHook and −LHHook configurations, respectively, which are two of the 20 different measures of the geometry of a disulfide bond. Analysis of the structures of the known allosteric disulfide bonds identified six other bonds that have these configurations, and they share some functional similarities with the plasminogen disulfides. This suggests that the −/+RHHook and −LHHook disulfides, along with the −RHStaple bond, are potential allosteric configurations.     

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₁.