Molecular simulation study of structural and dynamic properties of mixed DPPC/DPPE bilayers

Protonation of palmitic acid embedded in DPPC lipid bilayers obscures detection of ripple phase by FTIR spectroscopy

The transformation of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid bilayers from the gel (Lβ') to the fluid (Lα) phase involves an intermediate ripple (Pβ') phase forming a few degrees below the main transition temperature (Tm). While the exact cause of bilayer rippling is still debated, the presence of amphiphilic molecules, pH, and lipid bilayer architecture are all known to influence (pre)transition behavior. In particular, fatty acid chains interact with hydrophobic lipid tails, while the carboxylic groups simultaneously participate in proton transfer with interfacial water in the polar lipid region which is controlled by the pH of the surrounding aqueous medium. The molecular-level variations in the DPPC ripple phase in the presence of 2% palmitic acid (PA) were studied at pH levels 4.0, 7.3, and 9.1, where PA is fully protonated, partially protonated, or fully deprotonated. Bilayer thermotropic behavior was investigated by differential scanning calorimetry (DSC) and Fourier-transform infrared (FTIR) spectroscopy which agreed in their characterization of (pre)transition at pH of 9.1, but not at pH 4.0 and especially not at 7.3. Owing to the different insertion depths of protonated and deprotonated PA, along with the ability of protonated PA to undergo flip-flop in the bilayer, these two forms of PA show a different hydration pattern in the interfacial water layer. Finally, these results demonstrated the hitherto undiscovered potential of FTIR spectroscopy in the detection of the events occurring at the surface of lipid bilayers that obscure the low-cooperativity phase transition explored in this work.

Membrane Composition Allows the Optimization of Berberine Encapsulation in Liposomes

Berberine (BBR) is a natural molecule with noteworthy pharmacological properties, including the prevention of antibiotic resistance in Gram-negative bacteria. However, its oral bioavailability is poor, thus resulting in an impaired absorption and efficacy in humans. In combination with other drugs, liposomes have been shown to enhance the availability of the drug, representing a smart delivery system to target tissues and reduce negative side effects. To date, there is a lack of studies on BBR and liposomes that enable the rationalization and molecular-based design of such formulations for future use in humans. In this work, the encapsulation of BBR into liposomes is proposed to overcome current limitations using a combination of experimental and computational assays to rationalize the membrane composition of liposomes that maximizes BBR encapsulation. First, the encapsulation efficiency was measured for several membrane compositions, revealing that it is enhanced by cholesteryl hemisuccinate and, to a lesser extent, by cholesterol. The physical basis of the BBR encapsulation efficiency and permeability was clarified using molecular dynamics simulation: using the lipid composition, one can tune the capability of membranes to attract, i.e., to adsorb, the molecules onto their surface. Overall, these findings suggest a rational strategy to maximize the encapsulation efficiency of liposomes by using negatively charged lipids, thus representing the basis for designing delivery systems for BBR, useful to treat, e.g., antibiotic resistance.

Optical and molecular features of negatively curved surfaces created by POPE lipids: A crucial role of the initial conditions

Membrane fusion is closely related to plasma membrane domains rich in cone-shaped phosphatidylethanolamine (PE) lipids that can reverse membrane curvature under certain conditions. The phase transition of PE-based lipid membranes from the lamellar fluid phase (Lα) to the inverse hexagonal phase (HII) is commonly taken as a general model in reconstructing the membrane fusion pathway, and whose structural features have been mostly described so far using structural and microscopic techniques. The aim of this paper is to decipher the optical and molecular features of Lβ → Lα and especially of Lα → HII transition of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) lipids at pH = 7.0 when they are initially prepared in the form of both multi- and unilamellar liposomes (MLVs and LUVs). The distinction between optical properties of MLS- and LUVs-derived HII phase, provided from turbidity-sensitive temperature-dependent UV-Vis spectra, was attributed to different formation mechanisms of HII phase. Most importantly, from FTIR spectroscopic data of POPE lipids in Lβ (15 °C), Lα (50 °C) and HII (85 °C) phases we identified the changes in molecular features of POPE lipids during phase transitions. Among the latter, by far the most significant is different hydration pattern of POPE lipids in MLVs- and LUVs-derived HII phase which extends from the polar-apolar interface all the way to the terminal amino group of the POPE lipid, along with the changes in the conformation of glycerol backbone as evidenced by the signature of α-methylene groups. Molecular dynamics simulations confirmed higher water penetration in HII phase and provided insight into hydrogen bonding patterns.

Molecular dynamics simulations reliably identify vibrational modes in far-IR spectra of phospholipids

The properties of self-assembled phospholipid membranes are of essential importance in biochemistry and physical chemistry, providing a platform for many cellular life functions. Far-infrared (far-IR) vibrational spectroscopy, on the other hand, is a highly information-rich method to characterize intermolecular interactions and collective behaviour of lipids that can help explain, e.g., chain packing, thermodynamic phase behaviour, and sequestration. However, reliable interpretation of the far-IR spectra is still lacking. Here we present a molecular dynamics (MD) based approach to simulate vibrational modes of individual lipids and in an ensemble. The results are a good match to synchrotron far-IR measurements and enable identification of the molecular motions corresponding to each vibrational mode, thus allowing the correct interpretation of membrane spectra with high accuracy and resolving the longstanding ambiguities in the literature in this regard. Our results demonstrate the feasibility of using MD simulations for interpreting far-IR spectra broadly, opening new avenues for practical use of this powerful method.

Molecular simulation on the interaction between trehalose and asymmetric lipid bilayer mimicking the membrane of human red blood cells

Trehalose is widely acknowledged for its ability to stabilize plasma membranes during dehydration. However, the exact mechanism by which trehalose interacts with lipid bilayers remains presently unclear. In this study, we conducted atomistic molecular dynamic simulations on asymmetric model bilayers that mimic the membrane of human red blood cells at various trehalose and water contents. We considered three different hydration levels mimicking the full hydration to desiccation scenarios. Results indicate that the asymmetric distribution of lipids did not significantly influence the computed structural characteristics at full and low hydration. At dehydration, however, the order parameter obtained from the symmetric bilayer is significantly higher compared to those obtained from asymmetric ones. Analysis of hydrogen bonds revealed that the protective ability of trehalose is well described by the water replacement hypothesis at full and low hydration, while at dehydration other interaction mechanisms associated with trehalose exclusion from the bilayer may involve. In addition, we found that trehalose exclusion is not attributed to sugar saturation but rather to the reduction in hydration levels. It can be concluded that the protective effect of trehalose is not only related to the hydration level of the bilayer, but also closely tied to the asymmetric distribution of lipids within each leaflet.

Susceptibility of Legionella gormanii Membrane-Derived Phospholipids to the Peptide Action of Antimicrobial LL-37-Langmuir Monolayer Studies

LL-37 is the only member of the cathelicidin-type host defense peptide family in humans. It exhibits broad-spectrum bactericidal activity, which represents a distinctive advantage for future therapeutic targets. The presence of choline in the growth medium for bacteria changes the composition and physicochemical properties of their membranes, which affects LL-37's activity as an antimicrobial agent. In this study, the effect of the LL-37 peptide on the phospholipid monolayers at the liquid-air interface imitating the membranes of Legionella gormanii bacteria was determined. The Langmuir monolayer technique was employed to prepare model membranes composed of individual classes of phospholipids-phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylglycerol (PG), cardiolipin (CL)-isolated from L. gormanii bacteria supplemented or non-supplemented with exogenous choline. Compression isotherms were obtained for the monolayers with or without the addition of the peptide to the subphase. Then, penetration tests were carried out for the phospholipid monolayers compressed to a surface pressure of 30 mN/m, followed by the insertion of the peptide into the subphase. Changes in the mean molecular area were observed over time. Our findings demonstrate the diversified effect of LL-37 on the phospholipid monolayers, depending on the bacteria growth conditions. The substantial changes in membrane properties due to its interactions with LL-37 enable us to propose a feasible mechanism of peptide action at a molecular level. This can be associated with the stable incorporation of the peptide inside the monolayer or with the disruption of the membrane leading to the removal (desorption) of molecules into the subphase. Understanding the role of antimicrobial peptides is crucial for the design and development of new strategies and routes for combating resistance to conventional antibiotics.

Monolayer/Bilayer Equilibrium of Phospholipids in Gel or Liquid States: Interfacial Adsorption via Monomer or Liposome Diffusion?

Phospholipids (PLs) are widely used in the pharma industry and a better understanding of their behavior under different conditions is helpful for applications such as their use as medical transporters. The transition temperature Tm affects the lipid conformation and the interfacial tension between perfluoroperhydrophenanthrene (PFP) and an aqueous suspension of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC), 1,2-distearoyl-sn-glycero-3-phosphatidylcholine (DSPC), as well as a mixture of these PLs with cholesterol. Interfacial tensions were measured with the Du Noüy ring at quasi-equilibrium; the area per molecule was calculated according to the Gibbsian approach and a time-dependent tension gradient. Results show that the time tε to reach quasi-equilibrium was shorter when the temperature was above Tm, indicating a faster adsorption process (tε,DPPC,36 °C = 48 h, tε,DPPC,48 °C = 24 h) for PL in the liquid crystalline state than in the gel state (T < Tm). In addition, concentration-dependent results of the interfacial tension revealed that above the respective Tm and at all concentrations c > 0.1 mM, the average minimum interfacial tension for DPPC and DSPC (14.1 mN/m and 15.3 mN/m) does not differ significantly between those two lipids. Equilibrium between monolayers and bilayers shows that for T < Tm, surface pressures ∏ ≈ 31 mN/m are reached while for T > Tm, ∏ ≈ 41 mN/m. Mixtures with cholesterol only reach ∏ ≤ 31 mN/m Tm, with no significant difference between the two PLs. The higher interfacial tension of the mixture indicates stabilization of the liposomal conformation in the aqueous phase by the addition of cholesterol. The high diffusion coefficients show that adsorption is mainly based on liposomes.

Lipid doping of the sponge (L3) mesophase

The polymorphism of lipid aggregates has long attracted detailed study due to the myriad factors that determine the final mesophase observed. This study is driven by the need to understand mesophase behaviour for a number of applications, such as drug delivery and membrane protein crystallography. In the case of the latter, the role of the so-called 'sponge' (L3) mesophase has been often noted, but not extensively studied by itself. The L3 mesophase can be formed in monoolein/water systems on the addition of butanediol to water, which partitions the headgroup region of the membrane, and decreases its elastic moduli. Like cubic mesophases, it is bicontinuous, but unlike them, has no long-range translational symmetry. In our present study, we show that the formation of the L3 phase can delicately depend on the addition of dopant lipids to the mesophase. While electrostatically neutral molecules similar in shape to monoolein (DOPE, cholesterol) have little effect on the general mesophase behaviour, others (DOPC, DDM) significantly reduce the composition at which it can form. Additionally, we show that by combining cholesterol with the anionic lipid DOPG, it is possible to form the largest stable L3 mesophases observed to date, with characteristic lengths over 220 Å.

Implications of different membrane compartmentalization models in particle-based in silico studies

Studying membrane dynamics is important to understand the cellular response to environmental stimuli. A decisive spatial characteristic of the plasma membrane is its compartmental structure created by the actin-based membrane-skeleton (fences) and anchored transmembrane proteins (pickets). Particle-based reaction-diffusion simulation of the membrane offers a suitable temporal and spatial resolution to analyse its spatially heterogeneous and stochastic dynamics. Fences have been modelled via hop probabilities, potentials or explicit picket fences. Our study analyses the different approaches' constraints and their impact on simulation results and performance. Each of the methods comes with its own constraints; the picket fences require small timesteps, potential fences might induce a bias in diffusion in crowded systems, and probabilistic fences, in addition to carefully scaling the probability with the timesteps, induce higher computational costs for each propagation step.

Physicochemical Characteristics of Model Membranes Composed of Legionella gormanii Lipids

Legionella gormanii is one of the species belonging to the genus Legionella, which causes atypical community-acquired pneumonia. The most important virulence factors that enable the bacteria to colonize the host organism are associated with the cell surface. Lipids building the cell envelope are crucial not only for the membrane integrity of L. gormanii but also by virtue of being a dynamic site of interactions between the pathogen and the metabolites supplied by its host. The utilization of exogenous choline by the Legionella species results in changes in the lipids' composition, which influences the physicochemical properties of the cell surface. The aim of this study was to characterize the interfacial properties of the phospholipids extracted from L. gormanii cultured with (PL+choline) and without exogenous choline (PL-choline). The Langmuir monolayer technique coupled with the surface potential (SPOT) sensor and the Brewster angle microscope (BAM) made it possible to prepare the lipid monomolecular films (model membranes) and study their properties at the liquid/air interface at 20 °C and 37 °C. The results indicate the effect of the choline addition to the bacterial medium on the properties of the L. gormanii phospholipid membranes. The differences were revealed in the organization of monolayers, their molecular packing and ordering, degree of condensation and changes in the components' miscibility. These findings are the basis for further research on the mechanisms of adaptation of this pathogen, which by changing the native composition and properties of lipids, bypasses the action of antimicrobial compounds and avoids the host immune attack.

The role of hydrophobic patches of de novo designed MSI-78 and VG16KRKP antimicrobial peptides on fragmenting model bilayer membranes

Antimicrobial peptides (AMPs) with cell membrane lysing capability are considered potential candidates for the development of the next generation of antibiotics. Designing novel AMPs requires an in-depth understanding of the mechanism of action of the peptides. In this work, we used various biophysical techniques including ³¹P solid-state NMR to examine the interaction of model membranes with amphipathic de novo-designed peptides. Two such peptides, MSI-78 and VG16KRKP, were designed with different hydrophobicity and positive charges. The model lipid membranes were constituted by mixing lipids of varying degrees of 'area per lipid' (APL), which directly affected the packing properties of the membrane. The observed emergence of the isotropic peak in ³¹P NMR spectra as a function of time is a consequence of the fragmentation of the membrane mediated by the peptide interaction. The factors such as the charges, overall hydrophilicity of the AMPs, as well as lipid membrane packing, contributed to the kinetics of membrane fragmentation. Furthermore, we anticipate the designed AMPs follow the carpet and toroidal pore mechanisms when lysing the cell membrane. This study highlights the significance of the effect of the overall charges and the hydrophobicity of the novel AMPs designed for antimicrobial activity.

The reconciliation between the experimental and calculated octanol-water partition coefficient of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine using atomistic molecular dynamics: an open question

The octanol-water partition coefficient of 1,2-dipalmitoyl-sn-glycero-3-phosphatidylcholine (DPPC) was investigated using atomistic molecular dynamics simulations via thermodynamic integration and multistate Bennett acceptance ratio methods. The GAFF and CHARMM36 force fields were used with six water models widely used in molecular dynamics simulations. The OPC4 water model provided the best agreement with the experimental octanol-water partition coefficient of DPPC using the two force fields. However, there is still plenty of room for improvement in water models with correct estimation of surface tension that uses better and suitable non-bonded interaction parameters between water-water and water-DPPC. The Gibbs free energy of transferring DPPC from octanol to water phase was calculated to be 19.8 ± 0.3 and 20.2 ± 0.3 kcal mol-1, giving a partition coefficient of 14.5 ± 0.4 and 14.8 ± 0.3 for the GAFF and CHARMM36 force fields, respectively. This study reinforces the importance of developing new water models that reproduce experimental surface tensions to reconcile the water-water and water-DPPC non-bonded interactions and the existing discrepancy between experimental measurements of amphiphilic molecules that are important in many areas of scientific applications and industry such as biophysics, surfactant, colloids, membranes, medicine, nanotechnology, and food and pharmaceutical industries, and so on. It raises two important open questions: Is the experimental octanol-water partition coefficient of DPPC reliable? Or is its calculation accurate using the OPC4 water model? With respect to the experimental measurements, there may be non-treated aspects such as the formation of aggregates in aqueous phase and limit of detection of the applied method. And, in the calculation, some effects are not possible to be considered in a correct way or viable time such as calculating quantum effects, sampling all conformations, considering phase transitions, and correctly evaluating the intermolecular forces to estimate an accurate surface tension.Communicated by Ramaswamy H. Sarma.

Study of the Rv1417 and Rv2617c Membrane Proteins and Their Interactions with Nicotine Derivatives as Potential Inhibitors of Erp Virulence-Associated Factor in Mycobacterium tuberculosis: An In Silico Approach

The increasing emergence of Mycobacterium tuberculosis (Mtb) strains resistant to traditional anti-tuberculosis drugs has alarmed health services worldwide. The search for new therapeutic targets and effective drugs that counteract the virulence and multiplication of Mtb represents a challenge for the scientific community. Several studies have considered the erp gene a possible therapeutic target in the last two decades, since its disruption negatively impacts Mtb multiplication. This gene encodes the exported repetitive protein (Erp), which is located in the cell wall of Mtb. In vitro studies have shown that the Erp protein interacts with two putative membrane proteins, Rv1417 and Rv2617c, and the impairment of their interactions can decrease Mtb replication. In this study, we present five nicotine analogs that can inhibit the formation of heterodimers and trimers between these proteins. Through DFT calculations, molecular dynamics, docking, and other advanced in silico techniques, we have analyzed the molecular complexes, and show the effect these compounds have on protein interactions. The results show that four of these analogs can be possible candidates to counteract the pathogenicity of Mtb. This study aims to combine research on the Erp protein as a therapeutic target in the search for new drugs that serve to create new therapies against tuberculosis disease.

Influence of DPPE surface undulations on melting temperature determination: UV/Vis spectroscopic and MD study

One of the most distinguished quantities that describes lipid main phase transition, i.e. the transition from the gel (L_β(')) to the fluid (L_α) phase, is its melting temperature (T_m). Because melting is accompanied by a large change in enthalpy the, L_β(') → L_α transition can be monitored by various calorimetric, structural and spectroscopic techniques and T_m should be the same regardless of the metric monitored or the technique employed. However, in the case of DPPE multilamellar aggregates there is a small but systematic deviation of T_m values determined by DSC and FTIR spectroscopy. The aim of this paper is to explain this discrepancy by combined UV/Vis spectroscopic and MD computational approach. Multivariate analysis performed on temperature-dependent UV/Vis spectra of DPPE suspensions demonstrated that at 55 ± 1 °C certain phenomenon causes a small but detectable change in suspension turbidity, whereas a dominant change in the latter is registered at 63.2 ± 0.4 °C that coincides with T_m value determined from DSC curve. If this effect should be ignored, the overall data give T_m value the same as FTIR spectra data (61.0 ± 0.4 °C). As the classical MD simulations suggest that about 10° below T_m certain undulations appear at the surface of DPPE bilayers, we concluded that certain discontinuities in curvature fluctuations arise at reported temperature which are to some extent coupled with lipid melting. Ultimately, such events and the associated changes in curvature affect T_m value measured by different techniques.

Plasmonic lipogels: driving co-assembly of composites with peptide-based gels for controlled drug release

Supramolecular short peptide-based gels are promising materials for the controlled release of drugs (e.g. chemotherapeutic drugs) owing to the biocompatibility and similarity to cell matrix. However, the drug encapsulation and control over its release, mainly the hydrophilic drugs, can be a cumbersome task. This can be overcome through encapsulation/compartmentalization of drugs in liposomes, which can also enable spatiotemporal control and enhanced drug release through a trigger, such as photothermia. Having this in mind, we explored the assembly of silica-coated gold nanoparticles and liposomes (storage units) with dehydropeptide-based hydrogels as a proof-of-concept to afford peptide-based NIR light-responsive lipogels. Several liposomes compositions were assessed that displayed influence on the final assembly properties by combining with silica-coated gold nanorods (∼106 nm). Gold nanospheres (∼11 nm) were used to study the preparation method, which revealed the importance of initially combine liposomes with nanoparticles and then the gelator solution to achieve a closer proximity of the nanoparticles to the liposomes. The control over a hydrophilic model drug, 5(6)-carboxyfluorescein, was only achieved by its encapsulation in liposomes, in which the presence of silica-coated nanorods further enabled the use of photothermia to induce the liposomes phase transition and stimulate the drug release. Further, both composites, the liposomes and silica-coated gold nanorods, induced a lower elastic modulus, but also provided an enhanced gelation kinetics. Hereby, this work advances fabrication strategies for the development of short peptide-based hydrogels towards on-demand, sustained and controlled release of hydrophilic drugs through photothermia under NIR light irradiation.

Homeoviscous Adaptation of the Lipid Membrane of a Soil Bacterium Surviving under Diurnal Temperature Variation: A Molecular Simulation Perspective

A recent experiment has reported the lipidome remodeling of a soil-based plant-associated bacterium Methylobacterium extorquens due to diurnal temperature variations. The key adaptation strategy is the headgroup-specific remodeling of the acyl chain. To understand the idiosyncratic adaptation at the molecular level, we simulate the model membrane of the same bacterium using the reported lipidome compositions at four different experimental temperatures. We investigate the temperature-dependent packing density and fluidity of the membrane, the constancy of which is key to the homeoviscous adaptation. The results show that complex lipidome remodeling approximately preserves membrane properties under heat and cold stress. The headgroup-specific remodeling of the acyl chain serves to fine-tune the packing density and fluidity of the membrane at different temperatures. While lipids with strongly interacting headgroups are more abundant at higher temperatures, the lipidome is more dominated by lipids with weaker interacting headgroups at lower temperatures. This adaptation alleviates lipid membrane disruption caused by heat and cold stress. This study provides a molecular picture of the homeoviscous adaptation of the realistic lipid membrane of a soil-based bacterium.

Novel In Silico Insights into Rv1417 and Rv2617c as Potential Protein Targets: The Importance of the Medium on the Structural Interactions with Exported Repetitive Protein (Erp) of Mycobacterium tuberculosis

Nowadays, tuberculosis is the second leading cause of death from a monopathogenic transmitted disease, only ahead of COVID-19. The role of exported repetitive protein (Erp) in the virulence of Mycobacterium tuberculosis has been extensively demonstrated. In vitro and in vivo assays have identified that Erp interacts with Rv1417 and Rv2617c proteins, forming putative transient molecular complexes prior to localization to the cell envelope. Although new insights into the interactions and functions of Erp have emerged over the years, knowledge about its structure and protein–protein interactions at the atomistic level has not been sufficiently explored. In this work, we have combined several in silico methodologies to gain new insights into the structural relationship between these proteins. Two system conditions were evaluated by MD simulations: Rv1417 and Rv2617c embedded in a lipid membrane and another with a semi-polar solvent to mimic the electrostatic conditions on the membrane surface. The Erp protein was simulated as an unanchored structure. Stabilized structures were docked, and complexes were evaluated to recognize the main residues involved in protein–protein interactions. Our results show the influence of the medium on the structural conformation of proteins. Globular conformations were favored under high polarity conditions and showed a higher energetic affinity in complex formation. Meanwhile, disordered conformations were favored under semi-polar conditions and an increase in the number of contacts between residues was observed. In addition, the electrostatic potential analysis showed remarkable changes in protein interactions due to the polarity of the medium, demonstrating the relevance of Erp protein in heterodimer formation. On the other hand, contact analysis showed that several C-terminal residues of Erp were involved in the protein interactions, which seems to contradict experimental observations; however, these complexes could be transient forms. The findings presented in this work are intended to open new perspectives in the studies of Erp protein molecular interactions and to improve the knowledge about its function and role in the virulence of Mycobacterium tuberculosis.

How Does Temperature Affect the Dynamics of SARS-CoV-2 M Proteins? Insights from Molecular Dynamics Simulations

Enveloped viruses, in general, have several transmembrane proteins and glycoproteins, which assist the virus in entry and attachment onto the host cells. These proteins also play a significant role in determining the shape and size of the newly formed virus particles. The lipid membrane and the embedded proteins affect each other in non-trivial ways during the course of the viral life cycle. Unraveling the nature of the protein-protein and protein-lipid interactions, under various environmental and physiological conditions, could therefore prove to be crucial in development of therapeutics. Here, we study the M protein of SARS-CoV-2 to understand the effect of temperature on the properties of the protein-membrane system. The membrane-embedded dimeric M proteins were studied using atomistic and coarse-grained molecular dynamics simulations at temperatures ranging between 10 and 50 °C. While temperature-induced fluctuations are expected to be monotonic, we observe a steady rise in the protein dynamics up to 40 °C, beyond which it surprisingly reverts back to the low-temperature behavior. Detailed investigation reveals disordering of the membrane lipids in the presence of the protein, which induces additional curvature around the transmembrane region. Coarse-grained simulations indicate temperature-dependent aggregation of M protein dimers. Our study clearly indicates that the dynamics of membrane lipids and integral M protein of SARS-CoV-2 enables it to better associate and aggregate only at a certain temperature range (i.e., ~ 30-40 °C). This can have important implications in the protein aggregation and subsequent viral budding/fission processes.

Interaction of Docetaxel with Phosphatidylcholine Membranes: A Combined Experimental and Computational Study

The antineoplastic drug Docetaxel is a second generation taxane which is used against a great variety of cancers. The drug is highly lipophilic and produces a great array of severe toxic effects that limit its therapeutic effectiveness. The study of the interaction between Docetaxel and membranes is very scarce, however, it is required in order to get clues in relation with its function, mechanism of toxicity and possibilities of new formulations. Using phosphatidylcholine biomimetic membranes, we examine the interaction of Docetaxel with the phospholipid bilayer combining an experimental study, employing a series of biophysical techniques like Differential Scanning Calorimetry, X-Ray Diffraction and Infrared Spectroscopy, and a Molecular Dynamics simulation. Our experimental results indicated that Docetaxel incorporated into DPPC bilayer perturbing the gel to liquid crystalline phase transition and giving rise to immiscibility when the amount of the drug is increased. The drug promotes the gel ripple phase, increasing the bilayer thickness in the fluid phase, and is also able to alter the hydrogen-bonding interactions in the interfacial region of the bilayer producing a dehydration effect. The results from computational simulation agree with the experimental ones and located the Docetaxel molecule forming small clusters in the region of the carbon 8 of the acyl chain palisade overlapping with the carbonyl region of the phospholipid. Our results support the idea that the anticancer drug is embedded into the phospholipid bilayer to a limited amount and produces structural perturbations which might affect the function of the membrane.

New spirit of an old technique: Characterization of lipid phase transitions via UV/Vis spectroscopy

One of the advantages of investigating lipid phase transitions by thermoanalytical techniques such as DSC is manifested in the proportionality of the signal strength on a DSC curve, attributed to a particular thermotropic event, and its cooperativity degree. Accordingly, the pretransition of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) is less noticeable than its main phase transition; as a matter of fact, when DSC measurements are performed at low heating rate, such low-cooperativity phase transition could go (almost) unnoticed. The aim of this work is to present temperature-dependent UV/Vis spectroscopy, based on a temperature-dependent change in DPPC suspension turbidity, as a technique applicable for determination of lipid phase transition temperatures. Multivariate analyzes of the acquired UV/Vis spectra show that phase transitions of the low-cooperativity degree, such as pretransitions, can be identified with the same certainty as transitions of a high-cooperativity degree.

Interactions between Beta-2-Glycoprotein-1 and Phospholipid Bilayer-A Molecular Dynamic Study

This study aims to investigate the interactions appearing when the beta-2-glycoprotein-1 binds to a lipid bilayer. The inter- and intra-molecular forces acting between the two macromolecular systems have been investigated using a molecular dynamics simulation method. The importance of water bridges has also been addressed. Additionally, the viscoelastic response of the bilayer has been studied. In detail, the (saturated-chain) 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and (unsaturated-chain) 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) bilayers have been chosen to test their behavior near the protein. Both of the lipids have a polar head but different chemical structures and are similar to the main phospholipids present in the synovial fluid. This study is meaningful for further explaining the worsening friction properties in articular cartilage, as the inactivation of phospholipid bilayers by beta-2-glycoprotein-1 is believed to be a cause of the destruction of cartilage in most rheumatic diseases and osteoarthritis. It was found that the protein binds stronger to the DPPC bilayer than to the POPE, but in both cases, it has the potential to change the local bilayer stability. Nevertheless, the binding forces are placed within a small area (only a few lipids contribute to the binding, creating many interactions). However, together, they are not stronger than the covalent bonds between C-O, thus, potentially, it is possible to push the lipids into the bilayer but detaching the lipids' heads from the tail is not possible. Additionally, the protein causes water displacement from the vicinity of the bilayer, and this may be a contributor to the instability of the bilayer (disrupting the water bridges needed for the stabilization of the bilayer, especially in the case of DPPC where the heads are not so well stabilized by H-bonds as they are in POPE). Moreover, it was found that the diffusivity of lipids in the DPPC bilayer bound to the protein is significantly different from the diffusivity of the ones which are not in contact with the protein. The POPE bilayer is stiffer due to intramolecular interactions, which are stronger than in the DPPC; thus, the viscous to elastic effects in the POPE case are more significant than in the case of the DPPC. It is, therefore, harder to destabilize the POPE bilayer than the DPPC one.

Dissecting the Molecular Mechanism of Colistin Resistance in mcr-1 Bacteria

Colistin or polymyxin B is the last resort antibiotic to treat infections of multidrug-resistant Gram-negative bacteria by disrupting their outer membranes. The recent emergence of Gram-negative bacteria that demonstrate colistin resistance, particularly plasmid-mediated mobile colistin resistance (mcr), poses a big challenge to the treatment of multidrug resistance infections. Using molecular dynamics simulations, we explore the mechanism of colistin resistance in a model lipid A bilayer mimicking the Gram-negative mcr-1 bacterial outer membrane. The simulation results reveal that the outer membrane of normal Gram-negative bacteria is stabilized by salt bridges between positively charged divalent ions and negatively charged phosphate groups of the membranes. In the presence of positively charged polymyxin B, these salt bridges are disrupted, and calcium is released into the aqueous phase, resulting in membrane disruption. In contrast, the lipid A in the outer membrane of mcr-1 bacteria has a novel modification, this being a covalently attached phosphoethanolamine group. This group enables the formation of a large number of hydrogen bonds between the amine and phosphate groups, resulting in an electrostatic net on the membrane. This extensive noncovalent electrostatic cross-linking between the lipid molecules collectively enhances the membrane stability and results in resistance to the action of cationic peptides such as polymyxin B. The simulation results shed new atomistic insights for understanding the mechanistic basis of colistin resistance and provide clues for the design of new membrane disruptors and permeabilizers to treat mcr-1 infections.

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.

Weak Fragment Crystallizable (Fc) Domain Interactions Drive the Dynamic Assembly of IgG Oligomers upon Antigen Recognition

Activation of membrane receptors through clustering is a common mechanism found in various biological systems. Spatial proximity of receptors may be transduced across the membrane to initiate signaling pathways or alternatively be recognized by peripheral proteins or immune cells to trigger external effector functions. Here we show how specific immunoglobulin G (IgG) binding induces clustering of monomeric target molecules in lipid membranes through Fc-Fc interactions. We visualize and characterize the dynamic IgG oligomerization process and the molecular interactions involved using high-speed atomic force microscopy, single-molecule force spectroscopy, and quartz crystal microbalance experiments. We found that the Fc-Fc interaction strength is precisely tuned to be weak enough to prevent IgG oligomerization in solution at physiological titers, but enabling IgG oligomerization when Fabs additionally bind to their cognate surface epitopes, a mechanism that ultimately targets IgG-mediated effector functions such as classical complement activation to antigenic membranes.

Adsorption of impurities in vegetable oil: A molecular modelling study

Here, the adsorption of impurity species from triglyceride solvent representing a model vegetable oil is studied using atomistic molecular dynamics simulations. We compare the adsorption of water, glycerol, oleic acid, monoolein, and two types of phospholipids on model silica adsorbents differing in their OH-group density, i.e. hydrogen bonding ability, quartz and cristobalite. We find that the species containing charged groups, phospholipids DOPC and DOPE, adsorb significantly stronger than the nonionic impurities. Secondary contribution to adsorption arises from hydrogen bonding capability of the impurity species, the silica surface, and also the triglyceride solvent: in general, more hydrogen bonding sites in impurity species leads to enhanced adsorption but hydrogen bonding with solvent competes for the available sites. Interestingly, adsorption is weaker on cristobalite even though it has a higher hydrogen bonding site density than quartz. This is because the hydrogen bonds can saturate each other on the adsorbent. The finding demonstrates that optimal adsorption response is obtained with intermediate adsorbent hydrogen bonding site densities. Additionally, we find that monoolein and oleic acid show a concentration driven adsorption response and reverse micelle like aggregate formation in bulk triglyceride solvent even in the absence of water. The findings offer insight into adsorption phenomena at inorganic adsorbent - apolar solvent interfaces and provide guidelines for enhanced design of adsorbent materials for example for vegetable oil purification.

Effect of dipole moment on amphiphile solubility and partition into liquid ordered and liquid disordered phases in lipid bilayers

Association of amphiphiles with biomembranes is important for their availability at specific locations in organisms and cells, being critical for their biological function. A prominent role is usually attributed to the hydrophobic effect, and to electrostatic interactions between charged amphiphiles and lipids. This work explores a closely related and complementary aspect, namely the contribution made by dipole moments to the strength of the interactions established. Two xanthene amphiphiles with opposite relative orientations of their dipole and amphiphilic moments have been selected (Rhodamine-C₁₄ and Carboxyfluorescein-C₁₄). The membranes studied have distinct lipid compositions, representing typical cell membrane pools, ranging from internal membranes to the outer and inner leaflet of the plasma membrane. A comprehensive study is reported, including the affinity of the amphiphiles for the different membranes, the stability of the amphiphiles as monomers and their tendency to form small clusters, as well as their transverse location in the membrane. The orientation of the amphiphile dipole moment, which determines whether its interaction with the membrane dipole potential is repulsive or attractive, is found to exert a large influence on the association of the amphiphile with ordered lipid membranes. These interactions are also responsible for the formation of small clusters or stabilization of amphiphile monomers in the membrane. The results obtained allow understanding the prevalence of protein lipidation at the N-terminal for efficient targeting to the plasma membrane, as well as the tendency of GPI-anchored proteins (usually lipidated at the C-terminal) to form small clusters in the membrane ordered domains.

Direct interaction of fibrinogen with lipid microparticles modulates clotting kinetics and clot structure

Extensive studies revealed the role of blood lipid microparticles (liposomes, microvesicles) in activation of coagulation cascade. The direct interaction of fibrinogen/fibrin with lipid surfaces and its consequence for hemostasis received much less attention. We observed pronounced changes in both clot morphology and kinetics of fibrin clotting in the presence of artificial liposomes. The evidence was obtained that lipid microparticles per se present a diffusion barrier to the three-dimensional fibril assembling and pose spatial restrictions for fiber elongation. On the other hand, fibrinogen adsorption results in its high local concentration on liposome surface that accelerates fibrin polymerization. Adsorption induces Fg secondary structure alterations which may contribute to the abnormal clot morphology. In dependence on lipid composition and size of microparticles, the interplay of all the outlined mechanisms determines functionally important changes of clot morphology. The obtained results contribute to the knowledge of clotting mechanisms in the presence of artificial and natural lipid microparticles.

Interplay of Nanoparticle Rigidity and Its Translocation Ability through Cell Membrane

Understanding the endocytic process of nanoparticles (NPs) with different mechanical rigidities is critical to develop effective drug delivery vectors. Here, we perform experiments, coarse-grained molecular dynamics simulations, and theoretical analyses to investigate the role of NPs' mechanical rigidity in the cellular endocytic process. Experiments based on two types of engineered Au NPs that have similar properties but different rigidities are performed in order to investigate their cellular uptake efficiencies, and it has been found that the more rigid NPs can achieve a higher cellular uptake efficiency. Simulation results confirm that rigid NPs can achieve full internalization by forming a complete double-layer endosome coating, while relatively soft NPs can only reach 40% surface coverage by membrane lipids. Simulation results capture an intriguing translocation of multiple NPs with different rigidities in a cooperative manner where the NPs' mechanical rigidities regulate their translocation efficiencies. We find that theoretically rigid NPs require less energy to overcome the energy barrier for membrane internalization than soft NPs do, which is in good agreement with experiment and simulation results. This synergetic study offers useful insight into the design principle of a general NP-based drug delivery vector as well as the promising biomedical application of NP-based medicine.

Contribution of headgroup and chain length of glycerophospholipids to thermal stability and permeability of liposomes loaded with calcein

Biological membranes are complex systems that are composed of lipids, proteins and carbohydrates. They are difficult to study, so it is established practice to use lipid vesicles that consist of closed 'shells' of phospholipid bilayers as model systems to study various functional and structural aspects of lipid organisation. To define the effects of the structural properties of lipid vesicles on their phase behaviour, we investigated their headgroup and chain length, and the chemical bonds by which their acyl chains are attached to the glycerol moiety of glycerophospholipid species, in terms of phase transition temperature, enthalpy change and calcein permeability. We used differential scanning calorimetry to measure the temperature and enthalpy changes of phase transition, and fluorescence to follow calcein release through the bilayer structure. Our data show that longer acyl chains increase the stability of the lipid bilayers, whereas higher salt concentrations decrease the thermal stability and widen the phase transitions of these lipid bilayers. We discuss the possible reasons for the observed phase transition behaviour.

Microcavity-Supported Lipid Bilayers; Evaluation of Drug-Lipid Membrane Interactions by Electrochemical Impedance and Fluorescence Correlation Spectroscopy

Many drugs have intracellular or membrane-associated targets, thus understanding their interaction with the cell membrane is of value in drug development. Cell-free tools used to predict membrane interactions should replicate the molecular organization of the membrane. Microcavity array-supported lipid bilayer (MSLB) platforms are versatile biophysical models of the cell membrane that combine liposome-like membrane fluidity with stability and addressability. We used an MSLB herein to interrogate drug-membrane interactions across seven drugs from different classes, including nonsteroidal anti-inflammatories: ibuprofen (Ibu) and diclofenac (Dic); antibiotics: rifampicin (Rif), levofloxacin (Levo), and pefloxacin (Pef); and bisphosphonates: alendronate (Ale) and clodronate (Clo). Fluorescence lifetime correlation spectroscopy (FLCS) and electrochemical impedance spectroscopy (EIS) were used to evaluate the impact of drug on 1,2-dioleyl- sn-glycerophosphocholine and binary bilayers over physiologically relevant drug concentrations. Although FLCS data revealed Ibu, Levo, Pef, Ale, and Clo had no impact on lipid lateral mobility, EIS, which is more sensitive to membrane structural change, indicated modest but significant decreases to membrane resistivity consistent with adsorption but weak penetration of drugs at the membrane. Ale and Clo, evaluated at pH 5.25, did not impact the impedance of the membrane except at concentrations exceeding 4 mM. Conversely, Dic and Rif dramatically altered bilayer fluidity, suggesting their translocation through the bilayer, and EIS data showed that resistivity of the membrane decreased substantially with increasing drug concentration. Capacitance changes to the bilayer in most cases were insignificant. Using a Langmuir-Freundlich model to fit the EIS data, we propose R_sat as an empirical value that reflects permeation. Overall, the data indicate that Ibu, Levo, and Pef adsorb at the interface of the lipid membrane but Dic and Rif interact strongly, permeating the membrane core modifying the water/ion permeability of the bilayer structure. These observations are discussed in the context of previously reported data on drug permeability and log P.

Molecular dynamics simulation study of the effect of halothane on mixed DPPC/DPPE phospholipid membranes

We report results of a molecular dynamics simulation study of the effect of one general anesthetic, halothane, on some properties of mixed DPPC/DPPE phospholipid membranes. This is a suitable model for the study of simple, two-phospholipid membrane systems. From the simulation runs, we determined several membrane properties for five different molecular proportions of DPPC/DPPE. The effect of halothane on the studied membrane properties (area per lipid molecule, density of membrane, order parameter, etc.) was rather small. The distribution of halothane is not uniform through the bilayer thickness. Instead, there is a maximum of anesthetic concentration around 1.2 nm from the center of the membrane. The anesthetic molecule is located close to the phospholipid headgroups. The position of the halothane density maximum depends slightly on the DPPC/DPPE molar proportion. Snapshots taken over the plane of the membrane, as well as calculated two-dimensional radial distribution functions show that the anesthetic has no preference for either phospholipid (DPPC or DPPE). Our results indicate that this anesthetic molecule has only small effects on DPPC/DPPE mixed membranes. In addition, halothane displays no preferential location around DPPC or DPPE. This is probably due to the hydrophobic nature of halothane and to the fact that the chosen phospholipids have the same hydrophobic tails.

Exploring Dynamics and Structure of Biomolecules, Cryoprotectants, and Water Using Molecular Dynamics Simulations: Implications for Biostabilization and Biopreservation

Successful stabilization and preservation of biological materials often utilize low temperatures and dehydration to arrest molecular motion. Cryoprotectants are routinely employed to help the biological entities survive the physicochemical and mechanical stresses induced by cold or dryness. Molecular interactions between biomolecules, cryoprotectants, and water fundamentally determine the outcomes of preservation. The optimization of assays using the empirical approach is often limited in structural and temporal resolution, whereas classical molecular dynamics simulations can provide a cost-effective glimpse into the atomic-level structure and interaction of individual molecules that dictate macroscopic behavior. Computational research on biomolecules, cryoprotectants, and water has provided invaluable insights into the development of new cryoprotectants and the optimization of preservation methods. We describe the rapidly evolving state of the art of molecular simulations of these complex systems, summarize the molecular-scale protective and stabilizing mechanisms, and discuss the challenges that motivate continued innovation in this field.

Phenylalanine intercalation parameters for liquid-disordered phase domains - a membrane model study

BACKGROUND

Propensity of phenylalanine (Phe) for nonpolar environments drives its intercalation into phospholipid membranes, which has been implicated in metabolic and neurological disorders. The knowledge of Phe intercalation parameters can be instrumental in understanding various membrane processes triggered by interactions with Phe, in particular the early events leading to the formation of nucleation/docking sites for the self-assembly of Phe amyloid fibrils at the membrane surface.

RESULTS

In this study, we used monolayers of phosphatidylethanolamine (DPPE) and phosphatidylcholine (DPPC) to mimic the membrane outer leaflet. Its initial interaction with Phe was modeled by injecting Phe into the aqueous phase underneath the monolayer. Constant pressure insertion assays augmented with epifluorescence microscopy imaging were used to monitor Phe intercalation. Our primary goal was to determine the Phe intercalation area, A Phe. Two values were obtained for A Phe, 33 ± 2 and 48 ± 3 Å2.

CONCLUSIONS

Phe appeared to discriminate between DPPE and DPPC packing, and use two modes of intercalation. The area of A Phe 33 ± 2 Å2 is consistent with a Phe monomer intercalating into membrane by inserting the phenyl ring nearly normal to the membrane surface. This mode has been found to operate in DPPE membranes. For DPPC membranes however, the value of A Phe = 48 ± 3 Å2 suggests that, from saline, Phe tends to intercalate as a larger species plausibly dragging along a counterion, Na+, in a Na+(Phe) complex.

Thermodynamic Characterization of Mixed Monolayers of a Novel Oxazolidine Derivative and Phospholipids

Oxazolidine derivatives (OxD) are five ring-membered compounds that contain at least one oxygen and nitrogen in their molecular structure. OxD are known due to several therapeutic activities such as anticancer and antibiotic properties. In this paper, we performed a thermodynamic analysis of the mixed films composed by dipalmitoylphosphatidylglycerol (DPPG), dipalmitoylphosphoethanolamine (DPPE), dipalmitoyl phosphatidylcholine (DPPC) or L-α phosphatidylcholine (PC) with a novel oxazolidine derivate (OxD). Relevant thermodynamic parameters such as excess areas (ΔAE), excess free energies (ΔG), and Gibbs free energy of mixing (AG^mix) were derived from the surface pressure data. The topographical analysis was performed using atomic force microscopy. Based on the calculated values of the thermodynamic parameters, we observed that the miscibility of the mixed films was directly dependent on their composition. DPPG/OxD and DPPE/OxD systems present the best-mixed character at low pressures at OxD molar fraction equivalent to 0.25.

Molecular dynamics modeling of Pseudomonas aeruginosa outer membranes

Pseudomonas aeruginosa is a common Gram-negative bacterium and opportunistic human pathogen. The distinctive structure of its outer membrane (OM) and outer membrane vesicles (OMVs) plays a fundamental role in bacterial virulence, colonization ability, and antibiotic resistance. To provide critical insights into OM and OMV functionality, we conducted an all-atom molecular dynamics study of asymmetric membranes that are biologically relevant to P. aeruginosa. We hybridized a GLYCAM06-based lipopolysaccharides force field with the Stockholm lipids force field (Slipids) to model bilayer membranes with Lipid A molecules in one leaflet and physiologically relevant phospholipid molecules in the other, including 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine (DPPE), 1,2-dioleoyl-sn-glycero-3-phosphoethanolamine (DOPE), 1,2-dipalmitoyl-sn-glycero-3-phosphoglycerol (DPPG), and 1,2-dioleoyl-sn-glycero-3-phosphoglycerol (DOPG). In particular, a membrane with phospholipid composition representing the P. aeruginosa OM was constructed and modeled by mixing the physiologically dominant components. The detailed structure of membranes was characterized by area per lipid, transmembrane mass and charge densities, radial distribution function (RDF), deuterium order parameter (SCD) of acyl chains, and inclination angles of phosphates and disaccharide in Lipid A. The membrane fluidity in equilibrium and the hydration of functional groups were probed and characterized quantitatively. The consistent properties of the Lipid A leaflets in different membranes demonstrate its compatibility with various phospholipids present in the P. aeruginosa OM. The more ordered acyl chains of Lipid A compared to the cytoplasmic cell membrane contribute to the low permeability of bacterial outer membrane. The findings of this computational investigation of P. aeruginosa OM will further the understanding of microbial pathogenesis and enable future study of OMV biogenesis.

Structure Dependence of Pyridine and Benzene Derivatives on Interactions with Model Membranes

Pyridine-based small-molecule drugs, vitamins, and cofactors are vital for many cellular processes, but little is known about their interactions with membrane interfaces. These specific membrane interactions of these small molecules or ions can assist in diffusion across membranes or reach a membrane-bound target. This study explores how minor differences in small molecules (isoniazid, benzhydrazide, isonicotinamide, nicotinamide, picolinamide, and benzamide) can affect their interactions with model membranes. Langmuir monolayer studies of dipalmitoylphosphatidylcholine (DPPC) or dipalmitoylphosphatidylethanolamine (DPPE), in the presence of the molecules listed, show that isoniazid and isonicotinamide affect the DPPE monolayer at lower concentrations than the DPPC monolayer, demonstrating a preference for one phospholipid over the other. The Langmuir monolayer studies also suggest that nitrogen content and stereochemistry of the small molecule can affect the phospholipid monolayers differently. To determine the molecular interactions of the simple N-containing aromatic pyridines with a membrane-like interface, ¹H one-dimensional NMR and ¹H-¹H two-dimensional NMR techniques were utilized to obtain information about the position and orientation of the molecules of interest within aerosol-OT (AOT) reverse micelles. These studies show that all six of the molecules reside near the AOT sulfonate headgroups and ester linkages in similar positions, but nicotinamide and picolinamide tilt at the water-AOT interface to varying degrees. Combined, these studies demonstrate that small structural changes of small N-containing molecules can affect their specific interactions with membrane-like interfaces and specificity toward different membrane components.

Adsorption of phospholipids at oil/water interfaces during emulsification is controlled by stress relaxation and diffusion

Adsorption of phosphatidylcholines at oil/water interfaces strongly deviates from spread monolayers at air/water surfaces. Understanding its nature and consequences could vastly improve applications in medical nanoemulsions and biotechnologies. Adsorption kinetics at interfaces of water with different oil phases were measured by profile analysis tensiometry. Adsorption kinetics for 2 different phospholipids, DPPC and POPC, as well as 2 organic phases, squalene and squalane, show that formation of interfacial monolayers is initially dominated by stress-relaxation in the first minutes. Diffusion only gradually contributes to a decrease in interfacial tension at later stages of time and higher film pressures. The results can be applied for the optimization of emulsification protocols using mechanical treatments. Emulsions using phospholipids with unsaturated fatty acids are dominated much more strongly by stress-relaxation and cover interfaces very fast compared to those with saturated fatty acids. In contrast, phospholipid layers consisting of saturated fatty acids converge faster towards the equilibrium than those with unsaturated fatty acids.

Computational Insight Into the Structural Organization of Full-Length Toll-Like Receptor 4 Dimer in a Model Phospholipid Bilayer

Toll-like receptors (TLRs) are a unique category of pattern recognition receptors that recognize distinct pathogenic components, often utilizing the same set of downstream adaptors. Specific molecular features of extracellular, transmembrane (TM), and cytoplasmic domains of TLRs are crucial for coordinating the complex, innate immune signaling pathway. Here, we constructed a full-length structural model of TLR4-a widely studied member of the interleukin-1 receptor/TLR superfamily-using homology modeling, protein-protein docking, and molecular dynamics simulations to understand the differential domain organization of TLR4 in a membrane-aqueous environment. Results showed that each functional domain of the membrane-bound TLR4 displayed several structural transitions that are biophysically essential for plasma membrane integration. Specifically, the extracellular and cytoplasmic domains were partially immersed in the upper and lower leaflets of the membrane bilayer. Meanwhile, TM domains tilted considerably to overcome the hydrophobic mismatch with the bilayer core. Our analysis indicates an alternate dimerization or a potential oligomerization interface of TLR4-TM. Moreover, the helical properties of an isolated TM dimer partly agree with that of the full-length receptor. Furthermore, membrane-absorbed or solvent-exposed surfaces of the toll/interleukin-1 receptor domain are consistent with previous X-ray crystallography and biochemical studies. Collectively, we provided a complete structural model of membrane-bound TLR4 that strengthens our current understanding of the complex mechanism of receptor activation and adaptor recruitment in the innate immune signaling pathway.

Simulations of simple Bovine and Homo sapiens outer cortex ocular lens membrane models with a majority concentration of cholesterol

The lipid composition of bovine and human ocular lens membranes has been probed, and a variety of lipids have been found including phosphatidylcholine (PC), phosphatidylethanolamine (PE), sphingomyelin (SM), and cholesterol (CHOL) with cholesterol being present in particularly high concentrations. In this study, we use the all-atom CHARMM36 force field to simulate binary, ternary, and quaternary mixtures as models of the ocular lens. High concentration of cholesterol, in combination with different and varying diversity of phospholipids (PL) and sphingolipids (SL), affect the structure of the ocular lens lipid bilayer. The following analyses were done for each simulation: surface area per lipid, component surface area per lipid, deuterium order parameters (S_CD), electron density profiles (EDP), membrane thickness, hydrogen bonding, radial distribution functions, clustering, and sterol tilt angle distribution. The S_CD show significant bilayer alignment and packing in cholesterol-rich bilayers. The EDP show the transition from liquid crystalline to liquid ordered with the addition of cholesterol. Hydrogen bonds in our systems show the tendency for intramolecular interactions between cholesterol and fully saturated lipid tails for less complex bilayers. But with an increased number of components in the bilayer, the acyl chain of the lipids becomes a less important characteristic, and the headgroup of the lipid becomes more significant. Overall, cholesterol is the driving force of membrane structure of the ocular lens membrane where interactions between cholesterol, PL, and SL determine structure and function of the biomembrane. The goal of this work is to develop a baseline for further study of more physiologically realistic ocular lens lipid membranes. This article is part of a Special Issue entitled: Emergence of Complex Behavior in Biomembranes edited by Marjorie Longo.

Three-Dimensional Microphase Separation and Synergistic Permeability in Stacked Lipid-Polymer Hybrid Membranes

We present new structures of soft-material thin films that augment the functionality of substrate-mediated delivery systems. A hybrid material composed of phospholipids and block copolymers adopts a multilayered membrane structure supported on a solid surface. The hybrid films comprise intentional intramembrane heterogeneities that register across multilayers. These stacked domains convey unprecedented enhancement and control of permeability of solutes across micrometer-thick films. Using grazing incidence X-ray scattering, phase contrast atomic force microscopy, and confocal microscopy, we observed that in each lamella, lipid and polymers partition unevenly within the membrane plane segregating into lipid- or polymer-rich domains. Interestingly, we found evidence that like-domains align in registry across multilayers, thereby making phase separation three-dimensional. Phase boundaries exist over extended length scales to compensate the height mismatch between lipid and polymer molecules. We show that microphase separation in hybrid films can be exploited to augment the capability of drug-eluting substrates. Lipid-polymer hybrid films loaded with paclitaxel show synergistic permeability of drug compared to single-component counterparts. We present a thorough structural study of stacked lipid-polymer hybrid membranes and propose that the presence of registered domains and domain boundaries impart enhanced drug release functionality. This work offers new perspectives in designing thin films for controlled delivery applications.

The folding, stability and function of lactose permease differ in their dependence on bilayer lipid composition

Lipids play key roles in Biology. Mechanical properties of the lipid bilayer influence their neighbouring membrane proteins, however it is unknown whether different membrane protein properties have the same dependence on membrane mechanics, or whether mechanics are tuned to specific protein processes of the protein. We study the influence of lipid lateral pressure and electrostatic effects on the in vitro reconstitution, folding, stability and function of a representative of the ubiquitous major facilitator transporter superfamily, lactose permease. Increasing the outward chain lateral pressure in the bilayer, through addition of lamellar phosphatidylethanolamine lipids, lowers lactose permease folding and reconstitution yields but stabilises the folded state. The presence of phosphatidylethanolamine is however required for correct folding and function. An increase in headgroup negative charge through the addition of phosphatidylglycerol lipids favours protein reconstitution but is detrimental to topology and function. Overall the in vitro folding, reconstitution, topology, stability and function of lactose permease are found to have different dependences on bilayer composition. A regime of lipid composition is found where all properties are favoured, even if suboptimal. This lays ground rules for rational control of membrane proteins in nanotechnology and synthetic biology by manipulating global bilayer properties to tune membrane protein behaviour.

Controlled release of nitric oxide from liposomes

We report the ability to readily tune NO release from N-diazeniumdiolate-encapsulated liposomal structures by altering the NO donor molecule structure and/or phospholipid composition (independently or in combination). While encapsulating more stable NO donors expectedly enhanced the NO release (up to 48 h) from the liposomes, the phospholipid headgroup surface area proved equally useful in controlling NO-release kinetics by influencing the water uptake and concomitant N-diazeniumdiolate NO donor breakdown (to NO). The potential therapeutic utility of the NO-releasing liposomes was further assessed in biological/proteinaceous fluids. The NO-release kinetics were similar in buffer and serum.