Structure of fully hydrated fluid phase lipid bilayers with monounsaturated chains

Insights into short chain polyethylene penetration of phospholipid bilayers via atomistic molecular dynamics simulations

The escalation of global plastic production, reaching an annual output of 400 million tons, has significantly intensified concerns regarding plastic waste management. This has been exacerbated by improper recycling and disposal practices, contributing to the impending crisis of plastic pollution. Predictions indicate that by 2025, the environment will bear the burden of over ten billion metric tons of accumulated plastic waste. This situation has led to the concerning release of microplastics and nanoplastics (NPs) into the environment as plastic materials degrade, thereby posing risks to both ecosystems and human health. Nanoparticle interactions with living organisms have garnered significant attention due to their potential to disrupt vital biological processes. Of particular interest are lipid membranes, acting as crucial gatekeepers, underscoring the importance of comprehending the intricate process of NP penetration. Molecular dynamics (MD) simulations serve as a robust tool, offering molecular-level insights into these intricate interactions. In this study, we leverage all-atom MD simulations to delve into the interactions between lipid bilayers and polyethylene (PETH) chains of varying lengths. The investigation spans diverse lipid bilayer compositions-ranging from pure POPC to POPC:DPPC mixtures-revealing how PETH accommodates itself, adopts extended conformations, and influences membrane structure and ordering. Significantly, while longer PETH chains demonstrate limited passive diffusion, their potential to penetrate bilayers over extended timescales emerges as a significant revelation. Overall, this research significantly advances our comprehension of NP-membrane interactions, shedding light on the potential environmental and health implications that lie ahead.

Membrane specificity of the human cholesterol transfer protein STARD4

STARD4 regulates cholesterol homeostasis by transferring cholesterol between the plasma membrane and endoplasmic reticulum. The STARD4 structure features a helix-grip fold surrounding a large hydrophobic cavity holding the sterol. Its access is controlled by a gate formed by the Ω1 and Ω4 loops and the C-terminal α-helix. Little is known about the mechanisms by which STARD4 binds to membranes and extracts/releases cholesterol. All available structures of STARD4 are without a bound sterol and display the same closed conformation of the gate. The cholesterol transfer activity of the mouse STARD4 is enhanced in the presence of anionic lipids, and in particular of phosphatidylinositol biphosphates (PIP2) for which two binding sites were proposed on the mouse STARD4 surface. Yet only one of these sites is conserved in human STARD4. We here report the results of a liposome microarray-based assay and microseconds-long molecular dynamics simulations of human STARD4 with complex lipid bilayers mimicking the composition of the donor and acceptor membranes. We show that the binding of apo form of human STARD4 is sensitive to the presence of PIP2 through two specific binding sites, one of which was not identified on mouse STARD4. We report two novel conformations of the gate in holo-STARD4: a yet-unobserved close conformation and an open conformation of Ω4 shedding light on the opening/closure mechanism needed for cholesterol uptake/release. Overall, the modulation of human STARD4 membrane-binding by lipid composition, and by the presence of the cargo supports the capacity of human STARD4 to achieve directed transfer between specific organelle membranes.

Small-angle X-ray and neutron scattering applied to lipid-based nanoparticles: Recent advancements across different length scales

Lipid-based nanoparticles (LNPs), ranging from nanovesicles to non-lamellar assemblies, have gained significant attention in recent years, as versatile carriers for delivering drugs, vaccines, and nutrients. Small-angle scattering methods, employing X-rays (SAXS) or neutrons (SANS), represent unique tools to unveil structure, dynamics, and interactions of such particles on different length scales, spanning from the nano to the molecular scale. This review explores the state-of-the-art on scattering methods applied to unveil the structure of lipid-based nanoparticles and their interactions with drugs and bioactive molecules, to inform their rational design and formulation for medical applications. We will focus on complementary information accessible with X-rays or neutrons, ranging from insights on the structure and colloidal processes at a nanoscale level (SAXS) to details on the lipid organization and molecular interactions of LNPs (SANS). In addition, we will review new opportunities offered by Time-resolved (TR)-SAXS and -SANS for the investigation of dynamic processes involving LNPs. These span from real-time monitoring of LNPs structural evolution in response to endogenous or external stimuli (TR-SANS), to the investigation of the kinetics of lipid diffusion and exchange upon interaction with biomolecules (TR-SANS). Finally, we will spotlight novel combinations of SAXS and SANS with complementary on-line techniques, recently enabled at Large Scale Facilities for X-rays and neutrons. This emerging technology enables synchronized multi-method investigation, offering exciting opportunities for the simultaneous characterization of the structure and chemical or mechanical properties of LNPs.

Melittin-Phospholipase A2 Synergism Is Mediated by Liquid-Liquid Miscibility Phase Transition in Giant Unilamellar Vesicles

The primary constituents of honeybee venom, melittin and phospholipase A2 (PLA2), display toxin synergism in which the PLA2 activity is significantly enhanced by the presence of melittin. It has been shown previously that this is accomplished by the disruption in lipid packing, which allows PLA2 to become processive on the membrane surface. In this work, we show that melittin is capable of driving miscibility phase transition in giant unilamellar vesicles (GUVs) and that it raises the miscibility transition temperature (Tmisc) in a concentration-dependent manner. The induced phase separation enhances the processivity of PLA2, particularly at its boundaries, where a substantial difference in domain thickness creates a membrane discontinuity. The catalytic action of PLA2, in response, induces changes in the membrane, rendering it more conducive to melittin binding. This, in turn, facilitates further lipid phase separation and eventual vesicle lysis. Overall, our results show that melittin has powerful membrane-altering capabilities that activate PLA2 in various membrane contexts. More broadly, they exemplify how this biochemical system actively modulates and capitalizes on the spatial distribution of membrane lipids to efficiently achieve its objectives.

A two-step mechanism for the binding of the HIV-1 MPER epitope by the 10E8 antibody onto biosensor-supported lipid bilayers

HIV-1 antibodies targeting the carboxy-terminal area of the membrane-proximal external region (ctMPER) are close to exerting viral pan-neutralization. Here, we reconstituted the ctMPER epitope as the N-terminal extremity of the Env glycoprotein transmembrane domain helix and immobilized it onto biosensor-supported lipid bilayers. We assessed the binding mechanism of anti-MPER antibody 10E8 through Surface Plasmon Resonance, and found, through equilibrium and kinetic binding analyses as a function of bilayer thickness, peptide length, and paratope mutations, that 10E8 engages first with the epitope peptide (encounter), limited by ctMPER helix accessibility at the membrane surface, and then inserts into the lipid bilayer assisted by favorable Fab-membrane interactions (docking). This mechanistic information may help in devising new strategies to develop more efficient MPER-targeting vaccines.

Molecular interactions with bilayer membrane stacks using neutron and X-ray diffraction

Lamellar unit cell reconstruction from neutron and X-ray diffraction data provides information about the disposition and position of molecules and molecular segments with respect to the bilayer. When supplemented with the judicious use of molecular deuteration, the technique probes the molecular interactions and conformations within the bilayer membrane and the water layer which constitute the crystallographic unit cell. The perspective is model independent, and potentially, with a higher degree of resolution than is available with other techniques. In the case of neutron diffraction the measurement consists of carefully normalised diffracted intensity under conditions of contrast variation of the water layer. The subsequent Fourier reconstruction of the unit cell is made using the phase information from variation of peak intensities with contrast. Although the phase problem is not as easily solved for the corresponding X-ray measurements, an intuitive approach can often suffice. Here we discuss the two complimentary techniques as probes of scattering length density profiles of a bilayer, and how such a perspective provides information about the location and orientation of molecules within or between lipid bilayers. Within the basic paradigm of lamellar phases this method has provided, for example, detailed insights into the location and interaction of cryoprotectants and stress proteins, of the mechanisms of actions of viral proteins, antimicrobial compounds and drugs, and the underlying structure of the stratum corneum. In this paper we review these techniques and provide examples of the systems that have been examined. We finish with a future outlook on the use of these techniques to improve our understanding of the interactions of membranes with biomolecules.

Oleic Acid Could Act as a Channel Blocker in the Inhibition of nAChR: Insights from Molecular Dynamics Simulations

The activation of the muscular nicotinic acetylcholine receptor (nAChR) produces the opening of the channel, with the consequent increase in the permeability of cations, triggering an excitatory signal. Free fatty acids (FFA) are known to modulate the activity of the receptor as noncompetitive antagonists, acting at the membrane-AChR interface. We present molecular dynamics simulations of a model of nAChR in a desensitized closed state embedded in a lipid bilayer in which distinct membrane phospholipids were replaced by two different monounsaturated FFA that differ in the position of a double bond. This allowed us to detect and describe that the cis-18:1ω-9 FFA were located at the interface between the transmembrane segments of α2 and γ subunits diffused into the channel lumen with the consequent potential ability to block the channel to the passage of ions.

Comparison of structural effects of cholesterol, lanosterol, and oxysterol on phospholipid (POPC) bilayers

Membrane sterols contribute to the function of biomembranes by regulating the physical properties of the lipid bilayers. Cholesterol, a typical mammalian sterol, is biosynthesized by oxidation of lanosterol. From a molecular evolutionary perspective, lanosterol is considered the ancestral molecule of cholesterol. Here, we studied whether cholesterol is superior to lanosterol in regulating the physical properties of the lipid bilayer in terms of the structural effect on model biomembranes composed of a phospholipid. For comparison, oxysterol, which is formed by oxidation of cholesterol, was also studied. The phospholipid used was 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), which is abundantly found in mammalian biomembranes, and 7β-hydroxycholesterol, which is highly cytotoxic, was used as the oxysterol. The apparent molecular volume was calculated from the mass density determined by the flotation method using H2O and D2O, and the bilayer thickness was determined by reconstructing the electron density distribution from X-ray diffraction data of the POPC/sterol mixtures at a sterol concentration of 30 mol%. The apparent occupied area at the bilayer surface was calculated from the above two structural data. The cholesterol system had the thickest bilayer thickness and the smallest occupied area of the three sterols studied here. This indicates that the POPC/cholesterol bilayer has a better barrier property than the other two systems. Compared to cholesterol, the effects of lanosterol and 7β-hydroxycholesterol on lipid bilayer properties can be interpreted as suboptimal for the function of mammalian biomembranes.

Intrinsic Lipid Curvature and Bilayer Elasticity as Regulators of Channel Function: A Comparative Single-Molecule Study

Perturbations in bilayer material properties (thickness, lipid intrinsic curvature and elastic moduli) modulate the free energy difference between different membrane protein conformations, thereby leading to changes in the conformational preferences of bilayer-spanning proteins. To further explore the relative importance of curvature and elasticity in determining the changes in bilayer properties that underlie the modulation of channel function, we investigated how the micelle-forming amphiphiles Triton X-100, reduced Triton X-100 and the HII lipid phase promoter capsaicin modulate the function of alamethicin and gramicidin channels. Whether the amphiphile-induced changes in intrinsic curvature were negative or positive, amphiphile addition increased gramicidin channel appearance rates and lifetimes and stabilized the higher conductance states in alamethicin channels. When the intrinsic curvature was modulated by altering phospholipid head group interactions, however, maneuvers that promote a negative-going curvature stabilized the higher conductance states in alamethicin channels but destabilized gramicidin channels. Using gramicidin channels of different lengths to probe for changes in bilayer elasticity, we found that amphiphile adsorption increases bilayer elasticity, whereas altering head group interactions does not. We draw the following conclusions: first, confirming previous studies, both alamethicin and gramicidin channels are modulated by changes in lipid bilayer material properties, the changes occurring in parallel yet differing dependent on the property that is being changed; second, isolated, negative-going changes in curvature stabilize the higher current levels in alamethicin channels and destabilize gramicidin channels; third, increases in bilayer elasticity stabilize the higher current levels in alamethicin channels and stabilize gramicidin channels; and fourth, the energetic consequences of changes in elasticity tend to dominate over changes in curvature.

Dimeric Tubulin Modifies Mechanical Properties of Lipid Bilayer, as Probed Using Gramicidin A Channel

Using the gramicidin A channel as a molecular probe, we show that tubulin binding to planar lipid membranes changes the channel kinetics-seen as an increase in the lifetime of the channel dimer-and thus points towards modification of the membrane's mechanical properties. The effect is more pronounced in the presence of non-lamellar lipids in the lipid mixture used for membrane formation. To interpret these findings, we propose that tubulin binding redistributes the lateral pressure of lipid packing along the membrane depth, making it closer to the profile expected for lamellar lipids. This redistribution happens because tubulin perturbs the lipid headgroup spacing to reach the membrane's hydrophobic core via its amphiphilic α-helical domain. Specifically, it increases the forces of repulsion between the lipid headgroups and reduces such forces in the hydrophobic region. We suggest that the effect is reciprocal, meaning that alterations in lipid bilayer mechanics caused by membrane remodeling during cell proliferation in disease and development may also modulate tubulin membrane binding, thus exerting regulatory functions. One of those functions includes the regulation of protein-protein interactions at the membrane surface, as exemplified by VDAC complexation with tubulin.

Antimicrobial peptide interactions with bacterial cell membranes

Antimicrobial peptides (AMPs) are potential alternatives for common antibiotics because of their greater activity and efficiency against a broad range of viruses, bacteria, fungi, and parasites. In this project, two antimicrobial peptides including magainin 2 and protegrin 1 with α-helix and β-sheet secondary structures were selected to investigate their interactions with different lipid bilayers such as 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoserine (POPS), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol (POPG), and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE), POPC/POPG (7:3), POPC/POPS (7:3), POPG/POPE(1:3), and POPG/POPE(3:1). The obtained structures of the AMPs illustrated that protegrin 1 cannot maintain its secondary structure in the solution phase in contrast to magainin 2. The head groups of the lipid units play a key role in the stability of the lipid bilayers. The head parts of the lipid membranes by increasing the internal H-bond contribute to membrane compactness. The POPG and POPS units inside the POPC/POPG and POPC/POPS membranes increase the order of the POPC units. The cationic residues of the AMPs form remarkable electrostatic interactions with the negatively charged membrane surfaces, which play a key role in the stabilization process of the peptide secondary structures. The Arg residues of protegrin 1 and the Gly1, Lys4, Lys10, Lys11, Lys14, and Glu19 of the magainin 2 have the most important roles in the complexation process. The values of Gibbs binding energies (ΔG) indicate that the complexation process between AMPs and different bacterial membranes is favorable from the thermodynamic viewpoint and AMPs could form stable complexes with the lipid bilayers. As a result of ΔG values, protegrin 1 forms a more stable complex with POPG/POPE(3:1), while the α-helix has more affinity to the POPG/POPE(1:3) bacterial membranes. Therefore, it can be considered that β-sheet and α-helix AMPs are more effective against gram-positive and gram-negative bacteria, respectively. The results of this study can provide useful details about the antimicrobial peptide interactions with the bacterial cell, which can be employed for designing new antimicrobial materials with greater efficiency.Communicated by Ramaswamy H. Sarma.

Exploring Cyclic Peptide Nanotube Stability Across Diverse Lipid Bilayers and Unveiling Water Transport Dynamics

Cyclic Peptide Nanotubes (CPNTs) have emerged as compelling candidates for various applications, particularly as nanochannels within lipid bilayers. In this study, the stability of two CPNTs, namely 8 × [(Cys-Gly-Met-Gly)2] and 8 × [(Gly-Leu)4], are comprehensively investigated across different lipid bilayers, including 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), a mixed model membrane (POPE/POPG), and a realistic yeast model membrane. The results demonstrate that both CPNTs maintain their tubular structures in all lipid bilayers, with [(Cys-Gly-Met-Gly)2] showing increased stability over an extended period in these lipid membranes. The insertion of CPNTs shows negligible impact on lipid bilayer properties, including area per lipid, volume per lipid, and bilayer thickness. The study demonstrates that the CPNT preserves its two-line water movement pattern within all the lipid membranes, reaffirming their potential as water channels. The MSD curves further reveal that the dynamics of water molecules inside the nanotube are similar for all the bilayer systems with minor differences that arise due to different lipid environments.

A comparative study of influenza A M2 protein conformations in DOPC/DOPS liposomes and in native E. coli membranes

We compared the conformations of the transmembrane domain (TMD) of influenza A M2 (IAM2) protein reconstituted at pH 7.4 in DOPC/DOPS bilayers to those in isolated E. coli membranes, having preserved its native proteins and lipids. IAM2 is a single-pass transmembrane protein known to assemble into homo-tetrameric proton channel. To represent this channel, we made a construct containing the IAM2's TMD region flanked by the juxtamembrane residues. The single cysteine substitute, L43C, of leucine located in the bilayer polar region was paramagnetically tagged with a methanethiosulfonate nitroxide label for the ESR (electron spin resonance) study. We compared the conformations of the spin-labeled IAM2 residing in DOPC/DOPS and native E. coli membranes using continuous-wave (CW) ESR and double electron-electron resonance (DEER) spectroscopy. The total protein-to-lipid molar ratio spanned the range from 1:230 to 1:10,400⩦ The CW ESR spectra corresponded to a nearly rigid limit spin label dynamics in both environments. In all cases, the DEER data were reconstructed into the distance distributions showing well-resolved peaks at 1.68 nm and 2.37 nm. The peak distance ratio was 1.41±0.2 and the amplitude ratio was 2:1. This is what one expects from four nitroxide spin-labels located at the corners of a square, indicative of an axially symmetric tetramer. Distance modeling of DEER data with molecular modeling software applied to the NMR molecular structures (PDB: 2L0J) confirmed the symmetry and closed state of the C-terminal exit pore of the IAM2 tetramer in agreement with the NMR model. Thus, we can conclude that IAM2 TMD has similar conformations in model and native E. coli membranes of comparable thickness and fluidity, notwithstanding the complexity of the E. coli membranes caused by their lipid diversity and the abundance of integral and peripheral membrane proteins.

MEDUSA: A cloud-based tool for the analysis of X-ray diffuse scattering to obtain the bending modulus from oriented membrane stacks

An important mechanical property of cells is their membrane bending modulus, κ. Here, we introduce MEDUSA (MEmbrane DiffUse Scattering Analysis), a cloud-based analysis tool to determine the bending modulus, κ, from the analysis of X-ray diffuse scattering. MEDUSA uses GPU (graphics processing unit) accelerated hardware and a parallelized algorithm to run the calculations efficiently in a few seconds. MEDUSA's graphical user interface allows the user to upload 2-dimensional data collected from different sources, perform background subtraction and distortion corrections, select regions of interest, run the fitting procedure and output the fitted parameters, the membranes' bending modulus κ, and compressional modulus B.

A Rigorous Framework for Calculating Protein-Protein Binding Affinities in Membranes

Calculating the binding free energy of integral transmembrane (TM) proteins is crucial for understanding the mechanisms by which they recognize one another and reversibly associate. The glycophorin A (GpA) homodimer, composed of two α-helical segments, has long served as a model system for studying TM protein reversible association. The present work establishes a methodological framework for calculating the binding affinity of the GpA homodimer in the heterogeneous environment of a membrane. Our investigation carefully considered a variety of protocols, including the appropriate choice of the force field, rigorous standardization reflecting the experimental conditions, sampling algorithm, anisotropic environment, and collective variables, to accurately describe GpA dimerization via molecular dynamics-based approaches. Specifically, two strategies were explored: (i) an unrestrained potential mean force (PMF) calculation, which merely enhances sampling along the separation of the two binding partners without any restraint, and (ii) a so-called "geometrical route", whereby the α-helices are progressively separated with imposed restraints on their orientational, positional, and conformational degrees of freedom to accelerate convergence. Our simulations reveal that the simplified, unrestrained PMF approach is inadequate for the description of GpA dimerization. Instead, the geometrical route, tailored specifically to GpA in a membrane environment, yields excellent agreement with experimental data within a reasonable computational time. A dimerization free energy of -10.7 kcal/mol is obtained, in fairly good agreement with available experimental data. The geometrical route further helps elucidate how environmental forces drive association before helical interactions stabilize it. Our simulations also brought to light a distinct, long-lived spatial arrangement that potentially serves as an intermediate state during dimer formation. The methodological advances in the generalized geometrical route provide a powerful tool for accurate and efficient binding-affinity calculations of intricate TM protein complexes in inhomogeneous environments.

Cancer and Autism: How PTEN Mutations Degrade Function at the Membrane and Isoform Expression in the Human Brain

Mutations causing loss of PTEN lipid phosphatase activity can promote cancer, benign tumors (PHTS), and neurodevelopmental disorders (NDDs). Exactly how they preferentially trigger distinct phenotypic outcomes has been puzzling. Here, we demonstrate that PTEN mutations differentially allosterically bias P loop dynamics and its connection to the catalytic site, affecting catalytic activity. NDD-related mutations are likely to sample conformations of the functional wild-type state, while sampled conformations for the strong, cancer-related driver mutation hotspots favor catalysis-primed conformations, suggesting that NDD mutations are likely to be weaker, and our large-scale simulations show why. Prenatal PTEN isoform expression data suggest exons 5 and 7, which harbor NDD mutations, as cancer-risk carriers. Since cancer requires more than a single mutation, our conformational and genomic analysis helps discover how same protein mutations can foster different clinical manifestations, articulates a role for co-occurring background latent driver mutations, and uncovers relationships of splicing isoform expression to life expectancy.

The Functionality of Membrane-Inserting Proteins and Peptides: Curvature Sensing, Generation, and Pore Formation

Proteins and peptides with hydrophobic and amphiphilic segments are responsible for many biological functions. The sensing and generation of membrane curvature are the functions of several protein domains or motifs. While some specific membrane proteins play an essential role in controlling the curvature of distinct intracellular membranes, others participate in various cellular processes such as clathrin-mediated endocytosis, where several proteins sort themselves at the neck of the membrane bud. A few membrane-inserting proteins form nanopores that permeate selective ions and water to cross the membrane. In addition, many natural and synthetic small peptides and protein toxins disrupt the membrane by inducing nonspecific pores in the membrane. The pore formation causes cell death through the uncontrolled exchange between interior and exterior cellular contents. In this article, we discuss the insertion depth and orientation of protein/peptide helices, and their role as a sensor and inducer of membrane curvature as well as a pore former in the membrane. We anticipate that this extensive review will assist biophysicists to gain insight into curvature sensing, generation, and pore formation by membrane insertion.

Understanding the Mechanical Properties of Ultradeformable Liposomes Using Molecular Dynamics Simulations

Improving drug delivery efficiency to solid tumor sites is a central challenge in anticancer therapeutic research. Our previous experimental study (Guo et al., Nat. Commun. 2018, 9, 130) showed that soft, elastic liposomes had increased uptake and accumulation in cancer cells and tumors in vitro and in vivo respectively, relative to rigid particles. As a first step toward understanding how liposomes' molecular structure and composition modulates their elasticity, we performed all-atom and coarse-grained classical molecular dynamics (MD) simulations of lipid bilayers formed by mixing a long-tailed unsaturated phospholipid with a short-tailed saturated lipid with the same headgroup. The former types of phospholipids considered were 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) and 1,2-dipalmitoleoyl-sn-glycero-3-phosphocholine (termed here DPMPC). The shorter saturated lipids examined were 1,2-diheptanoyl-sn-glycero-3-phosphocholine (DHPC), 1,2-didecanoyl-sn-glycero-3-phosphocholine (DDPC), 1,2-dilauroyl-sn-glycero-3-phosphocholine (DLPC), and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). Several lipid concentrations and surface tensions were considered. Our results show that DOPC or DPMPC systems having 25-35 mol % of the shortest lipids DHPC or DDPC are the least rigid, having area compressibility moduli KA that are ∼10% smaller than the values observed in pure DOPC or DPMPC bilayers. These results agree with experimental measurements of the stretching modulus and lysis tension in liposomes with the same compositions. These mixed systems also have lower areas per lipid and form more uneven x-y interfaces with water, the tails of both primary and secondary lipids are more disordered, and the terminal methyl groups in the tails of the long lipid DOPC or DPMPC wriggle more in the vertical direction, compared to pure DOPC or DPMPC bilayers or their mixtures with the longer saturated lipid DLPC or DMPC. These observations confirm our hypothesis that adding increasing concentrations of the short unsaturated lipid DHPC or DDPC to DOPC or DPMPC bilayers alters lipid packing and thus makes the resulting liposomes more elastic and less rigid. No formation of lipid nanodomains was noted in our simulations, and no clear trends were observed in the lateral diffusivities of the lipids as the concentration, type of secondary lipid, and surface tension were varied.

Rotation of the c-Ring Promotes the Curvature Sorting of Monomeric ATP Synthases

ATP synthases are proteins that catalyse the formation of ATP through the rotatory movement of their membrane-spanning subunit. In mitochondria, ATP synthases are found to arrange as dimers at the high-curved edges of cristae. Here, a direct link is explored between the rotatory movement of ATP synthases and their preference for curved membranes. An active curvature sorting of ATP synthases in lipid nanotubes pulled from giant vesicles is found. Coarse-grained simulations confirm the curvature-seeking behaviour of rotating ATP synthases, promoting reversible and frequent protein-protein contacts. The formation of transient protein dimers relies on the membrane-mediated attractive interaction of the order of 1.5 kB T produced by a hydrophobic mismatch upon protein rotation. Transient dimers are sustained by a conic-like arrangement characterized by a wedge angle of θ ≈ 50°, producing a dynamic coupling between protein shape and membrane curvature. The results suggest a new role of the rotational movement of ATP synthases for their dynamic self-assembly in biological membranes.

Passive Translocation of Phospholipids in Asymmetric Model Membranes: Solid-State 1H NMR Characterization of Flip-Flop Kinetics Using Deuterated Sphingomyelin and Phosphatidylcholine

Although lateral and inter-leaflet lipid-lipid interactions in cell membranes play roles in maintaining asymmetric lipid bilayers, the molecular basis of these interactions is largely unknown. Here, we established a method to determine the distribution ratio of phospholipids between the outer and inner leaflets of asymmetric large unilamellar vesicles (aLUVs). The trimethylammonium group, (CH3)3N+, in the choline headgroup of N-palmitoyl-sphingomyelin (PSM) and 1,2-dioleoyl-sn-glycero-3-phosphatidylcholine (DOPC) gave rise to a relatively sharp signal in magic-angle spinning solid-state 1H NMR (MAS-ss-1H NMR). PSM and DOPC have the same headgroup structure, but one phospholipid was selectively observed by deuterating the trimethylammonium group of the other phospholipid. The addition of Pr3+ to the medium surrounding aLUVs selectively shifted the chemical shift of the (CH3)3N+ group in the outer leaflet from that in the inner leaflet, which allowed estimation of the inter-leaflet distribution ratio of the unlabeled lipid in aLUVs. Using this method, we evaluated the translocation of PSM and DOPC between the outer and inner leaflets of the cholesterol-containing aLUVs, with PSM and DOPC mostly distributed in the outer and inner leaflets, respectively, immediately after aLUV preparation; their flip and flop rates were approximately 2.7 and 6.4 × 10-6 s-1, respectively. During the passive symmetrization of aLUVs, the lipid translocation rate was decreased due to changes in the membrane order, probably through the formation of the registered liquid-ordered domains. Comparison of the result with that of symmetric LUVs revealed that lipid asymmetry may not significantly affect the lipid translocation rates, while the lateral lipid-lipid interaction may be a dominant factor in lipid translocation under these conditions. These findings highlight the importance of considering the effects of lateral lipid interactions within the same leaflet on lipid flip-flop rates when evaluating the asymmetry of phospholipids in the cell membrane.

Location of Polyelectrolytes in Swollen Lipid Oligobilayers

Osteoarthritis is caused by degeneration of the cartilage, which covers the bone ends of the joints and is decorated with an oligolamellar phospholipid (PL) bilayer. The gap between the bone ends is filled with synovial fluid mainly containing hyaluronic acid (HA). HA and PLs are supposed to reduce friction and protect the cartilage from wear in joint movement. However, a detailed understanding of the molecular mechanisms of joint lubrication is still missing. Previously, we found that aqueous solutions of HA and poly(allylamine hydrochloride) (PAH), the latter serving as a polymeric analogue to HA, adsorb onto the headgroups of surface-bound 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) oligobilayers and significantly enhance their stability with respect to shear forces, typically occurring in joint movement. We now investigated the precise location of PAH chains across the lipid films in neutron reflectivity measurements, as bridging of the oligobilayers by polyelectrolytes (PEs) might be the cause for their improved mechanical stability. In a first set of experiments, we used hydrogenated PAH and chain-deuterated DMPC (DMPC-d54) to improve the contrast between the lipids and potentially intruding PAH. However, due to difficulties in distinguishing between incorporation of water and PAH, penetration into the lipid chain region could hardly be proven quantitatively. Therefore, we designed a more elaborate experiment based on mixed films of DMPC-d54 and hydrogenated DMPC, which is insensitive to water penetration into the films. Beside facilitating a detailed structural characterization of the oligolamellar system, this elaborate approach showed that PAH adsorbs to the DMPC heads and penetrates the lipid tail strata. No PAH was found in the lipid head strata, which excludes bridging of several lipid bilayers by the PE chains. The data are consistent with the assumption that PAH bridges are formed between the headgroups of two adjacent bilayers and contribute to the enhanced mechanical stability.

Permeability of antioxidants through a lipid bilayer model with coarse-grained simulations

Oxidative stress caused by pollution and lifestyle changes causes an excess of free radicals that react chemically with cell constituents leading to irreversible damage. There are molecules known as antioxidants that reduce the levels of free radicals. Some pigments of fruits and vegetables known as anthocyanins have antioxidant properties. Their interaction with the cell membrane becomes a crucial step in studying these substances. In this research, molecular dynamics simulations, particularly, coarse-grained molecular dynamics (CGMD) were used. This technique aims to replace functional groups with corresponding beads that represent their level of polarity and affinities to other chemical groups. Also, umbrella sampling was carried out to obtain free energy profiles that describe well the orientation and location of antioxidants in a membrane considering Trolox, Cyanidin, Gallic Acid, and Resveratrol molecules to study the structural effects they cause on it. It was concluded in this study that an antioxidant when crossing the membrane does not cause either damage to the structural properties or the loss of packing and stratification of phospholipids. it was also observed that the most reactive part of the molecules could easily approach area A prone to lipid oxidation, which can describe the antioxidant capacity of these molecules.Communicated by Ramaswamy H. Sarma.

Cristae formation is a mechanical buckling event controlled by the inner membrane lipidome

Cristae are high curvature structures in the inner mitochondrial membrane (IMM) that are crucial for ATP production. While cristae-shaping proteins have been defined, analogous mechanisms for lipids have yet to be elucidated. Here we combine experimental lipidome dissection with multi-scale modeling to investigate how lipid interactions dictate IMM morphology and ATP generation. When modulating phospholipid (PL) saturation in engineered yeast strains, we observed a surprisingly abrupt breakpoint in IMM topology driven by a continuous loss of ATP synthase organization at cristae ridges. We found that cardiolipin (CL) specifically buffers the IMM against curvature loss, an effect that is independent of ATP synthase dimerization. To explain this interaction, we developed a continuum model for cristae tubule formation that integrates both lipid and protein-mediated curvatures. The model highlighted a snapthrough instability, which drives IMM collapse upon small changes in membrane properties. We also showed that CL is essential in low oxygen conditions that promote PL saturation. These results demonstrate that the mechanical function of CL is dependent on the surrounding lipid and protein components of the IMM.

Surfactant-Mediated Structural Modulations to Planar, Amphiphilic Multilamellar Stacks

The hydrophobic effect, a ubiquitous process in biology, is a primary thermodynamic driver of amphiphilic self-assembly. It leads to the formation of unique morphologies including two highly important classes of lamellar and micellar mesophases. The interactions between these two types of structures and their involved components have garnered significant interest because of their importance in key biochemical technologies related to the isolation, purification, and reconstitution of membrane proteins. This work investigates the structural organization of mixtures of the lamellar-forming phospholipid 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and two zwitterionic micelle-forming surfactants, being n-dodecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (Zwittergent 3-12 or DDAPS) and 1-oleoyl-2-hydroxy-sn-glycero-3-phosphocholine (O-Lyso-PC), when assembled by water vapor hydration with X-ray diffraction measurements, brightfield optical microscopy, wide-field fluorescence microscopy, and atomic force microscopy. The results reveal that multilamellar mesophases of these mixtures can be assembled across a wide range of POPC to surfactant (POPC:surfactant) concentration ratios, including ratios far surpassing the classical detergent-saturation limit of POPC bilayers without significant morphological disruptions to the lamellar motif. The mixed mesophases generally decreased in lamellar spacing (D) and headgroup-to-headgroup distance (Dhh) with a higher concentration of the doped surfactant, but trends in water layer thickness (Dw) between each bilayer in the stack are highly variable. Further structural characteristics including mesophase topography, bilayer thickness, and lamellar rupture force were revealed by atomic force microscopy (AFM), exhibiting homogeneous multilamellar stacks with no significant physical differences with changes in the surfactant concentration within the mesophases. Taken together, the outcomes present the assembly of unanticipated and highly unique mixed mesophases with varied structural trends from the involved surfactant and lipidic components. Modulations in their structural properties can be attributed to the surfactant's chemical specificity in relation to POPC, such as the headgroup hydration and the hydrophobic chain tail mismatch. Taken together, our results illustrate how specific chemical complexities of surfactant-lipid interactions can alter the morphologies of mixed mesophases and thereby alter the kinetic pathways by which surfactants dissolve lipid mesophases in bulk aqueous solutions.

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.

Spontaneous Formation of Ultrasmall Unilamellar Vesicles in Mixtures of an Amphiphilic Drug and a Phospholipid

We have observed ultrasmall unilamellar vesicles, with diameters of less than 20 nm, in mixtures of the tricyclic antidepressant drug amitriptyline hydrochloride (AMT) and the unsaturated zwitterionic phospholipid 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) in physiological saline solution. The size and shape of spontaneously formed self-assembled aggregates have been characterized using complementary techniques, i.e., small-angle neutron and X-ray scattering (SANS and SAXS) and cryo-transmission electron microscopy (cryo-TEM). We observe rodlike mixed micelles in more concentrated samples that grow considerably in length upon dilution, and a transition from micelles to vesicles is observed as the concentration approaches the critical micelle concentration of AMT. Unlike the micelles, the spontaneously formed vesicles decrease in size with each step of dilution, and ultrasmall unilamellar vesicles, with diameters as small as about 15 nm, were observed at the lowest concentrations. The spontaneously formed ultrasmall unilamellar vesicles maintain their size for as long we have investigated them (i.e., several months). To the best of our knowledge, such small vesicles have never before been reported to form spontaneously in a biocompatible phospholipid-based system. Most interestingly, the size of the vesicles was observed to be strongly dependent on the chemical structure of the phospholipid, and in mixtures of AMT and the phospholipid 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC), the vesicles were observed to be considerably larger in size. The self-assembly behavior in the phospholipid-drug surfactant system in many ways resembles the formation of equilibrium micelles and vesicles in mixed anionic/cationic surfactant systems.

The effect of metalation on antimicrobial piscidins imbedded in normal and oxidized lipid bilayers

Metalation of the N-terminal Amino Terminal Cu(ii)- and Ni(ii)-binding (ATCUN) motif may enhance the antimicrobial properties of piscidins. Molecular dynamics simulations of free and nickelated piscidins 1 and 3 (P1 and P3) were performed in 3 : 1 POPC/POPG and 2.6 : 1 : 0.4 POPC/POPG/aldo-PC bilayers (POPC, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine: POPG, 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoglycerol; aldo-PC, 1-palmitoyl-2-(9'-oxo-nonanoyl)-sn-glycero-3-phosphocholine) bilayer models. Nickel(ii) binding decreases the conformation dynamics of the ATCUN motif and lowers the charge of the N-terminus to allow it to embed deeper in the bilayer without significantly changing the overall depth due to interactions of the charged half-helix of the peptide with the headgroups. Phe1⋯Ni2+ cation-π and Phe2-Phe1 CH-π interactions contribute to a small fraction of structures within the nickelated P1 simulations and may partially protect a bound metal from metal-centered chemical activity. The substitution of Phe2 for Ile2 in P3 sterically blocks conformations with cation-π interactions offering less protection to the metal. This difference between metalated P1 and P3 may indicate a mechanism by which peptide sequence can influence antimicrobial properties. Any loss of bilayer integrity due to chain reversal of the oxidized phospholipid chains of aldo-PC may be enhanced in the presence of metalated piscidins.

Label-free quantification of protein binding to lipid vesicles using transparent waveguide evanescent-field scattering microscopy with liquid control

Recent innovations in microscopy techniques are paving the way for label-free studies of single nanoscopic biological entities such as viruses, lipid-nanoparticle drug carriers, and even proteins. One such technique is waveguide evanescent-field microscopy, which offers a relatively simple, yet sensitive, way of achieving label-free light scattering-based imaging of nanoparticles on surfaces. Herein, we extend the application of this technique by incorporating microfluidic liquid control and adapting the design for use with inverted microscopes by fabricating a waveguide on a transparent substrate. We furthermore formulate analytical models describing scattering and fluorescence intensities from single spherical and shell-like objects interacting with evanescent fields. The models are then applied to analyze scattering and fluorescence intensities from adsorbed polystyrene beads and to temporally resolve cholera-toxin B (CTB) binding to individual surface-immobilized glycosphingolipid GM1 containing vesicles. We also propose a self-consistent means to quantify the thickness of the CTB layer, revealing that protein-binding to individual vesicles can be characterized with sub-nm precision in a time-resolved manner.

In Search of a Molecular View of Peptide-Lipid Interactions in Membranes

Lipid bilayer membranes are often represented as a continuous nonpolar slab with a certain thickness bounded by two more polar interfaces. Phenomena such as peptide binding to the membrane surface, folding, insertion, translocation, and diffusion are typically interpreted on the basis of this view. In this Perspective, I argue that this membrane representation as a hydrophobic continuum solvent is not adequate to understand peptide-lipid interactions. Lipids are not small compared to membrane-active peptides: their sizes are similar. Therefore, peptide diffusion needs to be understood in terms of free volume, not classical continuum mechanics; peptide solubility or partitioning in membranes cannot be interpreted in terms of hydrophobic mismatch between membrane thickness and peptide length; peptide folding and translocation, often involving cationic peptides, can only be understood if realizing that lipids adapt to the presence of peptides and the membrane may undergo considerable lipid redistribution in the process. In all of those instances, the detailed molecular interactions between the peptide residues and the lipid components are essential to understand the mechanisms involved.

Unraveling the Transport Properties of RONS across Nitro-Oxidized Membranes

The potential of cold atmospheric plasma (CAP) in biomedical applications has received significant interest, due to its ability to generate reactive oxygen and nitrogen species (RONS). Upon exposure to living cells, CAP triggers alterations in various cellular components, such as the cell membrane. However, the permeation of RONS across nitrated and oxidized membranes remains understudied. To address this gap, we conducted molecular dynamics simulations, to investigate the permeation capabilities of RONS across modified cell membranes. This computational study investigated the translocation processes of less hydrophilic and hydrophilic RONS across the phospholipid bilayer (PLB), with various degrees of oxidation and nitration, and elucidated the impact of RONS on PLB permeability. The simulation results showed that less hydrophilic species, i.e., NO, NO₂, N₂O₄, and O₃, have a higher penetration ability through nitro-oxidized PLB compared to hydrophilic RONS, i.e., HNO₃, s-cis-HONO, s-trans-HONO, H₂O₂, HO₂, and OH. In particular, nitro-oxidation of PLB, induced by, e.g., cold atmospheric plasma, has minimal impact on the penetration of free energy barriers of less hydrophilic species, while it lowers these barriers for hydrophilic RONS, thereby enhancing their translocation across nitro-oxidized PLB. This research contributes to a better understanding of the translocation abilities of RONS in the field of plasma biomedical applications and highlights the need for further analysis of their role in intracellular signaling pathways.

Membrane plasticity induced by myo-inositol derived archaeal lipids: chemical synthesis and biophysical characterization

Archaeal membrane lipids have specific structures that allow Archaea to withstand extreme conditions of temperature and pressure. In order to understand the molecular parameters that govern such resistance, the synthesis of 1,2-di-O-phytanyl-sn-glycero-3-phosphoinositol (DoPhPI), an archaeal lipid derived from myo-inositol, is reported. Benzyl protected myo-inositol was first prepared and then transformed to phosphodiester derivatives using a phosphoramidite based-coupling reaction with archaeol. Aqueous dispersions of DoPhPI alone or mixed with DoPhPC can be extruded and form small unilamellar vesicles, as detected by DLS. Neutron, SAXS, and solid-state NMR demonstrated that the water dispersions could form a lamellar phase at room temperature that then evolves into cubic and hexagonal phases with increasing temperature. Phytanyl chains were also found to impart remarkable and nearly constant dynamics to the bilayer over wide temperature ranges. All these new properties of archaeal lipids are proposed as providers of plasticity and thus means for the archaeal membrane to resist extreme conditions.

Using novel click chemistry algorithm to design D3R inhibitors as blood-brain barrier permeants

Dopamine receptor D3 (D3R) has gained attention as a promising therapeutic target for neurological disorders. In this study, an innovative in silico click reaction strategy was employed to identify potential D3R binders. The ligand template, 1-phenyl-4-[4-(1H-1,2,3-triazol-5-yl)butyl]piperazine, with substitution at the 1,2,3-triazole ring, served as the starting point. Generated compounds underwent filtration based on their brain-to-blood concentration ratio (logBB), leading to the identification of 1-{4-[1-(decahydronaphthalen-1-yl)-1H-1,2,3-triazol-5-yl]butyl}-4-phenylpiperazine as the most promising candidate, displaying superior D3R affinity and blood-brain barrier (BBB) permeability compared to the reference ligand, eticlopride. Molecular dynamics simulations further supported these findings. This study presents a novel hit for designing D3R ligands and establishes a workflow utilizing in silico click chemistry to screen compounds with BBB permeability. The proposed click reaction-based algorithm holds significant potential as a valuable tool in the development of effective antipsychotic compounds.

Extracellular vesicle-matrix interactions

The extracellular matrix in microenvironments harbors a variety of signals to control cellular functions and the materiality of tissues. Most efforts to synthetically reconstitute the matrix by biomaterial design have focused on decoupling cell-secreted and polymer-based cues. Cells package molecules into nanoscale lipid membrane-bound extracellular vesicles and secrete them. Thus, extracellular vesicles inherently interact with the meshwork of the extracellular matrix. In this Review, we discuss various aspects of extracellular vesicle-matrix interactions. Cells receive feedback from the extracellular matrix and leverage intracellular processes to control the biogenesis of extracellular vesicles. Once secreted, various biomolecular and biophysical factors determine whether extracellular vesicles are locally incorporated into the matrix or transported out of the matrix to be taken up by other cells or deposited into tissues at a distal location. These insights can be utilized to develop engineered biomaterials where EV release and retention can be precisely controlled in host tissue to elicit various biological and therapeutic outcomes.

Coarse-grained modeling of polystyrene-modified CNTs and their interactions with lipid bilayers

In the present work, we describe Martini3 coarse-grained models of polystyrene and carboxyl-terminated polystyrene functionalized carbon nanotubes (CNTs) and investigate their interactions with lipid bilayers with and without cholesterol (CHOL) using molecular dynamics simulations. By changing the polystyrene chain length and grafting density at the end ring of the CNTs at two different nanotube concentrations, we observe the translocation of nanoparticles as well as changes in the lipid bilayer properties. Our results show that all developed models passively diffuse into the membranes without causing any damage to the membrane integrity, although high concentrations of CNTs induce structural and elastic changes in lipid bilayers. In the presence of CHOL, increasing CNT concentration results in decreased rates of CHOL transmembrane motions. On the other hand, CNTs are prone to lipid and polystyrene blockage, which affects their equilibrated configurations, and tilting behavior within the membranes. Hence, we demonstrate that polystyrene-functionalized CNTs are promising drug-carrier agents. However, polystyrene chain length and grafting density are important factors to consider to enhance the efficiency of drug delivery.

Structural factors governing binding of curvature-sensing peptides to bacterial extracellular vesicles covered with hydrophilic polysaccharide chains

Extracellular vesicles (EVs) have attracted an attention as important targets in the fields of biology and medical science because they contain physiologically active molecules. Curvature-sensing peptides are currently used as novel tools for marker-independent EV detection techniques. A structure-activity correlation study demonstrated that the α-helicity of the peptides is prominently involved in peptide binding to vesicles. However, whether a flexible structure changing from a random coil to an α-helix upon binding to vesicles or a restricted α-helical structure is an important factor in the detection of biogenic vesicles is still unclear. To address this issue, we compared the binding affinities of stapled and unstapled peptides for bacterial EVs with different surface polysaccharide chains. We found that unstapled peptides showed similar binding affinities for bacterial EVs regardless of surface polysaccharide chains, whereas stapled peptides showed substantially decreased binding affinities for bacterial EVs covered with capsular polysaccharides. This is probably because curvature-sensing peptides must pass through the layer of hydrophilic polysaccharide chains prior to binding to the hydrophobic membrane surface. While stapled peptides with restricted structures cannot easily pass through the layer of polysaccharide chains, unstapled peptides with flexible structures can easily approach the membrane surface. Therefore, we concluded that the structural flexibility of curvature-sensing peptides is a key factor for governing the highly sensitive detection of bacterial EVs.

Drude Polarizable Lipid Force Field with Explicit Treatment of Long-Range Dispersion: Parametrization and Validation for Saturated and Monounsaturated Zwitterionic Lipids

Accurate empirical force fields of lipid molecules are a critical component of molecular dynamics simulation studies aimed at investigating properties of monolayers, bilayers, micelles, vesicles, and liposomes, as well as heterogeneous systems, such as protein-membrane complexes, bacterial cell walls, and more. While the majority of lipid force field-based simulations have been performed using pairwise-additive nonpolarizable models, advances have been made in the development of the polarizable force field based on the classical Drude oscillator model. In the present study, we undertake further optimization of the Drude lipid force field, termed Drude2023, including improved treatment of the phosphate and glycerol linker region of PC and PE headgroups, additional optimization of the alkene group in monounsaturated lipids, and inclusion of long-range Lennard-Jones interactions using the particle-mesh Ewald method. Initial optimization targeted quantum mechanical (QM) data on small model compounds representative of the linker region. Subsequent optimization targeted QM data on larger model compounds, experimental data, and dihedral potentials of mean force from the CHARMM36 additive lipid force field using a parameter reweighting protocol. The use of both experimental and QM target data during the reweighting protocol is shown to produce physically reasonable parameters that reproduce a collection of experimental observables. Target data for optimization included surface area/lipid for DPPC, DSPC, DMPC, and DLPC bilayers and nuclear magnetic resonance (NMR) order parameters for DPPC bilayers. Validation data include prediction of membrane thickness, scattering form factors, electrostatic potential profiles, compressibility moduli, surface area per lipid, water permeability, NMR T1 relaxation times, diffusion constants, and monolayer surface tensions for a variety of saturated and unsaturated lipid mono- and bilayers. Overall, the agreement with experimental data is quite good, though the results are less satisfactory for the NMR T1 relaxation times for carbons near the ester groups. Notable improvements compared to the additive C36 force field were obtained for membrane dipole potentials, lipid diffusion coefficients, and water permeability with the exception of monounsaturated lipid bilayers. It is anticipated that the optimized polarizable Drude2023 force field will help generate more accurate molecular simulations of pure bilayers and heterogeneous systems containing membranes, advancing our understanding of the role of electronic polarization in these systems.

HylleraasMD: A Domain Decomposition-Based Hybrid Particle-Field Software for Multiscale Simulations of Soft Matter

We present HylleraasMD (HyMD), a comprehensive implementation of the recently proposed Hamiltonian formulation of hybrid particle-field molecular dynamics. The methodology is based on a tunable, grid-independent length-scale of coarse graining, obtained by filtering particle densities in reciprocal space. This enables systematic convergence of energies and forces by grid refinement, also eliminating nonphysical force aliasing. Separating the time integration of fast modes associated with internal molecular motion from slow modes associated with their density fields, we enable the first time-reversible, energy-conserving hybrid particle-field simulations. HyMD comprises the optional use of explicit electrostatics, which, in this formalism, corresponds to the long-range potential in particle-mesh Ewald. We demonstrate the ability of HyMD to perform simulations in the microcanonical and canonical ensembles with a series of test cases, comprising lipid bilayers and vesicles, surfactant micelles, and polypeptide chains, comparing our results to established literature. An on-the-fly increase of the characteristic coarse-grain length significantly speeds up dynamics, accelerating self-diffusion and leading to expedited aggregation. Exploiting this acceleration, we find that the time scales involved in the self-assembly of polymeric structures can lie in the tens to hundreds of picoseconds instead of the multimicrosecond regime observed with comparable coarse-grained models.

Novel Helical Trp- and Arg-Rich Antimicrobial Peptides Locate Near Membrane Surfaces and Rigidify Lipid Model Membranes

Antibiotics are losing effectiveness as bacteria become resistant to conventional drugs. To find new alternatives, antimicrobial peptides (AMPs) are rationally designed with different lengths, charges, hydrophobicities (H), and hydrophobic moments (μH), containing only three types of amino acids: arginine, tryptophan, and valine. Six AMPs with low minimum inhibitory concentrations (MICs) and <25% toxicity to mammalian cells are selected for biophysical studies. Their secondary structures are determined using circular dichroism (CD), which finds that the % α-helicity of AMPs depends on composition of the lipid model membranes (LMMs): gram-negative (G(-)) inner membrane (IM) >gram-positive (G(+)) > Euk33 (eukaryotic with 33 mol% cholesterol). The two most effective peptides, E2-35 (16 amino acid [AA] residues) and E2-05 (22 AAs), are predominantly helical in G(-) IM and G(+) LMMs. AMP/membrane interactions such as membrane elasticity, chain order parameter, and location of the peptides in the membrane are investigated by low-angle and wide-angle X-ray diffuse scattering (XDS). It is found that headgroup location correlates with efficacy and toxicity. The membrane bending modulus KC displays nonmonotonic changes due to increasing concentrations of E2-35 and E2-05 in G(-) and G(+) LMMs, suggesting a bacterial killing mechanism where domain formation causes ion and water leakage.

Modeling of flexible membrane-bound biomolecular complexes for solution small-angle scattering

Recent advances in protein expression protocols, sample handling, and experimental set up of small-angle scattering experiments have allowed users of the technique to structurally investigate biomolecules of growing complexity and structural disorder. Notable examples include intrinsically disordered proteins, multi-domain proteins and membrane proteins in suitable carrier systems. Here, we outline a modeling scheme for calculating the scattering profiles from such complex samples. This kind of modeling is necessary for structural information to be refined from the corresponding data. The scheme bases itself on a hybrid of classical form factor based modeling and the well-known spherical harmonics-based formulation of small-angle scattering amplitudes. Our framework can account for flexible domains alongside other structurally elaborate components of the molecular system in question. We demonstrate the utility of this modeling scheme through a recent example of a structural model of the growth hormone receptor membrane protein in a phospholipid bilayer nanodisc which is refined against experimental SAXS data. Additionally we investigate how the scattering profiles from the complex would appear under different scattering contrasts. For each contrast situation we discuss what structural information is contained and the related consequences for modeling of the data.

Unraveling the Aquaporin-3 Inhibitory Effect of Rottlerin by Experimental and Computational Approaches

The natural polyphenolic compound Rottlerin (RoT) showed anticancer properties in a variety of human cancers through the inhibition of several target molecules implicated in tumorigenesis, revealing its potential as an anticancer agent. Aquaporins (AQPs) are found overexpressed in different types of cancers and have recently emerged as promising pharmacological targets. Increasing evidence suggests that the water/glycerol channel aquaporin-3 (AQP3) plays a key role in cancer and metastasis. Here, we report the ability of RoT to inhibit human AQP3 activity with an IC₅₀ in the micromolar range (22.8 ± 5.82 µM for water and 6.7 ± 2.97 µM for glycerol permeability inhibition). Moreover, we have used molecular docking and molecular dynamics simulations to understand the structural determinants of RoT that explain its ability to inhibit AQP3. Our results show that RoT blocks AQP3-glycerol permeation by establishing strong and stable interactions at the extracellular region of AQP3 pores interacting with residues essential for glycerol permeation. Altogether, our multidisciplinary approach unveiled RoT as an anticancer drug against tumors where AQP3 is highly expressed providing new information to aquaporin research that may boost future drug design.

Location of dopamine in lipid bilayers and its relevance to neuromodulator function

Dopamine (DA) is a neurotransmitter that also acts as a neuromodulator, with both functions being essential to brain function. Here, we present the first experimental measurement of DA location in lipid bilayers using x-ray diffuse scattering, solid-state deuterium NMR, and electron paramagnetic resonance. We find that the association of DA with lipid headgroups as seen in electron density profiles leads to an increase of intermembrane repulsion most likely due to electrostatic charging. DA location in the lipid headgroup region also leads to an increase of the cross-sectional area per lipid without affecting the bending rigidity significantly. The order parameters measured by solid-state deuterium NMR decrease in the presence of DA for the acyl chains of PC and PS lipids, consistent with an increase in the area per lipid due to DA. Most importantly, these results support the hypothesis that three-dimensional diffusion of DA to target membranes could be followed by relatively more efficient two-dimensional diffusion to receptors within those membranes.

Structure of the gel phase of diC22:1PC lipid bilayers determined by x-ray diffraction

High-resolution x-ray data are reported for the ordered phases of long-chain di-monounsaturated C22:1 phosphocholine lipid bilayers. Similar to PC lipids that have saturated chains, diC22:1PC has a subgel phase and a gel phase, but dissimilarly, we find no ripple phase. Our quantitative focus is on the structure of the gel phase. We have recorded 17 lamellar orders, indicating a very well-ordered structure. Fitting to a model provides the phases of the orders. The Fourier construction of the electron density profile has two well-defined headgroup peaks and a very sharp and deep methyl trough. The wide-angle scattering exhibits two Bragg rods that provide the area per molecule. They have an intensity pattern quite different than that of lipids with saturated chains. Models of chain packing indicate that ground state chain configurations are tilted primarily toward next nearest neighbors with an angle that is also consistent with the modeling of the electron density profile. Wide-angle modeling also indicates broken mirror symmetry between the monolayers. Our wide-angle results and our electron density profile together leads to the hypothesis that the sn-1 and sn-2 chains have equivalent penetration depths in contrast to the gel phase structure of lipids with saturated hydrocarbon chains.

Exhaustive in silico design and screening of novel antipsychotic compounds with improved pharmacodynamics and blood-brain barrier permeation properties

Antipsychotic drugs or neuroleptics are widely used in the treatment of psychosis as a manifestation of schizophrenia and bipolar disorder. However, their effectiveness largely depends on the blood-brain barrier (BBB) permeation (pharmacokinetics) and drug-receptor pharmacodynamics. Therefore, in this study, we developed and implemented the in silico pipeline to design novel compounds (n = 260) as leads using the standard drug scaffolds with improved PK/PD properties from the standard scaffolds. As a result, the best candidates (n = 3) were evaluated in molecular docking to interact with serotonin and dopamine receptors. Finally, haloperidol (HAL) derivative (1-(4-fluorophenyl)-4-(4-hydroxy-4-{4-[(2-phenyl-1,3-thiazol-4-yl)methyl]phenyl}piperidin-1-yl)butan-1-one) was identified as a "magic shotgun" lead compound with better affinity to the 5-HT2A, 5-HT1D, D2, D3, and 5-HT1B receptors than the control molecule. Additionally, this hit substance was predicted to possess similar BBB permeation properties and much lower toxicological profiles in comparison to HAL. Overall, the proposed rational drug design platform for novel antipsychotic drugs based on the BBB permeation and receptor binding might be an invaluable asset for a medicinal chemist or translational pharmacologist.Communicated by Ramaswamy H. Sarma.

Effect of Cholesterol on the Structure and Composition of Glyco-DIBMA Lipid Particles

Different amphiphilic co-polymers have been introduced to produce polymer-lipid particles with nanodisc structure composed of an inner lipid bilayer and polymer chains self-assembled as an outer belt. These particles can be used to stabilize membrane proteins in solution and enable their characterization by means of biophysical methods, including small-angle X-ray scattering (SAXS). Some of these co-polymers have also been used to directly extract membrane proteins together with their associated lipids from native membranes. Styrene/maleic acid and diisobutylene/maleic acid are among the most commonly used co-polymers for producing polymer-lipid particles, named SMALPs and DIBMALPs, respectively. Recently, a new co-polymer, named Glyco-DIBMA, was produced by partial amidation of DIBMA with the amino sugar N-methyl-d-glucosamine. Polymer-lipid particles produced with Glyco-DIBMA, named Glyco-DIBMALPs, exhibit improved structural properties and stability compared to those of SMALPs and DIBMALPs while retaining the capability of directly extracting membrane proteins from native membranes. Here, we characterize the structure and lipid composition of Glyco-DIBMALPs produced with either 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) or 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC). Glyco-DIBMALPs were also prepared with mixtures of either POPC or DMPC and cholesterol at different mole fractions. We estimated the lipid content in the Glyco-DIBMALPs and determined the particle structure and morphology by SAXS. We show that the Glyco-DIBMALPs are nanodisc-like particles whose size and shape depend on the polymer/lipid ratio. This is relevant for designing nanodisc particles with a tunable diameter according to the size of the membrane protein to be incorporated. We also report that the addition of >20 mol % cholesterol strongly perturbed the formation of Glyco-DIBMALPs. Altogether, we describe a detailed characterization of the Glyco-DIBMALPs, which provides relevant inputs for future application of these particles in the biophysical investigation of membrane proteins.

Confinement effect on hydrolysis in small lipid vesicles

In living organisms most chemical reactions take place within the confines of lipid-membrane bound compartments, while confinement within the bounds of a lipid membrane is thought to be a key step in abiogenesis. In previous work we demonstrated that confinement in the aqueous cavity of a lipid vesicle affords protection against hydrolysis, a phenomenon that we term here confinement effect (C _e) and that we attributed to the interaction with the lipid membrane. Here, we show that both the size and the shape of the cavity of the vesicle modulate the C _e. We link this observation to the packing of the lipid following changes in membrane curvature, and formulate a mathematical model that relates the C _e to the radius of a spherical vesicle and the packing parameter of the lipids. These results suggest that the shape of the compartment where a molecule is located plays a major role in controlling the chemical reactivity of non-enzymatic reactions. Moreover, the mathematical treatment we propose offers a useful tool for the design of vesicles with predictable reaction rates of the confined molecules, e.g., drug delivery vesicles with confined prodrugs. The results also show that a crude form of signal transduction, devoid of complex biological machinery, can be achieved by any external stimuli that drastically changes the structure of the membrane, like the osmotic shocks used in the present work.

A skin pharmacokinetics study of permeation enhancers: The root cause of dynamic enhancement effect on in vivo drug permeation

Skin pharmacokinetics (SPK) of permeation enhancers can answer the question of why enhancement effects different at the kinetic level. Herein, SPK of permeation enhancers were classified into two categories, namely, lateral elimination (elimination to surrounding stratum corneum (SC)) and longitudinal elimination (elimination to deep epidermal (EP)). They were evaluated with a specific parameter for permeation enhancers, diffusion ratio (DR_SC-EP), according to results of tissue-distribution test, molecular dynamic (MD) simulation, and confocal laser scanning microscopy (CLSM). The linear relationship between k_e-enahcer and Δ C_max-drug (R² = 0.92), MRT_enhancer and Δ T_max-drug (R² = 0.97), AUC_t-enhancer and Δ AUC_t-drug (R² = 0.90) suggesting that SPK of permeation enhancers precisely controlled dynamic process of drug permeation in vivo. The molecular mechanisms of the dynamic effect of SPK process on drug transdermal behaviors were characterized by modulated-temperature differential scanning calorimetry (MTDSC), dielectric spectroscopy, small-angle X-ray scattering (SAXS), solid-state NMR. Permeation enhancers with high molecular weight (M.W.) and high polar surface area (P.S.A.) had good compatibility and strong interaction strength with SC, leading their lateral-elimination behavior, causing their low DR_SC-EP and resulting in low k_e-enhancer, long MRT_enhancer, and large AUC_t-enhancer. Consequently, skin barrier can be rapidly opened fast and to a great extent. In summary, compared with SPK of permeation enhancers with longitudinal elimination, SPK of permeation enhancers with lateral elimination can enable more sustainable and greater drug permeation. The information about SPK of permeation enhancers offered a criterion to estimate its permeation-enhancement effect on the drug and its subsequent application in transdermal formulations.

Incorporation and localisation of alkanes in a protomembrane model by neutron diffraction

Protomembranes at the origin of life were likely composed of short-chain lipids, readily available on the early Earth. Membranes formed by such lipids are less stable and more permeable under extreme conditions, so a novel membrane architecture was suggested to validate the accuracy of this assumption. The model membrane includes the presence of a layer of alkanes in the mid-plane of the protomembrane in between the two monolayer leaflets and lying perpendicular to the lipid acyl chains. Here, we investigated such a possibility experimentally for membranes formed by the short-chain phospholipid 1,2-didecanoyl-sn-glycero-3-phophocholine, including or not the alkanes eicosane, squalane or triacontane by means of neutron membrane diffraction and contrast variation. We found strong indications for incorporation of two of the three alkanes in the membrane mid-plane through the determination of neutron scattering length density profiles with hydrogenated vs deuterated alkanes and membrane swelling at various relative humidities indicating a slightly increased bilayer thickness when the alkanes are incorporated into the bilayers. The selectivity of the incorporation points out the role of the length of the n-alkanes with respect to the capacity of the membrane to incorporate them.

Coarse-Grained Molecular Dynamics of pH-Sensitive Lipids

pH-sensitive lipids represent a class of lipids that can be protonated and destabilized in acidic environments, as they become positively charged in response to low-pH conditions. They can be incorporated into lipidic nanoparticles such as liposomes, which are able to change their properties and allow specific drug delivery at the acidic conditions encountered in some pathological microenvironments. In this work, we used coarse-grained molecular-dynamic simulations to study the stability of neutral and charged lipid bilayers containing POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and various kinds of ISUCA ((F)2-(imidazol-1-yl)succinic acid)-derived lipids, which can act as pH-sensitive molecules. In order to explore such systems, we used a MARTINI-derived forcefield, previously parameterized using all-atom simulation results. We calculated the average area per lipid, the second-rank order parameter and the lipid diffusion coefficient of both lipid bilayers made of pure components and mixtures of lipids in different proportions, under neutral or acidic conditions. The results show that the use of ISUCA-derived lipids disturbs the lipid bilayer structure, with the effect being particularly marked under acidic conditions. Although more-in depth studies on these systems must be carried out, these initial results are encouraging and the lipids designed in this research could be a good basis for developing new pH-sensitive liposomes.

The Effect of Cholesterol in SOPC Lipid Bilayers at Low Temperatures

We study the behavior of lipid bilayers composed of SOPC (1-stearoyl-2-oleoyl-sn-glycero-3-phosphocholine) with different concentrations of cholesterol, ranging from 10 mol% to 50 mol% at 273 K. To this end, we carry out extensive atomistic molecular dynamic simulations with the aid of the Slipid force field aiming at computing basic bilayer parameters, as well as thermodynamic properties and structural characteristics. The obtained results are compared to available relevant experimental data and the outcome of atomistic simulations performed on bilayers composed of analogous phospholipids. Our results show a good quantitative, as well as qualitative, agreement with the main trends associated with the concentration increase in cholesterol. Moreover, it comes out that a change in the behavior of the bilayer is brought about at a concentration of about 30 mol% cholesterol. At this very concentration, some of the bilayer properties are found to exhibit a saturation and a significant long-range ordering of the lipid molecules in the membrane shows up.

How Does the Study MD of pH-Dependent Exposure of Nanoparticles Affect Cellular Uptake of Anticancer Drugs?

The lack of knowledge about the uptake of NPs by biological cells poses a significant problem for drug delivery. For this reason, designing an appropriate model is the main challenge for modelers. To address this problem, molecular modeling studies that can describe the mechanism of cellular uptake of drug-loaded nanoparticles have been conducted in recent decades. In this context, we developed three different models for the amphipathic nature of drug-loaded nanoparticles (MTX-SS-γ-PGA), whose cellular uptake mechanism was predicted by molecular dynamics studies. Many factors affect nanoparticle uptake, including nanoparticle physicochemical properties, protein-particle interactions, and subsequent agglomeration, diffusion, and sedimentation. Therefore, the scientific community needs to understand how these factors can be controlled and the NP uptake of nanoparticles. Based on these considerations, in this study, we investigated for the first time the effects of the selected physicochemical properties of the anticancer drug methotrexate (MTX) grafted with hydrophilic-γ-polyglutamic acid (MTX-SS-γ-PGA) on its cellular uptake at different pH values. To answer this question, we developed three theoretical models describing drug-loaded nanoparticles (MTX-SS-γ-PGA) at three different pH values, such as (1) pH 7.0 (the so-called neutral pH model), (2) pH 6.4 (the so-called tumor pH model), and (3) pH 2.0 (the so-called stomach pH model). Exceptionally, the electron density profile shows that the tumor model interacts more strongly with the head groups of the lipid bilayer than the other models due to charge fluctuations. Hydrogen bonding and RDF analyses provide information about the solution of the NPs with water and their interaction with the lipid bilayer. Finally, dipole moment and HOMO-LUMO analysis showed the free energy of the solution in the water phase and chemical reactivity, which are particularly useful for determining the cellular uptake of the NPs. The proposed study provides fundamental insights into molecular dynamics (MD) that will allow researchers to determine the influence of pH, structure, charge, and energetics of NPs on the cellular uptake of anticancer drugs. We believe that our current study will be useful in developing a new model for drug delivery to cancer cells with a much more efficient and less time-consuming model.