Search

The development of reliable and extensible molecular mechanics (MM) force fields-fast, empirical models characterizing the potential energy surface of molecular systems-is indispensable for biomolecular simulation and computer-aided drug design. Here, we introduce a generalized and extensible machine-learned MM force field, espaloma-0.3, and an end-to-end differentiable framework using graph neural networks to overcome the limitations of traditional rule-based methods. Trained in a single GPU-day to fit a large and diverse quantum chemical dataset of over 1.1 M energy and force calculations, espaloma-0.3 reproduces quantum chemical energetic properties of chemical domains highly relevant to drug discovery, including small molecules, peptides, and nucleic acids. Moreover, this force field maintains the quantum chemical energy-minimized geometries of small molecules and preserves the condensed phase properties of peptides and folded proteins, self-consistently parametrizing proteins and ligands to produce stable simulations leading to highly accurate predictions of binding free energies. This methodology demonstrates significant promise as a path forward for systematically building more accurate force fields that are easily extensible to new chemical domains of interest., Competing Interests: J. D. C. is a current member of the Scientific Advisory Board of OpenEye Scientific Software, Redesign Science, Ventus Therapeutics, and Interline Therapeutics, and has equity interests in Redesign Science and Interline Therapeutics. The Chodera laboratory receives or has received funding from multiple sources, including the National Institutes of Health, the National Science Foundation, the Parker Institute for Cancer Immunotherapy, Relay Therapeutics, Entasis Therapeutics, Silicon Therapeutics, EMD Serono (Merck KGaA), AstraZeneca, Vir Biotechnology, Bayer, XtalPi, Interline Therapeutics, the Molecular Sciences Software Institute, the Starr Cancer Consortium, the Open Force Field Consortium, Cycle for Survival, a Louis V. Gerstner Young Investigator Award, and the Sloan Kettering Institute. A complete funding history for the Chodera lab can be found at http://choderalab.org/funding. Y. W. has limited financial interests in Flagship Pioneering, Inc. and its subsidiaries. M. R. S. is an Open Science Fellow with Psivant Sciences and consults for Relay Therapeutics. D. L. M. serves on the scientific advisory boards of Anagenex and OpenEye Scientific Software, Cadence Molecular Sciences, and is an Open Science Fellow with Psivant., (This journal is © The Royal Society of Chemistry.)

Isothermal titration calorimetry (ITC) is a powerful biophysical tool to characterize energetic profiles of biomacromolecular interactions without any alteration of the underlying chemical structures. In this protocol, we describe procedures for performing, analyzing, and interpreting ITC data obtained from a cooperative riboswitch-ligand interaction., (© 2023. This is a U.S. government work and not under copyright protection in the U.S.; foreign copyright protection may apply.)

The essential pre-mRNA splicing factor U2AF2 (also called U2AF65) identifies polypyrimidine (Py) tract signals of nascent transcripts, despite length and sequence variations. Previous studies have shown that the U2AF2 RNA recognition motifs (RRM1 and RRM2) preferentially bind uridine-rich RNAs. Nonetheless, the specificity of the RRM1/RRM2 interface for the central Py tract nucleotide has yet to be investigated. We addressed this question by determining crystal structures of U2AF2 bound to a cytidine, guanosine, or adenosine at the central position of the Py tract, and compared U2AF2-bound uridine structures. Local movements of the RNA site accommodated the different nucleotides, whereas the polypeptide backbone remained similar among the structures. Accordingly, molecular dynamics simulations revealed flexible conformations of the central, U2AF2-bound nucleotide. The RNA binding affinities and splicing efficiencies of structure-guided mutants demonstrated that U2AF2 tolerates nucleotide substitutions at the central position of the Py tract. Moreover, enhanced UV-crosslinking and immunoprecipitation of endogenous U2AF2 in human erythroleukemia cells showed uridine-sensitive binding sites, with lower sequence conservation at the central nucleotide positions of otherwise uridine-rich, U2AF2-bound splice sites. Altogether, these results highlight the importance of RNA flexibility for protein recognition and take a step towards relating splice site motifs to pre-mRNA splicing efficiencies., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Riboswitches are structured non-coding RNAs often located upstream of essential genes in bacterial messenger RNAs. Such RNAs regulate expression of downstream genes by recognizing a specific cellular effector. Although nearly 50 riboswitch classes are known, only a handful recognize multiple effectors. Here, we report the 2.60-Å resolution co-crystal structure of a class I type I preQ 1 -sensing riboswitch that reveals two effectors stacked atop one another in a single binding pocket. These effectors bind with positive cooperativity in vitro and both molecules are necessary for gene regulation in bacterial cells. Stacked effector recognition appears to be a hallmark of the largest subgroup of preQ 1 riboswitches, including those from pathogens such as Neisseria gonorrhoeae. We postulate that binding to stacked effectors arose in the RNA World to closely position two substrates for RNA-mediated catalysis. These findings expand known effector recognition capabilities of riboswitches and have implications for antimicrobial development., (© 2022. The Author(s).)

RNA recognition by proteins is central to biology. Here we demonstrate the existence of a recurrent structural motif, the "arginine fork", that codifies arginine readout of cognate backbone and guanine nucleobase interactions in a variety of protein-RNA complexes derived from viruses, metabolic enzymes, and ribosomes. Nearly 30 years ago, a theoretical arginine fork model was posited to account for the specificity between the HIV-1 Tat protein and TAR RNA. This model predicted that a single arginine should form four complementary contacts with nearby phosphates, yielding a two-pronged backbone readout. Recent high-resolution structures of TAR-protein complexes have unveiled new details, including ( i ) arginine interactions with the phosphate backbone and the major-groove edge of guanine and ( ii ) simultaneous cation-π contacts between the guanidinium group and flanking nucleobases. These findings prompted us to search for arginine forks within experimental protein-RNA structures retrieved from the Protein Data Bank. The results revealed four distinct classes of arginine forks that we have defined using a rigorous but flexible nomenclature. Examples are presented in the context of ribosomal and nonribosomal interfaces with analysis of arginine dihedral angles and structural (suite) classification of RNA targets. When arginine fork chemical recognition principles were applied to existing structures with unusual arginine-guanine recognition, we found that the arginine fork geometry was more consistent with the experimental data, suggesting the utility of fork classifications to improve structural models. Software to analyze arginine-RNA interactions has been made available to the community.

Riboswitches are structured RNA motifs that recognize metabolites to alter the conformations of downstream sequences, leading to gene regulation. To investigate this molecular framework, we determined crystal structures of a preQ1-I riboswitch in effector-free and bound states at 2.00 Å and 2.65 Å-resolution. Both pseudoknots exhibited the elusive L2 loop, which displayed distinct conformations. Conversely, the Shine-Dalgarno sequence (SDS) in the S2 helix of each structure remained unbroken. The expectation that the effector-free state should expose the SDS prompted us to conduct solution experiments to delineate environmental changes to specific nucleobases in response to preQ1. We then used nudged elastic band computational methods to derive conformational-change pathways linking the crystallographically-determined effector-free and bound-state structures. Pathways featured: (i) unstacking and unpairing of L2 and S2 nucleobases without preQ1-exposing the SDS for translation and (ii) stacking and pairing L2 and S2 nucleobases with preQ1-sequestering the SDS. Our results reveal how preQ1 binding reorganizes L2 into a nucleobase-stacking spine that sequesters the SDS, linking effector recognition to biological function. The generality of stacking spines as conduits for effector-dependent, interdomain communication is discussed in light of their existence in adenine riboswitches, as well as the turnip yellow mosaic virus ribosome sensor., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Natural and lab-evolved proteins often recognize their RNA partners with exquisite affinity. Structural analysis of such complexes can offer valuable insight into sequence-selective recognition that can be exploited to alter biological function. Here, we describe the structure of a lab-evolved RNA recognition motif (RRM) bound to the HIV-1 trans-activation response (TAR) RNA element at 1.80 Å-resolution. The complex reveals a trio of arginines in an evolved β2-β3 loop penetrating deeply into the major groove to read conserved guanines while simultaneously forming cation-π and salt-bridge contacts. The observation that the evolved RRM engages TAR within a double-stranded stem is atypical compared to most RRMs. Mutagenesis, thermodynamic analysis and molecular dynamics validate the atypical binding mode and quantify molecular contributions that support the exceptionally tight binding of the TAR-protein complex (KD,App of 2.5 ± 0.1 nM). These findings led to the hypothesis that the β2-β3 loop can function as a standalone TAR-recognition module. Indeed, short constrained peptides comprising the β2-β3 loop still bind TAR (KD,App of 1.8 ± 0.5 μM) and significantly weaken TAR-dependent transcription. Our results provide a detailed understanding of TAR molecular recognition and reveal that a lab-evolved protein can be reduced to a minimal RNA-binding peptide.

A clear consensus concerning the mechanisms of intracellular secretory vesicle trafficking in astrocytes is lacking in the physiological literature. A good characterization of vesicle trafficking that may assist researchers in achieving that goal is the trajectory angle, defined as the angle between the trajectory of a vesicle and a line radial to the cell's nucleus. In this study, we provide a precise definition of the trajectory angle, describe and compare two methods for its calculation in terms of measureable trafficking parameters, and give recommendations for the appropriate use of each method. We investigated the trafficking of vesicles containing excitatory amino acid transporter 2 (EAAT2) fluorescently tagged with enhanced green fluorescent protein (EGFP) to quantify and validate the precision of each method. The motion of fluorescent puncta--taken to represent vesicles containing EAAT2-EGFP--was found to be typical of secretory vesicle trafficking. An exact method for calculating the trajectory angle of these puncta produced no error but required large computation time. An approximate method reduced the requisite computation time but produced an error depending on the inverse of the ratio of the punctum's initial distance from the nucleus centroid to its maximal displacement. Fitting this dependence to a power function allowed us to establish an exclusion distance from the centroid, beyond which the approximate method is less likely to produce an error above an acceptable 5%. We recommend that the exact method be used to calculate the trajectory angle for puncta closer to the nucleus centroid than this exclusion distance.

Astrocytes provide a principal pathway for glutamate uptake in the mammalian brain, a task accomplished by the powerful action of excitatory amino acid transporters (EAAT) 1 and 2. These transporters are synthesized within the endoplasmic reticulum and are then trafficked to the plasma membrane. The characteristics of their intracellular traffic within astrocytes have not been investigated. We monitored the trafficking of secretory vesicles laden with the recombinant fluorescent protein chimera of EAAT2 in cultured astrocytes. Such vesicles appeared as fluorescent puncta, and their trafficking parameters were obtained using original algorithms, which we describe here in detail. We determined the maximal displacement, average instantaneous speed, and trajectory angle of individual puncta/vesicles, with angles near 0° indicating radial movement directly away from or toward the nucleus and angles near 90° indicating tangential movement. Analysis of these trafficking parameters demonstrated that trafficking of EAAT2-laden vesicles has typical characteristics expected of the trafficking of secretory vesicles in cultured astrocytes. The distribution of trajectory angles for directional vesicles, i.e. those with a maximal displacement greater than 1 μm within the 40-s time-lapse imaging, was found to be unimodal, with angles near 0° being the most prominent (mode 7°). These measurements are in good agreement with previous measurements of trajectory angles of similar trafficking vesicles carrying cannabinoid receptor 1, evidencing the validity and robustness of our analytical approach and algorithms.