Protein-Nucleic Acid Interactions for RNA Polymerase II Elongation Factors by Molecular Dynamics Simulations

On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules

RNA molecules play a vital role in biological processes within the cell, with significant implications for science and medicine. Notably, the biological functions exerted by specific RNA molecules are often linked to the RNA conformational ensemble. However, the experimental characterization of such three-dimensional RNA structures is challenged by the structural heterogeneity of RNA and by its multiple dynamic interactions with binding partners such as small molecules, proteins, and metal ions. Consequently, our current understanding of the structure-function relationship of RNA molecules is still limited. In this context, we highlight molecular dynamics (MD) simulations as a powerful tool to complement experimental efforts on RNAs. Despite the recognized limitations of current force fields for RNA MD simulations, examining the dynamics of selected RNAs has provided valuable functional insights into their structures.

Insights into the DNA and RNA Interactions of Human Topoisomerase III Beta Using Molecular Dynamics Simulations

Human topoisomerase III beta (hTOP3B) is the only topoisomerase in the human cell that can act on both DNA and RNA substrates. Recent findings have emphasized the physiological importance of hTOP3B and consolidated it as a valuable drug target for antiviral and anticancer therapeutics. Although type IA topoisomerases of different organisms have been studied over the years, the step-by-step interaction of hTOP3B and nucleic acid substrates is still not well understood. Due to the lack of hTOP3B-RNA structures as well as DNA/RNA covalent complexes, computational investigations have been limited. In our study, we utilized molecular dynamics (MD) simulations to study the interactions between hTOP3B and nucleic acids to get a closer look into the residues that play a role in binding DNA or RNA and facilitate catalysis, along with the differences and similarities when hTOP3B interacts with DNA compared to RNA. For this, we generated multiple models of hTOP3B complexed with DNA and RNA sequences using the hTOP3B crystal structure and 8-mer single-stranded DNA and RNA sequences. These models include both covalent and noncovalent complexes, which are then subjected to MD simulations and analyzed. Our findings highlight the complexes' stability, sequence preference, and interactions of the binding pocket residues with different nucleotides. Our work demonstrates that hTOP3B forms stable complexes with both DNA and RNA and provides a better understanding of the enzyme's interaction with different nucleic acid substrate sequences.

Structural impacts of two disease-linked ADAR1 mutants: a molecular dynamics study

Adenosine deaminases acting on RNA (ADARs) are pivotal RNA-editing enzymes responsible for converting adenosine to inosine within double-stranded RNA (dsRNA). Dysregulation of ADAR1 editing activity, often arising from genetic mutations, has been linked to elevated interferon levels and the onset of autoinflammatory diseases. However, understanding the molecular underpinnings of this dysregulation is impeded by the lack of an experimentally determined structure for the ADAR1 deaminase domain. In this computational study, we utilized homology modeling and the AlphaFold2 to construct structural models of the ADAR1 deaminase domain in wild-type and two pathogenic variants, R892H and Y1112F, to decipher the structural impact on the reduced deaminase activity. Our findings illuminate the critical role of structural complementarity between the ADAR1 deaminase domain and dsRNA in enzyme-substrate recognition. That is, the relative position of E1008 and K1120 must be maintained so that they can insert into the minor and major grooves of the substrate dsRNA, respectively, facilitating the flipping-out of adenosine to be accommodated within a cavity surrounding E912. Both amino acid replacements studied, R892H at the orthosteric site and Y1112F at the allosteric site, alter K1120 position and ultimately hinder substrate RNA binding.

Binding of small molecule inhibitors to RNA polymerase-Spt5 complex impacts RNA and DNA stability

Spt5 is an elongation factor that associates with RNA polymerase II (Pol II) during transcription and has important functions in promoter-proximal pausing and elongation processivity. Spt5 was also recognized for its roles in the transcription of expanded-repeat genes that are related to neurodegenerative diseases. Recently, a set of Spt5-Pol II small molecule inhibitors (SPIs) were reported, which selectively inhibit mutant huntingtin gene transcription. Inhibition mechanisms as well as interaction sites of these SPIs with Pol II and Spt5 are not entirely known. In this study, we predicted the binding sites of three selected SPIs at the Pol II-Spt5 interface by docking and molecular dynamics simulations. Two molecules out of three demonstrated strong binding with Spt5 and Pol II, while the other molecule was more loosely bound and sampled multiple binding sites. Strongly bound SPIs indirectly affected RNA and DNA dynamics at the exit site as DNA became more flexible while RNA was stabilized by increased interactions with Spt5. Our results suggest that the transcription inhibition mechanism induced by SPIs can be related to Spt5-nucleic acid interactions, which were altered to some extent with strong binding of SPIs.

RNA-protein complexes and force field polarizability

Molecular dynamic (MD) simulations offer a way to study biomolecular interactions and their dynamics at the atomistic level. There are only a few studies of RNA-protein complexes in MD simulations, and here we wanted to study how force fields differ when simulating RNA-protein complexes: 1) argonaute 2 with bound guide RNA and a target RNA, 2) CasPhi-2 bound to CRISPR RNA and 3) Retinoic acid-inducible gene I C268F variant in complex with double-stranded RNA. We tested three non-polarizable force fields: Amber protein force fields ff14SB and ff19SB with RNA force field OL3, and the all-atom OPLS4 force field. Due to the highly charged and polar nature of RNA, we also tested the polarizable AMOEBA force field and the ff19SB and OL3 force fields with a polarizable water model O3P. Our results show that the non-polarizable force fields lead to compact and stable complexes. The polarizability in the force field or in the water model allows significantly more movement from the complex, but in some cases, this results in the disintegration of the complex structure, especially if the protein contains longer loop regions. Thus, one should be cautious when running long-scale simulations with polarizability. As a conclusion, all the tested force fields can be used to simulate RNA-protein complexes and the choice of the optimal force field depends on the studied system and research question.

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

C-B-A Test of DNA Force Fields

The DNA duplex may be locally strongly bent in complexes with proteins, for example, with polymerases or in a nucleosome. At such bends, the DNA helix is locally in the noncanonical forms A (with a narrow major groove and a large amount of north sugars) or C (with a narrow minor groove and a large share of BII phosphates). To model the formation of such complexes by molecular dynamics methods, the force field is required to reproduce these conformational transitions for a naked DNA. We analyzed the available experimental data on the B-C and B-A transitions under the conditions easily implemented in modeling: in an aqueous NaCl solution. We selected six DNA duplexes which conformations at different salt concentrations are known reliably enough. At low salt concentrations, poly(GC) and poly(A) are in the B-form, classical and slightly shifted to the A-form, respectively. The duplexes ATAT and GGTATACC have a strong and salt concentration dependent bias toward the A-form. The polymers poly(AC) and poly(G) take the C- and A-forms, respectively, at high salt concentrations. The reproduction of the behavior of these oligomers can serve as a test for the balance of interactions between the base stacking and the conformational flexibility of the sugar-phosphate backbone in a DNA force field. We tested the AMBER bsc1 and CHARMM36 force fields and their hybrids, and we failed to reproduce the experiment. In all the force fields, the salt concentration dependence is very weak. The known B-philicity of the AMBER force field proved to result from the B-philicity of its excessively strong base stacking. In the CHARMM force field, the B-form is a result of a fragile balance between the A-philic base stacking (especially for G:C pairs) and the C-philic backbone. Finally, we analyzed some recent simulations of the LacI-, SOX-4-, and Sac7d-DNA complex formation in the framework of the AMBER force field.

Molecular Dynamics Simulations of Protein RNA Complexes by Using an Advanced Electrostatic Model

Protein-RNA interactions are integral to the biological functions of RNA. It is well recognized that molecular dynamics (MD) simulations of protein-RNA complexes are more challenging than those of each component. The difficulty arises from the strong electrostatic interactions and the delicate balance between various types of physical forces at the interface. Previously, MD simulations of protein-RNA complexes have predominantly employed fixed-charge force fields. Although force field modifications have been developed to address problems identified in the simulations, some protein-RNA structures are still hard to reproduce by simulations. Here, we present MD simulations of two representative protein-RNA complexes using the AMOEBA polarizable force field. The van der Waals parameters were refined to reproduce accurate quantum-mechanical data of base-base and base-amino acid interactions. It was found that the refined parameters produced a more stable hydrogen-bond network in the interface. One of the complexes remained stable during the short simulations, whereas it could quickly break down in previous simulations using fixed-charge force fields. There was reversible breaking and formation of hydrogen bonds that are observed in the crystal structure, which may indicate the difference in solution and crystal structures. While further improvement and validation of the force fields are still needed, this work demonstrates that polarizable force fields are promising for the study of protein-RNA complexes.