1
|
Ogbonna E, Paul A, Farahat AA, Terrell JR, Mineva E, Ogbonna V, Boykin DW, Wilson WD. X-ray Structure Characterization of the Selective Recognition of AT Base Pair Sequences. ACS BIO & MED CHEM AU 2023; 3:335-348. [PMID: 37599788 PMCID: PMC10436263 DOI: 10.1021/acsbiomedchemau.3c00002] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/11/2023] [Revised: 03/23/2023] [Accepted: 03/24/2023] [Indexed: 08/22/2023]
Abstract
The rational design of small molecules that target specific DNA sequences is a promising strategy to modulate gene expression. This report focuses on a diamidinobenzimidazole compound, whose selective binding to the minor groove of AT DNA sequences holds broad significance in the molecular recognition of AT-rich human promoter sequences. The objective of this study is to provide a more detailed and systematized understanding, at an atomic level, of the molecular recognition mechanism of different AT-specific sequences by a rationally designed minor groove binder. The specialized method of X-ray crystallography was utilized to investigate how the sequence-dependent recognition properties in general, A-tract, and alternating AT sequences affect the binding of diamidinobenzimidazole in the DNA minor groove. While general and A-tract AT sequences give a narrower minor groove, the alternating AT sequences intrinsically have a wider minor groove which typically constricts upon binding. A strong and direct hydrogen bond between the N-H of the benzimidazole and an H-bond acceptor atom in the minor groove is essential for DNA recognition in all sequences described. In addition, the diamidine compound specifically utilizes an interfacial water molecule for its DNA binding. DNA complexes of AATT and AAAAAA recognition sites show that the diamidine compound can bind in two possible orientations with a preference for water-assisted hydrogen bonding at either cationic end. The complex structures of AAATTT, ATAT, ATATAT, and AAAA are bound in a singular orientation. Analysis of the helical parameters shows a minor groove expansion of about 1 Å across all the nonalternating DNA complexes. The results from this systematic approach will convey a greater understanding of the specific recognition of a diverse array of AT-rich sequences by small molecules and more insight into the design of small molecules with enhanced specificity to AT and mixed DNA sequences.
Collapse
Affiliation(s)
- Edwin
N. Ogbonna
- Department
of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30303-3083, United States
| | - Ananya Paul
- Department
of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30303-3083, United States
| | - Abdelbasset A. Farahat
- Department
of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30303-3083, United States
- Department
of Pharmaceutical Organic Chemistry, Faculty of Pharmacy, Mansoura University, Mansoura 35516, Egypt
- Master
of Pharmaceutical Sciences Program, California
North State University, 9700 W Taron Dr., Elk Grove, California 95757, United States
| | - J. Ross Terrell
- Department
of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30303-3083, United States
| | - Ekaterina Mineva
- Department
of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30303-3083, United States
| | - Victor Ogbonna
- Department
of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30303-3083, United States
| | - David W Boykin
- Department
of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30303-3083, United States
| | - W. David Wilson
- Department
of Chemistry and Center for Diagnostics and Therapeutics, Georgia State University, Atlanta, Georgia 30303-3083, United States
| |
Collapse
|
2
|
Bogomolov A, Filonov S, Chadaeva I, Rasskazov D, Khandaev B, Zolotareva K, Kazachek A, Oshchepkov D, Ivanisenko VA, Demenkov P, Podkolodnyy N, Kondratyuk E, Ponomarenko P, Podkolodnaya O, Mustafin Z, Savinkova L, Kolchanov N, Tverdokhleb N, Ponomarenko M. Candidate SNP Markers Significantly Altering the Affinity of TATA-Binding Protein for the Promoters of Human Hub Genes for Atherogenesis, Atherosclerosis and Atheroprotection. Int J Mol Sci 2023; 24:ijms24109010. [PMID: 37240358 DOI: 10.3390/ijms24109010] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2023] [Revised: 05/13/2023] [Accepted: 05/17/2023] [Indexed: 05/28/2023] Open
Abstract
Atherosclerosis is a systemic disease in which focal lesions in arteries promote the build-up of lipoproteins and cholesterol they are transporting. The development of atheroma (atherogenesis) narrows blood vessels, reduces the blood supply and leads to cardiovascular diseases. According to the World Health Organization (WHO), cardiovascular diseases are the leading cause of death, which has been especially boosted since the COVID-19 pandemic. There is a variety of contributors to atherosclerosis, including lifestyle factors and genetic predisposition. Antioxidant diets and recreational exercises act as atheroprotectors and can retard atherogenesis. The search for molecular markers of atherogenesis and atheroprotection for predictive, preventive and personalized medicine appears to be the most promising direction for the study of atherosclerosis. In this work, we have analyzed 1068 human genes associated with atherogenesis, atherosclerosis and atheroprotection. The hub genes regulating these processes have been found to be the most ancient. In silico analysis of all 5112 SNPs in their promoters has revealed 330 candidate SNP markers, which statistically significantly change the affinity of the TATA-binding protein (TBP) for these promoters. These molecular markers have made us confident that natural selection acts against underexpression of the hub genes for atherogenesis, atherosclerosis and atheroprotection. At the same time, upregulation of the one for atheroprotection promotes human health.
Collapse
Affiliation(s)
- Anton Bogomolov
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Sergey Filonov
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
- The Natural Sciences Department, Novosibirsk State University, Novosibirsk 630090, Russia
| | - Irina Chadaeva
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Dmitry Rasskazov
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Bato Khandaev
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
- The Natural Sciences Department, Novosibirsk State University, Novosibirsk 630090, Russia
| | - Karina Zolotareva
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
- The Natural Sciences Department, Novosibirsk State University, Novosibirsk 630090, Russia
| | - Anna Kazachek
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
- The Natural Sciences Department, Novosibirsk State University, Novosibirsk 630090, Russia
| | - Dmitry Oshchepkov
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Vladimir A Ivanisenko
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Pavel Demenkov
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Nikolay Podkolodnyy
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
- Institute of Computational Mathematics and Mathematical Geophysics, Novosibirsk 630090, Russia
| | - Ekaterina Kondratyuk
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Petr Ponomarenko
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Olga Podkolodnaya
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Zakhar Mustafin
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Ludmila Savinkova
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Nikolay Kolchanov
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Natalya Tverdokhleb
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| | - Mikhail Ponomarenko
- Institute of Cytology and Genetics, Siberian Branch of Russian Academy of Sciences (SB RAS), Novosibirsk 630090, Russia
| |
Collapse
|
3
|
Izadi ME, Sabzyan H. Reactive Molecular Dynamics Simulation of the Structural Damages of the B-DNA Induced by the Oxidation/Nitration of Guanine. J Phys Chem B 2022; 126:10347-10359. [PMID: 36448964 DOI: 10.1021/acs.jpcb.2c05151] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
Reactive molecular dynamics simulations (RMD) have been carried out to investigate structural alterations of the dodecamer double-strand B-DNA due to the oxidation/nitration modifications introduced to its guanine bases, including 8-oxoguanine, 8-nitroguanine, and 5-guanidino-4-nitroimidazole, considering two distribution patterns. These modifications may arise in the case of cancer treatment using oxidative/nitrosative reactive nitrogen species as anticancer agents. Results show that these mutations affect structural characteristics of the B-DNA dodecamer in the order 8-nitroguanine > 5-guanidino-4-nitroimidazole ≫ 8-oxoguanine. For instance, the base-pair per turn for these modified B-DNA are changed respectively to 9.79, 10.88 and 10.58 from 10.51 in the native defect-free B-DNA, which is compatible with the experimental value of 10.10. In addition, these mutations allow more water molecules to diffuse into the dodecamer structure and consequently increase the possibility of the penetration of reactive and nonreactive species toward constituting nucleic base-pairs. The largest variation of the B-DNA structure is observed for the mutated B-DNA with 8-nitroguanine modifications applied to its separated CG base-pairs along the dodecamer chain. The structural changes introduced by these nitro-/oxo-modified guanine bases can be considered as a critical step in the damage of the DNA structure and alterations of its function.
Collapse
Affiliation(s)
| | - Hassan Sabzyan
- Department of Chemistry, University of Isfahan, Isfahan81746-873441, I. R. Iran
| |
Collapse
|
4
|
Gordon-Kim C, Rha A, Poppitz GA, Smith-Carpenter J, Luu R, Roberson AB, Conklin R, Blake A, Lynn DG. Polyanion order controls liquid-to-solid phase transition in peptide/nucleic acid co-assembly. Front Mol Biosci 2022; 9:991728. [PMID: 36452451 PMCID: PMC9702359 DOI: 10.3389/fmolb.2022.991728] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Accepted: 10/25/2022] [Indexed: 01/06/2024] Open
Abstract
The Central Dogma highlights the mutualistic functions of protein and nucleic acid biopolymers, and this synergy appears prominently in the membraneless organelles widely distributed throughout prokaryotic and eukaryotic organisms alike. Ribonucleoprotein granules (RNPs), which are complex coacervates of RNA with proteins, are a prime example of these membranelles organelles and underly multiple essential cellular functions. Inspired by the highly dynamic character of these organelles and the recent studies that ATP both inhibits and templates phase separation of the fused in sarcoma (FUS) protein implicated in several neurodegenerative diseases, we explored the RNA templated ordering of a single motif of the Aβ peptide of Alzheimer's disease. We now know that this strong cross-β propensity motif alone assembles through a liquid-like coacervate phase that can be externally templated to form distinct supramolecular assemblies. Now we provide evidence that structured phosphates, ranging from complex structures like double stranded and quadraplex DNA to simple trimetaphosphate, differentially impact the liquid to solid phase transition necessary for paracrystalline assembly. The results from this simple model illustrate the potential of ordered environmental templates in the transition to potentially irreversible pathogenic assemblies and provides insight into the ordering dynamics necessary for creating functional synthetic polymer co-assemblies.
Collapse
Affiliation(s)
| | - Allisandra Rha
- Children’s Health of Orange County, Research Institute, Orange, CA, United States
| | - George A. Poppitz
- Department of Chemistry, Emory University, Atlanta, GA, United States
| | | | - Regina Luu
- Department of Chemistry, Emory University, Atlanta, GA, United States
| | | | - Russell Conklin
- Department of Chemistry, Emory University, Atlanta, GA, United States
| | - Alexis Blake
- Department of Chemistry, Emory University, Atlanta, GA, United States
| | - David G. Lynn
- Department of Chemistry, Emory University, Atlanta, GA, United States
- Department of Biology, Emory University, Atlanta, GA, United States
| |
Collapse
|
5
|
Georgakopoulou C, Thomos D, Tsolis T, Ypsilantis K, Plakatouras JC, Kordias D, Magklara A, Kouderis C, Kalampounias AG, Garoufis A. Synthesis, characterization, interactions with the DNA duplex dodecamer d(5'-CGCGAATTCGCG-3') 2 and cytotoxicity of binuclear η 6-arene-Ru(II) complexes. Dalton Trans 2022; 51:13808-13825. [PMID: 36039685 DOI: 10.1039/d2dt02304k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The novel binuclear η6-arene-Ru(II) complexes with the general formula {[(η6-cym)Ru(L)]2(μ-BL)}(PF6)4, and their corresponding water soluble {[(η6-cym)Ru(L)]2(μ-BL)}Cl4, where cym = p-cymene, L = 2,2'-bipyridine (bpy) and 1,10-phenanthroline (phen), BL = 4,4'-bipyridine (BL-1), 1,2-bis(4-pyridyl)ethane (BL-2) and 1,3-bis(4-pyridyl)propane (BL-3), were synthesized and characterized. The structure of {[(η6-cym)Ru(phen)]2(μ-BL-1)}(PF6)4 was determined by X-ray single crystal methods. The interaction of {[(η6-cym)Ru(phen)]2(μ-BL-i)}Cl4 (i = 1, 2, 3; (4), (5) and (6) correspondingly) with the DNA duplex d(5'-CGCGAATTCGCG-3')2 was studied by means of NMR techniques and fluorescence titrations. The results show that complex (4) binds with a Kb = 12.133 × 103 M-1 through both intercalation and groove binding, while (5) and (6) are groove binders (Kb = 2.333 × 103 M-1 and Kb = 3.336 × 103 M-1 correspondingly). Comparison with the mononuclear complex [(η6-cym)Ru(phen)(py)]2+ reveals that it binds to the d(5'-CGCGAATTCGCG-3')2 with a Kb value two orders of magnitude lower than (4) (Kb = 0.158 × 103 M-1), indicating that for the binuclear complexes both ruthenium moieties participate in the binding. The complexes were found to be cytotoxic against the A2780 and A2780 res. cancer cell line with a selectivity index (SI) in the range of 3.0-5.9.
Collapse
Affiliation(s)
| | - Dimitrios Thomos
- Department of Chemistry, University of Ioannina, GR-45110 Ioannina, Greece.
| | - Theodoros Tsolis
- Department of Chemistry, University of Ioannina, GR-45110 Ioannina, Greece.
| | | | - John C Plakatouras
- Department of Chemistry, University of Ioannina, GR-45110 Ioannina, Greece. .,University Research Center of Ioannina (URCI), Institute of Materials Science and Computing, Ioannina, Greece
| | - Dimitris Kordias
- Biomedical Research Institute-Foundation for Research and Technology, 45110 Ioannina, Greece.,Laboratory of Clinical Chemistry, Faculty of Medicine, University of Ioannina, 45110 Ioannina, Greece
| | - Angeliki Magklara
- Biomedical Research Institute-Foundation for Research and Technology, 45110 Ioannina, Greece.,Laboratory of Clinical Chemistry, Faculty of Medicine, University of Ioannina, 45110 Ioannina, Greece.,Institute of Biosciences, University Research Center of Ioannina (U.R.C.I.), Ioannina, Greece
| | | | - Angelos G Kalampounias
- Department of Chemistry, University of Ioannina, GR-45110 Ioannina, Greece. .,University Research Center of Ioannina (URCI), Institute of Materials Science and Computing, Ioannina, Greece
| | - Achilleas Garoufis
- Department of Chemistry, University of Ioannina, GR-45110 Ioannina, Greece. .,University Research Center of Ioannina (URCI), Institute of Materials Science and Computing, Ioannina, Greece
| |
Collapse
|
6
|
Kumar Bandela A, Sadihov‐Hanoch H, Cohen‐Luria R, Gordon C, Blake A, Poppitz G, Lynn DG, Ashkenasy G. The Systems Chemistry of Nucleic‐acid‐Peptide Networks. Isr J Chem 2022. [DOI: 10.1002/ijch.202200030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Anil Kumar Bandela
- Department of Chemistry Ben-Gurion University of the Negev Beer-Sheva 84105 Israel
| | - Hava Sadihov‐Hanoch
- Department of Chemistry Ben-Gurion University of the Negev Beer-Sheva 84105 Israel
| | - Rivka Cohen‐Luria
- Department of Chemistry Ben-Gurion University of the Negev Beer-Sheva 84105 Israel
| | - Christella Gordon
- Chemistry and Biology Emory University 1521 Dickey Drive NE Atlanta GA 30322 USA
| | - Alexis Blake
- Chemistry and Biology Emory University 1521 Dickey Drive NE Atlanta GA 30322 USA
| | - George Poppitz
- Chemistry and Biology Emory University 1521 Dickey Drive NE Atlanta GA 30322 USA
| | - David G. Lynn
- Chemistry and Biology Emory University 1521 Dickey Drive NE Atlanta GA 30322 USA
| | - Gonen Ashkenasy
- Department of Chemistry Ben-Gurion University of the Negev Beer-Sheva 84105 Israel
| |
Collapse
|
7
|
Grzechowiak M, Ruszkowska A, Sliwiak J, Urbanowicz A, Jaskolski M, Ruszkowski M. New aspects of DNA recognition by group II WRKY transcription factor revealed by structural and functional study of AtWRKY18 DNA binding domain. Int J Biol Macromol 2022; 213:589-601. [PMID: 35660042 DOI: 10.1016/j.ijbiomac.2022.05.186] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2022] [Revised: 05/18/2022] [Accepted: 05/29/2022] [Indexed: 01/28/2023]
Abstract
WRKY transcription factors (TFs) constitute one of the largest families of plant TFs. Based on the organization of domains and motifs, WRKY TFs are divided into three Groups (I-III). The WRKY subgroup IIa includes three representatives in A. thaliana, AtWRKY18, AtWRKY40, and AtWRKY60, that participate in biotic and abiotic stress responses. Here we present crystal structures of the DNA binding domain (DBD) of AtWRKY18 alone and in the complex with a DNA duplex containing the WRKY-recognition sequence, W-box. Subgroup IIa WRKY TFs are known to form homo and heterodimers. Our data suggest that the dimerization interface of the full-length AtWRKY18 involves contacts between the DBD subunits. DNA binding experiments and structural analysis point out novel aspects of DNA recognition by WRKY TFs. In particular, AtWRKY18-DBD preferentially binds an overlapping tandem of W-boxes accompanied by a quasi-W-box motif. The binding of DNA deforms the B-type double helix, which suggests that the DNA fragment must be prone to form a specific structure. This can explain why despite the short W-box consensus, WRKY TFs can precisely control gene expression. Finally, this first experimental structure of a Group II WRKY TF allowed us to compare Group I-III representatives.
Collapse
Affiliation(s)
- Marta Grzechowiak
- Department of Structural Biology of Eukaryotes, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Poznan 61-704, Poland
| | - Agnieszka Ruszkowska
- Department of Structural Chemistry and Biology of Nucleic Acids, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Poznan 61-704, Poland
| | - Joanna Sliwiak
- Laboratory of Protein Engineering, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Poznan 61-704, Poland
| | - Anna Urbanowicz
- Laboratory of Protein Engineering, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Poznan 61-704, Poland
| | - Mariusz Jaskolski
- Department of Structural Biology of Eukaryotes, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Poznan 61-704, Poland; Department of Crystallography, Faculty of Chemistry, A. Mickiewicz University, Poznan 61-614, Poland
| | - Milosz Ruszkowski
- Department of Structural Biology of Eukaryotes, Institute of Bioorganic Chemistry, Polish Academy of Sciences, Poznan 61-704, Poland.
| |
Collapse
|
8
|
Abstract
B-DNA, the informational molecule for life on earth, appears to contain ratios structured around the irrational number 1.618…, often known as the “golden ratio”. This occurs in the ratio of the length:width of one turn of the helix; the ratio of the spacing of the two helices; and in the axial structure of the molecule which has ten-fold rotational symmetry. That this occurs in the information-carrying molecule for life is unexpected, and suggests the action of some process. What this process might be is unclear, but it is central to any understanding of the formation of DNA, and so life.
Collapse
|
9
|
Dohnalová H, Lankaš F. Deciphering the mechanical properties of
B‐DNA
duplex. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1575] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Hana Dohnalová
- Department of Informatics and Chemistry University of Chemistry and Technology Prague Praha 6 Czech Republic
| | - Filip Lankaš
- Department of Informatics and Chemistry University of Chemistry and Technology Prague Praha 6 Czech Republic
| |
Collapse
|
10
|
Mishra A, Siwach P, Misra P, Dhiman S, Pandey AK, Srivastava P, Jayaram B. Intron exon boundary junctions in human genome have in-built unique structural and energetic signals. Nucleic Acids Res 2021; 49:2674-2683. [PMID: 33621338 PMCID: PMC7969029 DOI: 10.1093/nar/gkab098] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/06/2019] [Revised: 01/21/2021] [Accepted: 02/22/2021] [Indexed: 11/13/2022] Open
Abstract
Precise identification of correct exon–intron boundaries is a prerequisite to analyze the location and structure of genes. The existing framework for genomic signals, delineating exon and introns in a genomic segment, seems insufficient, predominantly due to poor sequence consensus as well as limitations of training on available experimental data sets. We present here a novel concept for characterizing exon–intron boundaries in genomic segments on the basis of structural and energetic properties. We analyzed boundary junctions on both sides of all the exons (3 28 368) of protein coding genes from human genome (GENCODE database) using 28 structural and three energy parameters. Study of sequence conservation at these sites shows very poor consensus. It is observed that DNA adopts a unique structural and energy state at the boundary junctions. Also, signals are somewhat different for housekeeping and tissue specific genes. Clustering of 31 parameters into four derived vectors gives some additional insights into the physical mechanisms involved in this biological process. Sites of structural and energy signals correlate well to the positions playing important roles in pre-mRNA splicing.
Collapse
Affiliation(s)
- Akhilesh Mishra
- Supercomputing Facility for Bioinformatics & Computational Biology, Indian Institute of Technology Delhi, India.,Kusuma School of Biological Sciences, Indian Institute of Technology, Delhi, India
| | - Priyanka Siwach
- Supercomputing Facility for Bioinformatics & Computational Biology, Indian Institute of Technology Delhi, India.,Department of Biotechnology, Chaudhary Devi Lal University, Sirsa, Haryana, India
| | - Pallavi Misra
- Supercomputing Facility for Bioinformatics & Computational Biology, Indian Institute of Technology Delhi, India
| | - Simran Dhiman
- Supercomputing Facility for Bioinformatics & Computational Biology, Indian Institute of Technology Delhi, India
| | | | - Parul Srivastava
- Supercomputing Facility for Bioinformatics & Computational Biology, Indian Institute of Technology Delhi, India
| | - B Jayaram
- Supercomputing Facility for Bioinformatics & Computational Biology, Indian Institute of Technology Delhi, India.,Kusuma School of Biological Sciences, Indian Institute of Technology, Delhi, India.,Department of Chemistry, Indian Institute of Technology, Delhi, India
| |
Collapse
|
11
|
Westwood MN, Ljunggren KD, Boyd B, Becker J, Dwyer TJ, Meints GA. Single-Base Lesions and Mismatches Alter the Backbone Conformational Dynamics in DNA. Biochemistry 2021; 60:873-885. [PMID: 33689312 DOI: 10.1021/acs.biochem.0c00784] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
DNA damage has been implicated in numerous human diseases, particularly cancer, and the aging process. Single-base lesions and mismatches in DNA can be cytotoxic or mutagenic and are recognized by a DNA glycosylase during the process of base excision repair. Altered local dynamics and conformational properties in damaged DNAs have previously been suggested to assist in recognition and specificity. Herein, we use solution nuclear magnetic resonance to quantify changes in BI-BII backbone conformational dynamics due to the presence of single-base lesions in DNA, including uracil, dihydrouracil, 1,N6-ethenoadenine, and T:G mismatches. Stepwise changes to the %BII and ΔG of the BI-BII dynamic equilibrium compared to those of unmodified sequences were observed. Additionally, the equilibrium skews toward endothermicity for the phosphates nearest the lesion/mismatched base pair. Finally, the phosphates with the greatest alterations correlate with those most relevant to the repair of enzyme binding. All of these results suggest local conformational rearrangement of the DNA backbone may play a role in lesion recognition by repair enzymes.
Collapse
Affiliation(s)
- M N Westwood
- Department of Chemistry, Missouri State University, 901 South National Avenue, Springfield, Missouri 65897, United States
| | - K D Ljunggren
- Department of Chemistry, Missouri State University, 901 South National Avenue, Springfield, Missouri 65897, United States
| | - Benjamin Boyd
- Department of Chemistry, Missouri State University, 901 South National Avenue, Springfield, Missouri 65897, United States
| | - Jaclyn Becker
- Department of Chemistry, Missouri State University, 901 South National Avenue, Springfield, Missouri 65897, United States
| | - Tammy J Dwyer
- Department of Chemistry and Biochemistry, University of San Diego, San Diego, California 92110, United States
| | - Gary A Meints
- Department of Chemistry, Missouri State University, 901 South National Avenue, Springfield, Missouri 65897, United States
| |
Collapse
|
12
|
Beyond the double helix: DNA structural diversity and the PDB. J Biol Chem 2021; 296:100553. [PMID: 33744292 PMCID: PMC8063756 DOI: 10.1016/j.jbc.2021.100553] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2020] [Revised: 01/15/2021] [Accepted: 03/16/2021] [Indexed: 12/11/2022] Open
Abstract
The determination of the double helical structure of DNA in 1953 remains the landmark event in the development of modern biological and biomedical science. This structure has also been the starting point for the determination of some 2000 DNA crystal structures in the subsequent 68 years. Their structural diversity has extended to the demonstration of sequence-dependent local structure in duplex DNA, to DNA bending in short and long sequences and in the DNA wound round the nucleosome, and to left-handed duplex DNAs. Beyond the double helix itself, in circumstances where DNA sequences are or can be induced to unwind from being duplex, a wide variety of topologies and forms can exist. Quadruplex structures, based on four-stranded cores of stacked G-quartets, are prevalent though not randomly distributed in the human and other genomes and can play roles in transcription, translation, and replication. Yet more complex folds can result in DNAs with extended tertiary structures and enzymatic/catalytic activity. The Protein Data Bank is the depository of all these structures, and the resource where structures can be critically examined and validated, as well as compared one with another to facilitate analysis of conformational and base morphology features. This review will briefly survey the major structural classes of DNAs and illustrate their significance, together with some examples of how the use of the Protein Data Bank by for example, data mining, has illuminated DNA structural concepts.
Collapse
|
13
|
Zerze GH, Stillinger FH, Debenedetti PG. Thermodynamics of DNA Hybridization from Atomistic Simulations. J Phys Chem B 2021; 125:771-779. [PMID: 33434025 DOI: 10.1021/acs.jpcb.0c09237] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
Studying DNA hybridization equilibrium at atomistic length scales, either via molecular dynamics (MD) or through commonly used advanced sampling approaches, is notoriously difficult. In this work, we describe an order-parameter-based advanced sampling technique to calculate the free energy of hybridization, and estimate the melting temperature of DNA oligomers at atomistic resolution. The free energy landscapes are reported as a function of a native-topology-based order parameter for the Drew-Dickerson dodecamer and for a range of DNA decamer sequences of different GC content. Our estimated melting temperatures match the experimental numbers within ±15 °C. As a test of the numerical reliability of the procedures employed, it was verified that the predicted free energy surfaces and melting temperatures of the d- and l-enantiomers of the Drew-Dickerson dodecamer were indistinguishable within numerical accuracy.
Collapse
Affiliation(s)
- Gül H Zerze
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| | - Frank H Stillinger
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Pablo G Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
14
|
Jawad B, Poudel L, Podgornik R, Steinmetz NF, Ching WY. Molecular mechanism and binding free energy of doxorubicin intercalation in DNA. Phys Chem Chem Phys 2019; 21:3877-3893. [PMID: 30702122 DOI: 10.1039/c8cp06776g] [Citation(s) in RCA: 54] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The intercalation process of binding doxorubicin (DOX) in DNA is studied by extensive MD simulations. Many molecular factors that control the binding affinity of DOX to DNA to form a stable complex are inspected and quantified by employing continuum solvation models for estimating the binding free energy. The modified MM-PB(GB)SA methodology provides a complete energetic profile of ΔGele, ΔGvDW, ΔGpolar, ΔGnon-polar, TΔStotal, ΔGdeform, ΔGcon, and ΔGion. To identify the sequence specificity of DOX, two different DNA sequences, d(CGATCG) or DNA1 and d(CGTACG) or DNA2, with one molecule (1 : 1 complex) or two molecule (2 : 1 complex) configurations of DOX were selected in this study. Our results show that the DNA deformation energy (ΔGdeform), the energy cost from translational and rotational entropic contributions (TΔStran+rot), the total electrostatic interactions (ΔGpolar-PB/GB + ΔGele) of incorporation, the intramolecular electrostatic interactions (ΔGele) and electrostatic polar solvation interactions (ΔGpolar-PB/GB) are all unfavorable to the binding of DOX to DNA. However, they are overcome by at least five favorable interactions: the van der Waals interactions (ΔGvDW), the non-polar solvation interaction (ΔGnon-polar), the vibrational entropic contribution (TΔSvib), and the standard concentration dependent free energies of DOX (ΔGcon) and the ionic solution (ΔGion). Specifically, the van der Waals interaction appears to be the major driving force to form a stable DOX-DNA complex. We also predict that DOX has stronger binding to DNA1 than DNA2. The DNA deformation penalty and entropy cost in the 2 : 1 complex are less than those in the 1 : 1 complex, thus they indicate that the 2 : 1 complex is more stable than the 1 : 1 complex. We have calculated the total binding free energy (BFE) (ΔGt-sim) using both MM-PBSA and MM-GBSA methods, which suggests a more stable DOX-DNA complex at lower ionic concentration. The calculated BFE from the modified MM-GBSA method for DOX-DNA1 and DOX-DNA2 in the 1 : 1 complex is -9.1 and -5.1 kcal mol-1 respectively. The same quantities from the modified MM-PBSA method are -12.74 and -8.35 kcal mol-1 respectively. The value of the total BFE ΔGt-sim in the 1 : 1 complex is in reasonable agreement with the experimental value of -7.7 ± 0.3 kcal mol-1.
Collapse
Affiliation(s)
- Bahaa Jawad
- Department of Physics and Astronomy, University of Missouri-Kansas City, Kansas City, Missouri 64110, USA.
| | | | | | | | | |
Collapse
|
15
|
Dou X, Meints GA, Sedaghat-Herati R. New Insights into the Interactions of a DNA Oligonucleotide with mPEGylated-PAMAM by Circular Dichroism and Solution NMR. J Phys Chem B 2019; 123:666-674. [PMID: 30562015 DOI: 10.1021/acs.jpcb.8b08517] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Dendrimers are well-defined, highly branched, synthetic three-dimensional molecules with a large number of reactive end groups. PAMAM dendrimers form stable complexes with DNA chemistries and constitute an important class of nonviral, cationic vectors in gene delivery. The aim of this study is to examine the interactions of a 12 bp DNA oligonucletide with PAMAM-G2 and mPEG- b-PAMAM-G3 having eight surface amine groups under physiological conditions, using constant DNA concentration but varying dendrimer concentration. 1D 31P NMR, 2D NOESY, and CD spectroscopic methods were employed to investigate the interactions between the dendrimer and the DNA. The CD experiments carried out with a constant DNA concentration of 25 μM and dendrimer concentrations from 0 to 100 μM indicated minimal change to the chirality of the DNA for both types of dendrimers. While the PAMAM-G2 dendrimer caused aggregation of the majority of the DNA, the 2D NMR data of the DNA with an mPEG- b-PAMAM-G3 dendrimer indicated general broadening of the 1D 31P peaks from the DNA phosphates, a small number of 1H chemical shift perturbations (CSPs), and reduction of specific 1H-1H NOE intensities. These data suggest there is minimal structural alteration of the DNA in the complex and indicate preferential binding of the dendrimer to the central AATT region of the DNA sequence. The results herein are the first such results demonstrating a soluble DNA complex with the mPEG- b-PAMAM-G3 dendrimer analyzed by multidimensional NMR.
Collapse
Affiliation(s)
- Xiaozheng Dou
- Department of Chemistry , Missouri State University , Springfield , Missouri 65897 , United States
| | - Gary A Meints
- Department of Chemistry , Missouri State University , Springfield , Missouri 65897 , United States
| | - Reza Sedaghat-Herati
- Department of Chemistry , Missouri State University , Springfield , Missouri 65897 , United States
| |
Collapse
|
16
|
Asanuma H, Murayama K, Kamiya Y, Kashida H. The DNA Duplex as an Aqueous One-Dimensional Soft Crystal Scaffold for Photochemistry. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2018. [DOI: 10.1246/bcsj.20180278] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Hiroyuki Asanuma
- Department of Biomolecular Engineering, Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8603, Japan
| | - Keiji Murayama
- Department of Biomolecular Engineering, Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8603, Japan
| | - Yukiko Kamiya
- Department of Biomolecular Engineering, Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8603, Japan
| | - Hiromu Kashida
- Department of Biomolecular Engineering, Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8603, Japan
| |
Collapse
|
17
|
Kashida H, Kawai H, Maruyama R, Kokubo Y, Araki Y, Wada T, Asanuma H. Quantitative evaluation of energy migration between identical chromophores enabled by breaking symmetry. Commun Chem 2018. [DOI: 10.1038/s42004-018-0093-0] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
|
18
|
Nishio T, Yoshikawa Y, Fukuda W, Umezawa N, Higuchi T, Fujiwara S, Imanaka T, Yoshikawa K. Branched-Chain Polyamine Found in Hyperthermophiles Induces Unique Temperature-Dependent Structural Changes in Genome-Size DNA. Chemphyschem 2018; 19:2299-2304. [PMID: 29931720 PMCID: PMC6175440 DOI: 10.1002/cphc.201800396] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2018] [Indexed: 12/31/2022]
Abstract
A pentavalent branched‐chain polyamine, N4‐bis(aminopropyl)spermidine 3(3)(3)4, is a unique polycation found in the hyperthermophilic archaeon Thermococcus kodakarensis, which grows at temperatures between 60 and 100 °C. We studied the effects of this branched‐chain polyamine on DNA structure at different temperatures up to 80 °C. Atomic force microscopic observation revealed that 3(3)(3)4 induces a mesh‐like structure on a large DNA (166 kbp) at 24 °C. With an increase in temperature, DNA molecules tend to unwind, and multiple nano‐loops with a diameter of 10–50 nm are generated along the DNA strand at 80 °C. These results were compared to those obtained with linear‐chain polyamines, homocaldopentamine 3334 and spermidine, the former of which is a structural isomer of 3(3)(3)4. These specific effects are expected to neatly concern with its role on high‐temperature preference in hyperthermophiles.
Collapse
Affiliation(s)
- Takashi Nishio
- Faculty of Life and Medical Sciences, Doshisha University, Kyotanabe, 610-0394, Japan
| | - Yuko Yoshikawa
- Faculty of Life and Medical Sciences, Doshisha University, Kyotanabe, 610-0394, Japan
| | - Wakao Fukuda
- School of Science and Technology, Kwansei-gakuin University, Sanda, 669-1337, Japan
| | - Naoki Umezawa
- Graduate School of Pharmaceutical Sciences, Nagoya City University, Nagoya, 467-8603, Japan
| | - Tsunehiko Higuchi
- Graduate School of Pharmaceutical Sciences, Nagoya City University, Nagoya, 467-8603, Japan
| | - Shinsuke Fujiwara
- School of Science and Technology, Kwansei-gakuin University, Sanda, 669-1337, Japan
| | - Tadayuki Imanaka
- Research Organization of Science and Technology, Ritsumeikan University, Kusatsu, 525-8577, Japan
| | - Kenichi Yoshikawa
- Faculty of Life and Medical Sciences, Doshisha University, Kyotanabe, 610-0394, Japan
| |
Collapse
|
19
|
Minguirbara A, Nsangou M. DFT study of geometrical and vibrational features of a 3',5'-deoxydisugar-monophosphate (dDSMP) DNA model in the presence of counterions and solvent. J Mol Model 2018. [PMID: 29516189 DOI: 10.1007/s00894-018-3629-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The B3LYP/6-31++G* theoretical level was used to study the influence of various hexahydrated monovalent (Li+, Na+, K+) and divalent (Mg2+) metal counterions in interaction with the charged PO2- group, on the geometrical and vibrational characteristics of the DNA fragments of 3',5'-dDSMP, represented by four conformers (g+g+, g+t, g-g- and g-t). All complexes were optimized through two solvation models [the explicit model (6H2O) and the hybrid model (6H2O/Continuum)]. The results obtained established that, in the hybrid model, counterions (Li+, Na+, K+, Mg2+) always remain in the bisector plane of the O1-P-O2 angle. When these counterions are explicitly hydrated, the smallest counterions (Li+, Na+) deviate from the bisector plane, while the largest counterions (K+ and Mg2+) always remain in the same plane. On the other hand, the present calculations reveal that the g+g+ conformer is the most stable in the presence of monovalent counterions, while conformers g+t and g-t are the most stable in the presence of the divalent counterion Mg2+. Finally, the hybrid solvation model seems to be in better agreement with the available crystallographic and spectroscopic (Raman) experiments than the explicit model. Indeed, the six conformational torsions of the C4'-C3'-O3'-PO-2-O5'-C5'-C4' segment of all complexes of the g-g- conformer in 6H2O/Continuum remain similar to the available experimental data of A- and B-DNA forms. The calculated wavenumbers of the g+g+ conformer in the presence of the monovalent counterion and of g-t conformer in presence of the divalent counterion in the hybrid model are in good agreement with the Raman experimental data of A- and B-DNA forms. In addition, the maximum deviation between the calculated wavenumbers in the 6H2O/Continuum for the g+g+ conformer and experimental value measured in an aqueous solution of the DMP-Na+ complex, is <1.07% for the PO2- (asymmetric and symmetric) stretching modes and <2.03% for the O5'-C5' and O3'-C3' stretching modes. Graphical abstract dDSMP-(OO)- Mg2+/6W/Continuum.
Collapse
Affiliation(s)
- Alain Minguirbara
- Department of Physics, Higher Teachers' Training College, University of Maroua, PO Box 46, Maroua, Cameroon
| | - Mama Nsangou
- Department of Physics, Higher Teachers' Training College, University of Maroua, PO Box 46, Maroua, Cameroon.
| |
Collapse
|
20
|
O' Lee DJ. Introducing a model of pairing based on base pair specific interactions between identical DNA sequences. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:075102. [PMID: 29219116 DOI: 10.1088/1361-648x/aaa043] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
At present, there have been suggested two types of physical mechanism that may facilitate preferential pairing between DNA molecules, with identical or similar base pair texts, without separation of base pairs. One mechanism solely relies on base pair specific patterns of helix distortion being the same on the two molecules, discussed extensively in the past. The other mechanism proposes that there are preferential interactions between base pairs of the same composition. We introduce a model, built on this second mechanism, where both thermal stretching and twisting fluctuations are included, as well as the base pair specific helix distortions. Firstly, we consider an approximation for weak pairing interactions, or short molecules. This yields a dependence of the energy on the square root of the molecular length, which could explain recent experimental data. However, analysis suggests that this approximation is no longer valid at large DNA lengths. In a second approximation, for long molecules, we define two adaptation lengths for twisting and stretching, over which the pairing interaction can limit the accumulation of helix disorder. When the pairing interaction is sufficiently strong, both adaptation lengths are finite; however, as we reduce pairing strength, the stretching adaptation length remains finite but the torsional one becomes infinite. This second state persists to arbitrarily weak values of the pairing strength; suggesting that, if the molecules are long enough, the pairing energy scales as length. To probe differences between the two pairing mechanisms, we also construct a model of similar form. However, now, pairing between identical sequences solely relies on the intrinsic helix distortion patterns. Between the two models, we see interesting qualitative differences. We discuss our findings, and suggest new work to distinguish between the two mechanisms.
Collapse
Affiliation(s)
- Dominic J O' Lee
- Department of Chemistry, Imperial College London, SW7 2AZ, London, United Kingdom
| |
Collapse
|
21
|
Kashida H, Kurihara A, Kawai H, Asanuma H. Orientation-dependent FRET system reveals differences in structures and flexibilities of nicked and gapped DNA duplexes. Nucleic Acids Res 2017; 45:e105. [PMID: 28369626 PMCID: PMC5499647 DOI: 10.1093/nar/gkx200] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2016] [Accepted: 03/17/2017] [Indexed: 12/18/2022] Open
Abstract
Differences in structures and flexibilities of DNA duplexes play important roles on recognition by DNA-binding proteins. We herein describe a novel method for structural analyses of DNA duplexes by using orientation dependence of Förster resonance energy transfer (FRET). We first analyzed canonical B-form duplex and correct structural parameters were obtained. The experimental FRET efficiencies were in excellent agreement with values theoretically calculated by using determined parameters. We then investigated DNA duplexes with nick and gaps, which are key intermediates in DNA repair systems. Effects of gap size on structures and flexibilities were successfully revealed. Since our method is facile and sensitive, it could be widely used to analyze DNA structures containing damages and non-natural molecules.
Collapse
Affiliation(s)
- Hiromu Kashida
- Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan.,PRESTO, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Ayako Kurihara
- Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Hayato Kawai
- Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| | - Hiroyuki Asanuma
- Graduate School of Engineering, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8603, Japan
| |
Collapse
|
22
|
Lee DJO. Statistical mechanical model for a closed loop plectoneme with weak helix specific forces. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2017; 29:145101. [PMID: 28251958 DOI: 10.1088/1361-648x/aa521c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We develop a statistical mechanical framework, based on a variational approximation, to describe closed loop plectonemes. This framework incorporates weak helix structure dependent forces into the determination of the free energy and average structure of a plectoneme. Notably, due to their chiral nature, helix structure dependent forces break the symmetry between left and right handed supercoiling. The theoretical approach, presented here, also provides a systematic way of enforcing the topological constraint of closed loop supercoiling in the variational approximation. At large plectoneme lengths, by considering correlation functions in an expansion in terms of the spatial mean twist density about its thermally averaged value, it can be argued that topological constraint may be approximated by replacing twist and writhe by their thermal averages. A Lagrange multiplier, containing the sum of average twist and writhe, can be added to the free energy to conveniently inforce this result. The average writhe can be calculated through the thermal average of the Gauss' integral in the variational approximation. Furthermore, this approach allows for a possible way to calculate finite size corrections due to the topological constraint. Using interaction energy terms from the mean-field Kornyshev-Leikin theory, for parameter values that correspond to weak helix dependent forces, we calculate the free energy, fluctuation magnitudes and mean geometric parameters for the plectoneme. We see a slight asymmetry, where interestingly, left handed supercoils have a looser structure than right handed ones, although with a lower free energy, unlike what the previous ground state calculations would suggest.
Collapse
Affiliation(s)
- Dominic J O' Lee
- Department of Chemistry, Imperial College London, SW7 2AZ, London, United Kingdom
| |
Collapse
|
23
|
Winger J, Bowman GD. The Sequence of Nucleosomal DNA Modulates Sliding by the Chd1 Chromatin Remodeler. J Mol Biol 2017; 429:808-822. [PMID: 28189426 DOI: 10.1016/j.jmb.2017.02.002] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2016] [Revised: 02/03/2017] [Accepted: 02/06/2017] [Indexed: 10/20/2022]
Abstract
Chromatin remodelers are ATP-dependent enzymes that are critical for reorganizing and repositioning nucleosomes in concert with many basic cellular processes. For the chromodomain helicase DNA-binding protein 1 (Chd1) remodeler, nucleosome sliding has been shown to depend on the DNA flanking the nucleosome, transcription factor binding at the nucleosome edge, and the presence of the histone H2A/H2B dimer on the entry side. Here, we report that Chd1 is also sensitive to the sequence of DNA within the nucleosome and slides nucleosomes made with the 601 Widom positioning sequence asymmetrically. Kinetic and equilibrium experiments show that poly(dA:dT) tracts perturb remodeling reactions if within one and a half helical turns of superhelix location 2 (SHL2), where the Chd1 ATPase engages nucleosomal DNA. These sequence-dependent effects do not rely on the Chd1 DNA-binding domain and are not due to differences in nucleosome affinity. Using site-specific cross-linking, we show that internal poly(dA:dT) tracts do not block the engagement of the ATPase motor with SHL2, yet they promote multiple translational positions of DNA with respect to both Chd1 and the histone core. We speculate that Chd1 senses the sequence-dependent response of DNA as the remodeler ATPase perturbs the duplex at SHL2. These results suggest that the sequence sensitivity of histones and remodelers occur at unique segments of DNA on the nucleosome, allowing them to work together or in opposition to determine nucleosome positions throughout the genome.
Collapse
Affiliation(s)
- Jessica Winger
- Department of Biology, Johns Hopkins University, Baltimore, MD 21218, USA
| | - Gregory D Bowman
- Department of Biology, Johns Hopkins University, Baltimore, MD 21218, USA; T.C. Jenkins Department of Biophysics, Johns Hopkins University, Baltimore, MD 21218, USA.
| |
Collapse
|
24
|
Amato J, Pagano A, Cosconati S, Amendola G, Fotticchia I, Iaccarino N, Marinello J, De Magis A, Capranico G, Novellino E, Pagano B, Randazzo A. Discovery of the first dual G-triplex/G-quadruplex stabilizing compound: a new opportunity in the targeting of G-rich DNA structures? Biochim Biophys Acta Gen Subj 2016; 1861:1271-1280. [PMID: 27836755 DOI: 10.1016/j.bbagen.2016.11.008] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2016] [Revised: 11/03/2016] [Accepted: 11/05/2016] [Indexed: 11/17/2022]
Abstract
BACKGROUND Guanine-rich DNA motifs can form non-canonical structures known as G-quadruplexes, whose role in tumorigenic processes makes them attractive drug-target candidates for cancer therapy. Recent studies revealed that the folding and unfolding pathways of G-quadruplexes proceed through a quite stable intermediate named G-triplex. METHODS Virtual screening was employed to identify a small set of putative G-triplex ligands. The G-triplex stabilizing properties of these compounds were analyzed by CD melting assay. DSC, non-denaturing gel electrophoresis, NMR and molecular modeling studies were performed to investigate the interaction between the selected compound 1 and G-rich DNA structures. Cytotoxic activity of 1 was evaluated by MTT cell proliferation assay. RESULTS The experiments led to the identification of a promising hit that was shown to bind preferentially to G-triplex and parallel-stranded G-quadruplexes over duplex and antiparallel G-quadruplexes. Molecular modeling results suggested a partial end-stacking of 1 to the external G-triad/G-tetrads as a binding mode. Biological assays showed that 1 is endowed with cytotoxic effect on human osteosarcoma cells. CONCLUSIONS A tandem application of virtual screening along with the experimental investigation was employed to discover a G-triplex-targeting ligand. Experiments revealed that the selected compound actually acts as a dual G-triplex/G-quadruplex stabilizer, thus stimulating further studies aimed at its optimization. GENERAL SIGNIFICANCE The discovery of molecules able to bind and stabilize G-triplex structures is highly appealing, but their transient state makes challenging their recognition. These findings suggest that the identification of ligands with dual G-triplex/G-quadruplex stabilizing properties may represent a new route for the design of anticancer agents targeting the G-rich DNA structures. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.
Collapse
Affiliation(s)
- Jussara Amato
- Department of Pharmacy, University of Naples Federico II, via D. Montesano 49, 80131 Naples, Italy
| | - Alessia Pagano
- Department of Pharmacy, University of Naples Federico II, via D. Montesano 49, 80131 Naples, Italy
| | - Sandro Cosconati
- DiSTABiF, Second University of Naples, via Vivaldi 43, 81100 Caserta, Italy
| | - Giorgio Amendola
- DiSTABiF, Second University of Naples, via Vivaldi 43, 81100 Caserta, Italy
| | - Iolanda Fotticchia
- Department of Pharmacy, University of Naples Federico II, via D. Montesano 49, 80131 Naples, Italy
| | - Nunzia Iaccarino
- Department of Pharmacy, University of Naples Federico II, via D. Montesano 49, 80131 Naples, Italy
| | - Jessica Marinello
- Department of Pharmacy and Biotechnology, University of Bologna, 40126 Bologna, Italy
| | - Alessio De Magis
- Department of Pharmacy and Biotechnology, University of Bologna, 40126 Bologna, Italy
| | - Giovanni Capranico
- Department of Pharmacy and Biotechnology, University of Bologna, 40126 Bologna, Italy
| | - Ettore Novellino
- Department of Pharmacy, University of Naples Federico II, via D. Montesano 49, 80131 Naples, Italy
| | - Bruno Pagano
- Department of Pharmacy, University of Naples Federico II, via D. Montesano 49, 80131 Naples, Italy
| | - Antonio Randazzo
- Department of Pharmacy, University of Naples Federico II, via D. Montesano 49, 80131 Naples, Italy.
| |
Collapse
|
25
|
|
26
|
Deiana M, Mettra B, Matczyszyn K, Piela K, Pitrat D, Olesiak-Banska J, Monnereau C, Andraud C, Samoc M. Interactions of a biocompatible water-soluble anthracenyl polymer derivative with double-stranded DNA. Phys Chem Chem Phys 2016; 17:30318-27. [PMID: 26506886 DOI: 10.1039/c5cp05381a] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We have studied the interaction of a polymeric water soluble anthracenyl derivative () with salmon testes DNA. The results from UV-Vis, fluorescence, Fourier transform infrared (FT-IR) and circular dichroism spectroscopies indicate that the groove binding process regulates the interaction between and DNA. The binding constants, calculated by absorption spectroscopy at 298, 304 and 310 K, were equal to 3.2 × 10(5) M(-1), 4.7 × 10(5) M(-1), and 6.6 × 10(5) M(-1) respectively, proving a relatively high affinity of for salmon testes DNA. Results of Hoechst 33258 displacement assays strongly support the groove binding mode of to DNA. The association stoichiometry of the :DNA adduct was found to be 1 for every 5 base pairs. FT-IR spectra, recorded at different /DNA molar ratios, indicate the involvement of the phosphate groups and adenine and thymine DNA bases in the association process. Thermodynamic results suggest that hydrophobic forces regulate the binding of with DNA without excluding some extent of involvement of van der Waals forces and hydrogen bonding arising due to surface binding between the hydrophilic polymeric arms of the ligand and the functional groups positioned on the edge of the groove. The resulting composite biomaterial could constitute a valuable candidate for future biological and/or photonic applications.
Collapse
Affiliation(s)
- Marco Deiana
- Advanced Materials Engineering and Modelling Group, Faculty of Chemistry, Wroclaw University of Technology, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland.
| | - Bastien Mettra
- Laboratoire de Chimie, CNRS UMR 5182, Ecole Normale Supérieure de Lyon, Université, Lyon 1, Lyon, France
| | - Katarzyna Matczyszyn
- Advanced Materials Engineering and Modelling Group, Faculty of Chemistry, Wroclaw University of Technology, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland.
| | - Katarzyna Piela
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wroclaw University of Technology, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland
| | - Delphine Pitrat
- Laboratoire de Chimie, CNRS UMR 5182, Ecole Normale Supérieure de Lyon, Université, Lyon 1, Lyon, France
| | - Joanna Olesiak-Banska
- Advanced Materials Engineering and Modelling Group, Faculty of Chemistry, Wroclaw University of Technology, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland.
| | - Cyrille Monnereau
- Laboratoire de Chimie, CNRS UMR 5182, Ecole Normale Supérieure de Lyon, Université, Lyon 1, Lyon, France
| | - Chantal Andraud
- Laboratoire de Chimie, CNRS UMR 5182, Ecole Normale Supérieure de Lyon, Université, Lyon 1, Lyon, France
| | - Marek Samoc
- Advanced Materials Engineering and Modelling Group, Faculty of Chemistry, Wroclaw University of Technology, Wyb. Wyspianskiego 27, 50-370 Wroclaw, Poland.
| |
Collapse
|
27
|
Mandal PK, Collie GW, Srivastava SC, Kauffmann B, Huc I. Structure elucidation of the Pribnow box consensus promoter sequence by racemic DNA crystallography. Nucleic Acids Res 2016; 44:5936-43. [PMID: 27137886 PMCID: PMC4937325 DOI: 10.1093/nar/gkw367] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2016] [Revised: 04/05/2016] [Accepted: 04/23/2016] [Indexed: 11/13/2022] Open
Abstract
It has previously been shown that the use of racemic mixtures of naturally chiral macromolecules such as protein and DNA can significantly aid the crystallogenesis process, thereby addressing one of the major bottlenecks to structure determination by X-ray crystallographic methods-that of crystal growth. Although previous studies have provided convincing evidence of the applicability of the racemic crystallization technique to DNA through the study of well-characterized DNA structures, we sought to apply this method to a historically challenging DNA sequence. For this purpose we chose a non-self-complementary DNA duplex containing the biologically-relevant Pribnow box consensus sequence 'TATAAT'. Four racemic crystal structures of this previously un-crystallizable DNA target are reported (with resolutions in the range of 1.65-2.3 Å), with further crystallographic studies and structural analysis providing insight into the racemic crystallization process as well as structural details of this highly pertinent DNA sequence.
Collapse
Affiliation(s)
- Pradeep K Mandal
- Université de Bordeaux, CBMN (UMR5248), Pessac 33600, France CNRS, CBMN (UMR5248), Pessac 33600, France
| | - Gavin W Collie
- Université de Bordeaux, CBMN (UMR5248), Pessac 33600, France CNRS, CBMN (UMR5248), Pessac 33600, France
| | | | - Brice Kauffmann
- Université de Bordeaux, Institut Européen de Chimie et Biologie (UMS3033), Pessac 33600, France CNRS, Institut Européen de Chimie et Biologie (UMS3033), Pessac 33600, France INSERM, Institut Européen de Chimie et Biologie (US001), Pessac 33600, France
| | - Ivan Huc
- Université de Bordeaux, CBMN (UMR5248), Pessac 33600, France CNRS, CBMN (UMR5248), Pessac 33600, France
| |
Collapse
|
28
|
O' Lee DJ, Danilowicz C, Rochester C, Kornyshev AA, Prentiss M. Evidence of protein-free homology recognition in magnetic bead force-extension experiments. Proc Math Phys Eng Sci 2016; 472:20160186. [PMID: 27493568 PMCID: PMC4971244 DOI: 10.1098/rspa.2016.0186] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
Earlier theoretical studies have proposed that the homology-dependent pairing of large tracts of dsDNA may be due to physical interactions between homologous regions. Such interactions could contribute to the sequence-dependent pairing of chromosome regions that may occur in the presence or the absence of double-strand breaks. Several experiments have indicated the recognition of homologous sequences in pure electrolytic solutions without proteins. Here, we report single-molecule force experiments with a designed 60 kb long dsDNA construct; one end attached to a solid surface and the other end to a magnetic bead. The 60 kb constructs contain two 10 kb long homologous tracts oriented head to head, so that their sequences match if the two tracts fold on each other. The distance between the bead and the surface is measured as a function of the force applied to the bead. At low forces, the construct molecules extend substantially less than normal, control dsDNA, indicating the existence of preferential interaction between the homologous regions. The force increase causes no abrupt but continuous unfolding of the paired homologous regions. Simple semi-phenomenological models of the unfolding mechanics are proposed, and their predictions are compared with the data.
Collapse
Affiliation(s)
- D J O' Lee
- Department of Chemistry , Imperial College London , London SW7 2AZ, UK
| | - C Danilowicz
- Department of Physics , Harvard University, Cambridge , MA 02138, USA
| | - C Rochester
- Department of Chemistry , Imperial College London , London SW7 2AZ, UK
| | - A A Kornyshev
- Department of Chemistry , Imperial College London , London SW7 2AZ, UK
| | - M Prentiss
- Department of Physics , Harvard University, Cambridge , MA 02138, USA
| |
Collapse
|
29
|
Thirugnanasambandam A, Karthik S, Artheswari G, Gautham N. DNA polymorphism in crystals: three stable conformations for the decadeoxynucleotide d(GCATGCATGC). ACTA CRYSTALLOGRAPHICA SECTION D-STRUCTURAL BIOLOGY 2016; 72:780-8. [PMID: 27303798 DOI: 10.1107/s2059798316006306] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Received: 01/27/2016] [Accepted: 04/14/2016] [Indexed: 11/11/2022]
Abstract
High-resolution structures of DNA fragments determined using X-ray crystallography or NMR have provided descriptions of a veritable alphabet of conformations. They have also shown that DNA is a flexible molecule, with some sequences capable of adopting two different structures. Here, the first example is presented of a DNA fragment that can assume three different and distinct conformations in crystals. The decanucleotide d(GCATGCATGC) was previously reported to assume a single-stranded double-fold structure. In one of the two crystal structures described here the decamer assumes both the double-fold conformation and, simultaneously, the more conventional B-type double-helical structure. In the other crystal the sequence assumes the A-type double-helical conformation. These results, taken together with CD spectra, which were recorded as the decamer was titrated against four metal ions and spermine, indicate that the molecule may exist as a mixed population of structures in solution. Small differences in the environmental conditions, such as the concentration of metal ion, may decide which of these crystallizes out. The results also support the idea that it may be possible for DNA to change its structure to suit the binding requirements of proteins or drugs.
Collapse
Affiliation(s)
| | - Selvam Karthik
- CAS in Crystallography and Biophysics, University of Madras, Guindy Campus, Chennai 600 025, India
| | - Gunanithi Artheswari
- CAS in Crystallography and Biophysics, University of Madras, Guindy Campus, Chennai 600 025, India
| | - Namasivayam Gautham
- CAS in Crystallography and Biophysics, University of Madras, Guindy Campus, Chennai 600 025, India
| |
Collapse
|
30
|
Mukherjee S, Kailasam S, Bansal M, Bhattacharyya D. Stacking interactions in RNA and DNA: Roll-slide energy hyperspace for ten unique dinucleotide steps. Biopolymers 2016; 103:134-47. [PMID: 25257334 DOI: 10.1002/bip.22566] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2014] [Revised: 09/17/2014] [Accepted: 09/17/2014] [Indexed: 12/15/2022]
Abstract
Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by ωB97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality.
Collapse
Affiliation(s)
- Sanchita Mukherjee
- Computational Science Division, Saha Institute of Nuclear Physics, Kolkata, 700064, West Bengal, India
| | | | | | | |
Collapse
|
31
|
Waters JT, Lu XJ, Galindo-Murillo R, Gumbart JC, Kim HD, Cheatham TE, Harvey SC. Transitions of Double-Stranded DNA Between the A- and B-Forms. J Phys Chem B 2016; 120:8449-56. [PMID: 27135262 DOI: 10.1021/acs.jpcb.6b02155] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The structure of double-stranded DNA (dsDNA) is sensitive to solvent conditions. In solution, B-DNA is the favored conformation under physiological conditions, while A-DNA is the form found under low water activity. The A-form is induced locally in some protein-DNA complexes, and repeated transitions between the B- and A-forms have been proposed to generate the forces used to drive dsDNA into viral capsids during genome packaging. Here, we report analyses on previous molecular dynamics (MD) simulations on B-DNA, along with new MD simulations on the transition from A-DNA to B-DNA in solution. We introduce the A-B Index (ABI), a new metric along the A-B continuum, to quantify our results. When A-DNA is placed in an equilibrated solution at physiological ionic strength, there is no energy barrier to the transition to the B-form, which begins within about 1 ns. The transition is essentially complete within 5 ns, although occasionally a stretch of a few base pairs will remain A-like for up to ∼10 ns. A comparison of four sequences with a range of predicted A-phobicities shows that more A-phobic sequences make the transition more rapidly than less A-phobic sequences. Simulations on dsDNA with a region of roughly one turn locked in the A-form allow us to characterize the A/B junction, which has an average bend angle of 20-30°. Fluctuations in this angle occur with characteristic times of about 10 ns.
Collapse
Affiliation(s)
- James T Waters
- School of Physics, Georgia Institute of Technology , Atlanta, Georgia 30332, United States
| | - Xiang-Jun Lu
- Department of Biological Sciences, Columbia University , New York, New York 10027, United States
| | - Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, L.S. Skaggs Pharmacy Institute, University of Utah , Salt Lake City, Utah 84112, United States
| | - James C Gumbart
- School of Physics, Georgia Institute of Technology , Atlanta, Georgia 30332, United States
| | - Harold D Kim
- School of Physics, Georgia Institute of Technology , Atlanta, Georgia 30332, United States
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, L.S. Skaggs Pharmacy Institute, University of Utah , Salt Lake City, Utah 84112, United States
| | - Stephen C Harvey
- Department of Biochemistry and Molecular Biophysics, Perelman School of Medicine, University of Pennsylvania , Philadelphia, Pennsylvania 19104, United States
| |
Collapse
|
32
|
Nikolova EN, Zhou H, Gottardo FL, Alvey HS, Kimsey IJ, Al-Hashimi HM. A historical account of Hoogsteen base-pairs in duplex DNA. Biopolymers 2016; 99:955-68. [PMID: 23818176 DOI: 10.1002/bip.22334] [Citation(s) in RCA: 59] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2013] [Accepted: 06/17/2013] [Indexed: 11/05/2022]
Abstract
In 1957, a unique pattern of hydrogen bonding between N3 and O4 on uracil and N7 and N6 on adenine was proposed to explain how poly(rU) strands can associate with poly(rA)-poly(rU) duplexes to form triplexes. Two years later, Karst Hoogsteen visualized such a noncanonical A-T base-pair through X-ray analysis of co-crystals containing 9-methyladenine and 1-methylthymine. Subsequent X-ray analyses of guanine and cytosine derivatives yielded the expected Watson-Crick base-pairing, but those of adenine and thymine (or uridine) did not yield Watson-Crick base-pairs, instead favoring "Hoogsteen" base-pairing. More than two decades ensued without experimental "proof" for A-T Watson-Crick base-pairs, while Hoogsteen base-pairs continued to surface in AT-rich sequences, closing base-pairs of apical loops, in structures of DNA bound to antibiotics and proteins, damaged and chemically modified DNA, and in polymerases that replicate DNA via Hoogsteen pairing. Recently, NMR studies have shown that base-pairs in duplex DNA exist as a dynamic equilibrium between Watson-Crick and Hoogsteen forms. There is now little doubt that Hoogsteen base-pairs exist in significant abundance in genomic DNA, where they can expand the structural and functional versatility of duplex DNA beyond that which can be achieved based only on Watson-Crick base-pairing. Here, we provide a historical account of the discovery and characterization of Hoogsteen base-pairs, hoping that this will inform future studies exploring the occurrence and functional importance of these alternative base-pairs.
Collapse
Affiliation(s)
- Evgenia N Nikolova
- Department of Chemistry & Biophysics, The University of Michigan, 930 North University Avenue, Ann Arbor, MI, 48109-1055; Integrative Structural & Computational Biology Department, The Scripps Research Institute, 10550 North Torrey Pines Road, La Jolla, CA, 92037
| | | | | | | | | | | |
Collapse
|
33
|
Saoji M, Paukstelis PJ. Sequence-dependent structural changes in a self-assembling DNA oligonucleotide. ACTA CRYSTALLOGRAPHICA. SECTION D, BIOLOGICAL CRYSTALLOGRAPHY 2015; 71:2471-8. [PMID: 26627654 PMCID: PMC4667286 DOI: 10.1107/s1399004715019598] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/09/2015] [Accepted: 10/15/2015] [Indexed: 11/10/2022]
Abstract
DNA has proved to be a remarkable molecule for the construction of sophisticated two-dimensional and three-dimensional architectures because of its programmability and structural predictability provided by complementary Watson-Crick base pairing. DNA oligonucleotides can, however, exhibit a great deal of local structural diversity. DNA conformation is strongly linked to both environmental conditions and the nucleobase identities inherent in the oligonucleotide sequence, but the exact relationship between sequence and local structure is not completely understood. This study examines how a single-nucleotide addition to a class of self-assembling DNA 13-mers leads to a significantly different overall structure under identical crystallization conditions. The DNA 13-mers self-assemble in the presence of Mg(2+) through a combination of Watson-Crick and noncanonical base-pairing interactions. The crystal structures described here show that all of the predicted Watson-Crick base pairs are present, with the major difference being a significant rearrangement of noncanonical base pairs. This includes the formation of a sheared A-G base pair, a junction of strands formed from base-triple interactions, and tertiary interactions that generate structural features similar to tandem sheared G-A base pairs. The adoption of this alternate noncanonical structure is dependent in part on the sequence in the Watson-Crick duplex region. These results provide important new insights into the sequence-structure relationship of short DNA oligonucleotides and demonstrate a unique interplay between Watson-Crick and noncanonical base pairs that is responsible for crystallization fate.
Collapse
Affiliation(s)
- Maithili Saoji
- Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742, USA
| | - Paul J. Paukstelis
- Department of Chemistry and Biochemistry, University of Maryland, College Park, MD 20742, USA
- Center for Biomolecular Structure and Organisation, Maryland NanoCenter, University of Maryland, College Park, MD 20742, USA
| |
Collapse
|
34
|
Gu C, Zhang J, Yang YI, Chen X, Ge H, Sun Y, Su X, Yang L, Xie S, Gao YQ. DNA Structural Correlation in Short and Long Ranges. J Phys Chem B 2015; 119:13980-90. [PMID: 26439165 DOI: 10.1021/acs.jpcb.5b06217] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Recent single-molecule measurements have revealed the DNA allostery in protein/DNA binding. MD simulations showed that this allosteric effect is associated with the deformation properties of DNA. In this study, we used MD simulations to further investigate the mechanism of DNA structural correlation, its dependence on DNA sequence, and the chemical modification of the bases. Besides a random sequence, poly d(AT) and poly d(GC) are also used as simpler model systems, which show the different bending and twisting flexibilities. The base-stacking interactions and the methyl group on the 5-carbon site of thymine causes local structures and flexibility to be very different for the two model systems, which further lead to obviously different tendencies of the conformational deformations, including the long-range allosteric effects.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | - Sunney Xie
- Department of Chemistry and Chemical Biology, Harvard University , Cambridge, Massachusetts 02138, United States
| | | |
Collapse
|
35
|
O' Lee DJ, Wynveen A, Albrecht T, Kornyshev AA. Which way up? Recognition of homologous DNA segments in parallel and antiparallel alignments. J Chem Phys 2015; 142:045101. [PMID: 25638008 DOI: 10.1063/1.4905291] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Homologous gene shuffling between DNA molecules promotes genetic diversity and is an important pathway for DNA repair. For this to occur, homologous genes need to find and recognize each other. However, despite its central role in homologous recombination, the mechanism of homology recognition has remained an unsolved puzzle of molecular biology. While specific proteins are known to play a role at later stages of recombination, an initial coarse grained recognition step has, however, been proposed. This relies on the sequence dependence of the DNA structural parameters, such as twist and rise, mediated by intermolecular interactions, in particular, electrostatic ones. In this proposed mechanism, sequences that have the same base pair text, or are homologous, have lower interaction energy than those sequences with uncorrelated base pair texts. The difference between the two energies is termed the "recognition energy." Here, we probe how the recognition energy changes when one DNA fragment slides past another, and consider, for the first time, homologous sequences in antiparallel alignment. This dependence on sliding is termed the "recognition well." We find there is a recognition well for anti-parallel, homologous DNA tracts, but only a very shallow one, so that their interaction will differ little from the interaction between two nonhomologous tracts. This fact may be utilized in single molecule experiments specially targeted to test the theory. As well as this, we test previous theoretical approximations in calculating the recognition well for parallel molecules against MC simulations and consider more rigorously the optimization of the orientations of the fragments about their long axes upon calculating these recognition energies. The more rigorous treatment affects the recognition energy a little, when the molecules are considered rigid. When torsional flexibility of the DNA molecules is introduced, we find excellent agreement between the analytical approximation and simulations.
Collapse
Affiliation(s)
- Dominic J O' Lee
- Department of Chemistry, Imperial College London, SW7 2AZ London, United Kingdom
| | - Aaron Wynveen
- School of Physics and Astronomy, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Tim Albrecht
- Department of Chemistry, Imperial College London, SW7 2AZ London, United Kingdom
| | - Alexei A Kornyshev
- Department of Chemistry, Imperial College London, SW7 2AZ London, United Kingdom
| |
Collapse
|
36
|
Srivastava P, Abou El Asrar R, Knies C, Abramov M, Froeyen M, Rozenski J, Rosemeyer H, Herdewijn P. Achiral, acyclic nucleic acids: synthesis and biophysical studies of a possible prebiotic polymer. Org Biomol Chem 2015; 13:9249-60. [DOI: 10.1039/c5ob00898k] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The search for prebiotic, nucleic acid precursors is, at its best, a speculative undertaking.
Collapse
Affiliation(s)
| | | | - C. Knies
- Organic Materials Chemistry and Bioorganic Chemistry
- Institute or Chemistry
- University of Osnabrück
- 49069 Osnabrück
- Germany
| | - M. Abramov
- Medicinal Chemistry
- KU Leuven
- B-3000 Leuven
- Belgium
| | - M. Froeyen
- Medicinal Chemistry
- KU Leuven
- B-3000 Leuven
- Belgium
| | - J. Rozenski
- Medicinal Chemistry
- KU Leuven
- B-3000 Leuven
- Belgium
| | - H. Rosemeyer
- Organic Materials Chemistry and Bioorganic Chemistry
- Institute or Chemistry
- University of Osnabrück
- 49069 Osnabrück
- Germany
| | - P. Herdewijn
- Medicinal Chemistry
- KU Leuven
- B-3000 Leuven
- Belgium
| |
Collapse
|
37
|
Kruse H, Šponer J. Towards biochemically relevant QM computations on nucleic acids: controlled electronic structure geometry optimization of nucleic acid structural motifs using penalty restraint functions. Phys Chem Chem Phys 2014; 17:1399-410. [PMID: 25427983 DOI: 10.1039/c4cp04680c] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Recent developments in dispersion-corrected density functional theory methods allow for the first time the description of large fragments of nucleic acids (hundreds of atoms) with an accuracy clearly surpassing the accuracy of common biomolecular force fields. Such calculations can significantly improve the description of the potential energy surface of nucleic acid molecules, which may be useful for studies of molecular interactions and conformational preferences of nucleic acids, as well as verification and parameterization of other methods. The first of such studies, however, demonstrated that successful applications of accurate QM calculations to larger nucleic acid building blocks are hampered by difficulties in obtaining geometries that are biochemically relevant and are not biased by non-native structural features. We present an approach that can greatly facilitate large-scale QM studies on nucleic acids, namely electronic structure geometry optimization of nucleic acid fragments utilizing a penalty function to restrain key internal coordinates with a specific focus on the torsional backbone angles. This work explores the viability of these restraint optimizations for DFT-D3, PM6-D3H and HF-3c optimizations on a set of examples (a UpA dinucleotide, a DNA G-quadruplex and a B-DNA fragment). Evaluation of different penalty function strengths reveals only a minor system-dependency and reasonable restraint values range from 0.01 to 0.05 Eh rad(-2) for the backbone torsions. Restraints are crucial to perform the QM calculations on biochemically relevant conformations in implicit solvation and gas phase geometry optimizations. The reasons for using restrained instead of constrained or unconstrained optimizations are explained and an open-source external optimizer is provided.
Collapse
Affiliation(s)
- Holger Kruse
- CEITEC - Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic.
| | | |
Collapse
|
38
|
Pingoud A, Wilson GG, Wende W. Type II restriction endonucleases--a historical perspective and more. Nucleic Acids Res 2014; 42:7489-527. [PMID: 24878924 PMCID: PMC4081073 DOI: 10.1093/nar/gku447] [Citation(s) in RCA: 169] [Impact Index Per Article: 16.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2014] [Revised: 05/02/2014] [Accepted: 05/07/2014] [Indexed: 12/17/2022] Open
Abstract
This article continues the series of Surveys and Summaries on restriction endonucleases (REases) begun this year in Nucleic Acids Research. Here we discuss 'Type II' REases, the kind used for DNA analysis and cloning. We focus on their biochemistry: what they are, what they do, and how they do it. Type II REases are produced by prokaryotes to combat bacteriophages. With extreme accuracy, each recognizes a particular sequence in double-stranded DNA and cleaves at a fixed position within or nearby. The discoveries of these enzymes in the 1970s, and of the uses to which they could be put, have since impacted every corner of the life sciences. They became the enabling tools of molecular biology, genetics and biotechnology, and made analysis at the most fundamental levels routine. Hundreds of different REases have been discovered and are available commercially. Their genes have been cloned, sequenced and overexpressed. Most have been characterized to some extent, but few have been studied in depth. Here, we describe the original discoveries in this field, and the properties of the first Type II REases investigated. We discuss the mechanisms of sequence recognition and catalysis, and the varied oligomeric modes in which Type II REases act. We describe the surprising heterogeneity revealed by comparisons of their sequences and structures.
Collapse
Affiliation(s)
- Alfred Pingoud
- Institute of Biochemistry, Justus-Liebig-University Giessen, Heinrich-Buff-Ring 58, D-35392 Giessen, Germany
| | - Geoffrey G Wilson
- New England Biolabs Inc., 240 County Road, Ipswich, MA 01938-2723, USA
| | - Wolfgang Wende
- Institute of Biochemistry, Justus-Liebig-University Giessen, Heinrich-Buff-Ring 58, D-35392 Giessen, Germany
| |
Collapse
|
39
|
Zhang F, Tsunoda M, Kikuchi Y, Wilkinson O, Millington CL, Margison GP, Williams DM, Takénaka A. O6-Carboxymethylguanine in DNA forms a sequence context-dependent wobble base-pair structure with thymine. ACTA ACUST UNITED AC 2014; 70:1669-79. [DOI: 10.1107/s1399004714006178] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2014] [Accepted: 03/19/2014] [Indexed: 11/10/2022]
Abstract
N-Nitrosation of glycine and its derivatives generates potent alkylating agents that can lead to the formation ofO6-carboxymethylguanine (O6-CMG) in DNA.O6-CMG has been identified in DNA derived from human colon tissue and its occurrence has been linked to diets high in red and processed meats, implying an association with the induction of colorectal cancer. By analogy toO6-methylguanine,O6-CMG is expected to be mutagenic, inducing G-to-A mutations that may be the molecular basis of increased cancer risk. Previously, the crystal structure of the DNA dodecamer d(CGCG[O6-CMG]ATTCGCG) has been reported, in whichO6-CMG forms a Watson–Crick-type pair with thymine similar to the canonical A:T pair. In order to further investigate the versatility ofO6-CMG in base-pair formation, the structure of the DNA dodecamer d(CGC[O6-CMG]AATTTGCG) containingO6-CMG at a different position has been determined by X-ray crystallography using four crystal forms obtained under conditions containing different solvent ions (Sr2+, Ba2+, Mg2+, K+or Na+) with and without Hoechst 33258. The most striking finding is that the pairing modes ofO6-CMG with T are quite different from those previously reported. In the present dodecamer, the T bases are displaced (wobbled) into the major groove to form a hydrogen bond between the thymine N3N—H and the carboxyl group ofO6-CMG. In addition, a water molecule is bridged through two hydrogen bonds between the thymine O2atom and the 2-amino group ofO6-CMG to stabilize the pairing. These interaction modes commonly occur in the four crystal forms, regardless of the differences in crystallization conditions. The previous and the present results show thatO6-CMG can form a base pair with T in two alternative modes: the Watson–Crick type and a high-wobble type, the nature of which may depend on the DNA-sequence context.
Collapse
|
40
|
Villani G. Coupling between hydrogen atoms transfer and stacking interaction in adenine-thymine/guanine-cytosine complexes: a theoretical study. J Phys Chem B 2014; 118:5439-52. [PMID: 24813562 DOI: 10.1021/jp502792r] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Four different complexes of two base pairs, an adenine-thymine and a guanine-cytosine one, have been studied in order to understand the modifications induced by the staking interaction between the two base pairs on the hydrogen atoms transfers between the bases in either base pair. The inclusion of these two kinds of interactions allows us to clarify if some properties, as the mechanism of hydrogen transfer, is exclusively a local effect of a base pair or can be modified by a more long-range interaction between the base pairs. The results on these four complexes are compared with those of the monomeric systems, the A-T and G-C base pair, and with those of the A-T and G-C dimers. The specificity of each complex and of each hydrogen bond has been analyzed.
Collapse
Affiliation(s)
- Giovanni Villani
- Istituto di Chimica dei Composti OrganoMetallici, ICCOM-UOS Pisa Area della Ricerca del CNR, Via G. Moruzzi, 1, I-56124 Pisa, Italy
| |
Collapse
|
41
|
Mukherjee S, Kailasam S, Bansal M, Bhattacharyya D. Energy hyperspace for stacking interaction inAU/AUdinucleotide step: Dispersion-corrected density functional theory study. Biopolymers 2013; 101:107-20. [DOI: 10.1002/bip.22289] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2013] [Revised: 05/13/2013] [Accepted: 05/13/2013] [Indexed: 12/20/2022]
Affiliation(s)
- Sanchita Mukherjee
- Biophysics Division; Saha Institute of Nuclear Physics; Kolkata 700064 India
| | - Senthilkumar Kailasam
- Molecular Biophysics Unit; Indian Institute of Science; Bangalore 560012 Karnataka India
| | - Manju Bansal
- Molecular Biophysics Unit; Indian Institute of Science; Bangalore 560012 Karnataka India
| | | |
Collapse
|
42
|
Zhang F, Tsunoda M, Suzuki K, Kikuchi Y, Wilkinson O, Millington CL, Margison GP, Williams DM, Czarina Morishita E, Takénaka A. Structures of DNA duplexes containing O6-carboxymethylguanine, a lesion associated with gastrointestinal cancer, reveal a mechanism for inducing pyrimidine transition mutations. Nucleic Acids Res 2013; 41:5524-32. [PMID: 23580550 PMCID: PMC3664797 DOI: 10.1093/nar/gkt198] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
N-nitrosation of glycine and its derivatives generates potent alkylating agents that can lead to the formation of O(6)-carboxymethylguanine (O(6)-CMG) in DNA. O(6)-CMG has been identified in DNA derived from human colon tissue, and its occurrence has been linked to diets high in red and processed meats. By analogy to O(6)-methylguanine, O(6)-CMG is expected to be highly mutagenic, inducing G to A mutations during DNA replication that can increase the risk of gastrointestinal and other cancers. Two crystal structures of DNA dodecamers d(CGCG[O(6)-CMG]ATTCGCG) and d(CGC[O(6)-CMG]AATTCGCG) in complex with Hoechst33258 reveal that each can form a self-complementary duplex to retain the B-form conformation. Electron density maps clearly show that O(6)-CMG forms a Watson-Crick-type pair with thymine similar to the canonical A:T pair, and it forms a reversed wobble pair with cytosine. In situ structural modeling suggests that a DNA polymerase can accept the Watson-Crick-type pair of O(6)-CMG with thymine, but might also accept the reversed wobble pair of O(6)-CMG with cytosine. Thus, O(6)-CMG would permit the mis-incorporation of dTTP during DNA replication. Alternatively, the triphosphate that would be formed by carboxymethylation of the nucleotide triphosphate pool d[O(6)-CMG]TP might compete with dATP incorporation opposite thymine in a DNA template.
Collapse
Affiliation(s)
- Fang Zhang
- Graduate School of Science and Engineering, Iwaki-Meisei University, Iwaki 970-8551, Japan
| | | | | | | | | | | | | | | | | | | |
Collapse
|
43
|
Villani G. Theoretical investigation of the coupling between hydrogen-atom transfer and stacking interaction in adenine-thymine dimers. Chemphyschem 2013; 14:1256-63. [PMID: 23494877 DOI: 10.1002/cphc.201200971] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2012] [Revised: 02/07/2013] [Indexed: 12/13/2022]
Abstract
Three different dimers of the adenine-thymine (A-T) base pair are studied to point out the changes of important properties (structure, atomic charge, energy and so on) induced by coupling between the movement of the atoms in the hydrogen bonds and the stacking interaction. The comparison of these results with those for the A-T monomer system explains the role of the stacking interaction in the hydrogen-atom transfer in this biologically important base pair. The results support the idea that this coupling depends on the exact dimer considered and is different for the N-N and N-O hydrogen bonds. In particular, the correlation between the hydrogen transfer and the stacking interaction is more relevant for the N-N bridge than for the N-O one. Also, the two different mechanisms of two-hydrogen transfer (step by step and concerted) can be modified by the stacking interaction between the base pairs.
Collapse
Affiliation(s)
- Giovanni Villani
- Istituto di Chimica dei Composti OrganoMetallici, ICCOM-UOS Pisa, Area della Ricerca del CNR, Via G. Moruzzi 1, 56124 Pisa, Italy.
| |
Collapse
|
44
|
Villani G. Theoretical investigation of the coupling between hydrogen atoms transfer and stacking interaction in guanine–cytosine dimers. Phys Chem Chem Phys 2013; 15:19242-52. [DOI: 10.1039/c3cp52855c] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
45
|
Abstract
Molecular docking, molecular mechanics, molecular dynamics and relaxation matrix simulation protocols have been extensively used to generate the structural details of ligand-receptor complexes in order to understand the binding interactions between the two entities. Experimental methods like NMR spectroscopy and X-ray crystallography are known to provide structural information about ligand-receptor complexes. In addition, fluorescence spectroscopy, circular dichroism (CD) spectroscopy and molecular docking have also been utilized to decode the phenomenon of the ligand-DNA interactions, with good correlation between experimental and computational results. The DNA binding affinity was demonstrated by analysing fluorescence spectral data. Structural rigidity of DNA upon ligand binding was identified by CD spectroscopy. Docking is carried out using the DNA-Dock program which results in the binding affinity data along with structural information like interatomic distances and H-bonding, etc. The complete structural analyses of various drug-DNA complexes have afforded results that indicate a specific DNA binding pattern of these ligands. It also exhibited that certain structural features of ligands can make a ligand to be AT- or GC-specific. It was also demonstrated that changing specificity from AT base pairs to GC base pairs further improved the DNA topoisomerase inhibiting activity in certain ligands. Thus, a specific molecular recognition signature encrypted in the structure of ligand can be decoded and can be effectively employed in designing more potent antiviral and antitumour agents.
Collapse
|
46
|
Beššeová I, Banáš P, Kührová P, Košinová P, Otyepka M, Šponer J. Simulations of A-RNA Duplexes. The Effect of Sequence, Solute Force Field, Water Model, and Salt Concentration. J Phys Chem B 2012; 116:9899-916. [DOI: 10.1021/jp3014817] [Citation(s) in RCA: 59] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Ivana Beššeová
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
| | - Pavel Banáš
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Petra Kührová
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Pavlína Košinová
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Michal Otyepka
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- Regional Centre of Advanced
Technologies and Materials, Department of Physical Chemistry, Faculty
of Science, Palacky University, tr. 17
listopadu 12, 771 46, Olomouc, Czech Republic
| | - Jiří Šponer
- Institute
of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska
135, 612 65 Brno, Czech Republic
- CEITEC - Central European Institute of Technology, Campus Bohunice, Kamenice
5, 625 00 Brno, Czech Republic
| |
Collapse
|
47
|
Banáš P, Mládek A, Otyepka M, Zgarbová M, Jurečka P, Svozil D, Lankaš F, Šponer J. Can We Accurately Describe the Structure of Adenine Tracts in B-DNA? Reference Quantum-Chemical Computations Reveal Overstabilization of Stacking by Molecular Mechanics. J Chem Theory Comput 2012; 8:2448-60. [PMID: 26588974 DOI: 10.1021/ct3001238] [Citation(s) in RCA: 58] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Sequence-dependent local variations of helical parameters, structure, and flexibility are crucial for molecular recognition processes involving B-DNA. A-tracts, i.e., stretches of several consecutive adenines in one strand that are in phase with the DNA helical repeat, mediate significant DNA bending. During the past few decades, there have been intense efforts to understand the sequence dependence of helical parameters in DNA. Molecular dynamics (MD) simulations can provide valuable insights into the molecular mechanism behind the relationship between sequence and structure. However, although recent improvements in empirical force fields have helped to capture many sequence-dependent B-DNA properties, several problems remain, such as underestimation of the helical twist and suspected underestimation of the propeller twist in A-tracts. Here, we employ reference quantum mechanical (QM) calculations, explicit solvent MD, and bioinformatics to analyze the underestimation of propeller twisting of A-tracts in simulations. Although we did not identify a straightforward explanation, we discovered two imbalances in the empirical force fields. The first was overestimation of stacking interactions accompanied by underestimation of base-pairing energy, which we attribute to anisotropic polarizabilities that are not reflected by the isotropic force fields. This may lead to overstacking with potentially important consequences for MD simulations of nucleic acids. The second observed imbalance was steric clash between A(N1) and T(N3) nitrogens of AT base pairs in force-field descriptions, resulting in overestimation of the AT pair stretch in MD simulations. We also substantially extend the available set of benchmark estimated CCSD(T)/CBS data for B-DNA base stacking and provide a code that allows the generation of diverse base-stacking geometries suitable for QM computations with predefined intra- and interbase pair parameters.
Collapse
Affiliation(s)
- Pavel Banáš
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Arnošt Mládek
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, tr. 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - Daniel Svozil
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic.,Laboratory of Informatics and Chemistry, Faculty of Chemical Technology, Institute of Chemical Technology, Technicka 5, 166 28 Prague, Czech Republic
| | - Filip Lankaš
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo nam. 6, 166 10 Prague, Czech Republic
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Kralovopolska 135, 612 65 Brno, Czech Republic.,CEITEC - Central European Institute of Technology, Masaryk University, Campus Bohunice, 625 00 Brno, Czech Republic
| |
Collapse
|
48
|
Mikołajczyk MM, Czyżnikowska Ż, Czeleń P, Bielecka U, Zaleśny R, Toman P, Bartkowiak W. Quantum chemical study of hole transfer coupling in nucleic acid base complexes containing 7-deazaadenine. Chem Phys Lett 2012. [DOI: 10.1016/j.cplett.2012.04.018] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
|
49
|
Bullwinkle TJ, Samorodnitsky D, Rosati RC, Koudelka GB. Determinants of bacteriophage 933W repressor DNA binding specificity. PLoS One 2012; 7:e34563. [PMID: 22509323 PMCID: PMC3317979 DOI: 10.1371/journal.pone.0034563] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2011] [Accepted: 03/06/2012] [Indexed: 11/22/2022] Open
Abstract
We reported previously that 933W repressor apparently does not cooperatively bind to adjacent sites on DNA and that the relative affinities of 933W repressor for its operators differ significantly from that of any other lambdoid bacteriophage. These findings indicate that the operational details of the lysis-lysogeny switch of bacteriophage 933W are unique among lambdoid bacteriophages. Since the functioning of the lysis-lysogeny switch in 933W bacteriophage uniquely and solely depends on the order of preference of 933W repressor for its operators, we examined the details of how 933W repressor recognizes its DNA sites. To identify the specificity determinants, we first created a molecular model of the 933W repressor-DNA complex and tested the predicted protein-DNA interactions. These results of these studies provide a picture of how 933W repressor recognizes its DNA sites. We also show that, opposite of what is normally observed for lambdoid phages, 933W operator sequences have evolved in such a way that the presence of the most commonly found base sequences at particular operator positions serves to decrease, rather than increase, the affinity of the protein for the site. This finding cautions against assuming that a consensus sequence derived from sequence analysis defines the optimal, highest affinity DNA binding site for a protein.
Collapse
Affiliation(s)
- Tammy J. Bullwinkle
- Department of Biological Sciences, University at Buffalo, Buffalo, New York, United States of America
| | - Daniel Samorodnitsky
- Department of Biological Sciences, University at Buffalo, Buffalo, New York, United States of America
| | - Rayna C. Rosati
- Department of Biological Sciences, University at Buffalo, Buffalo, New York, United States of America
| | - Gerald B. Koudelka
- Department of Biological Sciences, University at Buffalo, Buffalo, New York, United States of America
- * E-mail:
| |
Collapse
|
50
|
OSTROVSKII VICTORE, KADYSHEVICH ELENAA. HYDRATE MODEL OF THE EQUILIBRIUM DNA–WATER SYSTEMS. INTERNATIONAL JOURNAL OF NANOSCIENCE 2012. [DOI: 10.1142/s0219581x02000103] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
A hydrate model for the DNA– H2O system is proposed. For the notions developed, available data on a tendency of H2O molecules to form structured hydrates containing atomic groups, molecules, or atoms housed within the structural cavities formed by H-bonded H2O molecules are used. It is shown that the large and small cavities of the hydrate structure II are in close geometric agreement with N-bases and deoxyribose and with other atomic groups of DNA molecules, respectively. On the basis of the model proposed and with the Watson–Crick base-pairing scheme, the number of pairs of N-bases per helix turn (11.25), the density (1.161 g/cm3), and the helix step (0.567 nm) for the quasi-equilibrium DNA (RNA)–water system are computed, the last value nearly coinciding with the well-known alpha-helix step in the protein secondary structure (5.44 nm). Assumptions on the phenomenology of some stages of mitosis are presented.
Collapse
Affiliation(s)
- VICTOR E. OSTROVSKII
- Karpov Institute of Physical Chemistry, ul. Vorontsovo Pole 10, Moscow, 105064, Russia
| | - ELENA A. KADYSHEVICH
- Karpov Institute of Physical Chemistry, ul. Vorontsovo Pole 10, Moscow, 105064, Russia
| |
Collapse
|