On the Calculation of Acyl Chain Order Parameters from Lipid Simulations

Simulating Bacterial Membrane Models at the Atomistic Level: A Force Field Comparison

Molecular dynamics (MD) simulations are currently an indispensable tool to understand both the dynamic and nanoscale organization of cell membrane models. A large number of quantitative parameters can be extracted from these simulations, but their reliability is determined by the quality of the employed force field and the simulation parameters. Much of the work on parametrizing and optimizing force fields for biomembrane modeling has been focused on homogeneous bilayers with a single phospholipid type. However, these may not perform effectively or could even be unsuitable for lipid mixtures commonly employed in membrane models. This work aims to fill this gap by comparing MD simulation results of several bacterial membrane models using different force fields and simulation parameters, namely, CHARMM36, Slipids, and GROMOS-CKP. Furthermore, the hydrogen isotope exchange (HIE) method, combined with GROMOS-CKP (GROMOS-H2Q), was also tested to check for the impact of this acceleration strategy on the performance of the force field. A common set of simulation parameters was employed for all of the force fields in addition to those corresponding to the original parametrization of each of them. Furthermore, new experimental order parameter values determined from NMR of several lipid mixtures are also reported to compare them with those determined from MD simulations. Our results reveal that most of the calculated physical properties of bacterial membrane models from MD simulations are substantially force field and lipid composition dependent. Some lipid mixtures exhibit nearly ideal behaviors, while the interaction of different lipid types in other mixtures is highly synergistic. None of the employed force fields seem to be clearly superior to the other three, each having its own strengths and weaknesses. Slipids are notably effective at replicating the order parameters for all acyl chains, including those in lipid mixtures, but they offer the least accurate results for headgroup parameters. Conversely, CHARMM provides almost perfect estimates for the order parameters of the headgroups but tends to overestimate those of the lipid tails. The GROMOS parametrizations deliver reasonable order parameters for entire lipid molecules, including multicomponent bilayers, although they do not reach the accuracy of Slipids for tails or CHARMM for headgroups. Importantly, GROMOS-H2Q stands out for its computational efficiency, being at least 3 times faster than GROMOS, which is already faster than both CHARMM and Slipids. In turn, GROMOS-H2Q yields much higher compressibilities compared to all other parametrizations.

Free Energy Analysis of Peptide-Induced Pore Formation in Lipid Membranes by Bridging Atomistic and Coarse-Grained Simulations

Antimicrobial peptides (AMPs) are attractive materials for combating the antimicrobial resistance crisis because they can kill target microbes by directly disrupting cell membranes. Although thousands of AMPs have been discovered, their molecular mechanisms of action are still poorly understood. One broad mechanism for membrane disruption is the formation of membrane-spanning hydrophilic pores which can be stabilized by AMPs. In this study, we use molecular dynamics simulations to investigate the thermodynamics of pore formation in model single-component lipid membranes in the presence of one of three AMPs: aurein 1.2, melittin and magainin 2. To overcome the general challenge of modeling long time scale membrane-related behaviors, including AMP binding, clustering, and pore formation, we develop a generalizable methodology for sampling AMP-induced pore formation. This approach involves the long equilibration of peptides around a pore created with a nucleation collective variable by performing coarse-grained simulations, then backmapping equilibrated AMP-membrane configurations to all-atom resolution. We then perform all-atom simulations to resolve free energy profiles for pore formation while accurately modeling the interplay of lipid-peptide-solvent interactions that dictate pore formation free energies. Using this approach, we quantify free energy barriers for pore formation without direct biases on peptides or whole lipids, allowing us to investigate mechanisms of pore formation for these 3 AMPs that are a consequence of unbiased peptide diffusion and clustering. Further analysis of simulation trajectories then relates variations in pore lining by AMPs, AMP-induced lipid disruptions, and salt bridges between AMPs to the observed pore formation free energies and corresponding mechanisms. This methodology and mechanistic analysis have the potential to generalize beyond the AMPs in this study to improve our understanding of pore formation by AMPs and related antimicrobial materials.

Pore Formation by Amyloid-like Peptides: Effects of the Nonpolar-Polar Sequence Pattern

One of the mechanisms accounting for the toxicity of amyloid peptides in diseases like Alzheimer's and Parkinson's is the formation of pores on the plasma membrane of neurons. Here, we perform unbiased all-atom simulations of the full membrane damaging pathway, which includes adsorption, aggregation, and perforation of the lipid bilayer accounting for pore-like structures. Simulations are performed using four peptides made with the same amino acids. Differences in the nonpolar-polar sequence pattern of these peptides prompt them to adsorb into the membrane with the extended conformations oriented either parallel [peptide labeled F1, Ac-(FKFE)2-NH2], perpendicular (F4, Ac-FFFFKKEE-NH2), or with an intermediate orientation (F2, Ac-FFKKFFEE-NH2, and F3, Ac-FFFKFEKE-NH2) in regard to the membrane surface. At the water-lipid interface, only F1 fully self-assembles into β-sheets, and F2 peptides partially fold into an α-helical structure. The β-sheets of F1 emerge as electrostatic interactions attract neighboring peptides to intermediate distances where nonpolar side chains can interact within the dry core of the bilayer. This complex interplay between electrostatic and nonpolar interactions is not observed for the other peptides. Although β-sheets of F1 peptides are mostly parallel to the membrane, some of their edges penetrate deep inside the bilayer, dragging water molecules with them. This precedes pore formation, which starts with the flow of two water layers through the membrane that expand into a stable cylindrical pore delimited by polar faces of β-sheets spanning both leaflets of the bilayer.

The phase behavior of skin-barrier lipids: A combined approach of experiments and simulations

Skin barrier function is localized in its outermost layer, the stratum corneum (SC), which is comprised of corneocyte cells embedded in an extracellular lipid matrix containing ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs). The unique structure and composition of this lipid matrix are important for skin barrier function. In this study, experiments and molecular dynamics simulation were combined to investigate the structural properties and phase behavior of mixtures containing nonhydroxy sphingosine CER (CER NS), CHOL, and FFA. X-ray scattering for mixtures with varying CHOL levels revealed the presence of the 5.4 nm short periodicity phase in the presence of CHOL. Bilayers in coarse-grained multilayer simulations of the same compositions contained domains with thicknesses of approximately 5.3 and 5.8 nm that are associated with elevated levels, respectively, of CER sphingosine chains with CHOL, and CER acyl chains with FFA chains. The prevalence of the thicker domain increased with decreasing CHOL content. This might correspond to a phase with ∼5.8 nm spacing observed by x-rays (other details unknown) in mixtures with lower CHOL content. Scissoring and stretching frequencies from Fourier transform infrared spectroscopy (FTIR) also indicate interaction between FFA and CER acyl chains and little interaction between CER acyl and CER sphingosine chains, which requires CER molecules to adopt a predominantly extended conformation. In the simulated systems, neighbor preferences of extended CER chains align more closely with the FTIR observations than those of CERs with hairpin ceramide chains. Both FTIR and atomistic simulations of reverse mapped multilayer membranes detect a hexagonal to fluid phase transition between 65 and 80°C. These results demonstrate the utility of a collaborative experimental and simulation effort in gaining a more comprehensive understanding of SC lipid membranes.

Lateral diffusion of ions near membrane surface

Biological membranes isolate living cells from their environment, while allowing selective molecular transport between the inner and outer realms. For example, Na+ and K+ permeability through ionic channels contributes to neural conduction. Whether the ionic currents arise directly from cations in the bulk, or from the interface, is currently unclear. There are only scant results concerning lateral diffusion of ions on aquated membrane surfaces (and strong belief that this occurs through binding to a diffusing lipid). We performed classical molecular dynamics (MD) simulations of monovalent ions, Na+, K+, and Cl-, near the surface of the zwitterionic palmitoyl-oleoyl-phosphatidylcholine (POPC) membrane. Realistic force-fields for lipids (Amber's Lipid17 and Lipid21) and water (TIP4P-Ew) are tested for the mass and charge densities and the electrostatic potential across the membrane. These calculations reveal that the chloride can bind to the choline moiety through an intervening water molecule by forming a CH⋯OH hydrogen bond, while cations bind to both the phosphatic and carbonyl oxygens of phosphatidylcholine moieties. Upon transitioning from the bulk to the interface, a cation sheds some of its hydration water, which are replaced by headgroup atoms. Notably, an interfacial cation can bind 1-4 headgroup atoms, which is a key to understanding its surface hopping mechanism. We find that cation binding to three headgroup atoms immobilizes it, while binding to four energizes it. Consequently, the lateral cation diffusion rate is only 15-25 times slower than in the bulk, and 4-5 times faster than lipid self-diffusion. K+ diffusion is notably more anomalous than Na+, switching from sub- to super-diffusion after about 2 ns.

Topological Learning Approach to Characterizing Biological Membranes

Biological membranes play key roles in cellular compartmentalization, structure, and its signaling pathways. At varying temperatures, individual membrane lipids sample from different configurations, a process that frequently leads to higher-order phase behavior and phenomena. Here, we present a persistent homology (PH)-based method for quantifying the structural features of individual and bulk lipids, providing local and contextual information on lipid tail organization. Our method leverages the mathematical machinery of algebraic topology and machine learning to infer temperature-dependent structural information on lipids from static coordinates. To train our model, we generated multiple molecular dynamics trajectories of dipalmitoyl-phosphatidylcholine membranes at varying temperatures. A fingerprint was then constructed for each set of lipid coordinates by PH filtration, in which interaction spheres were grown around the lipid atoms while tracking their intersections. The sphere filtration formed a simplicial complex that captures enduring key topological features of the configuration landscape using homology, yielding persistence data. Following fingerprint extraction for physiologically relevant temperatures, the persistence data were used to train an attention-based neural network for assignment of effective temperature values to selected membrane regions. Our persistence homology-based method captures the local structural effects, via effective temperature, of lipids adjacent to other membrane constituents, e.g., sterols and proteins. This topological learning approach can predict lipid effective temperatures from static coordinates across multiple spatial resolutions. The tool, called MembTDA, can be accessed at https://github.com/hyunp2/Memb-TDA.

Interaction of MRI Contrast Agent [Gd(DOTA)]- with Lipid Membranes: A Molecular Dynamics Study

Contrast agents are important imaging probes in clinical MRI, allowing the identification of anatomic changes that otherwise would not be possible. Intensive research on the development of new contrast agents is being made to image specific pathological markers or sense local biochemical changes. The most widely used MRI contrast agents are based on gadolinium(III) complexes. Due to their very high charge density, they have low permeability through tight biological barriers such as the blood-brain barrier, hampering their application in the diagnosis of neurological disorders. In this study, we explore the interaction between the widely used contrast agent [Gd(DOTA)]- (Dotarem) and POPC lipid bilayers by means of molecular dynamics simulations. This metal complex is a standard reference where several chemical modifications have been introduced to improve key properties such as bioavailability and targeting. The simulations unveil detailed insights into the agent's interaction with the lipid bilayer, offering perspectives beyond experimental methods. Various properties, including the impact on global and local bilayer properties, were analyzed. As expected, the results indicate a low partition coefficient (KP) and high permeation barrier for this reference compound. Nevertheless, favorable interactions are established with the membrane leading to moderately long residence times. While coordination of one inner-sphere water molecule is maintained for the membrane-associated chelate, the physical-chemical attributes of [Gd(DOTA)]- as a MRI contrast agent are affected. Namely, increases in the rotational correlation times and in the residence time of the inner-sphere water are observed, with the former expected to significantly increase the water proton relaxivity. This work establishes a reference framework for the use of simulations to guide the rational design of new contrast agents with improved relaxivity and bioavailability and for the development of liposome-based formulations for use as imaging probes or theranostic agents.

Insights into Transfer of Supramolecular Doxorubicin/Congo Red Aggregates through Phospholipid Membranes

Doxorubicin (DOX) is a commonly used chemotherapeutic drug, from the anthracycline class, which is genotoxic to neoplastic cells via a DNA intercalation mechanism. It is effective and universal; however, it also causes numerous side effects. The most serious of them are cardiotoxicity and a decrease in the number of myeloid cells. For this reason, targeted DOX delivery systems are desirable, since they would allow lowering the drug dose and therefore limiting systemic side effects. Recently, synthetic dyes, in particular Congo red (CR), have been proposed as possible DOX carriers. CR is a planar molecule, built of a central biphenyl moiety and two substituted naphthalene rings, connected with diazo bonds. In water, it forms elongated ribbon-shaped supramolecular structures, which are able to selectively interact with immune complexes. In our previous studies, we have shown that CR aggregates can intercalate DOX molecules. In this way, they preclude DOX precipitation in water solutions and increase its uptake by MCF7 breast cancer cells. In the present work, we further explore the interactions between DOX, CR, and their aggregates (CR/DOX) with phospholipid membranes. In addition to neutral molecules, the protonated doxorubicin form, DXP, is also studied. Molecular dynamics simulations are employed to study the transfer of CR, DOX, DXP, and their aggregates through POPC bilayers. Interactions of CR, DOX, and CR/DOX with model monolayers are studied with Langmuir trough measurements. This study shows that CR may support the transfer of doxorubicin molecules into the bilayer. Both electrostatic and van der Waals interactions with lipids are important in this respect. The former promote the initial stages of the insertion process, the latter keep guest molecules inside the bilayer.

Long-Chain Lipids Facilitate Insertion of Large Nanoparticles into Membranes of Small Unilamellar Vesicles

Insertion of hydrophobic nanoparticles into phospholipid bilayers is limited to small particles that can incorporate into a hydrophobic membrane core between two lipid leaflets. Incorporation of nanoparticles above this size limit requires the development of challenging surface engineering methodologies. In principle, increasing the long-chain lipid component in the lipid mixture should facilitate incorporation of larger nanoparticles. Here, we explore the effect of incorporating very long phospholipids (C24:1) into small unilamellar vesicles on the membrane insertion efficiency of hydrophobic nanoparticles that are 5-11 nm in diameter. To this end, we improve an existing vesicle preparation protocol and utilized cryogenic electron microscopy imaging to examine the mode of interaction and evaluate the insertion efficiency of membrane-inserted nanoparticles. We also perform classical coarse-grained molecular dynamics simulations to identify changes in lipid membrane structural properties that may increase insertion efficiency. Our results indicate that long-chain lipids increase the insertion efficiency by preferentially accumulating near membrane-inserted nanoparticles to reduce the thermodynamically unfavorable disruption of the membrane.

Non-affine deformation analysis and 3D packing defects: A new way to probe membrane heterogeneity in molecular simulations

Here, we discuss a new framework developed over the last 5 years in our group to probe nanoscale membrane heterogeneity. The framework is based on the idea of characterizing lateral heterogeneity through non-affine deformation (NAD) measurements, transverse heterogeneity through three dimensional (3D) lipid packing defects, and using these approaches to formalize the seemingly trivial correlation between lateral organization and lipid packing in biological membranes. We find that measurements from NAD analysis, a prescription which is borrowed from Physics of glasses and granular material, can faithfully distinguish between liquid-ordered and disordered phases in membranes at molecular length scales and, can also be used to identify phase boundaries with high precision. Concomitantly, 3D-packing defects can not only distinguish between the two co-existing fluid phases based on their molecular scale packing (or membrane free volume), but also provide a route to connect the membrane domains to their functionality, such as exploring the molecular origins of inter-leaflet domain registration and peptide partitioning. The correlation between lateral membrane order and transverse packing presents novel molecular design-level features that can explain functions such as protein/peptide partitioning and small-molecule permeation dynamics in complex and heterogeneous membranes with high-fidelity. The framework allows us to explore the nature of lateral organization and molecular packing as a manifestation of intricate molecular interactions among a chemically rich variety of lipids and other molecules in a membrane with complex membrane composition and asymmetry across leaflets.

Studying the Effects of Dissolved Noble Gases and High Hydrostatic Pressure on the Spherical DOPC Bilayer Using Molecular Dynamic Simulations

Fine-grained molecular dynamics simulations have been conducted to depict lipid objects enclosed in water and interacting with a series of noble gases dissolved in the medium. The simple point-charge (SPC) water system, featuring a boundary composed of 1,2-Dioleoyl-sn-glycero-3-phosphocholine (DOPC) molecules, maintained stability throughout the simulation under standard conditions. This allowed for the accurate modeling of the effects of hydrostatic pressure at an ambient pressure of 25 bar. The chosen pressure references the 240 m depth of seawater: the horizon frequently used by commercial divers, who comprise the primary patient population of the neurological complication of inert gas narcosis and the consequences of high-pressure neurological syndrome. To quantify and validate the neurological effects of noble gases and discriminate them from high hydrostatic pressure, we reduced the dissolved gas molar concentration to 1.5%, three times smaller than what we previously tested for the planar bilayer (3.5%). The nucleation and growth of xenon, argon and neon nanobubbles proved consistent with the data from the planar bilayer simulations. On the other hand, hyperbaric helium induces only a residual distorting effect on the liposome, with no significant condensed gas fraction observed within the hydrophobic core. The bubbles were distributed over a large volume-both in the bulk solvent and in the lipid phase-thereby causing substantial membrane distortion. This finding serves as evidence of the validity of the multisite distortion hypothesis for the neurological effect of inert gases at high pressure.

Functional regulation of aquaporin dynamics by lipid bilayer composition

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics are examined. We demonstrate that SoPIP2;1's structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

Digestion of lipid micelles leads to increased membrane permeability

Lipid-based drug carriers are an attractive option to solubilise poorly water soluble therapeutics. Previously, we reported that the digestion of a short tail PC lipid (2C6PC) by the PLA2 enzyme has a significant effect on the structure and stability of the micelles it forms. Here, we studied the interactions of micelles of varying composition representing various degrees of digestion with a model ordered (70 mol% DPPC & 30 mol% cholesterol) and disordered (100% DOPC) lipid membrane. Micelles of all compositions disassociated when interacting with the two different membranes. As the percentage of digestion products (C6FA and C6LYSO) in the micelle increased, the disassociation occurred more rapidly. The C6FA inserts preferentially into both membranes. We find that all micelle components increase the area per lipid, increase the disorder and decrease the thickness of the membranes, and the 2C6PC lipid molecules have the most significant impact. Additionally, there is an increase in permeation of water into the membrane that accompanies the insertion of C6FA into the DOPC membranes. We show that the natural digestion of lipid micelles result in molecular species that can enhance the permeability of lipid membranes that in turn result in an enhanced delivery of drugs.

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.

Lipid and cholesterols modulate the dynamics of SARS-CoV-2 viral ion channel ORF3a and its pathogenic variants

SARS-CoV-2 accessory protein, ORF3a is a putative ion channel which immensely contributes to viral pathogenicity by modulating host immune responses and virus-host interactions. Relatively high expression of ORF3a in diseased individuals and implication with inflammasome activation, apoptosis and autophagy inhibition, ratifies as an effective target for developing vaccines and therapeutics. Herein, we present the elusive dynamics of ORF3a-dimeric state using all-atoms molecular dynamics (MD) simulations at μ-seconds scale in a heterogeneous lipid-mimetic system in multiple replicates. Additionally, we also explore the effect of non-synonymous pathogenic mutations on ORF3a ion channel activity and viral pathogenicity in different SARS-CoV-2 variants using various structure-based protein stability (ΔΔG) tools and computational saturation mutagenesis. Our study ascertains the role of phosphatidylcholines and cholesterol in modulating the structure of ORF3a, which perturbs the size and flexibility of the polar cavity that allows permeation of large cations. Discrete trend in ion channel pore radius and area per lipid arises the premise that presence of lipids might also affect the overall conformation of ORF3a. MD structural-ensembles, in some replicates rationalize the crucial role of TM2 in maintaining the native structure of ORF3a. We also infer that loss of structural stability primarily grounds for pathogenicity in more than half of the pathogenic variants of ORF3a. Overall, the effect of mutation on alteration of ion permeability of ORF3a, proposed in this study brings mechanistic insights into variant consequences on viral membrane proteins of SARS-CoV-2, which can be utilized for the development of novel therapeutics to treat COVID-19 and other coronavirus diseases.

Probing the Hydrophobic Region of a Lipid Bilayer at Specific Depths Using Vibrational Spectroscopy

A novel spectroscopic approach for studying the flexibility and mobility in the hydrophobic interior of lipid bilayers at specific depths is proposed. A set of test compounds featuring an azido moiety and a cyano or carboxylic acid moiety, connected by an alkyl chain of different lengths, was synthesized. FTIR data and molecular dynamics calculations indicated that the test compounds in a bilayer are oriented so that the cyano or carboxylic acid moiety is located in the lipid head-group region, while the azido group stays inside the bilayer at the depth determined by its alkyl chain length. We found that the asymmetric stretching mode of the azido group (νN3) can serve as a reporter of the membrane interior dynamics. FTIR and two-dimensional infrared (2DIR) studies were performed at different temperatures, ranging from 22 to 45 °C, covering the Lβ-Lα phase transition temperature of dipalmitoylphosphatidylcholine (∼41 °C). The width of the νN3 peak was found to be very sensitive to the phase transition and to the temperature in general. We introduced an order parameter, SN3, which characterizes restrictions to motion inside the bilayer. 2DIR spectra of νN3 showed different extents of inhomogeneity at different depths in the bilayer, with the smallest inhomogeneity in the middle of the leaflet. The spectral diffusion dynamics of the N3 peak was found to be dependent on the depth of the N3 group location in the bilayer. The obtained results enhance our understanding of the bilayer dynamics and can be extended to investigate membranes with more complex compositions.

Topological Learning Approach to Characterizing Biological Membranes

Biological membranes play key roles in cellular compartmentalization, structure, and its signaling pathways. At varying temperatures, individual membrane lipids sample from different configurations, a process that frequently leads to higher-order phase behavior and phenomena. Here we present a persistent homology-based method for quantifying the structural features of individual and bulk lipids, providing local and contextual information on lipid tail organization. Our method leverages the mathematical machinery of algebraic topology and machine learning to infer temperature-dependent structural information of lipids from static coordinates. To train our model, we generated multiple molecular dynamics trajectories of DPPC membranes at varying temperatures. A fingerprint was then constructed for each set of lipid coordinates by a persistent homology filtration, in which interactions spheres were grown around the lipid atoms while tracking their intersections. The sphere filtration formed a simplicial complex that captures enduring key topological features of the configuration landscape, using homology, yielding persistence data. Following fingerprint extraction for physiologically relevant temperatures, the persistence data were used to train an attention-based neural network for assignment of effective temperature values to selected membrane regions. Our persistence homology-based method captures the local structural effects, via effective temperature, of lipids adjacent to other membrane constituents, e.g. sterols and proteins. This topological learning approach can predict lipid effective temperatures from static coordinates across multiple spatial resolutions. The tool, called MembTDA, can be accessed at https://github.com/hyunp2/Memb-TDA.

Hop, Skip, and Jump: Hydrogen Molecular Transport through Amorphous Polyethylene Matrices Studied via Molecular Dynamics Simulations

In the pursuit of advancing and diversifying energy technologies for a more sustainable future, the possibilities of hydrogen (H2) usage will broaden, as will our understanding of its containment materials. Polyethylene (PE) has a vast assortment of uses and applications, which are growing with demands for alternative energy possibilities. One use of PE liner is as a prime candidate for nonmetallic piping and pressurized type IV storage devices. Such applications require PE to effectively prevent H2 transport through containment systems. To study the molecular transport mechanism of hydrogen through polymeric barriers, a system containing hydrogen molecules absorbed within amorphous PE is modeled here using molecular dynamics simulations. The simulations are conducted within a range of temperatures that span the glass transition temperature of amorphous PE. The simulated PE displays bulk density, radius of gyration, and self-diffusion coefficient that are consistent with experimental data. The simulated trajectories are interrogated to study the movement of the guest gas molecules. The results show that the diffusion coefficients increase with temperature, as expected. However, the mobility of the PE chains is found to affect the mobility of absorbed H2 molecules to a much lower extent than it affects that of CH4 molecules because of the much smaller size of the former than of the latter guest. From a molecular perspective, a "hopping" mechanism is observed, according to which H2 molecules hop between one vacant free volume space to another within the polymer matrix, in combination with longer, straight, undisturbed "jumps" or "skips" along directions aligned with regions of ordered PE chains. This suggests that the orientation of the crystallites within the semicrystalline PE matrix affects the H2 containment. Implications of these findings toward PE usage as containment material are discussed.

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.

Integration of an LPAR1 Antagonist into Liposomes Enhances Their Internalization and Tumor Accumulation in an Animal Model of Human Metastatic Breast Cancer

Lysophosphatidic acid receptor 1 (LPAR1) is elevated in breast cancer. The deregulation of LPAR1, including the function and level of expression, is linked to cancer initiation, progression, and metastasis. LPAR1 antagonists, AM095 or Ki16425, may be effective therapeutic molecules, yet their limited water solubility hinders in vivo delivery. In this study, we report on the synthesis of two liposomal formulations incorporating AM095 or Ki16425, embedded within the lipid bilayer, as targeted nanocarriers for metastatic breast cancer (MBC). The data show that the Ki16425 liposomal formulation exhibited a 50% increase in internalization by MBC mouse epithelial cells (4T1) and a 100% increase in tumor accumulation in a mouse model of MBC compared with that of a blank liposomal formulation (control). At the same time, normal mouse epithelial cells (EpH-4Ev) internalized the Ki16425 liposomal formulation 25% lesser than the control formulation. Molecular dynamics simulations show that the integration of AM095 or Ki16425 modified the physical and mechanical properties of the lipid bilayer, making it more flexible in these liposomal formulations compared with liposomes without drug. The incorporation of an LPAR1 antagonist within a liposomal drug delivery system represents a viable therapeutic approach for targeting the LPA-LPAR1 axis, which may hinder the progression of MBC.

Hydrophobic Matching Dictates over the Linear Rule of Mixtures in Binary Lipid Membranes

Biological membranes feature heterogeneous mixtures of lipids with different head and tail characteristics. Their biophysical properties are dictated by the intimate interaction among different constituent lipids. Previous studies suggest that the membrane area-per-lipid (APL) deviates from the linear rule of mixtures (LRM) for binary lipid membranes, but the underlying mechanism remains elusive. Our molecular dynamics (MD) simulations of binary lipid membranes consisting of lipids with different tail characteristics reveal a competitive mechanism whereby lipids tend to deform each other to minimize the hydrophobic mismatch between their tails. Depending on the relative tail lengths and saturation levels, this may result in an either positive or negative deviation of APL from the LRM. As lipid packing plays an essential role in membrane fusion and peptide-membrane binding, our findings may help guide the selection of lipids for the effective rational design of nanoliposomes and membrane-targeting peptides.

Lipid composition modulates interactions of p7 viroporin during membrane insertion

Viral proteins interact with lipid membranes during various stages in the viral life cycle to propagate infection. p7 is an ion channel forming protein of Hepatitis C virus (HCV) that participates in viral assembly. Studies show that it has close ties to lipid metabolism in the cell and anionic phosphatidylserine (PS) lipids are suggested to be key for its permeabilizing function, but the mechanism of its interaction with the lipid environment is largely unknown. To begin unraveling the molecular processes of the protein, we evaluated the impact of lipid environment on the binding and insertion mechanism of p7 prior to channel formation and viral assembly using molecular dynamics simulations. It is seen that p7 is sensitive to its lipid environment and results in different remodeling patterns in membranes. Helix 1 (H1) is especially important for peptide insertion, with deeper entry taking place when the membrane contains phosphatidylserine (PS). Helix 2 (H2) and the adjacent loop connecting to Helix 3 (H3) prompts recruitment of phosphatidylethanolamine (PE) lipids to the protein binding site in membrane models with lower surface charge. This work provides perspectives on the interplay between protein-lipid dynamics and membrane composition, and insights on membrane reorganization in mechanisms of disease.

Functional Regulation of Aquaporin Dynamics by Lipid Bilayer Composition

With the diversity of lipid-protein interactions, any observed membrane protein dynamics or functions directly depend on the lipid bilayer selection. However, the implications of lipid bilayer choice are seldom considered unless characteristic lipid-protein interactions have been previously reported. Using molecular dynamics simulation, we characterize the effects of membrane embedding on plant aquaporin SoPIP2;1, which has no reported high-affinity lipid interactions. The regulatory impacts of a realistic lipid bilayer, and nine different homogeneous bilayers, on varying SoPIP2;1 dynamics were examined. We demonstrate that SoPIP2;1s structure, thermodynamics, kinetics, and water transport are altered as a function of each membrane construct's ensemble properties. Notably, the realistic bilayer provides stabilization of non-functional SoPIP2;1 metastable states. Hydrophobic mismatch and lipid order parameter calculations further explain how lipid ensemble properties manipulate SoPIP2;1 behavior. Our results illustrate the importance of careful bilayer selection when studying membrane proteins. To this end, we advise cautionary measures when performing membrane protein molecular dynamics simulations.

Characterizing the Interactions of Cell-Membrane-Disrupting Peptides with Lipid-Functionalized Single-Walled Carbon Nanotubes

Lipid-functionalized single-walled carbon nanotubes (SWNTs) have garnered significant interest for their potential use in a wide range of biomedical applications. In this work, we used molecular dynamics simulations to study the equilibrium properties of SWNTs surrounded by the phosphatidylcholine (POPC) corona phase and their interactions with three cell membrane disruptor peptides: colistin, TAT peptide, and crotamine-derived peptide. Our results show that SWNTs favor asymmetrical positioning within the POPC corona, so that one side of the SWNT, covered by the thinnest part of the corona, comes in contact with charged and polar functional groups of POPC and water. We also observed that colistin and TAT insert deeply into the POPC corona, while crotamine-derived peptide only adsorbs to the corona surface. In separate simulations, we show that three examined peptides exhibit similar insertion and adsorption behaviors when interacting with POPC bilayers, confirming that peptide-induced perturbations to POPC in conjugates and bilayers are similar in nature and magnitude. Furthermore, we observed correlations between the peptide-induced structural perturbations and the near-infrared emission of the lipid-functionalized SWNTs, which suggest that the optical signal of the conjugates transduces the morphological changes in the lipid corona. Overall, our findings indicate that lipid-functionalized SWNTs could serve as simplified cell membrane model systems for prescreening of new antimicrobial compounds that disrupt cell membranes.

Insights into Early Steps of Decanoic Acid Self-Assemblies under Prebiotic Temperatures Using Molecular Dynamics Simulations

The origin of life possibly required processes in confined systems that facilitated simple chemical reactions and other more complex reactions impossible to achieve under the condition of infinite dilution. In this context, the self-assembly of micelles or vesicles derived from prebiotic amphiphilic molecules is a cornerstone in the chemical evolution pathway. A prime example of these building blocks is decanoic acid, a short-chain fatty acid capable of self-assembling under ambient conditions. This study explored a simplified system made of decanoic acids under temperatures ranging from 0 °C to 110 °C to replicate prebiotic conditions. The study revealed the first point of aggregation of decanoic acid into vesicles and examined the insertion of a prebiotic-like peptide in a primitive bilayer. The information gathered from this research provides critical insights into molecule interactions with primitive membranes, allowing us to understand the first nanometric compartments needed to trigger further reactions that were essential for the origin of life.

Fluorescent Probes cis- and trans-Parinaric Acids in Fluid and Gel Lipid Bilayers: A Molecular Dynamics Study

Fluorescence probes are indispensable tools in biochemical and biophysical membrane studies. Most of them possess extrinsic fluorophores, which often constitute a source of uncertainty and potential perturbation to the host system. In this regard, the few available intrinsically fluorescent membrane probes acquire increased importance. Among them, cis- and trans-parinaric acids (c-PnA and t-PnA, respectively) stand out as probes of membrane order and dynamics. These two compounds are long-chained fatty acids, differing solely in the configurations of two double bonds of their conjugated tetraene fluorophore. In this work, we employed all-atom and coarse-grained molecular dynamics simulations to study the behavior of c-PnA and t-PnA in lipid bilayers of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC), representative of the liquid disordered and solid ordered lipid phases, respectively. All-atom simulations indicate that the two probes show similar location and orientation in the simulated systems, with the carboxylate facing the water/lipid interface and the tail spanning the membrane leaflet. The two probes establish interactions with the solvent and lipids to a similar degree in POPC. However, the almost linear t-PnA molecules have tighter lipid packing around them, especially in DPPC, where they also interact more with positively charged lipid choline groups. Probably for these reasons, while both probes show similar partition (assessed from computed free energy profiles across bilayers) to POPC, t-PnA clearly partitions more extensively than c-PnA to the gel phase. t-PnA also displays more hindered fluorophore rotation, especially in DPPC. Our results agree very well with experimental fluorescence data from the literature and allow deeper understanding of the behavior of these two reporters of membrane organization.

Characterizing the Interactions of Cell Membrane-Disrupting Peptides with Lipid-Functionalized Single-Walled Carbon Nanotube Systems for Antimicrobial Screening

Lipid-functionalized single-walled carbon nanotubes (SWNTs) have garnered significant interest for their potential use in a wide range of biomedical applications. In this work, we used molecular dynamics simulations to study the equilibrium properties of SWNTs surrounded by the phosphatidylcholine (POPC) corona phase, and their interactions with three cell membrane disruptor peptides: colistin, TAT peptide, and crotamine-derived peptide. Our results show that SWNTs favor asymmetrical positioning within the POPC corona, so that one side of the SWNT, covered by the thinnest part of the corona, comes in contact with charged and polar functional groups of POPC and water. We also observed that colistin and TAT insert deeply into POPC corona, while crotamine-derived peptide only adsorbs to the corona surface. Compared to crotamine-derived peptide, colistin and TAT also induce larger perturbations in the thinnest region of the corona, by allowing more water molecules to directly contact the SWNT surface. In separate simulations, we show that three examined peptides exhibit similar insertion and adsorption behaviors when interacting with POPC bilayers, confirming that peptide-induced perturbations to POPC in conjugates and bilayers are similar in nature and magnitude. Furthermore, we observed correlations between the peptide-induced structural perturbations and the near-infrared emission of the lipid-functionalized SWNTs, which suggest that the optical signal of the conjugates transduces the morphological changes in the lipid corona. Overall, our findings indicate that lipid-functionalized SWNTs could serve as simplified cell membrane model systems for pre-screening of new antimicrobial compounds that disrupt cell membranes.

Dehydration does not affect lipid-based hydration lubrication

Phosphatidylcholine (PC) lipid bilayers at surfaces massively reduce sliding friction, via the hydration lubrication mechanism acting at their highly-hydrated phosphocholine headgroups, a central paradigm of biological lubrication, particularly at articular cartilage surfaces where low friction is crucial for joint well-being. Nanotribological measurements probed the effect on such lubrication of dehydration by dimethyl sulfoxide (DMSO), known to strongly dehydrate the phosphocholine headgroups of such PC bilayers, i.e. reduce the thickness of the inter-bilayer water layer, and thus expected to substantially degrade the hydration lubrication. Remarkably, and unexpectedly, we found that the dehydration has little effect on the friction. We used several approaches, including atomic force microscopy, small- and wide-angle X-ray scattering and all-atom molecular dynamics simulations to elucidate this. Our results show that while DMSO clearly removes hydration water from the lipid head-groups, this is offset by both higher areal head-group density and by rigidity-enhancement of the lipid bilayers, both of which act to reduce frictional dissipation. This sheds strong light on the robustness of lipid-based hydration lubrication in biological systems, despite the ubiquitous presence of bio-osmolytes which compete for hydration water.

Molecular dynamics simulations to study the role of biphenylalanine in promoting the antibacterial activity of ultrashort peptides

The hydrophobic amino acid biphenylalanine (B) plays a key role in the antibacterial activity of ultrashort peptides. In this study, the interactions of tetrapeptide BRBR-NH₂ (BRBR) and pentapeptide BRBRB-NH₂ (BRBRB) with dioleoylphosphatidylcholine/dioleoylphosphatidylglycerol (DOPC/DOPG) mixed model membrane were studied by molecular dynamics simulation to assess the role of biphenylalanine in promoting the antibacterial activity of ultrashort peptides. At low peptide concentrations, both peptides presented amphiphilic conformations; residues B of the pentapeptide approached the membrane faster than those of the tetrapeptide and made more contacts with the membrane; BRBRB exhibited stronger membrane affinity than BRBR. However, due to the low peptide concentrations, the effects of these two peptides on the membrane were not significantly different. At high peptide concentrations, the strong affinity of BRBRB made it have more interaction with membrane than BRBR and most residues B of BRBRB inserted into the membrane; BRBRB was more prone to aggregation and caused the membrane more disordered and thinner than BRBR. Hydrophobic residues often act as anchors in the antibacterial activity of ultrashort antimicrobial peptides. Adding a hydrophobic residue B to the C-terminal of BRBR could improve the ability of the peptide to "grasp" the membrane. At high peptide concentrations, the addition of residue B might enhance the antibacterial activity of the peptide. Thus, our results will be helpful in designing efficient antibacterial drugs.

Molecular simulation of lignin-related aromatic compound permeation through gram-negative bacterial outer membranes

Lignin, an abundant aromatic heteropolymer in secondary plant cell walls, is the single largest source of renewable aromatics in the biosphere. Leveraging this resource for renewable bioproducts through targeted microbial action depends on lignin fragment uptake by microbial hosts and subsequent enzymatic action to obtain the desired product. Recent computational work has emphasized that bacterial inner membranes are permeable to many aromatic compounds expected from lignin depolymerization processes. In this study, we expand on these findings through simulations for 42 lignin-related compounds across a gram-negative bacterial outer membrane model. Unbiased simulation trajectories indicate that spontaneous crossing for the full outer membrane is relatively rare at molecular simulation timescales, primarily due to preferential membrane partitioning and slow diffusion within the lipopolysaccharide layer within the outer membrane. Membrane partitioning and permeability coefficients were determined through replica exchange umbrella sampling simulations to overcome sampling limitations. We find that the glycosylated lipopolysaccharides found in the outer membrane increase the permeation barrier to many lignin-related compounds, particularly the most hydrophobic compounds. However, the effect is relatively modest; at industrially relevant concentrations, uncharged lignin-related compounds will readily diffuse across the outer membrane without the need for specific porins. Together, our results provide insight into the permeability of the bacterial outer membrane for assessing lignin fragment uptake and the future production of renewable bioproducts.

Magic mushroom extracts in lipid membranes

The active hallucinogen of magic mushrooms, psilocin, is being repurposed to treat nicotine addiction and treatment-resistant depression. Psilocin belongs to the tryptamine class of psychedelic compounds which include the hormone serotonin. It is believed that psilocin exerts its effect by binding to the serotonin 5-HT_2A receptor. However, recent in-vivo evidence suggests that psilocin may employ a different mechanism to exert its effects. Membrane-mediated receptor desensitization of neurotransmitter receptors is one such mechanism. We compare the impact of the neutral and charged versions of psilocin and serotonin on the properties of zwitterionic and anionic lipid membranes using molecular dynamics simulations and calorimetry. Both compounds partition to the lipid interface and induce membrane thinning. The tertiary amine in psilocin, as opposed to the primary amine in serotonin, limits psilocin's impact on the membrane although more psilocin partitions into the membrane than serotonin. Calorimetry corroborates that both compounds induce a classical melting point depression like anesthetics do. Our results also lend support to a membrane-mediated receptor-binding mechanism for both psilocin and serotonin and provide physical insights into subtle chemical changes that can alter the membrane-binding of psychedelic compounds.

Using molecular simulation to understand the skin barrier

Skin's effectiveness as a barrier to permeation of water and other chemicals rests almost entirely in the outermost layer of the epidermis, the stratum corneum (SC), which consists of layers of corneocytes surrounded by highly organized lipid lamellae. As the only continuous path through the SC, transdermal permeation necessarily involves diffusion through these lipid layers. The role of the SC as a protective barrier is supported by its exceptional lipid composition consisting of ceramides (CERs), cholesterol (CHOL), and free fatty acids (FFAs) and the complete absence of phospholipids, which are present in most biological membranes. Molecular simulation, which provides molecular level detail of lipid configurations that can be connected with barrier function, has become a popular tool for studying SC lipid systems. We review this ever-increasing body of literature with the goals of (1) enabling the experimental skin community to understand, interpret and use the information generated from the simulations, (2) providing simulation experts with a solid background in the chemistry of SC lipids including the composition, structure and organization, and barrier function, and (3) presenting a state of the art picture of the field of SC lipid simulations, highlighting the difficulties and best practices for studying these systems, to encourage the generation of robust reproducible studies in the future. This review describes molecular simulation methodology and then critically examines results derived from simulations using atomistic and then coarse-grained models.

Generation and Computational Characterization of a Complex Staphylococcus aureus Lipid Bilayer

Studies indicate a crucial cell membrane role in the antibiotic resistance of Staphylococcus aureus. To simulate its membrane structure and dynamics, a complex molecular-scale computational representation of the S. aureus lipid bilayer was developed. Phospholipid types and their amounts were optimized by reverse Monte Carlo to represent characterization data from the literature, leading to 19 different phospholipid types that combine three headgroups [phosphatidylglycerol, lysyl-phosphatidylglycerol (LPG), and cardiolipin] and 10 tails, including iso- and anteiso-branched saturated chains. The averaged lipid bilayer thickness was 36.7 Å, and area per headgroup was 67.8 Ų. Phosphorus and nitrogen density profiles showed that LPG headgroups tended to be bent and oriented more parallel to the bilayer plane. The water density profile showed that small amounts reached the membrane center. Carbon density profiles indicated hydrophobic interactions for all lipids in the middle of the bilayer. Bond vector order parameters along each tail demonstrated different C-H ordering even within distinct lipids of the same type; however, all tails followed similar trends in average order parameter. These complex simulations further revealed bilayer insights beyond those attainable with monodisperse, unbranched lipids. Longer tails often extended into the opposite leaflet. Carbon at and beyond a branch showed significantly decreased ordering compared to carbon in unbranched tails; this feature arose in every branched lipid. Diverse tail lengths distributed these disordered methyl groups throughout the middle third of the bilayer. Distributions in mobility and ordering reveal diverse properties that cannot be obtained with monodisperse lipids.

Coarse-grained description of monounsaturated peroxidized phospholipidbilayers

Lipid peroxides result from a reaction between cis-unsaturated lipid chains and singlet oxygen molecules leading to the addition of a peroxide OOH side group next to the acyl-chain double bond. It is now established that HP-POPC (hydroperoxidized POPC) molecules form stable, thin and laterally expanded bilayers. The difference in structural organisation arises from the hydrophilic character of the OOH side group that has strong affinity with the water interface region, leading to significant reorganisation of the bilayer. In this Article we describe a coarse-grained (CG) model of POPC and DOPC lipid peroxides within the framework of the Martini CG force-field (v2.2), derived from experimental data. We then discuss extensively the predicted structure and the influence of hydration and show how shifting the position of the unsaturated bonds along the chain changes the structure. Finally, we provide electron and neutron scattering length density profiles of the simulated bilayers.

Probing the Link between Pancratistatin and Mitochondrial Apoptosis through Changes in the Membrane Dynamics on the Nanoscale

Pancratistatin (PST) is a natural antiviral alkaloid that has demonstrated specificity toward cancerous cells and explicitly targets the mitochondria. PST initiates apoptosis while leaving healthy, noncancerous cells unscathed. However, the manner by which PST induces apoptosis remains elusive and impedes the advancement of PST as a natural anticancer therapeutic agent. Herein, we use neutron spin-echo (NSE) spectroscopy, molecular dynamics (MD) simulations, and supporting small angle scattering techniques to study PST's effect on membrane dynamics using biologically representative model membranes. Our data suggests that PST stiffens the inner mitochondrial membrane (IMM) by being preferentially associated with cardiolipin, which would lead to the relocation and release of cytochrome c. Second, PST has an ordering effect on the lipids and disrupts their distribution within the IMM, which would interfere with the maintenance and functionality of the active forms of proteins in the electron transport chain. These previously unreported findings implicate PST's effect on mitochondrial apoptosis.

Alteration of lipid bilayer mechanics by volatile anesthetics: Insights from μs-long molecular dynamics simulations

Very few drugs in clinical practice feature the chemical diversity, narrow therapeutic window, unique route of administration, and reversible cognitive effects of volatile anesthetics. The correlation between their hydrophobicity and their potency and the increasing amount of evidence suggesting that anesthetics exert their action on transmembrane proteins, justifies the investigation of their effects on phospholipid bilayers at the molecular level, given the strong functional and structural link between transmembrane proteins and the surrounding lipid matrix. Molecular dynamics simulations of a model lipid bilayer in the presence of ethylene, desflurane, methoxyflurane, and the nonimmobilizer 1,2-dichlorohexafluorocyclobutane (also called F6 or 2N) at different concentrations highlight the structural consequences of VA partitioning in the lipid phase, with a decrease of lipid order and bilayer thickness, an increase in overall lipid lateral mobility and area-per-lipid, and a marked reduction in the mechanical stiffness of the membrane, that strongly correlates with the compounds' hydrophobicity.

Analysis of the Conformational Dynamics of the Soluble and Membrane-Bound Forms of Interleukin-1 Receptor Type-1: Insights into Linker Flexibility and Domain Rearrangements

Interleukin-1 receptor type 1 (IL-1R1) is a key player in inflammation and immune responses. This receptor regulates IL-1 activity in two forms: as a membrane-bound form and as a soluble ectodomain. The details and differences between the conformational dynamics of the membrane-bound and the soluble IL-1R1 ectodomains (ECDs) remain largely elusive. Here, we study and compare the structural dynamics of the soluble and membrane-bound IL-1R1-ECDs using molecular dynamics (MD) simulations, focusing on the flexible interdomain linker of the ECD, as well as the spatial rearrangements between the Ig-like domains of the ECD. To explore the membrane-bound conformations, a full-length IL-1R1 structural model was developed and subjected to classical equilibrium MD. Comparative analysis of multiple MD trajectories of the soluble and the membrane-bound IL-1R1-ECDs reveals that (i) as somewhat expected, the extent of the visited “open-to-closed” transitional states differs significantly between the soluble and membrane-bound forms; (ii) the soluble form presents open-closed transitions, sampling a wider rotational motion between the Ig-like domains of the ECD, visiting closed and “twisted” conformations in higher extent, whereas the membrane-bound form is characterized by more conformationally restricted states; (iii) interestingly, the backbone dihedral angles of residues Glu202, Glu203 and Asn204, located in the flexible linker, display the highest variations during the transition between discrete conformational states detected in IL-1R1, thus appearing to work as the “central wheel of a clock’s movement”. The simulations and analyses presented in this contribution offer a deeper insight into the structure and dynamics of IL-1R1, which may be explored in a drug discovery setting.

PE and PET oligomers' interplay with membrane bilayers

The prevalence of microplastic pollution in nature and foodstuffs is fairly well identified. However, studies of micro- or nanoplastics' cell membrane permeation and health effects in humans are lacking. Our study focuses on examining the interactions of polyethylene (PE) and polyethylene terephthalate (PET) with bilayer membranes. We have performed molecular dynamics simulations to study how plastic oligomers behave in bilayers. In addition, we have studied membrane permeation of PE and Bis(2-hydroxyethyl) terephthalate (BHET), a type of PET monomer, with Parallel Artificial Membrane Permeability Assay (PAMPA). As a result, in simulations the molecules exhibited different movements and preferred locations in membrane. PAMPA studies suggested similar preferences in membrane, especially for PE plastic. Our results suggest that passive diffusion could be an important transport mechanism into cells for some small plastic oligomers. Both molecular dynamics simulations and PAMPA have potential for micro- and nanoplastics research.

Drug Meets Monolayer: Understanding the Interactions of Sterol Drugs with Models of the Lung Surfactant Monolayer Using Molecular Dynamics Simulations

The lung surfactant monolayer (LSM) is a thin layer of lipids and proteins that forms the air/water interface of the alveoli. The primary function of the LSM is to reduce the surface tension at the air/water interface during breathing. The LSM also forms the main biological barrier for any inhaled particles, including drugs, to treat lung diseases. Elucidating the mechanism by which these drugs bind to and absorb into the LSM requires a molecular-level understanding of any drug-induced changes to the morphology, structure, and phase changes of the LSM.Molecular dynamics simulations have been used extensively to study the structure and dynamics of the LSM. The monolayer is usually simulated in at least two states: the compressed state, mimicking exhalation, and the expanded state, mimicking inhalation. In this chapter, we provide detailed instructions on how to set up, run, and analyze coarse-grained MD simulations to study the concentration-dependent effect of a sterol drug on the LSM, both in the expanded and compressed state.

Molecular Organization of a Raft-like Domain in a Polyunsaturated Phospholipid Bilayer: A Supervised Machine Learning Analysis of Molecular Dynamics Simulations

Numerous health benefits are associated with omega-3 polyunsaturated fatty acids (n-3 PUFA) consumed in fish oils. An understanding of the mechanism remains elusive. The plasma membrane as a site of action is the focus in this study. With large-scale all-atom MD simulations run on a model membrane (1050 lipid molecules), we observed the evolution over time (6 μs) of a circular (raft-like) domain composed of N-palmitoylsphingomyelin (PSM) and cholesterol embedded into a surrounding (non-raft) patch composed of polyunsaturated 1-palmitoyl-2-docosahexaenoylphosphatylcholine (PDPC) (1:1:1 mol). A supervised machine learning algorithm was developed to characterize the migration of each lipid based on molecular conformation and the local environment. PDPC molecules were seen to infiltrate the ordered raft-like domain in a small amount, while a small concentration of PSM and cholesterol molecules was seen to migrate into the disordered non-raft region. Enclosing the raft-like domain, a narrow (∼2 nm in width) interfacial zone composed of PDPC, PSM, and cholesterol that buffers the substantial difference in order (ΔS_CD ≈ 0.12) between raft-like and non-raft environments was seen to form. Our results suggest that n-3 PUFA regulate the architecture of lipid rafts enriched in sphingolipids and cholesterol with a minimal effect on order within their interior in membranes.

The role of size and nature in nanoparticle binding to a model lung membrane: an atomistic study

Understanding the uptake of nanoparticles (NPs) by different types of cellular membranes plays a pivotal role in the design of NPs for medical applications and in avoiding adverse effects that result in nanotoxicity. Yet, the role of key design parameters, such as the bare NP material, NP size and surface reactivity, and the nature of NP coatings, in membrane remodelling and uptake mechanisms is still very poorly understood, particularly towards the lower range of NP dimensions that are beyond the experimental imaging resolution. The same can be said about the role of a particular membrane composition. Here, we systematically employ biased and unbiased molecular dynamics simulations to calculate the binding energy for three bare materials (Ag/SiO₂/TiO₂) and three NP sizes (1/3/5 nm diameter) with a representative lung surfactant membrane, and to study their binding kinetics. The calculated binding energies show that irrespective of size, Ag nanoparticles bind very strongly to the bilayer, while the NPs made of SiO₂ or TiO₂ experience very low to no binding. The unbiased simulations provide insight into how the NPs and membrane affect each other in terms of the solvent-accessible surface area (SASA) of the NPs and the defect types and fluidity of the membrane. Using these systematic fine-grained results in coarsening procedures will pave the way for simulations considering NP sizes that are well beyond the membrane thickness, i.e. closer to experimental dimensions, for which different binding characteristics and more significant membrane remodelling are expected.

Interaction of Human β Defensin Type 3 (hBD-3) with Different PIP2-Containing Membranes, a Molecular Dynamics Simulation Study

Human β defensin type 3 (hBD-3) is a cysteine-rich small antibacterial peptide. It belongs to the human innate immune system. hBD-3 has six cysteine residues, which form three pairs of disulfide bonds, and those bonds break in the reducing condition. It is known that hBD-3 can interact with bacterial membrane, and even eukaryotic cell membrane, which has a low concentration of phosphatidylinositol 4,5-bisphosphate (PIP2) lipids. PIP2 is a vital component in cell membranes and has been found to play important roles during antimicrobial peptide (AMP) interaction with membranes. To understand the functional mechanism of hBD-3 interacting with PIP2-containing membranes, the binding structures of hBD-3 on 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayers mixed with 10% of PIP2 were predicted using two kinds of methods. The first one is by placing the hBD-3 monomer in different orientations above the POPC + 10%PIP2 membrane to set up five different initial simulation systems and performing long-term simulations on each to predict the most stable binding structure. It was found that hBD-3 analogue binds on the mixed lipid membrane on the two loop regions. The second method is by running long-term simulations on one or nine hBD-3 dimers binding on POPC mixed with 10%PIP2 lipid bilayer starting from the solid-state NMR (ssNMR)-suggested orientation. The dimer dissociated, and the most stable binding of hBD-3 in wild-type on the mixed membrane is also through the two loop regions, which agrees with the prediction from both the first method and the lipid self-assembly result. The PIP2 lipids can form long-lasting hydrogen bonds with positively charged residues such as Arg and Lys on hBD-3, thus forming clusters with hBD-3. As a comparison, hBD-3 dimers binding with a combined bilayer having 1,2-palmitoyl-oleoyl-sn-glycero-3-phosphoserine (POPS) on the upper and POPC on the lower leaflets and the combined POPS + POPC bilayer mixing with 10%PIP2 were also studied. The long-term simulation result shows that hBD-3 can bind with the heads of negatively charged POPS and PIP2 lipids and form hydrogen bonds. The stable binding sites of hBD-3 on PIP2 or POPS mixed bilayers are still on the two loop regions. On the combined POPS + POPC mixed with 10%PIP2 bilayer, the binding of hBD-3 with PIP2 lipids became less stable and fewer because of the competition of binding with the POPS lipids. Besides that, binding with hBD-3 can decrease the membrane thickness of the POPC + PIP2, POPS + POPC, and POPS + POPC + PIP2 bilayers and make POPS and PIP2 lipids more flexible based on the order parameter calculations. Our results supply molecular insight on AMP binding with different membranes and can help understand the functional mechanism of hBD-3 disrupting PIP2-containing membranes.

Unravelling the Effects of Cholesterol on the Second-Order Nonlinear Optical Responses of Di-8-ANEPPS Dye Embedded in Phosphatidylcholine Lipid Bilayers

Cholesterol is known for its role in maintaining the correct fluidity and rigidity of the animals cell membranes and thus their functions. Assessing the content and the role of cholesterol in lipid bilayers is therefore of crucial importance for a deeper understanding and control of membrane functioning. In this computational work, we investigate bilayers built from three types of glycerophospholipid phosphatidylcholine (PC) lipids, namely dipalmitoylphosphatidylcholine (DPPC), 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), and dioleoylphosphatidylcholine (DOPC), and containing different amounts of cholesterol by analyzing the second-harmonic generation (SHG) nonlinear optical (NLO) response of a probe molecule, di-8-ANEPPS, inserted into the membranes. This molecular property presents the advantage to be specific to interfacial regions such as lipid bilayers. To unravel these effects, Molecular Dynamics (MD) simulations have been performed on both DPPC and DOPC lipids by varying the cholesterol mole fraction (from 0 to 0.66), while POPC was only considered as a pure bilayer. In the case of the structural properties of the bilayers, all the analyses converge toward the same conclusion: as the mole fraction of cholesterol increases, the systems become more rigid, confirming the condensing effect of cholesterol. In addition, the chromophore is progressively more aligned with respect to the normal to the bilayer. On the contrary, addition of unsaturation disorders the lipid bilayers, with barely no impact on the alignment of the chromophore. Then, using the frames obtained from the MD simulations, the first hyperpolarizability β of the dye in its environment has been computed at the TDDFT level. On the one hand, the addition of cholesterol induces a progressive increase of the diagonal component the β tensor parallel to the bilayer normal. On the other hand, larger β values have been calculated for the unsaturated than for the saturated lipid systems. In summary, this study illustrates the relationship between the composition and structure of the bilayers and the NLO responses of the embedded dye.

Modeling lipid layers of atopic skin and observation of changes in lipid layer properties with changes in ceramide content

BACKGROUND

Studies have shown that there is a high correlation between atopic dermatitis and decrease in ceramide content in the lipid bilayer of skin. Moreover, it has been shown that the reduction in ceramide content in the stratum corneum is unique to atopic dermatitis, indicating that there are particular structural differences between the lipid bilayers of normal and atopic skin.

AIM

This study aimed to compare the lipid bilayer of the atopic skin with that of the healthy skin and to establish a structural model of the lipid bilayer for atopy.

METHODS

Molecular dynamics simulations were performed using NAMD 2.8. Models of lipid bilayers of normal skin and atopic skin, and a model of lipid bilayer containing only ceramide were built with CHARMM-GUI. The thickness, area occupied per lipid, and alignment of lipids were compared among the three models. Potential mean force (PMF) of the sodium laureth sulfate (SLES) on lipid bilayers was calculated to predict the affinity between SLES and lipid bilayers.

RESULTS

Potential mean force calculations showed that the lipid bilayer of atopic skin was able to absorb the surfactant more easily than that of normal skin.

CONCLUSIONS

When the ceramide ratio is low, the thickness of lipid bilayer is reduced and its structure is weakened. Other structural differences between the lipid layers of normal and atopic skin included increased area per lipid and poor alignment of lipids. Further, the atopy lipid bilayer model was found to absorb more SLES than the normal skin lipid bilayer model.

Unsaturated fatty acids augment protein transport via the SecA:SecYEG translocon

The translocon SecYEG and the associated ATPase SecA form the primary protein secretion system in the cytoplasmic membrane of bacteria. The secretion is essentially dependent on the surrounding lipids, but the mechanistic understanding of their role in SecA : SecYEG activity is sparse. Here, we reveal that the unsaturated fatty acids (UFAs) of the membrane phospholipids, including tetraoleoyl-cardiolipin, stimulate SecA : SecYEG-mediated protein translocation up to ten-fold. Biophysical analysis and molecular dynamics simulations show that UFAs increase the area per lipid and cause loose packing of lipid head groups, where the N-terminal amphipathic helix of SecA docks. While UFAs do not affect the translocon folding, they promote SecA binding to the membrane, and the effect is enhanced up to fivefold at elevated ionic strength. Tight SecA : lipid interactions convert into the augmented translocation. Our results identify the fatty acid structure as a notable factor in SecA : SecYEG activity, which may be crucial for protein secretion in bacteria, which actively change their membrane composition in response to their habitat.

How stereochemistry of lipid components can affect lipid organization and the route of liposome internalization into cells

Though liposome-based drugs are in clinical use, the mechanism of cell internalization of liposomes is yet an object of controversy. The present experimental investigation, carried out on human glioblastoma cells, indicated different internalization routes for two diastereomeric liposomes. Molecular dynamics simulations of the lipid bilayers of the two formulations indicated that the different stereochemistry of a lipid component controls some parameters such as area per lipid molecule and fluidity of lipid membranes, surface potential and water organization at the lipid/water interface, all of which affect the interaction with biomolecules and cell components.

General Protocol for Constructing Molecular Models of Nanodiscs

Nanodisc technology is increasingly being applied for structural and biophysical studies of membrane proteins. In this work, we present a general protocol for constructing molecular models of nanodiscs for molecular dynamics simulations. The protocol is written in python and based on geometric equations, making it fast and easy to modify, enabling automation and customization of nanodiscs in silico. The novelty being the ability to construct any membrane scaffold protein (MSP) variant fast and easy given only an input sequence. We validated and tested the protocol by simulating seven different nanodiscs of various sizes and with different membrane scaffold proteins, both circularized and noncircularized. The structural and biophysical properties were analyzed and shown to be in good agreement with previously reported experimental data and simulation studies.

Long-Chain Fatty Acids in Bones and Their Link to Submicroscopic Vascularization Network: NMR Assignment and Relaxation Studies under Magic Angle Spinning Conditions in Intramuscular Bones of Atlantic Herring Fish

The long-lasting proton signals in bones are identified as long-chain fatty acids, including saturated, mono-, and di-unsaturated fatty acids, with direct nuclear magnetic resonance evidence. We used intramuscular bones from Atlantic Herring fish to avoid interference from lipid-rich marrows. The key is to recognize that these signals are from mobile phase materials and study them with J-coupled correlation spectroscopies under magic angle spinning conditions. We kept extensive ¹H-spin-echo records that allowed us to examine the effect of magic angle spinning on the transverse relaxation time of water and lipids over time. While it is impossible to distinguish based on chemical shifts, the relaxation data suggest that the signals are more consistent with the interpretation of phospholipid membranes than triglycerides in lipid droplets. In particular, the simultaneous T₂ changes in water and lipids suggest that the centrifugal impact of magic angle spinning alters the lipid's structure in very tight spaces. Additional evidence of phospholipid membranes came from the choline-γ resonance at 3.2 ppm in fresh samples, which disappears with magic angle spinning. Thus, the fatty acid signals are at least partially from membrane bilayer structures, and we propose that they are linked to the submicroscopic vascularization channels similar to the dense canaliculi network in mammalian bones. Our detection of phospholipids from bones depended critically on two factors: (1) the elimination of the overwhelming triglyceride signals from marrows and (2) the preservation of water that biomembranes require. The relaxation data reveal aspects of lipid fluidity that have not been elucidated by previous order parameter studies on model membranes. Relaxation times have long been considered difficult to interpret. A robust and renewed understanding may be beneficial.

The effect of acidic pH on the adsorption and lytic activity of the peptides Polybia-MP1 and its histidine-containing analog in anionic lipid membrane: a biophysical study by molecular dynamics and spectroscopy

Antimicrobial peptides (AMPs) are part of the innate immune system of many species. AMPs are short sequences rich in charged and non-polar residues. They act on the lipid phase of the plasma membrane without requiring membrane receptors. Polybia-MP1 (MP1), extracted from a native wasp, is a broad-spectrum bactericide, an inhibitor of cancer cell proliferation being non-hemolytic and non-cytotoxic. MP1 mechanism of action and its adsorption mode is not yet completely known. Its adsorption to lipid bilayer and lytic activity is most likely dependent on the ionization state of its two acidic and three basic residues and consequently on the bulk pH. Here we investigated the effect of bulk acidic (pH 5.5) and neutral pH (7.4) solution on the adsorption, insertion, and lytic activity of MP1 and its analog H-MP1 to anionic (7POPC:3POPG) model membrane. H-MP1 is a synthetic analog of MP1 with lysines replaced by histidines. Bulk pH changes could modulate this peptide efficiency. The combination of different experimental techniques and molecular dynamics (MD) simulations showed that the adsorption, insertion, and lytic activity of H-MP1 are highly sensitive to bulk pH in opposition to MP1. The atomistic details, provided by MD simulations, showed peptides contact their N-termini to the bilayer before the insertion and then lay parallel to the bilayer. Their hydrophobic faces inserted into the acyl chain phase disturb the lipid-packing.