1
|
Kim H, Pak Y. Isomerization Pathways of a Mismatched Base Pair of A:8OG in Free Duplex DNA. J Chem Inf Model 2024; 64:4511-4517. [PMID: 38767002 DOI: 10.1021/acs.jcim.4c00563] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
The A:8OG base pair (bp) is the outcome of DNA replication of the mismatched C:8OG bp. A high A:8OG bp population increases the C/G to A/T transversion mutation, which is responsible for various diseases. MutY is an important enzyme in the error-proof cycle and reverts A:8OG to C:8OG bp by cleaving adenine from the A:8OG bp. Several X-ray crystallography studies have determined the structure of MutY during the lesion scanning and lesion recognition stages. Interestingly, glycosidic bond (χ) angles of A:8OG bp in those two lesion recognition structures were found to differ, which implies that χ-torsion isomerization should occur during the lesion recognition process. In this study, as a first step to understanding this isomerization process, we characterized the intrinsic dynamic features of A:8OG in free DNAs by a free energy landscape simulation at the all-atom level. In this study, four isomerization states were assigned in the order of abundance: Aanti:8OGsyn > Aanti:8OGanti > Asyn:8OGanti ≈ Asyn:8OGsyn. Of these bp states, only 8OG in Asyn:8OGanti was located in the extrahelical space, whereas the purine bases (A and 8OG) in the other bp states remained inside the DNA helix. Also, free energy landscapes showed that the isomerization processes connecting these four bp states proceeded mostly in the intrahelical space via successive single glycosidic bond rotations of either A or 8OG.
Collapse
Affiliation(s)
- Hyeonjun Kim
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 46241, South Korea
| | - Youngshang Pak
- Department of Chemistry and Institute of Functional Materials, Pusan National University, Busan 46241, South Korea
| |
Collapse
|
2
|
Davila-Calderon J, Li ML, Penumutchu SR, Haddad C, Malcolm L, King J, Hargrove AE, Brewer G, Tolbert BS. Enterovirus evolution reveals the mechanism of an RNA-targeted antiviral and determinants of viral replication. SCIENCE ADVANCES 2024; 10:eadg3060. [PMID: 38363831 PMCID: PMC10871541 DOI: 10.1126/sciadv.adg3060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/14/2022] [Accepted: 01/16/2024] [Indexed: 02/18/2024]
Abstract
Selective pressures on viruses provide opportunities to establish target site specificity and mechanisms of antivirals. Enterovirus (EV)-A71 with resistant mutations in the stem loop (SL) II internal ribosome entry site (IRES) (SLIIresist) were selected at low doses of the antiviral dimethylamiloride (DMA)-135. The EV-A71 mutants were resistant to DMA-135 at concentrations that inhibit replication of wild-type virus. EV-A71 IRES structures harboring resistant mutations induced efficient expression of Luciferase messenger RNA in the presence of noncytotoxic doses of DMA-135. Nuclear magnetic resonance indicates that the mutations change the structure of SLII at the binding site of DMA-135 and at the surface recognized by the host protein AU-rich element/poly(U)-binding/degradation factor 1 (AUF1). Biophysical studies of complexes formed between AUF1, DMA-135, and either SLII or SLIIresist show that DMA-135 stabilizes a ternary complex with AUF1-SLII but not AUF1-SLIIresist. This work demonstrates how viral evolution elucidates the (DMA-135)-RNA binding site specificity in cells and provides insights into the viral pathways inhibited by the antiviral.
Collapse
Affiliation(s)
| | - Mei-Ling Li
- Department of Biochemistry and Molecular Biology, Rutgers Robert Wood Johnson Medical School, Piscataway, NJ, USA
| | | | - Christina Haddad
- Department of Chemistry, Case Western Reserve University, Cleveland, OH, USA
| | - Linzy Malcolm
- Department of Chemistry, Case Western Reserve University, Cleveland, OH, USA
| | - Josephine King
- Department of Chemistry, Case Western Reserve University, Cleveland, OH, USA
| | | | - Gary Brewer
- Department of Biochemistry and Molecular Biology, Rutgers Robert Wood Johnson Medical School, Piscataway, NJ, USA
| | - Blanton S. Tolbert
- Department of Chemistry, Case Western Reserve University, Cleveland, OH, USA
- Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA, USA
- Howard Hughes Medical Institute, Chevy Chase, MD, USA
| |
Collapse
|
3
|
Stone S, Ray D, Andricioaei I. Force-Field-Dependent DNA Breathing Dynamics: A Case Study of Hoogsteen Base Pairing in A6-DNA. J Chem Inf Model 2022; 62:6749-6761. [PMID: 36049242 PMCID: PMC9795553 DOI: 10.1021/acs.jcim.2c00519] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The Hoogsteen (HG) base pairing conformation, commonly observed in damaged and mutated DNA helices, facilitates DNA repair and DNA recognition. The free energy difference between HG and Watson-Crick (WC) base pairs has been computed in previous studies. However, the mechanism of the conformational transition is not well understood. A detailed understanding of the process of WC to HG base pair transition can provide a deeper understanding of DNA repair and recognition. In an earlier study, we explored the free energy landscape for this process using extensive computer simulation with the CHARMM36 force field. In this work, we study the impact of force field models in describing the WC to HG base pairing transition using meta-eABF enhanced sampling, quasi-harmonic entropy calculation, and nonbonded energy analysis. The secondary structures of both base pairing forms and the topology of the free energy landscapes were consistent over different force field models, although the relative free energy, entropy, and the interaction energies tend to vary. The relative stability of the WC and HG conformations is dictated by a delicate balance between the enthalpic stabilization and the reduced entropy of the structurally rigid HG structure. These findings highlight the impact that subtleties in force field models can have on accurately modeling DNA base pair dynamics and should stimulate further computational investigations into other dynamically important motions in DNA.
Collapse
Affiliation(s)
- Sharon
Emily Stone
- Department
of Chemistry, University of California Irvine, Irvine, California 92697, United States
| | - Dhiman Ray
- Department
of Chemistry, University of California Irvine, Irvine, California 92697, United States
| | - Ioan Andricioaei
- Department
of Chemistry, University of California Irvine, Irvine, California 92697, United States,Department
of Physics and Astronomy, University of
California Irvine, Irvine, California 92697, United States,
| |
Collapse
|
4
|
Geronimo I, De Vivo M. Alchemical Free-Energy Calculations of Watson-Crick and Hoogsteen Base Pairing Interconversion in DNA. J Chem Theory Comput 2022; 18:6966-6973. [PMID: 36201305 DOI: 10.1021/acs.jctc.2c00848] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Hoogsteen (HG) base pairs have a transient nature and can be structurally similar to Watson-Crick (WC) base pairs, making their occurrence and thermodynamic stability difficult to determine experimentally. Herein, we employed the restrain-free-energy perturbation-release (R-FEP-R) method to calculate the relative free energy of the WC and HG base pairing modes in isolated and bound DNA systems and predict the glycosyl torsion conformational preference of purine bases. Notably, this method does not require prior knowledge of the transition pathway between the two end states. Remarkably, relatively fast convergence was reached, with results in excellent agreement with experimental data for all the examined DNA systems. The R-REP-R method successfully determined the stability of HG base pairing and more generally, the conformational preference of purine bases, in these systems. Therefore, this computational approach can help to understand the dynamic equilibrium between the WC and HG base pairing modes in DNA.
Collapse
Affiliation(s)
- Inacrist Geronimo
- Laboratory of Molecular Modelling & Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, Genoa 16163, Italy
| | - Marco De Vivo
- Laboratory of Molecular Modelling & Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, Genoa 16163, Italy
| |
Collapse
|
5
|
Galindo-Murillo R, Cheatham TE. Transient Hoogsteen Base Pairs Observed in Unbiased Molecular Dynamics Simulations of DNA. J Phys Chem Lett 2022; 13:6283-6287. [PMID: 35775742 PMCID: PMC9466961 DOI: 10.1021/acs.jpclett.2c01348] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Duplex DNA is modeled as canonical B-DNA displaying the characteristic Watson-Crick base pairs. A less common and short-lived pairing of the nucleobases is the Hoogsteen (HG) conformation. The low population of the HG base pairs (<1%) necessitates extended sampling times in order to analyze through unbiased molecular dynamics (MD) simulations. We have discovered that with extended sampling times using multiple independent copies of an 18-mer sequence, the MD trajectories reproduce the expected and transient HG base pairing. Consistent with experiment, the percentage of the HG events are within the range of ∼0.1-1.0% over the combined aggregate sampling time of more than 3.6 ms. We present the reliability of the current AMBER set of nucleic acid force fields and tools to accurately simulate naturally occurring base-pairing and opening events without any bias or restraints. The mechanism consists of base pair fraying, flipping of the purine, and reformation with HG base pairs.
Collapse
|
6
|
Sequence dependence of transient Hoogsteen base pairing in DNA. PLoS Comput Biol 2022; 18:e1010113. [PMID: 35617357 PMCID: PMC9177043 DOI: 10.1371/journal.pcbi.1010113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2021] [Revised: 06/08/2022] [Accepted: 04/19/2022] [Indexed: 11/19/2022] Open
Abstract
Hoogsteen (HG) base pairing is characterized by a 180° rotation of the purine base with respect to the Watson-Crick-Franklin (WCF) motif. Recently, it has been found that both conformations coexist in a dynamical equilibrium and that several biological functions require HG pairs. This relevance has motivated experimental and computational investigations of the base-pairing transition. However, a systematic simulation of sequence variations has remained out of reach. Here, we employ advanced path-based methods to perform unprecedented free-energy calculations. Our methodology enables us to study the different mechanisms of purine rotation, either remaining inside or after flipping outside of the double helix. We study seven different sequences, which are neighbor variations of a well-studied A⋅T pair in A6-DNA. We observe the known effect of A⋅T steps favoring HG stability, and find evidence of triple-hydrogen-bonded neighbors hindering the inside transition. More importantly, we identify a dominant factor: the direction of the A rotation, with the 6-ring pointing either towards the longer or shorter segment of the chain, respectively relating to a lower or higher barrier. This highlights the role of DNA’s relative flexibility as a modulator of the WCF/HG dynamic equilibrium. Additionally, we provide a robust methodology for future HG proclivity studies. Recently, an alternative DNA base-pairing conformation, known as Hoogsteen (HG), has been found to coexist with the well-known Watson-Crick-Franklin (WCF) pairing. Several experimental and computational studies have focused on this heterogeneity, as it is involved in various recognition and replication processes. The WCF-to-HG transition mechanisms consist of a ±180° rotation of the purine base, occurring either inside of the double helix or while flipping temporarily outside of it. Even though molecular dynamics simulations can provide fine details about the transition pathways and their free-energy barriers, the computational cost has limited most studies to focus on only one particular chain (A6-DNA). Here, we investigate the sequence dependence of the base-pairing transition, by systematically varying the direct neighbors of a transitioning A⋅T pair; probing inside and outside pathways in seven distinct systems. We discover that triple-hydrogen-bonded neighboring base-pairs hinder the inside rotation mechanism, due to the reduced flexibility needed for internal base rotation. Across all sequences, we confirm that outside transitions have a lower free-energy barrier. Most importantly, we observe that the direction of the A rotation, with the A 6-ring pointing either towards the long or short end of the modelled DNA chain, has a determinant effect on the height of the free-energy barrier. These results point to a critical role of DNA’s small- and medium-scale flexibility in modulating the proclivity of HG base pairs; providing a handle that might be employed by several biological mechanisms.
Collapse
|
7
|
Ray D, Stone SE, Andricioaei I. Markovian Weighted Ensemble Milestoning (M-WEM): Long-Time Kinetics from Short Trajectories. J Chem Theory Comput 2021; 18:79-95. [PMID: 34910499 DOI: 10.1021/acs.jctc.1c00803] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We introduce a rare-event sampling scheme, named Markovian Weighted Ensemble Milestoning (M-WEM), which inlays a weighted ensemble framework within a Markovian milestoning theory to efficiently calculate thermodynamic and kinetic properties of long-time-scale biomolecular processes from short atomistic molecular dynamics simulations. M-WEM is tested on the Müller-Brown potential model, the conformational switching in alanine dipeptide, and the millisecond time-scale protein-ligand unbinding in a trypsin-benzamidine complex. Not only can M-WEM predict the kinetics of these processes with quantitative accuracy but it also allows for a scheme to reconstruct a multidimensional free-energy landscape along additional degrees of freedom, which are not part of the milestoning progress coordinate. For the ligand-receptor system, the experimental residence time, association and dissociation kinetics, and binding free energy could be reproduced using M-WEM within a simulation time of a few hundreds of nanoseconds, which is a fraction of the computational cost of other currently available methods, and close to 4 orders of magnitude less than the experimental residence time. Due to the high accuracy and low computational cost, the M-WEM approach can find potential applications in kinetics and free-energy-based computational drug design.
Collapse
Affiliation(s)
- Dhiman Ray
- Department of Chemistry, University of California Irvine, Irvine, California 92697, United States
| | - Sharon Emily Stone
- Department of Chemistry, University of California Irvine, Irvine, California 92697, United States
| | - Ioan Andricioaei
- Department of Chemistry, University of California Irvine, Irvine, California 92697, United States.,Department of Physics and Astronomy, University of California Irvine, Irvine, California 92697, United States
| |
Collapse
|