1
|
Ali Z, Kukhta T, Trant JF, Sharma P. An Atlas of the base inter-RNA stacks involved in bacterial translation. Biophys Chem 2024; 305:107144. [PMID: 38061282 DOI: 10.1016/j.bpc.2023.107144] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2023] [Accepted: 11/28/2023] [Indexed: 01/03/2024]
Abstract
Nucleobase-specific noncovalent interactions play a crucial role in translation. Herein, we provide a comprehensive analysis of the stacks between different RNA components in the crystal structures of the bacterial ribosome caught at different translation stages. Analysis of tRNA||rRNA stacks reveals distinct behaviour; both the A-and E-site tRNAs exhibit unique stacking patterns with 23S rRNA bases, while P-site tRNAs stack with 16S rRNA bases. Furthermore, E-site stacks exhibit diverse face orientations and ring topologies-rare for inter-chain RNA interactions-with higher average interaction energies than A or P-site stacks. This suggests that stacking may be essential for stabilizing tRNA progression through the E-site. Additionally, mRNA||rRNA stacks reveal other geometries, which depend on the tRNA binding site, whereas 16S rRNA||23S rRNA stacks highlight the importance of specific bases in maintaining the integrity of the translational complex by linking the two rRNAs. Furthermore, tRNA||mRNA stacks exhibit distinct geometries and energetics at the E-site, indicating their significance during tRNA translocation and elimination. Overall, both A and E-sites display a more diverse distribution of inter-RNA stacks compared to the P-site. Stacking interactions in the active ribosome are not simply accidental byproducts of biochemistry but are likely invoked to compensate and support the integrity and dynamics of translation.
Collapse
MESH Headings
- RNA, Ribosomal, 16S/chemistry
- RNA, Ribosomal, 16S/genetics
- RNA, Ribosomal, 16S/metabolism
- RNA, Ribosomal, 23S/chemistry
- RNA, Ribosomal, 23S/genetics
- RNA, Ribosomal, 23S/metabolism
- Ribosomes/chemistry
- RNA, Transfer/chemistry
- RNA, Transfer/genetics
- RNA, Transfer/metabolism
- RNA, Messenger/genetics
- RNA, Messenger/analysis
- RNA, Messenger/metabolism
- Nucleic Acid Conformation
Collapse
Affiliation(s)
- Zakir Ali
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Teagan Kukhta
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Ave., Windsor, ON N9B 3P4, Canada
| | - John F Trant
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Ave., Windsor, ON N9B 3P4, Canada; We-Spark Health Institute, 401 Sunset Ave., Windsor, ON N9B 3P4, Canada; Binary Star Research Services, LaSalle, ON N9J 3X8, Canada.
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India; Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Ave., Windsor, ON N9B 3P4, Canada.
| |
Collapse
|
2
|
Rutskoy B, Ozerov G, Bezrukov D. The Role of Bond Functions in Describing Intermolecular Electron Correlation for Van der Waals Dimers: A Study of (CH 4) 2 and Ne 2. Int J Mol Sci 2024; 25:1472. [PMID: 38338750 PMCID: PMC10855067 DOI: 10.3390/ijms25031472] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2023] [Revised: 01/19/2024] [Accepted: 01/22/2024] [Indexed: 02/12/2024] Open
Abstract
We present a study of the intermolecular interactions in van der Waals complexes of methane and neon dimers within the framework of the CCSD method. This approach was implemented and applied to calculate and examine the behavior of the contracted two-particle reduced density matrix (2-RDM). It was demonstrated that the region near the minimum of the two-particle density matrix correlation part, corresponding to the primary bulk of the Coulomb hole contribution, exerts a significant influence on the dispersion interaction energetics of the studied systems. As a result, the bond functions approach was applied to improve the convergence performance for the intermolecular correlation energy results with respect to the size of the atomic basis. For this, substantial acceleration was achieved by introducing an auxiliary basis of bond functions centered on the minima of the 2-RDM. For both methane and neon dimers, this general conclusion was confirmed with a series of CCSD calculations for the 2-RDM and the correlation energies.
Collapse
Affiliation(s)
- Bogdan Rutskoy
- National Research Centre “Kurchatov Institute”, Moscow 123182, Russia;
- Institute of Nuclear Physics and Technology, National Research Nuclear University “MEPhI” (Moscow Engineering Physics Institute), Moscow 115409, Russia
- Chemistry Department, M.V. Lomonosov Moscow State University, Moscow 119991, Russia;
| | - Georgiy Ozerov
- Chemistry Department, M.V. Lomonosov Moscow State University, Moscow 119991, Russia;
| | - Dmitry Bezrukov
- Chemistry Department, M.V. Lomonosov Moscow State University, Moscow 119991, Russia;
| |
Collapse
|
3
|
Ali Z, Kaur S, Kukhta T, Abu-Saleh AAAA, Jhunjhunwala A, Mitra A, Trant JF, Sharma P. Structural Mapping of the Base Stacks Containing Post-transcriptionally Modified Bases in RNA. J Phys Chem B 2023. [PMID: 37369074 DOI: 10.1021/acs.jpcb.3c02268] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/29/2023]
Abstract
Post-transcriptionally modified bases play vital roles in many biochemical processes involving RNA. Analysis of the non-covalent interactions associated with these bases in RNA is crucial for providing a more complete understanding of the RNA structure and function; however, the characterization of these interactions remains understudied. To address this limitation, we present a comprehensive analysis of base stacks involving all crystallographic occurrences of the most biologically relevant modified bases in a large dataset of high-resolution RNA crystal structures. This is accompanied by a geometrical classification of the stacking contacts using our established tools. Coupled with quantum chemical calculations and an analysis of the specific structural context of these stacks, this provides a map of the stacking conformations available to modified bases in RNA. Overall, our analysis is expected to facilitate structural research on altered RNA bases.
Collapse
Affiliation(s)
- Zakir Ali
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Sarabjeet Kaur
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
- Surface Chemistry and Catalysis: Characterisation and Application Team (COK-KAT), Leuven (Arenberg) Celestijnenlaan 200f─Box 2461, 3001 Leuven, Belgium
| | - Teagan Kukhta
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
| | - Abd Al-Aziz A Abu-Saleh
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
- Binary Star Research Services, LaSalle, Ontario N9J 3X8, Canada
| | - Ayush Jhunjhunwala
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology (IIIT-H), Gachibowli, Hyderabad, Telangana 500032, India
| | - Abhijit Mitra
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology (IIIT-H), Gachibowli, Hyderabad, Telangana 500032, India
| | - John F Trant
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
- Binary Star Research Services, LaSalle, Ontario N9J 3X8, Canada
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, Ontario N9B 3P4, Canada
| |
Collapse
|
4
|
Ali Z, Goyal A, Jhunjhunwala A, Mitra A, Trant JF, Sharma P. Structural and Energetic Features of Base-Base Stacking Contacts in RNA. J Chem Inf Model 2023; 63:655-669. [PMID: 36635230 DOI: 10.1021/acs.jcim.2c01116] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
Nucleobase π-π stacking is one of the crucial organizing interactions within three-dimensional (3D) RNA architectures. Characterizing the structural variability of these contacts in RNA crystal structures will help delineate their subtleties and their role in determining function. This analysis of different stacking geometries found in RNA X-ray crystal structures is the largest such survey to date; coupled with quantum-mechanical calculations on typical representatives of each possible stacking arrangement, we determined the distribution of stacking interaction energies. A total of 1,735,481 stacking contacts, spanning 359 of the 384 theoretically possible distinct stacking geometries, were identified. Our analysis reveals preferential occurrences of specific consecutive stacking arrangements in certain regions of RNA architectures. Quantum chemical calculations suggest that 88 of the 359 contacts possess intrinsically stable stacking geometries, whereas the remaining stacks require the RNA backbone or surrounding macromolecular environment to force their formation and maintain their stability. Our systematic analysis of π-π stacks in RNA highlights trends in the occurrence and localization of these noncovalent interactions and may help better understand the structural intricacies of functional RNA-based molecular architectures.
Collapse
Affiliation(s)
- Zakir Ali
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh160014, India
| | - Ambika Goyal
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh160014, India
| | - Ayush Jhunjhunwala
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad, Gachibowli, Hyderabad, Telangana500032, India
| | - Abhijit Mitra
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad, Gachibowli, Hyderabad, Telangana500032, India
| | - John F Trant
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, OntarioN9B 3P4, Canada
- Binary Star Research Services, LaSalle, OntarioN9J 3X8, Canada
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh160014, India
- Department of Chemistry and Biochemistry, University of Windsor, 401 Sunset Avenue, Windsor, OntarioN9B 3P4, Canada
| |
Collapse
|
5
|
De Nicola A, Correa A, Bracco S, Perego J, Sozzani P, Comotti A, Milano G. Collective dynamics of molecular rotors in periodic mesoporous organosilica: a combined solid-state 2H-NMR and molecular dynamics simulation study. Phys Chem Chem Phys 2022; 24:666-673. [PMID: 34904981 DOI: 10.1039/d1cp05013c] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Molecular rotors offer a platform to realize controlled dynamics and modulate the functions of solids. The motional mechanisms in arrays of rotors have not been explored in depth. Crystal-like porous organosilicas, comprising p-phenylene rotators pivoted onto a siloxane scaffold, were modelled using molecular dynamics (MD) simulations. Long simulations, on a microsecond scale, allowed to follow the reorientation statistics of rotor collections and single out group configurations and frequency distributions as a function of temperature. The motions observed in the MD simulations support a multiple-site model for rotor reorientations. Computed motional frequencies revealed a complex rotatory phenomenon combining an ultra-fast libration motion (oscillation up to 30°) with a slow and fast 180° flip reorientation. Adopting a multiple-site model provides a more accurate simulation of the 2H-NMR spectra and a rationalization of their temperature dependence. In particular, rotators endowed with distinct rates could be explained by the presence of slower rings locked in a T-shaped conformation.
Collapse
Affiliation(s)
- Antonio De Nicola
- Scuola Superiore Meridionale, Largo San Marcellino 10, 80138 Napoli, Italy.
| | - Andrea Correa
- Department of Chemical Sciences, University of Naples Federico II, Via Cintia, Complesso Monte S. Angelo, 80126, Napoli, Italy
| | - Silvia Bracco
- Department of Materials Science, University of Milano - Bicocca, Via R. Cozzi 55, 20125, Milano, Italy.
| | - Jacopo Perego
- Department of Materials Science, University of Milano - Bicocca, Via R. Cozzi 55, 20125, Milano, Italy.
| | - Piero Sozzani
- Department of Materials Science, University of Milano - Bicocca, Via R. Cozzi 55, 20125, Milano, Italy.
| | - Angiolina Comotti
- Department of Materials Science, University of Milano - Bicocca, Via R. Cozzi 55, 20125, Milano, Italy.
| | - Giuseppe Milano
- Department of Chemical, Materials and Production Engineering, University of Naples Federico II, Piazzale Tecchio 80, 80125 Naploli, Italy
| |
Collapse
|
6
|
Kaur S, Grover P, Wetmore SD, Sharma P. Role of Stacking Interactions in the Stability of Primitive Genetics: A Quantum Chemical View. J Chem Inf Model 2021; 61:4321-4330. [PMID: 34491053 DOI: 10.1021/acs.jcim.1c00276] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The origin of genetic material on earth is an age-old, entangled mystery that lacks a unanimous explanation. Recent studies have suggested that noncanonical bases such as barbituric acid (BA), melamine (MM), cyanuric acid (CA), and 2,4,6-triaminopyrimidine (TAP) may have undergone molecular selection within the "prebiotic soup" to spontaneously form supramolecular assemblies, which then covalently assembled into an RNA-like polymer (preRNA). However, information on the role of intrinsic interactions of these candidate heterocycles in their molecular selection as the components of preRNA, and the subsequent transition from preRNA to RNA, is currently missing in the literature. To fill this gap in our knowledge on the origin and evolution of primitive genetics, the present work employs density functional theory (B3LYP-D3) to evaluate and compare the stacking propensities of dimers containing prebiotic noncanonical (BA, MM, CA, and TAP) and/or canonical RNA bases (A, C, G, and U). Our detailed analysis of the variation in stacking strength with respect to four characteristic geometrical parameters between the monomers [i.e., the vertical distance, the angle of rotation, and (two) displacements in the x and y directions] reveals that stacking between nonidentical bases is preferred over identical bases for both prebiotic-prebiotic and canonical-canonical dimers. This not only underscores the similarity between the fundamental chemical properties of preRNA and RNA constituents but also supports the likelihood of the evolution of modern (RNA) genetics from primitive (preRNA) genetics. Furthermore, greater average stacking stabilization of canonical dimers than that of dimers containing one canonical and one preRNA nucleobase (by ∼5 kJ mol-1) or dimers solely containing preRNA nucleobases (by ∼12 kJ mol-1) indicates that enhanced stacking is an important factor that may have spurred the evolution of preRNA to an intermediate informational polymer to RNA. More importantly, our study identifies the central roles of CA, BA, and TAP in stacking stabilization within the preRNA and of BA in stacking interactions within the intermediate polymers and suggests that these heterocycles may have played distinct roles in various stages during the evolution from preRNA to RNA. Overall, our results highlight the significance of stacking interactions in the selection of nucleobase components of preRNA.
Collapse
Affiliation(s)
- Sarabjeet Kaur
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| | - Payal Grover
- Department of Chemistry, Dayanand Anglo-Vedic (DAV) College, Sector 10, Chandigarh 16011, India
| | - Stacey D Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, Alberta T1K 3M4, Canada
| | - Purshotam Sharma
- Computational Biochemistry Laboratory, Department of Chemistry and Centre for Advanced Studies in Chemistry, Panjab University, Chandigarh 160014, India
| |
Collapse
|
7
|
Deb I, Wong H, Tacubao C, Frank AT. Quantum Mechanics Helps Uncover Atypical Recognition Features in the Flavin Mononucleotide Riboswitch. J Phys Chem B 2021; 125:8342-8350. [PMID: 34310879 DOI: 10.1021/acs.jpcb.1c02702] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Estimating the binding energies of small molecules to RNA could help uncover their molecular recognition characteristics and be used to rationally design RNA-targeting chemical probes. Here, we leveraged the ability of the fragment molecular orbital (FMO) method to provide detailed pairwise energetic information to examine the interactions between the aptamer domain of the flavin mononucleotide (FMN)-responsive riboswitch and small-molecule ligands. After developing an efficient protocol for executing high-level FMO calculations on RNA-ligand complexes, we applied our protocol to nine FMN-aptamer-ligand complexes. We then used the results to identify "hot-spots" within the aptamer and decomposed pairwise interactions between the hot-spot residues and the ligands. Interestingly, we found that several of these hot-spot residues interact with the ligands via atypical CH···O hydrogen bonds and anion-π contacts, as well as (face-to-edge) T-shaped π-π interactions. We envision that our results should pave the way for the wider and more prominent use of FMO calculations to study structure-energy relationships in diverse RNA-ligand systems, which in turn may provide a basis for dissecting the molecular recognition characteristics of RNAs.
Collapse
Affiliation(s)
- Indrajit Deb
- Biophysics Program, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Hazel Wong
- Biophysics Program, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Colleen Tacubao
- Biophysics Program, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Aaron T Frank
- Biophysics Program, University of Michigan, Ann Arbor, Michigan 48109, United States.,Chemistry Department, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
8
|
Lin M, Malik FK, Guo JT. A comparative study of protein-ssDNA interactions. NAR Genom Bioinform 2021; 3:lqab006. [PMID: 33655206 PMCID: PMC7902235 DOI: 10.1093/nargab/lqab006] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2020] [Revised: 11/24/2020] [Accepted: 01/26/2021] [Indexed: 12/18/2022] Open
Abstract
Single-stranded DNA-binding proteins (SSBs) play crucial roles in DNA replication, recombination and repair, and serve as key players in the maintenance of genomic stability. While a number of SSBs bind single-stranded DNA (ssDNA) non-specifically, the others recognize and bind specific ssDNA sequences. The mechanisms underlying this binding discrepancy, however, are largely unknown. Here, we present a comparative study of protein-ssDNA interactions by annotating specific and non-specific SSBs and comparing structural features such as DNA-binding propensities and secondary structure types of residues in SSB-ssDNA interactions, protein-ssDNA hydrogen bonding and π-π interactions between specific and non-specific SSBs. Our results suggest that protein side chain-DNA base hydrogen bonds are the major contributors to protein-ssDNA binding specificity, while π-π interactions may mainly contribute to binding affinity. We also found the enrichment of aspartate in the specific SSBs, a key feature in specific protein-double-stranded DNA (dsDNA) interactions as reported in our previous study. In addition, no significant differences between specific and non-specific groups with respect of conformational changes upon ssDNA binding were found, suggesting that the flexibility of SSBs plays a lesser role than that of dsDNA-binding proteins in conferring binding specificity.
Collapse
Affiliation(s)
- Maoxuan Lin
- Department of Bioinformatics and Genomics, The University of North Carolina at Charlotte, Charlotte, NC 28223, USA
| | - Fareeha K Malik
- Department of Bioinformatics and Genomics, The University of North Carolina at Charlotte, Charlotte, NC 28223, USA
- Research Center of Modeling and Simulation, National University of Science and Technology, Islamabad, 44000, Pakistan
| | - Jun-tao Guo
- Department of Bioinformatics and Genomics, The University of North Carolina at Charlotte, Charlotte, NC 28223, USA
| |
Collapse
|
9
|
Mondal P. In silico decryption of serotonin-receptor binding: local non-covalent interactions and long-range conformational changes. RSC Adv 2020; 10:37995-38003. [PMID: 35515142 PMCID: PMC9057181 DOI: 10.1039/d0ra05559j] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2020] [Revised: 03/08/2021] [Accepted: 09/15/2020] [Indexed: 01/10/2023] Open
Abstract
Serotonin-receptor binding is the key step in the process behind serotonin functionality, including several psychological and physiological behaviours. This study is focused on identifying the main non-covalent interactions controlling the stability of serotonin-receptor complexes as well as the main conformational changes in the receptor due to serotonin-receptor binding using classical molecular dynamics simulations and quantum chemical calculations. A qualitative analysis based on two order parameters ((i) the centre of mass distance and (ii) the angle between the surface normals of each aromatic residue and serotonin in the binding site) on the serotonin-receptor complex trajectory suggests that the T-type stacking interaction is predominant in the binding site. Quantum chemical calculations of the stacking interaction energy provide the quantitative contributions of important aromatic residues to the stabilization of the complex. Furthermore, a three body stacking interaction (named 'L'-type) was observed and likely contributes to the stability of the complex. Direct and water-mediated hydrogen bonding between the residues in the binding site and serotonin contributes to the complex stability. Principal component analysis of the molecular dynamics simulation trajectory of the serotonin-receptor complex and the apo-receptor in water indicates that the whole receptor is significantly stabilized due to serotonin binding. An analysis based on the dynamic cross correlation function reflects the strong correlation between trans-membrane (TM)3, TM5, TM6 (containing residues responsible for the stacking interaction and hydrogen bonding) and mini-G0 which may participate in signal transduction leading to the functionality of serotonin.
Collapse
Affiliation(s)
- Padmabati Mondal
- Department of Chemistry and Center for Atomic, Molecular and Optical Sciences and Technologies, Indian Institute of Science Education and Research (IISER) Tirupati Tirupati 517507 Andhra Pradesh India +91 877 2500 926
| |
Collapse
|
10
|
Kaur R, Nikkel DJ, Wetmore SD. Computational studies of DNA repair: Insights into the function of monofunctional DNA glycosylases in the base excision repair pathway. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1471] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Rajwinder Kaur
- Department of Chemistry and Biochemistry University of Lethbridge Lethbridge Alberta Canada
| | - Dylan J. Nikkel
- Department of Chemistry and Biochemistry University of Lethbridge Lethbridge Alberta Canada
| | - Stacey D. Wetmore
- Department of Chemistry and Biochemistry University of Lethbridge Lethbridge Alberta Canada
| |
Collapse
|
11
|
Retureau R, Oguey C, Mauffret O, Hartmann B. Structural Explorations of NCp7-Nucleic Acid Complexes Give Keys to Decipher the Binding Process. J Mol Biol 2019; 431:1966-1980. [PMID: 30876916 DOI: 10.1016/j.jmb.2019.03.002] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2019] [Revised: 02/25/2019] [Accepted: 03/02/2019] [Indexed: 02/06/2023]
Abstract
A comprehensive view of all the structural aspects related to NCp7 is essential to understand how this protein, crucial in many steps of the HIV-1 cycle, binds and anneals nucleic acids (NAs), mainly thanks to two zinc fingers, ZF1 and ZF2. Here, we inspected the structural properties of the available experimental models of NCp7 bound to either DNA or RNA molecules, or free of ligand. Our analyses included the characterization of the relative positioning of ZF1 and ZF2, accessibility measurements and the exhaustive, quantitative mapping of the contacts between amino acids and nucleotides by a recent tessellation method, VLDM. This approach unveiled the intimate connection between NA binding process and the conformations explored by the free protein. It also provided new insights into the functional specializations of ZF1 and ZF2. The larger accessibility of ZF2 in free NCp7 and the consistency of the ZF2/NA interface in different models and conditions give ZF2 the lead of the binding process. ZF1 contributes to stabilize the complexes through various organizations of the ZF1/NA interface. This work outcome is a global binding scheme of NCp7 to DNA and RNA, and an example of how protein-NA complexes are stabilized.
Collapse
Affiliation(s)
- Romain Retureau
- LBPA, UMR 8113, ENS Paris-Saclay-CNRS, 61 avenue du Président Wilson, 94235 Cachan cedex, France
| | - Christophe Oguey
- LPTM, CNRS UMR 8089, Université de Cergy-Pontoise, 2 avenue Adolphe Chauvin, 95031 Cergy-Pontoise, France
| | - Olivier Mauffret
- LBPA, UMR 8113, ENS Paris-Saclay-CNRS, 61 avenue du Président Wilson, 94235 Cachan cedex, France.
| | - Brigitte Hartmann
- LBPA, UMR 8113, ENS Paris-Saclay-CNRS, 61 avenue du Président Wilson, 94235 Cachan cedex, France.
| |
Collapse
|
12
|
Punihaole D, Workman RJ, Upadhyay S, Bruggen CV, Schmitz AJ, Reineke TM, Frontiera RR. New Insights into Quinine-DNA Binding Using Raman Spectroscopy and Molecular Dynamics Simulations. J Phys Chem B 2018; 122:9840-9851. [PMID: 30336027 PMCID: PMC6425490 DOI: 10.1021/acs.jpcb.8b05795] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Quinine's ability to bind DNA and potentially inhibit transcription and translation has been examined as a mode of action for its antimalarial activity. UV absorption and fluorescence-based studies have lacked the chemical specificity to develop an unambiguous molecular-level picture of the binding interaction. To address this, we use Raman spectroscopy and molecular dynamics (MD) to investigate quinine-DNA interactions. We demonstrate that quinine's strongest Raman band in the fingerprint region, which derives from a symmetric stretching mode of the quinoline ring, is highly sensitive to the local chemical environment and pH. The frequency shifts observed for this mode in solvents of varying polarity can be explained in terms of the Stark effect using a simple Onsager solvation model, indicating that the vibration reports on the local electrostatic environment. However, specific chemical interactions between the quinoline ring and its environment, such as hydrogen bonding and π-stacking, perturb the frequency of this mode in a more complicated but predictable manner. We use this vibration as a spectroscopic probe to investigate the binding interaction between quinine and DNA. We find that, when the quinoline ring is protonated, quinine weakly intercalates into DNA by forming π-stacking interactions with the base pairs. The Raman spectra indicate that quinine can intercalate into DNA with a ratio reaching up to roughly one molecule per 25 base pairs. Our results are confirmed by MD simulations, which also show that the quinoline ring adopts a t-shaped π-stacking geometry with the DNA base pairs, whereas the quinuclidine head group weakly interacts with the phosphate backbone in the minor groove. We expect that the spectral correlations determined here will enable future studies to probe quinine's antimalarial activities, such as disrupting hemozoin biocrystallization, which is hypothesized to be, among other things, one of its primary modes of action against Plasmodium parasites.
Collapse
Affiliation(s)
- David Punihaole
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Riley J. Workman
- Department of Chemistry and Biochemistry, Center for Computational Sciences, Duquesne University, Pittsburgh, Pennsylvania 15282, United States
| | - Shiv Upadhyay
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Craig Van Bruggen
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Andrew J. Schmitz
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Theresa M. Reineke
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Renee R. Frontiera
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
13
|
Jacobs M, Greff Da Silveira L, Prampolini G, Livotto PR, Cacelli I. Interaction Energy Landscapes of Aromatic Heterocycles through a Reliable yet Affordable Computational Approach. J Chem Theory Comput 2018; 14:543-556. [DOI: 10.1021/acs.jctc.7b00602] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Affiliation(s)
- Matheus Jacobs
- Instituto
de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Leandro Greff Da Silveira
- Instituto
de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
- Departamento
de Ciências Exatas e da Terra, Universidade Regional Integrada do Alto Uruguay e da Missões (URI), Avenida Assis Brasil 709, CEP 98400-00 Frederico Westphalen, Brazil
| | - Giacomo Prampolini
- Istituto di Chimica
dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Paolo Roberto Livotto
- Instituto
de Química, Universidade Federal do Rio Grande do Sul, Avenida Bento Gonçalves 9500, CEP 91501-970 Porto Alegre, Brazil
| | - Ivo Cacelli
- Istituto di Chimica
dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi
3, I-56124 Pisa, Italy
| |
Collapse
|
14
|
Chawla M, Chermak E, Zhang Q, Bujnicki JM, Oliva R, Cavallo L. Occurrence and stability of lone pair-π stacking interactions between ribose and nucleobases in functional RNAs. Nucleic Acids Res 2017; 45:11019-11032. [PMID: 28977572 PMCID: PMC5737201 DOI: 10.1093/nar/gkx757] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2017] [Accepted: 08/17/2017] [Indexed: 12/13/2022] Open
Abstract
The specific folding pattern and function of RNA molecules lies in various weak interactions, in addition to the strong base-base pairing and stacking. One of these relatively weak interactions, characterized by the stacking of the O4' atom of a ribose on top of the heterocycle ring of a nucleobase, has been known to occur but has largely been ignored in the description of RNA structures. We identified 2015 ribose-base stacking interactions in a high-resolution set of non-redundant RNA crystal structures. They are widespread in structured RNA molecules and are located in structural motifs other than regular stems. Over 50% of them involve an adenine, as we found ribose-adenine contacts to be recurring elements in A-minor motifs. Fewer than 50% of the interactions involve a ribose and a base of neighboring residues, while approximately 30% of them involve a ribose and a nucleobase at least four residues apart. Some of them establish inter-domain or inter-molecular contacts and often implicate functionally relevant nucleotides. In vacuo ribose-nucleobase stacking interaction energies were calculated by quantum mechanics methods. Finally, we found that lone pair-π stacking interactions also occur between ribose and aromatic amino acids in RNA-protein complexes.
Collapse
Affiliation(s)
- Mohit Chawla
- King Abdullah University of Science and Technology (KAUST), Physical Sciences and Engineering Division, Thuwal 23955-6900, Saudi Arabia.,Kaust Catalysis Center, Thuwal 23955-6900, Saudi Arabia
| | - Edrisse Chermak
- King Abdullah University of Science and Technology (KAUST), Physical Sciences and Engineering Division, Thuwal 23955-6900, Saudi Arabia.,Kaust Catalysis Center, Thuwal 23955-6900, Saudi Arabia
| | - Qingyun Zhang
- King Abdullah University of Science and Technology (KAUST), Physical Sciences and Engineering Division, Thuwal 23955-6900, Saudi Arabia
| | - Janusz M Bujnicki
- Laboratory of Bioinformatics and Protein Engineering, International Institute of Molecular and Cell Biology in Warsaw, ul. Ks. Trojdena 4, 02-109 Warsaw, Poland.,Laboratory of Bioinformatics, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Adam Mickiewicz University, Umultowska 89, 61-614 Poznan, Poland
| | - Romina Oliva
- Department of Sciences and Technologies, University Parthenope of Naples, Centro Direzionale Isola C4, I-80143 Naples, Italy.,King Abdullah University of Science and Technology (KAUST), Computational Bioscience Research Center (CBRC), Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division, Thuwal 23955-6900, Saudi Arabia
| | - Luigi Cavallo
- King Abdullah University of Science and Technology (KAUST), Physical Sciences and Engineering Division, Thuwal 23955-6900, Saudi Arabia.,Kaust Catalysis Center, Thuwal 23955-6900, Saudi Arabia
| |
Collapse
|
15
|
Park J, Kim H, Kim S, Lee D, Kim MS, Shin DH. Crystal structure of D-glycero-Β-D-manno-heptose-1-phosphate adenylyltransferase fromBurkholderia pseudomallei. Proteins 2017; 86:124-131. [DOI: 10.1002/prot.25398] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2017] [Revised: 10/02/2017] [Accepted: 10/04/2017] [Indexed: 12/17/2022]
Affiliation(s)
- Jimin Park
- College of Pharmacy; Ewha W. University; Seoul Republic of Korea
| | - Hyojin Kim
- College of Pharmacy; Ewha W. University; Seoul Republic of Korea
| | - Suwon Kim
- College of Pharmacy; Ewha W. University; Seoul Republic of Korea
| | - Daeun Lee
- College of Pharmacy; Ewha W. University; Seoul Republic of Korea
| | - Mi-Sun Kim
- College of Pharmacy; Ewha W. University; Seoul Republic of Korea
| | - Dong Hae Shin
- College of Pharmacy; Ewha W. University; Seoul Republic of Korea
| |
Collapse
|
16
|
Gao XC, Hao Q, Wang CS. Improved Polarizable Dipole–Dipole Interaction Model for Hydrogen Bonding, Stacking, T-Shaped, and X–H···π Interactions. J Chem Theory Comput 2017; 13:2730-2741. [DOI: 10.1021/acs.jctc.6b00936] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Affiliation(s)
- Xi-Chan Gao
- School of Chemistry and Chemical
Engineering, Liaoning Normal University, Dalian 116029, P. R. China
| | - Qiang Hao
- School of Chemistry and Chemical
Engineering, Liaoning Normal University, Dalian 116029, P. R. China
| | - Chang-Sheng Wang
- School of Chemistry and Chemical
Engineering, Liaoning Normal University, Dalian 116029, P. R. China
| |
Collapse
|
17
|
An electrochemical aptasensor for multiplex antibiotics detection based on metal ions doped nanoscale MOFs as signal tracers and RecJ f exonuclease-assisted targets recycling amplification. Talanta 2016; 161:867-874. [PMID: 27769495 DOI: 10.1016/j.talanta.2016.09.051] [Citation(s) in RCA: 59] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2016] [Revised: 09/17/2016] [Accepted: 09/21/2016] [Indexed: 12/13/2022]
Abstract
An ultrasensitive electrochemical aptasensor for simultaneous detection of oxytetracycline (OTC) and kanamycin (KAN) has been developed based on metal ions doped metal organic frame materials (MOFs) as signal tracers and RecJf exonuclease-catalyzed targets recycling amplification. The aptasensor consists of capture beads (the anti-single-stranded DNA Antibody, as anti-ssDNA Ab, labeled on Dynabeads) and nanoscale MOF (NMOF) based signal tracers (simplified as Apts-MNM, the NMOF labeled with metal ions and the aptamers). Particularly, the MOF (UiO-66-NH2), with large internal surface areas, ultrahigh porosity and abundant amine groups in the pores, was employed as substrates to carry plenty of metal ions (Pb2+ or Cd2+) and label aptamers of OTC or KAN. Thus, the aptasensor is formed by the specific recognition between anti-ssDNA Ab and aptamers. In the presence of targets (OTC and KAN), aptamers prefer to form targets-Apts-MNM complexes in lieu of anti-ssDNA Ab-aptamer complexes, which results in the dissociation of Apts-MNM from capture beads. With the employment of RecJf exonuclease, targets-Apts-MNM in supernatant was digested into mononucleotides and liberated the target, which can further participate in the next reaction cycling to produce more signal tracers. After magnetic separation, the enhanced square wave voltammetry (SWV) signals were produced from signal tracers. The aptasensor exhibited a linear correlation in the range from 0.5pM to 50nM, with detection limits of 0.18pM and 0.15pM (S/N=3) toward OTC and KAN respectively. This strategy provides specificity and sensitive approach for multiplex antibiotics detection and has promising applications in food analysis.
Collapse
|
18
|
Wilson KA, Holland DJ, Wetmore SD. Topology of RNA-protein nucleobase-amino acid π-π interactions and comparison to analogous DNA-protein π-π contacts. RNA (NEW YORK, N.Y.) 2016; 22:696-708. [PMID: 26979279 PMCID: PMC4836644 DOI: 10.1261/rna.054924.115] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/18/2015] [Accepted: 02/13/2016] [Indexed: 06/05/2023]
Abstract
The present work analyzed 120 high-resolution X-ray crystal structures and identified 335 RNA-protein π-interactions (154 nonredundant) between a nucleobase and aromatic (W, H, F, or Y) or acyclic (R, E, or D) π-containing amino acid. Each contact was critically analyzed (including using a visual inspection protocol) to determine the most prevalent composition, structure, and strength of π-interactions at RNA-protein interfaces. These contacts most commonly involve F and U, with U:F interactions comprising one-fifth of the total number of contacts found. Furthermore, the RNA and protein π-systems adopt many different relative orientations, although there is a preference for more parallel (stacked) arrangements. Due to the variation in structure, the strength of the intermolecular forces between the RNA and protein components (as determined from accurate quantum chemical calculations) exhibits a significant range, with most of the contacts providing significant stability to the associated RNA-protein complex (up to -65 kJ mol(-1)). Comparison to the analogous DNA-protein π-interactions emphasizes differences in RNA- and DNA-protein π-interactions at the molecular level, including the greater abundance of RNA contacts and the involvement of different nucleobase/amino acid residues. Overall, our results provide a clearer picture of the molecular basis of nucleic acid-protein binding and underscore the important role of these contacts in biology, including the significant contribution of π-π interactions to the stability of nucleic acid-protein complexes. Nevertheless, more work is still needed in this area in order to further appreciate the properties and roles of RNA nucleobase-amino acid π-interactions in nature.
Collapse
Affiliation(s)
- Katie A Wilson
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, Alberta T1K 3M4, Canada
| | - Devany J Holland
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, Alberta T1K 3M4, Canada
| | - Stacey D Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, Lethbridge, Alberta T1K 3M4, Canada
| |
Collapse
|
19
|
Wang W, Liu J, Sun L. Surface shapes and surrounding environment analysis of single- and double-stranded DNA-binding proteins in protein-DNA interface. Proteins 2016; 84:979-89. [PMID: 27038080 DOI: 10.1002/prot.25045] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2016] [Revised: 03/15/2016] [Accepted: 03/25/2016] [Indexed: 11/12/2022]
Abstract
Protein-DNA bindings are critical to many biological processes. However, the structural mechanisms underlying these interactions are not fully understood. Here, we analyzed the residues shape (peak, flat, or valley) and the surrounding environment of double-stranded DNA-binding proteins (DSBs) and single-stranded DNA-binding proteins (SSBs) in protein-DNA interfaces. In the results, we found that the interface shapes, hydrogen bonds, and the surrounding environment present significant differences between the two kinds of proteins. Built on the investigation results, we constructed a random forest (RF) classifier to distinguish DSBs and SSBs with satisfying performance. In conclusion, we present a novel methodology to characterize protein interfaces, which will deepen our understanding of the specificity of proteins binding to ssDNA (single-stranded DNA) or dsDNA (double-stranded DNA). Proteins 2016; 84:979-989. © 2016 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Wei Wang
- Department of Computer Science and Technology, College of Computer and Information Engineering, Henan Normal University, Xinxiang, 453007, China.,Laboratory of Computation Intelligence and Information Processing, Engineering Technology Research Center for Computing Intelligence and Data Mining, Henan Province, China
| | - Juan Liu
- Institute of Computer Software, School of Computer, Wuhan University, Wuhan, 430072, China
| | - Lin Sun
- Department of Computer Science and Technology, College of Computer and Information Engineering, Henan Normal University, Xinxiang, 453007, China.,Laboratory of Computation Intelligence and Information Processing, Engineering Technology Research Center for Computing Intelligence and Data Mining, Henan Province, China
| |
Collapse
|
20
|
Synthesis of Diphenyl Pyridazinone-based flexible system for conformational studies through weak noncovalent interactions: Application in DNA binding. J CHEM SCI 2016. [DOI: 10.1007/s12039-016-1059-3] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
|
21
|
Chehel Amirani M, Tang T. A QM:MM model for the interaction of DNA nucleotides with carbon nanotubes. Phys Chem Chem Phys 2015; 17:7564-75. [PMID: 25708519 DOI: 10.1039/c4cp05222f] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Hybrid materials formed by DNA and carbon nanotubes (CNTs) have shown very interesting properties, but their simulation in solution using quantum mechanical approaches is still a challenge in the computational chemistry community. In this paper, we developed a QM:MM model to study the interactions between charged DNA nucleotides and carbon nanotubes in solution. All four types of DNA nucleotides were taken to interact with two CNTs of similar diameter but different chiralities: (4,4) and (7,0). The nucleotides and CNTs were treated at the QM level, while added water and neutralizing ions were modeled at the MM level. ONIOM simulations were performed at the (M06-2X/6-31G(d):Amber) level for the hybrids, as well as for individually solvated CNTs and nucleotides, which allowed us to evaluate the energy of binding. Our binding energy (BE) values range from 146.60 to 503.43 kJ mol(-1), indicating strong physisorption of nucleotides on CNTs. The relatively large BE, compared with past studies on nucleobase-CNT binding in a vacuum, could be due to the larger size of nucleotides compared with nucleobases, the charges on the nucleotides, and the inclusion of solution which causes the release of water molecules upon hybridization.
Collapse
|
22
|
An Y, Bloom JWG, Wheeler SE. Quantifying the π-Stacking Interactions in Nitroarene Binding Sites of Proteins. J Phys Chem B 2015; 119:14441-50. [DOI: 10.1021/acs.jpcb.5b08126] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Yi An
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
| | - Jacob W. G. Bloom
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
| | - Steven E. Wheeler
- Department of Chemistry, Texas A&M University, College Station, Texas 77842, United States
| |
Collapse
|
23
|
Wilson KA, Wells RA, Abendong MN, Anderson CB, Kung RW, Wetmore SD. Landscape of π-π and sugar-π contacts in DNA-protein interactions. J Biomol Struct Dyn 2015; 34:184-200. [PMID: 25723403 DOI: 10.1080/07391102.2015.1013157] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
There were 1765 contacts identified between DNA nucleobases or deoxyribose and cyclic (W, H, F, Y) or acyclic (R, E, D) amino acids in 672 X-ray structures of DNA-protein complexes. In this first study to compare π-interactions between the cyclic and acyclic amino acids, visual inspection was used to categorize amino acid interactions as nucleobase π-π (according to biological edge) or deoxyribose sugar-π (according to sugar edge). Overall, 54% of contacts are nucleobase π-π interactions, which involve all amino acids, but are more common for Y, F, and R, and involve all DNA nucleobases with similar frequencies. Among binding arrangements, cyclic amino acids prefer more planar (stacked) π-systems than the acyclic counterparts. Although sugar-π interactions were only previously identified with the cyclic amino acids and were found to be less common (38%) than nucleobase-cyclic amino acid contacts, sugar-π interactions are more common than nucleobase π-π contacts for the acyclic series (61% of contacts). Similar to DNA-protein π-π interactions, sugar-π contacts most frequently involve Y and R, although all amino acids adopt many binding orientations relative to deoxyribose. These DNA-protein π-interactions stabilize biological systems, by up to approximately -40 kJ mol(-1) for neutral nucleobase or sugar-amino acid interactions, but up to approximately -95 kJ mol(-1) for positively or negatively charged contacts. The high frequency and strength, despite variation in structure and composition, of these π-interactions point to an important function in biological systems.
Collapse
Affiliation(s)
- Katie A Wilson
- a Department of Chemistry and Biochemistry , University of Lethbridge , 4401 University Drive West, Lethbridge , AB T1K 3M4 , Canada
| | - Rachael A Wells
- a Department of Chemistry and Biochemistry , University of Lethbridge , 4401 University Drive West, Lethbridge , AB T1K 3M4 , Canada
| | - Minette N Abendong
- a Department of Chemistry and Biochemistry , University of Lethbridge , 4401 University Drive West, Lethbridge , AB T1K 3M4 , Canada
| | - Colin B Anderson
- a Department of Chemistry and Biochemistry , University of Lethbridge , 4401 University Drive West, Lethbridge , AB T1K 3M4 , Canada
| | - Ryan W Kung
- a Department of Chemistry and Biochemistry , University of Lethbridge , 4401 University Drive West, Lethbridge , AB T1K 3M4 , Canada
| | - Stacey D Wetmore
- a Department of Chemistry and Biochemistry , University of Lethbridge , 4401 University Drive West, Lethbridge , AB T1K 3M4 , Canada
| |
Collapse
|
24
|
Wilson KA, Wetmore SD. A Survey of DNA–Protein π–Interactions: A Comparison of Natural Occurrences and Structures, and Computationally Predicted Structures and Strengths. CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2015. [DOI: 10.1007/978-3-319-14163-3_17] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
|
25
|
Hussain HB, Wilson KA, Wetmore SD. Serine and Cysteine π-Interactions in Nature: A Comparison of the Frequency, Structure, and Stability of Contacts Involving Oxygen and Sulfur. Aust J Chem 2015. [DOI: 10.1071/ch14598] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Despite many DNA–protein π-interactions in high-resolution crystal structures, only four X–H···π or X···π interactions were found between serine (Ser) or cysteine (Cys) and DNA nucleobase π-systems in over 100 DNA–protein complexes (where X = O for Ser and X = S for Cys). Nevertheless, 126 non-covalent contacts occur between Ser or Cys and the aromatic amino acids in many binding arrangements within proteins. Furthermore, Ser and Cys protein–protein π-interactions occur with similar frequencies and strengths. Most importantly, due to the great stability that can be provided to biological macromolecules (up to –20 kJ mol–1 for neutral π-systems or –40 kJ mol–1 for cationic π-systems), Ser and Cys π-interactions should be considered when analyzing protein stability and function.
Collapse
|
26
|
Barone V, Cacelli I, Ferretti A, Prampolini G, Villani G. Proton and Electron Transfer Mechanisms in the Formation of Neutral and Charged Quinhydrone-Like Complexes: A Multilayered Computational Study. J Chem Theory Comput 2014; 10:4883-95. [DOI: 10.1021/ct500778u] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Affiliation(s)
- Vincenzo Barone
- Scuola Normale Superiore, Piazza
dei Cavalieri 7, I-56126 Pisa, Italy
| | - Ivo Cacelli
- Dipartimento
di Chimica e Chimica Industriale, Università di Pisa, Via Risorgimento
35, I-56126 Pisa, Italy
| | - Alessandro Ferretti
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Giacomo Prampolini
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Giovanni Villani
- Istituto di Chimica dei Composti OrganoMetallici (ICCOM-CNR), Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| |
Collapse
|
27
|
Chehel Amirani M, Tang T. Binding of nucleobases with graphene and carbon nanotube: a review of computational studies. J Biomol Struct Dyn 2014; 33:1567-97. [DOI: 10.1080/07391102.2014.954315] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
28
|
Wilson KA, Wetmore SD. Complex Conformational Heterogeneity of the Highly Flexible O6-Benzyl-guanine DNA Adduct. Chem Res Toxicol 2014; 27:1310-25. [PMID: 24941023 DOI: 10.1021/tx500178x] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Katie A. Wilson
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta, Canada T1K 3M4
| | - Stacey D. Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, Alberta, Canada T1K 3M4
| |
Collapse
|
29
|
Wilson KA, Kellie JL, Wetmore SD. DNA-protein π-interactions in nature: abundance, structure, composition and strength of contacts between aromatic amino acids and DNA nucleobases or deoxyribose sugar. Nucleic Acids Res 2014; 42:6726-41. [PMID: 24744240 PMCID: PMC4041443 DOI: 10.1093/nar/gku269] [Citation(s) in RCA: 108] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
Four hundred twenty-eight high-resolution DNA-protein complexes were chosen for a bioinformatics study. Although 164 crystal structures (38% of those searched) contained no interactions, 574 discrete π-contacts between the aromatic amino acids and the DNA nucleobases or deoxyribose were identified using strict criteria, including visual inspection. The abundance and structure of the interactions were determined by unequivocally classifying the contacts as either π-π stacking, π-π T-shaped or sugar-π contacts. Three hundred forty-four nucleobase-amino acid π-π contacts (60% of all interactions identified) were identified in 175 of the crystal structures searched. Unprecedented in the literature, 230 DNA-protein sugar-π contacts (40% of all interactions identified) were identified in 137 crystal structures, which involve C-H···π and/or lone-pair···π interactions, contain any amino acid and can be classified according to sugar atoms involved. Both π-π and sugar-π interactions display a range of relative monomer orientations and therefore interaction energies (up to -50 (-70) kJ mol(-1) for neutral (charged) interactions as determined using quantum chemical calculations). In general, DNA-protein π-interactions are more prevalent than perhaps currently accepted and the role of such interactions in many biological processes may yet to be uncovered.
Collapse
Affiliation(s)
- Katie A Wilson
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB, T1K 3M4, Canada
| | - Jennifer L Kellie
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB, T1K 3M4, Canada
| | - Stacey D Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB, T1K 3M4, Canada
| |
Collapse
|
30
|
An Y, Raju RK, Lu T, Wheeler SE. Aromatic interactions modulate the 5'-base selectivity of the DNA-binding autoantibody ED-10. J Phys Chem B 2014; 118:5653-9. [PMID: 24802982 DOI: 10.1021/jp502069a] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We present detailed computational analyses of the binding of four dinucleotides to a highly sequence-selective single-stranded DNA (ssDNA) binding antibody (ED-10) and selected point mutants. Anti-DNA antibodies are central to the pathogenesis of systemic lupus erythematosus (SLE), and a more complete understanding of the mode of binding of DNA and other ligands will be necessary to elucidate the role of anti-DNA antibodies in the kidney inflammation associated with SLE. Classical molecular mechanics based molecular dynamics simulations and density functional theory (DFT) computations were applied to pinpoint the origin of selectivity for the 5'-nucleotide. In particular, the strength of interactions between each nucleotide and the surrounding residues were computed using MMGBSA as well as DFT applied to a cluster model of the binding site. The results agree qualitatively with experimental binding free energies, and indicate that π-stacking, CH/π, NH/π, and hydrogen-bonding interactions all contribute to 5'-base selectivity in ED-10. Most importantly, the selectivity for dTdC over dAdC arises primarily from differences in the strength of π-stacking and XH/π interactions with the surrounding aromatic residues; hydrogen bonds play little role. These data suggest that a key Tyr residue, which is not present in other anti-DNA antibodies, plays a key role in the 5'-base selectivity, while we predict that the mutation of a single Trp residue can tune the selectivity for dTdC over dAdC.
Collapse
Affiliation(s)
- Yi An
- Department of Chemistry, Texas A&M University , College Station, Texas 77843, United States
| | | | | | | |
Collapse
|
31
|
Barone V, Cacelli I, Crescenzi O, d'Ischia M, Ferretti A, Prampolini G, Villani G. Unraveling the interplay of different contributions to the stability of the quinhydrone dimer. RSC Adv 2014. [DOI: 10.1039/c3ra46191b] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
|
32
|
Ruzziconi R, Bellachioma G, Ciancaleoni G, Lepri S, Superchi S, Zanasi R, Monaco G. Cationic half-sandwich quinolinophaneoxazoline-based (η6-p-cymene)ruthenium(ii) complexes exhibiting different chirality types: synthesis and structural determination by complementary spectroscopic methods. Dalton Trans 2014; 43:1636-50. [DOI: 10.1039/c3dt52291a] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
33
|
Copeland KL, Pellock SJ, Cox JR, Cafiero ML, Tschumper GS. Examination of tyrosine/adenine stacking interactions in protein complexes. J Phys Chem B 2013; 117:14001-8. [PMID: 24171662 DOI: 10.1021/jp408027j] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
The π-stacking interactions between tyrosine amino acid side chains and adenine-bearing ligands are examined. Crystalline protein structures from the protein data bank (PDB) exhibiting face-to-face tyrosine/adenine arrangements were used to construct 20 unique 4-methylphenol/N9-methyladenine (p-cresol/9MeA) model systems. Full geometry optimization of the 20 crystal structures with the M06-2X density functional theory method identified 11 unique low-energy conformations. CCSD(T) complete basis set (CBS) limit interaction energies were estimated for all of the structures to determine the magnitude of the interaction between the two ring systems. CCSD(T) computations with double-ζ basis sets (e.g., 6-31G*(0.25) and aug-cc-pVDZ) indicate that the MP2 method overbinds by as much as 3.07 kcal mol(-1) for the crystal structures and 3.90 kcal mol(-1) for the optimized structures. In the 20 crystal structures, the estimated CCSD(T) CBS limit interaction energy ranges from -4.00 to -6.83 kcal mol(-1), with an average interaction energy of -5.47 kcal mol(-1), values remarkably similar to the corresponding data for phenylalanine/adenine stacking interactions. Geometry optimization significantly increases the interaction energies of the p-cresol/9MeA model systems. The average estimated CCSD(T) CBS limit interaction energy of the 11 optimized structures is 3.23 kcal mol(-1) larger than that for the 20 crystal structures.
Collapse
Affiliation(s)
- Kari L Copeland
- Department of Chemistry and Biochemistry, University of Mississippi , University, Mississippi 38677, United States
| | | | | | | | | |
Collapse
|
34
|
Wells RA, Kellie JL, Wetmore SD. Significant strength of charged DNA-protein π-π interactions: a preliminary study of cytosine. J Phys Chem B 2013; 117:10462-74. [PMID: 23991905 DOI: 10.1021/jp406829d] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The present work characterized the preferred gas-phase structure and optimum interaction energy of both parallel stacked and perpendicular T-shaped dimers between cytosine (C), as a representative nucleobase, and aspartic/glutamic acid (DE), aspartate/glutamate (DE(-)) or arginine (R(+)), using detailed M06-2X/6-31+G(d,p) potential energy surface scans as a function of the relative monomer orientation. Through comparison to previous literature on the π-π interactions between the DNA nucleobases and the aromatic amino acid residues, this work will allow for comparisons between DNA-protein interactions involving aromatic and acyclic R-side chains, as well as comparisons of the relative geometric dependence and magnitude of π-π (C:DE), πcation-π (C:R(+)), and πanion-π (C:DE(-)) interactions. Our results show that the preferred relative monomer orientation is highly dependent on the monomer composition and charge, and is dictated by electrostatic-driven interactions. More importantly, for the first time, we report that the π-π interactions between cytosine and (neutral) aspartic/glutamic acid are up to approximately -40 kJ mol(-1), while the πcation-π or πanion-π interactions between cytosine and arginine or aspartate/glutamate are up to approximately -90 and -99 kJ mol(-1), respectively. An extensive investigation of the effects of the computational methodology implemented, including comparisons to detailed CCSD(T)/CBS potential energy surfaces and interaction energies, supports the use of M06-2X, as well as ωB97X-D, to study DNA-protein π-π interactions of varying composition and charge. Most importantly, the CCSD(T)/CBS results verify the strong nature of these DNA-protein π-π interactions, as well as the unique nature of the πcation-π and πanion-π counterparts. Therefore, our results emphasize that a wide variety of different types of noncovalent interactions between both cyclic and acyclic π-containing components can significantly contribute to the stability of DNA-protein complexes and likely play a larger role in biology than currently accepted.
Collapse
Affiliation(s)
- Rachael A Wells
- Department of Chemistry and Biochemistry, University of Lethbridge , 4401 University Drive West, Lethbridge, Alberta, Canada T1K 3M4
| | | | | |
Collapse
|
35
|
Badelin VG, Tyunina EY, Mezhevoi IN, Tarasova GN. Thermodynamic characteristics of the interaction between nicotinic acid and phenylalanine in an aqueous buffer solution at 298 K. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2013. [DOI: 10.1134/s0036024413070078] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
36
|
Kellie JL, Wetmore SD. Selecting DFT methods for use in optimizations of enzyme active sites: applications to ONIOM treatments of DNA glycosylases. CAN J CHEM 2013. [DOI: 10.1139/cjc-2012-0506] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
When using a hybrid methodology to treat an enzymatic reaction, many factors contribute to selecting the method for the high-level region, which can be complicated by the presence of dispersion-driven interactions such as π–π stacking. In addition, the proper treatment of the reaction center often requires a large number of heavy atoms to be included in the high-level region, precluding the use of ab initio methods such as MP2 as well as large basis sets, in the optimization step. In the present work, popular DFT methods were tested to identify an appropriate functional for treating the high-level region in ONIOM optimizations of reactions catalyzed by nonmetalloenzymes. Eight different DFT methods (B3LYP, B97-2, MPW1K, MPWB1K, BB1K, B1B95, M06-2X, and ωB97X-D) in combination with four double-ζ quality Pople basis sets were tested for their ability to optimize noncovalent interactions (hydrogen bonding and π–π) and characterize reactions (proton transfer, SN2 hydrolysis, and unimolecular cleavage). Although the primary focus of this study is accurate structure determination, energetics were also examined at both the optimization level of theory, and with triple-ζ quality basis set and select (M06-2X or ωB97X-D) methods. If dispersion-driven interactions exist within the active site, then MPWB1K/6-31G(d,p) or M06-2X/6-31+G(d,p) are recommended for the optimization step with subsequent triple-ζ quality single-point energies. However, since dispersion-corrected functionals (M06-2X and ωB97X-D) generally require diffuse functions to yield appropriate geometries, the possible size of the high-level region is greatly limited with these methods. In contrast, if the model is large enough to recover steric constraints on π–π interactions, then B3LYP with a small basis set performs comparatively well for the optimization step and is significantly less computationally expensive. Interestingly, the functionals that afford the best geometries often do not yield the best energetics, which emphasizes the importance of structural benchmark studies.
Collapse
Affiliation(s)
- Jennifer L. Kellie
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB T1K 3M4, Canada
| | - Stacey D. Wetmore
- Department of Chemistry and Biochemistry, University of Lethbridge, 4401 University Drive West, Lethbridge, AB T1K 3M4, Canada
| |
Collapse
|
37
|
Witts RN, Hopson EC, Koballa DE, Van Boening TA, Hopkins NH, Patterson EV, Nagan MC. Backbone-base interactions critical to quantum stabilization of transfer RNA anticodon structure. J Phys Chem B 2013; 117:7489-97. [PMID: 23742318 DOI: 10.1021/jp400084p] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Transfer RNA (tRNA) anticodons adopt a highly ordered 3'-stack without significant base overlap. Density functional theory at the M06-2X/6-31+G(d,p) level in combination with natural bond orbital analysis was utilized to calculate the intramolecular interactions within the tRNA anticodon that are responsible for stabilizing the stair-stepped conformation. Ten tRNA X-ray crystal structures were obtained from the PDB databank and were trimmed to include only the anticodon bases. Hydrogenic positions were added and optimized for the structures in the stair-stepped conformation. The sugar-phosphate backbone has been retained for these calculations, revealing the role it plays in RNA structural stability. It was found that electrostatic interactions between the sugar-phosphate backbone and the base provide the most stability, rather than the traditionally studied interbase stacking. Base-stacking interactions, though present, were weak and inconsistent. Aqueous solvation was found to have little effect on the intramolecular interactions.
Collapse
Affiliation(s)
- Rachel N Witts
- Department of Chemistry, Truman State University, 100 East Normal, Kirksville, Missouri 63501, USA
| | | | | | | | | | | | | |
Collapse
|
38
|
Riley KE, Hobza P. On the importance and origin of aromatic interactions in chemistry and biodisciplines. Acc Chem Res 2013; 46:927-36. [PMID: 22872015 DOI: 10.1021/ar300083h] [Citation(s) in RCA: 168] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Aromatic systems contain both σ- and π-electrons, which in turn constitute σ- and π-molecular orbitals (MOs). In discussing the properties of these systems, researchers typically refer to the highest occupied and lowest unoccupied MOs, which are π MOs. The characteristic properties of aromatic systems, such as their low ionization potentials and electron affinities, high polarizabilities and stabilities, and small band gaps (in spectroscopy called the N → V1 space), can easily be explained based on their electronic structure. These one-electron properties point to characteristic features of how aromatic systems interact with each other. Unlike hydrogen bonding systems, which primarily interact through electrostatic forces, complexes containing aromatic systems, especially aromatic stacked pairs, are predominantly stabilized by dispersion attraction. The stabilization energy in the benzene dimer is rather small (~2.5 kcal/mol) but strengthens with heteroatom substitution. The stacked interaction of aromatic nucleic acid bases is greater than 10 kcal/mol, and for the most stable stacked pair, guanine and cytosine, it reaches approximately 17 kcal/mol. Although these values do not equal the planar H-bonded interactions of these bases (~29 kcal/mol), stacking in DNA is more frequent than H-bonding and, unlike H-bonding, is not significantly weakened when passing from the gas phase to a water environment. Consequently, the stacking of aromatic systems represents the leading stabilization energy contribution in biomacromolecules and in related nanosystems. Therefore stacking (dispersion) interactions predominantly determine the double helical structure of DNA, which underlies its storage and transfer of genetic information. Similarly, dispersion is the dominant contributor to attractive interactions involving aromatic amino acids within the hydrophobic core of a protein, which is critical for folding. Therefore, understanding the nature of aromatic interactions, which depend greatly on quantum mechanical (QM) calculations, is of key importance in biomolecular science. This Account shows that accurate binding energies for aromatic complexes should be based on computations made at the (estimated) CCSD(T)/complete basis set limit (CBS) level of theory. This method is the least computationally intensive one that can give accurate stabilization energies for all common classes of noncovalent interactions (aromatic-aromatic, H-bonding, ionic, halogen bonding, charge-transfer, etc.). These results allow for direct comparison of binding energies between different interaction types. Conclusions based on lower-level QM calculations should be considered with care.
Collapse
Affiliation(s)
- Kevin E. Riley
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
| | - Pavel Hobza
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, 166 10 Prague, Czech Republic
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 771 46 Olomouc, Czech Republic
| |
Collapse
|
39
|
Yang J, Waller MP. JACOB: a dynamic database for computational chemistry benchmarking. J Chem Inf Model 2012; 52:3255-62. [PMID: 23157388 DOI: 10.1021/ci300374g] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
JACOB (just a collection of benchmarks) is a database that contains four diverse benchmark studies, which in-turn included 72 data sets, with a total of 122,356 individual results. The database is constructed upon a dynamic web framework that allows users to retrieve data from the database via predefined categories. Additional flexibility is made available via user-defined text-based queries. Requested sets of results are then automatically presented as bar graphs, with parameters of the graphs being controllable via the URL. JACOB is currently available at www.wallerlab.org/jacob.
Collapse
Affiliation(s)
- Jack Yang
- Theoretische Organische Chemie, Organisch-Chemisches Institut der Westfälische Wilhelms Universität Münster, Corrensstraße 40, 48149 Münster, Germany
| | | |
Collapse
|
40
|
Ramachandran S, Temple B, Alexandrova AN, Chaney SG, Dokholyan NV. Recognition of platinum-DNA adducts by HMGB1a. Biochemistry 2012; 51:7608-17. [PMID: 22950413 DOI: 10.1021/bi3008577] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Cisplatin (CP) and oxaliplatin (OX), platinum-based drugs used widely in chemotherapy, form adducts on intrastrand guanines (5'GG) in genomic DNA. DNA damage recognition proteins, transcription factors, mismatch repair proteins, and DNA polymerases discriminate between CP- and OX-GG DNA adducts, which could partly account for differences in the efficacy, toxicity, and mutagenicity of CP and OX. In addition, differential recognition of CP- and OX-GG adducts is highly dependent on the sequence context of the Pt-GG adduct. In particular, DNA binding protein domain HMGB1a binds to CP-GG DNA adducts with up to 53-fold greater affinity than to OX-GG adducts in the TGGA sequence context but shows much smaller differences in binding in the AGGC or TGGT sequence contexts. Here, simulations of the HMGB1a-Pt-DNA complex in the three sequence contexts revealed a higher number of interface contacts for the CP-DNA complex in the TGGA sequence context than in the OX-DNA complex. However, the number of interface contacts was similar in the TGGT and AGGC sequence contexts. The higher number of interface contacts in the CP-TGGA sequence context corresponded to a larger roll of the Pt-GG base pair step. Furthermore, geometric analysis of stacking of phenylalanine 37 in HMGB1a (Phe37) with the platinated guanines revealed more favorable stacking modes correlated with a larger roll of the Pt-GG base pair step in the TGGA sequence context. These data are consistent with our previous molecular dynamics simulations showing that the CP-TGGA complex was able to sample larger roll angles than the OX-TGGA complex or either CP- or OX-DNA complexes in the AGGC or TGGT sequences. We infer that the high binding affinity of HMGB1a for CP-TGGA is due to the greater flexibility of CP-TGGA compared to OX-TGGA and other Pt-DNA adducts. This increased flexibility is reflected in the ability of CP-TGGA to sample larger roll angles, which allows for a higher number of interface contacts between the Pt-DNA adduct and HMGB1a.
Collapse
Affiliation(s)
- Srinivas Ramachandran
- Department of Biochemistry and Biophysics, University of North Carolina, Chapel Hill, NC 27599, USA
| | | | | | | | | |
Collapse
|
41
|
|
42
|
Kellie JL, Navarro-Whyte L, Carvey MT, Wetmore SD. Combined effects of π-π stacking and hydrogen bonding on the (N1) acidity of uracil and hydrolysis of 2'-deoxyuridine. J Phys Chem B 2012; 116:2622-32. [PMID: 22296509 DOI: 10.1021/jp2121627] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
M06-2X/6-31+G(d,p) is used to study the simultaneous effects of π-π stacking interactions with phenylalanine (modeled as benzene) and hydrogen bonding with small molecules (HF, H(2)O, and NH(3)) on the N1 acidity of uracil and the hydrolytic deglycosylation of 2'-deoxyuridine (dU) (facilitated by fully (OH(-)) or partially (HCOO(-)···H(2)O) activated water). When phenylalanine is complexed with isolated uracil, the proton affinity of all acceptor sites significantly increases (by up to 28 kJ mol(-1)), while the N1 acidity slightly decreases (by ~6 kJ mol(-1)). When small molecules are hydrogen bound to uracil, addition of the phenylalanine ring can increase or decrease the acidity of uracil depending on the number and nature (acidity) of the molecules bound. Furthermore, a strong correlation between the effects of π-π stacking on the acidity of U and the dU deglycosylation reaction energetics is found, where the hydrolysis barrier can increase or decrease depending on the nature and number of small molecules bound, the nucleophile considered (which dictates the negative charge on U in the transition state), and the polarity of the (bulk) environment. These findings emphasize that the catalytic (or anticatalytic) role of the active-site aromatic amino acid residues is highly dependent on the situation under consideration. In the case of uracil-DNA glycosylase (UNG), which catalyzes the hydrolytic excision of uracil from DNA, the type of discrete hydrogen-bonding interactions with U, the nature of the nucleophile, and the anticipated weak, nonpolar environment in the active site suggest that phenylalanine will be slightly anticatalytic in the chemical step, and therefore experimentally observed contributions to catalysis may entirely result from associated structural changes that occur prior to deglycosylation.
Collapse
Affiliation(s)
- Jennifer L Kellie
- Department of Chemistry & Biochemistry, University of Lethbridge, 4401 University Drive, Lethbridge, Alberta T1K 3M4, Canada
| | | | | | | |
Collapse
|
43
|
Cacelli I, Cimoli A, Livotto PR, Prampolini G. An automated approach for the parameterization of accurate intermolecular force-fields: pyridine as a case study. J Comput Chem 2012; 33:1055-67. [PMID: 22410966 DOI: 10.1002/jcc.22937] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2011] [Revised: 12/23/2011] [Accepted: 12/27/2011] [Indexed: 11/09/2022]
Abstract
An automated protocol is proposed and validated, which integrates accurate quantum mechanical calculations with classical numerical simulations. Intermolecular force fields, (FF) suitable for molecular dynamics (MD) and Monte Carlo simulations, are parameterized through a novel iterative approach, fully based on quantum mechanical data, which has been automated and coded into the PICKY software, here presented. The whole procedure is tested and validated for pyridine, whose bulk phase, described through MD simulations performed with the specifically parameterized FF, is characterized by computing several of its thermodynamic, structural, and transport properties, comparing them with their experimental counterparts.
Collapse
Affiliation(s)
- Ivo Cacelli
- Dipartimento di Chimica e Chimica Industriale, Università degli Studi di Pisa, via Risorgimento 35, I-56126 Pisa, Italy
| | | | | | | |
Collapse
|
44
|
Rutledge LR, Wetmore SD. A computational proposal for the experimentally observed discriminatory behavior of hypoxanthine, a weak universal nucleobase. Phys Chem Chem Phys 2012; 14:2743-53. [PMID: 22270716 DOI: 10.1039/c2cp23600a] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
A computational model composed of six nucleobases was used to investigate why hypoxanthine does not yield duplexes of equal stability when paired opposite each of the natural DNA nucleobases. The magnitudes of all nearest-neighbor interactions in a DNA helix were calculated, including hydrogen-bonding, intra- and interstrand stacking interactions, as well as 1-3 intrastrand stacking interactions. Although the stacking interactions in DNA relevant arrangements are significant and account for at least one third of the total stabilization energy in our nucleobase complexes, the trends in the magnitude of the stacking interactions cannot explain the relative experimental melting temperatures previously reported in the literature. Furthermore, although the total hydrogen-bonding interactions explain why hypoxanthine preferentially pairs with cytosine, the experimental trend for the remaining nucleobases (A, T, G) is not explained. In fact, the calculated pairing preference of hypoxanthine matches that determined experimentally only when the sum of all types of nearest-neighbor interactions is considered. This finding highlights a strong correlation between the relative magnitude of the total nucleobase-nucleobase interactions and measured melting temperatures for DNA strands containing hypoxanthine despite the potential role of other factors (including hydration, temperature, sugar-phosphate backbone). By considering a large range of sequence combinations, we reveal that the binding preference of hypoxanthine is strongly dependent on the nucleobase sequence, which may explain the varied ability of hypoxanthine to universally bind to the natural bases. As a result, we propose that future work should closely examine the interplay between the dominant nucleobase-nucleobase interactions and the overall strand stability to fully understand how sequence context affects the universal binding properties of modified bases and to aid the design of new molecules with ambiguous pairing properties.
Collapse
Affiliation(s)
- Lesley R Rutledge
- Department of Chemistry & Biochemistry, University of Lethbridge, 4401 University Drive, Lethbridge, Alberta, Canada T1K 3M4
| | | |
Collapse
|
45
|
|
46
|
JISSY AK, RAMANA JHV, DATTA AYAN. π-Stacking interactions between G-quartets and circulenes: A computational study. J CHEM SCI 2011. [DOI: 10.1007/s12039-011-0155-7] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
47
|
Marshall MS, Burns LA, Sherrill CD. Basis set convergence of the coupled-cluster correction, δMP2CCSD(T): Best practices for benchmarking non-covalent interactions and the attendant revision of the S22, NBC10, HBC6, and HSG databases. J Chem Phys 2011; 135:194102. [DOI: 10.1063/1.3659142] [Citation(s) in RCA: 263] [Impact Index Per Article: 20.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
|
48
|
Butchosa C, Simon S, Voityuk AA. Conformational dependence of the electronic coupling for hole transfer between adenine and tryptophan. COMPUT THEOR CHEM 2011. [DOI: 10.1016/j.comptc.2011.04.025] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
49
|
Rutledge LR, Wetmore SD. Modeling the chemical step utilized by human alkyladenine DNA glycosylase: a concerted mechanism AIDS in selectively excising damaged purines. J Am Chem Soc 2011; 133:16258-69. [PMID: 21877721 DOI: 10.1021/ja207181c] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Human alkyladenine DNA glycosylase (AAG) initiates the repair of a wide variety of (neutral or cationic) alkylated and deaminated purines by flipping damaged nucleotides out of the DNA helix and catalyzing the hydrolytic N-glycosidic bond cleavage. Unfortunately, the limited number of studies on the catalytic pathway has left many unanswered questions about the hydrolysis mechanism. Therefore, detailed ONIOM(M06-2X/6-31G(d):AMBER) reaction potential energy surface scans are used to gain the first atomistic perspective of the repair pathway used by AAG. The lowest barrier for neutral 1,N(6)-ethenoadenine (εA) and cationic N(3)-methyladenine (3MeA) excision corresponds to a concerted (A(N)D(N)) mechanism, where our calculated ΔG(‡) = 87.3 kJ mol(-1) for εA cleavage is consistent with recent kinetic data. The use of a concerted mechanism supports previous speculations that AAG uses a nonspecific strategy to excise both neutral (εA) and cationic (3MeA) lesions. We find that AAG uses nonspecific active site DNA-protein π-π interactions to catalyze the removal of inherently more difficult to excise neutral lesions, and strongly bind to cationic lesions, which comes at the expense of raising the excision barrier for cationic substrates. Although proton transfer from the recently proposed general acid (protein-bound water) to neutral substrates does not occur, hydrogen-bond donation lowers the catalytic barrier, which clarifies the role of a general acid in the excision of neutral lesions. Finally, our work shows that the natural base adenine (A) is further inserted into the AAG active site than the damaged substrates, which results in the loss of a hydrogen bond with Y127 and misaligns the general base (E125) and water nucleophile to lead to poor nucleophile activation. Therefore, our work proposes how AAG discriminates against the natural purines in the chemical step and may also explain why some damaged pyrimidines are bound but are not excised by this enzyme.
Collapse
Affiliation(s)
- Lesley R Rutledge
- Department of Chemistry and Biochemistry, University of Lethbridge, Alberta T1K 3M4, Canada
| | | |
Collapse
|
50
|
Leavens FMV, Churchill CDM, Wang S, Wetmore SD. Evaluating how discrete water molecules affect protein-DNA π-π and π(+)-π stacking and T-shaped interactions: the case of histidine-adenine dimers. J Phys Chem B 2011; 115:10990-1003. [PMID: 21809837 DOI: 10.1021/jp205424z] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Changes in the magnitude of (M06-2X/6-31+G(d,p)) π-π stacking and T-shaped (nucleobase-edge and amino acid-edge) interactions between (neutral or protonated) histidine (His) and adenine (A) dimers upon microsolvation with up to four discrete water molecules were determined. A variety of histidine-water interactions were considered including conventional (N-H···O, N···H-O, C-H···O) hydrogen bonding and nonconventional (X-H···π (neutral His) or lone-pair···π (protonated His)) contacts. Overall, the effects of discrete His-H(2)O interactions on the neutral histidine-adenine π-π contacts are negligible (<3 kJ mol(-1) or 15%) regardless of the type of water binding, the number of water molecules bound, or the His-A dimer (stacked or (amino acid- or nucleobase-edge) T-shaped) configuration. This suggests that previously reported gas-phase binding strengths for a variety of neutral amino acid-nucleobase dimers are likely relevant for a wide variety of (microsolvated) environments. In contrast, the presence of water decreases the histidine-adenine π(+)-π interaction by up to 15 kJ mol(-1) (or 30%) for all water binding modes and orientations, as well as different stacked and T-shaped His(+)-A dimers. Regardless of the larger effect of discrete histidine-water interactions on the magnitude of the π(+)-π compared with π-π interactions, the π(+)-π binding strengths remain substantially larger than the corresponding π-π contacts. These findings emphasize the distinct nature of π(+)-π and π-π interactions and suggest that π(+)-π contacts can provide significant stabilization in biological systems relative to π-π contacts under many different environmental conditions.
Collapse
Affiliation(s)
- Fern M V Leavens
- Department of Chemistry & Biochemistry, University of Lethbridge, 4401 University Drive, Lethbridge, Alberta, Canada, T1K 3M4
| | | | | | | |
Collapse
|