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 application, we undertook the structural and dynamical analysis of the most abundant and conserved calcium-binding protein, namely, Calmodulin (CaM). CaM is constituted by two calcium-binding motifs called EF-hands, which in 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 nearnative 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 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 showed the capacity to handle a number of structurally and dynamically challenging situations including metal ion coordination, unbiased conformational sampling, and specific protein-peptide recognition.

Celiac disease (CD) is a chronic enteropathy elicited by a Th1 response to gluten peptides in the small intestine of genetically susceptible individuals. However, it remains unclear what drives the induction of inflammatory responses of this kind against harmless antigens in food. In a recent work, we have shown that the p31-43 peptide (p31-43) from α-gliadin can induce an innate immune response in the intestine and that this may initiate pathological adaptive immunity. The receptors and mechanisms responsible for the induction of innate immunity by p31-43 are unknown and here we present evidence that this may reflect conformational changes in the peptide that allow it to activate the NLRP3 inflammasome. Administration of p31-43, but not scrambled or inverted peptides, to normal mice induced enteropathy in the proximal small intestine, associated with increased production of type I interferon and mature IL-1β. P31-43 showed a sequence-specific spontaneous ability to form structured oligomers and aggregates in vitro and induced activation of the ASC speck complex. In parallel, the enteropathy induced by p31-43 in vivo did not occur in the absence of NLRP3 or caspase 1 and was inhibited by administration of the caspase 1 inhibitor Ac-YVAD-cmk. Collectively, these findings show that p31-43 gliadin has an intrinsic propensity to form oligomers which trigger the NLRP3 inflammasome and that this pathway is required for intestinal inflammation and pathology when p31-43 is administered orally to mice. This innate activation of the inflammasome may have important implications in the initial stages of CD pathogenesis.

Trypanothione synthetase (TryS) produces N 1 ,N 8-bis(glutathionyl)spermidine (or trypanothione) at the expense of ATP. Trypanothione is a metabolite unique and essential for survival and drug-resistance of trypanosomatid parasites. In this study, we report the mechanistic and biological characterisation of optimised N 5-substituted paullone analogues with anti-TryS activity. Several of the new derivatives retained submicromolar IC 50 against leishmanial TryS. The binding mode to TryS of the most potent paullones has been revealed by means of kinetic, biophysical and molecular modelling approaches. A subset of analogues showed an improved potency (EC 50 0.5-10 mM) and selectivity (20-35) against the clinically relevant stage of Leishmania braziliensis (mucocutaneous leishmaniasis) and L. infantum (visceral leishmaniasis). For a selected derivative, the mode of action involved intracellular depletion of trypanothione. Our findings shed light on the molecular interaction of TryS with rationally designed inhibitors and disclose a new set of compounds with on-target activity against different Leishmania species.

The capability to handle highly heterogeneous molecular assemblies in a consistent manner is among the greatest challenges faced when deriving simulation parameters. This is particularly the case for coarse-grained simulations in which chemical functional groups are lumped into effective interaction centers for which transferability between different chemical environments is not guaranteed. Here we introduce the parameterization of a set of CG phospholipids compatible with the latest version of the SIRAH force field for proteins. The newly introduced lipid species include different acylic chain lengths, partial unsaturation, as well as polar and acidic head groups that show a very good reproduction of structural membrane determinants, as areas per lipid, thickness, order parameter, etc., and their dependence with temperature. Simulation of membrane proteins showed unprecedented accuracy in the unbiased description of the thickness-dependent membrane-protein orientation in systems where this information is experimentally available (namely, the SarcoEndoplasmic Reticulum Calcium -SERCA-pump and its regulator Phospholamban). The interactions that lead to this faithful reproduction can be traced down to single amino acid-lipid interaction level and show full agreement with biochemical data present in the literature. Finally, the present parameterization is implemented in the GROMACS and AMBER simulation packages facilitating its use to a wide portion of the Biocomputing community..

Celiac disease (CeD) is a highly prevalent chronic immune-mediated enteropathy developed in genetically predisposed individuals after ingestion of a group of wheat proteins (called gliadins and glutenins). The 13mer agliadin peptide, p31-43, induces proinflammatory responses, observed by in vitro assays and animal models, that may contribute to innate immune mechanisms of CeD pathogenesis. Since a cellular receptor for p31-43 has not been identified, this raises the question of whether this peptide could mediate different biological effects. In this work, we aimed to characterize the p31-43 secondary structure by different biophysical and in silico techniques. By dynamic light scattering and using an oligomer/fibril-sensitive fluorescent probe, we showed the presence of oligomers of this peptide in solution. Furthermore, atomic force microscopy analysis showed p31-43 oligomers with different height distribution. Also, peptide concentration had a very strong influence on peptide self-organization process. Oligomers gradually increased their size at lower concentration. Whereas, at higher ones, oligomers increased their complexity, forming branched structures. By CD, we observed that p31-43 self-organized in a polyproline II conformation in equilibrium with b-sheets-like structures, whose pH remained stable in the range of 3-8. In addition, these findings were supported by molecular dynamics simulation. The formation of p31-43 nanostructures with increased b-sheet structure may help to explain the molecular etiopathogenesis in the induction of proinflammatory effects and subsequent damage at the intestinal mucosa in CeD.