1
|
Trewhella J, Vachette P, Larsen AH. Benchmarking predictive methods for small-angle X-ray scattering from atomic coordinates of proteins using maximum likelihood consensus data. IUCRJ 2024; 11:762-779. [PMID: 38989800 PMCID: PMC11364021 DOI: 10.1107/s205225252400486x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Accepted: 05/23/2024] [Indexed: 07/12/2024]
Abstract
Stimulated by informal conversations at the XVII International Small Angle Scattering (SAS) conference (Traverse City, 2017), an international team of experts undertook a round-robin exercise to produce a large dataset from proteins under standard solution conditions. These data were used to generate consensus SAS profiles for xylose isomerase, urate oxidase, xylanase, lysozyme and ribonuclease A. Here, we apply a new protocol using maximum likelihood with a larger number of the contributed datasets to generate improved consensus profiles. We investigate the fits of these profiles to predicted profiles from atomic coordinates that incorporate different models to account for the contribution to the scattering of water molecules of hydration surrounding proteins in solution. Programs using an implicit, shell-type hydration layer generally optimize fits to experimental data with the aid of two parameters that adjust the volume of the bulk solvent excluded by the protein and the contrast of the hydration layer. For these models, we found the error-weighted residual differences between the model and the experiment generally reflected the subsidiary maxima and minima in the consensus profiles that are determined by the size of the protein plus the hydration layer. By comparison, all-atom solute and solvent molecular dynamics (MD) simulations are without the benefit of adjustable parameters and, nonetheless, they yielded at least equally good fits with residual differences that are less reflective of the structure in the consensus profile. Further, where MD simulations accounted for the precise solvent composition of the experiment, specifically the inclusion of ions, the modelled radius of gyration values were significantly closer to the experiment. The power of adjustable parameters to mask real differences between a model and the structure present in solution is demonstrated by the results for the conformationally dynamic ribonuclease A and calculations with pseudo-experimental data. This study shows that, while methods invoking an implicit hydration layer have the unequivocal advantage of speed, care is needed to understand the influence of the adjustable parameters. All-atom solute and solvent MD simulations are slower but are less susceptible to false positives, and can account for thermal fluctuations in atomic positions, and more accurately represent the water molecules of hydration that contribute to the scattering profile.
Collapse
Affiliation(s)
- Jill Trewhella
- School of Life and Environmental SciencesUniversity of SydneyNSW2006Australia
| | - Patrice Vachette
- Institute for Integrative Biology of the Cell (12BC)Université Paris-Saclay, CEA, CNRSGif-sur-YvetteParis91198France
| | - Andreas Haahr Larsen
- Department of NeuroscienceUniversity of CopenhagenBlegdamsvej 32200CopenhagenDenmark
| |
Collapse
|
2
|
Muscat S, Martino G, Manigrasso J, Marcia M, De Vivo M. On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules. J Chem Theory Comput 2024. [PMID: 39150960 DOI: 10.1021/acs.jctc.4c00773] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/18/2024]
Abstract
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.
Collapse
Affiliation(s)
- Stefano Muscat
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Gianfranco Martino
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Jacopo Manigrasso
- Medicinal Chemistry, Research and Early Development, Cardiovascular, Renal and Metabolism (CVRM), BioPharmaceuticals R&D, AstraZeneca, 431 50 Mölndal, Sweden
| | - Marco Marcia
- European Molecular Biology Laboratory Grenoble, 71 Avenue des Martyrs, 38042 Grenoble, France
- Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, 751 23 Uppsala, Sweden
- Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| | - Marco De Vivo
- Laboratory of Molecular Modelling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genoa, Italy
| |
Collapse
|
3
|
Sun LZ, Qian JL, Cai P, Xu X. Mutual effects between single-stranded DNA conformation and Na +-Mg 2+ ion competition in mixed salt solutions. Phys Chem Chem Phys 2022; 24:20867-20881. [PMID: 36043348 DOI: 10.1039/d2cp02737b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The ion-dependence of single-stranded DNA (ssDNA) conformational changes has attracted growing attention because of its biological and technological importance. Although single-species ion effects have been extensively explored, it is challenging to study the ssDNA conformational properties under mixed monovalent/divalent ion conditions due to the complications of ssDNA flexibility and ion-ion competition. In this study, we apply Langevin dynamics simulations to investigate mixed Na+/Mg2+ ion-dependent ssDNA conformations. The ssDNA structure is described using a coarse-grained model, in which the phosphate, base, and sugar of each nucleotide are represented by three different beads. A novel improvement in our simulation model is that mixed-salt-related electrostatic interactions are computed via combining Manning counterion condensation (MCC) theory with the Monte Carlo tightly bound ion (MCTBI) model. Based on this MCC-MCTBI combination, we report new empirical functions to describe the ion-concentration-dependent and ssDNA conformation/structure-dependent electrostatic effects. The calculation results relating to the ion binding properties and the simulation results relating to the ssDNA conformational properties are validated against experimental results. In addition, our simulation results suggest a quantitative relationship between the ssDNA conformation and Na+-Mg2+ competition; this in turn reveals their mutual impact in the ion atmosphere.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China.
| | - Jun-Lin Qian
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China.
| | - Pinggen Cai
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China.
| | - Xiaojun Xu
- Institute of Bioinformatics and Medical Engineering, Jiangsu University of Technology, Changzhou, 213001, China
| |
Collapse
|
4
|
Predicting solution scattering patterns with explicit-solvent molecular simulations. Methods Enzymol 2022; 677:433-456. [DOI: 10.1016/bs.mie.2022.08.035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
5
|
Bernetti M, Hall KB, Bussi G. Reweighting of molecular simulations with explicit-solvent SAXS restraints elucidates ion-dependent RNA ensembles. Nucleic Acids Res 2021; 49:e84. [PMID: 34107023 PMCID: PMC8373061 DOI: 10.1093/nar/gkab459] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2021] [Revised: 05/07/2021] [Accepted: 05/16/2021] [Indexed: 01/03/2023] Open
Abstract
Small-angle X-ray scattering (SAXS) experiments are increasingly used to probe RNA structure. A number of forward models that relate measured SAXS intensities and structural features, and that are suitable to model either explicit-solvent effects or solute dynamics, have been proposed in the past years. Here, we introduce an approach that integrates atomistic molecular dynamics simulations and SAXS experiments to reconstruct RNA structural ensembles while simultaneously accounting for both RNA conformational dynamics and explicit-solvent effects. Our protocol exploits SAXS pure-solute forward models and enhanced sampling methods to sample an heterogenous ensemble of structures, with no information towards the experiments provided on-the-fly. The generated structural ensemble is then reweighted through the maximum entropy principle so as to match reference SAXS experimental data at multiple ionic conditions. Importantly, accurate explicit-solvent forward models are used at this reweighting stage. We apply this framework to the GTPase-associated center, a relevant RNA molecule involved in protein translation, in order to elucidate its ion-dependent conformational ensembles. We show that (a) both solvent and dynamics are crucial to reproduce experimental SAXS data and (b) the resulting dynamical ensembles contain an ion-dependent fraction of extended structures.
Collapse
Affiliation(s)
- Mattia Bernetti
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy
| | - Kathleen B Hall
- Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, MO 63110, USA
| | - Giovanni Bussi
- Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, Trieste 34136, Italy
| |
Collapse
|
6
|
Abstract
The biological functions of RNA range from gene regulation through catalysis and depend critically on its structure and flexibility. Conformational variations of flexible, non-base-paired components, including RNA hinges, bulges, or single-stranded tails, are well documented. Recent work has also identified variations in the structure of ubiquitous, base-paired duplexes found in almost all functional RNAs. Duplexes anchor the structures of folded RNAs, and their surface features are recognized by partner molecules. To date, no consistent picture has been obtained that describes the range of conformations assumed by RNA duplexes. Here, we apply wide angle, solution X-ray scattering (WAXS) to quantify these variations, by sampling length scales characteristic of helical geometries under different solution conditions. To identify the radius, helical rise, twist, and length of dsRNA helices, we exploit molecular dynamics generated structures, explicit solvent models, and ensemble optimization methods. Our results quantify the substantial and salt-dependent deviations of double-stranded (ds) RNA duplexes from the assumed canonical A-form conformation. Recent experiments underscore the need to properly describe the structures of RNA duplexes when interpreting the salt dependence of RNA conformations.
Collapse
Affiliation(s)
- Yen-Lin Chen
- School of Applied and Engineering Physics, Cornell University, Ithaca, NY 14853
| | - Lois Pollack
- School of Applied and Engineering Physics, Cornell University, Ithaca, NY 14853
| |
Collapse
|
7
|
Berntsson O, Rodriguez R, Henry L, Panman MR, Hughes AJ, Einholz C, Weber S, Ihalainen JA, Henning R, Kosheleva I, Schleicher E, Westenhoff S. Photoactivation of Drosophila melanogaster cryptochrome through sequential conformational transitions. SCIENCE ADVANCES 2019; 5:eaaw1531. [PMID: 31328161 PMCID: PMC6636987 DOI: 10.1126/sciadv.aaw1531] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/22/2018] [Accepted: 06/13/2019] [Indexed: 05/27/2023]
Abstract
Cryptochromes are blue-light photoreceptor proteins, which provide input to circadian clocks. The cryptochrome from Drosophila melanogaster (DmCry) modulates the degradation of Timeless and itself. It is unclear how light absorption by the chromophore and the subsequent redox reactions trigger these events. Here, we use nano- to millisecond time-resolved x-ray solution scattering to reveal the light-activated conformational changes in DmCry and the related (6-4) photolyase. DmCry undergoes a series of structural changes, culminating in the release of the carboxyl-terminal tail (CTT). The photolyase has a simpler structural response. We find that the CTT release in DmCry depends on pH. Mutation of a conserved histidine, important for the biochemical activity of DmCry, does not affect transduction of the structural signal to the CTT. Instead, molecular dynamics simulations suggest that it stabilizes the CTT in the resting-state conformation. Our structural photocycle unravels the first molecular events of signal transduction in an animal cryptochrome.
Collapse
Affiliation(s)
- Oskar Berntsson
- Department of Chemistry and Molecular Biology, University of Gothenburg, 40530 Gothenburg, Sweden
- MAX IV Laboratory, Lund University, 224 84 Lund, Sweden
| | - Ryan Rodriguez
- Institute of Physical Chemistry, Albert-Ludwigs-Universität Freiburg, 79104 Freiburg, Germany
| | - Léocadie Henry
- Department of Chemistry and Molecular Biology, University of Gothenburg, 40530 Gothenburg, Sweden
| | - Matthijs R. Panman
- Department of Chemistry and Molecular Biology, University of Gothenburg, 40530 Gothenburg, Sweden
| | - Ashley J. Hughes
- Department of Chemistry and Molecular Biology, University of Gothenburg, 40530 Gothenburg, Sweden
| | - Christopher Einholz
- Institute of Physical Chemistry, Albert-Ludwigs-Universität Freiburg, 79104 Freiburg, Germany
| | - Stefan Weber
- Institute of Physical Chemistry, Albert-Ludwigs-Universität Freiburg, 79104 Freiburg, Germany
| | - Janne A. Ihalainen
- Nanoscience Center, Department of Biological and Environmental Sciences, University of Jyväskylä, 40014 Jyväskylä, Finland
| | - Robert Henning
- Center for Advanced Radiation Sources, The University of Chicago, Chicago, IL 60637, USA
| | - Irina Kosheleva
- Center for Advanced Radiation Sources, The University of Chicago, Chicago, IL 60637, USA
| | - Erik Schleicher
- Institute of Physical Chemistry, Albert-Ludwigs-Universität Freiburg, 79104 Freiburg, Germany
| | - Sebastian Westenhoff
- Department of Chemistry and Molecular Biology, University of Gothenburg, 40530 Gothenburg, Sweden
| |
Collapse
|