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.

In silico induced effect of N-ε-lysine acetylation on microtubule stability and subsequent interaction of microtubule-associated proteins.

Plant systems have been considered valuable models for addressing fundamental questions of microtubule (MT) organization due to their considerable practical utility. Protein acetylation is a very common protein modification, and therate of acetylation can be modulated in cells in different biological states, and these changes can be detected at a molecular level. Here, we focused on K40, K112, and K394 residues as putative acetylation sites, which were shown to exist in both plants and mammals. Such residual effect of acetylation causes critical but unclear effect on MT stability. In turn, it was shown that acetylation indirectly affects the probability of interaction with different MAPs (Microtubule-associated proteins). In a multiscale study using an all-atom force field to reproduce several lattice-forming elements found on the surface the microtubule, we assembled a fragment of a plant microtubule composed of nine tubulins and used it as a model object along with the existing human complex. Triplets of tubulins assembled in a lattice cell were then simulated for both human and plant protein complexes, using a coarse-grained force field. We then analyzed the trajectories and identified some critical deformations of the MAP interaction surface. The initial coordinates were used to investigate the structural scenario in which autophagy-related protein 8 (ATG8) was able to interact with the MT fragment.

The Mechanism of Antimicrobial Small-Cationic Peptides from Coarse-Grained Simulations.

Antimicrobial cationic peptides (AMPs) are excellent candidates for use as therapeutic antimicrobial agents. Among them, short peptides possessing sequences of 9-11 amino acids have some advantages over long-sequence peptides. However, one of the main limitations of short peptides is that their mechanism of action at the molecular level is not well-known. In this article, we report a model based on multiscale molecular dynamics simulations of short peptides interacting with vesicles containing palmitoyl-oleoyl-phosphatidylglycerol (POPG)/palmitoyl-oleoyl-phosphatidylethanolamine (POPE). Simulations using this approach have allowed us to understand the different behaviors of peptides with antimicrobial activity with respect to those that do not produce this effect. We found remarkable agreement with a series of experimental results directly supporting our model. Moreover, these results allow us to understand the mechanism of action at the molecular level of these short peptides. Our simulations suggest that mechanical inhomogeneities appear in the membrane, promoting membrane rupture when a threshold concentration of peptides adsorbed on the membrane is achieved. These results explain the high structural demand for these peptides to maintain a delicate balance between the affinity for the bilayer surface, a low peptide-peptide repulsion (in order to reach the threshold concentration), and an acceptable tendency to penetrate into the bilayer. This mechanism is different from those proposed for peptides with long amino acid sequences. Such information is very useful from the medicinal chemistry point of view for the design of new small antimicrobial peptides.

AmberTools.

AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.

Substrate Recognition Properties from an Intermediate Structural State of the UreA Transporter.

Through a combination of comparative modeling, site-directed and classical random mutagenesis approaches, we previously identified critical residues for binding, recognition, and translocation of urea, and its inhibition by 2-thiourea and acetamide in the Aspergillus nidulans urea transporter, UreA. To deepen the structural characterization of UreA, we employed the artificial intelligence (AI) based AlphaFold2 (AF2) program. In this analysis, the resulting AF2 models lacked inward- and outward-facing cavities, suggesting a structural intermediate state of UreA. Moreover, the orientation of the W82, W84, N279, and T282 side chains showed a large variability, which in the case of W82 and W84, may operate as a gating mechanism in the ligand pathway. To test this hypothesis non-conservative and conservative substitutions of these amino acids were introduced, and binding and transport assessed for urea and its toxic analogue 2-thiourea, as well as binding of the structural analogue acetamide. As a result, residues W82, W84, N279, and T282 were implicated in substrate identification, selection, and translocation. Using molecular docking with Autodock Vina with flexible side chains, we corroborated the AF2 theoretical intermediate model, showing a remarkable correlation between docking scores and experimental affinities determined in wild-type and UreA mutants. The combination of AI-based modeling with classical docking, validated by comprehensive mutational analysis at the binding region, would suggest an unforeseen option to determine structural level details on a challenging family of proteins.

Wrapping Up Viruses at Multiscale Resolution: Optimizing PACKMOL and SIRAH Execution for Simulating the Zika Virus.

Simulating huge biomolecular complexes of million atoms at relevant biological time scales is becoming accessible to the broad scientific community. That proves to be crucial for urgent responses against emergent diseases in real time. Yet, there are still issues to sort regarding the system setup so that molecular dynamics (MD) simulations can be run in a simple and standard way. Here, we introduce an optimized pipeline for building and simulating enveloped virus-like particles (VLP). First, the membrane packing problem is tackled with new features and optimized options in PACKMOL. This allows preparing accurate membrane models of thousands of lipids in the context of a VLP within a few hours using a single CPU. Then, the assembly of the VLP system is done within the multiscale framework of the coarse-grained SIRAH force field. Finally, the equilibration protocol provides a system ready for production MD simulations within a few days on broadly accessible GPU resources. The pipeline is applied to study the Zika virus as a test case for large biomolecular systems. The VLP stabilizes at approximately 0.5 µs of MD simulation, reproducing correlations greater than 0.90 against experimental density maps from cryo-electron microscopy. Detailed structural analysis of the protein envelope also shows very good agreement in root-mean-square deviations and B-factors with the experimental data. The level of details attained shows for the first time a possible role for anionic phospholipids in stabilizing the envelope. Combining an efficient and reliable setup procedure with an accurate coarse-grained force field provides a valuable pipeline for simulating arbitrary viral systems or subcellular compartments, paving the way toward whole-cell simulations.

Post-Translational Modifications at the Coarse-Grained Level with the SIRAH Force Field.

Post-translational modifications (PTMs) on proteins significantly enlarge the physicochemical diversity present in biological macromolecules, altering function, localization, and interactions. Despite their critical role in regulating cellular processes, theoretical methods are not yet fully capable of coping with this diversity. These limitations are particularly more marked for coarse-grained (CG) models, in which comprehensive and self-consistent parametrizations are less frequent. Here we present a set of topologies and interaction parameters for the most common PTMs, fully compatible with the SIRAH force field. The PTMs introduced here reach the same level of structural description of the existing SIRAH force field, expanding the chemical spectrum with promising applications in dynamical protein-protein interactions in large and complex cellular environments.

Coarse-Grained Parameters for Divalent Cations within the SIRAH Force Field.

Although molecular dynamics simulations allow for the study of interactions among virtually all biomolecular entities, metal ions still pose significant challenges in achieving an accurate structural and dynamical description of many biological assemblies, particularly to coarse-grained (CG) models. Although the reduced computational cost of CG methods often makes them the technique of choice for the study of large biomolecular systems, the parameterization of metal ions is still very crude or not available for the vast majority of CG force fields. Here, we show that incorporating statistical data retrieved from the Protein Data Bank (PDB) to set specific Lennard-Jones interactions can produce structurally accurate CG molecular dynamics simulations using the SIRAH force field. We provide a set of interaction parameters for calcium, magnesium, and zinc ions, which cover more than 80% of the metal-bound structures reported in the PDB. Simulations performed on several proteins and DNA systems show that it is possible to preclude the use of topological constraints by modifying specific Lennard-Jones interactions.

Dimerization confers increased stability to nucleases in 5' halves from glycine and glutamic acid tRNAs.

We have previously shown that 5' halves from tRNAGlyGCC and tRNAGluCUC are the most enriched small RNAs in the extracellular space of human cell lines, and especially in the non-vesicular fraction. Extracellular RNAs are believed to require protection by either encapsulation in vesicles or ribonucleoprotein complex formation. However, deproteinization of non-vesicular tRNA halves does not affect their retention in size-exclusion chromatography. Thus, we considered alternative explanations for their extracellular stability. In-silico analysis of the sequence of these tRNA-derived fragments showed that tRNAGly 5' halves can form homodimers or heterodimers with tRNAGlu 5' halves. This capacity is virtually unique to glycine tRNAs. By analyzing synthetic oligonucleotides by size exclusion chromatography, we provide evidence that dimerization is possible in vitro. tRNA halves with single point substitutions preventing dimerization are degraded faster both in controlled nuclease digestion assays and after transfection in cells, showing that dimerization can stabilize tRNA halves against the action of cellular nucleases. Finally, we give evidence supporting dimerization of endogenous tRNAGlyGCC 5' halves inside cells. Considering recent reports have shown that 5' tRNA halves from Ala and Cys can form tetramers, our results highlight RNA intermolecular structures as a new layer of complexity in the biology of tRNA-derived fragments.

Mechanistic and biological characterisation of novel N5-substituted paullones targeting the biosynthesis of trypanothione in Leishmania.

Trypanothione synthetase (TryS) produces N1,N8-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 N5-substituted paullone analogues with anti-TryS activity. Several of the new derivatives retained submicromolar IC50 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 (EC50 0.5-10 µM) 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.

CGMD Platform: Integrated Web Servers for the Preparation, Running, and Analysis of Coarse-Grained Molecular Dynamics Simulations.

Advances in coarse-grained molecular dynamics (CGMD) simulations have extended the use of computational studies on biological macromolecules and their complexes, as well as the interactions of membrane protein and lipid complexes at a reduced level of representation, allowing longer and larger molecular dynamics simulations. Here, we present a computational platform dedicated to the preparation, running, and analysis of CGMD simulations. The platform is built on a completely revisited version of our Martini coarsE gRained MembrAne proteIn Dynamics (MERMAID) web server, and it integrates this with other three dedicated services. In its current version, the platform expands the existing implementation of the Martini force field for membrane proteins to also allow the simulation of soluble proteins using the Martini and the SIRAH force fields. Moreover, it offers an automated protocol for carrying out the backmapping of the coarse-grained description of the system into an atomistic one.

Imaging cAMP nanodomains in the heart.

3'-5'-cyclic adenosine monophosphate (cAMP) is a ubiquitous second messenger that modulates multiple cellular functions. It is now well established that cAMP can mediate a plethora of functional effects via a complex system of local regulatory mechanisms that result in compartmentalized signalling. The use of fluorescent probes to monitor cAMP in intact, living cells have been instrumental in furthering our appreciation of this ancestral and ubiquitous pathway and unexpected details of the nano-architecture of the cAMP signalling network are starting to emerge. Recent evidence shows that sympathetic control of cardiac contraction and relaxation is achieved via generation of multiple, distinct pools of cAMP that lead to differential phosphorylation of target proteins localized only tens of nanometres apart. The specific local control at these nanodomains is enabled by a distinct signalosome where effectors, targets, and regulators of the cAMP signal are clustered. In this review, we focus on recent advances using targeted fluorescent reporters for cAMP and how they have contributed to our current understanding of nanodomain cAMP signalling in the heart. We briefly discuss how this information can be exploited to design novel therapies and we highlight some of the questions that remain unanswered.

Multiscale modelization in a small virus: Mechanism of proton channeling and its role in triggering capsid disassembly.

In this work, we assess a previously advanced hypothesis that predicts the existence of ion channels in the capsid of small and non-enveloped icosahedral viruses. With this purpose we examine Triatoma Virus (TrV) as a case study. This virus has a stable capsid under highly acidic conditions but disassembles and releases the genome in alkaline environments. Our calculations range from a subtle sub-atomic proton interchange to the dismantling of a large-scale system representing several million of atoms. Our results provide structure-based explanations for the three roles played by the capsid to enable genome release. First, we observe, for the first time, the formation of a hydrophobic gate in the cavity along the five-fold axis of the wild-type virus capsid, which can be disrupted by an ion located in the pore. Second, the channel enables protons to permeate the capsid through a unidirectional Grotthuss-like mechanism, which is the most likely process through which the capsid senses pH. Finally, assuming that the proton leak promotes a charge imbalance in the interior of the capsid, we model an internal pressure that forces shell cracking using coarse-grained simulations. Although qualitatively, this last step could represent the mechanism of capsid opening that allows RNA release. All of our calculations are in agreement with current experimental data obtained using TrV and describe a cascade of events that could explain the destabilization and disassembly of similar icosahedral viruses.

Exploring DNA dynamics within oligonucleosomes with coarse-grained simulations: SIRAH force field extension for protein-DNA complexes.

Describing the regulation of chromatin segments by protein recognition events constitute a major goal in biology and biotechnology. Despite astonishing experimental developments, achieving nearly atomistic spatial/temporal resolution on such macromolecular systems remains a big challenge owing to the intrinsic flexibility of large biological assemblies. Although computer simulations have become a reliable complement to experimental techniques, computational cost limits their routine applications to relatively small systems. However, the development of accurate and cost-effective coarse-grained (CG) models helps to bridge the gap between molecular dynamics simulations and biologically relevant scales. Performing an exhaustive search on a set of well-resolved crystallographic protein-DNA complexes, we introduced improvements on the CG SIRAH force field to describe protein-DNA interfaces. Modifications were validated against a set of non redundant structures and applied to the simulation of the longest DNA segment in complex with proteins that has been crystallized to date, i.e. a tetranucleosome. Multimicrosecond simulation of this small chromatin segment evidences a large mobility of the external DNA filaments, which is consistent with results from FRET experiments in solution. Moreover, we found that the sub-microsecond dynamics of DNA is strongly modulated by the quaternary structure, partially overcoming the intrinsic dynamics dictated by the primary structure.

Dissecting the conformational determinants of chitosan and chitlac oligomers.

Chitosan and its highly hydrophilic 1-deoxy-lactit-1-yl derivative (Chitlac) are polysaccharides with increasing biomedical applications. Aimed to unravel their conformational properties we have performed a series of molecular dynamics simulations of Chitosan/Chitlac decamers, exploring different degrees of substitution (DS) of lactitol side chains. At low DS, two conformational regions with different populations are visited, while for DS ≥ 20% the oligomers remain mostly linear and only one main region of the glycosidic angles is sampled. These conformers are (locally) characterized by extended helical "propensities". Helical conformations 32 and 21, by far the most abundant, only develop in the main region. The accessible conformational space is clearly enlarged at high ionic strength, evidencing also a new region accessible to the glycosidic angles, with short and frequent interchange between regions. Simulations of neutral decamers share these features, pointing to a central role of electrostatic repulsion between charged moieties. These interactions seem to determine the conformational behavior of the chitosan backbone, with no evident influence of H-bond interactions. Finally, it is also shown that increasing temperature only slightly enlarges the available conformational space, but certainly without signs of a temperature-induced conformational transition.

Revisiting the human polypeptide GalNAc-T1 and T13 paralogs.

Polypeptide GalNAc-transferases (GalNAc-Ts) constitute a family of 20 human glycosyltransferases (comprising 9 subfamilies), which initiate mucin-type O-glycosylation. The O-glycoproteome is thought to be differentially regulated via the different substrate specificities and expression patterns of each GalNAc-T isoforms. Here, we present a comprehensive in vitro analysis of the peptide substrate specificity of GalNAc-T13, showing that it essentially overlaps with the ubiquitous expressed GalNAc-T1 isoform found in the same subfamily as T13. We have also identified and partially characterized nine splice variants of GalNAc-T13, which add further complexity to the GalNAc-T family. Two variants with changes in their lectin domains were characterized by in vitro glycosylation assays, and one (Δ39Ex9) was inactive while the second one (Ex10b) had essentially unaltered activity. We used reverse transcription-polymerase chain reaction analysis of human neuroblastoma cell lines, normal brain and a small panel of neuroblastoma tumors to demonstrate that several splice variants (Ex10b, ΔEx9, ΔEx2-7 and ΔEx6/8-39bpEx9) were highly expressed in tumor cell lines compared with normal brain, although the functional implications remain to be unveiled. In summary, the GalNAc-T13 isoform is predicted to function similarly to GalNAc-T1 against peptide substrates in vivo, in contrast to a prior report, but is unique by being selectively expressed in the brain.

The ataxia related G1107D mutation of the plasma membrane Ca2+ ATPase isoform 3 affects its interplay with calmodulin and the autoinhibition process.

The plasma membrane Ca2+ ATPases (PMCA pumps) have a long, cytosolic C-terminal regulatory region where a calmodulin-binding domain (CaM-BD) is located. Under basal conditions (low Ca2+), the C-terminal tail of the pump interacts with autoinhibitory sites proximal to the active center of the enzyme. In activating conditions (i.e., high Ca2+), Ca2+-bound CaM displaces the C-terminal tail from the autoinhibitory sites, restoring activity. We have recently identified a G1107D replacement within the CaM-BD of isoform 3 of the PMCA pump in a family affected by X-linked congenital cerebellar ataxia. Here, we investigate the effects of the G1107D replacement on the interplay of the mutated CaM-BD with both CaM and the pump core, by combining computational, biochemical and functional approaches. We provide evidence that the affinity of the isolated mutated CaM-BD for CaM is significantly reduced with respect to the wild type (wt) counterpart, and that the ability of CaM to activate the pump in vitro is thus decreased. Multiscale simulations support the conclusions on the detrimental effect of the mutation, indicating reduced stability of the CaM binding. We further show that the G1107D replacement impairs the autoinhibition mechanism of the PMCA3 pump as well, as the introduction of a negative charge perturbs the contacts between the CaM-BD and the pump core. Thus, the mutation affects both the ability of the pump to optimally transport Ca2+ in the activated state, and the autoinhibition mechanism in its resting state.

SIRAH tools: mapping, backmapping and visualization of coarse-grained models.

UNLABELLED: Coarse-grained (CG) models reduce the cost of molecular dynamics simulations keeping the essence of molecular interactions. Still, the diversity of CG representations (sizes, connectivity, naming, etc.) hampers the handling and visualization of such models. SIRAH Tools comprises a set of utilities to convert all-atoms coordinates to arbitrary residue-based CG schemes, write GROMACS' topological information at any resolution into PSF format and a VMD plugin to visualize, analyze and retrieve pseudo-atomistic information from CG trajectories performed with the SIRAH force field. These tools facilitate the use of intricate CG force fields outside the small developer's community. AVAILABILITY AND IMPLEMENTATION: Different utilities of SIRAH Tools are written in Perl, Tcl, or R. Documentation and source codes are freely distributed at http://www.sirahff.com CONTACT: : mmachado@pasteur.edu.uy or spantano@pasteur.edu.uy.

Comparative analysis reveals amino acids critical for anticancer activity of peptide CIGB-552.

Because of resistance development by cancer cells against current anticancer drugs, there is a considerable interest in developing novel antitumor agents. We have previously demonstrated that CIGB-552, a novel cell-penetrating synthetic peptide, was effective in reducing tumor size and increasing lifespan in tumor-bearing mice. Studies of protein-peptide interactions have shown that COMMD1 protein is a major mediator of CIGB-552 antitumor activity. Furthermore, a typical serine-protease degradation pattern for CIGB-552 in BALB/c mice serum was identified, yielding peptides which differ from CIGB-552 in size and physical properties. In the present study, we show the results obtained from a comparative analysis between CIGB-552 and its main metabolites regarding physicochemical properties, cellular internalization, and their capability to elicit apoptosis in MCF-7 cells. None of the analyzed metabolites proved to be as effective as CIGB-552 in promoting apoptosis in MCF-7. Taking into account these results, it seemed important to examine their cell-penetrating capacity and interaction with COMMD1. We show that internalization, a lipid binding-dependent process, is impaired as well as metabolite-COMMD1 interaction, key component of the apoptotic mechanism. Altogether, our results suggest that features conferred by the amino acid sequence are decisive for CIGB-552 biological activity, turning it into the minimal functional unit. Copyright © 2016 European Peptide Society and John Wiley & Sons, Ltd.