Despite significant advances in the performance and reliability of biomolecular simulation approaches, non- helical nucleic acid structures are proving difficult to fully and accurately model with currently available biomolecular force fields.
This research aims to assess, validate, and improve force fields for nucleic acids, and also to better understand the true conformational ensemble of model DNA and RNA systems by best fits to NMR experiment using Maximum Entropy methods. This includes characterization of not only the dominant conformations, but also excited states or low population states. Model systems include RNA dinucleotides, RNA tetranucleotides, RNA tetraloops, DNA mini-dumbbells, and NMR-derived NMR structures that are known to populate multiple structures or excited states. Methods employed include state-of-the-art multi-dimensional replica-exchange molecular dynamics (M-REMD) simulations and various NMR approaches to provide more insight into new tetranucleotides both in structure (NOE, J-coupling, etc.) and for excited states via NMR relaxation and NMR relaxation dispersion experiments. The M-REMD code will be extended to allow asynchronous ensemble instances (for greater efficiency and queue backfill) and adaptivity and steering with on-the-fly analyses to assess convergence and where more/fewer ensemble instances are needed. Various approaches to force field improvement include surrogate methods for re-weighting converged MD trajectory ensembles and project free energy surfaces with Multistate Bennett Acceptance Ratio methods as a function of force field parameter change, where possible, via parameter scanning, via fits to high level base-base interaction energies from high level QM, and small model compound fits to liquid densities and other methods. In addition, the group will work with the Open Force Field Initiative to adapt, assess, and validate their small molecule force fields to nucleic acids. Using M-REMD methods the team has proven the ability to fully sample the conformational ensemble of the proposed model systems with large-scale computation on GPUs with multiple different force fields. The group has considerable experience in large-scale simulation of DNA and RNA and a proven track record of collaboration and dissemination of research findings and results. Development of better methods for simulation of the structure, dynamics and interactions of nucleic acids provides the means to approach drug- ability and the design of novel therapeutics, and also to provide greater insight into the role of structure, dynamics and conformational change in function which at a basic science level provides unique capabilities that could have considerable health relevance if the methods are made to function correctly and accurately.

Public Health Relevance

By assessing, validating and improving biomolecular simulations approaches and their accuracy towards modeling nucleic acids, we can gain greater insight into nucleic acid structure, dynamics, function, and interactions with other biomolecules. Ultimately such approaches could be applied to better understand how to modulate function, understand biological processes mediated by nucleic acids, and provide for better capabilities for the design of drugs to target nucleic acids. For example, such approaches could provide new insight into anti- viral therapeutics, riboswitches, and greater insight into health related processes involving nucleic acids at the atomistic level.

Agency
National Institute of Health (NIH)
Institute
National Institute of General Medical Sciences (NIGMS)
Type
Research Project (R01)
Project #
2R01GM081411-06
Application #
9887485
Study Section
Macromolecular Structure and Function D Study Section (MSFD)
Program Officer
Lyster, Peter
Project Start
2008-02-01
Project End
2023-08-31
Budget Start
2019-09-23
Budget End
2020-08-31
Support Year
6
Fiscal Year
2019
Total Cost
Indirect Cost
Name
University of Utah
Department
Pharmacology
Type
Schools of Pharmacy
DUNS #
009095365
City
Salt Lake City
State
UT
Country
United States
Zip Code
84112
Galindo-Murillo, Rodrigo; Davis, Darrell R; Cheatham 3rd, Thomas E (2016) Probing the influence of hypermodified residues within the tRNA3(Lys) anticodon stem loop interacting with the A-loop primer sequence from HIV-1. Biochim Biophys Acta 1860:607-17
Waters, James T; Lu, Xiang-Jun; Galindo-Murillo, Rodrigo et al. (2016) Transitions of Double-Stranded DNA Between the A- and B-Forms. J Phys Chem B 120:8449-56
Galindo-Murillo, Rodrigo; Roe, Daniel R; Cheatham 3rd, Thomas E (2015) Convergence and reproducibility in molecular dynamics simulations of the DNA duplex d(GCACGAACGAACGAACGC). Biochim Biophys Acta 1850:1041-1058
Galindo-Murillo, Rodrigo; GarcĂ­a-Ramos, Juan Carlos; Ruiz-Azuara, Lena et al. (2015) Intercalation processes of copper complexes in DNA. Nucleic Acids Res 43:5364-76
Henriksen, Niel M; Hayatshahi, Hamed S; Davis, Darrell R et al. (2014) Structural and energetic analysis of 2-aminobenzimidazole inhibitors in complex with the hepatitis C virus IRES RNA using molecular dynamics simulations. J Chem Inf Model 54:1758-72
Pasi, Marco; Maddocks, John H; Beveridge, David et al. (2014) ?ABC: a systematic microsecond molecular dynamics study of tetranucleotide sequence effects in B-DNA. Nucleic Acids Res 42:12272-83
Roe, Daniel R; Bergonzo, Christina; Cheatham 3rd, Thomas E (2014) Evaluation of enhanced sampling provided by accelerated molecular dynamics with Hamiltonian replica exchange methods. J Phys Chem B 118:3543-52
Thibault, Julien C; Cheatham 3rd, Thomas E; Facelli, Julio C (2014) iBIOMES Lite: summarizing biomolecular simulation data in limited settings. J Chem Inf Model 54:1810-9
Galindo-Murillo, Rodrigo; Roe, Daniel R; Cheatham 3rd, Thomas E (2014) On the absence of intrahelical DNA dynamics on the ?s to ms timescale. Nat Commun 5:5152
Bergonzo, Christina; Galindo-Murillo, Rodrigo; Cheatham 3rd, Thomas E (2014) Molecular modeling of nucleic Acid structure: electrostatics and solvation. Curr Protoc Nucleic Acid Chem 55:7.9.1-27

Showing the most recent 10 out of 29 publications