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

Dynamic interchange between two protonation states is characteristic of active sites of cholinesterases

Cholinesterases are well-known and widely studied enzymes crucial to human health and involved in neurology, Alzheimer's, and lipid metabolism. The protonation pattern of active sites of cholinesterases influences all the chemical processes within, including reaction, covalent inhibition by nerve agents, and reactivation. Despite its significance, our comprehension of the fine structure of cholinesterases remains limited. In this study, we employed enhanced-sampling quantum-mechanical/molecular-mechanical calculations to show that cholinesterases predominantly operate as dynamic mixtures of two protonation states. The proton transfer between two non-catalytic glutamate residues follows the Grotthuss mechanism facilitated by a mediator water molecule. We show that this uncovered complexity of active sites presents a challenge for classical molecular dynamics simulations and calls for special treatment. The calculated proton transfer barrier of 1.65 kcal/mol initiates a discussion on the potential existence of two coupled low-barrier hydrogen bonds in the inhibited form of butyrylcholinesterase. These findings expand our understanding of structural features expressed by highly evolved enzymes and guide future advances in cholinesterase-related protein and drug design studies.

Effects of the Y432S Cancer-Associated Variant on the Reaction Mechanism of Human DNA Polymerase κ

Human DNA polymerases are vital for genetic information management. Their function involves catalyzing the synthesis of DNA strands with unparalleled accuracy, which ensures the fidelity and stability of the human genomic blueprint. Several disease-associated mutations and their functional impact on DNA polymerases have been reported. One particular polymerase, human DNA polymerase kappa (Pol κ), has been reported to be susceptible to several cancer-associated mutations. The Y432S mutation in Pol κ, associated with various cancers, is of interest due to its impact on polymerization activity and markedly reduced thermal stability. Here, we have used computational simulations to investigate the functional consequences of the Y432S using classical molecular dynamics (MD) and coupled quantum mechanics/molecular mechanics (QM/MM) methods. Our findings suggest that Y432S induces structural alterations in domains responsible for nucleotide addition and ternary complex stabilization while retaining structural features consistent with possible catalysis in the active site. Calculations of the minimum energy path associated with the reaction mechanism of the wild type (WT) and Y432S Pol κ indicate that, while both enzymes are catalytically competent (in terms of energetics and the active site's geometries), the cancer mutation results in an endoergic reaction and an increase in the catalytic barrier. Interactions with a third magnesium ion and environmental effects on nonbonded interactions, particularly involving key residues, contribute to the kinetic and thermodynamic distinctions between the WT and mutant during the catalytic reaction. The energetics and electronic findings suggest that active site residues favor the catalytic reaction with dCTP3- over dCTP4-.

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.

Probing the general base for DNA polymerization in telomerase: a molecular dynamics investigation

Telomerase is an RNA-dependent DNA polymerase that plays a role in the maintenance of the 3' end of the eukaryotic chromosome, known as a telomere, by catalyzing the DNA polymerization reaction in cancer and embryonic stem cells. The detailed molecular details of the DNA polymerization by telomerase, especially the general base for deprotonating the terminal 3'-hydroxyl, which triggers the chemical reaction, remain elusive. We conducted a computational investigation using hybrid quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations to probe the detailed mechanism of the reaction. Our simulations started with the telomerase:RNA:DNA:dNTP ternary complex, and by using enhanced sampling QM/MM MD simulations, we probed the general base involved directly in the polymerization. We report the participation of an aspartate (Asp344) coordinated to Mg and an active site water molecule, jointly acting as a base during nucleic acid addition. The Asp344 residue remains transiently protonated during the course of the reaction, and later it deprotonates by transferring its proton to the water at the end of the reaction.

Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.

Machines on Genes through the Computational Microscope

Macromolecular machines acting on genes are at the core of life's fundamental processes, including DNA replication and repair, gene transcription and regulation, chromatin packaging, RNA splicing, and genome editing. Here, we report the increasing role of computational biophysics in characterizing the mechanisms of "machines on genes", focusing on innovative applications of computational methods and their integration with structural and biophysical experiments. We showcase how state-of-the-art computational methods, including classical and ab initio molecular dynamics to enhanced sampling techniques, and coarse-grained approaches are used for understanding and exploring gene machines for real-world applications. As this review unfolds, advanced computational methods describe the biophysical function that is unseen through experimental techniques, accomplishing the power of the "computational microscope", an expression coined by Klaus Schulten to highlight the extraordinary capability of computer simulations. Pushing the frontiers of computational biophysics toward a pragmatic representation of large multimegadalton biomolecular complexes is instrumental in bridging the gap between experimentally obtained macroscopic observables and the molecular principles playing at the microscopic level. This understanding will help harness molecular machines for medical, pharmaceutical, and biotechnological purposes.

Challenges in Protein QM/MM Simulations with Intra-Backbone Link Atoms

Hybrid quantum mechanical/molecular mechanical (QM/MM) simulations fuel discoveries in many fields of science including computational biochemistry and enzymology. Development of more convenient tools leads to an increase in the number of works in which mechanical insights into enzymes' mode of operation are obtained. Most commonly, these tools feature hydrogen-capping (link atom) approach to provide coupling between QM and MM subsystems across a covalent bond. Extensive studies were conducted to provide a solid foundation for the correctness of such an approach when a bond to a nonpolar MM atom is considered. However, not every task may be accomplished this way. Certain scenarios of using QM/MM in computational enzymology encourage or even necessitate the incorporation of backbone atoms into the QM region. Two out of three backbone atoms are polar, and in QM/MM with electrostatic embedding, a neighboring link atom will be hyperpolarized. Several schemes to mitigate this effect were previously proposed alongside a rigorous assessment of quantitative effects on model systems. However, it was not clear whether they may translate into qualitatively different results and how link atom hyperpolarization may manifest itself in a real-life enzymological scenario. Here, we show that the consequences of such an artifact may be severe and may completely overturn the conclusions drawn from the simulations. Our case advocates for the use of charge redistribution schemes whenever intra-backbone QM/MM boundaries are considered. Moreover, we addressed how different boundary types and charge redistribution schemes influence backbone dynamics. We showed that the results are heavily dependent on which boundary MM terms are retained, with charge alteration being of secondary importance. In the worst case, only three intra-backbone boundaries may be used with relative confidence in the adequacy of resulting simulations, irrespective of the hyperpolarization mitigation scheme. Thus, advances in the field are certainly needed to fuel new discoveries. As of now, we believe that issues raised in this work might encourage authors in the field to report what boundaries, boundary MM terms, and charge redistribution schemes they are using, so their results may be correctly interpreted.

Between Protein Fold and Nucleophile Identity: Multiscale Modeling of the TEV Protease Enzyme-Substrate Complex

The cysteine protease from the tobacco etch virus (TEVp) is a well-known and widely utilized enzyme. TEVp's chymotrypsin-like fold is generally associated with serine catalytic triads that differ in terms of a reaction mechanism from the most well-studied papain-like cysteine proteases. The question of what dominates the TEVp mechanism, nucleophile identity, or structural composition has never been previously addressed. Here, we use enhanced sampling multiscale modeling to uncover that TEVp combines the features of two worlds in such a way that potentially hampers its activity. We show that TEVp cysteine is strictly in the anionic form in a free enzyme similar to papain. Peptide binding shifts the equilibrium toward the nucleophile's protonated form, characteristic of chymotrypsin-like proteases, although the cysteinyl anion form is still present and interconversion is rapid. This way cysteine protonation generates enzyme states that are a diversion from the most effective course of action, with only 13.2% of Michaelis complex sub-states able to initiate the reaction. As a result, we propose an updated view on the reaction mechanism catalyzed by TEVp. We also demonstrate that AlphaFold is able to construct protease-substrate complexes with high accuracy. We propose that our findings open a way for its industrious use in enzymological tasks. Unique features of TEVp discovered in this work open a discussion on the evolutionary history and trade-offs of optimizing serine triad-associated folds to cysteine as a nucleophile.

Current state of knowledge of human DNA polymerase eta protein structure and disease-causing mutations

POLη, encoded by the POLH gene, is a crucial protein for replicating damaged DNA and the most studied specialized translesion synthesis polymerases. Mutations in POLη are associated with cancer and the human syndrome xeroderma pigmentosum variant, which is characterized by extreme photosensitivity and an increased likelihood of developing skin cancers. The myriad of structural information about POLη is vast, covering dozens of different mutants, numerous crucial residues, domains, and posttranslational modifications that are essential for protein function within cells. Since POLη is key vital enzyme for cell survival, and mutations in this protein are related to aggressive diseases, understanding its structure is crucial for biomedical sciences, primarily due to its similarities with other Y-family polymerases and its potential as a targeted therapy-drug for tumors. This work provides an up-to-date review on structural aspects of the human POLη: from basic knowledge about critical residues and protein domains to its mutant variants, posttranslational modifications, and our current understanding of therapeutic molecules that target POLη. Thus, this review provides lessons about POLη's structure and gathers critical discussions and hypotheses that may contribute to understanding this protein's vital roles within the cells.

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.

Factors That Determine the Variation of Equilibrium and Kinetic Properties of QM/MM Enzyme Simulations: QM Region, Conformation, and Boundary Condition

To analyze the impact of various technical details on the results of quantum mechanical (QM)/molecular mechanical (MM) enzyme simulations, including the QM region size, catechol-O-methyltransferase (COMT) is studied as a model system using an approximate QM/MM method (DFTB3/CHARMM). The results show that key equilibrium and kinetic properties for methyl transfer in COMT exhibit limited variations with respect to the size of the QM region, which ranges from ∼100 to ∼500 atoms in this study. With extensive sampling, local and global structural characteristics of the enzyme are largely conserved across the studied QM regions, while the nature of the transition state (e.g., secondary kinetic isotope effect) and reaction exergonicity are largely maintained. Deviations in the free energy profile with different QM region sizes are similar in magnitude to those observed with changes in other simulation protocols, such as different initial enzyme conformations and boundary conditions. Electronic structural properties, such as the covariance matrix of residual charge fluctuations, appear to exhibit rather long-range correlations, especially when the peptide backbone is included in the QM region; this observation holds when a range-separated DFT approach is used as the QM region, suggesting that delocalization error is unlikely the origin. Overall, the analyses suggest that multiple simulation details determine the results of QM/MM enzyme simulations with comparable contributions.

Probing the Suitability of Different Ca2+ Parameters for Long Simulations of Diisopropyl Fluorophosphatase

Organophosphate hydrolases are promising as potential biotherapeutic agents to treat poisoning with pesticides or nerve gases. However, these enzymes often need to be further engineered in order to become useful in practice. One example of such enhancement is the alteration of enantioselectivity of diisopropyl fluorophosphatase (DFPase). Molecular modeling techniques offer a unique opportunity to address this task rationally by providing a physical description of the substrate-binding process. However, DFPase is a metalloenzyme, and correct modeling of metal cations is a challenging task generally coming with a tradeoff between simulation speed and accuracy. Here, we probe several molecular mechanical parameter combinations for their ability to empower long simulations needed to achieve a quantitative description of substrate binding. We demonstrate that a combination of the Amber19sb force field with the recently developed 12-6 Ca2+ models allows us to both correctly model DFPase and obtain new insights into the DFP binding process.

Multiple deprotonation paths of the nucleophile 3'-OH in the DNA synthesis reaction

DNA synthesis by polymerases is essential for life. Deprotonation of the nucleophile 3'-OH is thought to be the obligatory first step in the DNA synthesis reaction. We have examined each entity surrounding the nucleophile 3'-OH in the reaction catalyzed by human DNA polymerase (Pol) η and delineated the deprotonation process by combining mutagenesis with steady-state kinetics, high-resolution structures of in crystallo reactions, and molecular dynamics simulations. The conserved S113 residue, which forms a hydrogen bond with the primer 3'-OH in the ground state, stabilizes the primer end in the active site. Mutation of S113 to alanine destabilizes primer binding and reduces the catalytic efficiency. Displacement of a water molecule that is hydrogen bonded to the 3'-OH using the 2'-OH of a ribonucleotide or 2'-F has little effect on catalysis. Moreover, combining the S113A mutation with 2'-F replacement, which removes two potential hydrogen acceptors of the 3'-OH, does not reduce the catalytic efficiency. We conclude that the proton can leave the O3' via alternative paths, supporting the hypothesis that binding of the third Mg²⁺ initiates the reaction by breaking the α-β phosphodiester bond of an incoming deoxyribonucleoside triphosphate (dNTP).

Cooperative protein-solvent tuning of proton transfer energetics: carbonic anhydrase as a case study

We investigate the coupling between the proton transfer (PT) energetics and the protein-solvent dynamics using the intra-molecular PT in wild type (wt) human carbonic anhydrase II and its ten-fold faster mutant Y7F/N67Q as a test case. We calculate the energy variation upon PT, and from that we also calculate the PT reaction free energy, making use of a hybrid quantum mechanics/molecular dynamics approach. In agreement with the experimental data, we obtain that the reaction free energy is basically the same in the two systems. Yet, we show that the instantaneous PT energy is on average lower in the mutant possibly contributing to the faster PT rate. Analysis of the contribution to the PT energetics of the solvent and of each protein residue, also not in the vicinity of the active site, provides evidence for electrostatic tuning of the PT energy arising from the combined effect of the solvent and the protein environment. These findings open up a way to the more general task of the rational design of mutants with either enhanced or reduced PT rate.

Biomolecular QM/MM Simulations: What Are Some of the "Burning Issues"?

QM/MM simulations have become an indispensable tool in many chemical and biochemical investigations. Considering the tremendous degree of success, including recognition by a 2013 Nobel Prize in Chemistry, are there still "burning challenges" in QM/MM methods, especially for biomolecular systems? In this short Perspective, we discuss several issues that we believe greatly impact the robustness and quantitative applicability of QM/MM simulations to many, if not all, biomolecules. We highlight these issues with observations and relevant advances from recent studies in our group and others in the field. Despite such limited scope, we hope the discussions are of general interest and will stimulate additional developments that help push the field forward in meaningful directions.

Mechanism of Type IA Topoisomerases

Topoisomerases in the type IA subfamily can catalyze change in topology for both DNA and RNA substrates. A type IA topoisomerase may have been present in a last universal common ancestor (LUCA) with an RNA genome. Type IA topoisomerases have since evolved to catalyze the resolution of topological barriers encountered by genomes that require the passing of nucleic acid strand(s) through a break on a single DNA or RNA strand. Here, based on available structural and biochemical data, we discuss how a type IA topoisomerase may recognize and bind single-stranded DNA or RNA to initiate its required catalytic function. Active site residues assist in the nucleophilic attack of a phosphodiester bond between two nucleotides to form a covalent intermediate with a 5'-phosphotyrosine linkage to the cleaved nucleic acid. A divalent ion interaction helps to position the 3'-hydroxyl group at the precise location required for the cleaved phosphodiester bond to be rejoined following the passage of another nucleic acid strand through the break. In addition to type IA topoisomerase structures observed by X-ray crystallography, we now have evidence from biophysical studies for the dynamic conformations that are required for type IA topoisomerases to catalyze the change in the topology of the nucleic acid substrates.

Single-step Replacement of an Unreactive C-H Bond by a C-S Bond Using Polysulfide as the Direct Sulfur Source in Anaerobic Ergothioneine Biosynthesis

Ergothioneine, a natural longevity vitamin and antioxidant, is a thiol-histidine derivative. Recently, two types of biosynthetic pathways were reported. In the aerobic ergothioneine biosynthesis, a non-heme iron enzyme incorporates a sulfoxide to an sp² C-H bond in trimethyl-histidine (hercynine) through oxidation reactions. In contrast, in the anaerobic ergothioneine biosynthetic pathway in a green sulfur bacterium, Chlorobium limicola, a rhodanese domain containing protein (EanB) directly replaces this unreactive hercynine C-H bond with a C-S bond. Herein, we demonstrate that polysulfide (HSS_nSR) is the direct sulfur-source in EanB-catalysis. After identifying EanB's substrates, X-ray crystallography of several intermediate states along with mass spectrometry results provide additional mechanistic details for this reaction. Further, quantum mechanics/molecular mechanics (QM/MM) calculations reveal that protonation of N_π of hercynine by Tyr353 with the assistance of Thr414 is a key activation step for the hercynine sp² C-H bond in this trans-sulfuration reaction.