Inferring Structural Ensembles of Flexible and Dynamic Macromolecules Using Bayesian, Maximum Entropy, and Minimal-Ensemble Refinement Methods

Accounting for fast vs slow exchange in single molecule FRET experiments reveals hidden conformational states

Proteins are dynamic systems whose structural preferences determine their function. Unfortunately, building atomically detailed models of protein structural ensembles remains challenging, limiting our understanding of the relationships between sequence, structure, and function. Combining single molecule Förster resonance energy transfer (smFRET) experiments with molecular dynamics simulations could provide experimentally grounded, all-atom models of a protein's structural ensemble. However, agreement between the two techniques is often insufficient to achieve this goal. Here, we explore whether accounting for important experimental details like averaging across structures sampled during a given smFRET measurement is responsible for this apparent discrepancy. We present an approach to account for this time-averaging by leveraging the kinetic information available from Markov state models of a protein's dynamics. This allows us to accurately assess which timescales are averaged during an experiment. We find this approach significantly improves agreement between simulations and experiments in proteins with varying degrees of dynamics, including the well-ordered protein T4 lysozyme, the partially disordered protein apolipoprotein E (ApoE), and a disordered amyloid protein (Aβ40). We find evidence for hidden states that are not apparent in smFRET experiments because of time averaging with other structures, akin to states in fast exchange in NMR, and evaluate different force fields. Finally, we show how remaining discrepancies between computations and experiments can be used to guide additional simulations and build structural models for states that were previously unaccounted for. We expect our approach will enable combining simulations and experiments to understand the link between sequence, structure, and function in many settings.

Simple Adjustment of Intranucleotide Base-Phosphate Interaction in the OL3 AMBER Force Field Improves RNA Simulations

Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field (ff). Here we suggest a correction for the widely used AMBER OL3 ff by adding a simple adjustment of the nonbonded parameters. The reparameterization of the Lennard-Jones potential for the -H8···O5'- and -H6···O5'- atom pairs addresses an intranucleotide steric clash occurring in the type 0 base-phosphate interaction (0BPh). The nonbonded fix (NBfix) modification of 0BPh interactions (NBfix0BPh modification) was tuned via a reweighting approach and subsequently tested using an extensive set of standard and enhanced sampling simulations of both unstructured and folded RNA motifs. The modification corrects minor but visible intranucleotide clash for the anti nucleobase conformation. We observed that structural ensembles of small RNA benchmark motifs simulated with the NBfix0BPh modification provide better agreement with experiments. No side effects of the modification were observed in standard simulations of larger structured RNA motifs. We suggest that the combination of OL3 RNA ff and NBfix0BPh modification is a viable option to improve RNA MD simulations.

Conformational ensemble of the full-length SARS-CoV-2 nucleocapsid (N) protein based on molecular simulations and SAXS data

The nucleocapsid protein of the SARS-CoV-2 virus comprises two RNA-binding domains and three regions that are intrinsically disordered. While the structures of the RNA-binding domains have been solved using protein crystallography and NMR, current knowledge of the conformations of the full-length nucleocapsid protein is rather limited. To fill in this knowledge gap, we combined coarse-grained molecular simulations with data from small-angle X-ray scattering (SAXS) experiments using the ensemble refinement of SAXS (EROS) method. Our results show that the dimer of the full-length nucleocapsid protein exhibits large conformational fluctuations with its radius of gyration ranging from about 4 to 8 nm. The RNA-binding domains do not make direct contacts. The disordered region that links these two domains comprises a hydrophobic α-helix which makes frequent and nonspecific contacts with the RNA-binding domains. Each of the intrinsically disordered regions adopts conformations that are locally compact, yet on average, much more extended than Gaussian chains of equivalent lengths. We offer a detailed picture of the conformational ensemble of the nucleocapsid protein dimer under near-physiological conditions, which will be important for understanding the nucleocapsid assembly process.

How to Determine Accurate Conformational Ensembles by Metadynamics Metainference: A Chignolin Study Case

The reliability and usefulness of molecular dynamics simulations of equilibrium processes rests on their statistical precision and their capability to generate conformational ensembles in agreement with available experimental knowledge. Metadynamics Metainference (M&M), coupling molecular dynamics with the enhanced sampling ability of Metadynamics and with the ability to integrate experimental information of Metainference, can in principle achieve both goals. Here we show that three different Metadynamics setups provide converged estimate of the populations of the three-states populated by a model peptide. Errors are estimated correctly by block averaging, but higher precision is obtained by performing independent replicates. One effect of Metadynamics is that of dramatically decreasing the number of effective frames resulting from the simulations and this is relevant for M&M where the number of replicas should be large enough to capture the conformational heterogeneity behind the experimental data. Our simulations allow also us to propose that monitoring the relative error associated with conformational averaging can help to determine the minimum number of replicas to be simulated in the context of M&M simulations. Altogether our data provides useful indication on how to generate sound conformational ensemble in agreement with experimental data.

Unraveling the NaCl Concentration Effect on the First Stages of α-Synuclein Aggregation

Intraneuronal aggregation of the intrinsically disordered protein α-synuclein is at the core of Parkinson's disease and related neurodegenerative disorders. Several reports show that the concentration of salts in the medium heavily affects its aggregation rate and fibril morphology, but a characterization of the individual monomeric conformations underlying these effects is still lacking. In this work, we have applied our α-synuclein-optimized coarse-grained molecular dynamics approach to decipher the structural features of the protein monomer under a range of NaCl concentrations (0.0-1.0 M). The results show that key intramolecular contacts between the terminal domains are lost at intermediate concentrations (leading to extended conformations likely to fibrillate), but recovered at high concentrations (leading to compact conformations likely to evolve toward amorphous aggregates). The pattern of direct interactions of the terminal α-synuclein domains with Na⁺ and Cl^- ions plays a key role in explaining this effect. Our results are consistent with a recent study reporting a fibrillation enhancement at moderate NaCl concentrations but an inhibition at higher concentrations. The present work will contribute to improving our understanding of the structural features of monomeric α-synuclein, determining its NaCl-induced fibrillation propensity and the molecular basis of synucleinopathies, necessary for the future development of disease-halting therapies.

Capturing the Conformational Ensemble of the Mixed Folded Polyglutamine Protein Ataxin-3

Ataxin-3 is a deubiquitinase involved in protein quality control and other essential cellular functions. It preferentially interacts with polyubiquitin chains of four or more units attached to proteins delivered to the ubiquitin-proteasome system. Ataxin-3 is composed of an N-terminal Josephin domain and a flexible C terminus that contains two or three ubiquitin-interacting motifs (UIMs) and a polyglutamine tract, which, when expanded beyond a threshold, leads to protein aggregation and misfolding and causes spinocerebellar ataxia type 3. The high-resolution structure of the Josephin domain is available, but the structural and dynamical heterogeneity of ataxin-3 has so far hindered the structural description of the full-length protein. Here, we characterize non-expanded and expanded variants of ataxin-3 in terms of conformational ensembles adopted by the proteins in solution by jointly using experimental data from nuclear magnetic resonance and small-angle X-ray scattering with coarse-grained simulations. Our results pave the way to a molecular understanding of polyubiquitin recognition.