Insertion of oxidized nucleotide triggers rapid DNA polymerase opening

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.

A Role for N6-Methyladenine in DNA Damage Repair

The leading cause of mutation due to oxidative damage is 8-oxo-2'-deoxyguanosine (8-oxoG) mispairing with adenine (Ade), which can occur in two ways. First, guanine of a G:C DNA base pair can be oxidized. If not repaired in time, DNA polymerases can mispair Ade with 8-oxoG in the template. This 8-oxoG:A can be repaired by enzymes that remove Ade opposite to template 8-oxoG, or 8-oxoG opposite to Cyt. Second, free 8-oxo-dGTP can be misincorporated by DNA polymerases into DNA opposite template Ade. However, there is no known repair activity that removes 8-oxoG opposite to template Ade. We suggest that a major role of N6-methyladenine in mammalian DNA is minimizing incorporation of 8-oxoG opposite to Ade by DNA polymerases following adduct formation.

Competitive Binding of Mg2+ and Na+ Ions to Nucleic Acids: From Helices to Tertiary Structures

Nucleic acids generally reside in cellular aqueous solutions with mixed divalent/monovalent ions, and the competitive binding of divalent and monovalent ions is critical to the structures of nucleic acids because of their polyanionic nature. In this work, we first proposed a general and effective method for simulating a nucleic acid in mixed divalent/monovalent ion solutions with desired bulk ion concentrations via molecular dynamics (MD) simulations and investigated the competitive binding of Mg²⁺/Na⁺ ions to various nucleic acids by all-atom MD simulations. The extensive MD-based examinations show that single MD simulations conducted using the proposed method can yield desired bulk divalent/monovalent ion concentrations for various nucleic acids, including RNA tertiary structures. Our comprehensive analyses show that the global binding of Mg²⁺/Na⁺ to a nucleic acid is mainly dependent on its structure compactness, as well as Mg²⁺/Na⁺ concentrations, rather than the specific structure of the nucleic acid. Specifically, the relative global binding of Mg²⁺ over Na⁺ is stronger for a nucleic acid with higher effective surface charge density and higher relative Mg²⁺/Na⁺ concentrations. Furthermore, the local binding of Mg²⁺/Na⁺ to a phosphate of a nucleic acid mainly depends on the local phosphate density in addition to Mg²⁺/Na⁺ concentrations.

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.

Computational Simulations of DNA Polymerases: Detailed Insights on Structure/Function/Mechanism from Native Proteins to Cancer Variants

Genetic information is vital in the cell cycle of DNA-based organisms. DNA polymerases (DNA Pols) are crucial players in transactions dealing with these processes. Therefore, the detailed understanding of the structure, function, and mechanism of these proteins has been the focus of significant effort. Computational simulations have been applied to investigate various facets of DNA polymerase structure and function. These simulations have provided significant insights over the years. This perspective presents the results of various computational studies that have been employed to research different aspects of DNA polymerases including detailed reaction mechanism investigation, mutagenicity of different metal cations, possible factors for fidelity synthesis, and discovery/functional characterization of cancer-related mutations on DNA polymerases.

Divalent Ion-Mediated DNA-DNA Interactions: A Comparative Study of Triplex and Duplex

Ion-mediated interaction between DNAs is essential for DNA condensation, and it is generally believed that monovalent and nonspecifically binding divalent cations cannot induce the aggregation of double-stranded (ds) DNAs. Interestingly, recent experiments found that alkaline earth metal ions such as Mg²⁺ can induce the aggregation of triple-stranded (ts) DNAs, although there is still a lack of deep understanding of the surprising findings at the microscopic level. In this work, we employed all-atom dynamic simulations to directly calculate the potentials of mean force (PMFs) between tsDNAs, between dsDNAs, and between tsDNA and dsDNA in Mg²⁺ solutions. Our calculations show that the PMF between tsDNAs is apparently attractive and becomes more strongly attractive at higher [Mg²⁺], although the PMF between dsDNAs cannot become apparently attractive even at high [Mg²⁺]. Our analyses show that Mg²⁺ internally binds into grooves and externally binds to phosphate groups for both tsDNA and dsDNA, whereas the external binding of Mg²⁺ is much stronger for tsDNA. Such stronger external binding of Mg²⁺ for tsDNA favors more apparent ion-bridging between helices than for dsDNA. Furthermore, our analyses illustrate that bridging ions, as a special part of external binding ions, are tightly and positively coupled to ion-mediated attraction between DNAs.

The m6A methylation perturbs the Hoogsteen pairing-guided incorporation of an oxidized nucleotide

This study describes the structural implications and properties of m⁶A in reducing the incorporation of an oxidized nucleotide into DNA.

Natural nucleic acid bases can form Watson–Crick (WC) or Hoogsteen (HG) base pairs. Importantly, 8-oxo-2′-deoxyguanosine (8-oxo-dG) in DNA or 8-oxo-dG 5′-triphosphate (8-oxo-dGTP) favors a syn conformation because of the steric repulsion between O8 and O4′ of the deoxyribose ring. 8-oxo-dGTP can be incorporated into DNA opposite the templating adenine (A) using HG pairing as the dominant mechanism. Both RNA and DNA can be methylated at the N6 position of A to form N⁶-methyladenine (m⁶A). It has been found that certain viral infections may trigger an increase in the production of both 8-oxo-dGTP and m⁶A. The current study aims to systematically explore the effects of m⁶A methylation on HG base pairs and the consequent nucleotide incorporation. Our thermodynamic melting study shows that the m⁶A·8-oxo-dG is significantly less stable than the A·8-oxo-dG base pair in the paired region of a DNA duplex. Moreover, we have used pre-steady-state kinetics to examine the incorporation of 8-oxo-dGTP opposite m⁶A relative to A by a variety of reverse transcriptase (RT) enzymes and DNA polymerase (DNA pol) enzymes such as the human immunodeficiency virus type 1 (HIV-1) RT and human DNA pol β. The results demonstrate that all of these enzymes incorporate 8-oxo-dGTP less efficiently opposite m⁶A relative to A. Considering the steric bulk of the purine–purine pair between 8-oxo-dG and A, m⁶A methylation may affect the HG pairing to a great extent. Hence, it will be unfavorable to incorporate 8-oxo-dGTP into the growing strand opposite m⁶A. Moreover, the impeded incorporation of 8-oxo-dGTP opposite m⁶A has been extended to determine m⁶A at pre-defined positions in human rRNA. Our study may provide new insights into the roles of m⁶A in reducing the mutagenic potential of cellular 8-oxo-dGTP.

Protocols for Molecular Dynamics Simulations of RNA Nanostructures

Molecular dynamics (MD) simulations have been used as one of the main research tools to study a wide range of biological systems and bridge the gap between X-ray crystallography or NMR structures and biological mechanism. In the field of RNA nanostructures, MD simulations have been used to fix steric clashes in computationally designed RNA nanostructures, characterize the dynamics, and investigate the interaction between RNA and other biomolecules such as delivery agents and membranes.In this chapter we present examples of computational protocols for molecular dynamics simulations in explicit and implicit solvent using the Amber Molecular Dynamics Package. We also show examples of post-simulation analysis steps and briefly mention selected tools beyond the Amber package. Limitations of the methods, tools, and protocols are also discussed. Most of the examples are illustrated for a small RNA duplex (helix), but the protocols are applicable to any nucleic acid structure, subject only to the computational speed and memory limitations of the hardware available to the user.