Search

Chemical probing experiments such as SHAPE are routinely used to probe RNA molecules. In this work, we use atomistic molecular dynamics simulations to test the hypothesis that binding of RNA with SHAPE reagents is affected by cooperative effects leading to an observed reactivity that is dependent on the reagent concentration. We develop a general technique that enables the calculation of the affinity for arbitrary molecules as a function of their concentration in the grand-canonical ensemble. Our simulations of an RNA structural motif suggest that, at the concentration typically used in SHAPE experiments, cooperative binding would lead to a measurable concentration-dependent reactivity. We also provide a qualitative validation of this statement by analyzing a new set of experiments collected at different reagent concentrations.

Biomolecular force fields have been traditionally derived based on a mixture of reference quantum chemistry data and experimental information obtained on small fragments. However, the possibility to run extensive molecular dynamics simulations on larger systems achieving ergodic sampling is paving the way to directly using such simulations along with solution experiments obtained on macromolecular systems. Recently, a number of methods have been introduced to automatize this approach. Here we review these methods, highlight their relationship with machine learning methods, and discuss the open challenges in the field., Comment: Supporting information included in ancillary files. The article has been accepted by J. Chem. Phys

RNA function crucially depends on its structure. Thermodynamic models currently used for secondary structure prediction rely on computing the partition function of folding ensembles, and can thus estimate minimum free-energy structures and ensemble populations. These models sometimes fail in identifying native structures unless complemented by auxiliary experimental data. Here, we build a set of models that combine thermodynamic parameters, chemical probing data (DMS, SHAPE), and co-evolutionary data (Direct Coupling Analysis, DCA) through a network that outputs perturbations to the ensemble free energy. Perturbations are trained to increase the ensemble populations of a representative set of known native RNA structures. In the chemical probing nodes of the network, a convolutional window combines neighboring reactivities, enlightening their structural information content and the contribution of local conformational ensembles. Regularization is used to limit overfitting and improve transferability. The most transferable model is selected through a cross-validation strategy that estimates the performance of models on systems on which they are not trained. With the selected model we obtain increased ensemble populations for native structures and more accurate predictions in an independent validation set. The flexibility of the approach allows the model to be easily retrained and adapted to incorporate arbitrary experimental information., Comment: Supporting Information are reported in a Jupyter notebook, in both format .pdf (read-only) and .ipynb (interactive and editable). Full data required to run the .ipynb notebook can be found at the GitHub repository reported in the notebook

The relation between elasticity and yielding is investigated in a model polymer solid by Molecular-Dynamics simulations. By changing the bending stiffness of the chain and the bond length, semicrystalline and disordered glassy polymers - both with bond disorder - as well as nematic glassy polymers with bond ordering are obtained. It is found that in systems with bond disorder the ratio Ty/G between the shear yield strength Ty and the shear modulus G is close to the universal value of the atomic metallic glasses. The increase of the local nematic order in glasses leads to the increase of the shear modulus and the decrease of the shear yield strength, as observed in experiments on nematic thermosets. A tentative explanation of the subsequent reduction of the ratio Ty/G in terms of the distributions of the per-monomer stress is offered.