QM/MM Calculations on Protein-RNA Complexes: Understanding Limitations of Classical MD Simulations and Search for Reliable Cost-Effective QM Methods

Logo poskytovatele

Varování

Publikace nespadá pod Filozofickou fakultu, ale pod Středoevropský technologický institut. Oficiální stránka publikace je na webu muni.cz.
Autoři

POKORNÁ Pavlína KRUSE Holger KREPL Miroslav ŠPONER Jiří

Rok publikování 2018
Druh Článek v odborném periodiku
Časopis / Zdroj Journal of Chemical Theory and Computation
Fakulta / Pracoviště MU

Středoevropský technologický institut

Citace
www https://pubs.acs.org/doi/10.1021/acs.jctc.8b00670
Doi http://dx.doi.org/10.1021/acs.jctc.8b00670
Klíčová slova SMALL NUCLEAR RIBONUCLEOPROTEIN; MOLECULAR-DYNAMICS SIMULATIONS; DENSITY-FUNCTIONAL THEORY; L-HISTIDINE LIGAND; RIBOZYME CATALYSIS; FORCE-FIELD; BACKBONE PARAMETERS; CRYSTAL-STRUCTURE; HAIRPIN RIBOZYME; BINDING DOMAIN
Popis Although atomistic explicit-solvent Molecular Dynamics (MD) is a popular tool to study protein-RNA recognition, satisfactory MD description of protein-RNA complexes is not always achieved. Unfortunately, it is often difficult to separate MD simulation instabilities primarily caused by the simple point-charge molecular mechanics (MM) force fields from problems related to the notorious uncertainties in the starting structures. Herein, we report a series of large-scale QM/MM calculations on the U1A protein-RNA complex. This experimentally well-characterized system has an intricate protein-RNA interface, which is very unstable in MD simulations. The QM/MM calculations identify several H-bonds poorly described by the MM method and thus indicate the sources of instabilities of the U1A interface in MD simulations. The results suggest that advanced QM/MM computations could be used to indirectly rationalize problems seen in MM-based MD simulations of protein-RNA complexes. As the most accurate QM method, we employ the computationally demanding meta-GGA density functional TPSS-D3(BJ)/def2-TZVP level of theory. Because considerably faster methods would be needed to extend sampling and to study even larger protein-RNA interfaces, a set of low-cost QM/MM methods is compared to the TPSS-D3(BJ)/def2-TZVP data. The PBEh-3c and B97-3c density functional composite methods appear to be suitable for protein-RNA interfaces. In contrast, HF-3c and the tight-binding Hamiltonians DFTB3-D3 and GFN-xTB perform unsatisfactorily and do not provide any advantage over the MM description. These conclusions are supported also by similar analysis of a simple HutP protein-RNA interface, which is well-described by MD with the exception of just one H-bond. Some other methodological aspects of QM/MM calculations on protein-RNA interfaces are discussed.
Související projekty:

Používáte starou verzi internetového prohlížeče. Doporučujeme aktualizovat Váš prohlížeč na nejnovější verzi.