Energies and 2 '-Hydroxyl Group Orientations of RNA Backbone
Conformations. Benchmark CCSD(T)/CBS Database, Electronic
Analysis, and Assessment of DFT Methods and MD Simulations

J 2014

Energies and 2 '-Hydroxyl Group Orientations of RNA Backbone Conformations. Benchmark CCSD(T)/CBS Database, Electronic Analysis, and Assessment of DFT Methods and MD Simulations

MLÁDEK, Arnošt, Pavel BANÁŠ, Petr JUREČKA, Michal OTYEPKA, Marie ZGARBOVÁ et. al.

Základní údaje

Originální název

Energies and 2 '-Hydroxyl Group Orientations of RNA Backbone Conformations. Benchmark CCSD(T)/CBS Database, Electronic Analysis, and Assessment of DFT Methods and MD Simulations

Autoři

MLÁDEK, Arnošt (203 Česká republika), Pavel BANÁŠ (203 Česká republika), Petr JUREČKA (203 Česká republika), Michal OTYEPKA (203 Česká republika), Marie ZGARBOVÁ (203 Česká republika) a Jiří ŠPONER (203 Česká republika, garant, domácí)

Vydání

Journal of Chemical Theory and Computation, Washington DC, American Chemical Society, 2014, 1549-9618

Další údaje

Jazyk

angličtina

Typ výsledku

Článek v odborném periodiku

Obor

10403 Physical chemistry

Stát vydavatele

Spojené státy

Utajení

není předmětem státního či obchodního tajemství

Odkazy

URL

Impakt faktor

Impact factor: 5.498

Kód RIV

RIV/00216224:14740/14:00075659

Organizační jednotka

Středoevropský technologický institut

DOI

http://dx.doi.org/10.1021/ct400837p

UT WoS

000330142400044

Klíčová slova anglicky

DENSITY-FUNCTIONAL THEORY; SUGAR-PHOSPHATE BACKBONE; QUANTUM-CHEMICAL CALCULATIONS; WATSON-CRICK/SUGAR-EDGE; GAUSSIAN-BASIS SETS; ACID BASE-PAIRS; GENERALIZED GRADIENT APPROXIMATION; CORRELATED MOLECULAR CALCULATIONS; NONEMPIRICAL AB-INITIO

Štítky

kontrola MP, MP, rivok

Příznaky

Mezinárodní význam, Recenzováno

Změněno: 10. 6. 2014 11:25, Martina Prášilová

Anotace

V originále

Sugar-phosphate backbone is an electronically complex molecular segment imparting RNA molecules high flexibility and architectonic heterogeneity necessary for their biological functions. The structural variability of RNA molecules is amplified by the presence of the 2'-hydroxyl group, capable of forming multitude of intra- and intermolecular interactions. Bioinformatics studies based on X-ray structure database revealed that RNA backbone samples at least 46 substates known as rotameric families. The present study provides a comprehensive analysis of RNA backbone conformational preferences and 2'-hydroxyl group orientations. First, we create a benchmark database of estimated CCSD(T)/CBS relative energies of all rotameric families and test performance of dispersion-corrected DFT-D3 methods and molecular mechanics in vacuum and in continuum solvent. The performance of the DFT-D3 methods is in general quite satisfactory. The B-LYP-D3 method provides the best trade-off between accuracy and computational demands. B3-LYP-D3 slightly outperforms the new PW6B95-D3 and MPW1B95-D3 and is the second most accurate density functional of the study. The best agreement with CCSD(T)/CBS is provided by DSD-B-LYP-D3 double-hybrid functional, although its large-scale applications may be limited by high computational costs. Molecular mechanics does not reproduce the fine energy differences between the RNA backbone substates. We also demonstrate that the differences in the magnitude of the hyperconjugation effect do not correlate with the energy ranking of the backbone conformations. Further, we investigated the 2'-hydroxyl group orientation preferences. For all families, we conducted a QM and MM hydroxyl group rigid scan in gas phase and solvent. We then carried out set of explicit solvent MD simulations of folded RNAs and analyze 2'-hydroxyl group orientations of different backbone families in MD. The solvent energy profiles determined primarily by the sugar pucker match well with the distribution data derived from the simulations. The QM and MM energy profiles predict the same 2'-hydroxyl group orientation preferences. Finally, we demonstrate that the high energy of unfavorable and rarely sampled 2'-hydroxyl group orientations can be attributed to clashes between occupied orbitals.

Návaznosti

ED1.1.00/02.0068, projekt VaV

Název: CEITEC - central european institute of technology

Citovat

MLÁDEK, Arnošt, Pavel BANÁŠ, Petr JUREČKA, Michal OTYEPKA, Marie ZGARBOVÁ a Jiří ŠPONER. Energies and 2 '-Hydroxyl Group Orientations of RNA Backbone Conformations. Benchmark CCSD(T)/CBS Database, Electronic Analysis, and Assessment of DFT Methods and MD Simulations. Journal of Chemical Theory and Computation. Washington DC: American Chemical Society, 2014, roč. 10, č. 1, s. 463-480. ISSN 1549-9618. Dostupné z: https://dx.doi.org/10.1021/ct400837p.

@article{1185382,
   author = {Mládek, Arnošt and Banáš, Pavel and Jurečka, Petr and Otyepka, Michal and Zgarbová, Marie and Šponer, Jiří},
   article_location = {Washington DC},
   article_number = {1},
   doi = {http://dx.doi.org/10.1021/ct400837p},
   keywords = {DENSITY-FUNCTIONAL THEORY; SUGAR-PHOSPHATE BACKBONE; QUANTUM-CHEMICAL CALCULATIONS; WATSON-CRICK/SUGAR-EDGE; GAUSSIAN-BASIS SETS; ACID BASE-PAIRS; GENERALIZED GRADIENT APPROXIMATION; CORRELATED MOLECULAR CALCULATIONS; NONEMPIRICAL AB-INITIO},
   language = {eng},
   issn = {1549-9618},
   journal = {Journal of Chemical Theory and Computation},
   title = {Energies and 2 '-Hydroxyl Group Orientations of RNA Backbone Conformations. Benchmark CCSD(T)/CBS Database, Electronic Analysis, and Assessment of DFT Methods and MD Simulations},
   url = {http://pubs.acs.org/doi/pdfplus/10.1021/ct400837p},
   volume = {10},
   year = {2014}
}

TY  - JOUR
ID  - 1185382
AU  - Mládek, Arnošt - Banáš, Pavel - Jurečka, Petr - Otyepka, Michal - Zgarbová, Marie - Šponer, Jiří
PY  - 2014
TI  - Energies and 2 '-Hydroxyl Group Orientations of RNA Backbone Conformations. Benchmark CCSD(T)/CBS Database, Electronic Analysis, and Assessment of DFT Methods and MD Simulations
JF  - Journal of Chemical Theory and Computation
VL  - 10
IS  - 1
SP  - 463-480
EP  - 463-480
PB  - American Chemical Society
SN  - 15499618
KW  - DENSITY-FUNCTIONAL THEORY
KW  - SUGAR-PHOSPHATE BACKBONE
KW  - QUANTUM-CHEMICAL CALCULATIONS
KW  - WATSON-CRICK/SUGAR-EDGE
KW  - GAUSSIAN-BASIS SETS
KW  - ACID BASE-PAIRS
KW  - GENERALIZED GRADIENT APPROXIMATION
KW  - CORRELATED MOLECULAR CALCULATIONS
KW  - NONEMPIRICAL AB-INITIO
UR  - http://pubs.acs.org/doi/pdfplus/10.1021/ct400837p
L2  - http://pubs.acs.org/doi/pdfplus/10.1021/ct400837p
N2  - Sugar-phosphate backbone is an electronically complex molecular segment imparting RNA molecules high flexibility and architectonic heterogeneity necessary for their biological functions. The structural variability of RNA molecules is amplified by the presence of the 2'-hydroxyl group, capable of forming multitude of intra- and intermolecular interactions. Bioinformatics studies based on X-ray structure database revealed that RNA backbone samples at least 46 substates known as rotameric families. The present study provides a comprehensive analysis of RNA backbone conformational preferences and 2'-hydroxyl group orientations. First, we create a benchmark database of estimated CCSD(T)/CBS relative energies of all rotameric families and test performance of dispersion-corrected DFT-D3 methods and molecular mechanics in vacuum and in continuum solvent. The performance of the DFT-D3 methods is in general quite satisfactory. The B-LYP-D3 method provides the best trade-off between accuracy and computational demands. B3-LYP-D3 slightly outperforms the new PW6B95-D3 and MPW1B95-D3 and is the second most accurate density functional of the study. The best agreement with CCSD(T)/CBS is provided by DSD-B-LYP-D3 double-hybrid functional, although its large-scale applications may be limited by high computational costs. Molecular mechanics does not reproduce the fine energy differences between the RNA backbone substates. We also demonstrate that the differences in the magnitude of the hyperconjugation effect do not correlate with the energy ranking of the backbone conformations. Further, we investigated the 2'-hydroxyl group orientation preferences. For all families, we conducted a QM and MM hydroxyl group rigid scan in gas phase and solvent. We then carried out set of explicit solvent MD simulations of folded RNAs and analyze 2'-hydroxyl group orientations of different backbone families in MD. The solvent energy profiles determined primarily by the sugar pucker match well with the distribution data derived from the simulations. The QM and MM energy profiles predict the same 2'-hydroxyl group orientation preferences. Finally, we demonstrate that the high energy of unfavorable and rarely sampled 2'-hydroxyl group orientations can be attributed to clashes between occupied orbitals.
ER  -

MLÁDEK, Arnošt, Pavel BANÁŠ, Petr JUREČKA, Michal OTYEPKA, Marie ZGARBOVÁ a Jiří ŠPONER. Energies and 2 '-Hydroxyl Group Orientations of RNA Backbone Conformations. Benchmark CCSD(T)/CBS Database, Electronic Analysis, and Assessment of DFT Methods and MD Simulations. \textit{Journal of Chemical Theory and Computation}. Washington DC: American Chemical Society, 2014, roč.~10, č.~1, s.~463-480. ISSN~1549-9618. Dostupné z: https://dx.doi.org/10.1021/ct400837p.

Podrobný výpis o publikaci