Short-Range Imbalances in the AMBER Lennard-Jones Potential for (Deoxy)Ribose center dot center dot center dot Nucleobase Lone-Pair center dot center dot center dot pi Contacts in Nucleic Acids

Kód výsledku v IS VaVaI

<a href="https://www.isvavai.cz/riv?ss=detail&h=RIV%2F68081707%3A_____%2F21%3A00555405" target="_blank" >RIV/68081707:_____/21:00555405 - isvavai.cz</a>

Výsledek na webu

<a href="https://pubs.acs.org/doi/10.1021/acs.jcim.1c01047" target="_blank" >https://pubs.acs.org/doi/10.1021/acs.jcim.1c01047</a>

DOI - Digital Object Identifier

<a href="http://dx.doi.org/10.1021/acs.jcim.1c01047" target="_blank" >10.1021/acs.jcim.1c01047</a>

Jazyk výsledku

angličtina

Název v původním jazyce

Short-Range Imbalances in the AMBER Lennard-Jones Potential for (Deoxy)Ribose center dot center dot center dot Nucleobase Lone-Pair center dot center dot center dot pi Contacts in Nucleic Acids

Popis výsledku v původním jazyce

The lone-pair center dot center dot center dot pi (lp center dot center dot center dot pi) (deoxy)ribose center dot center dot center dot nucleobase stacking is a recurring interaction in Z-DNA and RNAs that is characterized by subvan der Waals lp center dot center dot center dot pi contacts (<3.0 A). It is a part of the structural signature of CpG Z-step motifs in Z-DNA and r(UNCG) tetraloops that are known to behave poorly in molecular dynamics (MD) simulations. Although the exact origin of the MD simulation issues remains unclear, a significant part of the problem might be due to an imbalanced description of nonbonded interactions, including the characteristic lp center dot center dot center dot pi stacking. To gain insights into the links between lp center dot center dot center dot pi stacking and MD, we present an in-depth comparison between accurate large-basis-set double-hybrid Kohn-Sham density functional theory calculations DSD-BLYP-D3/ma-def2-QZVPP (DHDF-D3) and data obtained with the nonbonded potential of the AMBER force field (AFF) for NpN Z-steps (N = G, A, C, and U). Among other differences, we found that the AFF overestimates the DHDF-D3 lp center dot center dot center dot pi distances by similar to 0.1-0.2 A, while the deviation between the DHDF-D3 and AFF descriptions sharply increases in the short-range region of the interaction. Based on atom-in-molecule polarizabilities and symmetry-adapted perturbation theory analysis, we inferred that the DHDF-D3 versus AFF differences partly originate in identical nucleobase carbon atom Lennard-Jones (LJ) parameters despite the presence/absence of connected electron-withdrawing groups that lead to different effective volumes or vdW radii. Thus, to precisely model the very short CpG lp center dot center dot center dot pi contact distances, we recommend revision of the nucleobase atom LJ parameters. Additionally, we suggest that the large discrepancy between DHDF-D3 and AFF short-range repulsive part of the interaction energy potential may significantly contribute to the poor performances of MD simulations of nucleic acid systems containing Z-steps. Understanding where, and if possible why, the point-charge-type effective potentials reach their limits is vital for developing next-generation FFs and for addressing specific issues in contemporary MD simulations.

Název v anglickém jazyce

Short-Range Imbalances in the AMBER Lennard-Jones Potential for (Deoxy)Ribose center dot center dot center dot Nucleobase Lone-Pair center dot center dot center dot pi Contacts in Nucleic Acids

Popis výsledku anglicky

The lone-pair center dot center dot center dot pi (lp center dot center dot center dot pi) (deoxy)ribose center dot center dot center dot nucleobase stacking is a recurring interaction in Z-DNA and RNAs that is characterized by subvan der Waals lp center dot center dot center dot pi contacts (<3.0 A). It is a part of the structural signature of CpG Z-step motifs in Z-DNA and r(UNCG) tetraloops that are known to behave poorly in molecular dynamics (MD) simulations. Although the exact origin of the MD simulation issues remains unclear, a significant part of the problem might be due to an imbalanced description of nonbonded interactions, including the characteristic lp center dot center dot center dot pi stacking. To gain insights into the links between lp center dot center dot center dot pi stacking and MD, we present an in-depth comparison between accurate large-basis-set double-hybrid Kohn-Sham density functional theory calculations DSD-BLYP-D3/ma-def2-QZVPP (DHDF-D3) and data obtained with the nonbonded potential of the AMBER force field (AFF) for NpN Z-steps (N = G, A, C, and U). Among other differences, we found that the AFF overestimates the DHDF-D3 lp center dot center dot center dot pi distances by similar to 0.1-0.2 A, while the deviation between the DHDF-D3 and AFF descriptions sharply increases in the short-range region of the interaction. Based on atom-in-molecule polarizabilities and symmetry-adapted perturbation theory analysis, we inferred that the DHDF-D3 versus AFF differences partly originate in identical nucleobase carbon atom Lennard-Jones (LJ) parameters despite the presence/absence of connected electron-withdrawing groups that lead to different effective volumes or vdW radii. Thus, to precisely model the very short CpG lp center dot center dot center dot pi contact distances, we recommend revision of the nucleobase atom LJ parameters. Additionally, we suggest that the large discrepancy between DHDF-D3 and AFF short-range repulsive part of the interaction energy potential may significantly contribute to the poor performances of MD simulations of nucleic acid systems containing Z-steps. Understanding where, and if possible why, the point-charge-type effective potentials reach their limits is vital for developing next-generation FFs and for addressing specific issues in contemporary MD simulations.

Druh

J_imp - Článek v periodiku v databázi Web of Science

CEP obor

—

OECD FORD obor

30104 - Pharmacology and pharmacy

Projekt

Výsledek vznikl pri realizaci vícero projektů. Více informací v záložce Projekty.

Návaznosti

P - Projekt vyzkumu a vyvoje financovany z verejnych zdroju (s odkazem do CEP)

Rok uplatnění

2021

Kód důvěrnosti údajů

S - Úplné a pravdivé údaje o projektu nepodléhají ochraně podle zvláštních právních předpisů

Název periodika

Journal of Chemical Information and Modeling

ISSN

1549-9596

e-ISSN

1549-960X

Svazek periodika

61

Číslo periodika v rámci svazku

11

Stát vydavatele periodika

US - Spojené státy americké

Počet stran výsledku

14

Strana od-do

5644-5657

Kód UT WoS článku

000757001900031

EID výsledku v databázi Scopus

—