GROMACS—the road ahead

Understanding Genetic Risks: Computational Exploration of Human β-Synuclein nsSNPs and their Potential Impact on Structural Alteration

Synucleins are pivotal in neurodegenerative conditions. Beta-synuclein (β-synuclein) is part of the synuclein protein family alongside alpha-synuclein (α-synuclein) and gamma-synuclein (γ-synuclein). These proteins, found mainly in brain tissue and cancers, are soluble and unstructured. β-synuclein shares significant similarity with α-synuclein, especially in their N-terminus, with a 90% match. However, their aggregation tendencies differ significantly. While α-synuclein aggregation is believed to be counteracted by β-synuclein, which occurs in conditions like Parkinson's disease, β-synuclein may counteract α-synuclein's toxic effects on the nervous system, offering potential treatment for neurodegenerative diseases. Under normal circumstances, β-synuclein may guard against disease by interacting with α-synuclein. Yet, in pathological environments with heightened levels or toxic substances, it might contribute to disease. Our research aims to explore potential harmful mutations in the β-synuclein using computational tools to predict their destabilizing impact on protein structure. Consensus analysis revealed rs1207608813 (A63P), rs1340051870 (S72F), and rs1581178262 (G36C) as deleterious. These findings highlight the intricate relationship between nsSNPs and protein function, shedding light on their potential implications in disease pathways. Understanding the structural consequences of nsSNPs is crucial for elucidating their role in pathogenesis and developing targeted therapeutic interventions. Our results offer a robust computational framework for identifying neurodegenerative disorder-related mutations from SNP datasets, potentially reducing the costs associated with experimental characterization.

Phospholipid tail asymmetry allows cellular adaptation to anoxic environments

Membrane biophysical properties are critical to cell fitness and depend on unsaturated phospholipid acyl tails. These can only be produced in aerobic environments since eukaryotic desaturases require molecular oxygen. This raises the question of how cells maintain bilayer properties in anoxic environments. Using advanced microscopy, molecular dynamics simulations, and lipidomics by mass spectrometry we demonstrated the existence of an alternative pathway to regulate membrane fluidity that exploits phospholipid acyl tail length asymmetry, replacing unsaturated species in the membrane lipidome. We show that the fission yeast, Schizosaccharomyces japonicus, which can grow in aerobic and anaerobic conditions, is capable of utilizing this strategy, whereas its sister species, the well-known model organism Schizosaccharomyces pombe, cannot. The incorporation of asymmetric-tailed phospholipids might be a general adaptation to hypoxic environmental niches.

On the Structure and Flip-flop of Free Docosahexaenoic Acid in a Model Human Brain Membrane

Among the omega-3 fatty acids, docosahexaenoic acid (DHA, sn22:6) is particularly vital in human brain cell membranes. There is considerable interest in DHA because low-level DHA has been associated with declined cognitive function and poor visual acuity. In this work, atomistic molecular dynamics simulations were used to investigate the effects of free protonated DHA (DHAP) in molar fractions of 0, 17, 30, and 38% in a realistic model of a healthy brain cell membrane comprising 26 lipid types. Numerous flip-flop events of DHAP were observed and categorized as successful or aborted. Novel use of the machine learning technique, density-based spatial clustering of applications with noise (DBSCAN), effectively identified flip-flop events by way of clustering. Our data show that increasing amounts of DHAP in the membrane disorder the bilayer packing, fluidize the membrane, and increase the rates of successful flip-flop from k = 0.2 μs^-1 (17% DHAP) to 0.8 μs^-1 (30% DHAP) and to 1.3 μs^-1 (38% DHAP). In addition, we also provided a detailed understanding of the flip-flop mechanism of DHAP across this complex membrane. Interestingly, we noted the role of hydrogen bonds in two distinct coordinated flip-flop phenomena between two DHAP molecules: double flip-flop and assisted flip-flop. Understanding the effects of various concentrations of DHAP on the dynamics within a lipid membrane and the resulting structural properties of the membrane are important when considering the use of DHAP as a dietary supplement or as a potential therapeutic.

Asymmetric glycerophospholipids impart distinctive biophysical properties to lipid bilayers

Phospholipids are a diverse group of biomolecules consisting of a hydrophilic headgroup and two hydrophobic acyl tails. The nature of the head and length and saturation of the acyl tails are important for defining the biophysical properties of lipid bilayers. It has recently been shown that the membranes of certain yeast species contain high levels of unusual asymmetric phospholipids consisting of one long and one medium-chain acyl moiety, a configuration not common in mammalian cells or other well-studied model yeast species. This raises the possibility that structurally asymmetric glycerophospholipids impart distinctive biophysical properties to the yeast membranes. Previously, it has been shown that lipids with asymmetric length tails form a mixed interdigitated gel phase and exhibit unusual endotherm behavior upon heating and cooling. Here, however, we address physiologically relevant temperature conditions and, using atomistic molecular dynamics simulations and environmentally sensitive fluorescent membrane probes, characterize key biophysical parameters (such as lipid packing, diffusion coefficient, membrane thickness, and area per lipid) in membranes composed of both length-asymmetric glycerophospholipids and ergosterol. Interestingly, we show that saturated but asymmetric glycerophospholipids maintain membrane lipid order across a wide range of temperatures. We also show that these asymmetric lipids can substiture of unsaturated symmetric lipids in the phase behaviour of ternary lipid bilayers. This may allow cells to maintain membrane fluidity, even in environments that lack oxygen, which is required for the synthesis of unsaturated lipids and sterols.

Laurdan and Di-4-ANEPPDHQ Influence the Properties of Lipid Membranes: A Classical Molecular Dynamics and Fluorescence Study

Environmentally sensitive (ES) dyes have been used for many decades to study the lipid order of cell membranes, as different lipid phases play a crucial role in a wide variety of cell processes. Yet, the understanding of how ES dyes behave, interact, and affect membranes at the atomistic scale is lacking, partially due to the lack of molecular dynamics (MD) models of these dyes. Here, we present ground- and excited-state MD models of commonly used ES dyes, Laurdan and di-4-ANEPPDHQ, and use MD simulations to study the behavior of these dyes in a disordered and an ordered membrane. We also investigate the effect that these two dyes have on the hydration and lipid order of the membranes, where we see a significant effect on the hydration of lipids proximal to the dyes. These findings are combined with experimental fluorescence experiments of ordered and disordered vesicles and live HeLa cells stained by the aforementioned dyes, where the generalized polarization (GP) values were measured at different concentrations of the dyes. We observe a small but significant decrease of GP at higher Laurdan concentrations in vesicles, while the same effect is not observed in cell membranes. The opposite effect is observed with di-4-ANEPPDHQ where no significant change in GP is seen for vesicles but a very substantial and significant decrease is seen in cell membranes. Together, our results show the profound effect that ES dyes have on membranes, and the presented MD models will be important for further understanding of these effects.

A potential for molecular simulation of compounds with linear moieties

The harmonic angle bending potential is used in many force fields for (bio)molecular simulation. The force associated with this potential is discontinuous at angles close to 180°, which can lead to numeric instabilities. Angle bending of linear groups, such as alkynes or nitriles, or linear molecules, such as carbon dioxide, can be treated by a simple harmonic potential if we describe the fluctuations as a deviation from a reference position of the central atom, the position of which is determined by the flanking atoms. The force constant for the linear angle potential can be derived analytically from the corresponding force constant in the traditional potential. The new potential is tested on the properties of alkynes, nitriles, and carbon dioxide. We find that the angles of the linear groups remain about 2° closer to 180° using the new potential. The bond and angle force constants for carbon dioxide were tuned to reproduce the experimentally determined frequencies. An interesting finding was that simulations of liquid carbon dioxide under pressure with the new flexible model were stable only when explicitly modeling the long-range Lennard-Jones (LJ) interactions due to the very long-range nature of the LJ interactions (>1.7 nm). In the other tested liquids, we find that a Lennard-Jones cutoff of 1.1 nm yields similar results as the particle mesh Ewald algorithm for LJ interactions. Algorithmic factors influencing the stability of liquid simulations are discussed as well. Finally, we demonstrate that the linear angle potential can be used in free energy perturbation calculations.

Time-Resolved Fluorescence Anisotropy of a Molecular Rotor Resolves Microscopic Viscosity Parameters in Complex Environments

Understanding viscosity in complex environments remains a largely unanswered question despite its importance in determining reaction rates in vivo. Here, time-resolved fluorescence anisotropy imaging (TR-FAIM) is combined with fluorescent molecular rotors (FMRs) to simultaneously determine two non-equivalent viscosity-related parameters in complex heterogeneous environments. The parameters, FMR rotational correlation time and lifetime, are extracted from fluorescence anisotropy decays, which in heterogeneous environments show dip-and-rise behavior due to multiple dye populations. Decays of this kind are found both in artificially constructed adiposomes and in live cell lipid droplet organelles. Molecular dynamics simulations are used to assign each population to nano-environments within the lipid systems. The less viscous population corresponds to the state showing an average 25° tilt to the lipid membrane normal, and the more viscous population to the state showing an average 55° tilt. This combined experimental and simulation approach enables a comprehensive description of the FMR probe behavior within viscous nano-environments in complex, biological systems.

Double-buffered, heterogeneous CPU + GPU integral digestion algorithm for single-excitation calculations involving a large number of excited states

The most widely used quantum-chemical models for excited states are single-excitation theories, a category that includes configuration interaction with single substitutions, time-dependent density functional theory, and also a recently developed ab initio exciton model. When a large number of excited states are desired, these calculations incur a significant bottleneck in the "digestion" step in which two-electron integrals are contracted with density or density-like matrices. We present an implementation that moves this step onto graphical processing units (GPUs), and introduce a double-buffer scheme that minimizes latency by computing integrals on the central processing units (CPUs) concurrently with their digestion on the GPUs. An automatic code generation scheme simplifies the implementation of high-performance GPU kernels. For the exciton model, which requires separate excited-state calculations on each electronically coupled chromophore, the heterogeneous implementation described here results in speedups of 2-6× versus a CPU-only implementation. For traditional time-dependent density functional theory calculations, we obtain speedups of up to 5× when a large number of excited states is computed. © 2018 Wiley Periodicals, Inc.

Two-Dimensional Folding of Polypeptides into Molecular Nanostructures at Surfaces

Herein we report the fabrication of molecular nanostructures on surfaces via two-dimensional (2D) folding of the nine amino acid peptide bradykinin. Soft-landing electrospray ion beam deposition in conjunction with high-resolution imaging by scanning tunneling microscopy is used to fabricate and investigate the molecular nanostructures. Subnanometer resolved images evidence the large conformational freedom of the molecules if thermal motion is inhibited and the formation of stable uniform dimers of only one specific conformation when diffusion can take place. Molecular dynamics modeling supported by density functional theory calculations give atomically precise insight into the induced-fit binding scheme when the folded dimer is formed. In the absence of solvent, we find a hierarchy of binding strength from polar to nonpolar, manifested in an inverted polar-nonpolar segregation which suppresses unspecific interactions at the rim of the nanostructure. The demonstrated 2D-folding scheme resembles many key properties of its native 3D counterpart and shows that functional, molecular nanostructures on surfaces fabricated by folding could be just as versatile and specific.

Emergent stem cell homeostasis in the C. elegans germline is revealed by hybrid modeling

The establishment of homeostasis among cell growth, differentiation, and apoptosis is of key importance for organogenesis. Stem cells respond to temporally and spatially regulated signals by switching from mitotic proliferation to asymmetric cell division and differentiation. Executable computer models of signaling pathways can accurately reproduce a wide range of biological phenomena by reducing detailed chemical kinetics to a discrete, finite form. Moreover, coordinated cell movements and physical cell-cell interactions are required for the formation of three-dimensional structures that are the building blocks of organs. To capture all these aspects, we have developed a hybrid executable/physical model describing stem cell proliferation, differentiation, and homeostasis in the Caenorhabditis elegans germline. Using this hybrid model, we are able to track cell lineages and dynamic cell movements during germ cell differentiation. We further show how apoptosis regulates germ cell homeostasis in the gonad, and propose a role for intercellular pressure in developmental control. Finally, we use the model to demonstrate how an executable model can be developed from the hybrid system, identifying a mechanism that ensures invariance in fate patterns in the presence of instability.

Identification of Distinct Conformations of the Angiotensin-II Type 1 Receptor Associated with the Gq/11 Protein Pathway and the β-Arrestin Pathway Using Molecular Dynamics Simulations

Biased signaling represents the ability of G protein-coupled receptors to engage distinct pathways with various efficacies depending on the ligand used or on mutations in the receptor. The angiotensin-II type 1 (AT1) receptor, a prototypical class A G protein-coupled receptor, can activate various effectors upon stimulation with the endogenous ligand angiotensin-II (AngII), including the Gq/11 protein and β-arrestins. It is believed that the activation of those two pathways can be associated with distinct conformations of the AT1 receptor. To verify this hypothesis, microseconds of molecular dynamics simulations were computed to explore the conformational landscape sampled by the WT-AT1 receptor, the N111G-AT1 receptor (constitutively active and biased for the Gq/11 pathway), and the D74N-AT1 receptor (biased for the β-arrestin1 and -2 pathways) in their apo-forms and in complex with AngII. The molecular dynamics simulations of the AngII-WT-AT1, N111G-AT1, and AngII-N111G-AT1 receptors revealed specific structural rearrangements compared with the initial and ground state of the receptor. Simulations of the D74N-AT1 receptor revealed that the mutation stabilizes the receptor in the initial ground state. The presence of AngII further stabilized the ground state of the D74N-AT1 receptor. The biased agonist [Sar(1),Ile(8)]AngII also showed a preference for the ground state of the WT-AT1 receptor compared with AngII. These results suggest that activation of the Gq/11 pathway is associated with a specific conformational transition stabilized by the agonist, whereas the activation of the β-arrestin pathway is linked to the stabilization of the ground state of the receptor.

A prospective overview of the essential requirements in molecular modeling for nanomedicine design

Nanotechnology has presented many new challenges and opportunities in the area of nanomedicine design. The issues related to nanoconjugation, nanosystem-mediated targeted drug delivery, transitional stability of nanovehicles, the integrity of drug transport, drug-delivery mechanisms and chemical structural design require a pre-estimated and determined course of assumptive actions with property and characteristic estimations for optimal nanomedicine design. Molecular modeling in nanomedicine encompasses these pre-estimations and predictions of pertinent design data via interactive computographic software. Recently, an increasing amount of research has been reported where specialized software is being developed and employed in an attempt to bridge the gap between drug discovery, materials science and biology. This review provides an assimilative and concise incursion into the current and future strategies of molecular-modeling applications in nanomedicine design and aims to describe the utilization of molecular models and theoretical-chemistry computographic techniques for expansive nanomedicine design and development.

Thermodynamic and solution state NMR characterization of the binding of secondary and conjugated bile acids to STARD5

STARD5 is a member of the STARD4 sub-family of START domain containing proteins specialized in the non-vesicular transport of lipids and sterols. We recently reported that STARD5 binds primary bile acids. Herein, we report on the biophysical and structural characterization of the binding of secondary and conjugated bile acids by STARD5 at physiological concentrations. We found that the absence of the 7α-OH group and its epimerization increase the affinity of secondary bile acids for STARD5. According to NMR titration and molecular modeling, the affinity depends mainly on the number and positions of the steroid ring hydroxyl groups and to a lesser extent on the presence or type of bile acid side-chain conjugation. Primary and secondary bile acids have different binding modes and display different positioning within the STARD5 binding pocket. The relative STARD5 affinity for the different bile acids studied is: DCA>LCA>CDCA>GDCA>TDCA>CA>UDCA. TCA and GCA do not bind significantly to STARD5. The impact of the ligand chemical structure on the thermodynamics of binding is discussed. The discovery of these new ligands suggests that STARD5 is involved in the cellular response elicited by bile acids and offers many entry points to decipher its physiological role.

Insights on the solubility of CO2 in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide from the microscopic point of view

Emissions of greenhouse gases due to human activities have been well documented as well as the effects on global warming resulting from it. Efforts to reduce greenhouse gases at the source are crucial to curb climate change, but due to insignificant economic incentives to reduce usage of fossil fuels, not a lot of progress has been made by this route. This necessitates additional measures to reduce the occurrence of greenhouse gases in the atmosphere. Here we used theoretical methods to study the solubility of carbon dioxide in ionic liquids (ILs) since sequestration of CO2 in ILs has been proposed as a possible technology for reducing the emissions of CO2 to the atmosphere. Ionic liquids form a class of solvents with melting temperatures below 100 °C and, due to very low vapor pressures, which are not volatile. We have performed molecular dynamics (MD) simulations of 1-ethyl-3-methylimidazolium (C2mim) bis(trifluoromethylsulfonyl)imide (Tf2N) and its mixtures with carbon dioxide in order to investigate the CO2 concentration effect on the CO2-cation and CO2-anion interactions. A systematic investigation of CO2 concentration effects on resulting equilibrium liquid structure, and the local environment of the ions is provided. The Quantum Theory of Atoms in Molecules (QTAIM) was used to determine the interaction energy for CO2-cation and CO2-anion complexes from uncorrelated structures derived from MD simulations. A spatial distribution function analysis demonstrates the specific interactions between CO2 and the ionic liquid. Our findings indicate that the total volume of the system increases with the CO2 concentration, with a molar volume of CO2 of about 0.038 L/mol, corresponding to liquid CO2 under a pressure of 100 bar. In other words, the IL effectively pressurizes the CO2 inside its matrix. The thermodynamics of CO2 solvation in C2 min-Tf2N were computed using free energy techniques, and the solubility of CO2 is found to be higher in this IL (-3.7 ± 1 kcal/mol) than in water (+0.2 kJ/mol), predominantly due to anion-CO2 interactions.

Novel Ultrathin Membranes Composed of Organic Ions

Until recently, construction of bilayers was an exclusive mission of nature. It requires careful choice of compounds, whose delicate interplay between head group attraction and chain repulsion engenders a truly unique balance over a narrow temperature range. We report the investigation of artificial bilayers composed of long-chained organic ions, such as dodecyltrimethylammonium (DMA(+)) and perfluorooctaonate (PFO(-)). Various ratios of DMA/PFO surfactants result in bilayers of different stability, thickness, area per molecule, and density profiles. In our quest for water filtration, we incorporated aquaporin protein into the DMA/PFO bilayer but did not observe sufficient stability of the system. We discuss further steps to utilize these surfactant bilayers as highly selective, salt-impermeable membranes.

GROMACS 4.5: a high-throughput and highly parallel open source molecular simulation toolkit

MOTIVATION

Molecular simulation has historically been a low-throughput technique, but faster computers and increasing amounts of genomic and structural data are changing this by enabling large-scale automated simulation of, for instance, many conformers or mutants of biomolecules with or without a range of ligands. At the same time, advances in performance and scaling now make it possible to model complex biomolecular interaction and function in a manner directly testable by experiment. These applications share a need for fast and efficient software that can be deployed on massive scale in clusters, web servers, distributed computing or cloud resources.

RESULTS

Here, we present a range of new simulation algorithms and features developed during the past 4 years, leading up to the GROMACS 4.5 software package. The software now automatically handles wide classes of biomolecules, such as proteins, nucleic acids and lipids, and comes with all commonly used force fields for these molecules built-in. GROMACS supports several implicit solvent models, as well as new free-energy algorithms, and the software now uses multithreading for efficient parallelization even on low-end systems, including windows-based workstations. Together with hand-tuned assembly kernels and state-of-the-art parallelization, this provides extremely high performance and cost efficiency for high-throughput as well as massively parallel simulations.

AVAILABILITY

GROMACS is an open source and free software available from http://www.gromacs.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Critical hydrogen bond formation for activation of the angiotensin II type 1 receptor

G protein-coupled receptors contain selectively important residues that play central roles in the conformational changes that occur during receptor activation. Asparagine 111 (N111(3.35)) is such a residue within the angiotensin II type 1 (AT(1)) receptor. Substitution of N111(3.35) for glycine leads to a constitutively active receptor, whereas substitution for tryptophan leads to an inactivable receptor. Here, we analyzed the AT(1) receptor and two mutants (N111G and N111W) by molecular dynamics simulations, which revealed a novel molecular switch involving the strictly conserved residue D74(2.50). Indeed, D74(2.50) forms a stable hydrogen bond (H-bond) with the residue in position 111(3.35) in the wild-type and the inactivable receptor. However, in the constitutively active mutant N111G-AT(1) receptor, residue D74 is reoriented to form a new H-bond with another strictly conserved residue, N46(1.50). When expressed in HEK293 cells, the mutant N46G-AT(1) receptor was poorly activable, although it retained a high binding affinity. Interestingly, the mutant N46G/N111G-AT(1) receptor was also inactivable. Molecular dynamics simulations also revealed the presence of a cluster of hydrophobic residues from transmembrane domains 2, 3, and 7 that appears to stabilize the inactive form of the receptor. Whereas this hydrophobic cluster and the H-bond between D74(2.50) and W111(3.35) are more stable in the inactivable N111W-AT(1) receptor, the mutant N111W/F77A-AT(1) receptor, designed to weaken the hydrophobic core, showed significant agonist-induced signaling. These results support the potential for the formation of an H-bond between residues D74(2.50) and N46(1.50) in the activation of the AT(1) receptor.

Force Field Benchmark of Organic Liquids: Density, Enthalpy of Vaporization, Heat Capacities, Surface Tension, Isothermal Compressibility, Volumetric Expansion Coefficient, and Dielectric Constant

The chemical composition of small organic molecules is often very similar to amino acid side chains or the bases in nucleic acids, and hence there is no a priori reason why a molecular mechanics force field could not describe both organic liquids and biomolecules with a single parameter set. Here, we devise a benchmark for force fields in order to test the ability of existing force fields to reproduce some key properties of organic liquids, namely, the density, enthalpy of vaporization, the surface tension, the heat capacity at constant volume and pressure, the isothermal compressibility, the volumetric expansion coefficient, and the static dielectric constant. Well over 1200 experimental measurements were used for comparison to the simulations of 146 organic liquids. Novel polynomial interpolations of the dielectric constant (32 molecules), heat capacity at constant pressure (three molecules), and the isothermal compressibility (53 molecules) as a function of the temperature have been made, based on experimental data, in order to be able to compare simulation results to them. To compute the heat capacities, we applied the two phase thermodynamics method (Lin et al. J. Chem. Phys.2003, 119, 11792), which allows one to compute thermodynamic properties on the basis of the density of states as derived from the velocity autocorrelation function. The method is implemented in a new utility within the GROMACS molecular simulation package, named g_dos, and a detailed exposé of the underlying equations is presented. The purpose of this work is to establish the state of the art of two popular force fields, OPLS/AA (all-atom optimized potential for liquid simulation) and GAFF (generalized Amber force field), to find common bottlenecks, i.e., particularly difficult molecules, and to serve as a reference point for future force field development. To make for a fair playing field, all molecules were evaluated with the same parameter settings, such as thermostats and barostats, treatment of electrostatic interactions, and system size (1000 molecules). The densities and enthalpy of vaporization from an independent data set based on simulations using the CHARMM General Force Field (CGenFF) presented by Vanommeslaeghe et al. (J. Comput. Chem.2010, 31, 671) are included for comparison. We find that, overall, the OPLS/AA force field performs somewhat better than GAFF, but there are significant issues with reproduction of the surface tension and dielectric constants for both force fields.