ANNaMo: Coarse-grained modeling for folding and assembly of RNA and DNA systems.

The folding of RNA and DNA strands plays crucial roles in biological systems and bionanotechnology. However, studying these processes with high-resolution numerical models is beyond current computational capabilities due to the timescales and system sizes involved. In this article, we present a new coarse-grained model for investigating the folding dynamics of nucleic acids. Our model represents three nucleotides with a patchy particle and is parameterized using well-established nearest-neighbor models. Thanks to the reduction of degrees of freedom and to a bond-swapping mechanism, our model allows for simulations at timescales and length scales that are currently inaccessible to more detailed models. To validate the performance of our model, we conducted extensive simulations of various systems: We examined the thermodynamics of DNA hairpins, capturing their stability and structural transitions, the folding of an MMTV pseudoknot, which is a complex RNA structure involved in viral replication, and also explored the folding of an RNA tile containing a k-type pseudoknot. Finally, we evaluated the performance of the new model in reproducing the melting temperatures of oligomers and the dependence on the toehold length of the displacement rate in toehold-mediated displacement processes, a key reaction used in molecular computing. All in all, the successful reproduction of experimental data and favorable comparisons with existing coarse-grained models validate the effectiveness of the new model.

A Primer on the oxDNA Model of DNA: When to Use it, How to Simulate it and How to Interpret the Results.

The oxDNA model of Deoxyribonucleic acid has been applied widely to systems in biology, biophysics and nanotechnology. It is currently available via two independent open source packages. Here we present a set of clearly documented exemplar simulations that simultaneously provide both an introduction to simulating the model, and a review of the model's fundamental properties. We outline how simulation results can be interpreted in terms of-and feed into our understanding of-less detailed models that operate at larger length scales, and provide guidance on whether simulating a system with oxDNA is worthwhile.

Modelling realistic microgels in an explicit solvent.

Thermoresponsive microgels are polymeric colloidal networks that can change their size in response to a temperature variation. This peculiar feature is driven by the nature of the solvent-polymer interactions, which triggers the so-called volume phase transition from a swollen to a collapsed state above a characteristic temperature. Recently, an advanced modelling protocol to assemble realistic, disordered microgels has been shown to reproduce experimental swelling behavior and form factors. In the original framework, the solvent was taken into account in an implicit way, condensing solvent-polymer interactions in an effective attraction between monomers. To go one step further, in this work we perform simulations of realistic microgels in an explicit solvent. We identify a suitable model which fully captures the main features of the implicit model and further provides information on the solvent uptake by the interior of the microgel network and on its role in the collapse kinetics. These results pave the way for addressing problems where solvent effects are dominant, such as the case of microgels at liquid-liquid interfaces.

Small-angle neutron scattering and molecular dynamics structural study of gelling DNA nanostars.

DNA oligomers with properly designed sequences self-assemble into well defined constructs. Here, we exploit this methodology to produce bulk quantities of tetravalent DNA nanostars (each one composed of 196 nucleotides) and to explore the structural signatures of their aggregation process. We report small-angle neutron scattering experiments focused on the evaluation of both the form factor and the temperature evolution of the scattered intensity at a nanostar concentration where the system forms a tetravalent equilibrium gel. We also perform molecular dynamics simulations of one isolated tetramer to evaluate the form factor numerically, without resorting to any approximate shape. The numerical form factor is found to be in very good agreement with the experimental one. Simulations predict an essentially temperature-independent form factor, offering the possibility to extract the effective structure factor and its evolution during the equilibrium gelation.