Mechanism of Peroxynitrite Interaction with Ferric M. tuberculosis Nitrobindin: A Computational Study.

Nitrobindins (Nbs) are all-ß-barrel heme proteins present along the evolutionary ladder. They display a highly solvent-exposed ferric heme group with the iron atom being coordinated by the proximal His residue and a water molecule at the distal position. Ferric nitrobindins (Nb(III)) play a role in the conversion of toxic peroxynitrite (ONOO-) to harmless nitrate, with the value of the second-order rate constant being similar to those of most heme proteins. The value of the second-order rate constant of Nbs increases as the pH decreases; this suggests that Nb(III) preferentially reacts with peroxynitrous acid (ONOOH), although ONOO- is more nucleophilic. In this work, we shed light on the molecular basis of the ONOO- and ONOOH reactivity of ferric Mycobacterium tuberculosis Nb (Mt-Nb(III)) by dissecting the ligand migration toward the active site, the water molecule release, and the ligand binding process by computer simulations. Classical molecular dynamics simulations were performed by employing a steered molecular dynamics approach and the Jarzynski equality to obtain ligand migration free energy profiles for both ONOO- and ONOOH. Our results indicate that ONOO- and ONOOH migration is almost unhindered, consistent with the exposed metal center of Mt-Nb(III). To further analyze the ligand binding process, we computed potential energy profiles for the displacement of the Fe(III)-coordinated water molecule using a hybrid QM/MM scheme at the DFT level and a nudged elastic band approach. These results indicate that ONOO- exhibits a much larger barrier for ligand displacement than ONOOH, suggesting that water displacement is assisted by protonation of the leaving group by the incoming ONOOH.

Binding mechanism of disulfide species to ferric hemeproteins: The case of metmyoglobin.

The interactions of the heme iron of hemeproteins with sulfide and disulfide compounds are of potential interest as physiological signaling processes. While the interaction with hydrogen sulfide has been described computationally and experimentally, the reaction with disulfide, and specifically the molecular mechanism for ligand binding has not been studied in detail. In this work, we study the association process for disulfane and its conjugate base disulfanide at different pH conditions. Additionally, by means of advanced sampling techniques based on multiple steered molecular dynamics, we provide free energy profiles for ligand migration for both acid/base species, showing a similar behavior to the previously reported for the related H2S/HS¯ pair. Finally, we studied the ligand interchange reaction (H2O/H2S, HS¯ and H2O/HSSH, HSS¯) by means of hybrid quantum mechanics-molecular mechanics calculations. We show that the anionic species are able to displace more efficiently the H2O bound to the iron, and that the H-bond network in the distal cavity can help the neutral species to perform the reaction. Altogether, we provide a molecular explanation for the experimental information and show that the global association process depends on a fine balance between the migration towards the active site and the ligand interchange reaction.

Direct Cysteine Desulfurase Activity Determination by NMR and the Study of the Functional Role of Key Structural Elements of Human NFS1.

The mitochondrial cysteine desulfurase NFS1 is an essential PLP-dependent enzyme involved in iron-sulfur cluster assembly. The enzyme catalyzes the desulfurization of the l-Cys substrate, producing a persulfide and l-Ala as products. In this study, we set the measurement of the product l-Ala by NMR in vitro by means of 1H NMR spectra acquisition. This methodology provided us with the possibility of monitoring the reaction in both fixed-time and real-time experiments, with high sensitivity and accuracy. By studying I452A, W454A, Q456A, and H457A NFS1 variants, we found that the C-terminal stretch (CTS) of the enzyme is critical for function. Specifically, mutation of the extremely conserved position W454 resulted in highly decreased activity. Additionally, we worked on two singular variants: "GGG" and C158A. In the former, the catalytic Cys-loop was altered by including two Gly residues to increase the flexibility of this loop. This variant had significantly impaired activity, indicating that the Cys-loop motions are fine-tuned in the wild-type enzyme. In turn, for C158A, we found an unanticipated increase in l-Cys desulfurase activity. Furthermore, we carried out molecular dynamics simulations of the supercomplex dedicated to iron-sulfur cluster biosynthesis, which includes NFS1, ACP, ISD11, ISCU2, and FXN subunits. We identified CTS as a key element that established interactions with ISCU2 and FXN concurrently; we found specific interactions that are established when FXN is present, reinforcing the idea that FXN not only forms part of the iron-sulfur cluster assembly site but also modulates the internal motions of ISCU2.

Best Practices on QM/MM Simulations of Biological Systems.

During the second half of the 20th century, following structural biology hallmark works on DNA and proteins, biochemists shifted their questions from "what does this molecule look like?" to "how does this process work?". Prompted by the theoretical and practical developments in computational chemistry, this led to the emergence of biomolecular simulations and, along with the 2013 Nobel Prize in Chemistry, to the development of hybrid QM/MM methods. QM/MM methods are necessary whenever the problem we want to address involves chemical reactivity and/or a change in the system's electronic structure, with archetypal examples being the studies of an enzyme's reaction mechanism and a metalloprotein's active site. In the last decades QM/MM methods have seen an increasing adoption driven by their incorporation in widely used biomolecular simulation software. However, properly setting up a QM/MM simulation is not an easy task, and several issues need to be properly addressed to obtain meaningful results. In the present work, we describe both the theoretical concepts and practical issues that need to be considered when performing QM/MM simulations. We start with a brief historical perspective on the development of these methods and describe when and why QM/MM methods are mandatory. Then we show how to properly select and analyze the performance of the QM level of theory, the QM system size, and the position and type of the boundaries. We show the relevance of performing prior QM model system (or QM cluster) calculations in a vacuum and how to use the corresponding results to adequately calibrate those derived from QM/MM. We also discuss how to prepare the starting structure and how to select an adequate simulation strategy, including those based on geometry optimizations as well as free energy methods. In particular, we focus on the determination of free energy profiles using multiple steered molecular dynamics (MSMD) combined with Jarzynski's equation. Finally, we describe the results for two illustrative and complementary examples: the reaction performed by chorismate mutase and the study of ligand binding to hemoglobins. Overall, we provide many practical recommendations (or shortcuts) together with important conceptualizations that we hope will encourage more and more researchers to incorporate QM/MM studies into their research projects.

Class III Peroxidases PRX01, PRX44, and PRX73 Control Root Hair Growth in Arabidopsis thaliana.

Root hair cells are important sensors of soil conditions. They grow towards and absorb water-soluble nutrients. This fast and oscillatory growth is mediated by continuous remodeling of the cell wall. Root hair cell walls contain polysaccharides and hydroxyproline-rich glycoproteins, including extensins (EXTs). Class-III peroxidases (PRXs) are secreted into the apoplastic space and are thought to trigger either cell wall loosening or polymerization of cell wall components, such as Tyr-mediated assembly of EXT networks (EXT-PRXs). The precise role of these EXT-PRXs is unknown. Using genetic, biochemical, and modeling approaches, we identified and characterized three root-hair-specific putative EXT-PRXs, PRX01, PRX44, and PRX73. prx01,44,73 triple mutation and PRX44 and PRX73 overexpression had opposite effects on root hair growth, peroxidase activity, and ROS production, with a clear impact on cell wall thickness. We use an EXT fluorescent reporter with contrasting levels of cell wall insolubilization in prx01,44,73 and PRX44-overexpressing background plants. In this study, we propose that PRX01, PRX44, and PRX73 control EXT-mediated cell wall properties during polar expansion of root hair cells.

Conformational stability, dynamics and function of human frataxin: Tryptophan side chain interplay.

In humans, the loss of frataxin results in Friedreich's Ataxia, a neurodegenerative disease, in which a deficit in the iron-sulfur cluster assembly is observed. In this work, we analyzed three frataxin variants in which one tryptophan was replaced by a glycine: W155G, W168G and W173G. As expected, given its localization in the assembly site, W155G was not able to activate the desulfurase activity of the supercomplex for iron-sulfur cluster assembly. In turn, W168G, which was significantly more unstable than W155G, was fully active. W173G, which was highly unstable as W168G, showed a significantly decreased activity, only slightly higher than W155G. As W168G and W173G were highly sensitive to proteolysis, we investigated the protein motions by molecular dynamic simulations. We observed that W173G may display altered motions at the Trp155 site. Furthermore, we revealed a H-bond network in which Trp155 takes part, involving residues Gln148, Asn151, Gln153 and Arg165. We suggest that this motion modulation that specifically alters the population of different Trp155 rotamers can be directly transferred to the assembly site, altering the dynamics of the ISCU His137 key residue. This hypothesis was also contrasted by means of molecular dynamic simulations of frataxin in the context of the complete supercomplex. We propose that the supercomplex requires very definite motions of Trp155 to consolidate the assembly site.

Control of distal lysine coordination in a monomeric hemoglobin: A role for heme peripheral interactions.

THB1 is a monomeric truncated hemoglobin (TrHb) found in the cytoplasm of the green alga Chlamydomonas reinhardtii. The canonical heme coordination scheme in hemoglobins is a proximal histidine ligand and an open distal site. In THB1, the latter site is occupied by Lys53, which is likely to facilitate Fe(II)/Fe(III) redox cycling but hinders dioxygen binding, two features inherent to the NO dioxygenase activity of the protein. TrHb surveys show that a lysine at a position aligning with Lys53 is an insufficient determinant of coordination, and in this study, we sought to identify factors controlling lysine affinity for the heme iron. We solved the "Lys-off" X-ray structure of THB1, represented by the cyanide adduct of the Fe(III) protein, and hypothesized that interactions that differ between the known "Lys-on" structure and the Lys-off structure participate in the control of Lys53 affinity for the heme iron. We applied an experimental approach (site-directed mutagenesis, heme modification, pH titrations in the Fe(III) and Fe(II) states) and a computational approach (MD simulations in the Fe(II) state) to assess the role of heme propionate-protein interactions, distal helix capping, and the composition of the distal pocket. All THB1 modifications resulted in a weakening of lysine affinity and affected the coupling between Lys53 proton binding and heme redox potential. The results supported the importance of specific heme peripheral interactions for the pH stability of iron coordination and the ability of the protein to undergo redox reactions.

Distal lysine (de)coordination in the algal hemoglobin THB1: A combined computer simulation and experimental study.

THB1 is a monomeric truncated hemoglobin from the green alga Chlamydomonas reinhardtii. In the absence of exogenous ligands and at neutral pH, the heme group of THB1 is coordinated by two protein residues, Lys53 and His77. THB1 is thought to function as a nitric oxide dioxygenase, and the distal binding of O2 requires the cleavage of the Fe-Lys53 bond accompanied by protonation and expulsion of the lysine from the heme cavity into the solvent. Nuclear magnetic resonance spectroscopy and crystallographic data have provided dynamic and structural insights of the process, but the details of the mechanism have not been fully elucidated. We applied a combination of computer simulations and site-directed mutagenesis experiments to shed light on this issue. Molecular dynamics simulations and hybrid quantum mechanics/molecular mechanics restrained optimizations were performed to explore the nature of the transition between the decoordinated and lysine-bound states of the ferrous heme in THB1. Lys49 and Arg52, which form ionic interactions with the heme propionates in the X-ray structure of lysine-bound THB1, were observed to assist in maintaining Lys53 inside the protein cavity and play a key role in the transition. Lys49Ala, Arg52Ala and Lys49Ala/Arg52Ala THB1 variants were prepared, and the consequences of the replacements on the Lys (de)coordination equilibrium were characterized experimentally for comparison with computational prediction. The results reinforced the dynamic role of protein-propionate interactions and strongly suggested that cleavage of the Fe-Lys53 bond and ensuing conformational rearrangement is facilitated by protonation of the amino group inside the distal cavity.

Heme-Thiolate Perturbation in Cystathionine ß-Synthase by Mercury Compounds.

Cystathionine ß-synthase (CBS) is an enzyme involved in sulfur metabolism that catalyzes the pyridoxal phosphate-dependent condensation of homocysteine with serine or cysteine to form cystathionine and water or hydrogen sulfide (H2S), respectively. CBS possesses a b-type heme coordinated by histidine and cysteine. Fe(III)-CBS is inert toward exogenous ligands, while Fe(II)-CBS is reactive. Both Fe(III)- and Fe(II)-CBS are sensitive to mercury compounds. In this study, we describe the kinetics of the reactions with mercuric chloride (HgCl2) and p-chloromercuribenzoic acid. These reactions were multiphasic and resulted in five-coordinate CBS lacking thiolate ligation, with six-coordinate species as intermediates. Computational QM/MM studies supported the feasibility of formation of species in which the thiolate is proximal to both the iron ion and the mercury compound. The reactions of Fe(II)-CBS were faster than those of Fe(III)-CBS. The observed rate constants of the first phase increased hyperbolically with concentration of the mercury compounds, with limiting values of 0.3-0.4 s-1 for Fe(III)-CBS and 40 ± 4 s-1 for Fe(II)-CBS. The data were interpreted in terms of alternative models of conformational selection or induced fit. Exposure of Fe(III)-CBS to HgCl2 led to heme release and activity loss. Our study reveals the complexity of the interactions between mercury compounds and CBS.

Relationship between activity and stability: Design and characterization of stable variants of human frataxin.

The relationships between conformational dynamics, stability and protein function are not obvious. Frataxin (FXN) is an essential protein that forms part of a supercomplex dedicated to the iron-sulfur (Fe-S) cluster assembly within the mitochondrial matrix. In humans, the loss of FXN expression or a decrease in its functionality results in Friedreich's Ataxia, a cardio-neurodegenerative disease. Recently, the way in which FXN interacts with the rest of the subunits of the supercomplex was uncovered. This opens a window to explore relationships between structural dynamics and function. In this study, we prepared a set of FXN variants spanning a broad range of conformational stabilities. Variants S160I, S160M and A204R were more stable than the wild-type and showed similar biological activity. Additionally, we prepared SILCAR, a variant that combines S160I, L203C and A204R mutations. SILCAR was 2.4 kcal mol-1 more stable and equally active. Some of the variants were significantly more resistant to proteolysis than the wild-type FXN. SILCAR showed the highest resistance, suggesting a more rigid structure. It was corroborated by means of molecular dynamics simulations. Relaxation dispersion NMR experiments comparing SILCAR and wild-type variants suggested similar internal motions in the microsecond to millisecond timescale. Instead, variant S157I showed higher denaturation resistance but a significant lower function, similarly to that observed for the FRDA variant N146K. We concluded that the contribution of particular side chains to the conformational stability of FXN might be highly subordinated to their impact on both the protein function and the stability of the functional supercomplex.

Structure of the Human ACP-ISD11 Heterodimer.

In recent years, the mammalian mitochondrial protein complex for iron-sulfur cluster assembly has been the focus of important studies. This is partly because of its high degree of relevance in cell metabolism and because mutations of the involved proteins are the cause of several human diseases. Cysteine desulfurase NFS1 is the key enzyme of the complex. At present, it is well-known that the active form of NFS1 is stabilized by the small protein ISD11. In this work, the structure of the human mitochondrial ACP-ISD11 heterodimer was determined at 2.0 Å resolution. ACP-ISD11 forms a cooperative unit stabilized by several ionic interactions, hydrogen bonds, and apolar interactions. The 4'-phosphopantetheine-acyl chain, which is covalently bound to ACP, interacts with several residues of ISD11, modulating together with ACP the foldability of ISD11. Recombinant human ACP-ISD11 was able to interact with the NFS1 desulfurase, thus yielding an active enzyme, and the NFS1/ACP-ISD11 core complex was activated by frataxin and ISCU proteins. Internal motions of ACP-ISD11 were studied by molecular dynamics simulations, showing the persistence of the interactions between both protein chains. The conformation of the dimer is similar to that found in the context of the (NFS1/ACP-ISD11)2 supercomplex core, which contains the Escherichia coli ACP instead of the human variant. This fact suggests a sequential mechanism for supercomplex consolidation, in which the ACP-ISD11 complex may fold independently and, after that, the NFS1 dimer would be stabilized.

Ligand Binding Rate Constants in Heme Proteins Using Markov State Models and Molecular Dynamics Simulations.

Computer simulation studies of the molecular basis for ligand migration in proteins allow the description of key events such as the transition between docking sites, displacement of existing ligands and solvent molecules, and open/closure of specific "gates", among others. In heme proteins, ligand migration from the solvent to the active site preludes the binding to the heme iron and triggers different functions. In this work, molecular dynamics simulations, a Markov State Model of migration and empirical kinetic equations are combined to study the migration of O2 and NO in two truncated hemoglobins of Mycobacterium tuberculosis (Mt-TrHbN and Mt-TrHbO). For Mt-TrHbN, we show that the difference in the association constant in the oxy and deoxy states relies mainly in the displacement of water molecules anchored in the distal cavity in the deoxy form. The results here provide a valuable approach to study ligand migration in globins.

Thermal Stability of Globins: Implications of Flexibility and Heme Coordination Studied by Molecular Dynamics Simulations.

Proteins are sensitive to temperature, and abrupt changes in the normal temperature conditions can have a profound impact on both structure and function, leading to protein unfolding. However, the adaptation of certain organisms to extreme conditions raises questions about the structural features that permit the structure and function of proteins to be preserved under these adverse conditions. To gain insight into the molecular basis of protein thermostability in the globin family, we have examined three representative examples: human neuroglobin, horse heart myoglobin, and Drosophila hemoglobin, which differ in their melting temperatures and coordination states of the heme iron in the absence of external ligands. In order to elucidate the possible mechanisms that govern the thermostability of these proteins, microsecond-scale classical molecular dynamics simulations were performed at different temperatures. Structural fluctuations and essential dynamics were analyzed, indicating that the flexibility of the CD region, which includes the two short C and D helixes and the connecting CD loop, is directly related to the thermostability. We observed that a larger inherent flexibility of the protein produces higher thermostability, probably concentrating the thermal fluctuations observed at high temperature in flexible regions, preventing unfolding. Globally, the results of this work improve our understanding of thermostability in the globin family.

Filling the Gaps to Solve the Extensin Puzzle.

Extensins (EXTs) are highly repetitive plant O-glycoproteins that require several post-translational modifications (PTMs) to become functional in plant cell walls. First, they are hydroxylated on contiguous proline residues; then they are O-glycosylated on hydroxyproline and serine. After secretion into the apoplast, O-glycosylated EXTs form a tridimensional network organized by inter- and intra-Tyr linkages. Recent studies have made significant progress in the identification of the enzymatic machinery required to process EXTs, which includes prolyl 4-hydroxylases, glycosyltransferases, papain-type cysteine endopeptidases, and peroxidases. EXTs are abundant in plant tissues and are particularly important in rapidly expanding root hairs and pollen tubes, which grow in a polar manner. Small changes in EXT PTMs affect fast-growing cells, although the molecular mechanisms underlying this regulation are unknown. In this review, we highlight recent advances in our understanding of EXT modifications throughout the secretory pathway, EXT assembly in cell walls, and possible sensing mechanisms involving the Catharanthus roseus cell surface sensor receptor-like kinases located at the interface between the apoplast and the cytoplasmic side of the plasma membrane.

Tertiary and quaternary structural basis of oxygen affinity in human hemoglobin as revealed by multiscale simulations.

Human hemoglobin (Hb) is a benchmark protein of structural biology that shaped our view of allosterism over 60 years ago, with the introduction of the MWC model based on Perutz structures of the oxy(R) and deoxy(T) states and the more recent Tertiary Two-State model that proposed the existence of individual subunit states -"r" and "t"-, whose structure is yet unknown. Cooperative oxygen binding is essential for Hb function, and despite decades of research there are still open questions related to how tertiary and quaternary changes regulate oxygen affinity. In the present work, we have determined the free energy profiles of oxygen migration and for HisE7 gate opening, with QM/MM calculations of the oxygen binding energy in order to address the influence of tertiary differences in the control of oxygen affinity. Our results show that in the α subunit the low to high affinity transition is achieved by a proximal effect that mostly affects oxygen dissociation and is the driving force of the allosteric transition, while in the ß subunit the affinity change results from a complex interplay of proximal and distal effects, including an increase in the HE7 gate opening, that as shown by free energy profiles promotes oxygen uptake.

Coarse-Grained Simulations of Heme Proteins: Validation and Study of Large Conformational Transitions.

Heme proteins are ubiquitous in nature and perform many diverse functions in all kingdoms of life. Many of these functions are related to large-scale conformational transitions and allosteric processes. Sampling of these large conformational changes is computationally very challenging. In this context, coarse-grain simulations emerge as an efficient approach to explore the conformational landscape. In this work, we present a coarse-grained model of the heme group and thoroughly validate this model in different benchmark examples, which include the monomeric heme proteins myoglobin and neuroglobin and the tetrameric human hemoglobin where we evaluated the method's ability to explore conformational changes (as the formation of hexacoordinated species) and allosteric transitions (as the well-known R → T transition). The obtained results are compared with atomistic molecular dynamics simulations. Overall, the results indicate that this approach conserves the essential dynamical information on different allosteric processes.

Structural Study of a Flexible Active Site Loop in Human Indoleamine 2,3-Dioxygenase and Its Functional Implications.

Human indoleamine 2,3-dioxygenase catalyzes the oxidative cleavage of tryptophan to N-formyl kynurenine, the initial and rate-limiting step in the kynurenine pathway. Additionally, this enzyme has been identified as a possible target for cancer therapy. A 20-amino acid protein segment (the JK loop), which connects the J and K helices, was not resolved in the reported hIDO crystal structure. Previous studies have shown that this loop undergoes structural rearrangement upon substrate binding. In this work, we apply a combination of replica exchange molecular dynamics simulations and site-directed mutagenesis experiments to characterize the structure and dynamics of this protein region. Our simulations show that the JK loop can be divided into two regions: the first region (JK loop(C)) displays specific and well-defined conformations and is within hydrogen bonding distance of the substrate, while the second region (JK loop(N)) is highly disordered and exposed to the solvent. The peculiar flexible nature of JK loop(N) suggests that it may function as a target for post-translational modifications and/or a mediator for protein-protein interactions. In contrast, hydrogen bonding interactions are observed between the substrate and Thr379 in the highly conserved "GTGG" motif of JK loop(C), thereby anchoring JK loop(C) in a closed conformation, which secures the appropriate substrate binding mode for catalysis. Site-directed mutagenesis experiments confirm the key role of this residue, highlighting the importance of the JK loop(C) conformation in regulating the enzymatic activity. Furthermore, the existence of the partially and totally open conformations in the substrate-free form suggests a role of JK loop(C) in controlling substrate and product dynamics.

CG2AA: backmapping protein coarse-grained structures.

UNLABELLED: Coarse grain (CG) models allow long-scale simulations with a much lower computational cost than that of all-atom simulations. However, the absence of atomistic detail impedes the analysis of specific atomic interactions that are determinant in most interesting biomolecular processes. In order to study these phenomena, it is necessary to reconstruct the atomistic structure from the CG representation. This structure can be analyzed by itself or be used as an onset for atomistic molecular dynamics simulations. In this work, we present a computer program that accurately reconstructs the atomistic structure from a CG model for proteins, using a simple geometrical algorithm. AVAILABILITY AND IMPLEMENTATION: The software is free and available online at http://www.ic.fcen.uba.ar/cg2aa/cg2aa.py SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online. CONTACT: lula@qi.fcen.uba.ar.

Engineered chimeras reveal the structural basis of hexacoordination in globins: a case study of neuroglobin and myoglobin.

BACKGROUND: Myoglobin (Mb) and neuroglobin (Ngb) are representative members of pentacoordinated and bis-histidyl, hexacoordinated globins. In spite of their low sequence identity, they show surprisingly similar three-dimensional folds. The ability of Ngb to form a hexacoordinated bis-histidyl complex with the distal HisE7 has a strong impact on ligand affinity. The factors governing such different behaviors have not been completely understood yet, even though they are extremely relevant to establish structure-function relationships within the globin superfamily. METHODS: In this work we generated chimeric proteins by swapping a previously identified regulatory segment - the CD region - and evaluated comparatively the structural and functional properties of the resulting proteins by molecular dynamics simulations, and spectroscopic and kinetic investigations. RESULTS: Our results show that chimeric proteins display heme coordination properties displaced towards those expected for the corresponding CD region. In particular, in the absence of exogenous ligands, chimeric Mb is found as a partially hexacoordinated bis-histidyl species, whereas chimeric Ngb shows a lower equilibrium constant for forming a hexacoordinated bis-histidyl species. CONCLUSIONS: While these results confirm the regulatory role of the CD region for bis-histidyl hexacoordination, they also suggest that additional sources contribute to fine tune the equilibrium. General significance Globins constitute a ubiquitous group of heme proteins widely found in all kingdoms of life. These findings raise challenging questions regarding the structure-function relationships in these proteins, as bis-histidyl hexacoordination emerges as a novel regulatory mechanism of the physiological function of globins.

Molecular dynamics simulations identify time scale of conformational changes responsible for conformational selection in molecular recognition of HIV-1 transactivation responsive RNA.

The HIV-1 Tat protein and several small molecules bind to HIV-1 transactivation responsive RNA (TAR) by selecting sparsely populated but pre-existing conformations. Thus, a complete characterization of TAR conformational ensemble and dynamics is crucial to understand this paradigmatic system and could facilitate the discovery of new antivirals targeting this essential regulatory element. We show here that molecular dynamics simulations can be effectively used toward this goal by bridging the gap between functionally relevant time scales that are inaccessible to current experimental techniques. Specifically, we have performed several independent microsecond long molecular simulations of TAR based on one of the most advanced force fields available for RNA, the parmbsc0 AMBER. Our simulations are first validated against available experimental data, yielding an excellent agreement with measured residual dipolar couplings and order parameter S(2). This contrast with previous molecular dynamics simulations (Salmon et al., J. Am. Chem. Soc. 2013 135, 5457-5466) based on the CHARMM36 force field, which could achieve only modest accord with the experimental RDC values. Next, we direct the computation toward characterizing the internal dynamics of TAR over the microsecond time scale. We show that the conformational fluctuations observed over this previously elusive time scale have a strong functionally oriented character in that they are primed to sustain and assist ligand binding.