Microsecond molecular dynamics simulations of lipid mixing

Nanoscale lipid domains determine the dynamic molecular portraits of mixed DOPC/DOPS bilayers in a fluid phase: A computational insight

Lateral heterogeneity, or mosaicity, is a fundamental property inherent to cell membranes that is crucial for their functioning. While microscopic inhomogeneities (e.g. rafts) are easily detected experimentally, lipid domains with nanoscale dimensions (nanoclusters of nanodomains, NDs) resist reliable characterization by instrumental methods. In such a case, important insight can be gained via computer modeling. Here, NDs composed of lipid's head groups in the mixed zwitterionic dioleoylphosphatidylcholine (DOPC) and negatively charged dioleoylphosphatidylserine (DOPS) bilayers were studied by molecular dynamics. A new algorithm has been developed to identify NDs. Unlike most similar methods, it implicitly considers the heterogeneous distribution of lipid head atomic density and does not require subjectively chosen parameters. In DOPS-rich membranes, lipids form more compact and stable NDs due to strong interlipid interactions. In DOPC-rich systems, NDs arise due to the "packing" effect of weakly bound lipid heads. The clustering picture is related to the physical properties of the bilayer surface: DOPS-rich systems show more pronounced surface heterogeneity of hydrophilic/hydrophobic regions compared to DOPC-rich ones. The results obtained are important for the effective quantitative characterization of the "dynamic molecular portrait" of a membrane surface - its "fingerprint" characterizing dynamical distribution of its physicochemical properties.

Insertion and Anchoring of HIV-1 Fusion Peptide into Complex Membrane Mimicking Human T-cell

A fundamental understanding of how HIV-1 envelope (Env) protein facilitates fusion is still lacking. The HIV-1 fusion peptide, consisting of 15 to 22 residues, is the N-terminus of the gp41 subunit of the Env protein. Further, this peptide, a promising vaccine candidate, initiates viral entry into target cells by inserting and anchoring into human immune cells. The influence of membrane lipid reorganization and the conformational changes of the fusion peptide during the membrane insertion and anchoring processes, which can significantly affect HIV-1 cell entry, remains largely unexplored due to the limitations of experimental measurements. In this work, we investigate the insertion of the fusion peptide into an immune cell membrane mimic through multiscale molecular dynamics simulations. We mimic the native T-cell by constructing a 9-lipid asymmetric membrane, along with geometrical restraints accounting for insertion in the context of gp41. To account for the slow timescale of lipid mixing while enabling conformational changes, we implement a protocol to go back and forth between atomistic and coarse-grained simulations. Our study provides a molecular understanding of the interactions between the HIV-1 fusion peptide and the T-cell membrane, highlighting the importance of conformational flexibility of fusion peptides and local lipid reorganization in stabilizing the anchoring of gp41 into the targeted host membrane during the early events of HIV-1 cell entry. Importantly, we identify a motif within the fusion peptide critical for fusion that can be further manipulated in future immunological studies.

Improved Highly Mobile Membrane Mimetic Model for Investigating Protein-Cholesterol Interactions

Cholesterol (CHL) plays an integral role in modulating the function and activity of various mammalian membrane proteins. Due to the slow dynamics of lipids, conventional computational studies of protein-CHL interactions rely on either long-time scale atomistic simulations or coarse-grained approximations to sample the process. A highly mobile membrane mimetic (HMMM) has been developed to enhance lipid diffusion and thus used to facilitate the investigation of lipid interactions with peripheral membrane proteins and, with customized in silico solvents to replace phospholipid tails, with integral membrane proteins. Here, we report an updated HMMM model that is able to include CHL, a nonphospholipid component of the membrane, henceforth called HMMM-CHL. To this end, we had to optimize the effect of the customized solvents on CHL behavior in the membrane. Furthermore, the new solvent is compatible with simulations using force-based switching protocols. In the HMMM-CHL, both improved CHL dynamics and accelerated lipid diffusion are integrated. To test the updated model, we have applied it to the characterization of protein-CHL interactions in two membrane protein systems, the human β2-adrenergic receptor (β2AR) and the mitochondrial voltage-dependent anion channel 1 (VDAC-1). Our HMMM-CHL simulations successfully identified CHL binding sites and captured detailed CHL interactions in excellent consistency with experimental data as well as other simulation results, indicating the utility of the improved model in applications where an enhanced sampling of protein-CHL interactions is desired.

Widom Line in Supercritical Water in Terms of Changes in Local Structure: Theoretical Perspective

Performing molecular dynamics simulations with the TIP4P/2005 water model along 9 isobars (from 175 to 375 bar) in the temperature range between 300 and 1100 K, we have found that the loci of the extrema in the rate of change of specific structural properties can be used to define purely structure-based Widom lines. We have examined several parameters that describe the local structure of water, such as the tetrahedral arrangement, nearest neighbor distance, local density around the molecules, and the size of the largest dense domain. The last two parameters were determined using the Voronoi polyhedral and density-based spatial clustering of applications with noise methods, respectively. By analyzing the moments of the associated distributions, we show that along a given isobar, the temperature at which we observe a maximum in the fluctuation, the rate of change of the average values, or in the skewness values unambiguously determines the Widom line that is in agreement with the experimentally detected, thermodynamic response function-based ones.

Structure of the Bacterial Cell Envelope and Interactions with Antimicrobials: Insights from Molecular Dynamics Simulations

Bacteria have evolved over 3 billion years, shaping our intrinsic and symbiotic coexistence with these single-celled organisms. With rising populations of drug-resistant strains, the search for novel antimicrobials is an ongoing area of research. Advances in high-performance computing platforms have led to a variety of molecular dynamics simulation strategies to study the interactions of antimicrobial molecules with different compartments of the bacterial cell envelope of both Gram-positive and Gram-negative species. In this review, we begin with a detailed description of the structural aspects of the bacterial cell envelope. Simulations concerned with the transport and associated free energy of small molecules and ions through the outer membrane, peptidoglycan, inner membrane and outer membrane porins are discussed. Since surfactants are widely used as antimicrobials, a section is devoted to the interactions of surfactants with the cell wall and inner membranes. The review ends with a discussion on antimicrobial peptides and the insights gained from the molecular simulations on the free energy of translocation. Challenges involved in developing accurate molecular models and coarse-grained strategies that provide a trade-off between atomic details with a gain in sampling time are highlighted. The need for efficient sampling strategies to obtain accurate free energies of translocation is also discussed. Molecular dynamics simulations have evolved as a powerful tool that can potentially be used to design and develop novel antimicrobials and strategies to effectively treat bacterial infections.

Efficient calculation of orientation-dependent lipid dynamics from membrane simulations

Molecular dynamics simulations of lipid membranes have become increasingly impactful in biophysics because they offer atomistic resolution of structural fluctuations in relation to their functional outputs. Yet quantitative characterization of multiscale processes is a formidable challenge due to the distribution of motions that evade analysis of discrete simulation data. Here we investigate the efficient calculation of CH bond relaxation rates from membrane simulations. Widely used computational approaches offer numerical simplicity but fall short of capturing crucial aspects of the orientation dependence of the dynamics. To circumvent this problem, we introduced a robust framework based on liquid crystal theory which considers explicitly the CH bond motions with respect to the director axis (bilayer normal). Analysis of the orientation dependence of the dynamics shows excellent agreement with experiment, illustrating how the ordering potential affects the calculated relaxation rates. Furthermore, a fit-based resampling of the autocorrelation function of the bond fluctuations validates the new approach for low-temporal resolution data. The recovered relaxation rates indicate that at short timescales, both with and without cholesterol, the local motions of CH bonds describe the bilayer microviscosity and resemble liquid hydrocarbons. Our results establish the critical role of the orientational anisotropy in analysis of membrane simulations, explain fundamental aspects of lipid dynamics, and provide guidelines for extracting information that can be compared to experimental data.

Orientation and Membrane Partition Free Energy of PeT-Based Voltage-Sensitive Dyes from Molecular Simulations

Voltage measurement via small-molecule fluorescent indicators is a valuable approach in deciphering complex dynamics in electrically excitable cells. However, our understanding of various physicochemical properties governing the performance of fluorescent voltage sensors based on the photoinduced electron transfer (PeT) mechanism remains incomplete. Here, through extensive molecular dynamics and free energy calculations, we systematically examine the orientation and membrane partition of three PeT-based voltage-sensing VoltageFluor (VF) dyes in different lipid environment. We show that the symmetry of the molecular scaffold and the net charge of the hydrophilic headgroup of a given VF dye dominate its orientation and membrane partition, respectively. Our work provides a mechanistic understanding of the physical properties contributing to the voltage sensitivity, signal-to-noise ratio, as well as membrane distribution of VF dyes and sheds light onto rational design principles of PeT-based fluorescent probes in general.

Shocker─A Molecular Dynamics Protocol and Tool for Accelerating and Analyzing the Effects of Osmotic Shocks

The process of osmosis, a fundamental phenomenon in life, drives water through a semipermeable membrane in response to a solute concentration gradient across this membrane. In vitro, osmotic shocks are often used to drive shape changes in lipid vesicles, for instance, to study fission events in the context of artificial cells. While experimental techniques provide a macroscopic picture of large-scale membrane remodeling processes, molecular dynamics (MD) simulations are a powerful tool to study membrane deformations at the molecular level. However, simulating an osmotic shock is a time-consuming process due to slow water diffusion across the membrane, making it practically impossible to examine its effects in classic MD simulations. In this article, we present Shocker, a Python-based MD tool for simulating the effects of an osmotic shock by selecting and relocating water particles across a membrane over the course of several pumping cycles. Although this method is primarily aimed at efficiently simulating volume changes in vesicles, it can also handle membrane tubes and double bilayer systems. Additionally, Shocker is force field-independent and compatible with both coarse-grained and all-atom systems. We demonstrate that our tool is applicable to simulate both hypertonic and hypotonic osmotic shocks for a range of vesicular and bilamellar setups, including complex multicomponent systems containing membrane proteins or crowded internal solutions.

Molecular mechanism of calcitriol enhances membrane water permeability

Helicobacter pylori (H. pylori) exhibits a unique membrane lipid composition, including dimyristoyl phosphatidylethanolamine (DMPE) and cholesterol, unlike other Gram-negative bacteria. Calcitriol has antimicrobial activity against H. pylori, but cholesterol enhances antibiotics resistance in H. pylori. This study explored the changes in membrane structure and the molecular mechanisms of cholesterol/calcitriol translocation using well-tempered metadynamics (WT-MetaD) simulations and microsecond conventional molecular dynamics (CMD) simulations. Calcitriol facilitated water transport across the membrane, while cholesterol had the opposite effect. The differing effects might result from the tail 25-hydroxyl group and a wider range of orientations of calcitriol in the DMPE/dimyristoyl phosphatidylglycerol (DMPG) (3:1) membrane. Calcitriol moves across the bilayer center without changing its orientation along the membrane Z-axis, becomes parallel to the membrane surface at the membrane-water interface, and then rotates approximately 90° in this interface. The translocation mechanism of calcitriol is quite different from the flip-flop of cholesterol. Moreover, calcitriol crossed from one layer to another more easily than cholesterol, causing successive perturbations to the hydrophobic core and increasing water permeation. These results improve our understanding of the relationship between cholesterol/calcitriol concentrations and the lipid bilayer structure and the role of lipid composition in water permeation.

How a few help all: cooperative crossing of lipid membranes by COSAN anions

Experimental results show that the presence of a concentration gradient of certain nano-ions (most notably cobaltabisdicarbollide ([o-COSAN]- anions), induce a current across intact artificial phospholipid bilayers in spite of the high Born free energy estimated for these ions. The mechanism underlying this observed translocation of nano-anions across membranes has yet to be determined. Here we show, using molecular dynamics simulations, that the permeation of [o-COSAN]- anions across a lipid bilayer proceeds in a cooperative manner. Single nano-ions can enter the bilayer but permeation is hampered by a free energy barrier of about 8kBT. The interaction between these nano-ions inside a leaflet induces a flip-flop translocation mechanism with the formation of transient, elongated structure inside the membrane. This cooperative flip-flop allows an efficient distribution of [o-COSAN]- anions in both leaflets of the bilayer. These results suggest the existence of a new mechanism for permeation of nano-ions across lipid membranes, relevant for those that have the appropriate self-assembly character.

Dynamic lipid interactions in the plasma membrane Na+,K+-ATPase

The function of ion-transporting Na+,K+-ATPases depends on the surrounding lipid environment in biological membranes. Two established lipid-interaction sites A and B within the transmembrane domain have been observed to induce protein activation and stabilization, respectively. In addition, lipid-mediated inhibition has been assigned to a site C, but with the exact location not experimentally confirmed. Also, possible effects on lipid interactions by disease mutants dwelling in the membrane-protein interface remain relatively uncharacterized. We simulated human Na+,K+-ATPase α1β1FXYD homology models in E1 and E2 states in an asymmetric, multicomponent plasma membrane to determine both wild-type and disease mutant lipid-protein interactions. The simulated wild-type lipid interactions at the established sites A and B were in agreement with experimental results thereby confirming the membrane-protein model system. The less well-characterized, proposed inhibitory site C was dominated by lipids lacking inhibitory properties. Instead, two sites hosting inhibitory lipids were identified at the extracellular side and also a cytoplasmic CHL-binding site that provide putative alternative locations of Na+,K+-ATPase inhibition. Three disease mutations, Leu302Arg, Glu840Arg and Met859Arg resided in the lipid-protein interface and caused drastic changes in the lipid interactions. The simulation results show that lipid interactions to the human Na+,K+-ATPase α1β1FXYD protein in the plasma membrane are highly state-dependent and can be disturbed by disease mutations located in the lipid interface, which can open up for new venues to understand genetic disorders.

Nucleoside Analog Reverse-Transcriptase Inhibitors in Membrane Environment: Molecular Dynamics Simulations

The behavior of four drugs from the family of nucleoside analog reverse-transcriptase inhibitors (zalcitabine, stavudine, didanosine, and apricitabine) in a membrane environment was traced using molecular dynamics simulations. The simulation models included bilayers and monolayers composed of POPC and POPG phospholipids. It was demonstrated that the drugs have a higher affinity towards POPG membranes than POPC membranes due to attractive long-range electrostatic interactions. The results obtained for monolayers were consistent with those obtained for bilayers. The drugs accumulated in the phospholipid polar headgroup region. Two adsorption modes were distinguished. They differed in the degree of penetration of the hydrophilic headgroup region. Hydrogen bonds between drug molecules and phospholipid heads were responsible for adsorption. It was shown that apricitabine penetrated the hydrophilic part of the POPC and POPG membranes more effectively than the other drugs. Van der Waals interactions between S atoms and lipids were responsible for this.

Utilizing Machine Learning to Greatly Expand the Range and Accuracy of Bottom-Up Coarse-Grained Models through Virtual Particles

Coarse-grained (CG) models parametrized using atomistic reference data, i.e., "bottom up" CG models, have proven useful in the study of biomolecules and other soft matter. However, the construction of highly accurate, low resolution CG models of biomolecules remains challenging. We demonstrate in this work how virtual particles, CG sites with no atomistic correspondence, can be incorporated into CG models within the context of relative entropy minimization (REM) as latent variables. The methodology presented, variational derivative relative entropy minimization (VD-REM), enables optimization of virtual particle interactions through a gradient descent algorithm aided by machine learning. We apply this methodology to the challenging case of a solvent-free CG model of a 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipid bilayer and demonstrate that introduction of virtual particles captures solvent-mediated behavior and higher-order correlations which REM alone cannot capture in a more standard CG model based only on the mapping of collections of atoms to the CG sites.

Measuring pico-Newton Forces with Lipid Anchors as Force Sensors in Molecular Dynamics Simulations

Binding forces between biomolecules are ubiquitous in nature but sometimes as weak as a few pico-Newtons (pN). In many cases, the binding partners are attached to biomembranes with the help of a lipid anchor. One important example are glycolipids that promote membrane adhesion through weak carbohydrate-carbohydrate binding between adjacent membranes. Here, we use molecular dynamics (MD) simulations to quantify the forces generated by bonds involving membrane-anchored molecules. We introduce a method in which the protrusion of the lipid anchors from the membrane acts as the force sensor. Our results with two different glycolipids reveal binding forces of up to 20 pN and corroborate the recent notion that carbohydrate-carbohydrate interactions are generic rather than specific.

Formation of lipid raft nanodomains in homogeneous ternary lipid mixture of POPC/DPSM/cholesterol: Theoretical insights

Lipid rafts, in biological membranes, are cholesterol-rich nanodomains that regulate many protein activities and cellular processes. Understanding the formation of the lipid-raft nanodomains helps us elucidate many complex interactions in the cell. In this study, the formation of lipid-raft nanodomains in a ternary palmitoyl-oleoyl-phosphatidylcholine/stearoyl-sphingomyelin/cholesterol (POPC/DPSM/Chol) lipid mixture, the most realistic surrogate model for biological membranes, has been successfully observed for the first time in-silico using microsecond timescale molecular dynamics simulations. The model reveals the formation of cholesterol-induced nanodomains with raft-like characteristics and their underlying mechanism: the cholesterol molecules segregate themselves into cholesterol nanodomains and then enrich the cholesterol-rich domain with sphingomyelin molecules to form a lipid raft thanks to the weak bonding of cholesterol with sphingomyelin. Besides, it is found that the increase in cholesterol concentration enhances the biophysical properties (e.g., bilayer thickness, area per lipid headgroup, and order parameter) of the lipid raft nanodomains. Such findings suggest that the POPC/DPSM/Chol bilayer is a suitable model to fundamentally extend the nanodomain evolution to investigate their lifetime and kinetics as well as to study protein-lipid interaction, protein-protein interaction, and selection of therapeutic molecules in the presence of lipid rafts.

The Structures of Heterogeneous Membranes and Their Interactions with an Anticancer Peptide: A Molecular Dynamics Study

We conduct molecular dynamics simulations of model heterogeneous membranes and their interactions with a 24-amino acid peptide-NAF-144-67. NAF-144-67 is an anticancer peptide that selectively permeates and kills malignant cells; it does not permeate normal cells. We examine three membranes with different binary mixtures of lipids, DOPC-DOPA, DOPC-DOPS, and DOPC-DOPE, with a single peptide embedded in each as models for the diversity of biological membranes. We illustrate that the peptide organization in the membrane depends on the types of nearby phospholipids and is influenced by the charge and size of the head groups. The present study sheds light on early events of permeation and the mechanisms by which an amphiphilic peptide crosses from an aqueous solution to a hydrophobic membrane. Understanding the translocation mechanism is likely to help the design of new permeants.

Understanding the Interactions between Triolein and Cosolvent Binary Mixtures Using Molecular Dynamics Simulations

Biodiesel is one of the emerging renewable sources of energy to replace fossil-fuel-based resources. It is produced by a transesterification reaction in which a triglyceride reacts with methanol in the presence of a catalyst. The reaction is slow because of the low solubility of methanol in triglycerides, which results in low concentrations of methanol available to react with triglyceride. To speed up the reaction, cosolvents are added to create a single phase which helps to improve the concentration of methanol in the triglyceride phase. In this study, molecular dynamics simulations are used to help understand the role of cosolvents in the solvation of triglyceride (triolein). Six binary mixtures of triolein/cosolvent were used to study the solvation of triolein at 298.15 K. Results of 100 ns simulations at constant temperature and pressure to simulate mixing experiments show that in the first 10 ns all the binary mixtures remain largely unmixed. However, for the cosolvents that are fully miscible with triolein, the partial densities across the simulation boxes show that the systems are fully mixed in the final 10 ns. Some solvents were found to interact strongly with the polar part of triolein, while others interacted with the aliphatic part. The radial distribution functions and clustering of the solvents around triolein were also used as indicators for solvation. The presence of cosolvents also influenced the conformation of triolein molecules. In the presence of solvents that solubilize it, triolein preferred a propeller conformation but took up a trident conformation when there is less or no solubilization. The results show that tetrahydrofuran is the best solvent at solubilizing triolein, followed by cyclopentyl methyl ether, diethyl ether, and hexane. With 1,4-dioxane, the solubility improves with an increase in temperature. The miscibility of a solvent in triolein is aided by its ability to interact with both the polar and nonpolar parts of triolein.

All-Atom Modeling of Complex Cellular Membranes

Cell membranes are composed of a variety of lipids and proteins where they interact with each other to fulfill their roles. The first step in modeling these interactions in molecular simulations is to have reliable mimetics of the membrane's lipid environment. This Feature Article presents our recent efforts to model complex cellular membranes using all-atom force fields. A short review of the CHARMM36 (C36) lipid force field and its recent update to incorporate the long-range dispersion is presented. Key examples of model membranes mimicking various species and organelles are given. These include single-celled organisms such as bacteria (E. coli., chlamydia, and P. aeruginosa) and yeast (plasma membrane, endoplasmic reticulum, and trans-Golgi network) and more advanced ones such as plants (soybean and Arabidopsis thaliana) and mammals (ocular lens, stratum corneum, and peripheral nerve myelin). Leaflet asymmetry in composition has also been applied to some of these models. With the increased lipid diversity in the C36 lipid FF, these complex models can better reflect the structural, mechanical, and dynamic properties of realistic membranes and open an opportunity to study biological processes involving other molecules.

Leaflet Asymmetry Modeling in the Lipid Composition of Escherichia coli Cytoplasmic Membranes

Lipid composition asymmetry between leaflets is important to cell function and plays a key role in the "positive inside" rule in transmembrane proteins. In this work, Escherichia coli inner plasma membrane models reflecting this asymmetry have been investigated at the early-log and stationary stages during the bacterial lifecycle using all-atom molecular dynamics simulations. The CHARMM36 lipid force field is used, and selected membrane properties are tested for variations between two leaflets and whole membranes. Our models include bacterial lipids with a cyclopropane moiety on the sn-2 acyl chain in the stationary membrane model. The PE/PG ratio for two leaflets reflects the "positive inside" rule of membrane proteins, set to 6.8 and 2.8 for the inner and outer leaflets of the two models, respectively. We are the first to model leaflet asymmetry in the lipid composition of E. coli cytoplasmic membranes and observe the effect on membrane properties in leaflets and whole membranes. Specifically, our results show that for the stationary phase bilayer, the surface area per lipid (SA/lipid) is larger, the thickness (2D_C and D_B) is smaller, the tilt angle is larger, the tilt modulus is smaller, and the deuterium order parameters (S_CD) of sn-1 and sn-2 tails are lower, compared to the early-log stage. Moreover, the stationary stage bilayer has a positive spontaneous curvature, while the early-log stage has a near flat spontaneous curvature. For leaflet asymmetry, the inner leaflet has a larger SA/lipid, a smaller thickness, a smaller elastic tilt modulus (a larger tilt angle), and lower S_CD, compared to the outer leaflet in both stages. Moreover, an asymmetric membrane involves a lipid tilt and a lateral extension, varying from a reference state of a pre-equilibrium membrane. This work encourages a more profound exploration of leaflet asymmetry in various other membrane models and how this might affect the structure and function of membrane-associated peptides and proteins.

Mechanistic Understanding from Molecular Dynamics in Pharmaceutical Research 2: Lipid Membrane in Drug Design

We review the use of molecular dynamics (MD) simulation as a drug design tool in the context of the role that the lipid membrane can play in drug action, i.e., the interaction between candidate drug molecules and lipid membranes. In the standard "lock and key" paradigm, only the interaction between the drug and a specific active site of a specific protein is considered; the environment in which the drug acts is, from a biophysical perspective, far more complex than this. The possible mechanisms though which a drug can be designed to tinker with physiological processes are significantly broader than merely fitting to a single active site of a single protein. In this paper, we focus on the role of the lipid membrane, arguably the most important element outside the proteins themselves, as a case study. We discuss work that has been carried out, using MD simulation, concerning the transfection of drugs through membranes that act as biological barriers in the path of the drugs, the behavior of drug molecules within membranes, how their collective behavior can affect the structure and properties of the membrane and, finally, the role lipid membranes, to which the vast majority of drug target proteins are associated, can play in mediating the interaction between drug and target protein. This review paper is the second in a two-part series covering MD simulation as a tool in pharmaceutical research; both are designed as pedagogical review papers aimed at both pharmaceutical scientists interested in exploring how the tool of MD simulation can be applied to their research and computational scientists interested in exploring the possibility of a pharmaceutical context for their research.

Design of Lipid-Based Nanocarriers via Cation Modulation of Ethanol-Interdigitated Lipid Membranes

Short-chain alcohols (i.e., ethanol) can induce membrane interdigitation in saturated-chain phosphatidylcholines (PCs). In this process, alcohol molecules intercalate between phosphate heads, increasing lateral separation and favoring hydrophobic interactions between opposing acyl chains, which interpenetrate forming an interdigitated phase. Unraveling mechanisms underlying the interactions between ethanol and model lipid membranes has implications for cell biology, biochemistry, and for the formulation of lipid-based nanocarriers. However, investigations of ethanol-lipid membrane systems have been carried out in deionized water, which limits their applicability. Here, using a combination of small- and wide-angle X-ray scattering, small-angle neutron scattering, and all-atom molecular dynamics simulations, we analyzed the effect of varying CaCl₂ and NaCl concentrations on ethanol-induced interdigitation. We observed that while ethanol addition leads to the interdigitation of bulk phase 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) bilayers in the presence of CaCl₂ and NaCl regardless of the salt concentration, the ethanol-induced interdigitation of vesicular DPPC depends on the choice of cation and its concentration. These findings unravel a key role for cations in the ethanol-induced interdigitation of lipid membranes in either bulk phase or vesicular form.

Rational design of hyperstable antibacterial peptides for food preservation

We describe the design of peptides with properties like thermostability, pH stability, and antibacterial activity against a few bacterial food pathogens. Insights obtained from classical structure-function analysis of natural peptides and their mutants through antimicrobial and enzymatic assays are used to rationally develop a set of peptides. pH and thermostability assays were performed to demonstrate robust antimicrobial activity post-treatment with high temperatures and at wide pH ranges. We have also investigated the mode of action of these hyperstable peptides using membrane permeability assays, electron microscopy, and molecular dynamics simulations. Notably, through mutational studies, we show that these peptides elicit their antibacterial action via both membrane destabilization and inhibition of intracellular trypsin-the two functions attributable to separate peptide segments. Finally, toxicity studies and food preservation assays demonstrate the safety and efficacy of the designed peptides for food preservation. Overall, the study provides a general 'blueprint' for the development of stable antimicrobial peptides (AMPs). Insights obtained from this work may also be combined with combinatorial methods in high-throughput studies for future development of antimicrobials for various applications.

Unsupervised Learning Methods for Molecular Simulation Data

Unsupervised learning is becoming an essential tool to analyze the increasingly large amounts of data produced by atomistic and molecular simulations, in material science, solid state physics, biophysics, and biochemistry. In this Review, we provide a comprehensive overview of the methods of unsupervised learning that have been most commonly used to investigate simulation data and indicate likely directions for further developments in the field. In particular, we discuss feature representation of molecular systems and present state-of-the-art algorithms of dimensionality reduction, density estimation, and clustering, and kinetic models. We divide our discussion into self-contained sections, each discussing a specific method. In each section, we briefly touch upon the mathematical and algorithmic foundations of the method, highlight its strengths and limitations, and describe the specific ways in which it has been used-or can be used-to analyze molecular simulation data.

The ecological roles of microbial lipopeptides: Where are we going?

Lipopeptides (LPs) are secondary metabolites produced by a diversity of bacteria and fungi. Their unique chemical structure comprises both a peptide and a lipid moiety. LPs are of major biotechnological interest owing to their emulsification, antitumor, immunomodulatory, and antimicrobial activities. To date, these versatile compounds have been applied across multiple industries, from pharmaceuticals through to food processing, cosmetics, agriculture, heavy metal, and hydrocarbon bioremediation. The variety of LP structures and the diversity of the environments from which LP-producing microorganisms have been isolated suggest important functions in their natural environment. However, our understanding of the ecological role of LPs is limited. In this review, the mode of action and the role of LPs in motility, antimicrobial activity, heavy metals removal and biofilm formation are addressed. We include discussion on the need to characterise LPs from a diversity of microorganisms, with a focus on taxa inhabiting 'extreme' environments. We introduce the use of computational target fishing and molecular dynamics simulations as powerful tools to investigate the process of interaction between LPs and cell membranes. Together, these advances will provide new understanding of the mechanism of action of novel LPs, providing greater insights into the roles of LPs in the natural environment.

Molecular dynamics simulations and experimental studies reveal differential permeability of withaferin-A and withanone across the model cell membrane

Poor bioavailability due to the inability to cross the cell membrane is one of the major reasons for the failure of a drug in clinical trials. We have used molecular dynamics simulations to predict the membrane permeability of natural drugs-withanolides (withaferin-A and withanone) that have similar structures but remarkably differ in their cytotoxicity. We found that whereas withaferin-A, could proficiently transverse through the model membrane, withanone showed weak permeability. The free energy profiles for the interaction of withanolides with the model bilayer membrane revealed that whereas the polar head group of the membrane caused high resistance for the passage of withanone, the interior of the membrane behaves similarly for both withanolides. The solvation analysis further revealed that the high solvation of terminal O5 oxygen of withaferin-A was the major driving force for its high permeability; it interacted with the phosphate group of the membrane that led to its smooth passage across the bilayer. The computational predictions were tested by raising and recruiting unique antibodies that react to withaferin-A and withanone. The time-lapsed analyses of control and treated cells demonstrated higher permeation of withaferin-A as compared to withanone. The concurrence between the computation and experimental results thus re-emphasised the use of computational methods for predicting permeability and hence bioavailability of natural drug compounds in the drug development process.

Role of polyunsaturated phospholipids in liquid-ordered and liquid-disordered phases

Polyunsaturated phospholipids play a strong repulsive role in the liquid-disordered phase but a weak role in the liquid-ordered phase.

Membrane permeability of small molecules from unbiased molecular dynamics simulations

Permeation of many small molecules through lipid bilayers can be directly observed in molecular dynamics simulations on the nano- and microsecond timescale. While unbiased simulations provide an unobstructed view of the permeation process, their feasibility for computing permeability coefficients depends on various factors that differ for each permeant. The present work studies three small molecules for which unbiased simulations of permeation are feasible within less than a microsecond, one hydrophobic (oxygen), one hydrophilic (water), and one amphiphilic (ethanol). Permeabilities are computed using two approaches: counting methods and a maximum-likelihood estimation for the inhomogeneous solubility diffusion (ISD) model. Counting methods yield nearly model-free estimates of the permeability for all three permeants. While the ISD-based approach is reasonable for oxygen, it lacks precision for water due to insufficient sampling and results in misleading estimates for ethanol due to invalid model assumptions. It is also demonstrated that simulations using a Langevin thermostat with collision frequencies of 1/ps and 5/ps yield oxygen permeabilities and diffusion constants that are lower than those using Nosé-Hoover by statistically significant margins. In contrast, permeabilities from trajectories generated with Nosé-Hoover and the microcanonical ensemble do not show statistically significant differences. As molecular simulations become more affordable and accurate, calculation of permeability for an expanding range of molecules will be feasible using unbiased simulations. The present work summarizes theoretical underpinnings, identifies pitfalls, and develops best practices for such simulations.

Tail-Oxidized Cholesterol Enhances Membrane Permeability for Small Solutes

Cholesterol renders mammalian cell membranes more compact by reducing the amount of voids in the membrane structure. Because of this, cholesterol is known to regulate the ability of cell membranes to prevent the permeation of water and water-soluble molecules through the membranes. Meanwhile, it is also known that even seemingly tiny modifications in the chemical structure of cholesterol can lead to notable changes in membrane properties. The question is, how significantly do these small changes in cholesterol structure affect the permeability barrier function of cell membranes? In this work, we applied fluorescence methods as well as atomistic molecular dynamics simulations to characterize changes in lipid membrane permeability induced by cholesterol oxidation. The studied 7β-hydroxycholesterol (7β-OH-chol) and 27-hydroxycholesterol (27-OH-chol) represent two distinct groups of oxysterols, namely, ring- and tail-oxidized cholesterols, respectively. Our previous research showed that the oxidation of the cholesterol tail has only a marginal effect on the structure of a lipid bilayer; however, oxidation was found to disturb membrane dynamics by introducing a mechanism that allows sterol molecules to move rapidly back and forth across the membrane-bobbing. Herein, we show that bobbing of 27-OH-chol accelerates fluorescence quenching of NBD-lipid probes in the inner leaflet of liposomes by dithionite added to the liposomal suspension. Systematic experiments using fluorescence quenching spectroscopy and microscopy led to the conclusion that the presence of 27-OH-chol increases membrane permeability to the dithionite anion. Atomistic molecular dynamics simulations demonstrated that 27-OH-chol also facilitates water transport across the membrane. The results support the view that oxysterol bobbing gives rise to successive perturbations to the hydrophobic core of the membrane, and these perturbations promote the permeation of water and small water-soluble molecules through a lipid bilayer. The observed impairment of permeability can have important consequences for eukaryotic organisms. The effects described for 27-OH-chol were not observed for 7β-OH-chol which represents ring-oxidized sterols.

EML webinar overview: Simulation-assisted discovery of membrane targeting nanomedicine

The COVID-19 pandemic has brought infectious diseases again to the forefront of global public health concerns. In this EML webinar (Gao, 2020), we discuss some recent work on simulation-assisted discovery of membrane targeting nanomedicine to counter increasing antimicrobial resistance and potential application of similar ideas to the current pandemic. A recent report led by the world health organization (WHO) warned that 10 million people worldwide could die of bacterial infections each year by 2050. To avert the crisis, membrane targeting antibiotics are drawing increasing attention due to their intrinsic advantage of low resistance development. In collaboration with a number of experimental groups, we show examples of simulation-assisted discovery of molecular agents capable of selectively penetrating and aggregating in bacterial lipid membranes, causing membrane permeability/rupture. Through systematic all-atom molecular dynamics simulations and free energy analysis, we demonstrate that the membrane activity of the molecular agents correlates with their ability to enter, perturb and permeabilize the lipid bilayers. Further study on different cell membranes demonstrates that the selectivity results from the presence of cholesterol in mammalian but not in bacterial membranes, as the cholesterol can condense the hydrophobic region of membrane, preventing the penetration of the molecular agents. Following the molecular penetration, we establish a continuum theory and derive the energetic driving force for the domain aggregation and pore growth on lipid membrane. We show that the energy barrier to membrane pore formation can be significantly lowered through molecular aggregation on a large domain with intrinsic curvature and a sharp interface. The theory is consistent with experimental observations and validated with coarse-grained molecular dynamics simulations of molecular domain aggregation leading to pore formation in a lipid membrane. The mechanistic modelling and simulation provide some fundamental principles on how molecular antimicrobials interact with bacterial membranes and damage them through domain aggregation and pore formation. For treating viral infections and cancer therapy, we discuss potential size- and lipid-type-based selectivity principles for developing membrane active nanomedicine. These studies suggest a general simulation-assisted platform to accelerate discovery and innovation in nanomedicine against infectious diseases. EML Webinar speakers are updated at https://imechanica.org/node/24132.

Effects of aquaporin-lipid molar ratio on the permeability of an aquaporin Z-phospholipid membrane system

Aquaporins are water-permeable membrane-channel proteins found in biological cell membranes that selectively exclude ions and large molecules and have high water permeability, which makes them promising candidates for water desalination systems. To effectively apply the properties of aquaporins in the desalination process, many studies have been conducted on aquaporin-lipid membrane systems using phospholipids, which are the main component of cell membranes. Many parametric studies have evaluated the permeability of such systems with various aquaporin types and lipid compositions. In this study, we performed molecular dynamics simulations for four cases with different protein-lipid molar ratios (1:50, 1:75, 1:100, and 1:150) between aquaporin Z and the phospholipids, and we propose a possibility of the existence of optimal protein-lipid molar ratio to maximize water permeability. Elucidating these simulation results from a structural viewpoint suggests that there is a relationship between the permeability and changes in the hydrophobic thickness of the lipid membrane adjacent to the aquaporin as a structural parameter. The results of this study can help optimize the design of an aquaporin-lipid membrane by considering its molar ratio at an early stage of development.

Assessing Barriers for Antimicrobial Penetration in Complex Asymmetric Bacterial Membranes: A Case Study with Thymol

The bacterial cell envelope is a complex multilayered structure evolved to protect bacteria in hostile environments. An understanding of the molecular basis for the interaction and transport of antibacterial therapeutics with the bacterial cell envelope will enable the development of drug molecules to combat bacterial infections and suppress the emergence of drug-resistant strains. Here we report the successful creation of an in vitro supported lipid bilayer (SLB) platform of the outer membrane (OM) of E. coli, an archetypical Gram-negative bacterium, containing the full smooth lipopolysaccharide (S-LPS) architecture of the membrane. Using this platform, we performed fluorescence correlation spectroscopy (FCS) in combination with molecular dynamics (MD) simulations to measure lipid diffusivities and provide molecular insights into the transport of natural antimicrobial agent thymol. Lipid diffusivities measured on symmetric supported lipid bilayers made up of inner membrane lipids show a distinct increase in the presence of thymol as also corroborated by MD simulations. However, lipid diffusivities in the asymmetric OM consisting of only S-LPS are invariant upon exposure to thymol. Increasing the phospholipid content in the LPS-containing outer leaflet improved the penetration toward thymol as reflected in slightly higher relative diffusivity changes in the inner leaflet when compared with the outer leaflet. Free-energy computations reveal the presence of a barrier (∼6 kT) only in the core-saccharide region of the OM for the translocation of thymol while the external O-antigen part is easily traversed. In contrast, thymol spontaneously inserts into the inner membrane. In addition to providing leaflet-resolved penetration barriers in bacterial membranes, we also assess the ability of small molecules to penetrate various membrane components. With rising bacterial resistance, our study opens up the possibility of screening potential antimicrobial drug candidates using these realistic model platforms for Gram-negative bacteria.

The Neutrally Charged Diarylurea Compound PQ401 Kills Antibiotic-Resistant and Antibiotic-Tolerant Staphylococcus aureus

Membrane-damaging antimicrobial agents have great potential to treat multidrug-resistant or multidrug-tolerant bacteria against which conventional antibiotics are not effective. However, their therapeutic applications are often hampered due to their low selectivity to bacterial over mammalian membranes or their potential for cross-resistance to a broad spectrum of cationic membrane-active antimicrobial agents. We discovered that the diarylurea derivative compound PQ401 has antimicrobial potency against multidrug-resistant and multidrug-tolerant Staphylococcus aureus. PQ401 selectively disrupts bacterial membrane lipid bilayers in comparison to mammalian membranes. Unlike cationic membrane-active antimicrobials, the neutral form of PQ401 rather than its cationic form exhibits maximum membrane activity. Overall, our results demonstrate that PQ401 could be a promising lead compound that overcomes the current limitations of membrane selectivity and cross-resistance. Also, this work provides deeper insight into the design and development of new noncharged membrane-targeting therapeutics to combat hard-to-cure bacterial infections.

Resistance or tolerance to traditional antibiotics is a challenging issue in antimicrobial chemotherapy. Moreover, traditional bactericidal antibiotics kill only actively growing bacterial cells, whereas nongrowing metabolically inactive cells are tolerant to and therefore “persist” in the presence of legacy antibiotics. Here, we report that the diarylurea derivative PQ401, previously characterized as an inhibitor of the insulin-like growth factor I receptor, kills both antibiotic-resistant and nongrowing antibiotic-tolerant methicillin-resistant Staphylococcus aureus (MRSA) by lipid bilayer disruption. PQ401 showed several beneficial properties as an antimicrobial lead compound, including rapid killing kinetics, low probability for resistance development, high selectivity to bacterial membranes compared to mammalian membranes, and synergism with gentamicin. In contrast to well-studied membrane-disrupting cationic antimicrobial low-molecular-weight compounds and peptides, molecular dynamic simulations supported by efficacy data demonstrate that the neutral form of PQ401 penetrates and subsequently embeds into bacterial lipid bilayers more effectively than the cationic form. Lastly, PQ401 showed efficacy in both the Caenorhabditis elegans and Galleria mellonella models of MRSA infection. These data suggest that PQ401 may be a lead candidate for repurposing as a membrane-active antimicrobial and has potential for further development as a human antibacterial therapeutic for difficult-to-treat infections caused by both drug-resistant and -tolerant S. aureus.

Carotenoids promote lateral packing and condensation of lipid membranes

Carotenoids are pigment molecules that protect biomembranes against degradation and may be involved in the formation of functional bacterial membrane microdomains. Little is known on whether different types of carotenoids have different effects on the membrane or if there is any concentration dependence of these effects. In this work, we present results from molecular dynamics simulations of phospholipid bilayers containing different amounts of either β-carotene or zeaxanthin. Both β-carotene and zeaxanthin show the ability to laterally condense the membrane lipids and reduce their inter-leaflet interactions. With increasing concentrations, both carotenoids increase the bilayer thickness and rigidity. The results reveal that carotenoids have similar effects to cholesterol on regulating the behavior of fluid-phase membranes, suggesting that they could function as sterol substitutes and confirming their potential role in the formation of functional membrane domains.

Divalent cations bind to phosphoinositides to induce ion and isomer specific propensities for nano-cluster initiation in bilayer membranes

We report all-atom molecular dynamics simulations of asymmetric bilayers containing phosphoinositides in the presence of monovalent and divalent cations. We have characterized the molecular mechanism by which these divalent cations interact with phosphoinositides. Ca²⁺ desolvates more readily, consistent with single-molecule calculations, and forms a network of ionic-like bonds that serve as a 'molecular glue' that allows a single ion to coordinate with up to three phosphatidylinositol-(4,5)-bisphosphate (PI(4, 5)P₂) lipids. The phosphatidylinositol-(3,5)-bisphosphate isomer shows no such effect and neither does PI(4, 5)P₂ in the presence of Mg²⁺. The resulting network of Ca²⁺-mediated lipid-lipid bonds grows to span the entire simulation space and therefore has implications for the lateral distribution of phosophoinositides in the bilayer. We observe context-specific differences in lipid diffusion rates, lipid surface densities and bilayer structure. The molecular-scale delineation of ion-lipid arrangements reported here provides insight into similar nanocluster formation induced by peripheral proteins to regulate the formation of functional signalling complexes on the membrane.

A Minimal Membrane Metal Transport System: Dynamics and Energetics of mer Proteins

The mer operon in bacteria encodes a set of proteins and enzymes that impart resistance to environmental mercury toxicity by importing Hg2+ and reducing it to volatile Hg(0). Because the reduction occurs in the cytoplasm, mercuric ions must first be transported across the cytoplasmic membrane by one of a few known transporters. MerF is the smallest of these, containing only two transmembrane helices and two pairs of vicinal cysteines that coordinate mercuric ions. In this work, we use molecular dynamics simulations to characterize the dynamics of MerF in its apo and Hg2+ -bound states. We find that the apo state positions one of the cysteine pairs closer to the periplasmic side of the membrane, while in the bound state the same pair approaches the cytoplasmic side. This finding is consistent with the functional requirement of accepting Hg2+ from the periplasmic space, sequestering it on acceptance, and transferring it to the cytoplasm. Conformational changes in the TM helices facilitate the functional interaction of the two cysteine pairs. Free-energy calculations provide a barrier of 16 kcal/mol for the association of the periplasmic Hg2+ -bound protein MerP with MerF and 7 kcal/mol for the subsequent association of MerF's two cysteine pairs. Despite the significant conformational changes required to move the binding site across the membrane, coarse-grained simulations of multiple copies of MerF support the expectation that it functions as a monomer. Our results demonstrate how conformational changes and binding thermodynamics could lead to such a small membrane protein acting as an ion transporter. Published 2019. This article is a U.S. Government work and is in the public domain in the USA.

A molecular dynamics study proposing the existence of statistical structural heterogeneity due to chain orientation in the POPC-cholesterol bilayer

We performed molecular dynamics simulations of a lipid bilayer consisting of POPC and cholesterol at temperatures from 283 to 308K and cholesterol concentrations from 0 to 50% mol/mol. The purpose of this study was to look for the existence of structural differences in the region delimited by these parameters and, in particular, in a region where coexistence of liquid disordered and liquid ordered phases has been proposed. Our interest in this range of concentration and temperature responds to the fact that polyene ionophore activity varies considerably along it. Two force fields, CHARMM36 and Slipids, were compared in order to determine the most suitable. Both force fields predict non-monotonic behaviors consistent with the existence of phase transitions. We found the presence of lateral structural heterogeneity, statistical in nature, in some of the bilayers occurring in this range of temperatures and sterol concentrations. This heterogeneity was produced by correlated ordering of the POPC tails and not due to cholesterol enrichment, and lasts for tens of nanoseconds. We relate these observations to the action of polyenes in these membranes.

Interaction of Antimicrobial Lipopeptides with Bacterial Lipid Bilayers

The resistance of pathogens to traditional antibiotics is currently a global issue of enormous concern. As the discovery and development of new antibiotics become increasingly challenging, synthetic antimicrobial lipopeptides (AMLPs) are now receiving renewed attention as a new class of antimicrobial agents. In contrast to traditional antibiotics, AMLPs act by physically disrupting the cell membrane (rather than targeting specific proteins), thus reducing the risk of inducing bacterial resistance. In this study, we use microsecond-timescale atomistic molecular dynamics simulations to quantify the interaction of a short AMLP (C16-KKK) with model bacterial lipid bilayers. In particular, we investigate how fundamental transmembrane properties change in relation to a range of lipopeptide concentrations. A number of structural, mechanical, and dynamical features are found to be significantly altered in a non-linear fashion. At 10 mol% concentration, lipopeptides have a condensing effect on bacterial bilayers, characterized by a decrease in the area per lipid and an increase in the bilayer order. Higher AMLP concentrations of 25 and 40 mol% destabilize the membrane by disrupting the bilayer core structure, inducing membrane thinning and water leakage. Important transmembrane properties such as the lateral pressure and dipole potential profiles are also affected. Potential implications on membrane function and associated proteins are discussed.

Antimicrobial peptide ROAD-1 triggers phase change in local membrane environment to execute its activity

Emergence of antibiotic-resistant pathogens has paved way for development of newer class of drugs that would not be susceptible to resistance. Antimicrobial peptides such as defensins that target the microbial membrane are promising candidates. ROAD-1 is an alpha-defensin present in the oral cavity of rhesus macaque and shares very high sequence similarity to human enteric defensin 5. In this study we have performed microsecond long all atom molecular dynamic simulations to understand the mechanism of action of ROAD-1. We find that ROAD-1 is able to adopt an energetically stable conformation predominantly stabilized by electrostatic interactions only in presence of bacterial membranes. In mammalian membrane even though it gets absorbed onto the bilayer, it is unable to adopt an equilibrium conformation. Binding of ROAD-1 to bilayer induces clustering of POPG molecules up to 15 Å around the peptide. POPG molecules show higher order parameters than the neighboring POPE implying coexistence of different phases. Analysis of binding free energy of ROAD-1-membrane complex indicates Arg1, Arg2, Arg7, and Arg25 to play key role in its antimicrobial activity. Unlike its homolog HD5, ROAD-1 is not observed to form a dimer. Our study gives insight into the membrane-bound conformation of ROAD-1 and its mechanism of action that can aid in designing defensin-based therapeutics. Graphical abstract Antimicrobial peptide ROAD-1 adopts a different membrane-bound conformation as compared with HD5 even though they belong to the same family implying a different mechanism of action.

Accelerating Membrane Simulations with Hydrogen Mass Repartitioning

The time step of atomistic molecular dynamics (MD) simulations is determined by the fastest motions in the system and is typically limited to 2 fs. An increasingly popular approach is to increase the mass of the hydrogen atoms to ∼3 amu and decrease the mass of the parent atom by an equivalent amount. This approach, known as hydrogen-mass repartitioning (HMR), permits time steps up to 4 fs with reasonable simulation stability. While HMR has been applied in many published studies to date, it has not been extensively tested for membrane-containing systems. Here, we compare the results of simulations of a variety of membranes and membrane-protein systems run using a 2 fs time step and a 4 fs time step with HMR. For pure membrane systems, we find almost no difference in structural properties, such as area-per-lipid, electron density profiles, and order parameters, although there are differences in kinetic properties such as the diffusion constant. Conductance through a porin in an applied field, partitioning of a small peptide, hydrogen-bond dynamics, and membrane mixing show very little dependence on HMR and the time step. We also tested a 9 Å cutoff as compared to the standard CHARMM cutoff of 12 Å, finding significant deviations in many properties tested. We conclude that HMR is a valid approach for membrane systems, but a 9 Å cutoff is not.

A selective membrane-targeting repurposed antibiotic with activity against persistent methicillin-resistant Staphylococcus aureus

Treatment of Staphylococcus aureus infections is complicated by the development of antibiotic tolerance, a consequence of the ability of S. aureus to enter into a nongrowing, dormant state in which the organisms are referred to as persisters. We report that the clinically approved anthelmintic agent bithionol kills methicillin-resistant S. aureus (MRSA) persister cells, which correlates with its ability to disrupt the integrity of Gram-positive bacterial membranes. Critically, bithionol exhibits significant selectivity for bacterial compared with mammalian cell membranes. All-atom molecular dynamics (MD) simulations demonstrate that the selectivity of bithionol for bacterial membranes correlates with its ability to penetrate and embed in bacterial-mimic lipid bilayers, but not in cholesterol-rich mammalian-mimic lipid bilayers. In addition to causing rapid membrane permeabilization, the insertion of bithionol increases membrane fluidity. By using bithionol and nTZDpa (another membrane-active antimicrobial agent), as well as analogs of these compounds, we show that the activity of membrane-active compounds against MRSA persisters positively correlates with their ability to increase membrane fluidity, thereby establishing an accurate biophysical indicator for estimating antipersister potency. Finally, we demonstrate that, in combination with gentamicin, bithionol effectively reduces bacterial burdens in a mouse model of chronic deep-seated MRSA infection. This work highlights the potential repurposing of bithionol as an antipersister therapeutic agent.

Nanotoxicity of Boron Nitride Nanosheet to Bacterial Membranes

Boron nitride (BN) nanosheet is a novel material with great potential in biomedical applications. A deep understanding of the basic interaction mechanisms between biosystems and foreign BN nanosheets can help to better clarify the potential risks of these nanomaterials and provide guidance on their safe design. In this paper, we show that BN nanosheets can cause degradation of bacterial cell membranes via experimental and simulation-based approaches. Our extensive molecular dynamics simulations results reveal that BN nanosheets cause toxicity to both bacterial outer and inner membranes in which hydrophobic effect plays an important role. The spontaneous lipid extraction by BN nanosheets is in agreement with the free-energy calculations. A liquid-to-gel phase transition is induced by the BN nanosheet in the outer model membrane of bacteria, indicating that the BN nanosheet may cause higher toxicity to the outer membrane than to the inner membrane. Our findings may offer new insights into the molecular basis of BN's cytotoxicity and antibacterial activity.

Physical properties of model biological lipid bilayers: insights from all-atom molecular dynamics simulations

The physical properties of lipid bilayers are sensitive to the specific type and composition of the lipids that make up the many different types of cell membranes. Studying model bilayers of representative heterogeneous compositions can provide key insights into membrane functionality. In this work, we use atomistic molecular dynamics simulations to characterize key properties in a number of bilayer membranes of varying composition. We first examine basic properties, such as lipid area, volume, and bilayer thickness, of simple, homogeneous bilayers comprising several lipid types, which are prevalent in biological membranes. Such lipids are then used in simulations of heterogeneous systems representative of bacterial, mammalian, and cancer membranes. Our analysis is especially focused on depth-dependent, transmembrane profiles; in particular, we calculate lateral pressure and dipole potential profiles, two fundamental properties which play key roles in a large number of biological functions.

Molecular-Level "Observations" of the Behavior of Gold Nanoparticles in Aqueous Solution and Interacting with a Lipid Bilayer Membrane

We use coarse-grained molecular dynamics simulations to "observe" details of interactions between ligand-covered gold nanoparticles and a lipid bilayer model membrane. In molecular dynamics simulations, one puts the individual atoms and groups of atoms of the physical system to be "observed" into a simulation box, specifies the forms of the potential energies of interactions between them (ultimately quantum based), and lets them individually move classically according to Newton's equations of motion, based on the forces arising from the assumed potential energy forms. The atoms that are chemically bonded to each other stay chemically bonded, following known potentials (force fields) that permit internal degrees of freedom (internal rotation, torsion, vibrations), and the interactions between nonbonded atoms are simplified to Lennard-Jones forms (in our case) and coulombic (where electrical charges are present) in which the parameters are previously optimized to reproduce thermodynamic properties or are based on quantum electronic calculations. The system is started out at a reasonable set of coordinates for all atoms or groups of atoms, and then permitted to develop according to the equations of motion, one small step (usually 10 fs time step) at a time, for millions of steps until the system is at a quasi-equilibrium (usually reached after hundreds of nanoseconds). We then let the system play out its motions further for many nanoseconds to observe the behavior, periodically taking snapshots (saving all positions and energies), and post-processing the snapshots to obtain various average descriptions of the system. Alkanethiols of various lengths serve as examples of hydrophobic ligands and methyl-terminated PEG with various numbers of monomer units serve as examples of hydrophilic ligands. Spherical gold particles of various diameters as well as gold nanorods form the core to which ligands are attached. The nanoparticles are characterized at the molecular level, especially the distributions of ligand configurations and their dependence on ligand length, and surface coverage. Self-assembly of the bilayer from an isotropic solution and observation of membrane properties that correspond well to experimental values validate the simulations. The mechanism of permeation of a gold NP coated with either a hydrophobic or a hydrophilic ligand, and its dependence on surface coverage, ligand length, core diameter, and core shape, is investigated. Lipid response such as lipid flip-flops, lipid extraction, and changes in order parameter of the lipid tails are examined in detail. The mechanism of permeation of a PEGylated nanorod is shown to occur by tilting, lying down, rotating, and straightening up. The nature of the information provided by molecular dynamics simulations permits understanding of the detailed behavior of gold nanoparticles interacting with lipid membranes which in turn helps to understand why some known systems work better than others and aids the design of new particles and improvement of methods for preparing existing ones.