Assessing the Accuracy of the SIRAH Force Field to Model DNA at Coarse Grain Level

SIRAH Late Harvest: Coarse-Grained Models for Protein Glycosylation

Glycans constitute one of the most complex families of biological molecules. Despite their crucial role in a plethora of biological processes, they remain largely uncharacterized because of their high complexity. Their intrinsic flexibility and the vast variability associated with the many combination possibilities have hampered their experimental determination. Although theoretical methods have proven to be a valid alternative to the study of glycans, the large size associated with polysaccharides, proteoglycans, and glycolipids poses significant challenges to a fully atomistic description of biologically relevant glycoconjugates. On the other hand, the exquisite dependence on hydrogen bonds to determine glycans' structure makes the development of simplified or coarse-grained (CG) representations extremely challenging. This is particularly the case when glycan representations are expected to be compatible with CG force fields that include several molecular types. We introduce a CG representation able to simulate a wide variety of polysaccharides and common glycosylation motifs in proteins, which is fully compatible with the CG SIRAH force field. Examples of application to N-glycosylated proteins, including antibody recognition and calcium-mediated glycan-protein interactions, highlight the versatility of the enlarged set of CG molecules provided by SIRAH.

The SIRAH force field: A suite for simulations of complex biological systems at the coarse-grained and multiscale levels

The different combinations of molecular dynamics simulations with coarse-grained representations have acquired considerable popularity among the scientific community. Especially in biocomputing, the significant speedup granted by simplified molecular models opened the possibility of increasing the diversity and complexity of macromolecular systems, providing realistic insights on large assemblies for more extended time windows. However, a holistic view of biological ensembles' structural and dynamic features requires a self-consistent force field, namely, a set of equations and parameters that describe the intra and intermolecular interactions among moieties of diverse chemical nature (i.e., nucleic and amino acids, lipids, solvent, ions, etc.). Nevertheless, examples of such force fields are scarce in the literature at the fully atomistic and coarse-grained levels. Moreover, the number of force fields capable of handling simultaneously different scales is restricted to a handful. Among those, the SIRAH force field, developed in our group, furnishes a set of topologies and tools that facilitate the setting up and running of molecular dynamics simulations at the coarse-grained and multiscale levels. SIRAH uses the same classical pairwise Hamiltonian function implemented in the most popular molecular dynamics software. In particular, it runs natively in AMBER and Gromacs engines, and porting it to other simulation packages is straightforward. This review describes the underlying philosophy behind the development of SIRAH over the years and across families of biological molecules, discussing current limitations and future implementations.

The automated optimisation of a coarse-grained force field using free energy data

Atomistic models provide a detailed representation of molecular systems, but are sometimes inadequate for simulations of large systems over long timescales. Coarse-grained models enable accelerated simulations by reducing the number of degrees of freedom, at the cost of reduced accuracy. New optimisation processes to parameterise these models could improve their quality and range of applicability. We present an automated approach for the optimisation of coarse-grained force fields, by reproducing free energy data derived from atomistic molecular simulations. To illustrate the approach, we implemented hydration free energy gradients as a new target for force field optimisation in ForceBalance and applied it successfully to optimise the un-charged side-chains and the protein backbone in the SIRAH protein coarse-grain force field. The optimised parameters closely reproduced hydration free energies of atomistic models and gave improved agreement with experiment.

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.

From quantum to subcellular scales: multi-scale simulation approaches and the SIRAH force field

Modern molecular and cellular biology profits from astonishing resolution structural methods, currently even reaching the whole cell level. This is encompassed by the development of computational methods providing a deep view into the structure and dynamics of molecular processes happening at very different scales in time and space. Linking such scales is of paramount importance when aiming at far-reaching biological questions. Computational methods at the interface between classical and coarse-grained resolutions are gaining momentum with several research groups dedicating important efforts to their development and tuning. An overview of such methods is addressed herein, with special emphasis on the SIRAH force field for coarse-grained and multi-scale simulations. Moreover, we provide proof of concept calculations on the implementation of a multi-scale simulation scheme including quantum calculations on a classical fine-grained/coarse-grained representation of double-stranded DNA. This opens the possibility to include the effect of large conformational fluctuations in chromatin segments on, for instance, the reactivity of particular base pairs within the same simulation framework.

The SIRAH 2.0 Force Field: Altius, Fortius, Citius

A new version of the coarse-grained (CG) SIRAH force field for proteins has been developed. Modifications to bonded and non-bonded interactions on the existing molecular topologies significantly ameliorate the structural description and flexibility of a non-redundant set of proteins. The SIRAH 2.0 force field has also been ported to the popular simulation package AMBER, which along with the former implementation in GROMACS expands significantly the potential range of users and performance of this CG force field on CPU/GPU codes. As a non-trivial example of its application, we undertook the structural and dynamical analysis of the most abundant and conserved calcium-binding protein, calmodulin (CaM). CaM is composed of two calcium-binding motifs called EF-hands, which in the presence of calcium specifically recognize a cognate peptide by embracing it. CG simulations of CaM bound to four calcium ions in the presence or absence of a binding peptide (holo and apo forms, respectively) resulted in good and stable ion coordination. The simulation of the holo form starting from an experimental structure sampled near-native conformations, retrieving quasi-atomistic precision. Removing the binding peptide enabled the EF-hands to perform large reciprocal movements, comparable to those observed in NMR structures. On the other hand, the isolated peptide starting from the helical conformation experienced spontaneous unfolding, in agreement with previous experimental data. However, repositioning the peptide in the neighborhood of one EF-hand not only prevented the peptide from unfolding but also drove CaM to a fully bound conformation, with both EF-hands embracing the cognate peptide, resembling the experimental holo structure. Therefore, SIRAH 2.0 shows the capacity to handle a number of structurally and dynamically challenging situations, including metal ion coordination, unbiased conformational sampling, and specific protein-peptide recognition.

Enantioselective Catalysis by Using Short, Structurally Defined DNA Hairpins as Scaffold for Hybrid Catalysts

A new type of DNA metal complex hybrid catalyst, which is based on single-stranded DNA oligonucleotides, is described. It was shown that oligonucleotides as short as 14 nucleotides that fold into hairpin structures are suitable as nucleic acid components for DNA hybrid catalysts. With these catalysts, excellent enantioinduction in asymmetric Diels-Alder reactions with selectivity values as high as 96 % enantiomeric excess (ee) can be achieved. Molecular dynamics simulations indicate that a rather flexible loop combined with a rigid stem region provides DNA scaffolds with these high selectivity values.

A Solvent-Mediated Coarse-Grained Model of DNA Derived with the Systematic Newton Inversion Method

We present a new class of coarse-grained (CG) force fields (FFs) for B-DNA with explicit ions suited for large-scale mesoscale simulations at microsecond-micrometer scale using a wide spectrum of particle simulation methods from molecular dynamics to dissipative particle dynamics. The effective solvent-mediated pairwise interactions making up the FFs are obtained by inverting radial distribution functions and other particle-particle distributions obtained from all-atom simulations of numbers of octadecamer DNA fragments from the Ascona B-DNA library. The inverse Monte Carlo (IMC) method, later known as Newton inversion (NI) (Lyubartsev, A. P.; Laaksonen, A. Phys. Rev. E, 1995, 52, 3730-3737), was used together with the iterative Boltzmann inversion (IBI) scheme to compute the effective CG potentials. We show that this systematic structure-based approach is capable of providing converged potentials that accurately reproduce the structural features of the underlying atomistic system within a few percents of relative difference. We also show that a simple one-site-per-nucleotide model with 10 intramolecular pair interaction potentials is able to reproduce key features of DNA, for example, the persistence length and its dependence on the ionic concentration, experimentally determined around 50 nm at physiological salt concentration.

Exploring LacI-DNA dynamics by multiscale simulations using the SIRAH force field

The lac repressor protein (LacI) together with its target regulatory sequence are a common model for studying DNA looping and its implications on transcriptional control in bacteria. Owing to the molecular size of this system, standard all-atom (AA) simulations are prohibitive for achieving relevant biological time scales. As an alternative, multiscale models, which combine AA descriptions at particular regions with coarse-grained (CG) representations of the remaining components, were used to address this computational challenge while preserving the relevant details of the system. In this work, we implement a new multiscale approach based on the SIRAH force field to gain deeper insights into the dynamics of the LacI-DNA system. Our methodology allows for a dual resolution treatment of the solute and solvent, explicitly representing the protein, DNA, and solvent environment without compromising the AA region. Starting from the P1 loop configuration in an undertwisted conformation, we were able to observe the transition to the more stable overtwisted state. Additionally, a detailed characterization of the conformational space sampled by the DNA loop was done. In agreement with experimental and theoretical evidence, we observed the transient formation of kinks at the loop, which were stabilized by the presence of counterions at the minor groove. We also show that the loop's intrinsic flexibility can account for reported FRET measurements and bent conformations required to bind the CAP transcription factor.