Simulating the fidelity and the three Mg mechanism of pol η and clarifying the validity of transition state theory in enzyme catalysis

Structural basis of transcription: RNA polymerase II substrate binding and metal coordination using a free-electron laser

Catalysis and translocation of multisubunit DNA-directed RNA polymerases underlie all cellular mRNA synthesis. RNA polymerase II (Pol II) synthesizes eukaryotic pre-mRNAs from a DNA template strand buried in its active site. Structural details of catalysis at near-atomic resolution and precise arrangement of key active site components have been elusive. Here, we present the free-electron laser (FEL) structures of a matched ATP-bound Pol II and the hyperactive Rpb1 T834P bridge helix (BH) mutant at the highest resolution to date. The radiation-damage-free FEL structures reveal the full active site interaction network, including the trigger loop (TL) in the closed conformation, bonafide occupancy of both site A and B Mg2+, and, more importantly, a putative third (site C) Mg2+ analogous to that described for some DNA polymerases but not observed previously for cellular RNA polymerases. Molecular dynamics (MD) simulations of the structures indicate that the third Mg2+ is coordinated and stabilized at its observed position. TL residues provide half of the substrate binding pocket while multiple TL/BH interactions induce conformational changes that could allow translocation upon substrate hydrolysis. Consistent with TL/BH communication, a FEL structure and MD simulations of the T834P mutant reveal rearrangement of some active site interactions supporting potential plasticity in active site function and long-distance effects on both the width of the central channel and TL conformation, likely underlying its increased elongation rate at the expense of fidelity.

Structures of kinetic intermediate states of HIV-1 reverse transcriptase DNA synthesis

Reverse transcription of the retroviral single-stranded RNA into double-stranded DNA is an integral step during HIV-1 replication, and reverse transcriptase (RT) is a primary target for antiviral therapy. Despite a wealth of structural information on RT, we lack critical insight into the intermediate kinetic states of DNA synthesis. Using catalytically active substrates, and a novel blot/diffusion cryo-electron microscopy approach, we captured 11 structures that define the substrate binding, reactant, transition and product states of dATP addition by RT at 1.9 to 2.4 Å resolution in the active site. Initial dATP binding to RT-template/primer complex involves a single Mg 2+ (site B), and promotes partial closure of the active site pocket by a large conformational change in the β3-β4 loop in the Fingers domain, and formation of a negatively charged pocket where a second "drifting" Mg 2+ can bind (site A). During the transition state, the α-phosphate oxygen from a previously unobserved dATP conformer aligns with the site A Mg 2+ and the primer 3'-OH for nucleophilic attack. In the product state, we captured two substrate conformations in the active site: 1) dATP that had yet to be incorporated into the nascent DNA, and 2) an incorporated dAMP with the pyrophosphate leaving group coordinated by metal B and stabilized through H- bonds in the active site of RT. This study provides insights into a fundamental chemical reaction that impacts polymerase fidelity, nucleoside inhibitor drug design, and mechanisms of drug resistance.

Structural basis of transcription: RNA Polymerase II substrate binding and metal coordination at 3.0 Å using a free-electron laser

Catalysis and translocation of multi-subunit DNA-directed RNA polymerases underlie all cellular mRNA synthesis. RNA polymerase II (Pol II) synthesizes eukaryotic pre-mRNAs from a DNA template strand buried in its active site. Structural details of catalysis at near atomic resolution and precise arrangement of key active site components have been elusive. Here we present the free electron laser (FEL) structure of a matched ATP-bound Pol II, revealing the full active site interaction network at the highest resolution to date, including the trigger loop (TL) in the closed conformation, bonafide occupancy of both site A and B Mg2+, and a putative third (site C) Mg2+ analogous to that described for some DNA polymerases but not observed previously for cellular RNA polymerases. Molecular dynamics (MD) simulations of the structure indicate that the third Mg2+ is coordinated and stabilized at its observed position. TL residues provide half of the substrate binding pocket while multiple TL/bridge helix (BH) interactions induce conformational changes that could propel translocation upon substrate hydrolysis. Consistent with TL/BH communication, a FEL structure and MD simulations of the hyperactive Rpb1 T834P bridge helix mutant reveals rearrangement of some active site interactions supporting potential plasticity in active site function and long-distance effects on both the width of the central channel and TL conformation, likely underlying its increased elongation rate at the expense of fidelity.

Screening, and Optimization of Fermentation Medium to Produce Secondary Metabolites from Bacillus amyloliquefaciens, for the Biocontrol of Early Leaf Spot Disease, and Growth Promoting Effects on Peanut (Arachis hypogaea L.)

A novel Bacillus amyloliquefaciens BAM strain, with novel fermentation nutrient mediums and compositions, could produce potent antifungal secondary metabolites, as the existing strains face resistance from fungus pathogens. In the current study, we introduced two novel nutrient mediums for the fermentation process, semolina and peanut root extract, as carbon and nitrogen sources in order to maximize the antifungal effects of B. amyloliquefaciens against Cercaspora arachidichola to control early leaf spot disease in peanuts. Based on a single-factor test and the central composite design of response surface methodology, the optimum fermentation medium for Bacillus amyloliquefaciens antagonistic substance was determined, containing 15 gm/L of semolina flour, 12.5 gm/L of beef extract, and 0.5 gm/L of magnesium sulfate, which inhibited the fungal growth by 91%. In vitro, antagonistic activity showed that the fermentation broth of B. amyloliquefaciens BAM with the optimized medium formulation had an inhibition rate of (92.62 ± 2.07)% on the growth of C. arachidichola. Disease control effects in pot experiments show that the pre-infection spray of B. amyloliquefaciens BAM broth had significant efficiency of (92.00 ± 3.79)% in comparison to post-infection spray. B. amyloliquefaciens BAM broth significantly promoted peanut plant growth and physiological parameters and reduced the biotic stress of C. archidechola. Studies revealed that B. amyloliquefaciens BAM with a novel fermentation formulation could be an ideal biocontrol and biofertilizer agent and help in early disease management of early leaf spots in peanuts.

On the Recognition of Natural Substrate CTP and Endogenous Inhibitor ddhCTP of SARS-CoV-2 RNA-Dependent RNA Polymerase: A Molecular Dynamics Study

The novel coronavirus SARS-CoV-2 is the causative agent of the COVID-19 outbreak that is affecting the entire planet. As the pandemic is still spreading worldwide, with multiple mutations of the virus, it is of interest and of help to employ computational methods for identifying potential inhibitors of the enzymes responsible for viral replication. Attractive antiviral nucleotide analogue RNA-dependent RNA polymerase (RdRp) chain terminator inhibitors are investigated with this purpose. This study, based on molecular dynamics (MD) simulations, addresses the important aspects of the incorporation of an endogenously synthesized nucleoside triphosphate, ddhCTP, in comparison with the natural nucleobase cytidine triphosphate (CTP) in RdRp. The ddhCTP species is the product of the viperin antiviral protein as part of the innate immune response. The absence of the ribose 3'-OH in ddhCTP could have important implications in its inhibitory mechanism of RdRp. We built an in silico model of the RNA strand embedded in RdRp using experimental methods, starting from the cryo-electron microscopy structure and exploiting the information obtained by spectrometry on the RNA sequence. We determined that the model was stable during the MD simulation time. The obtained results provide deeper insights into the incorporation of nucleoside triphosphates, whose molecular mechanism by the RdRp active site still remains elusive.

Medium optimization to analyze the protein composition of Bacillus pumilus HR10 antagonizing Sphaeropsis sapinea

A previous study found that a biocontrol bacterium, Bacillus pumilus HR10, inhibited the Sphaeropsis shoot blight disease of pine, and the fermentation broth of HR10 strain contained protein antifungal substances. The optimal formulation of the fermentation medium for the antagonistic substance of B. pumilus HR10 was finally obtained by single-factor test, Packett-Burman test, steepest ascent test and Box-Behnken Design (BBD) response surface test, and the best formulation of the fermentation medium for the antagonistic substance of B. pumilus HR10 was 12 g/L corn meal, 15 g/L beef extract and 13 g/L magnesium sulfate, with a predicted bacterial inhibition rate of 89%. The fermentation filtrate of B. pumilus HR10 cultured with the optimized medium formulation was verified to have an inhibition rate of (87.04 ± 3.2) % on the growth of Sphaeropsis sapinea by three replicate tests. The antagonistic crude protein of B. pumilus HR10 were further isolated and identified using HiTrap Capto Q strong Ion-Exchange Chromatography and LC-MS-MS, and it was speculated that glycoside hydrolase (Ghy), beta-glucanase (Beta), arabinogalactan endonuclease β-1,4-galactanase (Arab), and immunosuppressant A (ImA) are proteins with antagonistic activity against S. sapinea in the B. pumilus HR10.

In crystallo observation of three metal ion promoted DNA polymerase misincorporation

Error-free replication of DNA is essential for life. Despite the proofreading capability of several polymerases, intrinsic polymerase fidelity is in general much higher than what base-pairing energies can provide. Although researchers have investigated this long-standing question with kinetics, structural determination, and computational simulations, the structural factors that dictate polymerase fidelity are not fully resolved. Time-resolved crystallography has elucidated correct nucleotide incorporation and established a three-metal-ion-dependent catalytic mechanism for polymerases. Using X-ray time-resolved crystallography, we visualize the complete DNA misincorporation process catalyzed by DNA polymerase η. The resulting molecular snapshots suggest primer 3´-OH alignment mediated by A-site metal ion binding is the key step in substrate discrimination. Moreover, we observe that C-site metal ion binding preceded the nucleotidyl transfer reaction and demonstrate that the C-site metal ion is strictly required for misincorporation. Our results highlight the essential but separate roles of the three metal ions in DNA synthesis.

By observing DNA polymerase misincorporation with time-resolved crystallography, the authors visualize three-metal ion dependent polymerase catalysis and identify A-site metal-mediated primer alignment as a key step in nucleotide discrimination.

Two-Metal-Ion Catalysis: Inhibition of DNA Polymerase Activity by a Third Divalent Metal Ion

Almost all DNA polymerases (pols) exhibit bell-shaped activity curves as a function of both pH and Mg²⁺ concentration. The pol activity is reduced when the pH deviates from the optimal value. When the pH is too low the concentration of a deprotonated general base (namely, the attacking 3′-hydroxyl of the 3′ terminal residue of the primer strand) is reduced exponentially. When the pH is too high the concentration of a protonated general acid (i.e., the leaving pyrophosphate group) is reduced. Similarly, the pol activity also decreases when the concentration of the divalent metal ions deviates from its optimal value: when it is too low, the binding of the two catalytic divalent metal ions required for the full activity is incomplete, and when it is too high a third divalent metal ion binds to pyrophosphate, keeping it in the replication complex longer and serving as a substrate for pyrophosphorylysis within the complex. Currently, there is a controversy about the role of the third metal ion which we will address in this review.

Controlled Trafficking of Multiple and Diverse Cations Prompts Nucleic Acid Hydrolysis

Recent in crystallo reaction intermediates have detailed how nucleic acid hydrolysis occurs in the RNA ribonuclease H1 (RNase H1), a fundamental metalloenzyme involved in maintaining the human genome. At odds with the previous characterization, these in crystallo structures unexpectedly captured multiple metal ions (K⁺ and Mg²⁺) transiently bound in the vicinity of the two-metal-ion active site. Using multi-microsecond all-atom molecular dynamics and free-energy simulations, we investigated the functional implications of the dynamic exchange of multiple K⁺ and Mg²⁺ ions at the RNase H1 reaction center. We found that such ions are timely positioned at non-overlapping locations near the active site, at different stages of catalysis, being crucial for both reactants' alignment and leaving group departure. We also found that this cation trafficking is tightly regulated by variations of the solution's ionic strength and is aided by two conserved second-shell residues, E188 and K196, suggesting a mechanism for the cations' recruitment during catalysis. These results indicate that controlled trafficking of multi-cation dynamics, opportunely prompted by second-shell residues, is functionally essential to the complex enzymatic machinery of the RNase H1. These findings revise the current knowledge on the RNase H1 catalysis and open new catalytic possibilities for other similar metalloenzymes including, but not limited to, CRISPR-Cas9, group II intron ribozyme and the human spliceosome.

Building better polymerases: Engineering the replication of expanded genetic alphabets

DNA polymerases are today used throughout scientific research, biotechnology, and medicine, in part for their ability to interact with unnatural forms of DNA created by synthetic biologists. Here especially, natural DNA polymerases often do not have the "performance specifications" needed for transformative technologies. This creates a need for science-guided rational (or semi-rational) engineering to identify variants that replicate unnatural base pairs (UBPs), unnatural backbones, tags, or other evolutionarily novel features of unnatural DNA. In this review, we provide a brief overview of the chemistry and properties of replicative DNA polymerases and their evolved variants, focusing on the Klenow fragment of Taq DNA polymerase (Klentaq). We describe comparative structural, enzymatic, and molecular dynamics studies of WT and Klentaq variants, complexed with natural or noncanonical substrates. Combining these methods provides insight into how specific amino acid substitutions distant from the active site in a Klentaq DNA polymerase variant (ZP Klentaq) contribute to its ability to replicate UBPs with improved efficiency compared with Klentaq. This approach can therefore serve to guide any future rational engineering of replicative DNA polymerases.

Selective Metal Ion Utilization Contributes to the Transformation of the Activity of Yeast Polymerase η from DNA Polymerization toward RNA Polymerization

Polymerase eta (Polη) is a translesion synthesis DNA polymerase directly linked to cancer development. It can bypass several DNA lesions thereby rescuing DNA damage-stalled replication complexes. We previously presented evidence implicating Saccharomyces cerevisiae Polη in transcription elongation, and identified its specific RNA extension and translesion RNA synthetic activities. However, RNA synthesis by Polη proved rather inefficient under conditions optimal for DNA synthesis. Searching for factors that could enhance its RNA synthetic activity, we have identified the divalent cation of manganese. Here, we show that manganese triggers drastic changes in the activity of Polη. Kinetics experiments indicate that manganese increases the efficiency of ribonucleoside incorporation into RNA by ~400–2000-fold opposite undamaged DNA, and ~3000 and ~6000-fold opposite TT dimer and 8oxoG, respectively. Importantly, preference for the correct base is maintained with manganese during RNA synthesis. In contrast, activity is strongly impaired, and base discrimination is almost lost during DNA synthesis by Polη with manganese. Moreover, Polη shows strong preference for manganese during RNA synthesis even at a 25-fold excess magnesium concentration. Based on this, we suggest that a new regulatory mechanism, selective metal cofactor utilization, modulates the specificity of Polη helping it to perform distinct activities needed for its separate functions during replication and transcription.

Cations in motion: QM/MM studies of the dynamic and electrostatic roles of H+ and Mg2+ ions in enzyme reactions

Here we discuss current trends in the simulations of enzymatic reactions focusing on phosphate catalysis. The mechanistic details of the proton transfers coupled to the phosphate cleavage is one of the key challenges in QM/MM calculations of these and other enzyme catalyzed reactions. The lack of experimental information offers both an opportunity for computations as well as often unresolved controversies. We discuss the example of small GTPases including the important human Ras protein. The high dimensionality and chemical complexity of these reactions demand carefully chosen computational techniques both in terms of the underlying quantum chemical theory and the sampling of the conformational ensemble. We also point out the important role of Mg²⁺ ions, and recent advances in their transient involvement in the catalytic mechanisms.

Characterization of six recombinant human RNase H2 bearing Aicardi-Goutiéres syndrome causing mutations

Mammalian RNase H2 is a heterotrimeric enzyme consisting of one catalytic subunit (A) and two accessory subunits (B and C). RNase H2 is involved in the removal of a single ribonucleotide embedded in genomic DNA and removal of RNA of RNA/DNA hybrids. In humans, mutation of the RNase H2 gene causes a severe neuroinflammatory disorder Aicardi-Goutières syndrome (AGS). Here, we examined the activity and stability of six recombinant human RNase H2 variants bearing one AGS-causing mutation, A-G37S (Gly37 in the A subunit is replaced with Ser), A-N212I, A-R291H, B-A177T, B-V185G, or C-R69W. The activity of A-G37S was 0.3-1% of that of the wild-type RNase H2 (WT), while those of other five variants were 51-120%. In circular dichroism measurement, the melting temperatures of variants were 50-53°C, lower than that of WT (56°C). These results suggested that A-G37S had decreased activity and stability than WT, while other five variants had decreased stability but retained activity. In gel filtration chromatography of the purified enzyme preparation, WT migrated as a heterotrimer, while A-R291H eluted in two separate peaks containing either the heterotrimer or only the A subunit, suggesting that some AGS-causing mutations affect the heterotrimer-forming stability of RNase H2.

Extensive free-energy simulations identify water as the base in nucleotide addition by DNA polymerase

Transphosphorylation of nucleotide triphosphates is the central reaction in DNA replication by DNA polymerase as well as many other biological processes. Despite its importance, the microscopic chemical mechanism of transphosphorylation of nucleotide triphosphates is, in most cases, unknown. Here we use extensive simulations of DNA polymerase η to test mechanistic hypotheses. We systematically survey the reactive space by calculating 2D free-energy surfaces for 10 different plausible mechanisms that have been proposed. We supplement these free-energy surfaces with calculations of pKa for a number of potentially acidic protons in different states relevant to the catalytic cycle. We find that among all of the conditions that we test, the smallest activation barrier occurs for a reaction where a Mg2+-coordinated water deprotonates the nucleophilic 3'-OH, and this deprotonation is concerted with the phosphoryl transfer. The presence of a third Mg2+ in the active site lowers the activation barrier for the water-as-base mechanism, as does protonation of the pyrophosphate leaving group, which is consistent with general acid catalysis. The results demonstrate the value of simulations, when used in conjunction with experimental data, to help establish a microscopic chemical mechanism in a complex environment.

A Transition-State Perspective on Y-Family DNA Polymerase η Fidelity in Comparison with X-Family DNA Polymerases λ and β

Deoxynucleotide misincorporation efficiencies can span a wide 10⁴-fold range, from ∼10^-2 to ∼10^-6, depending principally on polymerase (pol) identity and DNA sequence context. We have addressed DNA pol fidelity mechanisms from a transition-state (TS) perspective using our "tool-kit" of dATP- and dGTP-β,γ substrate analogues in which the pyrophosphate leaving group (p K_a4 = 8.9) has been replaced by a series of bisphosphonates covering a broad acidity range spanning p K_a4 values from 7.8 (CF₂) to 12.3 [C(CH₃)₂]. Here, we have used a linear free energy relationship (LFER) analysis, in the form of a Brønsted plot of log( k_pol) versus p K_a4, for Y-family error-prone pol η and X-family pols λ and β to determine the extent to which different electrostatic active site environments alter k_pol values. The apparent chemical rate constant ( k_pol) is the rate-determining step for the three pols. The pols each exhibit a distinct catalytic signature that differs for formation of right (A·T) and wrong (G·T) incorporations observed as changes in slopes and displacements of the Brønsted lines, in relation to a reference LFER. Common to this signature among all three pols is a split linear pattern in which the analogues containing two halogens show k_pol values that are systematically lower than would be predicted from their p K_a4 values measured in aqueous solution. We discuss how metal ions and active site amino acids are responsible for causing "effective" p K_a4 values that differ for dihalo and non-dihalo substrates as well as for individual R and S stereoisomers for CHF and CHCl.

Uncovering a unique approach for damaged DNA replication: A computational investigation of a mutagenic tobacco-derived thymine lesion

4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone is a potent nicotine carcinogen that leads to many DNA lesions, the most persistent being the O2-[4-oxo-4-(3-pyridyl)butyl]thymine adduct (POB-T). Although the experimental mutagenic profile for the minor groove POB-T lesion has been previously reported, the findings are puzzling in terms of the human polymerases involved. Specifically, while pol κ typically replicates minor groove adducts, in vivo studies indicate pol η replicates POB-T despite being known for processing major groove adducts. Our multiscale modeling approach reveals that the canonical (anti) glycosidic orientation of POB-T can fit in the pol κ active site, but only a unique (syn) POB-T conformation is accommodated by pol η. These distinct binding orientations rationalize the differential in vitro mutagenic spectra based on the preferential stabilization of dGTP and dTTP opposite the lesion for pol κ and η, respectively. Overall, by uncovering the first evidence for the replication of a damaged pyrimidine in the syn glycosidic orientation, the current work provides the insight necessary to clarify a discrepancy in the DNA replication literature, expand the biological role of the critical human pol η, and understand the mutational signature in human cancers associated with tobacco exposure.

Crystallographic evidence for two-metal-ion catalysis in human pol η

Extensive evidence exists that DNA polymerases use two metal ions to catalyze the phosphoryl transfer reaction. Recently, competing evidence emerged, suggesting that a third metal ion, known as MnC, may be involved in catalysis. The binding of MnC was observed in crystal structures of the replication complexes of human polymerase (pol) η, pol β, and pol μ. Its occupancy (q_MnC ) in the pol η replication complexes exhibited a strong correlation with the occupancy of the formed product pyrophosphate (q_PPi ), i.e., q_MnC ∝ q_PPi . However, a key piece of information was missing that is needed to distinguish between two possible sequences of events: (i) the chemical reaction occurs first with only two meal ions, followed by the binding of MnC in a "catch-the-product" mode; and (ii) MnC binds first, followed by the chemical reaction with all three metal ions in a "push-the-reaction-forward" mode. Both mechanisms can lead to a strong correlation between q_MnC and q_PPi . However, q_MnC ≤ q_PPi in the first scenario, whereas q_MnC ≥ q_PPi in the second. In this study, an analysis of crystallographic data published recently for pol η complexes shows that the formation of the product pyrophosphate definitely precedes the binding of MnC. Therefore, just like all other DNA polymerases, human pol η employs a two-metal-ion catalytic mechanism. Rather than help to catalyze the reaction, MnC stabilizes the formed product, which remains trapped inside the crystals, before it slowly diffuses out.

Exploring the Role of the Third Active Site Metal Ion in DNA Polymerase η with QM/MM Free Energy Simulations

The enzyme human DNA polymerase η (Pol η) is critical for bypassing lesions during DNA replication. In addition to the two Mg²⁺ ions aligning the active site, experiments suggest that a third Mg²⁺ ion could play an essential catalytic role. Herein the role of this third metal ion is investigated with quantum mechanical/molecular mechanical (QM/MM) free energy simulations of the phosphoryl transfer reaction and a proposed self-activating proton transfer from the incoming nucleotide to the pyrophosphate leaving group. The simulations with only two metal ions in the active site support a sequential mechanism, with phosphoryl transfer followed by relatively fast proton transfer. The simulations with three metal ions in the active site suggest that the third metal ion may play a catalytic role through electrostatic interactions with the leaving group. These electrostatic interactions stabilize the product, making the phosphoryl transfer reaction more thermodynamically favorable with a lower free energy barrier relative to the activated state corresponding to the deprotonated 3'OH nucleophile, and also inhibit the subsequent proton transfer. The possibility that Mg²⁺-bound hydroxide acts as the base deprotonating the 3'OH nucleophile is also explored.

Structural Articulation of Biochemical Reactions Using Restrained Geometries and Topology Switching

A strategy named "restrained geometries and topology switching" (RGATS) is presented to obtain detailed trajectories for complex biochemical reactions using molecular mechanics (MM) methods. It enables prediction of realistic dynamical pathways for chemical reactions, especially for accurately characterizing the structural adjustments of highly complex environments to any proximal biochemical reaction. It can be used to generate reactive conformations, model stepwise or concerted reactions in complex environments, and probe the influence of changes in the environment. Its ability to take reactively nonoptimal conformations and generate favorable starting conformations for a biochemical reaction is illustrated for a proton transfer between two model compounds. Its ability to study concerted reactions in explicit solvent is illustrated using proton transfers between an ammonium ion and two conserved histidines in an ammonia transporter channel embedded in a lipid membrane. Its ability to characterize the changes induced by subtle differences in the active site environment is illustrated using nucleotide addition by a DNA polymerase in the presence of two versus three Mg²⁺ ions. RGATS can be employed within any MM program and requires no additional software implementation. This allows the full assortment of computational methods implemented in all available MM programs to be used to tackle virtually any question about biochemical reactions that is answerable without using a quantum mechanical (QM) model. It can also be applied to generate reasonable starting structures for more detailed and expensive QM or QM/MM methods. In particular, this strategy enables rapid prediction of reactant, intermediary, or product state structures in any macromolecular context, with the only requirement being that the structure in any one of these states is either known or can be accurately modeled.