Lipopolysaccharides (LPS) are the main constituent of the outer bacterial membrane of Gram-negative bacteria. Lipid-A is the structural region of LPS that interacts with the innate immune system and induces inflammatory responses. It is formed by a phosphorylated β-d-glucosaminyl-(1→6)-α-N-glucosamine disaccharide backbone containing ester-linked and amide-linked long-chain fatty acids, which may vary in length and number depending on the bacterial strains and the environment. Phenotypical variation (i.e., number of acyl chains), cation type, and temperature influence the phase transition, aggregate structure, and endotoxic activity of Lipid-A. We have applied an extension of the GROMOS force field 45a4 carbohydrate parameter set to investigate the behavior of hexa- and pentaacylated Lipid-A of Pseudomonas aeruginosa at two temperatures (300 and 328 K) and in the presence of mono- and divalent cations (represented by Ca(2+) and Na(+), respectively) through molecular dynamics simulations. The distinct phase of Lipid-A aggregates was characterized by structural properties, deuterium order parameters, the molecular shape of the lipid units (conical versus cylindrical), and molecular packing. Our results show that Na(+) ions induce a transition from the lamellar to nonlamellar phase. In contrast, the bilayer integrity is maintained in the presence of Ca(2+) ions. Through these findings, we present microscopic insights on the influence of different cations on the molecular behavior of Lipid-A associated with the lamellar to nonlamellar transition.
The supra-molecular structure of LPS aggregates governs outer membrane permeability and activation of the host immune response during Gram-negative bacterial infections. Molecular dynamics simulations unveil at atomic resolution the subtle balance between cation hydration and cross-linking ability in modulating phase transitions of LPS membranes.
Curvature is an intrinsic feature of biological membranes underlying vital cellular processes such as endocytosis, membrane fusion–fission, trafficking, and remodeling. The continuous expansion of the spatiotemporal scales accessible to computational simulations nowadays makes possible quasi-atomistic molecular dynamics simulations of these processes. In despite of that, computation of the shapes and curvatures associated with the dynamics of biological membranes remains challenging. For this reason, the effect of curvature is often neglected in the analysis of quantities essential for the accurate description of membrane properties (e.g., area and volume per lipid, density profiles, membrane thickness). We propose an algorithm for surface assessment via grid evaluation (SuAVE) that relies on the application of a radial base function to interpolate points scattered across an interface of any shape. This enables the representation of the chemical interface as fully differentiable so that related geometrical properties can be calculated through the straightforward employment of well-established differential geometry techniques. Hence, the effect of different types or degrees of curvature can be accurately taken into account in the calculations of structural properties of any interfaces regardless of chemical composition, asymmetry, and level of atom coarseness. The main functionalities implemented in SuAVE are featured for a number of tetraacylated and hexaacylated Lipid-A membranes of distinct curvatures and a surfactant micelle. We show that the properties calculated for moderately to highly curved membranes differ significantly between curvature-dependent and -independent algorithms. The SuAVE software is freely available from .
Glutathione S-transferases (GSTs) are enzymes capable of metabolizing cytotoxic compounds. The enzyme AgGSTE2, member of epsilon class GSTs (GSTE), is the most important GST conferring resistance to dichloro-diphenyl-trichloroethane (DDT) in Anopheles gambiae. We have investigated the conformational dynamics of three GSTE variants (GSTE2, GSTE2-I114T/F120L, GSTE5) from A. gambiae. Large-scale motions of helices H2 and H4 and conformational transition of the C-terminal governs the opening of the G-site and is expected to affect substrate binding and product release. This structural rearrangement places Glu116 (Glu120 in GSTE5) close of the thiol group of the tripeptide glutathione (GSH) cofactor, making this residue a candidate to act as a base in the activation of DDT. The structural rearrangement is noticeable for AgGSTE2-F120L, which has been shown to confer increased DDT-resistance. The other variants exhibit a more subtle rearrangement. These findings corroborate the hypothesis that the increase of the conformational dynamics of GST Epsilon class isoforms from A. gambiae promotes higher DDTase activity.Keywords: molecular dynamics simulation, evolutional constraint, positive and negative selection, metabolic resistance, malaria vector IntroductionGlutathione transferases (EC 2.5.1.18) are highly promiscuous proteins where broad functional promiscuity co-exists with highly conserved structural fold. Glutathione S-transferases (GSTs) constitute a large family of cytosolic and membrane-bound proteins that catalyze the nucleophilic addition of the tripeptide glutathione (GSH) to a variety of electrophilic toxins and drugs (xenobiotics), leading to their excretion. 1,2 Lately, several other activities have also been associated with GSTs, including steroid and leukotriene biosynthesis, peroxide degradation, doublebond cis-trans isomerization, dehydroascorbate reduction, Michael addition, and noncatalytic ligand binding and transport activity. 3 These polymorphic enzymes belong to supergene families whose members are generated by punctual mutations, gene duplication and alternative splicing. 4 The GST family is subdivided into several classes accordingly to its occurrence in different taxa.The number of isoforms per class varies widely, ranging from one to forty. 4 A single GST isoform is capable of conjugating glutathione to several hydrophobic substrates. Such promiscuity when coupled with the large number of isozymes generates a large range of potential substrates.Insects exhibit at least six classes of GSTs (Sigma, Omega, Theta, Zeta, Delta and Epsilon) with the first four classes present in almost all living organisms. 4,5 The Epsilon and Delta classes are arthropod specific and some members of these classes have been associated to insecticide resistance in culicid vectors. 5,6 The metabolism of toxic compounds in insects is conducted by a series of enzymes from different phases and GSTs act in phase II. GSTs metabolize these compounds in two ways: one has been cited above (through GSH conjugation) and the oth...
Molecular dynamics (MD) simulations represent an essential tool in the toolbox of modern chemistry, enabling the prediction of experimental observables for a variety of chemical systems and processes and majorly impacting the study of biological membranes. However, the chemical diversity of complex lipids beyond phospholipids brings new challenges to well-established protocols used in MD simulations of soft matter and requires continuous assessment to ensure simulation reproducibility and minimize unphysical behavior. Lipopolysaccharides (LPS) are highly charged glycolipids whose aggregation in a lamellar arrangement requires the binding of numerous cations to oppositely charged groups deep inside the membrane. The delicate balance between the fully hydrated carbohydrate region and the smaller hydrophobic core makes LPS membranes very sensitive to the choice of equilibration protocol. In this work, we show that the protocol successfully used to equilibrate phospholipid bilayers when applied to complex lipopolysaccharide membranes occasionally leads to a small expansion of the simulation box very early in the equilibration phase. Although the use of a barostat algorithm controls the system dimension and particle distances according to the target pressure, fluctuation in the fleeting pressure occasionally enables a few water molecules to trickle into the hydrophobic region of the membrane, with spurious solvent buildup. We show that this effect stems from the initial steps of NPT equilibration, where initial pressure can be fairly high. This can be solved with the use of a stepwise-thermalization NVT/NPT protocol, as demonstrated for atomistic MD simulations of LPS/DPPE and lipid-A membranes in the presence of different salts using an extension of the GROMOS forcefield within the GROMACS software. This equilibration protocol should be standard procedure for the generation of consistent structural ensembles of charged glycolipids starting from atomic coordinates not previously pre-equilibrated. Although different ways to deal with this issue can be envisioned, we investigated one alternative that could be readily available in major MD engines with general users in mind.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.