Pentavalent Organo-Vanadates as Transition State Analogues for Phosphoryl Transfer Reactions

Medicinal applications of vanadium complexes with Schiff bases

Many transition metal complexes have been explored for their therapeutic properties after the discovery of cisplatin. Schiff bases have an efficient complexation tendency with the transition metals and several medicinal properties have been reported. However, fewer studies have reported the medicinal utility of vanadium and its Schiff base complexes. This paper provides a comprehensive overview of vanadium complexes with Schiff bases along with their mechanistic insight. Vanadium complexes in + 4 and + 5 oxidation states have exhibited well-defined geometry and found to be thermodynamically stable. The studies have reported the G0/G1 phase cell cycle arrest and decreased delta psi m, inducing mitochondrial membrane depolarization in cancer cell lines along with the alterations in the metabolism of the cancer cells upon dosing with the vanadium complexes. Cancer cell invasion and growth are also found to be markedly reduced by peroxo complexes of vanadium. The studies included in the review paper have been taken from leading indexing databases and focus was laid on recent reports in literature. The biological potential of vanadium complexes of Schiff bases opens new horizons for future interdisciplinary studies and investigation focussed on understanding the biochemistry of these complexes, along with designing new complexes which have better bioavailability, solubility and low or non-toxicity.

Tipiracil binds to uridine site and inhibits Nsp15 endoribonuclease NendoU from SARS-CoV-2

SARS-CoV-2 Nsp15 is a uridine-specific endoribonuclease with C-terminal catalytic domain belonging to the EndoU family that is highly conserved in coronaviruses. As endoribonuclease activity seems to be responsible for the interference with the innate immune response, Nsp15 emerges as an attractive target for therapeutic intervention. Here we report the first structures with bound nucleotides and show how the enzyme specifically recognizes uridine moiety. In addition to a uridine site we present evidence for a second base binding site that can accommodate any base. The structure with a transition state analog, uridine vanadate, confirms interactions key to catalytic mechanisms. In the presence of manganese ions, the enzyme cleaves unpaired RNAs. This acquired knowledge was instrumental in identifying Tipiracil, an FDA approved drug that is used in the treatment of colorectal cancer, as a potential anti-COVID-19 drug. Using crystallography, biochemical, and whole-cell assays, we demonstrate that Tipiracil inhibits SARS-CoV-2 Nsp15 by interacting with the uridine binding pocket in the enzyme’s active site. Our findings provide new insights for the development of uracil scaffold-based drugs.

Youngchang Kim, Jacek Wower, and colleagues explore the sequence specificity, metal ion dependence and catalytic mechanism of the Nsp15 endoribonuclease NendoU from SARS-CoV-2. The authors also solve five new crystal structures of the enzyme in complex with 5’UMP, 3’UMP, 5’cGpU, uridine 2′,3′-vanadate (transition state analog) and Tipiracil (uracil mimic), and demonstrate that Tipiracil inhibits SARS-CoV-2 Nsp15 by interacting with the uridine binding pocket in the enzyme’s active site.

Tungstate as a Transition State Analog for Catalysis by Alkaline Phosphatase

The catalytic mechanisms underlying Escherichia coli alkaline phosphatase's (AP) remarkable rate enhancement have been probed extensively. Past work indicated that whereas the serine nucleophile (Ser102) electrostatically repels the product phosphate, another oxyanion, tungstate, binds more strongly in the presence of Ser102. These results predict a covalent bond between the serine nucleophile and tungstate, a model that we test herein. The crystal structure of tungstate-bound alkaline phosphatase provides evidence for a covalent adduct model and further shows that the ligand adopts trigonal bipyramidal geometry, which is infrequently observed for tungstate in small molecules and other active sites but mirrors the geometry of the presumed phosphoryl transfer transition state. The AP active site is known to stabilize another oxyanion, vanadate, in trigonal bipyramidal geometry, but the extent to which binding of either ligand reproduces the energetics of the transition state cannot be deduced from structural inspection alone. To test for transition state analog behavior, we determined the relationship between catalytic activity and affinity for tungstate and vanadate for a series of 20 AP variants. Affinity and activity were highly correlated for tungstate (r(2) = 0.89) but not vanadate (r(2) = 0.23), indicating that the tungstate•AP complex may better mimic this enzyme's transition state properties. The results herein suggest that tungstate will be a valuable tool for further dissecting AP catalysis and may prove helpful in mechanistic studies of other phosphoryl transfer enzymes.

QM/MM simulation (B3LYP) of the RNase A cleavage-transesterification reaction supports a triester A(N) + D(N) associative mechanism with an O2' H internal proton transfer

The mechanism of the backbone cleavage-transesterification step of the RNase A enzyme remains controversial even after 60 years of study. We report quantum mechanics/molecule mechanics (QM/MM) free energy calculations for two optimized reaction paths based on an analysis of all structural data and identified by a search for reaction coordinates using a reliable quantum chemistry method (B3LYP), equilibrated structural optimizations, and free energy estimations. Both paths are initiated by nucleophilic attack of the ribose O2' oxygen on the neighboring diester phosphate bond, and both reach the same product state (PS) (a O3'-O2' cyclic phosphate and a O5' hydroxyl terminated fragment). Path 1, resembles the widely accepted dianionic transition-state (TS) general acid (His119)/base (His12) classical mechanism. However, this path has a barrier (25 kcal/mol) higher than that of the rate-limiting hydrolysis step and a very loose TS. In Path 2, the proton initially coordinating the O2' migrates to the nonbridging O1P in the initial reaction path rather than directly to the general base resulting in a triester (substrate as base) AN + DN mechanism with a monoanionic weakly stable intermediate. The structures in the transition region are associative with low barriers (TS1 10, TS2 7.5 kcal/mol). The Path 2 mechanism is consistent with the many results from enzyme and buffer catalyzed and uncatalyzed analog reactions and leads to a PS consistent with the reactive state for the following hydrolysis step. The differences between the consistently estimated barriers in Path 1 and 2 lead to a 10(11) difference in rate strongly supporting the less accepted triester mechanism.

Fundamental challenges in mechanistic enzymology: progress toward understanding the rate enhancements of enzymes

Enzymes are remarkable catalysts that lie at the heart of biology, accelerating chemical reactions to an astounding extent with extraordinary specificity. Enormous progress in understanding the chemical basis of enzymatic transformations and the basic mechanisms underlying rate enhancements over the past decades is apparent. Nevertheless, it has been difficult to achieve a quantitative understanding of how the underlying mechanisms account for the energetics of catalysis, because of the complexity of enzyme systems and the absence of underlying energetic additivity. We review case studies from our own work that illustrate the power of precisely defined and clearly articulated questions when dealing with such complex and multifaceted systems, and we also use this approach to evaluate our current ability to design enzymes. We close by highlighting a series of questions that help frame some of what remains to be understood, and we encourage the reader to define additional questions and directions that will deepen and broaden our understanding of enzymes and their catalysis.

Metavanadate at the active site of the phosphatase VHZ

Vanadate is a potent modulator of a number of biological processes and has been shown by crystal structures and NMR spectroscopy to interact with numerous enzymes. Although these effects often occur under conditions where oligomeric forms dominate, the crystal structures and NMR data suggest that the inhibitory form is usually monomeric orthovanadate, a particularly good inhibitor of phosphatases because of its ability to form stable trigonal-bipyramidal complexes. We performed a computational analysis of a 1.14 Å structure of the phosphatase VHZ in complex with an unusual metavanadate species and compared it with two classical trigonal-bipyramidal vanadate-phosphatase complexes. The results support extensive delocalized bonding to the apical ligands in the classical structures. In contrast, in the VHZ metavanadate complex, the central, planar VO(3)(-) moiety has only one apical ligand, the nucleophilic Cys95, and a gap in electron density between V and S. A computational analysis showed that the V-S interaction is primarily ionic. A mechanism is proposed to explain the formation of metavanadate in the active site from a dimeric vanadate species that previous crystallographic evidence has shown to be able to bind to the active sites of phosphatases related to VHZ. Together, the results show that the interaction of vanadate with biological systems is not solely reliant upon the prior formation of a particular inhibitory form in solution. The catalytic properties of an enzyme may act upon the oligomeric forms primarily present in solution to generate species such as the metavanadate ion observed in the VHZ structure.

High-resolution analysis of Zn(2+) coordination in the alkaline phosphatase superfamily by EXAFS and x-ray crystallography

Comparisons among evolutionarily related enzymes offer opportunities to reveal how structural differences produce different catalytic activities. Two structurally related enzymes, Escherichia coli alkaline phosphatase (AP) and Xanthomonas axonopodis nucleotide pyrophosphatase/phosphodiesterase (NPP), have nearly identical binuclear Zn(2+) catalytic centers but show tremendous differential specificity for hydrolysis of phosphate monoesters or phosphate diesters. To determine if there are differences in Zn(2+) coordination in the two enzymes that might contribute to catalytic specificity, we analyzed both x-ray absorption spectroscopic and x-ray crystallographic data. We report a 1.29-Å crystal structure of AP with bound phosphate, allowing evaluation of interactions at the AP metal site with high resolution. To make systematic comparisons between AP and NPP, we measured zinc extended x-ray absorption fine structure for AP and NPP in the free-enzyme forms, with AMP and inorganic phosphate ground-state analogs and with vanadate transition-state analogs. These studies yielded average zinc-ligand distances in AP and NPP free-enzyme forms and ground-state analog forms that were identical within error, suggesting little difference in metal ion coordination among these forms. Upon binding of vanadate to both enzymes, small increases in average metal-ligand distances were observed, consistent with an increased coordination number. Slightly longer increases were observed in NPP relative to AP, which could arise from subtle rearrangements of the active site or differences in the geometry of the bound vanadyl species. Overall, the results suggest that the binuclear Zn(2+) catalytic site remains very similar between AP and NPP during the course of a reaction cycle.

Crystal structure of Onconase at 1.1 Å resolution--insights into substrate binding and collective motion

Onconase® (ONC) is an amphibian member of the pancreatic ribonuclease superfamily that is selectively toxic to tumor cells. It is a much less efficient enzyme than the archetypal ribonuclease A and, in an attempt to gain further insight, we report the first atomic resolution crystal structure of ONC, determined in complex with sulfate ions at 100 K. The electron density map is of a quality sufficient to reveal significant nonplanarity in several peptide bonds. The majority of active site residues are very well defined, with the exceptions being Lys31 from the catalytic triad and Lys33 from the B₁ subsite, which are relatively mobile but rigidify upon nucleotide binding. Cryocooling causes a compaction of the unit cell and the protein contained within. This is principally the result of an inward movement of one of the lobes of the enzyme (lobe 2), which also narrows the active site cleft. Binding a nucleotide in place of sulfate is associated with an approximately perpendicular movement of lobe 2 and has little further effect on the cleft width. Aspects of this deformation are present in the principal axes of anisotropy extracted from C^α atomic displacement parameters, indicating its intrinsic nature. The three lowest-frequency modes of ONC motion predicted by an anisotropic network model are compaction/expansion variations in which lobe 2 is the prime mover. Two of these have high similarity to the cryocooling response and imply that the essential ‘breathing’ motion of ribonuclease A is conserved in ONC. Instead, shifts in conformational equilibria may contribute to the reduced ribonucleolytic activity of ONC.

Prioritization of charge over geometry in transition state analogues of a dual specificity protein kinase

The direct observation of a transition state analogue (TSA) complex for tyrosine phosphorylation by a signaling kinase has been achieved using (19)F NMR analysis of MEK6 in complex with tetrafluoroaluminate (AlF(4)(-)), ADP, and p38α MAP kinase (acceptor residue: Tyr182). Solvent-induced isotope shifts and chemical shifts for the AlF(4)(-) moiety indicate that two fluorine atoms are coordinated by the two catalytic magnesium ions of the kinase active site, while the two remaining fluorides are liganded by protein residues only. An equivalent, yet distinct, AlF(4)(-) complex involving the alternative acceptor residue in p38α (Thr180) is only observed when the Tyr182 is mutated to phenylalanine. The formation of octahedral AlF(4)(-) species for both acceptor residues, rather than the trigonal bipyramidal AlF(3)(0) previously identified in the only other metal fluoride complex with a protein kinase, shows the requirement of MEK6 for a TSA that is isoelectronic with the migrating phosphoryl group. This requirement has hitherto only been demonstrated for proteins having a single catalytic magnesium ion.

o-Iminobenzosemiquinonate and o-imino-p-methylbenzosemiquinonate anion radicals coupled VO2+ stabilization

The diamagnetic VO(2+)-iminobenzosemiquinonate anion radical (L(R)(IS)(•-), R = H, Me) complexes, (L(-))(VO(2+))(L(R)(IS)(•-)): (L(1)(-))(VO(2+))(L(H)(IS)(•-))•3/2MeOH (1•3/2MeOH), (L(2)(-))(VO(2+))(L(H)(IS)(•-)) (2), and (L(2)(-))(VO(2+))(L(Me)(IS)(•-))•1/2 L(Me)(AP) (3•1/2 L(Me)(AP)), incorporating tridentate monoanionic NNO-donor ligands {L = L(1)(-) or L(2)(-), L(1)H = (2-[(phenylpyridin-2-yl-methylene)amino]phenol; L(2)H = 1-(2-pyridylazo)-2-naphthol; L(H)(IS)(•-) = o-iminobenzosemiquinonate anion radical; L(Me)(IS)(•-) = o-imino-p-methylbenzosemiquinonate anion radical; and L(Me)(AP) = o-amino-p-methylphenol} have been isolated and characterized by elemental analyses, IR, mass, NMR, and UV-vis spectra, including the single-crystal X-ray structure determinations of 1•3/2MeOH and 3•1/2 L(Me)(AP). Complexes 1•3/2MeOH, 2, and 3•1/2 L(Me)(AP) absorb strongly in the visible region because of intraligand (IL) and ligand-to-metal charge transfers (LMCT). 1•3/2MeOH is luminescent (λ(ext), 333 nm; λ(em), 522, 553 nm) in frozen dichloromethane-toluene glass at 77 K due to π(diimine→)π(diimine)* transition. The V-O(phenolato) (cis to the V═O) lengths, 1.940(2) and 1.984(2) Å, respectively, in 1•3/2MeOH and 3•1/2 L(Me)(AP) are consistent with the VO(2+) description. The V-O(iminosemiquinonate) (trans to the V═O) lengths, 2.1324(19) in 1•3/2MeOH and 2.083(2) Å in 3•1/2 L(Me)(AP), are expectedly ∼0.20 Å longer due to the trans influence of the V═O bond. Because of the stronger affinity of the paramagnetic VO(2+) ion to the L(H)(IS)(•-) or L(Me)(IS)(•-), the V-N(iminosemiquinonate) lengths, 1.908(2) and 1.921(2) Å, respectively, in 1•3/2MeOH and 3•1/2 L(Me)(AP), are unexpectedly shorter. Density functional theory (DFT) calculations using B3LYP, B3PW91, and PBE1PBE functionals on 1 and 2 have established that the closed shell singlet (CSS) solutions (VO(3+)-amidophenolato (L(R)(AP)(2-)) coordination) of these complexes are unstable with respect to triplet perturbations. But BS (1,1) M(s) = 0 (VO(2+)-iminobenzosemiquinonate anion radical (L(R)(IS)(•-)) coordination) solutions of these species are stable and reproduce the experimental bond parameters well. Spin density distributions of one electron oxidized cations are consistent with the [(L(-))(VO(2+))(L(R)(IQ))](+) descriptions [VO(2+)-o-iminobenzoquinone (L(R)(IQ)) coordination], and one electron reduced anions are consistent with the [(L(•2-))(VO(3+))(L(R)(AP)(2-))](-) descriptions [VO(3+)-amidophenolato (L(R)(AP)(2-)) coordination], incorporating the diimine anion radical (L(1)(•2-)) or azo anion radical (L(2)(3-)). Although, cations and anions are not isolable, but electro-and spectro-electrochemical experiments have shown that 3(+) and 3(-) ions are more stable than 1(+), 2(+) and 1(-), 2(-) ions. In all cases, the reductions occur with simultaneous two electron transfer, may be due to formation of coupled diimine/azo anion radical-VO(2+) species as in [(L(•2-))(VO(2+))(L(R)(AP)(2-))](2-).

Molecular dynamics simulation of bovine pancreatic ribonuclease A-CpA and transition state-like complexes

The mechanisms of enzymes are intimately connected with their overall structure and dynamics in solution. Experimentally, it is considerably challenging to provide detailed atomic level information about the conformational events that occur at different stages along the chemical reaction path. Here, theoretical tools may offer new potential insights that complement those obtained from experiments that may not yield an unambiguous mechanistic interpretation. In this study, we apply molecular dynamics simulations of bovine pancreatic ribonuclease A, an archetype ribonuclease, to study the conformational dynamics, structural relaxation, and differential solvation that occur at discrete stages of the transesterification and cleavage reaction. Simulations were performed with explicit solvation with rigorous electrostatics and utilize recently developed molecular mechanical force field parameters for transphosphorylation and hydrolysis transition state analogues. Herein, we present results for the enzyme complexed with the dinucleotide substrate cytidilyl-3',5'-adenosine (CpA) in the reactant, and transphosphorylation and hydrolysis transition states. A detailed analysis of active site structures and hydrogen-bond patterns is presented and compared. The integrity of the overall backbone structure is preserved in the simulations and supports a mechanism whereby His12 stabilizes accumulating negative charge at the transition states through hydrogen-bond donation to the nonbridge oxygens. Lys41 is shown to be highly versatile along the reaction coordinate and can aid in the stabilization of the dianionic transition state, while being poised to act as a general acid catalyst in the hydrolysis step.

Transition state analogue structures of human phosphoglycerate kinase establish the importance of charge balance in catalysis

Transition state analogue (TSA) complexes formed by phosphoglycerate kinase (PGK) have been used to test the hypothesis that balancing of charge within the transition state dominates enzyme-catalyzed phosphoryl transfer. High-resolution structures of trifluoromagnesate (MgF(3)(-)) and tetrafluoroaluminate (AlF(4)(-)) complexes of PGK have been determined using X-ray crystallography and (19)F-based NMR methods, revealing the nature of the catalytically relevant state of this archetypal metabolic kinase. Importantly, the side chain of K219, which coordinates the alpha-phosphate group in previous ground state structures, is sequestered into coordinating the metal fluoride, thereby creating a charge environment complementary to the transferring phosphoryl group. In line with the dominance of charge balance in transition state organization, the substitution K219A induces a corresponding reduction in charge in the bound aluminum fluoride species, which changes to a trifluoroaluminate (AlF(3)(0)) complex. The AlF(3)(0) moiety retains the octahedral geometry observed within AlF(4)(-) TSA complexes, which endorses the proposal that some of the widely reported trigonal AlF(3)(0) complexes of phosphoryl transfer enzymes may have been misassigned and in reality contain MgF(3)(-).

Kinetic analysis of beta-phosphoglucomutase and its inhibition by magnesium fluoride

The isomerization of beta-glucose-1-phosphate (betaG1P) to beta-glucose-6-phosphate (G6P) catalyzed by beta-phosphoglucomutase (betaPGM) has been examined using steady- and presteady-state kinetic analysis. In the presence of low concentrations of beta-glucose-1,6-bisphosphate (betaG16BP), the reaction proceeds through a Ping Pong Bi Bi mechanism with substrate inhibition (kcat = 65 s(-1), K(betaG1P) = 15 microM, K(betaG16BP) = 0.7 microM, Ki = 122 microM). If alphaG16BP is used as a cofactor, more complex kinetic behavior is observed, but the nonlinear progress curves can be fit to reveal further catalytic parameters (kcat = 74 s(-1), K(betaG1P) = 15 microM, K(betaG16BP) = 0.8 microM, Ki = 122 microM, K(alphaG16BP) = 91 microM for productive binding, K(alphaG16BP) = 21 microM for unproductive binding). These data reveal that variations in the substrate structure affect transition-state affinity (approximately 140,000-fold in terms of rate acceleration) substantially more than ground-state binding (110-fold in terms of binding affinity). When fluoride and magnesium ions are present, time-dependent inhibition of the betaPGM is observed. The concentration dependence of the parameters obtained from fitting these progress curves shows that a betaG1P x MgF3(-) x betaPGM inhibitory complex is formed under the reaction conditions. The overall stability constant for this complex is approximately 2 x 10(-16) M(5) and suggests an affinity of the MgF3(-) moiety to this transition-state analogue (TSA) of < or = 70 nM. The detailed kinetic analysis shows how a special type of TSA that does not exist in solution is assembled in the active site of an enzyme. Further experiments show that under the conditions of previous structural studies, phosphorylated glucose only persists when bound to the enzyme as the TSA. The preference for TSA formation when fluoride is present, and the hydrolysis of substrates when it is not, rules out the formation of a stable pentavalent phosphorane intermediate in the active site of betaPGM.

Transition state analogues for nucleotidyl transfer reactions: Structure and stability of pentavalent vanadate and phosphate ester dianions

The structures and energy of phosphate dimethyl ester and vanadate dimethyl ester have been calculated using B3LYP/TZVP density functional quantum chemical methods and polarized continuum (PCM) and Langevin dipoles solvation models. These calculations were carried out to obtain fundamental information on the ability of vanadate esters to function as transition state analogues for the nucleotidyl transfer reaction catalyzed by DNA polymerases. Base-catalyzed methanolysis of the phosphate and vanadate dimethyl esters were the model reactions examined in this study. The structures of the phosphate and vanadate dimethyl esters and pentavalent intermediates in aqueous solution were optimized and evaluated at the PCM/B3LYP/TZVP level. The three-dimensional free energy surfaces for the studied reactions were determined at the PCM/B3LYP/TZVP//B3LYP/TZVP level. Comparison with experimental structural data obtained from the Cambridge Structural Database and with the observed kinetics of phosphate diester hydrolysis demonstrated that the level of theory chosen for these studies was appropriate. The results showed that structurally and electrostatically the vanadate dimethylester and a five-coordinate nearly trigonal bipyramidal intermediate were reasonable analogues for the parent phosphorus systems. Despite these similarities in structure, the energetics of the two systems were different, and the transition states of the two model reactions were found on different areas of the potential energy surface. When the binding energy of a transition state-DNA polymerase complex was extrapolated to a transition state analogue-DNA polymerase complex, the formation of a simple dianionic pentavalent vanadate ester adduct in the enzyme active site was not found to be sufficiently favorable. This finding suggests that additional stabilization of this adduct is needed before this type of transition state analogue will be likely to yield stable adducts with this class of enzymes. New possible candidates for such complexes are suggested.

Enzymatic and Structural Characterisation of Amphinase, a Novel Cytotoxic Ribonuclease from Rana pipiens Oocytes

Besides Onconase (ONC) and its V11/N20/R103-variant, oocytes of the Northern Leopard frog (Rana pipiens) contain another homologue of ribonuclease A, which we named Amphinase (Amph). Four variants (Amph-1-4) were isolated and sequenced, each 114 amino acid residues in length and N-glycosylated at two positions. Sequence identities (a) among the variants and (b) versus ONC are 86.8-99.1% and 38.2-40.0%, respectively. When compared with other amphibian ribonucleases, a typical pattern of cysteine residues is evident but the N-terminal pyroglutamate residue is replaced by a six-residue extension. Amph variants have relatively weak ribonucleolytic activity that is insensitive to human ribonuclease inhibitor protein (RI). Values of k(cat)/K(M) with hypersensitive fluorogenic substrates are 10(4) and 10(2)-fold lower than the maximum values exhibited by ribonuclease A and ONC, respectively, and there is little cytosine/uracil or adenine/guanine discrimination at the B(1) or B(2) subsites, respectively. Amph variants have cytotoxic activity toward A-253 carcinoma cells that requires intact ribonucleolytic activity. The glycan component has little or no influence over single-stranded RNA cleavage, RI evasion or cytotoxicity. The crystal structures of natural and recombinant Amph-2 (determined at 1.8 and 1.9 A resolution, respectively) reveal that the N terminus is unlikely to play a catalytic role (but an unusual alpha2-beta1 loop may do so) and the B(2) subsite is rudimentary. At the active site, structural features that may contribute to the enzyme's low ribonucleolytic activity are the fixture of Lys14 in an obstructive position, the accompanying ejection of Lys42, and a lack of constraints on the conformations of Lys42 and His107.

A comparison of vanadate to a 2'-5' linkage at the active site of a small ribozyme suggests a role for water in transition-state stabilization

The potential for water to participate in RNA catalyzed reactions has been the topic of several recent studies. Here, we report crystals of a minimal, hinged hairpin ribozyme in complex with the transition-state analog vanadate at 2.05 A resolution. Waters are present in the active site and are discussed in light of existing views of catalytic strategies employed by the hairpin ribozyme. A second structure harboring a 2',5'-phosphodiester linkage at the site of cleavage was also solved at 2.35 A resolution and corroborates the assignment of active site waters in the structure containing vanadate. A comparison of the two structures reveals that the 2',5' structure adopts a conformation that resembles the reaction intermediate in terms of (1) the positioning of its nonbridging oxygens and (2) the covalent attachment of the 2'-O nucleophile with the scissile G+1 phosphorus. The 2',5'-linked structure was then overlaid with scissile bonds of other small ribozymes including the glmS metabolite-sensing riboswitch and the hammerhead ribozyme, and suggests the potential of the 2',5' linkage to elicit a reaction-intermediate conformation without the need to form metalloenzyme complexes. The hairpin ribozyme structures presented here also suggest how water molecules bound at each of the nonbridging oxygens of G+1 may electrostatically stabilize the transition state in a manner that supplements nucleobase functional groups. Such coordination has not been reported for small ribozymes, but is consistent with the structures of protein enzymes. Overall, this work establishes significant parallels between the RNA and protein enzyme worlds.

Crystal structure of tryptophanyl-tRNA synthetase complexed with adenosine-5' tetraphosphate: evidence for distributed use of catalytic binding energy in amino acid activation by class I aminoacyl-tRNA synthetases

Tryptophanyl-tRNA synthetase (TrpRS) is a functionally dimeric ligase, which specifically couples hydrolysis of ATP to AMP and pyrophosphate to the formation of an ester bond between tryptophan and the cognate tRNA. TrpRS from Bacillus stearothermophilus binds the ATP analogue, adenosine-5' tetraphosphate (AQP) competitively with ATP during pyrophosphate exchange. Estimates of binding affinity from this competitive inhibition and from isothermal titration calorimetry show that AQP binds 200 times more tightly than ATP both under conditions of induced-fit, where binding is coupled to an unfavorable conformational change, and under exchange conditions, where there is no conformational change. These binding data provide an indirect experimental measurement of +3.0 kcal/mol for the conformational free energy change associated with induced-fit assembly of the active site. Thermodynamic parameters derived from the calorimetry reveal very modest enthalpic changes, consistent with binding driven largely by a favorable entropy change. The 2.5 A structure of the TrpRS:AQP complex, determined de novo by X-ray crystallography, resembles that of the previously described, pre-transition state TrpRS:ATP complexes. The anticodon-binding domain untwists relative to the Rossmann-fold domain by 20% of the way toward the orientation observed for the Products complex. An unexpected tetraphosphate conformation allows the gamma and deltad phosphate groups to occupy positions equivalent to those occupied by the beta and gamma phosphates of ATP. The beta-phosphate effects a 1.11 A extension that relocates the alpha-phosphate toward the tryptophan carboxylate while the PPi mimic moves deeper into the KMSKS loop. This configuration improves interactions between enzyme and nucleotide significantly and uniformly in the adenosine and PPi binding subsites. A new hydrogen bond forms between S194 from the class I KMSKS signature sequence and the PPi mimic. These complementary thermodynamic and structural data are all consistent with the conclusion that the tetraphosphate mimics a transition-state in which the KMSKS loop develops increasingly tight bonds to the PPi leaving group, weakening linkage to the Palpha as it is relocated by an energetically favorable domain movement. Consistent with extensive mutational data on Tyrosyl-tRNA synthetase, this aspect of the mechanism develops high transition-state affinity for the adenosine and pyrophosphate moieties, which move significantly, relative to one another, during the catalytic step.

Genetic selection for critical residues in ribonucleases

Homologous mammalian proteins were subjected to an exhaustive search for residues that are critical to their structure/function. Error-prone polymerase chain reactions were used to generate random mutations in the genes of bovine pancreatic ribonuclease (RNase A) and human angiogenin, and a genetic selection based on the intrinsic cytotoxicity of ribonucleolytic activity was used to isolate inactive variants. Twenty-three of the 124 residues in RNase A were found to be intolerant to substitution with at least one particular amino acid. Twenty-nine of the 123 residues in angiogenin were likewise intolerant. In both RNase A and angiogenin, only six residues appeared to be wholly intolerant to substitution: two histidine residues involved in general acid/base catalysis and four cysteine residues that form two disulfide bonds. With few exceptions, the remaining critical residues were buried in the hydrophobic core of the proteins. Most of these residues were found to tolerate only conservative substitutions. The importance of a particular residue as revealed by this genetic selection correlated with its sequence conservation, though several non-conserved residues were found to be critical for protein structure/function. Despite voluminous research on RNase A, the importance of many residues identified herein was unknown, and those can now serve as targets for future work. Moreover, a comparison of the critical residues in RNase A and human angiogenin, which share only 35% amino acid sequence identity, provides a unique perspective on the molecular evolution of the RNase A superfamily, as well as an impetus for applying this methodology to other ribonucleases.

Syntheses of Silatranyl- and Germatranyluridines

[reaction: see text] Silatranyluridine 1 and germatranyluridine 2 have been prepared in five steps from oxazolinouridine 3 in 27 and 29% yields, respectively. These compounds are novel transition-state analogues (TSAs) for RNA hydrolysis and offer a number of advantages over traditional vanadium- or rhenium-based TSAs. Germatrane 2 is completely stable in D(2)O at room temperature, and the half-life of silatrane 1 in D(2)O was found to be >7 days by (1)H NMR spectroscopy.

The power of vanadate in crystallographic investigations of phosphoryl transfer enzymes

The formation of transition state mimics of phosphoryl transfer reactions with the metal oxoanion vanadate is a powerful technique in macromolecular crystallography. The tendency of vanadate to form pentacovalent complexes exhibiting trigonal bipyramidal geometry makes this compound a close approximation of the transition state for such reactions. In many cases, vanadate complexes provide the most accurate visualization of the transition state that can be reasonably achieved. A survey of the Protein Data Bank reveals that a relatively small number of structures (39, representing 23 unique proteins) include vanadate, yet these structures represent four of the six E.C. categories of enzymes, and were obtained in crystals with pH values ranging from 5.0 to 7.8. Vanadate has additional advantages over other compounds such as aluminum fluoride, beryllium fluoride and nitrate used for visualization of transition state mimics in that vanadate readily forms covalent bonds with a variety of ligands and has produced a wider variety of transition state mimics. Given the hundreds of crystal structures that have been solved for phosphoryl transfer enzymes, it is surprising that vanadate has not been used more frequently for visualization of transition state analogs. We propose that an opportunity exists for vanadate to become a more commonly utilized component of the macromolecular crystallographer's toolbox.

Contribution of Active-Site Residues to the Function of Onconase, a Ribonuclease with Antitumoral Activity

Onconase (ONC), a homologue of ribonuclease A (RNase A), is in clinical trials for the treatment of cancer. ONC possesses a conserved active-site catalytic triad, which is composed of His10, Lys31, and His97. The three-dimensional structure of ONC suggests that two additional residues, Lys9 and an N-terminal lactam formed from a glutamine residue (Pca1), could also contribute to catalysis. To determine the role of Pca1, Lys9, and Lys31 in the function of ONC, site-directed mutagenesis was used to replace each with alanine. Values of k(cat)/K(M) for the variants were determined with a novel fluorogenic substrate, which was designed to match the nucleobase specificity of ONC and gives the highest known k(cat)/K(M) value for the enzyme. The K9A and K31A variants display 10(3)-fold lower k(cat)/K(M) values than the wild-type enzyme, and a K9A/K31A double variant suffers a >10(4)-fold decrease in catalytic activity. In addition, replacing Lys9 or Lys31 eliminates the antitumoral activity of ONC. The side chains of Pca1 and Lys9 form a hydrogen bond in crystalline ONC. Replacing Pca1 with an alanine residue lowers the catalytic activity of ONC by 20-fold. Yet, replacing Pca1 in the K9A variant enzyme does not further reduce catalytic activity, revealing that the function of the N-terminal pyroglutamate residue is to secure Lys9. The thermodynamic cycle derived from k(cat)/K(M) values indicates that the Pca1...Lys9 hydrogen bond contributes 2.0 kcal/mol to the stabilization of the rate-limiting transition state during catalysis. Finally, binding isotherms with a substrate analogue indicate that Lys9 and Lys31 contribute little to substrate binding and that the low intrinsic catalytic activity of ONC originates largely from the low affinity of the enzyme for its substrate. These findings could assist the further development of ONC as a cancer chemotherapeutic.

Temperature dependence of the backbone dynamics of ribonuclease A in the ground state and bound to the inhibitor 5'-phosphothymidine (3'-5')pyrophosphate adenosine 3'-phosphate

The interaction of the dinucleotide inhibitor 5'-phosphothymidine(3',5')pyrophosphate adenosine 3'-phosphate (pTppAp) with bovine pancreatic ribonuclease A (RNase A) was characterized by calorimetry and solution NMR spectroscopy. Calorimetric data show that binding of pTppAp to RNase A is exothermic (DeltaH = -60.1 +/- 4.1 kJ/mol) with a dissociation constant of 16 nM at 298 K. At this temperature, the binding results in an entropy loss (TDeltaS = -16.8 +/- 7.3 kJ/mol) that is more favorable than that with the product analogue, 2'-CMP (TDeltaS = -31.3 +/- 0.9 kJ/mol). Temperature-dependent calorimetric experiments give a DeltaC(p) for ligand binding of -230 +/- 100 J/mol K. Binding of pTppAp results in noticeable effects on the backbone amide chemical shifts and dynamics. Amide backbone (15)N NMR spin-relaxation studies were performed on both apo RNase A and RNase A/pTppAp as a function of temperature. At each temperature, the model-free-determined order parameters, S(2), were significantly higher for RNase A/pTppAp than for the apo enzyme indicating a decrease in the conformational entropy of the protein upon ligand binding. Furthermore, the magnitude of this difference varies along the amino acid sequence specifically locating the entropic changes. The temperature dependence of S(2) at each residue enabled assessment of the local heat capacity changes (DeltaC(p)) from ligand binding. In an overall, average sense, DeltaC(p) for the protein backbone, determined from the NMR dynamics measurements, did not differ between apo RNase A and RNase A/pTppAp indicating that backbone dynamics contribute little to DeltaC(p) for protein-ligand interactions in this system. However, residue-by-residue comparison of the temperature-dependent change in entropy (DeltaS(B)) between free and bound forms reveals nonzero contributions to DeltaC(p) at individual sites. The balance of positive and negative changes reveals a redistribution of energetics upon binding. Furthermore, experiment and semiempirical estimates suggest that a large negative DeltaC(p) should accompany binding of pTppAp, and we conclude that this contribution must arise from factors other than amide backbone dynamics.

Transition state stabilization by a catalytic RNA

The hairpin ribozyme catalyzes sequence-specific cleavage of RNA through transesterification of the scissile phosphate. Vanadate has previously been used as a transition state mimic of protein enzymes that catalyze the same reaction. Comparison of the 2.2 angstrom resolution structure of a vanadate-hairpin ribozyme complex with structures of precursor and product complexes reveals a rigid active site that makes more hydrogen bonds to the transition state than to the precursor or product. Because of the paucity of RNA functional groups capable of general acid-base or electrostatic catalysis, transition state stabilization is likely to be an important catalytic strategy for ribozymes.

Is the PTPase-vanadate complex a true transition state analogue?

Vanadate can often bind to phosphoryl transfer enzymes to form a trigonal-bipyramidal structure at the active site. The enzyme-vanadate dissociation constants in these enzymes are much lower than those for phosphate. Therefore, enzyme-bound vanadate moieties are often considered as transition state analogues. To test whether the enzyme-vanadate complex is a true transition state analogue beyond the simple geometry and binding affinity arguments and whether the bond orders of the VO bonds in the complex approach those of the PO bonds in the transition state, the binding properties of vanadate in the Yersinia protein-tyrosine phosphatase (PTPase) and its T410A, D356N, W354A, R409K, and D356A mutants have been studied by steady-state kinetic measurements and by difference Raman measurements. The results of the kinetic measurements show no correlation between K(I) and kcat or kcat/K(m) in these mutants. In addition, our analysis of the Raman data shows that the bond order change of the nonbridging V--O bonds in the vanadate complexes does not correlate with the kinetic parameters in a number of PTPase variants as predicted by the transition state binding paradigm. Furthermore, the ionization state of the bound vanadate moiety is not invariant across the PTPase variants studied, and the average bond order of the nonbridging V--O bonds decreased by 0.06-0.07 valence unit in the wild type and all of the mutant PTPases, either in dianionic or in monoanionic form. Thus the complex would resemble an associative transition state, contrary to the previously determined dissociative structure of the transition state. Therefore, it is concluded that vanadate is not a true transition state analogue for the PTPase reactions.

Mechanism of inhibition of the beta-lactamase of Enterobacter cloacae P99 by 1:1 complexes of vanadate with hydroxamic acids

The class C beta-lactamase of Enterobacter cloacae P99 is competitively inhibited by low concentrations of 1:1 complexes of vanadate and hydroxamic acids. Structure-activity studies indicated that the hydroxamic acid functional group was essential to this inhibition. Both aryl and alkyl hydroxamic acids form inhibitory ternary complexes with vanadate and the enzyme, although, in certain cases of the latter, the inhibition may not be seen because of the low formation constants of the vanadate-hydroxamic acid complex. After all of the vanadate species present in solution had been taken into account, "real" K(i) values for the vanadate complexes could be determined. The K(i) value of the best of the inhibitors that were investigated, the 1:1 complex of vanadate with 4-nitrobenzohydroxamic acid, was 0.48 microM. Kinetics studies showed that the association and dissociation rate constants of this complex with the enzyme were 1.48 x 10(6) s(-1) M(-1) and 0.73 s(-1), respectively; the magnitude of the latter indicates covalent interaction of the complex with the enzyme. (51)V NMR and UV-vis spectra suggest that the structure of the vanadate complex bound to the enzyme may be very similar to that in solution. A (13)C NMR spectrum of the enzyme complex with 4-nitrobenzo[(13)C]hydroxamic acid and vanadate yields a coordination-induced shift (CIS) of 7.74 ppm. This is significantly larger than that of the vanadate complex in free solution (3.62 ppm), suggesting either, somewhat contrary to the (51)V and UV-vis spectra, greater interaction between vanadium and the hydroxamate carbonyl oxygen in the enzyme complex than in free solution or, more likely, polarization of the hydroxamate by interaction, e.g., hydrogen bonding, with the enzyme. Molecular modeling indicates that a pentacoordinated vanadate complex may well be able to snugly occupy the enzyme active site; Asn 152 is suitably placed to hydrogen bond to the hydroxamic acid oxygen atom. The experimental results are in accord with a model whereby the vanadate-hydroxamate-enzyme complex is a moderately good analogue of the transition state of the reaction of the beta-lactamase with phosphonate inhibitors.

Quantitative analysis of the effect of salt concentration on enzymatic catalysis

Like pH, salt concentration can have a dramatic effect on enzymatic catalysis. Here, a general equation is derived for the quantitative analysis of salt-rate profiles: k(cat)/K(M) = (k(cat)/K(M))(MAX)/[1+([Na+]/K[Na+])(n')], where (k(cat)/K(M))(MAX) is the physical limit of k(cat)/K(M), K(Na+) is the salt concentration at which k(cat)/K(M) = (k(cat)/K(M))(MAX)/2, and -n' is the slope of the linear region in a plot of log(k(cat)/K(M)) versus log [Na+]. The value of n' is of special utility, as it reflects the contribution of Coulombic interactions to the uniform binding of the bound states. This equation was used to analyze salt effects on catalysis by ribonuclease A (RNase A), which is a cationic enzyme that catalyzes the cleavage of an anionic substrate, RNA, with k(cat)/K(M) values that can exceed 10(9) M(-1) s(-1). Lys7, Arg10, and Lys66 comprise enzymic subsites that are remote from the active site. Replacing Lys7, Arg10, and Lys66 with alanine decreases the charge on the enzyme as well as the value of n'. Likewise, decreasing the number of phosphoryl groups in the substrate decreases the value of n'. Replacing Lys41, a key active-site residue, with arginine creates a catalyst that is limited by the chemical conversion of substrate to product. This change increases the value of n', as expected for a catalyst that is more sensitive to changes in the binding of the chemical transition state. Hence, the quantitative analysis of salt-rate profiles can provide valuable insight into the role of Coulombic interactions in enzymatic catalysis.

Contribution of the active site histidine residues of ribonuclease A to nucleic acid binding

His12 and His119 are critical for catalysis of RNA cleavage by ribonuclease A (RNase A). Substitution of either residue with an alanine decreases the value of k(cat)/K(M) by more than 10(4)-fold. His12 and His119 are proximal to the scissile phosphoryl group of an RNA substrate in enzyme-substrate complexes. Here, the role of these active site histidines in RNA binding was investigated by monitoring the effect of mutagenesis and pH on the stability of enzyme-nucleic acid complexes. X-ray diffraction analysis of the H12A and H119A variants at a resolution of 1.7 and 1.8 A, respectively, shows that the amino acid substitutions do not perturb the overall structure of the variants. Isothermal titration calorimetric studies on the complexation of wild-type RNase A and the variants with 3'-UMP at pH 6.0 show that His12 and His119 contribute 1.4 and 1.1 kcal/mol to complex stability, respectively. Determination of the stability of the complex of wild-type RNase A and 6-carboxyfluorescein approximately d(AUAA) at varying pHs by fluorescence anisotropy shows that the stability increases by 2.4 kcal/mol as the pH decreases from 8.0 to 4.0. At pH 4.0, replacing His12 with an alanine residue decreases the stability of the complex with 6-carboxyfluorescein approximately d(AUAA) by 2.3 kcal/mol. Together, these structural and thermodynamic data provide the first thorough analysis of the contribution of histidine residues to nucleic acid binding.

Mechanism of reaction of acyl phosph(on)ates with the beta-lactamase of Enterobacter cloacae P99

A series of acyl phosph(on)ates has been prepared to more closely examine the details of the interactions of this class of molecule with beta-lactamases. In general, they were found to react with the class C beta-lactamase of Enterobacter cloacae P99 in two ways, by acylation and by phosphylation. The acyl-enzymes generated by the former reaction were transiently stable with half-lives of between 3 and 45 s, of comparable lifetime therefore to those generated by the inhibitory beta-lactams cefotaxime, cefuroxime, and cefoxitin. On the other hand, phosphylation led to a completely inactive enzyme. In general, the second-order rate constants for acylation (k(cat)/K(m)) were larger than for phosphylation (k(i)). As expected on chemical grounds, phosphylation was found to be relatively more effective for the phosphonates than the phosphates. The acyl phosphates were much more effective acylating agents however. The acylation reaction was found to be enhanced by hydrophobic substituents in both the acyl and leaving group moieties. Thus, the most reactive compound in this series was benzo[b]thiophene-2-carbonyl 2'-naphthyl phosphate with a K(m) value of 0.15 microM and a k(cat) of 0.2 s(-1); k(cat)/K(m) is therefore 1.3 x 10(6) s(-1) M(-1), making this compound the most specific acyclic acylation reagent for this beta-lactamase yet described. Significant substrate inhibition by this compound suggested that further binding regions may be available for exploitation in inhibitor design. A linear free energy analysis showed that the transition states for acylation of the beta-lactamase by aroyl phosphates are analogues of the corresponding aryl boronic acid adducts. Molecular modeling suggested that the aroyl phosphates react with the P99 beta-lactamase with the aroyl group in the side chain/acyl group site of normal substrates and the phosphate in the leaving group site. In this orientation, the phosphate leaving group interacts strongly with Lys 315.

Bridging the gap between proteins and nucleic acids: a metal-independent RNAseA mimic with two protein-like functionalities

Two synthetically modified nucleoside triphosphate analogues (adenosine modified with an imidazole and uridine modified with a cationic amine) are enzymatically polymerized in tandem along a degenerate DNA library for the combinatorial selection of an RNAse A mimic. The selected activity is consistent with both electrostatic and general acid/base catalysis at physiological pH in the absence of divalent metal cations. The simultaneous use of two modified nucleotides to enrich the catalytic repertoire of DNA-based catalysts has never before been demonstrated and evidence of general acid/base catalysis at pH 7.4 for a DNAzyme has never been previously observed in the absence of a divalent metal cation or added cofactor. This work illustrates how the incorporation of protein-like functionalities in nucleic acids can bridge the gap between proteins and oligonucleotides underscoring the potential for using nucleic acid scaffolds in the development of new materials and improved catalysts for use in chemistry and medicine.