A Practical Guide to All-Atom and Coarse-Grained Molecular Dynamics Simulations Using Amber and Gromacs: A Case Study of Disulfide-Bond Impact on the Intrinsically Disordered Amyloid Beta.

Intrinsically disordered proteins (IDPs) pose challenges to conventional experimental techniques due to their large-scale conformational fluctuations and transient structural elements. This work presents computational methods for studying IDPs at various resolutions using the Amber and Gromacs packages with both all-atom (Amber ff19SB with the OPC water model) and coarse-grained (Martini 3 and SIRAH) approaches. The effectiveness of these methodologies is demonstrated by examining the monomeric form of amyloid-ß (Aß42), an IDP, with and without disulfide bonds at different resolutions. Our results clearly show that the addition of a disulfide bond decreases the ß-content of Aß42; however, it increases the tendency of the monomeric Aß42 to form fibril-like conformations, explaining the various aggregation rates observed in experiments. Moreover, analysis of the monomeric Aß42 compactness, secondary structure content, and comparison between calculated and experimental chemical shifts demonstrates that all three methods provide a reasonable choice to study IDPs; however, coarse-grained approaches may lack some atomistic details, such as secondary structure recognition, due to the simplifications used. In general, this study not only explains the role of disulfide bonds in Aß42 but also provides a step-by-step protocol for setting up, conducting, and analyzing molecular dynamics (MD) simulations, which is adaptable for studying other biomacromolecules, including folded and disordered proteins and peptides.

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.

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.

Investigation of the Structure of Full-Length Tau Proteins with Coarse-Grained and All-Atom Molecular Dynamics Simulations.

Tau proteins not only have many important biological functions but also are associated with several neurodegenerative diseases, such as Parkinson's disease and Alzheimer's disease (AD). However, it is still a challenge to identify the atomic structure of full-length tau proteins due to their lengthy and disordered characteristics and the factor that there are no crystal structures of full-length tau proteins available. We performed multi- and large-scale molecular dynamics simulations of the full-length tau monomer (the 2N4R isoform and 441 residues) in aqueous solution under biological conditions with coarse-grained and all-atom force fields. The obtained atomic structures produced radii of gyration and chemical shifts that are in excellent agreement with those of experiment. The generated monomer structure ensemble would be very useful for further studying the oligomerization mechanism and discovering tau oligomerization inhibitors, which are important events in AD drug development.

Early stages in Aß1-42 spontaneous aggregation: An unbiased dataset from coarse-grained molecular dynamics simulations.

The small soluble aggregates of Aß1-42 are broadly documented as potential targets for the development of new compounds with the capacity to inhibit the early stages of Alzheimer´s disease. Nevertheless, Aß1-42 peptides show an intrinsically disordered character with a high propensity for aggregation, which complicates the identification of conserved structural patterns. Because of this, experimental techniques find substantial difficulties in the characterization of such soluble oligomers. Theoretical techniques, such as molecular dynamics (MD) simulations, provide a possible workaround for this problem. However, the computational cost associated with comprehensively sampling the vast conformational space accessible to these peptides might become prohibitive. In this sense, coarse-grained (CG) simulations can effectively overcome that hurdle at a fraction of the computational cost. In this dataset, we furnish an extensive collection of Aß1-42 peptides in dimeric conformation generated with the SIRAH force field for CG MD simulation. It comprises 25 independent trajectories in .xtc (gromacs) format of Aß1-42 couples of peptides that evolve towards dimeric states along eleven µs-long unbiased simulations. Thanks to the backmapping capabilities of our force field, pseudo atomistic coordinates can be straightforwardly recovered from MD trajectories reported here and analyzed with popular molecular editing programs. This set of simulations performed at room conditions and physiological salt concentrations may furnish a complete collection of inter-peptide interfaces that can be used in high-throughput docking or as new starting states for peptide oligomerization seeding of Aß1-42 dimerization.

The stressed life of a lipid in the Zika virus membrane.

Protein-lipid interactions modulate a plethora of physiopathologic processes and have been the subject of countless studies. However, these kinds of interactions in the context of viral envelopes have remained relatively unexplored, partially because the intrinsically small dimensions of the molecular systems escape to the current resolution of experimental techniques. However, coarse-grained and multiscale simulations may fill that niche, providing nearly atomistic resolution at an affordable computational price. Here we use multiscale simulations to characterize the lipid-protein interactions in the envelope of the Zika Virus, a prominent member of the Flavivirus genus. Comparisons between the viral envelope and simpler molecular systems indicate that the viral membrane is under extreme pressures and asymmetric forces. Furthermore, the dense net of protein-protein contacts established by the envelope proteins creates poorly solvated regions that destabilize the external leaflet leading to a decoupled dynamics between both membrane layers. These findings lead to the idea that the Flaviviral membrane may store a significant amount of elastic energy, playing an active role in the membrane fusion process.

Investigation of the pH-dependent aggregation mechanisms of GCSF using low resolution protein characterization techniques and advanced molecular dynamics simulations.

Granulocyte-colony stimulating factor (GCSF) is a widely used therapeutic protein to treat neutropenia. GCSF has an increased propensity to aggregate if the pH is increased above 5.0. Although GCSF is very well experimentally characterized, the exact pH-dependent aggregation mechanism of GCSF is still under debate. This study aimed to model the complex pH-dependent aggregation behavior of GCSF using state-of-the-art simulation techniques. The conformational stability of GCSF was investigated by performing metadynamics simulations, while the protein-protein interactions were investigated using coarse-grained (CG) simulations of multiple GCSF monomers. The CG simulations were directly compared with small-angle X-ray (SAXS) data. The metadynamics simulations demonstrated that the orientations of Trp residues in GCSF are dependent on pH. The conformational change of Trp residues is due to the loss of Trp-His interactions at the physiological pH, which in turn may increase protein flexibility. The helical structure of GCSF was not affected by the pH conditions of the simulations. Our CG simulations indicate that at pH 4.0, the colloidal stability may be more important than the conformational stability of GCSF. The electrostatic potential surface and CG simulations suggested that the basic residues are mainly responsible for colloidal stability as deprotonation of these residues causes a reduction of the highly positively charged electrostatic barrier close to the aggregation-prone long loop regions.

CUTie2: The Attack of the Cyclic Nucleotide Sensor Clones.

The detection of small molecules in living cells using genetically encoded FRET sensors has revolutionized our understanding of signaling pathways at the sub-cellular level. However, engineering fluorescent proteins and specific binding domains to create new sensors remains challenging because of the difficulties associated with the large size of the polypeptides involved, and their intrinsically huge conformational variability. Indeed, FRET sensors' design still relies on vague structural notions, and trial and error combinations of linkers and protein modules. We recently designed a FRET sensor for the second messenger cAMP named CUTie (Cyclic nucleotide Universal Tag for imaging experiments), which granted sub-micrometer resolution in living cells. Here we apply a combination of sequence/structure analysis to produce a new-generation FRET sensor for the second messenger cGMP based on Protein kinase G I (PKGI), which we named CUTie2. Coarse-grained molecular dynamics simulations achieved an exhaustive sampling of the relevant spatio-temporal coordinates providing a quasi-quantitative prediction of the FRET efficiency, as confirmed by in vitro experiments. Moreover, biochemical characterization showed that the cGMP binding module maintains virtually the same affinity and selectivity for its ligand thant the full-length protein. The computational approach proposed here is easily generalizable to other allosteric protein modules, providing a cost effective-strategy for the custom design of FRET sensors.

A homogeneous dataset of polyglutamine and glutamine rich aggregating peptides simulations.

This dataset contains a collection of molecular dynamics (MD) simulations of polyglutamine (polyQ) and glutamine-rich (Q-rich) peptides in the multi-microsecond timescale. Primary data from coarse-grained simulations performed using the SIRAH force field has been processed to provide fully atomistic coordinates. The dataset encloses MD trajectories of polyQs of 4 (Q4), 11 (Q11), and 36 (Q36) amino acids long. In the case of Q11, simulations in presence of Q5 and QEQQQ peptides, which modulate aggregation, are also included. The dataset also comprises MD trajectories of the gliadin related p31-43 peptide, and Insulin's C-peptide at pH=7 and pH=3.2, which constitute examples of Q-rich and Q-poor aggregating peptides. The dataset grants molecular insights on the role of glutamines in spontaneous and unbiased ab-initio aggregation of a series of peptides using a homogeneous set of simulations [1]. The trajectory files are provided in Protein Data Bank (PDB) format containing the Cartesian coordinates of all heavy atoms in the aggregating peptides. Further analyses of the trajectories can be performed directly using any molecular visualization/analysis software suites.

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.

Cues to Opening Mechanisms From in Silico Electric Field Excitation of Cx26 Hemichannel and in Vitro Mutagenesis Studies in HeLa Transfectans.

Connexin channels play numerous essential roles in virtually every organ by mediating solute exchange between adjacent cells, or between cytoplasm and extracellular milieu. Our understanding of the structure-function relationship of connexin channels relies on X-ray crystallographic data for human connexin 26 (hCx26) intercellular gap junction channels. Comparison of experimental data and molecular dynamics simulations suggests that the published structures represent neither fully-open nor closed configurations. To facilitate the search for alternative stable configurations, we developed a coarse grained (CG) molecular model of the hCx26 hemichannel and studied its responses to external electric fields. When challenged by a field of 0.06 V/nm, the hemichannel relaxed toward a novel configuration characterized by a widened pore and an increased bending of the second transmembrane helix (TM2) at the level of the conserved Pro87. A point mutation that inhibited such transition in our simulations impeded hemichannel opening in electrophysiology and dye uptake experiments conducted on HeLa tranfectants. These results suggest that the hCx26 hemichannel uses a global degree of freedom to transit between different configuration states, which may be shared among the whole connexin family.