Dynamics of 5R-Tg Base Flipping in DNA Duplexes Based on Simulations─Agreement with Experiments and Beyond.
Ontology highlight
ABSTRACT:
SUBMITTER: Wang SD
PROVIDER: S-EPMC8790752 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:DNA base flipping is a fundamental theme in DNA biophysics. The dynamics for a B-DNA base to spontaneously flip out of the double helix has significant implications in various DNA-protein interactions but are still poorly understood. The spontaneous base-flipping rate obtained previously via the imino proton exchange assay is most likely the rate of base wobbling instead of flipping. Using the diffusion-decelerated fluorescence correlation spectroscopy together with molecular dynamics simulations, we show that a base of a single mismatched base pair (T-G, T-T, or T-C) in a double-stranded DNA can spontaneously flip out of the DNA duplex. The extrahelical lifetimes are on the order of 10 ms, whereas the intrahelical lifetimes range from 0.3 to 20 s depending on the stability of the base pairs. These findings provide detailed understanding on the dynamics of DNA base flipping and lay down foundation to fully understand how exactly the repair proteins search and locate the target mismatched base among a vast excess of matched DNA bases.
| S-EPMC4050552 | biostudies-literature
Project description:DNA base modifications and mutations are observed in all genomes throughout the kingdoms of life. Proteins involved in their establishment and removal were shown to use a base flipping mechanism to access their substrates. To better understand how proteins flip DNA bases to modify or remove them, we optimized and developed a pipeline of methods to step-by-step detect the process starting with protein-DNA interaction, base flipping itself and the ensuing DNA base modification or excision. As methylcytosine is the best-studied DNA modification, here we focus on the process of writing, modifying and reading this DNA base. Using multicolor electrophoretic mobility shift assays, we show that the methylcytosine modifier Tet1 exhibits little DNA sequence specificity with only a slight preference for methylated CpG containing DNA. A combination of chloroacetaldehyde treatment and high-resolution melting temperature analysis allowed us to detect base flipping induced by the methylcytosine modifier Tet1 as well as the methylcytosine writer M.HpaII. Finally, we show that high-resolution melting temperature analysis can be used to detect the activity of glycosylases, methyltransferases and dioxigenases on DNA substrates. Taken together, this DNA base flipping analytical pipeline (BaFAP) provide a complete toolbox for the fast and sensitive analysis of proteins that bind, flip and modify or excise DNA bases.
| S-EPMC6994035 | biostudies-literature
Project description:Individual DNA bases are known to be able to flip out of the helical stack, providing enzymes with access to the genetic information otherwise hidden inside the helix. Consequently, base flipping is a necessary first step to many more complex biological processes such as DNA transcription or replication. Much remains unknown about this elementary step, despite a wealth of experimental and theoretical studies. From the theoretical point of view, the involved timescale of milliseconds or longer requires the use of enhanced sampling techniques. In contrast to previous theoretical studies employing umbrella sampling along a predefined flipping coordinate, this study attempts to induce flipping without prior knowledge of the pathway, using information from a molecular dynamics simulation of a B-DNA fragment and the conformational flooding method. The relevance to base flipping of the principal components of the simulation is assayed, and a combination of modes optimally related to the flipping of the base through either helical groove is derived for each of the two bases of the central guanine-cytosine basepair. By applying an artificial flooding potential along these collective coordinates, the flipping mechanism is accelerated to within the scope of molecular dynamics simulations. The associated free energy surface is found to feature local minima corresponding to partially flipped states, particularly relevant to flipping in isolated DNA; further transitions from these minima to the fully flipped conformation are accelerated by additional flooding potentials. The associated free energy profiles feature similar barrier heights for both bases and pathways; the flipped state beyond is a broad and rugged attraction basin, only a few kcal/mol higher in energy than the closed conformation. This result diverges from previous works but echoes some aspects of recent experimental findings, justifying the need for novel approaches to this difficult problem: this contribution represents a first step in this direction. Important structural factors involved in flipping, both local (sugar-phosphate backbone dihedral angles) and global (helical axis bend), are also identified.
| S-EPMC1913169 | biostudies-literature
Project description:On the basis of non-covalent bond interactions in nucleic acids, we synthesized the deoxyadenosine derivatives tethering a phenyl group (X) and a naphthyl group (Z) by an amide linker, which mimic a Watson-Crick base pair. Circular dichroism spectra indicated that the duplexes containing X and Z formed a similar conformation regardless of the opposite nucleotide species (A, G, C, T and an abasic site analogue F), which was not observed for the natural duplexes. The values among the natural duplexes containing the A/A, A/G, A/C, A/T and A/F pairs differed by 5.2 kcal mol(-1) while that among the duplexes containing X or Z in place of the adenine differed by only 1.9 or 2.8 kcal mol(-1), respectively. Fluorescence quenching experiments confirmed that 2-amino purine opposite X adopted an unstacked conformation. The structural and thermodynamic analyses suggest that the aromatic hydrocarbon group of X and Z intercalates into a double helix, resulting in the opposite nucleotide base flipping into an unstacked position regardless of the nucleotide species. This observation implies that modifications at the aromatic hydrocarbon group and the amide linker may expand the application of the base pair-mimic nucleosides for molecular biology and biotechnology.
| S-EPMC1316115 | biostudies-literature
Project description:The crystal structure of Rad4/Rad23, the yeast homolog of the human nucleotide excision repair (NER) lesion recognition factor XPC-RAD23B ( Min , J. H. and Pavletich , N. P. ( 2007 ) Nature 449 , 570 - 575 ) reveals that the lesion-partner base is flipped out of the helix and binds to amino acids of the protein. This suggests the hypothesis that the flipping of this partner base must overcome a free energy barrier, which constitutes one element contributing to changes in the thermodynamic properties induced by the DNA damage and sensed by the recognition protein. We explored this hypothesis by computing complete flipping free energy profiles for two lesions derived from the procarcinogenic polycyclic aromatic hydrocarbons (PAHs), dibenzo[a,l]pyrene (DB[a,l]P) and benzo[a]pyrene (B[a]P), R-trans-anti-DB[a,l]P-N(6)-dA (R-DB[a,l]P-dA) and R-trans-anti-B[a]P-N(6)-dA (R-B[a]P-dA), and the corresponding unmodified duplex. The DB[a,l]P and B[a]P adducts differ in number and organization of their aromatic rings. We integrate these results with prior profiles for the R-trans-anti-DB[a,l]P-dG adduct ( Zheng , H. et al. ( 2010 ) Chem. Res. Toxicol. 23 , 1868 - 1870 ). All adopt conformational themes involving intercalation of the PAH aromatic ring system into the DNA duplex; however, R-DB[a,l]P-dA and R-B[a]P-dA intercalate from the major groove, while R-DB[a,l]P-dG intercalates from the minor groove. These structural differences produce different computed van der Waals stacking interaction energies between the flipping partner base with the lesion aromatic ring system and adjacent bases; we find that the better the stacking, the higher the relative flipping free energy barrier and hence lower flipping probability. The better relative NER susceptibilities correlate with greater ease of flipping in these three differently intercalated lesions. In addition to partner base flipping, the Rad4/Rad23 crystal structure shows that a protein-β-hairpin, BHD3, intrudes from the major groove side between the DNA strands at the lesion site. We present a molecular modeling study for the R-DB[a,l]P-dG lesion in Rad4/Rad23 showing BHD3 β-hairpin intrusion with lesion eviction, and we hypothesize that lesion steric effects play a role in the recognition of intercalated adducts.
| S-EPMC3865170 | biostudies-literature
Project description:Locked nucleic acid (LNA) is a chemically modified nucleic acid with its sugar ring locked in an RNA-like (C3'-endo) conformation. LNAs show extraordinary thermal stabilities when hybridized with DNA, RNA or LNA itself. We performed molecular dynamics simulations on five isosequential duplexes (LNA-DNA, LNA-LNA, LNA-RNA, RNA-DNA and RNA-RNA) in order to characterize their structure, dynamics and hydration. Structurally, the LNA-DNA and LNA-RNA duplexes are found to be similar to regular RNA-DNA and RNA-RNA duplexes, whereas the LNA-LNA duplex is found to have its helix partly unwound and does not resemble RNA-RNA duplex in a number of properties. Duplexes with an LNA strand have on average longer interstrand phosphate distances compared to RNA-DNA and RNA-RNA duplexes. Furthermore, intrastrand phosphate distances in LNA strands are found to be shorter than in DNA and slightly shorter than in RNA. In case of induced sugar puckering, LNA is found to tune the sugar puckers in partner DNA strand toward C3'-endo conformations more efficiently than RNA. The LNA-LNA duplex has lesser backbone flexibility compared to the RNA-RNA duplex. Finally, LNA is less hydrated compared to DNA or RNA but is found to have a well-organized water structure.
| S-EPMC2275159 | biostudies-literature
Project description:Mixed double helices formed by RNA and DNA strands, commonly referred to as hybrid duplexes or hybrids, are essential in biological processes like transcription and reverse transcription. They are also important for their applications in CRISPR gene editing and nanotechnology. Yet, despite their significance, the hybrid duplexes have been seldom modeled by atomistic molecular dynamics methodology, and there is no benchmark study systematically assessing the force-field performance. Here, we present an extensive benchmark study of polypurine tract (PPT) and Dickerson-Drew dodecamer hybrid duplexes using contemporary and commonly utilized pairwise additive and polarizable nucleic acid force fields. Our findings indicate that none of the available force-field choices accurately reproduces all the characteristic structural details of the hybrid duplexes. The AMBER force fields are unable to populate the C3'-endo (north) pucker of the DNA strand and underestimate inclination. The CHARMM force field accurately describes the C3'-endo pucker and inclination but shows base pair instability. The polarizable force fields struggle with accurately reproducing the helical parameters. Some force-field combinations even demonstrate a discernible conflict between the RNA and DNA parameters. In this work, we offer a candid assessment of the force-field performance for mixed DNA/RNA duplexes. We provide guidance on selecting utilizable force-field combinations and also highlight potential pitfalls and best practices for obtaining optimal performance.
| S-EPMC11325551 | biostudies-literature
Project description:Gramicidin A (gA) channels provide an ideal system to test molecular dynamics (MD) simulations of membrane proteins. The peptide backbone lines a cation-selective pore, and due to the small channel size, the average structure and extent of fluctuations of all atoms in the peptide will influence ion permeation. This raises the question of how well molecular mechanical force fields used in MD simulations and potential of mean force (PMF) calculations can predict structure and dynamics as well as ion permeation. To address this question, we undertook a comparative study of nuclear magnetic resonance (NMR) observables predicted by fully atomistic MD simulations on a gA dimer embedded in a sodium dodecyl sulfate (SDS) micelle with measurements of the gA dimer backbone and tryptophan side chain dynamics using solution-state (15)N NMR on gA dimers in SDS micelles (Vostrikov, V. V.; Gu, H.; Ingólfsson, H. I.; Hinton, J. F.; Andersen, O. S.; Roux, B.; Koeppe, R. E., II. J. Phys. Chem. B2011, DOI 10.1021/jp200906y , accompanying article). This comparison enables us to examine the robustness of the MD simulations done using different force fields as well as their ability to predict important features of the gA channel. We find that MD is able to predict NMR observables, including the generalized order parameters (S(2)), the (15)N spin-lattice (T(1)) and spin-spin (T(2)) relaxation times, and the (1)H-(15)N nuclear Overhauser effect (NOE), with remarkable accuracy. To examine further how differences in the force fields can affect the channel conductance, we calculated the PMF for K(+) and Na(+) permeation through a gA channel in a dimyristoylphosphatidylcholine (DMPC) bilayer. In this case, we find that MD is less successful in quantitatively predicting the single-channel conductance.
| S-EPMC3107394 | biostudies-literature
Project description:DNA mismatch recognition and repair is vital for preserving the fidelity of the genome. Conserved across prokaryotes and eukaryotes, MutS is the primary protein that is responsible for recognizing a variety of DNA mismatches. From molecular dynamics simulations of the Escherichia coli MutS-DNA complex, we describe significant conformational dynamics in the DNA surrounding a G·T mismatch that involves weakening of the basepair hydrogen bonding in the basepair adjacent to the mismatch and, in one simulation, complete base opening via the major groove. The energetics of base flipping was further examined with Hamiltonian replica exchange free energy calculations revealing a stable flipped-out state with an initial barrier of ~2 kcal/mol. Furthermore, we observe changes in the local DNA structure as well as in the MutS structure that appear to be correlated with base flipping. Our results suggest a role of base flipping as part of the repair initiation mechanism most likely leading to sliding-clamp formation.
| S-EPMC3207177 | biostudies-literature
Project description:Twenty AdoMet-dependent methyltransferases (MTases) have been characterized structurally by X-ray crystallography and NMR. These include seven DNA MTases, five RNA MTases, four protein MTases and four small molecule MTases acting on the carbon, oxygen or nitrogen atoms of their substrates. The MTases share a common core structure of a mixed seven-stranded beta-sheet (6 downward arrow 7 upward arrow 5 downward arrow 4 downward arrow 1 downward arrow 2 downward arrow 3 downward arrow) referred to as an 'AdoMet-dependent MTase fold', with the exception of a protein arginine MTase which contains a compact consensus fold lacking the antiparallel hairpin strands (6 downward arrow 7 upward arrow). The consensus fold is useful to identify hypothetical MTases during structural proteomics efforts on unannotated proteins. The same core structure works for very different classes of MTase including those that act on substrates differing in size from small molecules (catechol or glycine) to macromolecules (DNA, RNA and protein). DNA MTases use a 'base flipping' mechanism to deliver a specific base within a DNA molecule into a typically concave catalytic pocket. Base flipping involves rotation of backbone bonds in double-stranded DNA to expose an out-of-stack nucleotide, which can then be a substrate for an enzyme-catalyzed chemical reaction. The phenomenon is fully established for DNA MTases and for DNA base excision repair enzymes, and is likely to prove general for enzymes that require access to unpaired, mismatched or damaged nucleotides within base-paired regions in DNA and RNA. Several newly discovered MTase families in eukaryotes (DNA 5mC MTases and protein arginine and lysine MTases) offer new challenges in the MTase field.
| S-EPMC55914 | biostudies-literature