Exploring the Molecular Dynamics of a Lipid-A Vesicle at the Atom Level: Morphology and Permeation Mechanism.

Lipid-A was previously shown to spontaneously aggregate into a vesicle via the hybrid particle field approach. We assess the validity of the proposed vesiculation mechanism by simulating the resulting lipid-A vesicle at the atom level. The spatial confinement imposed by the vesicle geometry on the conformation and packing of lipid-A induces significant heterogeneity of physical properties in the inner and outer leaflets. It also induces tighter molecular packing and lower acyl chain order compared to the lamellar arrangement. Around 5% of water molecules passively permeates the vesicle membrane inward and outward. The permeation is facilitated by interactions with water molecules that are transported across the membrane by a network of electrostatic interactions with the hydrogen bond donors/acceptors in the N-acetylglucosamine ring and upper region of the acyl chains of lipid-A. The permeation process takes place at low rates but still at higher frequencies than observed for the lamellar arrangement of lipid-A. These findings not only substantiate the proposed lipid-A vesiculation mechanism but also reveal the complex structural dynamics of an important nonlamellar arrangement of lipid-A.

Surface Assessment via Grid Evaluation (SuAVE) for Every Surface Curvature and Cavity Shape.

The surface assessment via grid evaluation (SuAVE) software was developed to account for the effect of curvature in the calculations of structural properties of chemical interfaces regardless of the chemical composition, asymmetry, and level of atom coarseness. It employs differential geometry techniques, enabling the representation of chemical surfaces as fully differentiable. In this article, we present novel developments of SuAVE to treat closed surfaces and complex cavity shapes. These developments expand the repertoire of curvature-dependent analyses already available in the previous version of SuAVE (e.g., area per lipid, density profiles, membrane thickness, deuterium-order parameters, volume per lipid, and surface curvature angle) to include new functionalities applicable to soft matter (e.g., sphericity, average radius, principal moment of inertia, and roundness) and crystalline porous materials (e.g., pore diameter, internal void volume, total area, and the total void volume of the unit cell structure). SuAVE can accurately handle chemical systems with high and low atom density as demonstrated for two distinct chemical systems: the lipid A vesicle and a set of selected metal-organic frameworks. The SuAVE software v2.0 is fully parallel and benefits from a compiler that supports OpenMP. SuAVE is freely available from https://github.com/SuAVE-Software/source and https://www.biomatsite.net/.

The tug of war between Al3+ and Na+ for order-disorder transitions in lipid-A membranes.

Cations play a critical role in the stability and morphology of lipid-A aggregates by neutralizing, hydrating and cross-linking these glycolipid molecules. Monophosphorylated lipid-A is the major immunostimulatory principle in commercially available adjuvants containing Al3+ such as adjuvant system 04 (AS04). The antagonist/agonist immunomodulatory properties of lipid-A are associated with chemical variations (e.g. the number of acyl chains and phosphate groups) and their aggregate arrangements (e.g. lamellar, nonlamellar or mixed). Therefore, the identification of the active form of lipid-A can provide valuable guidance in the development of vaccine adjuvants capable of boosting the immune system with decreased reactogenicity. Although the effect of mono and divalent cations on the structural polymorphism and endotoxicity of LPS has been previously investigated, much less is known about the effect of trivalent cations. We have investigated the effect of NaCl and AlCl3 salt solutions on the structural dynamics and stability of mono and diphosphorylated lipid-A membranes via atomistic MD simulations. The Al3+ ion exerts two major effects on the structural dynamics of lipid-A membranes. It acts as an efficient cross-linker of mono or diphosphorylated lipid-A molecules, thus stabilizing the lamellar arrangement of these glycolipids. It also alters the lipid-A packing and membrane fluidity, inducing disorder → order structural transitions of the membrane. This effect is promptly reversed upon the addition of NaCl solution, which promotes a nearly threefold increase in the amount of water in the carbohydrate moiety of the Al3+-containing lipid-A membranes. The exchange dynamics and residence times of cation-coordinated water molecules in these membranes provide insights into the molecular mechanism for the Na+-induced transition from a densely packed ordered phase to a disordered one. Al3+ counter-ions favor ordered lamellar aggregates, which has been previously associated with the lack of endotoxic activity and cytokine-inducing action. The resulting microscopic understanding of the structure and dynamics of lipid-A aggregates in the presence of Al3+ and Na+ salts can provide valuable guidance in the development of vaccine adjuvants capable of boosting the immune system with decreased reactogenicity.

Aggregation of Lipid A Variants: A Hybrid Particle-Field Model.

Lipid A is one of the three components of bacterial lipopolysaccharides constituting the outer membrane of Gram-negative bacteria, and is recognized to have an important biological role in the inflammatory response of mammalians. Its biological activity is modulated by the number of acyl-chains that are present in the lipid and by the dielectric medium, i.e., the type of counter-ions, through electrostatic interactions. In this paper, we report on a coarse-grained model of chemical variants of Lipid A based on the hybrid particle-field/molecular dynamics approach (hPF-MD). In particular, we investigate the stability of Lipid A bilayers for two different hexa- and tetra-acylated structures. Comparing particle density profiles along bilayer cross-sections, we find good agreement between the hPF-MD model and reference all-atom simulation for both chemical variants of Lipid A. hPF-MD models of constituted bilayers composed by hexa-acylated Lipid A in water are stable within the simulation time. We further validate our model by verifying that the phase behavior of Lipid A/counterion/water mixtures is correctly reproduced. In particular, hPF-MD simulations predict the correct self-assembly of different lamellar and micellar phases from an initially random distribution of Lipid A molecules with counterions in water. Finally, it is possible to observe the spontaneous formation and stability of Lipid A vesicles by fusion of micellar aggregates.

Out of Sight, Out of Mind: The Effect of the Equilibration Protocol on the Structural Ensembles of Charged Glycolipid Bilayers.

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.

Rotational Profiler: A Fast, Automated, and Interactive Server to Derive Torsional Dihedral Potentials for Classical Molecular Simulations.

Rotational Profiler provides an analytical algorithm to compute sets of classical torsional dihedral parameters by fitting an empirical energy profile to a reference one that can be obtained experimentally or by quantum-mechanical methods. The resulting profiles are compatible with the functional forms in the most widely used biomolecular force fields (e.g., GROMOS, AMBER, OPLS, and CHARMM). The linear least-squares regression method is used to generate sets of parameters that best satisfy the fitting. Rotational Profiler is free to use, analytical, and force field/package independent. The formalism is herein described, and its usage, in an interactive and automated manner, is made available as a Web server at http://rotprof.lncc.br.

SuAVE: A Tool for Analyzing Curvature-Dependent Properties in Chemical Interfaces.

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 www.biomatsite.net/suave-software .

Conformational Dynamics and Responsiveness of Weak and Strong Polyelectrolyte Brushes: Atomistic Simulations of Poly(dimethyl aminoethyl methacrylate) and Poly(2-(methacryloyloxy)ethyl trimethylammonium chloride).

The complex solution behavior of polymer brushes is key to control their properties, including for biomedical applications and catalysis. The swelling behavior of poly(dimethyl aminoethyl methacrylate) (PDMAEMA) and poly(2-(methacryloyloxy)ethyl trimethylammonium chloride) (PMETAC) in response to changes in pH, solvent, and salt types has been investigated using atomistic molecular dynamics simulations. PDMAEMA and PMETAC have been selected as canonical models for weak and strong polyelectrolytes whose complex conformational behavior is particularly challenging for the development and validation of atomistic models. The GROMOS-derived atomic parameters reproduce the experimental swelling coefficients obtained from ellipsometry measurements for brushes of 5-15 nm thickness. The present atomistic models capture the protonated morphology of PDMAEMA, the swollen and collapsed conformations of PDMAEMA and PMETAC in good and bad solvents, and the salt-selective response of PMETAC. The modular nature of the molecular models allows for the simple extension of atomic parameters to a variety of polymers or copolymers.

Polymyxin Binding to the Bacterial Outer Membrane Reveals Cation Displacement and Increasing Membrane Curvature in Susceptible but Not in Resistant Lipopolysaccharide Chemotypes.

Lipid-A is the causative agent of Gram-negative sepsis and is responsible for an increasingly high mortality rate among hospitalized patients. Compounds that bind Lipid-A can limit this inflammatory process. The cationic antimicrobial peptide polymyxin B (Pmx-B) is one of the simplest molecules capable of selectively binding to Lipid-A and may serve as a model for further development of Lipid-A binding agents. Gram-negative bacteria resistance to Pmx-B relies on the upregulation of a number of regulatory systems, which promote chemical modifications of the lipopolysaccharide (LPS) structure and leads to major changes in the physical-chemical properties of the outer membrane. A detailed understanding of how the chemical structure of the LPS modulates macroscopic properties of the outer membrane is paramount for the design and optimization of novel drugs targeting clinically relevant strains. We have performed a systematic investigation of Pmx-B binding to outer membrane models composed of distinct LPS chemotypes experimentally shown to be either resistant or susceptible to the peptide. Molecular dynamics simulations were carried out for Pmx-B bound to the penta- and hexa-acylated forms of Lipid-A (more susceptible) and Lipid-A modified with 4-amino-4-deoxy-l-arabinose (resistant) as well as the penta-acylated form of LPS Re (less susceptible). The present simulations show that upon binding to the bacterial outer membrane surface, Pmx-B promotes cation displacement and structural changes in membrane curvature and integrity as a function of the LPS chemotype susceptibility or resistance to the antimicrobial peptide.