Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.     

Magnesium dependence of the amplified conformational switch in the trans-acting hepatitis delta virus ribozyme

The ability of divalent metal ions to participate in both structure formation and catalytic chemistry of RNA enzymes (ribozymes) has made it difficult to separate their cause and effect in ribozyme function. For example, the recently solved crystal structures of precursor and product forms of the cis-cleaving genomic hepatitis delta virus (HDV) ribozyme show a divalent metal ion bound in the active site that is released upon catalysis due to an RNA conformational change. This conformational switch is associated with a repositioning of the catalytically involved base C75 in the active-site cleft, thus controlling catalysis. These findings confirm previous data from fluorescence resonance energy transfer (FRET) on a trans-acting form of the HDV ribozyme that found a global conformational change to accompany catalysis. Here, we further test the conformational switch model by measuring the Mg(2+) dependence of the global conformational change of the trans-acting HDV ribozyme, using circular dichroism and time-resolved FRET as complementary probes of secondary and tertiary structure formation, respectively. We observe significant differences in both structure and Mg(2+) affinity of the precursor and product forms, in the presence and absence of 300 mM Na(+) background. The precursor shortens while the product extends with increasing Mg(2+) concentration, essentially amplifying the structural differences observed in the crystal structures. In addition, the precursor has an approximately 2-fold and approximately 13-fold lower Mg(2+) affinity than the product in secondary and tertiary structure formation, respectively. We also have compared the C75 wild-type with the catalytically inactive C75U mutant and find significant differences in global structure and Mg(2+) affinity for both their precursor and product forms. Significantly, the Mg(2+) affinity of the C75 wild-type is 1.7-2.1-fold lower than that of the C75U mutant, in accord with the notion that C75 is essential for a catalytic conformational change that leads to a decrease in the local divalent metal ion affinity and release of a catalytic metal. Thus, a consistent picture emerges in which divalent metal ions and RNA functional groups are intimately intertwined in affecting structural dynamics and catalysis in the HDV ribozyme.     

Folding of the hammerhead ribozyme: pyrrolo-cytosine fluorescence separates core folding from global folding and reveals a pH-dependent conformational change

The catalytic activity of the hammerhead ribozyme is limited by its ability to fold into the native tertiary structure. Analysis of folding has been hampered by a lack of assays that can independently monitor the environment of nucleobases throughout the ribozyme-substrate complex in real time. Here, we report the development and application of a new folding assay in which we use pyrrolo-cytosine (pyC) fluorescence to (1) probe active-site formation, (2) examine the ability of peripheral ribozyme domains to support native folding, (3) identify a pH-dependent conformational change within the ribozyme, and (4) explore its influence on the equilibrium between the folded and unfolded core of the hammerhead ribozyme. We conclude that the natural ribozyme folds in two distinct noncooperative steps and the pH-dependent correlation between core folding and activity is linked to formation of the G8-C3 base pair.     

Ground-state coordination of a catalytic metal to the scissile phosphate of a tertiary-stabilized Hammerhead ribozyme

Although the Hammerhead ribozyme (HHRz) has long been used as a model system in the field of ribozyme enzymology, several details of its mechanism are still not well understood. In particular, significant questions remain concerning the disposition and role of catalytic metals in the HHRz. Previous metal-rescue experiments using a "minimal" HHRz resulted in prediction of a catalytic metal that is bound in the A9/G10.1 site in the ground state of the reaction and that bridges to the scissile phosphate further along the reaction pathway. "Native" or extended HHRz constructs contain tertiary contacts that stabilize a more compact structure at moderate ionic strength. We performed Cd(2+) rescue experiments on an extended HHRz from Schistosoma mansoni using stereo-pure scissile phosphorothioate-substituted substrates in order to determine whether a metal ion makes contact with the scissile phosphate in the ground state or further along the reaction coordinate. Inhibition in Ca(2+)/Mg(2+) and rescue by thiophilic Cd(2+) was specific for the R(p)-S stereoisomer of the scissile phosphate. The affinity of the rescuing Cd(2+), measured in two different ionic strength backgrounds, increased fourfold to 17-fold when the pro-R(p) oxygen is replaced by sulfur. These data support a model in which the rescuing metal ion makes a ground-state interaction with the scissile phosphate in the native HHRz. The resulting model for Mg(2+) activation in the HHRz places a metal ion in contact with the scissile phosphate, where it may provide ground-state electrostatic activation of the substrate.     

Local conformational changes in the catalytic core of the trans-acting hepatitis delta virus ribozyme accompany catalysis

The hepatitis delta virus (HDV) is a human pathogen and satellite RNA of the hepatitis B virus. It utilizes a self-cleaving catalytic RNA motif to process multimeric intermediates in the double-rolling circle replication of its genome. Previous kinetic analyses have suggested that a particular cytosine residue (C(75)) with a pK(a) close to neutrality acts as a general acid or base in cleavage chemistry. The crystal structure of the product form of a cis-acting HDV ribozyme shows this residue positioned close to the 5'-OH leaving group of the reaction by a trefoil turn in the RNA backbone. By modifying G(76) of the trefoil turn of a synthetic trans-cleaving HDV ribozyme to the fluorescent 2-aminopurine (AP), we can directly monitor local conformational changes in the catalytic core. In the ribozyme-substrate complex (precursor), AP fluorescence is strongly quenched, suggesting that AP(76) is stacked with other bases and that the trefoil turn is not formed. In contrast, formation of the product complex upon substrate cleavage or direct product binding results in a significant increase in fluorescence, consistent with AP(76) becoming unstacked and solvent-exposed as evidenced in the trefoil turn. Using AP fluorescence and fluorescence resonance energy transfer (FRET) in concert, we demonstrate that this local conformational change in the trefoil turn is kinetically coincidental with a previously observed global structural change of the ribozyme. Our data show that, at least in the trans-acting HDV ribozyme, C(75) becomes positioned for reaction chemistry only along the trajectory from precursor to product.     

The folding pathway of the genomic hepatitis delta virus ribozyme is dominated by slow folding of the pseudoknots

Hepatitis delta virus (HDV) replicates by a double rolling-circle mechanism that requires self-cleavage by closely related genomic and antigenomic versions of a ribozyme. We have previously shown that the uncleaved genomic ribozyme is subject to a variety of alternative (Alt) pairings. Sequence upstream of the ribozyme can regulate self-cleavage activity by formation of an Alt 1 ribozyme-containing structure that severely inhibits self-cleavage, or a P(-1) self-structure that permits rapid self-cleavage. Here, we test three other alternative pairings: Alt P1, Alt 2, and Alt 3. Alt P1 and Alt 3 contain primarily ribozyme-ribozyme interactions, while Alt 2 involves ribozyme-flanking sequence interaction. A number of single and double mutant ribozymes were prepared to increase or decrease the stability of the alternative pairings, and rates of self-cleavage were determined. Results of these experiments were consistent with the existence of the proposed alternative pairings and their ability to inhibit the overall rate of native ribozyme folding. Local misfolds are treated as internal equilibrium constants in a binding polynomial that modulates the intrinsic rate of cleavage. This model of equilibrium effects of misfolds should be general and apply to other ribozymes. All of the alternative pairings sequester a pseudoknot-forming strand. Folding of ribozymes containing Alt P1 and Alt 2 was accelerated by urea as long as the native ribozyme fold was sufficiently stable, while folding of mutants in which both of these alternative pairings had been removed were not stimulated by urea. This behavior suggests that the pseudoknots form by capture of an unfolded or appropriately rearranged alternative pairing by its complementary native strand. Fast-folding mutants were prepared by either weakening alternative pairings or by strengthening native pairings. A kinetic model was developed that accommodates these features and explains the position of the rate-limiting step for the G11C mutant. Implications of these results for structural and dynamic studies of the uncleaved HDV ribozyme are discussed.     

An active site rearrangement within the Tetrahymena group I ribozyme releases nonproductive interactions and allows formation of catalytic interactions

Biological catalysis hinges on the precise structural integrity of an active site that binds and transforms its substrates and meeting this requirement presents a unique challenge for RNA enzymes. Functional RNAs, including ribozymes, fold into their active conformations within rugged energy landscapes that often contain misfolded conformers. Here we uncover and characterize one such “off-pathway” species within an active site after overall folding of the ribozyme is complete. The Tetrahymena group I ribozyme (E) catalyzes cleavage of an oligonucleotide substrate (S) by an exogenous guanosine (G) cofactor. We tested whether specific catalytic interactions with G are present in the preceding E•S•G and E•G ground-state complexes. We monitored interactions with G via the effects of 2′- and 3′-deoxy (–H) and −amino (–NH₂) substitutions on G binding. These and prior results reveal that G is bound in an inactive configuration within E•G, with the nucleophilic 3′-OH making a nonproductive interaction with an active site metal ion termed M_A and with the adjacent 2′-OH making no interaction. Upon S binding, a rearrangement occurs that allows both –OH groups to contact a different active site metal ion, termed M_C, to make what are likely to be their catalytic interactions. The reactive phosphoryl group on S promotes this change, presumably by repositioning the metal ions with respect to G. This conformational transition demonstrates local rearrangements within an otherwise folded RNA, underscoring RNA''s difficulty in specifying a unique conformation and highlighting Nature''s potential to use local transitions of RNA in complex function.     

Assessment of metal-assisted nucleophile activation in the hepatitis delta virus ribozyme from molecular simulation and 3D-RISM

The hepatitis delta virus ribozyme is an efficient catalyst of RNA 2′-O-transphosphorylation and has emerged as a key experimental system for identifying and characterizing fundamental features of RNA catalysis. Recent structural and biochemical data have led to a proposed mechanistic model whereby an active site Mg²⁺ ion facilitates deprotonation of the O2′ nucleophile, and a protonated cytosine residue (C75) acts as an acid to donate a proton to the O5′ leaving group as noted in a previous study. This model assumes that the active site Mg²⁺ ion forms an inner-sphere coordination with the O2′ nucleophile and a nonbridging oxygen of the scissile phosphate. These contacts, however, are not fully resolved in the crystal structure, and biochemical data are not able to unambiguously exclude other mechanistic models. In order to explore the feasibility of this model, we exhaustively mapped the free energy surfaces with different active site ion occupancies via quantum mechanical/molecular mechanical (QM/MM) simulations. We further incorporate a three-dimensional reference interaction site model for the solvated ion atmosphere that allows these calculations to consider not only the rate associated with the chemical steps, but also the probability of observing the system in the presumed active state with the Mg²⁺ ion bound. The QM/MM results predict that a pathway involving metal-assisted nucleophile activation is feasible based on the rate-controlling transition state barrier departing from the presumed metal-bound active state. However, QM/MM results for a similar pathway in the absence of Mg²⁺ are not consistent with experimental data, suggesting that a structural model in which the crystallographically determined Mg²⁺ is simply replaced with Na⁺ is likely incorrect. It should be emphasized, however, that these results hinge upon the assumption of the validity of the presumed Mg²⁺-bound starting state, which has not yet been definitively verified experimentally, nor explored in depth computationally. Thus, further experimental and theoretical study is needed such that a consensus view of the catalytic mechanism emerges.     

Conservation of helical structure contributes to functional metal ion interactions in the catalytic domain of ribonuclease P RNA

Like protein enzymes, catalytic RNAs contain conserved structure motifs important for function. A universal feature of the catalytic domain of ribonuclease P RNA is a bulged-helix motif within the P1-P4 helix junction. Here, we show that changes in bulged nucleotide identity and position within helix P4 affect both catalysis and substrate binding, while a subset of the mutations resulted only in catalytic defects. We find that the proximity of the bulge to sites of metal ion coordination in P4 is important for catalysis; moving the bulge distal to these sites and deleting it had similarly large effects, while moving it proximal to these sites had only a moderate effect on catalysis. To test whether the effects of the mutations are linked to metal ion interactions, we used terbium-dependent cleavage of the phosphate backbone to probe metal ion-binding sites in the wild-type and mutant ribozymes. We detect cleavages at specific sites within the catalytic domain, including helix P4 and J3/4, which have previously been shown to participate directly in metal ion interactions. Mutations introduced into P4 cause local changes in the terbium cleavage pattern due to alternate metal ion-binding configurations with the helix. In addition, a bulge deletion mutation results in a 100-fold decrease in the single turnover cleavage rate constant at saturating magnesium levels, and a reduced affinity for magnesium ions important for catalysis. In light of the alternate terbium cleavage pattern in P4 caused by bulge deletion, this decreased ability to utilize magnesium ions for catalysis appears to be due to localized structural changes in the ribozyme's catalytic core that weaken metal ion interactions in P4 and J3/4. The information reported here, therefore, provides evidence that the universal conservation of the P4 structure is based in part on optimization of metal ion interactions important for catalysis.     

M42 aminopeptidase catalytic site: the structural and functional role of a strictly conserved aspartate residue

The M42 aminopeptidases are a family of dinuclear aminopeptidases widely distributed in Prokaryotes. They are potentially associated to the proteasome, achieving complete peptide destruction. Their most peculiar characteristic is their quaternary structure, a tetrahedron-shaped particle made of twelve subunits. The catalytic site of M42 aminopeptidases is defined by seven conserved residues. Five of them are involved in metal ion binding which is important to maintain both the activity and the oligomeric state. The sixth conserved residue, a glutamate, is the catalytic base deprotonating the water molecule during peptide bond hydrolysis. The seventh residue is an aspartate whose function remains poorly understood. This aspartate residue, however, must have a critical role as it is strictly conserved in all MH clan enzymes. It forms some kind of catalytic triad with the histidine residue and the metal ion of the M2 binding site. We assess its role in TmPep1050, an M42 aminopeptidase of Thermotoga maritima, through a mutational approach. Asp-62 was substituted with alanine, asparagine, or glutamate residue. The Asp-62 substitutions completely abolished TmPep1050 activity and impeded dodecamer formation. They also interfered with metal ion binding as only one cobalt ion is bound per subunit instead of two. The structure of Asp62Ala variant was solved at 1.5 Å showing how the substitution has an impact on the active site fold. We propose a structural role for Asp-62, helping to stabilize a crucial loop in the active site and to position correctly the catalytic base and a metal ion ligand of the M1 site.     

A pH-dependent conformational change, rather than the chemical step, appears to be rate-limiting in the hammerhead ribozyme cleavage reaction.

We have investigated the chemical basis for a previously observed 7.8 A conformational change in the hammerhead ribozyme that positions the substrate for in-line attack. We have found that the conformational change can only be observed at or above pH 8.5 (in the presence of Co(2+)) and requires the presence of an ionizable 2'-OH at the cleavage site, and note that this observed apparent pK(a) of 8.5 for the conformational change is within experimental error (+/-0.5) of the previously reported apparent kinetic pK(a) of 8.5 for the hammerhead ribozyme in the presence of Co(2+). We have solved two crystal structures of hammerhead ribozymes having 2'-OCH(3) or 2'-F substitutions at the cleavage site and have found that these will not undergo a conformational change equivalent to that observed for the hammerhead ribozyme having an unmodified attacking nucleophile under otherwise identical conditions. We have also characterized the kinetics of cleavage in the crystal. In addition to verifying that the particular sequence of RNA that we crystallized cleaves faster in the crystal than in solution, we also find that the extent of cleavage in the crystal is complete, unlike in solution where this and most other hammerhead ribozyme substrates are cleaved only to about 70 % completion. The initial cleavage rate in the crystal obeys the expected log-linear relation between cleavage-rate and pH with a slope of 0.7, as has been observed for other hammerhead ribozyme sequences in solution, indicating that in both the crystal and in solution the pH-dependent step is rate-limiting. However, the cleavage rate in the crystal is biphasic, with the most dramatic distinction between initial (slower) and final (faster) phases appearing at pH 6.0. The initial phase corresponds to the pH-dependent cleavage rate observed in solution, but the second, faster phase is roughly pH-independent and closely parallels the cleavage rate observed at pH 8 (0.4/minute). This result is particularly remarkable because it entails that the rapidly cleaving phase at pH 6 is comparable to the cleavage rate for the fastest cleaving hammerhead ribozymes at pH 6. Based upon these observations, we conclude that the pH-dependent conformational change is the rate-determining step under standard conditions for the hammerhead ribozyme self-cleavage reaction, and that an ionizable 2'-proton at cleavage site is required for this conformational change. We further hypothesize that deprotonation of the cleavage-site 2'-oxygen drives this conformational change.     

Influence of decreased solvent permittivity on the structure and magnesium(II)-binding properties of the catalytic domain 5 of a group II intron ribozyme

Although it is well known that the so-called “equivalent solution” or “effective” solvent permittivity (dielectric constant) in proteins and nucleic acids is lower than in bulk water, this fact is commonly neglected in (bioinorganic) studies of such compounds. Using domain 5 of the group II intron ribozyme Sc.ai5γ, we describe here the influence of 1,4-dioxane-d₈ on the structure and magnesium(II)-binding properties of this catalytic domain. Applying one- and two-dimensional NMR, we observe distinct structural changes in the functionally important bulge region following a decrease in solvent permittivity. Concomitantly, an increase by a factor of 1.5 in the affinity of Mg²⁺ towards the individual-binding sites in the catalytic core domain is observed upon addition of 1,4-dioxane-d₈. This has led to the detection of a new metal ion coordination site near the GU wobble pair in the catalytic triad. Our results show that solvent permittivity is an important factor in the formation of intrinsic RNA structures and affects their metal ion-binding properties. Hence, solvent permittivity should be taken into account in future studies.     

Subtle but variable conformational rearrangements in the replication cycle of Sulfolobus solfataricus P2 DNA polymerase IV (Dpo4) may accommodate lesion bypass

The possible conformational changes of DNA polymerase IV (Dpo4) before and after the nucleotidyl-transfer reaction are investigated at the atomic level by dynamics simulations to gain insight into the mechanism of low-fidelity polymerases and identify slow and possibly critical steps. The absence of significant conformational changes in Dpo4 before chemistry when the incoming nucleotide is removed supports the notion that the "induced-fit" mechanism employed to interpret fidelity in some replicative and repair DNA polymerases does not exist in Dpo4. However, significant correlated movements in the little finger and finger domains, as well as DNA sliding and subtle catalytic-residue rearrangements, occur after the chemical reaction when both active-site metal ions are released. Subsequently, Dpo4's little finger grips the DNA through two arginine residues and pushes it forward. These metal ion correlated movements may define subtle, and possibly characteristic, conformational adjustments that operate in some Y-family polymerase members in lieu of the prominent subdomain motions required for catalytic cycling in other DNA polymerases like polymerase beta. Such subtle changes do not easily provide a tight fit for correct incoming substrates as in higher-fidelity polymerases, but introduce in low-fidelity polymerases different fidelity checks as well as the variable conformational-mobility potential required to bypass different lesions.     

Crystal structure of glycerophosphodiester phosphodiesterase (GDPD) from Thermoanaerobacter tengcongensis, a metal ion-dependent enzyme: insight into the catalytic mechanism

Glycerophosphodiester phosphodiesterase (GDPD; EC 3.1.4.46) catalyzes the hydrolysis of a glycerophosphodiester to an alcohol and glycerol 3-phosphate in glycerol metabolism. It has an important role in the synthesis of a variety of products that participate in many biochemical pathways. We report the crystal structure of the Thermoanaerobacter tengcongensis GDPD (ttGDPD) at 1.91 A resolution, with a calcium ion and glycerol as a substrate mimic coordinated at this calcium ion (PDB entry 2pz0). The ttGDPD dimer with an intermolecular disulfide bridge and two hydrogen bonds is considered as the potential functional unit. We used site-directed mutagenesis to characterize ttGDPD as a metal ion-dependent enzyme, identified a cluster of residues involved in substrate binding and the catalytic reaction, and we propose a possible general acid-base catalytic mechanism for ttGDPD. Superposing the active site with the homologous structure GDPD from Agrobacterium tumefaciens (PDB entry 1zcc), which binds a sulfate ion in the active site, the sulfate ion can represent the phosphate moiety of the substrate, simulating the binding mode of the true substrate of GDPD.     

Opening of the ADP-bound active site in the ABC transporter ATPase dimer: evidence for a constant contact, alternating sites model for the catalytic cycle

ABC transporters are ubiquitous, ATP-dependent transmembrane pumps. The mechanism by which ATP hydrolysis in the nucleotide-binding domain (NBD) effects conformational changes in the transmembrane domain that lead to allocrite translocation remains largely unknown. A possible aspect of this mechanism was suggested by previous molecular dynamics simulations of the MJ0796 NBD dimer, which revealed a novel, nucleotide-dependent intrasubunit conformational change involving the relative rotation of the helical and catalytic subdomains. Here, we find that in four of five simulations of the ADP/ATP-bound dimer, the relative rotation of the helical and catalytic subdomains in the ADP-bound monomer results in opening of the ADP-bound active site, probably sufficient or close to sufficient to allow nucleotide exchange. We also observe that in all five simulations of the ADP/ATP-bound dimer, the intimate contact of the LSGGQ signature sequence with the ATP gamma-phosphate is weakened by the intrasubunit conformational change within the ADP-bound monomer. We discuss how these results support a constant contact model for the function of the NBD dimer in contrast to switch models, in which the NBDs are proposed to fully disassociate during the catalytic cycle.     

Molecular mechanics simulations of a conformational rearrangement of D-xylose in the active site of D-xylose isomerase.

A proposed reaction mechanism for the enzyme D-xylose isomerase involves the ring opening of the cyclic substrate with a subsequent conformational rearrangement to an extended open-chain form. Restrained energy minimization was used to simulate the rearrangement. In the ring-opening step, the substrate energy function was gradually altered from a cyclic to an open-chain form, with energy minimization after each change. The protein/sugar contact energy did not increase significantly during the process, showing that there was no steric hindrance to ring opening. The conformational rearrangement involves an alteration in the coordination of the substrate to metal ion [1], which was induced by gradually changing restraints on metal/ligand distances. By allowing varying amounts of flexibility in the protein and examining a simplified model system, the interactions of the sugar with metal ion [1] and its immediate ligands were found to be the most important contributors to the energy barrier for the change. Only small changes in the positions of protein atoms were required. The energy barrier to the rearrangement was estimated to be less than the Arrhenius activation energy for the enzymatic reaction. This is in accordance with experimental indications that the isomerization step is rate determining.     

Structural and biochemical characterization of a new Mg(2+) binding site near Tyr94 in the restriction endonuclease PvuII

We have determined the crystal structure of the PvuII endonuclease in the presence of Mg(2+). According to the structural data, divalent metal ion binding in the PvuII subunits is highly asymmetric. The PvuII-Mg(2+) complex has two distinct metal ion binding sites, one in each monomer. One site is formed by the catalytic residues Asp58 and Glu68, and has extensive similarities to a catalytically important site found in all structurally examined restriction endonucleases. The other binding site is located in the other monomer, in the immediate vicinity of the hydroxyl group of Tyr94; it has no analogy to metal ion binding sites found so far in restriction endonucleases. To assign the number of metal ions involved and to better understand the role of Mg(2+) binding to Tyr94 for the function of PvuII, we have exchanged Tyr94 by Phe and characterized the metal ion dependence of DNA cleavage of wild-type PvuII and the Y94F variant. Wild-type PvuII cleaves both strands of the DNA in a concerted reaction. Mg(2+) binding, as measured by the Mg(2+) dependence of DNA cleavage, occurs with a Hill coefficient of 4, meaning that at least two metal ions are bound to each subunit in a cooperative fashion upon formation of the active complex. Quenched-flow experiments show that DNA cleavage occurs about tenfold faster if Mg(2+) is pre-incubated with enzyme or DNA than if preformed enzyme-DNA complexes are mixed with Mg(2+). These results show that Mg(2+) cannot easily enter the active center of the preformed enzyme-DNA complex, but that for fast cleavage the metal ions must already be bound to the apoenzyme and carried with the enzyme into the enzyme-DNA complex. The Y94F variant, in contrast to wild-type PvuII, does not cleave DNA in a concerted manner and metal ion binding occurs with a Hill coefficient of 1. These results indicate that removal of the Mg(2+) binding site at Tyr94 completely disrupts the cooperativity in DNA cleavage. Moreover, in quenched-flow experiments Y94F cleaves DNA about ten times more slowly than wild-type PvuII, regardless of the order of mixing. From these results we conclude that wild-type PvuII cleaves DNA in a fast and concerted reaction, because the Mg(2+) required for catalysis are already bound at the enzyme, one of them at Tyr94. We suggest that this Mg(2+) is shifted to the active center during binding of a specific DNA substrate. These results, for the first time, shed light on the pathway by which metal ions as essential cofactors enter the catalytic center of restriction endonucleases.     

Metal ions and phosphate binding in the H-N-H motif: crystal structures of the nuclease domain of ColE7/Im7 in complex with a phosphate ion and different divalent metal ions

H-N-H is a motif found in the nuclease domain of a subfamily of bacteria toxins, including colicin E7, that are capable of cleaving DNA nonspecifically. This H-N-H motif has also been identified in a subfamily of homing endonucleases, which cleave DNA site specifically. To better understand the role of metal ions in the H-N-H motif during DNA hydrolysis, we crystallized the nuclease domain of colicin E7 (nuclease-ColE7) in complex with its inhibitor Im7 in two different crystal forms, and we resolved the structures of EDTA-treated, Zn(2+)-bound and Mn(2+)-bound complexes in the presence of phosphate ions at resolutions of 2.6 A to 2.0 A. This study offers the first determination of the structure of a metal-free and substrate-free enzyme in the H-N-H family. The H-N-H motif contains two antiparallel beta-strands linked to a C-terminal alpha-helix, with a divalent metal ion located in the center. Here we show that the metal-binding sites in the center of the H-N-H motif, for the EDTA-treated and Mg(2+)-soaked complex crystals, were occupied by water molecules, indicating that an alkaline earth metal ion does not reside in the same position as a transition metal ion in the H-N-H motif. However, a Zn(2+) or Mn(2+) ions were observed in the center of the H-N-H motif in cases of Zn(2+) or Mn(2+)-soaked crystals, as confirmed in anomalous difference maps. A phosphate ion was found to bridge between the divalent transition metal ion and His545. Based on these structures and structural comparisons with other nucleases, we suggest a functional role for the divalent transition metal ion in the H-N-H motif in stabilizing the phosphoanion in the transition state during hydrolysis.     

A 3D structural model and dynamics of hepatitis C virus NS3/4A protease (genotype 4a,strain ED43) suggest conformational instability of the catalytic triad: implications in catalysis and drug resistivity

Egypt has the highest prevalence of hepatitis C virus (HCV) infection worldwide with a frequency of 15%. More than 90% of these infections are due to genotype 4, and the subtype 4a (HCV-4a) predominates. Moreover, due to the increased mobility of people, HCV-4a has recently spread to several European countries. The protease domain of the HCV nonstructural protein 3 (NS3) has been targeted for inhibition by several drugs. This approach has had marked success in inhibiting genotype 1 (HCV-1), the predominant genotype in the USA, Europe, and Japan. However, HCV-4a was found to resist inhibition by a number of these drugs, and little progress has been made to understand the structural basis of its drug resistivity. As a step forward, we sequenced the NS3 HCV-4a protease gene (strain ED43) and subsequently built a 3D structural model threaded through a template crystal structure of HCV-1b NS3 protease. The model protease, HCV-4a, shares 83% sequence identity with the template protease, HCV-1b, and has nearly identical rigid structural features. Molecular dynamics simulations predict similar overall dynamics of the two proteases. However, local dynamics and 4D analysis of the interactions between the catalytic triad residues (His57, Asp81, and Ser139) indicate conformational instability of the catalytic site in HCV-4a NS3 protease. These results suggest that the divergent dynamics behavior, more than the rigid structure, could be related to the altered catalytic activity and drug resistivity seen in HCV-4a.     

Functional coupling between a distal interaction and the cleavage site in bacterial RNase-P-RNA-mediated cleavage

Bacterial RNase P consists of one protein and one RNA [RNase P RNA (RPR)]. RPR can process tRNA precursors correctly in the absence of the protein. Here we have used model hairpin loop substrates corresponding to the acceptor, T-stem, and T-loop of a precursor tRNA to study the importance of the T-loop structure in RPR-alone reaction. T-stem/loop (TSL) interacts with a region in RPR [TSL binding site (TBS)], forming TSL/TBS interaction. Altering the T-loop structure affects both cleavage site selection and rate of cleavage at the correct site + 1 and at the alternative site − 1. The magnitude of variation depended on the structures of the T-loop and the TBS region, with as much as a 150-fold reduction in the rate of cleavage at + 1. Interestingly, for one T-loop structure mutant, no difference in the rate at − 1 was detected compared to cleavage of the substrate with an unchanged T-loop, indicating that, in this case, the altered T-loop structure primarily influences events required for efficient cleavage at the correct site + 1. We also provide data supporting a functional link between a productive TSL/TBS interaction and events at the cleavage site. Collectively, our findings emphasize the interplay between separate regions upon formation of a productive RPR substrate that leads to efficient and accurate cleavage. These new data provide support for an induced-fit mechanism in bacterial RPR-mediated cleavage at the correct site + 1.