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

Exploring How System Dimensions and Periodic Boundary Conditions Influence the Molecular Dynamics Simulation of A6H Peptide Self-Assembly Nanostructures

This work presents a study on the effects of periodic boundary conditions (PBC) on the energetic/structural properties and hydrogen bond dynamics (HB) using molecular dynamics (MD) simulations of peptide membranes composed of alanine and histidine. Our results highlight that simulations using small surface areas for the peptide membrane may result in nonconvergent values for membrane properties, which are only observed in regions simulated at a certain distance from the PBCs. Specifically, regarding hydrogen bonds, a property pervasive in peptide membranes, our findings indicate a significant increase in the lifetime of these interactions, reaching values ∼19% higher when observed in structures free from PBCs. For peptide mobility in these nanomembranes, our results compare regions simulated directly under the influence of PBCs with regions free from these conditions, emphasizing greater mobility of amino acid psi/phi angles in the latter model.

Interfacial Glucose to Regulate Hydrated Lipid Bilayer Properties: Influence of Concentrations

A series of atomistic molecular dynamics (MD) simulations were carried out with a hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer with the variation of glucose concentrations from 0 to 30 wt % in the presence of 0.3 M NaCl. The study suggested that although the thickness of the lipid bilayer dropped significantly with the increase in glucose concentration, it expanded laterally at high glucose levels due to the intercalation of glucose between the headgroups of adjacent lipids. We adopted the surface assessment via the grid evaluation method to compute the deviation of the bilayer's key structural features for the different amounts of glucose present. This suggested that the accumulation of glucose molecules near the headgroups influences the local lipid bilayer undulation and crimping of the lipid tails. We find that the area compressibility modulus increases with the glucose level, causing enhanced bilayer rigidity arising from the slow lateral diffusion of lipids. The restricted lipid motion at high glucose concentrations controls the sustainability of the curved bilayer surface. Calculations revealed that certain orientations of CO → of interfacial glucose with the PN → of lipid headgroups are preferred, which helps the glucose to form direct hydrogen bonds (HBs) with the lipid headgroups. Such lipid-glucose (LG) HBs relax slowly at low glucose concentrations and exhibit a higher lifetime, whereas fast structural relaxation of LG HBs with a shorter lifetime was noticed at a higher glucose level. In contrast, lipid-water (LW) HBs exhibited a higher lifetime at a higher glucose level, which gradually decreased with the glucose level lowering. The study interprets that the glucose concentration-driven LW and LG interactions are mutually inclusive. Our detailed analysis will exemplify small saccharide concentration-driven membrane stabilizing efficiency, which is, in general, helpful for drug delivery study.

Vesicle protrusion induced by antimicrobial peptides suggests common carpet mechanism for short antimicrobial peptides

Short-cationic alpha-helical antimicrobial peptides (SCHAMPs) are promising candidates to combat the growing global threat of antimicrobial resistance. They are short-sequenced, selective against bacteria, and have rapid action by destroying membranes. A full understanding of their mechanism of action will provide key information to design more potent and selective SCHAMPs. Molecular Dynamics (MD) simulations are invaluable tools that provide detailed insights into the peptide-membrane interaction at the atomic- and meso-scale level. We use atomistic and coarse-grained MD to look into the exact steps that four promising SCHAMPs-BP100, Decoralin, Neurokinin-1, and Temporin L-take when they interact with membranes. Following experimental set-ups, we explored the effects of SCHAMPs on anionic membranes and vesicles at multiple peptide concentrations. Our results showed all four peptides shared similar binding steps, initially binding to the membrane through electrostatic interactions and then flipping on their axes, dehydrating, and inserting their hydrophobic moieties into the membrane core. At higher concentrations, fully alpha-helical peptides induced membrane budding and protrusions. Our results suggest the carpet mode of action is fit for the description of SCHAMPs lysis activity and discuss the importance of large hydrophobic residues in SCHAMPs design and activity.

Analyzing Lipid Membrane Defects via a Coarse-Grained to Triangulated Surface Map: The Role of Lipid Order and Local Curvature in Molecular Binding

Lipid packing defects are known to serve as quantitative indicators for protein binding to lipid membranes. This paper presents a protocol for mapping molecular lipid detail onto a triangulated continuum leaflet representation. Besides establishing the desired forward counterpart to the existing inverse TS2CG map, this coarse-grained to triangulated surface (CG2TS) map enables straightforward extraction of the defect characteristics for any membrane geometry found in nature. We have applied our protocol to investigate the role of local curvature and varying lipid packing on the defect constant π. We find that the defect size is greatly influenced by both factors, arguing strongly against the usual assignment of a single defect constant in the case of more realistic membrane conditions. An important discovery is that lipids in the gel phase produce larger defects, or a higher π, in domains of high (local) curvature than the same lipid in a liquid phase of any curvature. This finding suggests that membranes featuring very ordered lipid packing can bind proteins via large defects in curved regions. Finally, we propose a route for estimating defect constants directly from the standard membrane properties. Identifying the precise role of composition, lipid (tail) order, and (local) curvature in defects for the irregular lipid structures that are (temporally) present in many biological processes is instrumental for obtaining fundamental insight as well as for a rational design of membrane binding targets.

The Mechanism of Antimicrobial Small-Cationic Peptides from Coarse-Grained Simulations

Antimicrobial cationic peptides (AMPs) are excellent candidates for use as therapeutic antimicrobial agents. Among them, short peptides possessing sequences of 9-11 amino acids have some advantages over long-sequence peptides. However, one of the main limitations of short peptides is that their mechanism of action at the molecular level is not well-known. In this article, we report a model based on multiscale molecular dynamics simulations of short peptides interacting with vesicles containing palmitoyl-oleoyl-phosphatidylglycerol (POPG)/palmitoyl-oleoyl-phosphatidylethanolamine (POPE). Simulations using this approach have allowed us to understand the different behaviors of peptides with antimicrobial activity with respect to those that do not produce this effect. We found remarkable agreement with a series of experimental results directly supporting our model. Moreover, these results allow us to understand the mechanism of action at the molecular level of these short peptides. Our simulations suggest that mechanical inhomogeneities appear in the membrane, promoting membrane rupture when a threshold concentration of peptides adsorbed on the membrane is achieved. These results explain the high structural demand for these peptides to maintain a delicate balance between the affinity for the bilayer surface, a low peptide-peptide repulsion (in order to reach the threshold concentration), and an acceptable tendency to penetrate into the bilayer. This mechanism is different from those proposed for peptides with long amino acid sequences. Such information is very useful from the medicinal chemistry point of view for the design of new small antimicrobial peptides.

Accelerating Lipid Flip-Flop at Low Concentrations: A General Mechanism for Membrane Binding Peptides

We report a physicochemical investigation of the lipid transport properties of model lipid membranes in the presence of the antimicrobial peptide indolicidin through comparisons of experimental SANS/SAXS scattering techniques to fully atomistic molecular dynamics simulations. In agreement with the experiment, we show that upon peripheral binding of the peptides, even at low concentrations, lipid flip-flop dynamics is greatly accelerated. Computer modeling elucidates the interplay between structural changes and lipid dynamics induced by peptides and proposes a mechanism for the mode of action of antimicrobial peptides, assessing the major role of entropy for the catalysis of the flipping events. The mechanism introduced here is universal for all peptides with preferential peripheral binding to the membrane as it does not depend on the specific amino acid sequence.

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.

A brief history of visualizing membrane systems in molecular dynamics simulations

Understanding lipid dynamics and function, from the level of single, isolated molecules to large assemblies, is more than ever an intensive area of research. The interactions of lipids with other molecules, particularly membrane proteins, are now extensively studied. With advances in the development of force fields for molecular dynamics simulations (MD) and increases in computational resources, the creation of realistic and complex membrane systems is now common. In this perspective, we will review four decades of the history of molecular dynamics simulations applied to membranes and lipids through the prism of molecular graphics.

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/.

Self-Assembly of Lipid Mixtures in Solutions: Structures, Dynamics Processes and Mechanical Properties

The self-assembly of lipid mixtures in aqueous solution was investigated by dissipative particle dynamics simulation. Two types of lipid molecules were modelled, where three mixed structures, i.e., the membrane, perforated membrane and vesicle, were determined in the self-assembly processes. Phase behaviour was investigated by using the phase diagrams based on the tail chain lengths for the two types of lipids. Several parameters, such as chain number and average radius of gyration, were employed to explore the structural formations of the membrane and perforated membrane in the dynamic processes. Interface tension was used to demonstrate the mechanical properties of the membrane and perforated membrane in the equilibrium state and dynamics processes. Results help us to understand the self-assembly mechanism of the biomolecule mixtures, which has a potential application for designing the lipid molecule-based bio-membranes in solutions.

Self-assembly of an in silico designed dipeptide derivative to obtain photo-responsive vesicles

Photo-responsive vesicles self-assembled from in silico designed peptide derivatives were investigated using coarse-grained molecular dynamics simulations.

Simulations reveal that antimicrobial BP100 induces local membrane thinning, slows lipid dynamics and favors water penetration

MD simulations reveal that BP100 peptide induces local membrane thinning and negative curvature, slows lipid dynamics and increases the water life time in the lipid hydrophobic core and transmembrane water transport in the direction of the peptide.

The Role of Key Amino Acids in the Antimicrobial Mechanism of a Bacteriocin Model Revealed by Molecular Simulations

The AS-48 bacteriocin is a potent antimicrobial polypeptide with enhanced stability due to its circular sequence of peptidic bonds. The mechanism of biological action is still not well understood in spite of both the elucidation of the molecular structure some years ago and several experiments performed that yielded valuable information about the AS-48 bacterial membrane poration activity. In this work, we present a computational study at an atomistic scale to analyze the membrane disruption mechanism. The process is based on the two-stage model: (1) peptide binding to the bilayer surface and (2) membrane poration due to the surface tension exerted by the peptide. Indeed, the induced membrane tension mechanism is able to explain stable formation of pores leading to membrane disruption. The atomistic detail obtained from the simulations allows one to envisage the contribution of the different amino acids during the poration process. Clustering of cationic residues and hydrophobic interactions between peptide and lipids seem to be essential ingredients in the process. GLU amino acids have shown to enhance the membrane disrupting ability of the bacteriocin. TRP24–TRP24 interactions make also an important contribution in the initial stages of the poration mechanism. The detailed atomistic information obtained from the simulations can serve to better understand bacteriocin structural characteristics to design more potent antimicrobial therapies.

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.

Impact of Increased Membrane Realism on Conformational Sampling of Proteins

The realism and accuracy of lipid bilayer simulations through molecular dynamics (MD) are heavily dependent on the lipid composition. While the field is pushing toward implementing more heterogeneous and realistic membrane compositions, a lack of high-resolution lipidomic data prevents some membrane protein systems from being modeled with the highest level of realism. Given the additional diversity of real-world cellular membranes and protein-lipid interactions, it is still not fully understood how altering membrane complexity affects modeled membrane protein functions or if it matters over long-timescale simulations. This is especially true for organisms whose membrane environments have little to no computational study, such as the plant plasma membrane. Tackling these issues in tandem, a generalized, realistic, and asymmetric plant plasma membrane with more than 10 different lipid species is constructed herein. Classical MD simulations of pure membrane constructs were performed to evaluate how altering the compositional complexity of the membrane impacted the plant membrane properties. The apo form of a plant sugar transporter, OsSWEET2b, was inserted into membrane models where lipid diversity was calculated in either a size-dependent or size-independent manner. An adaptive sampling simulation regime validated by Markov-state models was performed to capture the gating dynamics of OsSWEET2b in each of these membrane constructs. In comparison to previous OsSWEET2b simulations performed in a pure POPC bilayer, we confirm that simulations performed within a native-like membrane composition alter the stabilization of apo OsSWEET2b conformational states by ∼1 kcal/mol. The free-energy barriers of intermediate conformational states decrease when realistic membrane complexity is simplified, albeit roughly within sampling error, suggesting that protein-specific responses to membranes differ due to altered packing caused by compositional fluctuations. This work serves as a case study where a more realistic bilayer composition makes unbiased conformational sampling easier to achieve than with simplified bilayers.

Lipid Oxidation: Role of Membrane Phase-Separated Domains

Lipid oxidation is associated with several inflammatory and neurodegenerative diseases, but many questions to unravel its effects on biomembranes are still open due to the complexity of the topic. For instance, recent studies indicated that phase-separated domains can have a significant effect on membrane function. It is reported that domain interfaces are "hot spots" for pore formation, but the underlying mechanisms and the effect of oxidation-induced phase separation on membranes remain elusive. Thus, to evaluate the permeability of the membrane coexisting of liquid-ordered (Lo) and liquid-disordered (Ld) domains, we performed atomistic molecular dynamics simulations. Specifically, we studied the membrane permeability of nonoxidized or oxidized homogeneous membranes (single-phase) and at the Lo/Ld domain interfaces of heterogeneous membranes, where the Ld domain is composed of either oxidized or nonoxidized lipids. Our simulation results reveal that the addition of only 1.5% of lipid aldehyde molecules at the Lo/Ld domain interfaces of heterogeneous membranes increases the membrane permeability, whereas their addition at homogeneous membranes does not have any effect. This study is of interest for a better understanding of cancer treatment methods based on oxidative stress (causing among others lipid oxidation), such as plasma medicine and photodynamic therapy.

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.