1
|
Calcagno F, Maryasin B, Garavelli M, Avagliano D, Rivalta I. Modeling solvent effects and convergence of 31P-NMR shielding calculations with COBRAMM. J Comput Chem 2024; 45:1562-1575. [PMID: 38514234 DOI: 10.1002/jcc.27338] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2023] [Revised: 02/15/2024] [Accepted: 02/22/2024] [Indexed: 03/23/2024]
Abstract
Solvent effects on 31P-NMR parameters for triphenylphosphine oxide and triphenylphosphine in chloroform have been extensively investigated by testing different solvation models. The solvent is described implicitly, mixed implicitly/explicitly, and using full explicit models. Polarizable continuum model (PCM), molecular dynamic (MD) simulations, and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations are used to disclose the effects of solute/solvent interactions and, more generally, the role of the embedding in NMR simulations. The results show the beneficial effect of carrying out QM/MM optimizations on top of geometries directly extracted from classical MD simulations, used to ensure representative conformational sampling. The nuclear shielding convergence has been tested against a different number of snapshots and with the inclusion of solvent shells into the QM region. An automated MD//QM/MM//GIAO protocol, implemented in the COBRAMM package, is here proposed and tested on trimethyl phosphite showing that our approach boosts the convergence of nuclear shielding satisfactorily. The present work aims to be a stepping-stone to assess proper QM/MM computational strategies in simulating chemical shifts in non-homogeneous systems like supramolecular and biological systems.
Collapse
Affiliation(s)
- Francesco Calcagno
- Department of Industrial Chemistry "Toso Montanari", University of Bologna, Bologna, Italy
- Center for Chemical Catalysis - C3, University of Bologna, Bologna, Italy
| | - Boris Maryasin
- Institute of Organic Chemistry, University of Vienna, Vienna, Austria
- Institute of Theoretical Chemistry, University of Vienna, Vienna, Austria
| | - Marco Garavelli
- Department of Industrial Chemistry "Toso Montanari", University of Bologna, Bologna, Italy
| | - Davide Avagliano
- Department of Industrial Chemistry "Toso Montanari", University of Bologna, Bologna, Italy
| | - Ivan Rivalta
- Department of Industrial Chemistry "Toso Montanari", University of Bologna, Bologna, Italy
- Center for Chemical Catalysis - C3, University of Bologna, Bologna, Italy
- ENSL, CNRS, Laboratoire de Chimie UMR 5182, Lyon, France
| |
Collapse
|
2
|
Szántó JK, Dietschreit JCB, Shein M, Schütz AK, Ochsenfeld C. Systematic QM/MM Study for Predicting 31P NMR Chemical Shifts of Adenosine Nucleotides in Solution and Stages of ATP Hydrolysis in a Protein Environment. J Chem Theory Comput 2024; 20:2433-2444. [PMID: 38497488 DOI: 10.1021/acs.jctc.3c01280] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2024]
Abstract
NMR (nuclear magnetic resonance) spectroscopy allows for important atomistic insights into the structure and dynamics of biological macromolecules; however, reliable assignments of experimental spectra are often difficult. Herein, quantum mechanical/molecular mechanical (QM/MM) calculations can provide crucial support. A major problem for the simulations is that experimental NMR signals are time-averaged over much longer time scales, and since computed chemical shifts are highly sensitive to local changes in the electronic and structural environment, sufficiently large averages over representative structural ensembles are essential. This entails high computational demands for reliable simulations. For NMR measurements in biological systems, a nucleus of major interest is 31P since it is both highly present (e.g., in nucleic acids) and easily observable. The focus of our present study is to develop a robust and computationally cost-efficient framework for simulating 31P NMR chemical shifts of nucleotides. We apply this scheme to study the different stages of the ATP hydrolysis reaction catalyzed by p97. Our methodology is based on MM molecular dynamics (MM-MD) sampling, followed by QM/MM structure optimizations and NMR calculations. Overall, our study is one of the most comprehensive QM-based 31P studies in a protein environment and the first to provide computed NMR chemical shifts for multiple nucleotide states in a protein environment. This study sheds light on a process that is challenging to probe experimentally and aims to bridge the gap between measured and calculated NMR spectroscopic properties.
Collapse
Affiliation(s)
- Judit Katalin Szántó
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
| | - Johannes C B Dietschreit
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Mikhail Shein
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 München, Germany
| | - Anne K Schütz
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 München, Germany
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, D-70569 Stuttgart, Germany
| |
Collapse
|
3
|
Hersh WH, Chan TY. Improving the accuracy of 31P NMR chemical shift calculations by use of scaling methods. Beilstein J Org Chem 2023; 19:36-56. [PMID: 36726479 PMCID: PMC9843238 DOI: 10.3762/bjoc.19.4] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2022] [Accepted: 12/14/2022] [Indexed: 01/12/2023] Open
Abstract
Calculation of 31P NMR chemical shifts for a series of tri- and tetracoordinate phosphorus compounds using several basis sets and density functional theory (DFT) functionals gave a modest fit to experimental chemical shifts, but an excellent linear fit when plotted against the experimental values. The resultant scaling methods were then applied to a variety of "large" compounds previously selected by Latypov et al. and a set of stereoisomeric and unusual compounds selected here. No one method was best for all structural types. For compounds that contain P-P bonds and P-C multiple bonds, the Latypov et al. method using the PBE0 functional was best (mean absolute deviation/root mean square deviation (MAD/RMSD) = 6.9/8.5 ppm and 6.6/8.2 ppm, respectively), but for the full set of compounds gave higher deviations (MAD/RMSD = 8.2/12.3 ppm), and failed by over 60 ppm for a three-membered phosphorus heterocycle. Use of the M06-2X functional for both the structural optimization and NMR chemical shift calculation was best overall for the compounds without P-C multiple bonds (MAD/RMSD = 5.4/7.1 ppm), but failed by 30-49 ppm for compounds having any P-C multiple-bond character. Failures of these magnitudes have not been reported previously for these widely used functionals. These failures were then used to screen a variety of recommended functionals, leading to better overall methods for calculation of these chemical shifts: optimization with the M06-2X functional and NMR calculation with the PBE0 or ωB97x-D functionals gave values for MAD/RMSD = 6.9/8.5 ppm and 6.8/9.1 ppm, respectively, over an experimental chemical shift range of -181 to 356 ppm. Due to the unexplained failures observed, we recommend use of more than one method when looking at novel structures.
Collapse
Affiliation(s)
- William H Hersh
- Department of Chemistry and Biochemistry, Queens College, Queens, NY 11367-1597, USA,Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, USA
| | - Tsz-Yeung Chan
- Department of Chemistry and Biochemistry, Queens College, Queens, NY 11367-1597, USA,Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, NY 10016, USA
| |
Collapse
|
4
|
Bakker MJ, Mládek A, Semrád H, Zapletal V, Pavlíková Přecechtělová J. Improving IDP theoretical chemical shift accuracy and efficiency through a combined MD/ADMA/DFT and machine learning approach. Phys Chem Chem Phys 2022; 24:27678-27692. [PMID: 36373847 DOI: 10.1039/d2cp01638a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
This work extends the multi-scale computational scheme for the quantum mechanics (QM) calculations of Nuclear Magnetic Resonance (NMR) chemical shifts (CSs) in proteins that lack a well-defined 3D structure. The scheme couples the sampling of an intrinsically disordered protein (IDP) by classical molecular dynamics (MD) with protein fragmentation using the adjustable density matrix assembler (ADMA) and density functional theory (DFT) calculations. In contrast to our early investigation on IDPs (Pavlíková Přecechtělová et al., J. Chem. Theory Comput., 2019, 15, 5642-5658) and the state-of-the art NMR calculations for structured proteins, a partial re-optimization was implemented on the raw MD geometries in vibrational normal mode coordinates to enhance the accuracy of the MD/ADMA/DFT computational scheme. In addition, machine-learning based cluster analysis was performed on the scheme to explore its potential in producing protein structure ensembles (CLUSTER ensembles) that yield accurate CSs at a reduced computational cost. The performance of the cluster-based calculations is validated against results obtained with conventional structural ensembles consisting of MD snapshots extracted from the MD trajectory at regular time intervals (REGULAR ensembles). CS calculations performed with the refined MD/ADMA/DFT framework employed the 6-311++G(d,p) basis set that outperformed IGLO-III calculations with the same density functional approximation (B3LYP) and both explicit and implicit solvation. The partial geometry optimization did not universally improve the agreement of computed CSs with the experiment but substantially decreased errors associated with the ensemble averaging. A CLUSTER ensemble with 50 structures yielded ensemble averages close to those obtained with a REGULAR ensemble consisting of 500 MD frames. The cluster based calculations thus required only a fraction of the computational time.
Collapse
Affiliation(s)
- Michael J Bakker
- Faculty of Pharmacy in Hradec Králové, Charles University, Akademika Heyrovského 1203/8, 500 05 Hradec Králové, Czech Republic.
| | - Arnošt Mládek
- Faculty of Pharmacy in Hradec Králové, Charles University, Akademika Heyrovského 1203/8, 500 05 Hradec Králové, Czech Republic.
| | - Hugo Semrád
- Faculty of Pharmacy in Hradec Králové, Charles University, Akademika Heyrovského 1203/8, 500 05 Hradec Králové, Czech Republic. .,Department of Chemistry, Faculty of Science, Masaryk University, Kotlářská 267/2, 611 37 Brno, Czech Republic
| | - Vojtěch Zapletal
- Faculty of Pharmacy in Hradec Králové, Charles University, Akademika Heyrovského 1203/8, 500 05 Hradec Králové, Czech Republic.
| | - Jana Pavlíková Přecechtělová
- Faculty of Pharmacy in Hradec Králové, Charles University, Akademika Heyrovského 1203/8, 500 05 Hradec Králové, Czech Republic.
| |
Collapse
|
5
|
Fukal J, Zgarbová M, Jurečka P, Šebera J, Sychrovský V. Probabilistic Interpretation of NMR J-Couplings Determines BI-BII State Equilibria in DNA. J Chem Theory Comput 2022; 18:6989-6999. [PMID: 36206364 DOI: 10.1021/acs.jctc.2c00733] [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
Interpretation of 3JP,H3' NMR scalar spin-spin coupling constants in DNA becomes more reliable by including distinct structural states such as BI and BII, using the weighted-static or, better still, the recently implemented adiabatic-MD (Ad-MD) method. The calculation method employs an adiabatic ("Ad") dependence of 3JP,H3' coupling on NMR-assigned torsion angle, ε, weighted by P(ε) probability distribution calculated by molecular dynamics (MD). Ad-MD calculations enable cross-validation of the bsc1, OL15, and OL21 force fields and various parametrizations of the Karplus equation describing the dependence of 3JP,H3' coupling on ε torsion (KE). The mean absolute deviation of Ad-MD 3JP,H3' couplings from the experimental values in Dickerson-Drew DNA is comparable to the scatter of 3JP,H3' couplings among four separate NMR experiments. A commonly accepted assumption of homogeneity of one kind of structure-dynamic state within DNA (BI or BII) is questionable because the principal characteristics of relevant P(ε) probabilities (shapes and positioning) vary with DNA sequence. The theory outlined in the present work sets limits to future reparameterization of MD force fields, as relevant to NMR data.
Collapse
Affiliation(s)
- Jiří Fukal
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic.,Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Marie Zgarbová
- Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Petr Jurečka
- Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Jakub Šebera
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic
| | - Vladimír Sychrovský
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Flemingovo náměstí 2, 166 10 Praha 6, Czech Republic.,Department of Electrotechnology, Electrical Engineering, Czech Technical University, Technická 2, 166 27 Praha 6, Czech Republic
| |
Collapse
|
6
|
Fukal J, Buděšínský M, Páv O, Jurečka P, Zgarbová M, Šebera J, Sychrovský V. The Ad-MD method to calculate NMR shift including effects due to conformational dynamics: The 31 P NMR shift in DNA. J Comput Chem 2022; 43:132-143. [PMID: 34729803 DOI: 10.1002/jcc.26778] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Revised: 10/12/2021] [Accepted: 10/12/2021] [Indexed: 11/11/2022]
Abstract
A method for averaging of NMR parameters by molecular dynamics (MD) has been derived from the method of statistical averaging in MD snapshots, benchmarked and applied to structurally dynamic interpretation of the 31 P NMR shift (δ31P ) in DNA phosphates. The method employs adiabatic dependence of an NMR parameter on selected geometric parameter(s) that is weighted by MD-calculated probability distribution(s) for the geometric parameter(s) (Ad-MD method). The usage of Ad-MD for polymers is computationally convenient when one pre-calculated structural dependence of an NMR parameter is employed for all chemically equivalent units differing only in dynamic behavior. The Ad-MD method is benchmarked against the statistical averaging method for δ31P in the model phosphates featuring distinctively different structures and dynamic behavior. The applicability of Ad-MD is illustrated by calculating 31 P NMR spectra in the Dickerson-Drew DNA dodecamer. δ31P was calculated with the B3LYP/IGLO-III/PCM(water) and the probability distributions for the torsion angles adjacent to the phosphorus atoms in the DNA phosphates were calculated using the OL15 force field.
Collapse
Affiliation(s)
- Jiří Fukal
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic.,Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
| | - Miloš Buděšínský
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic
| | - Ondřej Páv
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic
| | - Petr Jurečka
- Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, Olomouc, Czech Republic
| | - Marie Zgarbová
- Department of Physical Chemistry, Faculty of Science, Palacký University Olomouc, Olomouc, Czech Republic
| | - Jakub Šebera
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic
| | - Vladimír Sychrovský
- Institute of Organic Chemistry and Biochemistry, Czech Academy of Sciences, Prague, Czech Republic.,Department of Electrotechnology, Faculty of Electrical Engineering, Czech Technical University, Prague, Czech Republic
| |
Collapse
|
7
|
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
|
8
|
Wilson AL, Outeiral C, Dowd SE, Doig AJ, Popelier PLA, Waltho JP, Almond A. Deconvolution of conformational exchange from Raman spectra of aqueous RNA nucleosides. Commun Chem 2020; 3:56. [PMID: 36703475 PMCID: PMC9814580 DOI: 10.1038/s42004-020-0298-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2019] [Accepted: 04/06/2020] [Indexed: 01/29/2023] Open
Abstract
Ribonucleic acids (RNAs) are key to the central dogma of molecular biology. While Raman spectroscopy holds great potential for studying RNA conformational dynamics, current computational Raman prediction and assignment methods are limited in terms of system size and inclusion of conformational exchange. Here, a framework is presented that predicts Raman spectra using mixtures of sub-spectra corresponding to major conformers calculated using classical and ab initio molecular dynamics. Experimental optimization allowed purines and pyrimidines to be characterized as predominantly syn and anti, respectively, and ribose into exchange between equivalent south and north populations. These measurements are in excellent agreement with Raman spectroscopy of ribonucleosides, and previous experimental and computational results. This framework provides a measure of ribonucleoside solution populations and conformational exchange in RNA subunits. It complements other experimental techniques and could be extended to other molecules, such as proteins and carbohydrates, enabling biological insights and providing a new analytical tool.
Collapse
Affiliation(s)
- Alex L. Wilson
- grid.5379.80000000121662407Manchester Institute of Biotechnology and Department of Chemistry, School of Natural Science, Faculty of Science and Engineering, The University of Manchester, M1 7DN Manchester, UK
| | - Carlos Outeiral
- grid.5379.80000000121662407Manchester Institute of Biotechnology and Department of Chemistry, School of Natural Science, Faculty of Science and Engineering, The University of Manchester, M1 7DN Manchester, UK
| | - Sarah E. Dowd
- grid.5379.80000000121662407Manchester Institute of Biotechnology and Department of Chemistry, School of Natural Science, Faculty of Science and Engineering, The University of Manchester, M1 7DN Manchester, UK
| | - Andrew J. Doig
- grid.5379.80000000121662407Division of Neuroscience and Experimental Psychology, Michael Smith Building, School of Biological Sciences, Faculty of Biology, Medicine and Health, The University of Manchester, M13 9PT Manchester, UK
| | - Paul L. A. Popelier
- grid.5379.80000000121662407Manchester Institute of Biotechnology and Department of Chemistry, School of Natural Science, Faculty of Science and Engineering, The University of Manchester, M1 7DN Manchester, UK
| | - Jonathan P. Waltho
- grid.5379.80000000121662407Manchester Institute of Biotechnology and Department of Chemistry, School of Natural Science, Faculty of Science and Engineering, The University of Manchester, M1 7DN Manchester, UK ,grid.11835.3e0000 0004 1936 9262Krebs Institute for Biomolecular Research, Department of Molecular Biology and Biotechnology, The University of Sheffield, S10 2TN Sheffield, UK
| | - Andrew Almond
- grid.5379.80000000121662407Manchester Institute of Biotechnology and Department of Chemistry, School of Natural Science, Faculty of Science and Engineering, The University of Manchester, M1 7DN Manchester, UK
| |
Collapse
|
9
|
Pavlíková Přecechtělová J, Mládek A, Zapletal V, Hritz J. Quantum Chemical Calculations of NMR Chemical Shifts in Phosphorylated Intrinsically Disordered Proteins. J Chem Theory Comput 2019; 15:5642-5658. [DOI: 10.1021/acs.jctc.8b00257] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Jana Pavlíková Přecechtělová
- Faculty of Pharmacy in Hradec Králové, Charles University, Akademika Heyrovského 1203, 500 05 Hradec Králové, Czech Republic
| | | | | | | |
Collapse
|
10
|
Fukal J, Páv O, Buděšínský M, Rosenberg I, Šebera J, Sychrovský V. Structural interpretation of the 31P NMR chemical shifts in thiophosphate and phosphate: key effects due to spin–orbit and explicit solvent. Phys Chem Chem Phys 2019; 21:9924-9934. [DOI: 10.1039/c9cp01460h] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Structural interpretation of the 31P NMR shifts measured in different molecules including thiophosphate or phosphate group was obtained by means of theoretical calculations including the effects of geometry, molecular dynamics, solvent, relativistic effects and the effect of NMR reference.
Collapse
Affiliation(s)
- J. Fukal
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- Praha 6
- Czech Republic
| | - O. Páv
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- Praha 6
- Czech Republic
| | - M. Buděšínský
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- Praha 6
- Czech Republic
| | - I. Rosenberg
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- Praha 6
- Czech Republic
| | - J. Šebera
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- Praha 6
- Czech Republic
| | - V. Sychrovský
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- Praha 6
- Czech Republic
- Department of Electrotechnology
| |
Collapse
|
11
|
Fukal J, Páv O, Buděšínský M, Šebera J, Sychrovský V. The benchmark of 31P NMR parameters in phosphate: a case study on structurally constrained and flexible phosphate. Phys Chem Chem Phys 2017; 19:31830-31841. [DOI: 10.1039/c7cp06969c] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
A benchmark for structural interpretation of the 31P NMR shift and the 2JP,C spin–spin coupling in the phosphate group was obtained by means of theoretical calculations and measurements in diethylphosphate and 5,5-dimethyl-2-hydroxy-1,3,2-dioxaphosphinane 2-oxide.
Collapse
Affiliation(s)
- Jiří Fukal
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- Praha 6
- Czech Republic
| | - Ondřej Páv
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- Praha 6
- Czech Republic
| | - Miloš Buděšínský
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- Praha 6
- Czech Republic
| | - Jakub Šebera
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- Praha 6
- Czech Republic
| | - Vladimír Sychrovský
- Institute of Organic Chemistry and Biochemistry
- Academy of Sciences of the Czech Republic
- Praha 6
- Czech Republic
- Department of Electrotechnology
| |
Collapse
|
12
|
Vícha J, Babinský M, Demo G, Otrusinová O, Jansen S, Pekárová B, Žídek L, Munzarová ML. The influence of Mg2+ coordination on 13C and 15N chemical shifts in CKI1RD protein domain from experiment and molecular dynamics/density functional theory calculations. Proteins 2016; 84:686-99. [PMID: 26879585 DOI: 10.1002/prot.25019] [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] [Received: 12/10/2015] [Revised: 01/29/2016] [Accepted: 02/10/2016] [Indexed: 11/09/2022]
Abstract
Sequence dependence of (13) C and (15) N chemical shifts in the receiver domain of CKI1 protein from Arabidopsis thaliana, CKI1RD , and its complexed form, CKI1RD •Mg(2+), was studied by means of MD/DFT calculations. MD simulations of a 20-ns production run length were performed. Nine explicitly hydrated structures of increasing complexity were explored, up to a 40-amino-acid structure. The size of the model necessary depended on the type of nucleus, the type of amino acid and its sequence neighbors, other spatially close amino acids, and the orientation of amino acid NH groups and their surface/interior position. Using models covering a 10 and a 15 Å environment of Mg(2+), a semi-quantitative agreement has been obtained between experiment and theory for the V67-I73 sequence. The influence of Mg(2+) binding was described better by the 15 Å as compared to the 10 Å model. Thirteen chemical shifts were analyzed in terms of the effect of Mg(2+) insertion and geometry preparation. The effect of geometry was significant and opposite in sign to the effect of Mg(2+) binding. The strongest individual effects were found for (15) N of D70, S74, and V68, where the electrostatics dominated; for (13) Cβ of D69 and (15) N of K76, where the influences were equal, and for (13) Cα of F72 and (13) Cβ of K76, where the geometry adjustment dominated. A partial correlation between dominant geometry influence and torsion angle shifts upon the coordination has been observed.
Collapse
Affiliation(s)
- Jan Vícha
- Central European Institute of Technology, Masaryk University, Kamenice 5, Brno, 62500, Czech Republic
| | - Martin Babinský
- Central European Institute of Technology, Masaryk University, Kamenice 5, Brno, 62500, Czech Republic
| | - Gabriel Demo
- Central European Institute of Technology, Masaryk University, Kamenice 5, Brno, 62500, Czech Republic
| | - Olga Otrusinová
- Central European Institute of Technology, Masaryk University, Kamenice 5, Brno, 62500, Czech Republic.,Faculty of Science, National Centre for Biomolecular Research, Masaryk University, Kamenice 5, Brno, 62500, Czech Republic
| | - Séverine Jansen
- Central European Institute of Technology, Masaryk University, Kamenice 5, Brno, 62500, Czech Republic
| | - Blanka Pekárová
- Central European Institute of Technology, Masaryk University, Kamenice 5, Brno, 62500, Czech Republic
| | - Lukáš Žídek
- Central European Institute of Technology, Masaryk University, Kamenice 5, Brno, 62500, Czech Republic.,Faculty of Science, National Centre for Biomolecular Research, Masaryk University, Kamenice 5, Brno, 62500, Czech Republic
| | - Markéta L Munzarová
- Faculty of Science, Department of Chemistry, Masaryk University, Kotlářská 2, Brno, 61137, Czech Republic
| |
Collapse
|
13
|
Imeddourene AB, Xu X, Zargarian L, Oguey C, Foloppe N, Mauffret O, Hartmann B. The intrinsic mechanics of B-DNA in solution characterized by NMR. Nucleic Acids Res 2016; 44:3432-47. [PMID: 26883628 PMCID: PMC4838374 DOI: 10.1093/nar/gkw084] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2015] [Accepted: 02/01/2016] [Indexed: 12/19/2022] Open
Abstract
Experimental characterization of the structural couplings in free B-DNA in solution has been elusive, because of subtle effects that are challenging to tackle. Here, the exploitation of the NMR measurements collected on four dodecamers containing a substantial set of dinucleotide sequences provides new, consistent correlations revealing the DNA intrinsic mechanics. The difference between two successive residual dipolar couplings (ΔRDCs) involving C6/8-H6/8, C3′-H3′ and C4′-H4′ vectors are correlated to the 31P chemical shifts (δP), which reflect the populations of the BI and BII backbone states. The δPs are also correlated to the internucleotide distances (Dinter) involving H6/8, H2′ and H2″ protons. Calculations of NMR quantities on high resolution X-ray structures and controlled models of DNA enable to interpret these couplings: the studied ΔRDCs depend mostly on roll, while Dinter are mainly sensitive to twist or slide. Overall, these relations demonstrate how δP measurements inform on key inter base parameters, in addition to probe the BI↔BII backbone equilibrium, and shed new light into coordinated motions of phosphate groups and bases in free B-DNA in solution. Inspection of the 5′ and 3′ ends of the dodecamers also supplies new information on the fraying events, otherwise neglected.
Collapse
Affiliation(s)
- Akli Ben Imeddourene
- Laboratoire de Biologie et Pharmacologie Appliquée, ENS Cachan, CNRS, Université Paris-Saclay, 61 avenue du Président Wilson, 94235 Cachan cedex, France Université Pierre et Marie Curie, 4 Place Jussieu, 75005 Paris, France
| | - Xiaoqian Xu
- Laboratoire de Biologie et Pharmacologie Appliquée, ENS Cachan, CNRS, Université Paris-Saclay, 61 avenue du Président Wilson, 94235 Cachan cedex, France Department of Life Sciences, East China Normal University, 200062 Shanghai, People's Republic of China
| | - Loussiné Zargarian
- Laboratoire de Biologie et Pharmacologie Appliquée, ENS Cachan, CNRS, Université Paris-Saclay, 61 avenue du Président Wilson, 94235 Cachan cedex, France
| | - Christophe Oguey
- Laboratoire de Physique Théorique et Modélisation, UMR 8089, CNRS, Université de Cergy-Pontoise, Cergy-Pontoise, France
| | | | - Olivier Mauffret
- Laboratoire de Biologie et Pharmacologie Appliquée, ENS Cachan, CNRS, Université Paris-Saclay, 61 avenue du Président Wilson, 94235 Cachan cedex, France
| | - Brigitte Hartmann
- Laboratoire de Biologie et Pharmacologie Appliquée, ENS Cachan, CNRS, Université Paris-Saclay, 61 avenue du Président Wilson, 94235 Cachan cedex, France
| |
Collapse
|
14
|
Ben Imeddourene A, Elbahnsi A, Guéroult M, Oguey C, Foloppe N, Hartmann B. Simulations Meet Experiment to Reveal New Insights into DNA Intrinsic Mechanics. PLoS Comput Biol 2015; 11:e1004631. [PMID: 26657165 PMCID: PMC4689557 DOI: 10.1371/journal.pcbi.1004631] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2015] [Accepted: 10/28/2015] [Indexed: 01/30/2023] Open
Abstract
The accurate prediction of the structure and dynamics of DNA remains a major challenge in computational biology due to the dearth of precise experimental information on DNA free in solution and limitations in the DNA force-fields underpinning the simulations. A new generation of force-fields has been developed to better represent the sequence-dependent B-DNA intrinsic mechanics, in particular with respect to the BI ↔ BII backbone equilibrium, which is essential to understand the B-DNA properties. Here, the performance of MD simulations with the newly updated force-fields Parmbsc0εζOLI and CHARMM36 was tested against a large ensemble of recent NMR data collected on four DNA dodecamers involved in nucleosome positioning. We find impressive progress towards a coherent, realistic representation of B-DNA in solution, despite residual shortcomings. This improved representation allows new and deeper interpretation of the experimental observables, including regarding the behavior of facing phosphate groups in complementary dinucleotides, and their modulation by the sequence. It also provides the opportunity to extensively revisit and refine the coupling between backbone states and inter base pair parameters, which emerges as a common theme across all the complementary dinucleotides. In sum, the global agreement between simulations and experiment reveals new aspects of intrinsic DNA mechanics, a key component of DNA-protein recognition.
Collapse
Affiliation(s)
- Akli Ben Imeddourene
- LBPA, CNRS, ENS Cachan, Université Paris-Saclay, Cachan, France
- Université Pierre et Marie Curie, Paris, France
| | - Ahmad Elbahnsi
- LBPA, CNRS, ENS Cachan, Université Paris-Saclay, Cachan, France
- LPTM, UMR 8089, Université de Cergy-Pontoise, Cergy-Pontoise, France
| | - Marc Guéroult
- UMR S665, INSERM, Université Paris Diderot, INTS, Paris, France
| | - Christophe Oguey
- LPTM, UMR 8089, Université de Cergy-Pontoise, Cergy-Pontoise, France
| | | | - Brigitte Hartmann
- LBPA, CNRS, ENS Cachan, Université Paris-Saclay, Cachan, France
- * E-mail: (NF); (BH)
| |
Collapse
|
15
|
Zgarbová M, Šponer J, Otyepka M, Cheatham TE, Galindo-Murillo R, Jurečka P. Refinement of the Sugar-Phosphate Backbone Torsion Beta for AMBER Force Fields Improves the Description of Z- and B-DNA. J Chem Theory Comput 2015; 11:5723-36. [PMID: 26588601 DOI: 10.1021/acs.jctc.5b00716] [Citation(s) in RCA: 336] [Impact Index Per Article: 37.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Z-DNA duplexes are a particularly complicated test case for current force fields. We performed a set of explicit solvent molecular dynamics (MD) simulations with various AMBER force field parametrizations including our recent refinements of the ε/ζ and glycosidic torsions. None of these force fields described the ZI/ZII and other backbone substates correctly, and all of them underpredicted the population of the important ZI substate. We show that this underprediction can be attributed to an inaccurate potential for the sugar-phosphate backbone torsion angle β. We suggest a refinement of this potential, β(OL1), which was derived using our recently introduced methodology that includes conformation-dependent solvation effects. The new potential significantly increases the stability of the dominant ZI backbone substate and improves the overall description of the Z-DNA backbone. It also has a positive (albeit small) impact on another important DNA form, the antiparallel guanine quadruplex (G-DNA), and improves the description of the canonical B-DNA backbone by increasing the population of BII backbone substates, providing a better agreement with experiment. We recommend using β(OL1) in combination with our previously introduced corrections, εζ(OL1) and χ(OL4), (the combination being named OL15) as a possible alternative to the current β torsion potential for more accurate modeling of nucleic acids.
Collapse
Affiliation(s)
- Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Jiří Šponer
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic.,Institute of Biophysics, Academy of Sciences of the Czech Republic , Královopolská 135, 612 65 Brno, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic
| | - Thomas E Cheatham
- Department of Medicinal Chemistry, University of Utah , 30 South 2000 East, Skaggs 105, Salt Lake City, Utah 84112, United States
| | - Rodrigo Galindo-Murillo
- Department of Medicinal Chemistry, University of Utah , 30 South 2000 East, Skaggs 105, Salt Lake City, Utah 84112, United States
| | - Petr Jurečka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University , 17. listopadu 12, 77146 Olomouc, Czech Republic
| |
Collapse
|
16
|
Zhou Q, Xia X, Luo Z, Liang H, Shakhnovich E. Searching the Sequence Space for Potent Aptamers Using SELEX in Silico. J Chem Theory Comput 2015; 11:5939-46. [PMID: 26642994 DOI: 10.1021/acs.jctc.5b00707] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
To isolate functional nucleic acids that bind to defined targets with high affinity and specificity, which are known as aptamers, the systematic evolution of ligands by exponential enrichment (SELEX) methodology has emerged as the preferred approach. Here, we propose a computational approach, SELEX in silico, that allows the sequence space to be more thoroughly explored regarding binding of a certain target. Our approach consists of two steps: (i) secondary structure-based sequence screening, which aims to collect the sequences that can form a desired RNA motif as an enhanced initial library, followed by (ii) sequence enrichment regarding target binding by molecular dynamics simulation-based virtual screening. Our SELEX in silico method provided a practical computational solution to three key problems in aptamer sequence searching: design of nucleic acid libraries, knowledge of sequence enrichment, and identification of potent aptamers. Six potent theophylline-binding aptamers, which were isolated by SELEX in silico from a sequence space containing 4(13) sequences, were experimentally verified to bind theophylline with high affinity: Kd ranging from 0.16 to 0.52 μM, compared with the dissociation constant of the original aptamer-theophylline, 0.32 μM. These results demonstrate the significant potential of SELEX in silico as a new method for aptamer discovery and optimization.
Collapse
Affiliation(s)
- Qingtong Zhou
- Department of Chemistry and Chemical Biology, Harvard University , Cambridge, Massachusetts 02138, United States
| | - Xiaole Xia
- Department of Chemistry and Chemical Biology, Harvard University , Cambridge, Massachusetts 02138, United States.,Key Laboratory of Industrial Biotechnology, Ministry of Education, School of Biotechnology, Jiangnan University , Wuxi, Jiangsu 214122, People's Republic of China
| | | | | | - Eugene Shakhnovich
- Department of Chemistry and Chemical Biology, Harvard University , Cambridge, Massachusetts 02138, United States
| |
Collapse
|
17
|
Andrushchenko V, Benda L, Páv O, Dračínský M, Bouř P. Vibrational Properties of the Phosphate Group Investigated by Molecular Dynamics and Density Functional Theory. J Phys Chem B 2015; 119:10682-92. [PMID: 26193890 DOI: 10.1021/acs.jpcb.5b05124] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The phosphate group (PO2(-)) is an important building block occurring in many components of living matter including nucleic acids. It provides distinct features in vibrational spectra and is useful as a local probe of NA conformation and interactions with the environment. For this purpose, it is desirable to explore in detail various factors influencing spectral shapes of characteristic phosphate vibrations. In the present study, effects of the solvent and conformational averaging are analyzed for simple model molecules, dimethylphosphate, ethylmethylphosphate, and ethylmethylthiophosphate. Infrared absorption (IR) and Raman spectra were measured and calculated using a combination of molecular dynamics (MD) and density functional theory (DFT). To fully understand the link between the structure and the spectra, the solvent has to be explicitly included in the computational modeling. The results indicate that vibrational properties of the phosphate moiety are very sensitive to its conformation and interactions with the aqueous environment indeed. Polarizable continuum solvent models without explicit water molecules provided significantly worse agreement with the experiment. The combined MD/DFT approach captures well spectral characteristics for the model systems and constitutes the most reliable basis for exploration of phosphate vibrational properties in biomolecular structural studies.
Collapse
Affiliation(s)
- Valery Andrushchenko
- Institute of Organic Chemistry and Biochemistry, Flemingovo náměstí 2, 16610, Prague, Czech Republic
| | - Ladislav Benda
- Institute of Organic Chemistry and Biochemistry, Flemingovo náměstí 2, 16610, Prague, Czech Republic
| | - Ondřej Páv
- Institute of Organic Chemistry and Biochemistry, Flemingovo náměstí 2, 16610, Prague, Czech Republic
| | - Martin Dračínský
- Institute of Organic Chemistry and Biochemistry, Flemingovo náměstí 2, 16610, Prague, Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry, Flemingovo náměstí 2, 16610, Prague, Czech Republic
| |
Collapse
|
18
|
Bruce CD, Ferrara MM, Manka JL, Davis ZS, Register J. Dynamic hydrogen bonding and DNA flexibility in minor groove binders: molecular dynamics simulation of the polyamide f-ImPyIm bound to the Mlu1 (MCB) sequence 5'-ACGCGT-3' in 2:1 motif. J Mol Recognit 2015; 28:325-37. [PMID: 25711379 DOI: 10.1002/jmr.2448] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2014] [Revised: 10/10/2014] [Accepted: 11/09/2014] [Indexed: 12/28/2022]
Abstract
Molecular dynamics simulations of the DNA 10-mer 5'-CCACGCGTGG-3' alone and complexed with the formamido-imidazole-pyrrole-imidazole (f-ImPyIm) polyamide minor groove binder in a 2:1 fashion were conducted for 50 ns using the pbsc0 parameters within the AMBER 12 software package. The change in DNA structure upon binding of f-ImPyIm was evaluated via minor groove width and depth, base pair parameters of Slide, Twist, Roll, Stretch, Stagger, Opening, Propeller, and x-displacement, dihedral angle distributions of ζ, ε, α, and γ determined using the Curves+ software program, and hydrogen bond formation. The dynamic hydrogen bonding between the f-ImPyIm and its cognate DNA sequence was compared to the static image used to predict sequence recognition by polyamide minor groove binders. Many of the predicted hydrogen bonds were present in less than 50% of the simulation; however, persistent hydrogen bonds between G5/15 and the formamido group of f-ImPyIm were observed. It was determined that the DNA is wider in the Complex than without the polyamide binder; however, there is flexibility in this particular sequence, even in the presence of the f-ImPyIm as evidenced by the range of minor groove widths the DNA exhibits and the dynamics of the hydrogen bonding that binds the two f-ImPyIm ions to the minor groove. The Complex consisting of the DNA and the 2 f-ImPyIm binders shows slight fraying of the 5' end of the 10-mer at the end of the simulation, but the portion of the oligomer responsible for recognition and binding is stable throughout the simulation. Several structural changes in the Complex indicate that minor groove binders may have a more active role in inhibiting transcription than just preventing binding of important transcription factors.
Collapse
Affiliation(s)
- Chrystal D Bruce
- Department of Chemistry, John Carroll University, 1 John Carroll Boulevard, University Heights, OH, 44118, USA
| | | | | | | | | |
Collapse
|
19
|
Abstract
Recent applications of solid-state NMR spectroscopy to studies of nucleic acids and their components.
Collapse
Affiliation(s)
- Martin Dračínský
- Institute of Organic Chemistry and Biochemistry
- Prague
- Czech Republic
| | | |
Collapse
|
20
|
Victora A, Möller HM, Exner TE. Accurate ab initio prediction of NMR chemical shifts of nucleic acids and nucleic acids/protein complexes. Nucleic Acids Res 2014; 42:e173. [PMID: 25404135 PMCID: PMC4267612 DOI: 10.1093/nar/gku1006] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
NMR chemical shift predictions based on empirical methods are nowadays indispensable tools during resonance assignment and 3D structure calculation of proteins. However, owing to the very limited statistical data basis, such methods are still in their infancy in the field of nucleic acids, especially when non-canonical structures and nucleic acid complexes are considered. Here, we present an ab initio approach for predicting proton chemical shifts of arbitrary nucleic acid structures based on state-of-the-art fragment-based quantum chemical calculations. We tested our prediction method on a diverse set of nucleic acid structures including double-stranded DNA, hairpins, DNA/protein complexes and chemically-modified DNA. Overall, our quantum chemical calculations yield highly/very accurate predictions with mean absolute deviations of 0.3–0.6 ppm and correlation coefficients (r2) usually above 0.9. This will allow for identifying misassignments and validating 3D structures. Furthermore, our calculations reveal that chemical shifts of protons involved in hydrogen bonding are predicted significantly less accurately. This is in part caused by insufficient inclusion of solvation effects. However, it also points toward shortcomings of current force fields used for structure determination of nucleic acids. Our quantum chemical calculations could therefore provide input for force field optimization.
Collapse
Affiliation(s)
- Andrea Victora
- Department of Chemistry and Zukunftskolleg, Universität Konstanz, 78457 Konstanz, Germany
| | - Heiko M Möller
- Institute of Chemistry, University of Potsdam, Karl-Liebknecht-Strasse 24-25, 14476 Potsdam OT Golm, Germany
| | - Thomas E Exner
- Department of Chemistry and Zukunftskolleg, Universität Konstanz, 78457 Konstanz, Germany Institute of Pharmacy, Eberhard Karls Universität Tübingen, Auf der Morgenstelle 8, 72076 Tübingen, Germany
| |
Collapse
|
21
|
Xu X, Ben Imeddourene A, Zargarian L, Foloppe N, Mauffret O, Hartmann B. NMR studies of DNA support the role of pre-existing minor groove variations in nucleosome indirect readout. Biochemistry 2014; 53:5601-12. [PMID: 25102280 DOI: 10.1021/bi500504y] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
We investigated how the intrinsic sequence-dependent properties probed via the phosphate linkages (BI ↔ BII equilibrium) influence the preferred shape of free DNA, and how this affects the nucleosome formation. First, this exploits NMR solution studies of four B-DNA dodecamers that together cover 39 base pairs of the 5' half of the sequence 601, of special interest for nucleosome formation. The results validate our previous prediction of a systematic, general sequence effect on the intrinsic backbone BII propensities. NMR provides new evidence that the backbone behavior is intimately coupled to the minor groove width. Second, application of the backbone behavior predictions to the full sequence 601 and other relevant sequences demonstrates that alternation of intrinsic low and high BII propensities, coupled to intrinsic narrow and wide minor grooves, largely coincides with the sinusoidal variations of the DNA minor groove width observed in crystallographic structures of the nucleosome. This correspondence is much poorer with low affinity sequences. Overall, the results indicate that nucleosome formation involves an indirect readout process implicating pre-existing DNA minor groove conformations. It also illustrates how the prediction of the intrinsic structural DNA behavior offers a powerful framework to gain explanatory insight on how proteins read DNA.
Collapse
Affiliation(s)
- Xiaoqian Xu
- LBPA, UMR 8113, ENS de Cachan CNRS , 61 avenue du Président Wilson, 94235 Cachan cedex, France
| | | | | | | | | | | |
Collapse
|
22
|
Dračínský M, Möller HM, Exner TE. Conformational Sampling by Ab Initio Molecular Dynamics Simulations Improves NMR Chemical Shift Predictions. J Chem Theory Comput 2013; 9:3806-15. [PMID: 26584127 DOI: 10.1021/ct400282h] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Car-Parrinello molecular dynamics simulations were performed for N-methyl acetamide as a small test system for amide groups in protein backbones, and NMR chemical shifts were calculated based on the generated ensemble. If conformational sampling and explicit solvent molecules are taken into account, excellent agreement between the calculated and experimental chemical shifts is obtained. These results represent a landmark improvement over calculations based on classical molecular dynamics (MD) simulations especially for amide protons, which are predicted too high-field shifted based on the latter ensembles. We were able to show that the better results are caused by the solute-solvents interactions forming shorter hydrogen bonds as well as by the internal degrees of freedom of the solute. Inspired by these results, we propose our approach as a new tool for the validation of force fields due to its power of identifying the structural reasons for discrepancies between the experimental and calculated data.
Collapse
Affiliation(s)
- Martin Dračínský
- Institute of Organic Chemistry and Biochemistry, Academy of Sciences , Flemingovo náměstí 2, 166 10 Prague, Czech Republic.,Department of Chemistry, Durham University , DH1 3LE Durham, United Kingdom
| | - Heiko M Möller
- Department of Chemistry, University of Konstanz , 78457 Konstanz, Germany
| | - Thomas E Exner
- Department of Chemistry, University of Konstanz , 78457 Konstanz, Germany.,Theoretical Medicinal Chemistry and Biophysics, Institute of Pharmacy, Eberhard Karls University Tübingen , Auf der Morgenstelle 8, 72076 Tübingen, Germany
| |
Collapse
|
23
|
Zgarbová M, Luque FJ, Šponer J, Cheatham TE, Otyepka M, Jurečka P. Toward Improved Description of DNA Backbone: Revisiting Epsilon and Zeta Torsion Force Field Parameters. J Chem Theory Comput 2013; 9:2339-2354. [PMID: 24058302 PMCID: PMC3775469 DOI: 10.1021/ct400154j] [Citation(s) in RCA: 221] [Impact Index Per Article: 20.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
We present a refinement of the backbone torsion parameters ε and ζ of the Cornell et al. AMBER force field for DNA simulations. The new parameters, denoted as εζOL1, were derived from quantum-mechanical calculations with inclusion of conformation-dependent solvation effects according to the recently reported methodology (J. Chem. Theory Comput. 2012, 7(9), 2886-2902). The performance of the refined parameters was analyzed by means of extended molecular dynamics (MD) simulations for several representative systems. The results showed that the εζOL1 refinement improves the backbone description of B-DNA double helices and G-DNA stem. In B-DNA simulations, we observed an average increase of the helical twist and narrowing of the major groove, thus achieving better agreement with X-ray and solution NMR data. The balance between populations of BI and BII backbone substates was shifted towards the BII state, in better agreement with ensemble-refined solution experimental results. Furthermore, the refined parameters decreased the backbone RMS deviations in B-DNA MD simulations. In the antiparallel guanine quadruplex (G-DNA) the εζOL1 modification improved the description of non-canonical α/γ backbone substates, which were shown to be coupled to the ε/ζ torsion potential. Thus, the refinement is suggested as a possible alternative to the current ε/ζ torsion potential, which may enable more accurate modeling of nucleic acids. However, long-term testing is recommended before its routine application in DNA simulations.
Collapse
Affiliation(s)
- Marie Zgarbová
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 17 listopadu 12, 771 46, Olomouc, Czech Republic
| | - F. Javier Luque
- Departament de Fisicoquímica and Institut de Biomedicina (IBUB), Facultad de Farmacia, Universidad de Barcelona, Campus Torribera, Prat de la Riba 171, Edifici Verdaguer, 080921 Santa Coloma de Gramanet, Spain
| | - Jiří Šponer
- Institute of Biophysics, Academy of Sciences of the Czech Republic, Královopolská 135, 612 65 Brno, Czech Republic
- CEITEC – Central European Institute of Technology, Campus Bohunice, Kamenice 5, 625 00 Brno, Czech Republic
| | - Thomas E. Cheatham
- Departments of Medicinal Chemistry, University of Utah, 2000 East 30 South Skaggs 105, Salt Lake City, UT 84112, USA
| | - Michal Otyepka
- Regional Centre of Advanced Technologies and Materials, Department of Physical Chemistry, Faculty of Science, Palacky University, 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, 17 listopadu 12, 771 46, Olomouc, Czech Republic
| |
Collapse
|