Active site dynamics and combined quantum mechanics/molecular mechanics (QM/MM) modelling of a HIV-1 reverse transcriptase/DNA/dTTP complex

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

Modeling the Enzymatic Mechanism of the SARS-CoV-2 RNA-Dependent RNA Polymerase by DFT/MM-MD: An Unusual Active Site Leading to High Replication Rates

Viral infection relies on the hijacking of cellular machineries to enforce the reproduction of the infecting virus and its subsequent diffusion. In this context, the replication of the viral genome is a key step performed by specific enzymes, i.e., polymerases. The replication of SARS-CoV-2, the causative agent of the COVID-19 pandemics, is based on the duplication of its RNA genome, an action performed by the viral RNA-dependent RNA polymerase. In this contribution, by using highly demanding DFT/MM-MD computations coupled to 2D-umbrella sampling techniques, we have determined the chemical mechanisms leading to the inclusion of a nucleotide in the nascent viral RNA strand. These results highlight the high efficiency of the polymerase, which lowers the activation free energy to less than 10 kcal/mol. Furthermore, the SARS-CoV-2 polymerase active site is slightly different from those usually found in other similar enzymes, and in particular, it lacks the possibility to enforce a proton shuttle via a nearby histidine. Our simulations show that this absence is partially compensated by lysine whose proton assists the reaction, opening up an alternative, but highly efficient, reactive channel. Our results present the first mechanistic resolution of SARS-CoV-2 genome replication at the DFT/MM-MD level and shed light on its unusual enzymatic reactivity paving the way for the future rational design of antivirals targeting emerging RNA viruses.

Computing Metal-Binding Proteins for Therapeutic Benefit

Over one third of biomolecules rely on metal ions to exert their cellular functions. Metal ions can play a structural role by stabilizing the structure of biomolecules, a functional role by promoting a wide variety of biochemical reactions, and a regulatory role by acting as messengers upon binding to proteins regulating cellular metal-homeostasis. These diverse roles in biology ascribe critical implications to metal-binding proteins in the onset of many diseases. Hence, it is of utmost importance to exhaustively unlock the different mechanistic facets of metal-binding proteins and to harness this knowledge to rationally devise novel therapeutic strategies to prevent or cure pathological states associated with metal-dependent cellular dysfunctions. In this compendium, we illustrate how the use of a computational arsenal based on docking, classical, and quantum-classical molecular dynamics simulations can contribute to extricate the minutiae of the catalytic, transport, and inhibition mechanisms of metal-binding proteins at the atomic level. This knowledge represents a fertile ground and an essential prerequisite for selectively targeting metal-binding proteins with small-molecule inhibitors aiming to (i) abrogate deregulated metal-dependent (mis)functions or (ii) leverage metal-dyshomeostasis to selectively trigger harmful cells death.

Mechanistic Insights on Human Phosphoglucomutase Revealed by Transition Path Sampling and Molecular Dynamics Calculations

Human α-phosphoglucomutase 1 (α-PGM) catalyzes the isomerization of glucose-1-phosphate into glucose-6-phosphate (G6P) through two sequential phosphoryl transfer steps with a glucose-1,6-bisphosphate (G16P) intermediate. Given that the release of G6P in the gluconeogenesis raises the glucose output levels, α-PGM represents a tempting pharmacological target for type 2 diabetes. Here, we provide the first theoretical study of the catalytic mechanism of human α-PGM. We performed transition-path sampling simulations to unveil the atomic details of the two catalytic chemical steps, which could be key for developing transition state (TS) analogue molecules with inhibitory properties. Our calculations revealed that both steps proceed through a concerted S_N 2-like mechanism, with a loose metaphosphate-like TS. Even though experimental data suggests that the two steps are identical, we observed noticeable differences: 1) the transition state ensemble has a well-defined TS region and a late TS for the second step, and 2) larger coordinated protein motions are required to reach the TS of the second step. We have identified key residues (Arg23, Ser117, His118, Lys389), and the Mg²⁺ ion that contribute in different ways to the reaction coordinate. Accelerated molecular dynamics simulations suggest that the G16P intermediate may reorient without leaving the enzymatic binding pocket, through significant conformational rearrangements of the G16P and of specific loop regions of the human α-PGM.

Quantum mechanics implementation in drug-design workflows: does it really help?

The pharmaceutical industry is progressively operating in an era where development costs are constantly under pressure, higher percentages of drugs are demanded, and the drug-discovery process is a trial-and-error run. The profit that flows in with the discovery of new drugs has always been the motivation for the industry to keep up the pace and keep abreast with the endless demand for medicines. The process of finding a molecule that binds to the target protein using in silico tools has made computational chemistry a valuable tool in drug discovery in both academic research and pharmaceutical industry. However, the complexity of many protein-ligand interactions challenges the accuracy and efficiency of the commonly used empirical methods. The usefulness of quantum mechanics (QM) in drug-protein interaction cannot be overemphasized; however, this approach has little significance in some empirical methods. In this review, we discuss recent developments in, and application of, QM to medically relevant biomolecules. We critically discuss the different types of QM-based methods and their proposed application to incorporating them into drug-design and -discovery workflows while trying to answer a critical question: are QM-based methods of real help in drug-design and -discovery research and industry?

Metadynamics supports molecular dynamics simulation-based binding affinities of eucalyptol and beta-cyclodextrin inclusion complexes

The development of various molecular dynamics methods enables the detailed investigation of association processes, like host–guest complexes, including their dynamics and, additionally, the release of the guest compound.

Role of R292K mutation in influenza H7N9 neuraminidase toward oseltamivir susceptibility: MD and MM/PB(GB)SA study

The H7N9 avian influenza virus is a novel re-assortment from at least four different strains of virus. Neuraminidase, which is a glycoprotein on the surface membrane, has been the target for drug treatment. However, some H7N9 strains that have been isolated from patient after drug treatment have a R292K mutation in neuraminidase. This substitution was found to facilitate drug resistance using protein- and virus- assays, in particular it gave a high resistance to the most commonly used drug, oseltamivir. The aim of this research is to understand the source of oseltamivir resistance using MD simulations and the MM/PB(GB)SA binding free energy approaches. Both methods can predict the reduced susceptibility of oseltamivir in good agreement to the IC ₅₀ binding energy, although MM/GBSA underestimates this prediction compared to the MM/PBSA calculation. Electrostatic interaction is the main contribution for oseltamivir binding in terms of both interaction and solvation. We found that the source of the drug resistance is a decrease in the binding interaction combined with the reduction of the dehydration penalty. The smaller K292 mutated residue has a larger binding pocket cavity compared to the wild-type resulting in the loss of drug carboxylate-K292 hydrogen bonding and an increased accessibility for water molecules around the K292 mutated residue. In addition, oseltamivir does not bind well to the R292K mutant complex as shown by the high degree of fluctuation in ligand RMSD during the simulation and the change in angular distribution of bulky side chain groups.

Investigating the mutation resistance of nonnucleoside inhibitors of HIV-RT using multiple microsecond atomistic simulations

Inhibiting HIV reverse transcriptase through the use of nonnucleoside reverse transcriptase inhibitors (NNRTIs) has become an essential component in drug regimens for the treatment of HIV. Older NNRTIs, such as nevirapine, are structurally rigid, exhibiting decreased inhibitory function on development of common mutations in the NNRTI-binding pocket, which is located around 10 Å from the catalytically active binding site. The newer generation of drugs, such as rilpivirine, are more flexible and resistant to binding pocket mutations but the mechanism by which they actually inhibit protein function and avoid mutations is not well-understood. To this end, we have performed 2-2.4 µs simulations with explicit solvent in an isobaric-isothermal ensemble of six different systems: apo wild-type, apo K103N/Y181C mutant, nevirapine-bound wild-type, nevirapine-bound mutant, rilpivirine-bound wild type, and rilpivirine-bound mutant. Analysis of protein conformations, principal components of motion, and mutual information between residues points to an inhibitory mechanism in which the primer grip stretches away from the catalytic triad of aspartic acids necessary for polymerization of HIV-encoding DNA, but is still unable to reveal a specific structural mechanism behind mutation resistance.

Combined quantum mechanics/molecular mechanics (QM/MM) methods in computational enzymology

Computational enzymology is a rapidly maturing field that is increasingly integral to understanding mechanisms of enzyme-catalyzed reactions and their practical applications. Combined quantum mechanics/molecular mechanics (QM/MM) methods are important in this field. By treating the reacting species with a quantum mechanical method (i.e., a method that calculates the electronic structure of the active site) and including the enzyme environment with simpler molecular mechanical methods, enzyme reactions can be modeled. Here, we review QM/MM methods and their application to enzyme-catalyzed reactions to investigate fundamental and practical problems in enzymology. A range of QM/MM methods is available, from cheaper and more approximate methods, which can be used for molecular dynamics simulations, to highly accurate electronic structure methods. We discuss how modeling of reactions using such methods can provide detailed insight into enzyme mechanisms and illustrate this by reviewing some recent applications. We outline some practical considerations for such simulations. Further, we highlight applications that show how QM/MM methods can contribute to the practical development and application of enzymology, e.g., in the interpretation and prediction of the effects of mutagenesis and in drug and catalyst design.

The structural role of Mg2+ ions in a class I RNA polymerase ribozyme: a molecular simulation study

According to the RNA world hypothesis, self-replicating ribozymes, storing the genetic information and being able to perform catalysis, were the constituents of the first living organisms. In particular, RNA polymerase ribozymes, similar to current proteinaceous enzymatic polymerases, may have been able to promote the synthesis of RNA strands in a primitive world. Polymerase catalysis is usually assisted by Mg(2+) ions, but it is not always trivial to find out experimentally the number of Mg(2+) ions placed in the active site as well as the identity and the number of their coordination ligands. Here, we addressed this issue in an artificial class I ligase ribozyme. On the basis of a recently solved crystal structure, we constructed computational models of reactant and product states of this ribozyme, considering monometallic and bimetallic species. Our models were relaxed by force field based molecular dynamics (MD) simulations and mixed quantum-classical (QM/MM) MD. The structural and dynamical properties of these models were consistent with experimental data and were validated by a comparison with the catalytic sites of proteinaceous DNA and RNA polymerases. Consistently with enzymatic polymerases, our results suggest that class I RNA ligases most probably contain two magnesium ions in the active site and they may, therefore, catalyze the junction of two RNA strands via "a two Mg(2+) ions" mechanism.

Correlation between electron localization and metal ion mutagenicity in DNA synthesis from QM/MM calculations

DNA polymerases require two divalent metal ions in the active site for catalysis. Mg(2+) has been confirmed to be the most probable cation utilized by most polymerases in vivo. Other metal ions are either potent mutagens or inhibitors. We used structural and topological analyses based on ab initio QM/MM calculations to study human DNA polymerase λ (Polλ) with different metals in the active site. Our results indicate a slightly longer O3'-Pα distance (∼3.6 Å) for most inhibitor cations compared to the natural and mutagenic metals (∼3.3-3.4 Å). Optimization with a larger basis set for the previously reported transition state (TS) structures (Cisneros et al., DNA Repair, 2008, 7, 1824.) gives barriers of 17.4 kcal mol(-1) and 15.1 kcal mol(-1) for the Mg(2+) and Mn(2+) catalyzed reactions respectively. Relying on the key relation between the topological signature of a metal cation and its selectivity within biological systems (de Courcy et al., J. Chem. Theor. Comput., 2010, 6, 1048.) we have performed electron localization function (ELF) topological analyses. These analyses show that all inhibitor and mutagenic metals considered, except Na(+), present a "split" of the outer-shell density of the metal. This "splitting" is not observed for the non-mutagenic Mg(2+) metal. Population and multipole analyses on the ELF basins reveal that the electronic dipolar and quadrupolar polarization is significantly different with Mg(2+) compared to all other cations. Our results shed light at the atomic level on the subtle differences between Mg(2+), mutagenic, and inhibitor metals in DNA polymerases. These results provide a correlation between the electronic distribution of the cations in the active site and the possible consequences on DNA synthesis.

Increased dispersion and solubility of carbon nanotubes noncovalently modified by the polysaccharide biopolymer, chitosan: MD simulations

In order to explain the solubility of carbon nanotubes (CNT), including single walled CNTs, wrapped with chitosan of a 60% degree of deacetylation, MD simulations were applied to represent three chitosan concentrations, using two pristine CNTs (pCNT-pCNT), and one and two CNTs wrapped (pCNT-cwCNT and cwCNT-cwCNT). The CNT aggregation was observed in pCNT-pCNT and pCNT-cwCNT due to van der Waals interactions between tube-tube aromatic rings, and inter-CNT bridging by chitosan, respectively. At higher chitosan concentrations, such that most to all of CNTs were wrapped with chitosan, charge-charge repulsion was found to separate robustly the cwCNTs and lead to a well dispersed solution.

Ring puckering: a metric for evaluating the accuracy of AM1, PM3, PM3CARB-1, and SCC-DFTB carbohydrate QM/MM simulations

The puckered conformations of furanose and pyranose carbohydrate rings are central to analyzing the action of enzymes on carbohydrates. Enzyme reaction mechanisms are generally inaccessible to experiments and so have become the focus of QM(semiempirical)/MM simulations. We show that the complete free energy of puckering is required to evaluate the accuracy of semiempirical methods used to study reactions involving carbohydrates. Interestingly, we find that reducing the free energy space to lower dimensions results in near meaningless minimum energy pathways. We analyze the furanose and pyranose free energy pucker surfaces and volumes using AM1, PM3, PM3CARB-1, and SCC-DFTB. A comparison with DFT optimized structures and a HF free energy surface reveals that SCC-DFTB provides the best semiempirical description of five- and six-membered carbohydrate ring deformation.

Combined QM/MM mechanistic study of the acylation process in furin complexed with the H5N1 avian influenza virus hemagglutinin's cleavage site

Combined quantum mechanical/molecular mechanical (QM/MM) techniques have been applied to investigate the detailed reaction mechanism of the first step of the acylation process by furin in which the cleavage site of the highly pathogenic avian influenza virus subtype H5N1 (HPH5) acts as its substrate. The energy profile shows a simultaneous mechanism, known as a concerted reaction, of the two subprocesses: the proton transfer from Ser368 to His194 and the nucleophilic attack on the carbonyl carbon of the scissile peptide of the HPH5 cleavage site with a formation of tetrahedral intermediate (INT). The calculated energy barrier for this reaction is 16.2 kcal.mol(-1) at QM/MM B3LYP/6-31+G*//PM3-CHARMM22 level of theory. Once the reaction proceeds, the ordering of the electrostatic stabilization by protein environment is of the enzyme-substrate < transition state < INT complexes. Asp153 was found to play the most important role in the enzymatic reaction by providing the highest degree of intermediate complex stabilization. In addition, the negatively charged carbonyl oxygen of INT is well stabilized by the oxyanion hole constructed by Asn295's carboxamide and Ser368's backbone.

Hydrolysis of aspartic acid phosphoramidate nucleotides: a comparative quantum chemical study

L-Aspartic acid has recently been found to be a good leaving group during HIV reverse transcriptase catalyzed incorporation of deoxyadenosine monophosphate (dAMP) in DNA. This showed that L-Asp is a good mimic for the pyrophosphate moiety of deoxyadenosine triphosphate. The present work explores the thermochemistry and mechanism for hydrolysis of several models for L-aspartic-dAMP using B3LYP/DGDZVP, MP2/6-311++G** and G3MP2 level of theory. The effect of the new compound is gradually investigated: starting from a simple methyl amine leaving group up to the aspartic acid leaving group. The enzymatic environment was mimicked by involving two Mg(2+) ions and some important active site residues in the reaction. All reactions are compared to the corresponding O-coupled leaving group, which is methanol for methyl amine and malic acid for aspartic acid. With methyl amine as a leaving group a tautomeric associative or tautomeric dissociative mechanism is preferred and the barrier is lower than the comparable mechanism with methanol as a leaving group. The calculations on the aspartic acid in the enzymatic environment show that qualitatively the mechanism is the same as for triphosphate but the barrier for hydrolysis by the associative mechanism is higher for L-aspartic-dAMP than for L-malic-dAMP and pyrophosphate.

QM/MM methods for biomolecular systems

Combined quantum-mechanics/molecular-mechanics (QM/MM) approaches have become the method of choice for modeling reactions in biomolecular systems. Quantum-mechanical (QM) methods are required for describing chemical reactions and other electronic processes, such as charge transfer or electronic excitation. However, QM methods are restricted to systems of up to a few hundred atoms. However, the size and conformational complexity of biopolymers calls for methods capable of treating up to several 100,000 atoms and allowing for simulations over time scales of tens of nanoseconds. This is achieved by highly efficient, force-field-based molecular mechanics (MM) methods. Thus to model large biomolecules the logical approach is to combine the two techniques and to use a QM method for the chemically active region (e.g., substrates and co-factors in an enzymatic reaction) and an MM treatment for the surroundings (e.g., protein and solvent). The resulting schemes are commonly referred to as combined or hybrid QM/MM methods. They enable the modeling of reactive biomolecular systems at a reasonable computational effort while providing the necessary accuracy.

Catalytic mechanism of human DNA polymerase lambda with Mg2+ and Mn2+ from ab initio quantum mechanical/molecular mechanical studies

DNA polymerases play a crucial role in the cell cycle due to their involvement in genome replication and repair. Understanding the reaction mechanism by which these polymerases carry out their function can provide insights into these processes. Recently, the crystal structures of human DNA polymerase lambda (Pollambda) have been reported both for pre- and post-catalytic complexes [García-Díaz et al., DNA Repair 3 (2007), 1333]. Here we employ the pre-catalytic complex as a starting structure for the determination of the catalytic mechanism of Pollambda using ab initio quantum mechanical/molecular mechanical methods. The reaction path has been calculated using Mg(2+) and Mn(2+) as the catalytic metals. In both cases the reaction proceeds through a two-step mechanism where the 3'-OH of the primer sugar ring is deprotonated by one of the conserved Asp residues (D490) in the active site before the incorporation of the nucleotide to the nascent DNA chain. A significant charge transfer is observed between both metals and some residues in the active site as the reaction proceeds. The optimized reactant and product structures agree with the reported crystal structures. In addition, the calculated reaction barriers for both metals are close to experimentally estimated barriers. Energy decomposition analysis to explain individual residue contributions suggests that several amino acids surrounding the active site are important for catalysis. Some of these residues, including R420, R488 and E529, have been implicated in catalysis by previous mutagenesis experiments on the homologous residues on Polbeta. Furthermore, Pollambda residues R420 and E529 found to be important from the energy decomposition analysis, are homologous to residues R183 and E295 in Polbeta, both of which are linked to cancer. In addition, residues R386, E391, K422 and K472 appear to have an important role in catalysis and could be a potential target for mutagenesis experiments. There is partial conservation of these residues across the Pol X family of DNA polymerases.

Computational enzymology: modelling the mechanisms of biological catalysts

Simulations and modelling [e.g. with combined QM/MM (quantum mechanics/molecular mechanics) methods] are increasingly important in investigations of enzyme-catalysed reaction mechanisms. Calculations offer the potential of uniquely detailed, atomic-level insight into the fundamental processes of biological catalysis. Highly accurate methods promise quantitative comparison with experiments, and reliable predictions of mechanisms, revolutionizing enzymology.

Hybrid quantum mechanical/molecular mechanical molecular dynamics simulations of HIV-1 integrase/inhibitor complexes

Human immunodeficiency virus (HIV)-1 integrase (IN) is an attractive target for development of acquired immunodeficiency syndrome chemotherapy. In this study, conventional and coupled quantum mechanical and molecular mechanical (QM/MM) molecular dynamics (MD) simulations of HIV-1 IN complexed with 5CITEP (IN-5CITEP) were carried out. In addition to differences in the bound position of 5CITEP, significant differences at the two levels of theory were observed in the metal coordination geometry and the areas involving residues 116-119 and 140-166. In the conventional MD simulation, the coordination of Mg(2+) was found to be a near-perfect octahedral geometry whereas a distorted octahedral complex was observed in QM/MM. All of the above reasons lead to a different pattern of protein-ligand salt link formation that was not observed in the classical MD simulation. Furthermore to provide a theoretical understanding of inhibition mechanisms of 5CITEP and its derivative (DKA), hybrid QM/MM MD simulations of the two complexes (IN-5CITEP and IN-DKA) have been performed. The results reveal that areas involving residues 60-68, 116-119, and 140-149 were substantially different among the two systems. The two systems show similar pattern of metal coordination geometry, i.e., a distorted octahedron. In IN-DKA, both OD1 and OD2 of Asp-64 coordinate the Mg(2+) in a monodentate fashion whereas only OD1 is chelated to the metal as observed in IN-5CITEP. The high potency of DKA as compared to 5CITEP is supported by a strong salt link formed between its carboxylate moiety and the ammonium group of Lys-159. Detailed comparisons between HIV-1 IN complexed with DKA and with 5CITEP provide information about ligand structure effects on protein-ligand interactions in particular with the Lys-159. This is useful for the design of new selective HIV-1 IN inhibitors.