DoGlycans-Tools for Preparing Carbohydrate Structures for Atomistic Simulations of Glycoproteins, Glycolipids, and Carbohydrate Polymers for GROMACS.
Ontology highlight
ABSTRACT:
SUBMITTER: Danne R
PROVIDER: S-EPMC5662928 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We describe a versatile method to enforce the rotation of subsets of atoms, e.g., a protein subunit, in molecular dynamics (MD) simulations. In particular, we introduce a "flexible axis" technique that allows realistic flexible adaptions of both the rotary subunit as well as the local rotation axis during the simulation. A variety of useful rotation potentials were implemented for the GROMACS 4.5 MD package. Application to the molecular motor F(1)-ATP synthase demonstrates the advantages of the flexible axis approach over the established fixed axis rotation technique.
| S-EPMC3091370 | biostudies-literature
Project description:In this Data in Brief article we provide a data package of GROMACS input files for atomistic molecular dynamics simulations of multicomponent, asymmetric lipid bilayers using the OPLS-AA force field. These data include 14 model bilayers composed of 8 different lipid molecules. The lipids present in these models are: cholesterol (CHOL), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylethanolamine (POPE), 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidyl-ethanolamine (SOPE), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (POPS), 1-stearoyl-2-oleoyl-sn-glycero-3-phosphatidylserine (SOPS), N-palmitoyl-D-erythro-sphingosyl-phosphatidylcholine (SM16), and N-lignoceroyl-D-erythro-sphingosyl-phosphatidylcholine (SM24). The bilayers׳ compositions are based on lipidomic studies of PC-3 prostate cancer cells and exosomes discussed in Llorente et al. (2013) [1], showing an increase in the section of long-tail lipid species (SOPS, SOPE, and SM24) in the exosomes. Former knowledge about lipid asymmetry in cell membranes was accounted for in the models, meaning that the model of the inner leaflet is composed of a mixture of PC, PS, PE, and cholesterol, while the extracellular leaflet is composed of SM, PC and cholesterol discussed in Van Meer et al. (2008) [2]. The provided data include lipids׳ topologies, equilibrated structures of asymmetric bilayers, all force field parameters, and input files with parameters describing simulation conditions (md.mdp). The data is associated with the research article "Interdigitation of Long-Chain Sphingomyelin Induces Coupling of Membrane Leaflets in a Cholesterol Dependent Manner" (Róg et al., 2016) [3].
| S-EPMC5063793 | biostudies-literature
Project description:Cholera toxin (CT) is the etiological agent of cholera. Here we report that multiple classes of fucosylated glycoconjugates function in CT binding and intoxication of intestinal epithelial cells. In Colo205 cells, knockout of B3GNT5, the enzyme required for synthesis of lacto- and neolacto-series glycosphingolipids (GSLs), reduces CT binding but sensitizes cells to intoxication. Overexpressing B3GNT5 to generate more fucosylated GSLs confers protection against intoxication, indicating that fucosylated GSLs act as decoy receptors for CT. Knockout (KO) of B3GALT5 causes increased production of fucosylated O-linked and N-linked glycoproteins, and leads to increased CT binding and intoxication. Knockout of B3GNT5 in B3GALT5 KO cells eliminates production of fucosylated GSLs but increases intoxication, identifying fucosylated glycoproteins as functional receptors for CT. These findings provide insight into molecular determinants regulating CT sensitivity of host cells.
2024-08-31 | GSE242156 | GEO
Project description:Sodium channels are chief proteins involved in electrical signaling in the nervous system, enabling critical functions like heartbeat and brain activity. New high-resolution X-ray structures for bacterial sodium channels have created an opportunity to see how these proteins operate at the molecular level. An important challenge to overcome is establishing relationships between the structures and functions of mammalian and bacterial channels. Bacterial sodium channels are known to exhibit the main structural features of their mammalian counterparts, as well as several key functional characteristics, including selective ion conduction, voltage-dependent gating, pore-based inactivation and modulation by local anesthetic, antiarrhythmic and antiepileptic drugs. Simulations have begun to shed light on each of these features in the past few years. Despite deviations in selectivity signatures for bacterial and mammalian channels, simulations have uncovered the nature of the multiion conduction mechanism associated with Na(+) binding to a high-field strength site established by charged glutamate side chains. Simulations demonstrated a surprising level of flexibility of the protein, showing that these side chains are active participants in the permeation process. They have also uncovered changes in protein structure, leading to asymmetrical collapses of the activation gate that have been proposed to correspond to inactivated structures. These observations offer the potential to examine the mechanisms of state-dependent drug activity, focusing on pore-blocking and pore-based slow inactivation in bacterial channels, without the complexities of inactivation on multiple timescales seen in eukaryotic channels. Simulations have provided molecular views of the interactions of drugs, consistent with sites predicted in mammalian channels, as well as a wealth of other sites as potential new drug targets. In this chapter, we survey the new insights into sodium channel function that have emerged from studies of simpler bacterial channels, which provide an excellent learning platform, and promising avenues for mechanistic discovery and pharmacological development.
| S-EPMC6591151 | biostudies-literature
Project description:Mixtures of long- and short-tail phosphatidylcholine lipids are known to self-assemble into a variety of aggregates combining flat bilayerlike and curved micellelike features, commonly called bicelles. Atomistic simulations of bilayer ribbons and perforated bilayers containing dimyristoylphosphatidylcholine (DMPC, di-C(14) tails) and dihexanoylphosphatidylcholine (DHPC, di-C(6) tails) have been carried out to investigate the partitioning of these components between flat and curved microenvironments and the stabilization of the bilayer edge by DHPC. To approach equilibrium partitioning of lipids on an achievable simulation timescale, configuration-bias Monte Carlo mutation moves were used to allow individual lipids to change tail length within a semigrand-canonical ensemble. Since acceptance probabilities for direct transitions between DMPC and DHPC were negligible, a third component with intermediate tail length (didecanoylphosphatidylcholine, di-C(10) tails) was included at a low concentration to serve as an intermediate for transitions between DMPC and DHPC. Strong enrichment of DHPC is seen at ribbon and pore edges, with an excess linear density of approximately 3 nm(-1). The simulation model yields estimates for the onset of edge stability with increasing bilayer DHPC content between 5% and 15% DHPC at 300 K and between 7% and 17% DHPC at 323 K, higher than experimental estimates. Local structure and composition at points of close contact between pores suggest a possible mechanism for effective attractions between pores, providing a rationalization for the tendency of bicelle mixtures to aggregate into perforated vesicles and perforated sheets.
| S-EPMC2884235 | biostudies-other
Project description:Cholera toxin (CT) is the etiological agent of cholera. Here we report that multiple classes of fucosylated glycoconjugates function in CT binding and intoxication of intestinal epithelial cells. In Colo205 cells, knockout of B3GNT5, the enzyme required for synthesis of lacto- and neolacto-series glycosphingolipids (GSLs), reduces CT binding but sensitizes cells to intoxication. Overexpressing B3GNT5 to generate more fucosylated GSLs confers protection against intoxication, indicating that fucosylated GSLs act as decoy receptors for CT. Knockout (KO) of B3GALT5 causes increased production of fucosylated O-linked and N-linked glycoproteins, and leads to increased CT binding and intoxication. Knockout of B3GNT5 in B3GALT5 KO cells eliminates production of fucosylated GSLs but increases intoxication, identifying fucosylated glycoproteins as functional receptors for CT. These findings provide insight into molecular determinants regulating CT sensitivity of host cells.
| S-EPMC10418270 | biostudies-literature
Project description:The molecular dynamics simulation package GROMACS runs efficiently on a wide variety of hardware from commodity workstations to high performance computing clusters. Hardware features are well-exploited with a combination of single instruction multiple data, multithreading, and message passing interface (MPI)-based single program multiple data/multiple program multiple data parallelism while graphics processing units (GPUs) can be used as accelerators to compute interactions off-loaded from the CPU. Here, we evaluate which hardware produces trajectories with GROMACS 4.6 or 5.0 in the most economical way. We have assembled and benchmarked compute nodes with various CPU/GPU combinations to identify optimal compositions in terms of raw trajectory production rate, performance-to-price ratio, energy efficiency, and several other criteria. Although hardware prices are naturally subject to trends and fluctuations, general tendencies are clearly visible. Adding any type of GPU significantly boosts a node's simulation performance. For inexpensive consumer-class GPUs this improvement equally reflects in the performance-to-price ratio. Although memory issues in consumer-class GPUs could pass unnoticed as these cards do not support error checking and correction memory, unreliable GPUs can be sorted out with memory checking tools. Apart from the obvious determinants for cost-efficiency like hardware expenses and raw performance, the energy consumption of a node is a major cost factor. Over the typical hardware lifetime until replacement of a few years, the costs for electrical power and cooling can become larger than the costs of the hardware itself. Taking that into account, nodes with a well-balanced ratio of CPU and consumer-class GPU resources produce the maximum amount of GROMACS trajectory over their lifetime.
| S-EPMC5042102 | biostudies-literature
Project description:In this data article we provide topologies and force field parameters files for molecular dynamics simulations of lipids in the OPLS-aa force field using the GROMACS package. This is the first systematic parameterization of lipid molecules in this force field. Topologies are provided for four phosphatidylcholines: saturated DPPC, mono-cis unsaturated POPC and DOPC, and mono-trans unsaturated PEPC. Parameterization of the phosphatidylcholines was achieved in two steps: first, we supplemented the OPLS force field parameters for DPPC with new parameters for torsion angles and van der Waals parameters for the carbon and hydrogen atoms in the acyl chains, as well as new partial atomic charges and parameters for torsion angles in the phosphatidylcholine and glycerol moieties [1]. Next, we derived parameters for the cis and trans double bonds and the neighboring them single bonds [2]. Additionally, we provide GROMACS input files with parameters describing simulation conditions (md.mdp), which are strongly recommended to be used with these lipids models. The data are associated with the research article "Cis and trans unsaturated phosphatidylcholine bilayers: a molecular dynamics simulation study" [2] and provided as supporting materials.
| S-EPMC4602361 | biostudies-literature
Project description:The celebrated Miller experiments reported on the spontaneous formation of amino acids from a mixture of simple molecules reacting under an electric discharge, giving birth to the research field of prebiotic chemistry. However, the chemical reactions involved in those experiments have never been studied at the atomic level. Here we report on, to our knowledge, the first ab initio computer simulations of Miller-like experiments in the condensed phase. Our study, based on the recent method of treatment of aqueous systems under electric fields and on metadynamics analysis of chemical reactions, shows that glycine spontaneously forms from mixtures of simple molecules once an electric field is switched on and identifies formic acid and formamide as key intermediate products of the early steps of the Miller reactions, and the crucible of formation of complex biological molecules.
| S-EPMC4183268 | biostudies-literature
Project description:The accuracy of atomistic biomolecular modeling and simulation studies depend on the accuracy of the input structures. Preparing these structures for an atomistic modeling task, such as molecular dynamics (MD) simulation, can involve the use of a variety of different tools for: correcting errors, adding missing atoms, filling valences with hydrogens, predicting pK values for titratable amino acids, assigning predefined partial charges and radii to all atoms, and generating force field parameter/topology files for MD. Identifying, installing and effectively using the appropriate tools for each of these tasks can be difficult for novice and time-consuming for experienced users. H++ (http://biophysics.cs.vt.edu/) is a free open-source web server that automates the above key steps in the preparation of biomolecular structures for molecular modeling and simulations. H++ also performs extensive error and consistency checking, providing error/warning messages together with the suggested corrections. In addition to numerous minor improvements, the latest version of H++ includes several new capabilities and options: fix erroneous (flipped) side chain conformations for HIS, GLN and ASN, include a ligand in the input structure, process nucleic acid structures and generate a solvent box with specified number of common ions for explicit solvent MD.
| S-EPMC3394296 | biostudies-literature