Magnesium-cationic dummy atom molecules enhance representation of DNA polymerase beta in molecular dynamics simulations: improved accuracy in studies of structural features and mutational effects

Human Polβ Natural Polymorphic Variants G118V and R149I Affects Substate Binding and Catalysis

DNA polymerase β (Polβ) expression is essential for the cell's response to DNA damage that occurs during natural cellular processes. Polβ is considered the main reparative DNA polymerase, whose role is to fill the DNA gaps arising in the base excision repair pathway. Mutations in Polβ can lead to cancer, neurodegenerative diseases, or premature aging. Many single-nucleotide polymorphisms have been identified in the POLB gene, but the consequences of these polymorphisms are not always clear. It is known that some polymorphic variants in the Polβ sequence reduce the efficiency of DNA repair, thereby raising the frequency of mutations in the genome. In the current work, we studied two polymorphic variants (G118V and R149I separately) of human Polβ that affect its DNA-binding region. It was found that each amino acid substitution alters Polβ's affinity for gapped DNA. Each polymorphic variant also weakens its binding affinity for dATP. The G118V variant was found to greatly affect Polβ's ability to fill gapped DNA and slowed the catalytic rate as compared to the wild-type enzyme. Thus, these polymorphic variants seem to decrease the ability of Polβ to maintain base excision repair efficiency.

Mechanism of Calcium Permeation in a Glutamate Receptor Ion Channel

The α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) are neurotransmitter-activated cation channels ubiquitously expressed in vertebrate brains. The regulation of calcium flux through the channel pore by RNA-editing is linked to synaptic plasticity while excessive calcium influx poses a risk for neurodegeneration. Unfortunately, the molecular mechanisms underlying this key process are mostly unknown. Here, we investigated calcium conduction in calcium-permeable AMPAR using Molecular Dynamics (MD) simulations with recently introduced multisite force-field parameters for Ca²⁺. Our calculations are consistent with experiment and explain the distinct calcium permeability in different RNA-edited forms of GluA2. For one of the identified metal binding sites, multiscale Quantum Mechanics/Molecular Mechanics (QM/MM) simulations further validated the results from MD and revealed small but reproducible charge transfer between the metal ion and its first solvation shell. In addition, the ion occupancy derived from MD simulations independently reproduced the Ca²⁺ binding profile in an X-ray structure of an NaK channel mimicking the AMPAR selectivity filter. This integrated study comprising X-ray crystallography, multisite MD, and multiscale QM/MM simulations provides unprecedented insights into Ca²⁺ permeation mechanisms in AMPARs, and paves the way for studying other biological processes in which Ca²⁺ plays a pivotal role.

A Polarizable Cationic Dummy Metal Ion Model

A novel locally polarizable multisite model based on the original cation dummy atom (CDA) model is described for molecular dynamics simulations of ions in condensed phases. Polarization effects are introduced by the electronegativity equalization model (EEM) method where charges on the metal ion and its dummy atoms can fluctuate to respond to the environment. This model includes explicit polarization and ion-induced interactions and can be coupled with nonpolarizable or polarizable water models, making it more transferable to simpler force fields. This approach allows us to enhance the original fixed charge CDA model where the charge distribution cannot adapt to the local solvent structure. To illustrate the new CDApol model, we examined properties of the Zn2+, Al3+, and Zr4+ ions in aqueous solution. The polarizable model and Lennard-Jones parameters were refined for octahedrally coordinated Zn2+, Al3+, and Zr4+ CDAs to reproduce thermodynamic and geometrical properties. Using this locally polarizable model, we were able to obtain the experimental hydration free energy, ion-oxygen distance, and coordination number coupled with the standard 12-6 Lennard-Jones model. This model can be used in myriad additional applications where local polarization and charge transfer is important.

Computational Modeling of Supramolecular Metallo-organic Cages-Challenges and Opportunities

Self-assembled metallo-organic cages have emerged as promising biomimetic platforms that can encapsulate whole substrates akin to an enzyme active site. Extensive experimental work has enabled access to a variety of structures, with a few notable examples showing catalytic behavior. However, computational investigations of metallo-organic cages are scarce, not least due to the challenges associated with their modeling and the lack of accurate and efficient protocols to evaluate these systems. In this review, we discuss key molecular principles governing the design of functional metallo-organic cages, from the assembly of building blocks through binding and catalysis. For each of these processes, computational protocols will be reviewed, considering their inherent strengths and weaknesses. We will demonstrate that while each approach may have its own specific pitfalls, they can be a powerful tool for rationalizing experimental observables and to guide synthetic efforts. To illustrate this point, we present several examples where modeling has helped to elucidate fundamental principles behind molecular recognition and reactivity. We highlight the importance of combining computational and experimental efforts to speed up supramolecular catalyst design while reducing time and resources.

"Bucket brigade" using lysine residues in RNA-dependent RNA polymerase of SARS-CoV-2

The RNA-dependent RNA polymerase (RdRp) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a promising drug target for coronavirus disease 2019 (COVID-19) because it plays the most important role in the replication of the RNA genome. Nucleotide analogs such as remdesivir and favipiravir are thought to interfere with the RNA replication by RdRp. More specifically, they are expected to compete with nucleoside triphosphates, such as ATP. However, the process in which these drug molecules and nucleoside triphosphates are taken up by RdRp remains unknown. In this study, we performed all-atom molecular dynamics simulations to clarify the recognition mechanism of RdRp for these drug molecules and ATP that were at a distance. The ligand recognition ability of RdRp decreased in the order of remdesivir, favipiravir, and ATP. We also identified six recognition paths. Three of them were commonly found in all ligands, and the remaining three paths were ligand-dependent ones. In the common two paths, it was observed that the multiple lysine residues of RdRp carried the ligands to the binding site like a "bucket brigade." In the remaining common path, the ligands directly reached the binding site. Our findings contribute to the understanding of the efficient ligand recognition by RdRp at the atomic level.

Developing and Assessing Nonbonded Dummy Models of Magnesium Ion with Different Hydration Free Energy References

A large diversity in the targeted hydration free energies (HFEs) during model parameterization of metal ions was reported in the literature with a difference by dozens of kcal/mol. Here, we developed a series of nonbonded dummy models of the Mg²⁺ ion targeting different HFE references in TIP3P water, followed by assessments of the designed models in the simulations of MgCl₂ solution and biological systems. Together with the comparison of existing models, we conclude that the difference in the targeted HFEs has a limited influence on the model performance, while the usability of these models differs from case to case. The feasibility of reproducing more properties of Mg²⁺ such as diffusion constants and water exchange rates using a nonbonded dummy model is demonstrated. Underestimated activity derivative and osmotic coefficient of MgCl₂ solutions in high concentration reveal a necessity for further optimization of ion-pair interactions. The developed dummy models are applicable to metal coordination with Asp, Glu, and His residues in metalloenzymes, and the performance in predicting monodentate or bidentate binding modes of Asp/Glu residues depends on the complexity of metal centers and the choice of protein force fields. When both the binding modes coexist, the nonbonded dummy models outperform point charge models, probably in need of considering polarization of metal-binding residues by, for instance, charge calibration in classical force fields. This work is valuable for the use and further development of magnesium ion models for simulations of metal-containing systems with good accuracy.

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.

Determining the atomic charge of calcium ion requires the information of its coordination geometry in an EF-hand motif

It is challenging to parameterize the force field for calcium ions (Ca²⁺) in calcium-binding proteins because of their unique coordination chemistry that involves the surrounding atoms required for stability. In this work, we observed a wide variation in Ca²⁺ binding loop conformations of the Ca²⁺-binding protein calmodulin, which adopts the most populated ternary structures determined from the molecular dynamics simulations, followed by ab initio quantum mechanical (QM) calculations on all 12 amino acids in the loop that coordinate Ca²⁺ in aqueous solution. Ca²⁺ charges were derived by fitting to the electrostatic potential in the context of a classical or polarizable force field (PFF). We discovered that the atomic radius of Ca²⁺ in conventional force fields is too large for the QM calculation to capture the variation in the coordination geometry of Ca²⁺ in its ionic form, leading to unphysical charges. Specifically, we found that the fitted atomic charges of Ca²⁺ in the context of PFF depend on the coordinating geometry of electronegative atoms from the amino acids in the loop. Although nearby water molecules do not influence the atomic charge of Ca²⁺, they are crucial for compensating for the coordination of Ca²⁺ due to the conformational flexibility in the EF-hand loop. Our method advances the development of force fields for metal ions and protein binding sites in dynamic environments.

Evaluating the Performance of a Non-Bonded Cu2+ Model Including Jahn-Teller Effect into the Binding of Tyrosinase Inhibitors

Tyrosinase (TYR) is a metalloenzyme classified as a type-3 copper protein, which is involved in the synthesis of melanin through a catalytic process beginning with the conversion of the amino acid l-Tyrosine (l-Tyr) to l-3,4-dihydroxyphenylalanine (l-DOPA). It plays an important role in the mechanism of melanogenesis in various organisms including mammals, plants, and fungi. Herein, we used a combination of computational molecular modeling techniques including molecular dynamic (MD) simulations and the linear interaction energy (LIE) model to evaluate the binding free energy of a set of analogs of kojic acid (KA) in complex with TYR. For the MD simulations, we used a dummy model including the description of the Jahn–Teller effect for Cu²⁺ ions in the active site of this enzyme. Our results show that the LIE model predicts the TYR binding affinities of the inhibitor in close agreement to experimental results. Overall, we demonstrate that the classical model provides a suitable description of the main interactions between analogs of KA and Cu²⁺ ions in the active site of TYR.

Novel phosphorylation-dependent regulation in an unstructured protein

This work explores how phosphorylation of an unstructured protein region in inhibitor-2 (I2) regulates protein phosphatase-1 (PP1) enzyme activity using molecular dynamics (MD). Free I2 is largely unstructured; however, when bound to PP1, three segments adopt a stable structure. In particular, an I2 helix (i-helix) blocks the PP1 active site and inhibits phosphatase activity. I2 phosphorylation in the PP1-I2 complex activates phosphatase activity without I2 dissociation. The I2 Thr74 regulatory phosphorylation site is in an unstructured domain in PP1-I2. PP1-I2 MD demonstrated that I2 phosphorylation promotes early steps of PP1-I2 activation in explicit solvent models. Moreover, phosphorylation-dependent activation occurred in PP1-I2 complexes derived from I2 orthologs with diverse sequences from human, yeast, worm, and protozoa. This system allowed exploration of features of the 73-residue unstructured human I2 domain critical for phosphorylation-dependent activation. These studies revealed that components of I2 unstructured domain are strategically positioned for phosphorylation responsiveness including a transient α-helix. There was no evidence that electrostatic interactions of I2 phosphothreonine74 influenced PP1-I2 activation. Instead, phosphorylation altered the conformation of residues around Thr74. Phosphorylation uncurled the distance between I2 residues Glu71 to Tyr76 to promote PP1-I2 activation, whereas reduced distances reduced activation. This I2 residue Glu71 to Tyr76 distance distribution, independently from Thr74 phosphorylation, controls I2 i-helix displacement from the PP1 active site leading to PP1-I2 activation.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Assessing the Performance of the Nonbonded Mg2+ Models in a Two-Metal-Dependent Ribonuclease

Magnesium ions (Mg²⁺), abundant in living cells, are essential for biomolecular structure, dynamics, and function. The biological importance of Mg²⁺ has motivated continuous development and improvement of various Mg²⁺ models for molecular dynamics (MD) simulations during the last decades. There are four types of nonbonded Mg²⁺ models: the point charge models based on a 12-6 or 12-6-4 type Lennard-Jones (LJ) potential, and the multisite models based on a 12-6 or 12-6-4 LJ potential. Here, we systematically assessed the performance of these four types of nonbonded Mg²⁺ models (21 models in total) in terms of maintaining a challenging intermediate state configuration captured in the structure of a prototypical two-metal-ion RNase H complex with an RNA/DNA hybrid. Our data demonstrate that the 12-6-4 multisite models, which account for charge-induced dipole interactions, perform the best in reproducing all the unique coordination modes in this intermediate state and maintaining the correct carboxylate denticity. Our benchmark work provides a useful guideline for MD simulations and structural refinement of Mg²⁺-containing biomolecular systems.

A Local Allosteric Network in Heat Shock Protein 70 (Hsp70) Links Inhibitor Binding to Enzyme Activity and Distal Protein-Protein Interactions

Allosteric inhibitors can be more difficult to optimize without an understanding of how their binding influences the conformational motions of the target. Here, we used an integrated computational and experimental approach to probe the molecular mechanism of an allosteric inhibitor of heat shock protein 70 (Hsp70). The anticancer compound, MKT-077, is known to bind a conserved site in members of the Hsp70 family, which favors the ADP-bound state and interferes with a protein-protein interaction (PPI) at long range. However, the binding site does not overlap with either the nucleotide-binding cleft or the PPI contact surface, so its mechanism is unclear. To this end, we modeled Hsp70's internal dynamics and studied how MKT-077 alters local sampling of its allosteric states. The results pointed to a set of concerted motions between five loops in Hsp70's nucleotide-binding domain (NBD), surrounding the MKT-077 binding site. To test this prediction, we mutated key residues and monitored chaperone activities in vitro. Together, the results indicate that MKT-077 interacts with loop222 to favor a pseudo-ADP bound conformer of Hsp70's NBD, even when ATP is present. We used this knowledge to synthesize an analog of MKT-077 that would better prevent motions of loop222 and confirmed that it had improved antiproliferative activity in breast cancer cells. These results provide an example of how to unlock and leverage the complex mechanisms of allosteric inhibitors.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Extending the Nonbonded Cationic Dummy Model to Account for Ion-Induced Dipole Interactions

Modeling metalloproteins often requires classical molecular dynamics (MD) simulations in order to capture their relevant motions, which in turn necessitates reliable descriptions of the metal centers involved. One of the most successful approaches to date is provided by the "cationic dummy model", where the positive charge of the metal ion is transferred toward dummy particles that are bonded to the central metal ion in a predefined coordination geometry. While this approach allows for ligand exchange, and captures the correct electrostatics as demonstrated for different divalent metal ions, current dummy models neglect ion-induced dipole interactions. In the present work, we resolve this weakness by taking advantage of the recently introduced 12-6-4 type Lennard-Jones potential to include ion-induced dipole interactions. We revise our previous dummy model for Mg²⁺ and demonstrate that the resulting model can simultaneously reproduce the experimental solvation free energy and metal-ligand distances without the need for artificial restraints or bonds. As ion-induced dipole interactions become particularly important for highly charged metal ions, we develop dummy models for the biologically relevant ions Al³⁺, Fe³⁺, and Cr³⁺. Finally, the effectiveness of our new models is demonstrated in MD simulations of several diverse (and highly challenging to simulate) metalloproteins.

Molecular simulations of self-assembly processes in metal-organic frameworks: Model dependence

Molecular simulation is a powerful tool for investigating microscopic behavior in various chemical systems, where the use of suitable models is critical to successfully reproduce the structural and dynamic properties of the real systems of interest. In this context, molecular dynamics simulation studies of self-assembly processes in metal-organic frameworks (MOFs), a well-known class of porous materials with interesting chemical and physical properties, are relatively challenging, where a reasonably accurate representation of metal-ligand interactions is anticipated to play an important role. In the current study, we both investigate the performance of some existing models and introduce and test new models to help explore the self-assembly in an archetypal Zn-carboxylate MOF system. To this end, the behavior of six different Zn-ion models, three solvent models, and two ligand models was examined and validated against key experimental structural parameters. To explore longer time scale ordering events during MOF self-assembly via explicit solvent simulations, it is necessary to identify a suitable combination of simplified model components representing metal ions, organic ligands, and solvent molecules. It was observed that an extended cationic dummy atom (ECDA) Zn-ion model combined with an all-atom carboxylate ligand model and a simple dipolar solvent model can reproduce characteristic experimental structures for the archetypal MOF system. The successful use of these models in extensive sets of molecular simulations, which provide key insights into the self-assembly mechanism of this archetypal MOF system occurring during the early stages of this process, has been very recently reported.

Development of a multisite model for Ni(II) ion in solution from thermodynamic and kinetic data

Force-field parameters are developed for a multisite model of Ni(II) ions to be used in molecular dynamics simulations combined to enhanced sampling methods. The performances of two charge-partitioning schemes are validated by taking into account structural, thermodynamic, and kinetic observables. One of the two models, featuring partial charges on the dummy atoms only, matches both Ni(II) free energy of solvation and water exchange rates. Such model is particularly suited to study complexation events at a fully dynamic description. © 2017 Wiley Periodicals, Inc.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Functional roles of magnesium binding to extracellular signal-regulated kinase 2 explored by molecular dynamics simulations and principal component analysis

Molecular dynamics (MD) simulations coupled with principal component (PC) analysis were carried out to study functional roles of Mg²⁺ binding to extracellular signal-regulated kinase 2 (ERK2). The results suggest that Mg²⁺ binding heavily decreases eigenvalue of the first principal component and totally inhibits motion strength of ERK2, which favors stabilization of ERK2 structure. Binding free energy predictions indicate that Mg²⁺ binding produces an important effect on binding ability of adenosine triphosphate (ATP) to ERK2 and strengthens the ATP binding. The calculations of residue-based free energy decomposition show that lack of Mg²⁺ weakens interactions between the hydrophobic rings of ATP and five residues I29, V37, A50, L105, and L154. Hydrogen bond analyses also prove that Mg²⁺ binding increases occupancies of hydrogen bonds formed between ATP and residues K52, Q103, D104, and M106. We expect that this study can provide a significant theoretical hint for designs of anticancer drugs targeting ERK2.

Development of Site-Specific Mg(2+)-RNA Force Field Parameters: A Dream or Reality? Guidelines from Combined Molecular Dynamics and Quantum Mechanics Simulations

The vital contribution of Mg²⁺ ions to RNA biology is challenging to dissect at the experimental level. This calls for the integrative support of atomistic simulations, which at the classical level are plagued by limited accuracy. Indeed, force fields intrinsically neglect nontrivial electronic effects that Mg²⁺ exerts on its surrounding ligands in varying RNA coordination environments. Here, we present a combined computational study based on classical molecular dynamics (MD) and Density Functional Theory (DFT) calculations, aimed at characterizing (i) the performance of five Mg²⁺ force field (FF) models in RNA systems and (ii) how charge transfer and polarization affect the binding of Mg²⁺ ions in different coordination motifs. As a result, a total of ∼2.5 μs MD simulations (100/200 ns for each run) for two prototypical Mg²⁺-dependent ribozymes showed remarkable differences in terms of populations of inner-sphere coordination site types. Most importantly, complementary DFT calculations unveiled that differences in charge transfer and polarization among recurrent Mg²⁺-RNA coordination motifs are surprisingly small. In particular, the charge of the Mg²⁺ ions substantially remains constant through different coordination sites, suggesting that the common philosophy of developing site-specific Mg²⁺ ion parameters is not in line with the physical origin of the Mg²⁺-RNA MD simulations inaccuracies. Overall, this study constitutes a guideline for an adept use of current Mg²⁺ models and provides novel insights for the rational development of next-generation Mg²⁺ FFs to be employed for atomistic simulations of RNA.

Clarifying binding difference of ATP and ADP to extracellular signal-regulated kinase 2 by using molecular dynamics simulations

Extracellular signal-regulated kinase 2 is a promising target for designs and development of anticancer drugs. Molecular dynamics simulations and molecular mechanics Poisson-Boltzmann method were applied to study binding difference of ADP and ATP to extracellular signal-regulated kinase 2. The results prove that the binding ability of ATP to extracellular signal-regulated kinase 2 is stronger than that of ADP. Principal component analysis performed by using molecular dynamics trajectories suggests that binding of ADP and ATP to extracellular signal-regulated kinase 2 change motion directions of two helices α1 and α2. Residue-based free energy decomposition method was adopted to calculate contributions of separate residues to associations of ADP and ATP with extracellular signal-regulated kinase 2. The results show that ADP and ATP produce strong CH-π interactions with five residues Ile29, Val37, Ala50, Leu105, and Leu154. In addition, five hydrogen bonding interactions of ADP and ATP with residues Lys52, Gln103, Asp104, and Met106 also stabilize bindings of ADP and ATP to extracellular signal-regulated kinase 2. Overall, the CH-π interactions of ATP with five residues Ile29, Val37, Ala50, Leu105, and Leu154 are stronger than ADP. This study is expected to contribute a significant theoretical hint for designs of anticancer drugs targeting extracellular signal-regulated kinase 2.

The control of the discrimination between dNTP and rNTP in DNA and RNA polymerase

Understanding the origin of discrimination between rNTP and dNTP by DNA/RNA polymerases is important both for gaining fundamental knowledge on the corresponding systems and for advancing the design of specific drugs. This work explores the nature of this discrimination by systematic calculations of the transition state (TS) binding energy in RB69 DNA polymerase (gp43) and T7 RNA polymerase. The calculations reproduce the observed trend, in particular when they included the water contribution obtained by the water flooding approach. Our detailed study confirms the idea that the discrimination is due to the steric interaction between the 2'OH and Tyr416 in DNA polymerase, while the electrostatic interaction is the source of the discrimination in RNA polymerase. Proteins 2016; 84:1616-1624. © 2016 Wiley Periodicals, Inc.

Divalent Ion Dependent Conformational Changes in an RNA Stem-Loop Observed by Molecular Dynamics

We compare the performance of five magnesium (Mg²⁺) ion models in simulations of an RNA stem loop which has an experimentally determined divalent ion dependent conformational shift. We show that despite their differences in parametrization and resulting van der Waals terms, including differences in the functional form of the nonbonded potential, when the RNA adopts its folded conformation, all models behave similarly across ten independent microsecond length simulations with each ion model. However, when the entire structure ensemble is accounted for, chelation of Mg²⁺ to RNA is seen in three of the five models, most egregiously and likely artifactual in simulations using a 12-6-4 model for the Lennard-Jones potential. Despite the simple nature of the fixed point-charge and van der Waals sphere models employed, and with the exception of the likely oversampled directed chelation of the 12-6-4 potential models, RNA–Mg²⁺ interactions via first shell water molecules are surprisingly well described by modern parameters, allowing us to observe the spontaneous conformational shift from Mg²⁺ free RNA to Mg²⁺ associated RNA structure in unrestrained molecular dynamics simulations.

Structural Studies of Trypanosoma brucei RNA Editing Ligases and Their Binding Partner Proteins

To study the mechanism of ligating nicked RNA strands, we conducted molecular dynamics simulations of Trypanosoma brucei RNA editing ligases L1 and L2 complexed with double-stranded RNA (dsRNA) fragments. In each resulting model, a Mg(2+) ion coordinates the 5'-PO4 of the nicked nucleotide and the 3'-OH of the terminal nucleotide for a nucleophilic reaction consistent with the postulated step 3 chemistry of the ligation mechanism. Moreover, coordination of the 3'-OH to the Mg(2+) ion may lower its pKa, thereby rendering it a more effective nucleophile as an oxyanion. Thus, Mg(2+) may play a twofold role: bringing the reactants into the proximity of each other and activating the nucleophile. We also conducted solvated interaction energy calculations to explore whether ligation specificities can be correlated to ligase-dsRNA binding affinity changes. The calculated dsRNA binding affinities are stronger for both L1 and L2 when the terminal nucleotide is changed from cytosine to guanine, in line with their experimentally measured ligation specificities. Because the ligation mechanism is also influenced by interactions of the ligase with partner proteins from the editosome subcomplex, we also modeled the structure of the RNA-bound L2 in complex with the oligonucleotide binding (OB) domain of largest editosome interacting protein A1. The resulting L2-dsRNA-A1 model, which is consistent with mutagenesis and binding data recorded to date, provides the first atomic-level glimpse of plausible interactions around the RNA ligation site in the presence of an OB domain presented in-trans to a nucleic acid ligase.

CHARMM Force-Fields with Modified Polyphosphate Parameters Allow Stable Simulation of the ATP-Bound Structure of Ca(2+)-ATPase

Adenosine triphosphate (ATP) is an indispensable energy source in cells. In a wide variety of biological phenomena like glycolysis, muscle contraction/relaxation, and active ion transport, chemical energy released from ATP hydrolysis is converted to mechanical forces to bring about large-scale conformational changes in proteins. Investigation of structure-function relationships in these proteins by molecular dynamics (MD) simulations requires modeling of ATP in solution and ATP bound to proteins with accurate force-field parameters. In this study, we derived new force-field parameters for the triphosphate moiety of ATP based on the high-precision quantum calculations of methyl triphosphate. We tested our new parameters on membrane-embedded sarcoplasmic reticulum Ca(2+)-ATPase and four soluble proteins. The ATP-bound structure of Ca(2+)-ATPase remains stable during MD simulations, contrary to the outcome in shorter simulations using original parameters. Similar results were obtained with the four ATP-bound soluble proteins. The new force-field parameters were also tested by investigating the range of conformations sampled during replica-exchange MD simulations of ATP in explicit water. Modified parameters allowed a much wider range of conformational sampling compared with the bias toward extended forms with original parameters. A diverse range of structures agrees with the broad distribution of ATP conformations in proteins deposited in the Protein Data Bank. These simulations suggest that the modified parameters will be useful in studies of ATP in solution and of the many ATP-utilizing proteins.

Development and Application of a Nonbonded Cu(2+) Model That Includes the Jahn-Teller Effect

Metal ions are both ubiquitous to and crucial in biology. In classical simulations, they are typically described as simple van der Waals spheres, making it difficult to provide reliable force field descriptions for them. An alternative is given by nonbonded dummy models, in which the central metal atom is surrounded by dummy particles that each carry a partial charge. While such dummy models already exist for other metal ions, none is available yet for Cu(2+) because of the challenge to reproduce the Jahn-Teller distortion. This challenge is addressed in the current study, where, for the first time, a dummy model including a Jahn-Teller effect is developed for Cu(2+). We successfully validate its usefulness by studying metal binding in two biological systems: the amyloid-β peptide and the mixed-metal enzyme superoxide dismutase. We believe that our parameters will be of significant value for the computational study of Cu(2+)-dependent biological systems using classical models.

Comparison of structural, thermodynamic, kinetic and mass transport properties of Mg(2+) ion models commonly used in biomolecular simulations

The prevalence of Mg(2+) ions in biology and their essential role in nucleic acid structure and function has motivated the development of various Mg(2+) ion models for use in molecular simulations. Currently, the most widely used models in biomolecular simulations represent a nonbonded metal ion as an ion-centered point charge surrounded by a nonelectrostatic pairwise potential that takes into account dispersion interactions and exchange effects that give rise to the ion's excluded volume. One strategy toward developing improved models for biomolecular simulations is to first identify a Mg(2+) model that is consistent with the simulation force fields that closely reproduces a range of properties in aqueous solution, and then, in a second step, balance the ion-water and ion-solute interactions by tuning parameters in a pairwise fashion where necessary. The present work addresses the first step in which we compare 17 different nonbonded single-site Mg(2+) ion models with respect to their ability to simultaneously reproduce structural, thermodynamic, kinetic and mass transport properties in aqueous solution. None of the models based on a 12-6 nonelectrostatic nonbonded potential was able to reproduce the experimental radial distribution function, solvation free energy, exchange barrier and diffusion constant. The models based on a 12-6-4 potential offered improvement, and one model in particular, in conjunction with the SPC/E water model, performed exceptionally well for all properties. The results reported here establish useful benchmark calculations for Mg(2+) ion models that provide insight into the origin of the behavior in aqueous solution, and may aid in the development of next-generation models that target specific binding sites in biomolecules.

Molecular dynamics study of the opening mechanism for DNA polymerase I

During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although X-ray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme:DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics.

QM/MM molecular dynamics studies of metal binding proteins

Mixed quantum-classical (quantum mechanical/molecular mechanical (QM/MM)) simulations have strongly contributed to providing insights into the understanding of several structural and mechanistic aspects of biological molecules. They played a particularly important role in metal binding proteins, where the electronic effects of transition metals have to be explicitly taken into account for the correct representation of the underlying biochemical process. In this review, after a brief description of the basic concepts of the QM/MM method, we provide an overview of its capabilities using selected examples taken from our work. Specifically, we will focus on heme peroxidases, metallo-β-lactamases, α-synuclein and ligase ribozymes to show how this approach is capable of describing the catalytic and/or structural role played by transition (Fe, Zn or Cu) and main group (Mg) metals. Applications will reveal how metal ions influence the formation and reduction of high redox intermediates in catalytic cycles and enhance drug metabolism, amyloidogenic aggregate formation and nucleic acid synthesis. In turn, it will become manifest that the protein frame directs and modulates the properties and reactivity of the metal ions.

Molecular dynamics study on conformational differences between dGMP and 8-oxo-dGMP: Effects of metal ions

The modified nucleotide base 7,8-dihydro-8-oxo-guanine (8-oxo-G) is one of the major sources of spontaneous mutagenesis. Nucleotide-sanitizing enzymes, such as the MutT homolog-1 (MTH1) and nudix-type motif 5 (NUDT5), selectively remove 8-oxo-G from the cellular pool of nucleotides. Previous studies showed that, although the syn conformation generally predominates in purine nucleotides with a bulky substituent at the 8-position, 8-oxo-dGMP binds to both MTH1 and NUDT5 in the anti conformation. This study was initiated to investigate the possibility that 8-oxo-dGMP itself may adopt the anti conformation. Molecular dynamics simulations of mononucleotides (dGMP, 8-oxo-dGMP) in aqueous solution were performed. 8-oxo-dGMP adopted the anti conformation as well as the syn conformation, and the proportion of adopting the anti conformation increased in the presence of metal ions. When 8-oxo-dGMP was in the anti conformation, a metal ion was located between the oxygen atom of phosphate and the oxygen atom at the 8-position of 8-oxo-G. The types of stable anti conformations of 8-oxo-dGMP differed, depending on the ionic radii and charges of coexisting ions. These data suggested a role for metal ions, other than as cofactors for the hydrolysis of the di- and tri-phosphate forms of mononucleotides; that the metal ions help retain the anti conformation of the N-glycosidic torsion angle of 8-oxo-dGMP to promote the binding between the 8-oxo-G deoxynucleotide and the nucleotide-sanitizing enzymes.

Force field independent metal parameters using a nonbonded dummy model

The cationic dummy atom approach provides a powerful nonbonded description for a range of alkaline-earth and transition-metal centers, capturing both structural and electrostatic effects. In this work we refine existing literature parameters for octahedrally coordinated Mn²⁺, Zn²⁺, Mg²⁺, and Ca²⁺, as well as providing new parameters for Ni²⁺, Co²⁺, and Fe²⁺. In all the cases, we are able to reproduce both M²⁺–O distances and experimental solvation free energies, which has not been achieved to date for transition metals using any other model. The parameters have also been tested using two different water models and show consistent performance. Therefore, our parameters are easily transferable to any force field that describes nonbonded interactions using Coulomb and Lennard-Jones potentials. Finally, we demonstrate the stability of our parameters in both the human and Escherichia coli variants of the enzyme glyoxalase I as showcase systems, as both enzymes are active with a range of transition metals. The parameters presented in this work provide a valuable resource for the molecular simulation community, as they extend the range of metal ions that can be studied using classical approaches, while also providing a starting point for subsequent parametrization of new metal centers.

Modeling catalytic promiscuity in the alkaline phosphatase superfamily

In recent years, it has become increasingly clear that promiscuity plays a key role in the evolution of new enzyme function. This finding has helped to elucidate fundamental aspects of molecular evolution. While there has been extensive experimental work on enzyme promiscuity, computational modeling of the chemical details of such promiscuity has traditionally fallen behind the advances in experimental studies, not least due to the nearly prohibitive computational cost involved in examining multiple substrates with multiple potential mechanisms and binding modes in atomic detail with a reasonable degree of accuracy. However, recent advances in both computational methodologies and power have allowed us to reach a stage in the field where we can start to overcome this problem, and molecular simulations can now provide accurate and efficient descriptions of complex biological systems with substantially less computational cost. This has led to significant advances in our understanding of enzyme function and evolution in a broader sense. Here, we will discuss currently available computational approaches that can allow us to probe the underlying molecular basis for enzyme specificity and selectivity, discussing the inherent strengths and weaknesses of each approach. As a case study, we will discuss recent computational work on different members of the alkaline phosphatase superfamily (AP) using a range of different approaches, showing the complementary insights they have provided. We have selected this particular superfamily, as it poses a number of significant challenges for theory, ranging from the complexity of the actual reaction mechanisms involved to the reliable modeling of the catalytic metal centers, as well as the very large system sizes. We will demonstrate that, through current advances in methodologies, computational tools can provide significant insight into the molecular basis for catalytic promiscuity, and, therefore, in turn, the mechanisms of protein functional evolution.

Substrate recognition by norovirus polymerase: microsecond molecular dynamics study

Molecular dynamics simulations of complexes between Norwalk virus RNA dependent RNA polymerase and its natural CTP and 2dCTP (both containing the O5'-C5'-C4'-O4' sequence of atoms bridging the triphosphate and sugar moiety) or modified coCTP (C5'-O5'-C4'-O4'), cocCTP (C5'-O5'-C4'-C4'') substrates were produced by means of CUDA programmable graphical processing units and the ACEMD software package. It enabled us to gain microsecond MD trajectories clearly showing that similar nucleoside triphosphates can bind surprisingly differently into the active site of the Norwalk virus RNA dependent RNA polymerase. It corresponds to their different modes of action (CTP-substrate, 2dCTP-poor substrate, coCTP-chain terminator, cocCTP-inhibitor). Moreover, extremely rare events-as repetitive pervasion of Arg182 into a potentially reaction promoting arrangement-were captured.

Conformational fluctuations of UreG, an intrinsically disordered enzyme

UreG proteins are small GTP binding (G) proteins that catalyze the hydrolysis of GTP necessary for the maturation of urease, a virulence factor in bacterial pathogenesis. UreG proteins are the first documented cases of intrinsically disordered enzymes. The comprehension of the dynamics of folding-unfolding events occurring in this protein could shed light on the enzymatic mechanism of UreG. Here, we used the recently developed replica exchange with solute tempering (REST2) computational methodology to explore the conformational space of UreG from Helicobacter pylori (HpUreG) and to identify its structural fluctuations. The same simulation and analysis protocol has been applied to HypB from Methanocaldococcus jannaschii (MjHypB), which is closely related to UreG in both sequence and function, even though it is not intrinsically disordered. A comparison of the two systems reveals that both HpUreG and MjHypB feature a substantial rigidity of the protein regions involved in catalysis, justifying its residual catalytic activity. On the other hand, HpUreG tends to unfold more than MjHypB in portions involved in protein-protein interactions with metallochaperones necessary for the formation of multiprotein complexes known to be involved in urease activation.

How calcium inhibits the magnesium-dependent kinase gsk3β: a molecular simulation study

Glycogen synthase kinase 3β (GSK3β) is a ubiquitous serine/threonine kinase that plays a pivotal role in many biological processes. GSK3β catalyzes the transfer of γ-phosphate of ATP to the unique substrate Ser/Thr residues with the assistance of two natural activating cofactors Mg(2+). Interestingly, the biological observation reveals that a non-native Ca(2+) ion can inhibit the GSK3β catalytic activity. Here, the inhibitory mechanism of GSK3β by the displacement of native Mg(2+) at site 1 by Ca(2+) was investigated by means of 80 ns comparative molecular dynamics (MD) simulations of the GSK3β···Mg(2+)-2/ATP/Mg(2+) -1 and GSK3β···Mg(2+)-2/ATP/Ca(2+)-1 systems. MD simulation results revealed that using the AMBER point charge model force field for Mg(2+) was more appropriate in the reproduction of the active site architectural characteristics of GSK3β than using the magnesium-cationic dummy atom model force field. Compared with the native Mg(2+) bound system, the misalignment of the critical triphosphate moiety of ATP, the erroneous coordination environments around the Mg(2+) ion at site 2, and the rupture of the key hydrogen bond between the invariant Lys85 and the ATP O(β2) atom in the Ca(2+) substituted system were observed in the MD simulation due to the Ca(2+) ion in active site in order to achieve its preferred sevenfold coordination geometry, which adequately abolish the enzymatic activity. The obtained results are valuable in understanding the possible mechanism by why Ca(2+) inhibits the GSK3β activity and also provide insights into the mechanism of Ca(2+) inhibition in other structurally related protein kinases.

Why nature really chose phosphate

Phosphoryl transfer plays key roles in signaling, energy transduction, protein synthesis, and maintaining the integrity of the genetic material. On the surface, it would appear to be a simple nucleophile displacement reaction. However, this simplicity is deceptive, as, even in aqueous solution, the low-lying d-orbitals on the phosphorus atom allow for eight distinct mechanistic possibilities, before even introducing the complexities of the enzyme catalyzed reactions. To further complicate matters, while powerful, traditional experimental techniques such as the use of linear free-energy relationships (LFER) or measuring isotope effects cannot make unique distinctions between different potential mechanisms. A quarter of a century has passed since Westheimer wrote his seminal review, 'Why Nature Chose Phosphate' (Science 235 (1987), 1173), and a lot has changed in the field since then. The present review revisits this biologically crucial issue, exploring both relevant enzymatic systems as well as the corresponding chemistry in aqueous solution, and demonstrating that the only way key questions in this field are likely to be resolved is through careful theoretical studies (which of course should be able to reproduce all relevant experimental data). Finally, we demonstrate that the reason that nature really chose phosphate is due to interplay between two counteracting effects: on the one hand, phosphates are negatively charged and the resulting charge-charge repulsion with the attacking nucleophile contributes to the very high barrier for hydrolysis, making phosphate esters among the most inert compounds known. However, biology is not only about reducing the barrier to unfavorable chemical reactions. That is, the same charge-charge repulsion that makes phosphate ester hydrolysis so unfavorable also makes it possible to regulate, by exploiting the electrostatics. This means that phosphate ester hydrolysis can not only be turned on, but also be turned off, by fine tuning the electrostatic environment and the present review demonstrates numerous examples where this is the case. Without this capacity for regulation, it would be impossible to have for instance a signaling or metabolic cascade, where the action of each participant is determined by the fine-tuned activity of the previous piece in the production line. This makes phosphate esters the ideal compounds to facilitate life as we know it.