Breakage of Hydrophobic Contacts Limits the Rate of Passive Lipid Exchange between Membranes

Improving the Computational Efficiency of the Adaptive Biasing Force Sampling by Leveraging the Telescopic-Solvation Scheme

The number of solvent molecules present in the system during molecular dynamics is the balancing act between the need to remove the boundary effects present in the system and the computational cost. Application of the telescopic-solvation box scheme during the estimation of the potential of mean force (PMF) can be advantageous in situations where the contribution of solvent far from the site of interest toward the whole PMF is negligible, as previously demonstrated in the case of adaptive steered molecular dynamics and umbrella sampling. This work explores the application of the telescopic-solvation box scheme during enhanced sampling by the stratified adaptive biasing force (ABF) family of methods, including ABF, extended ABF, well-tempered-metadynamics extended ABF, and multiwalker extended ABF. During this scheme, the number of water molecules differed in each stratified window, whose number depended on the value of the collective variable being sampled in that window. Two systems were used to verify the viability of the telescopic scheme: unfolding (Ala)10 peptide in water and insertion of α-tocopherol in a bilayer membrane. In the first system, the 1D and 2D PMFs obtained by the telescopic-solvation scheme matched well with the benchmark PMFs estimated with a standard solvation box. The minimal energy path connecting the α-helical and extended conformational states revealed that the unfolding process of (Ala)10 in water involved multiple closely spaced metastable states. As for the second system, the PMF, equilibrium location of α-tocopherol, and the free energy associated with the desorption and flipping of α-tocopherol obtained within the scope of the telescopic-solvation box scheme agreed with their standard solvation box values. Enhanced sampling with ABF and its variants in conjunction with the telescopic-solvation scheme results in a similar quality of the estimated PMF compared to sampling with a standard solvation box, albeit with reduced computational cost.

Elucidating Daptomycin's Antibacterial Efficacy: Insights into the Tripartite Complex with Lipid II and Phospholipids in Bacterial Septum Membrane

This study elucidated the mechanism of formation of a tripartite complex containing daptomycin (Dap), lipid II, and phospholipid phosphatidylglycerol in the bacterial septum membrane, which was previously reported as the cause of the antibacterial action of Dap against gram-positive bacteria via molecular dynamics and enhanced sampling methods. Others have suggested that this transient complex ushers in the inhibition of cell wall synthesis by obstructing the downstream polymerization and cross-linking processes involving lipid II, which is absent in the presence of cardiolipin lipid in the membrane. In this work, we observed that the complex was stabilized by Ca2+-mediated electrostatic interactions between Dap and lipid head groups, hydrophobic interaction, hydrogen bonds, and salt bridges between the lipopeptide and lipids and was associated with Dap concentration-dependent membrane depolarization, thinning of the bilayer, and increased lipid tail disorder. Residues Orn6 and Kyn13, along with the DXDG motif, made simultaneous contact with constituent lipids, hence playing a crucial role in the formation of the complex. Incorporating cardiolipin into the membrane model led to its competitively displacing lipid II away from the Dap, reducing the lifetime of the complex and the nonexistence of lipid tail disorder and membrane depolarization. No evidence of water permeation inside the membrane hydrophobic interior was noted in all of the systems studied. Additionally, it was shown that using hydrophobic contacts between Dap and lipids as collective variables for enhanced sampling gave rise to a free energy barrier for the translocation of the lipopeptide. A better understanding of Dap's antibacterial mechanism, as studied through this work, will help develop lipopeptide-based antibiotics for rising Dap-resistant bacteria.

Ligand Lipophilicity Determines Molecular Mechanisms of Nanoparticle Adsorption to Lipid Bilayers

The interactions of ligand-functionalized nanoparticles with the cell membrane affect cellular uptake, cytotoxicity, and related behaviors, but relating these interactions to ligand properties remains challenging. In this work, we perform coarse-grained molecular dynamics simulations to study how the adsorption of ligand-functionalized cationic gold nanoparticles (NPs) to a single-component lipid bilayer (as a model cell membrane) is influenced by ligand end group lipophilicity. A set of 2 nm diameter NPs, each coated with a monolayer of organic ligands that differ only in their end groups, was simulated to mimic NPs recently studied experimentally. Metadynamics calculations were performed to determine key features of the free energy landscape for adsorption as a function of the distance of the NP from the bilayer and the number of NP-lipid contacts. These simulations revealed that NP adsorption is thermodynamically favorable for all NPs due to the extraction of lipids from the bilayer and into the NP monolayer. To resolve ligand-dependent differences in adsorption behavior, string method calculations were performed to compute minimum free energy pathways for adsorption. These calculations revealed a surprising nonmonotonic dependence of the free energy barrier for adsorption on ligand end group lipophilicity. Large free energy barriers are predicted for the least lipophilic end groups because favorable NP-lipid contacts are initiated only through the unfavorable protrusion of lipid tail groups out of the bilayer. The smallest free energy barriers are predicted for end groups of intermediate lipophilicity which promote NP-lipid contacts by intercalating within the bilayer. Unexpectedly, large free energy barriers are also predicted for the most lipophilic end groups which remain sequestered within the ligand monolayer rather than intercalating within the bilayer. These trends are broadly in agreement with past experimental measurements and reveal how subtle variations in ligand lipophilicity dictate adsorption mechanisms and associated kinetics by influencing the interplay of lipid-ligand interactions.

Passive permeability controls synthesis for the allelochemical sorgoleone in sorghum root exudate

Competition for soil nutrients and water with other plants foster competition within the biosphere for access to these limited resources. The roots for the common grain sorghum produce multiple small molecules that are released via root exudates into the soil to compete with other plants. Sorgoleone is one such compound, which suppresses weed growth near sorghum by acting as a quinone analog and interferes with photosynthesis. Since sorghum also grows photosynthetically, and may be susceptible to sorgoleone action if present in tissues above ground, it is essential to exude sorgoleone efficiently. However, since the P450 enzymes that synthesize sorgoleone are intracellular, the release mechanism for sorgoleone remain unclear. In this study, we conducted an in silico assessment for sorgoleone and its precursors to passively permeate biological membranes. To facilitate accurate simulation, CHARMM parameters were newly optimized for sorgoleone and its precursors. These parameters were used to conduct 1 μs of unbiased molecular dynamics simulations to compare the permeability of sorgoleone with its precursors molecules. We find that interleaflet transfer is maximized for sorgoleone, suggesting that the precursor molecules may remain in the same leaflet for access by biosynthetic P450 enzymes. Since no sorgoleone was extracted during unbiased simulations, we compute a permeability coefficient using the inhomogeneous solubility diffusion model. The requisite free energy and diffusivity profiles for sorgoleone through a sorghum membrane model were determined through Replica Exchange Umbrella Sampling (REUS) simulations. The REUS calculations highlight that any soluble sorgoleone would quickly insert into a lipid bilayer, and would readily transit. When sorgoleone forms aggregates in root exudate as indicated by our equilibrium simulations, aggregate formation would lower the effective concentration in aqueous solution, creating a concentration gradient that would facilitate passive transport. This suggests that sorgoleone synthesis occurs within sorghum root cells and that sorgoleone is exuded by permeating through the cell membrane without the need for a transport protein once the extracellular sorgoleone aggregate is formed.

Telescoping-Solvation-Box Protocol-Based Umbrella Sampling Method for Potential Mean Force Estimation

Previously, it was shown that the telescoping box scheme, in combination with adaptive steered molecule dynamics (ASMD), can be used to estimate the potential of mean force (PMF) with a decrease in computational cost associated with large solvation boxes. Since ASMD reduces to umbrella sampling (US) in the limit of infinitely slow pulling velocity, a hypothesis was made that the telescoping box scheme can be extended to include the US method. The hypothesis was tested using the unfolding pathway of a polyalanine peptide in a water box and translocation of α-tocopherol through the human membrane. Two different approaches were tried: telescoping US (TELUS), in which the number of solvent molecules was linearly coupled to the reaction coordinate, and Block-TELUS, which was a compromise between the fixed solvation box of the US and the window-wise variable solvation box of TELUS. In the case of polyalanine peptide in a water box, both approaches gave overlapping potential of mean force (PMF) concerning the benchmark US-PMF. Window-wise comparison of properties like root-mean-square inner product, Ramachandran plot, α-helix content, and hydrogen bond formation was used to verify that both approaches captured the same dynamics as the US method. Applying the Block-TELUS protocol in the system with diffusing α-tocopherol through the bilayer resulted in overlapping PMF to its US benchmark. A comparison between the window-wise orientation of the chromanol headgroup also produced almost identical results. This study concluded that with the careful application of telescoping solvation boxes, a less computationally expensive US could be performed for systems where the effect of distant solvent molecules on the configurational space sampled in the window depends on the value of the reaction coordinate.

Lipophilic Anchors that Embed Bioconjugates in Bilayer Membranes: A Review

A wide range of biomaterials and engineered cell surfaces are composed of bioconjugates embedded in liposome membranes, surface-immobilized bilayers, or the plasma membranes of living cells. This review article summarizes the various ways that Nature anchors integral and peripheral proteins in a cell membrane and describes the strategies devised by chemical biologists to label a membrane protein in living cells. Also discussed are modern synthetic and semisynthetic methods to produce lipidated proteins. Subsequent sections describe methods to anchor a three-component synthetic construct that is composed of a lipophilic membrane anchor, hydrophilic linker, and exposed functional component. The surface exposed payload can be a fluorophore, aptamer, oligonucleotide, polypeptide, peptide nucleic acid, polysaccharide, branched dendrimer, or linear polymer. Hydrocarbon chains are commonly used as the membrane anchor, and a general experimental trend is that a two chain lipid anchor has higher membrane affinity than a cholesteryl or single chain lipid anchor. Amphiphilic fluorescent dyes are effective molecular probes for cell membrane imaging and a zwitterionic linker between the fluorophore and the lipid anchor promotes high persistence in the plasma membrane of living cells. A relatively new advance is the development of switchable membrane anchors as molecular tools for fundamental studies or as technology platforms for applied biomaterials.

Ceramide-1-phosphate transfer protein enhances lipid transport by disrupting hydrophobic lipid-membrane contacts

Cellular distributions of the sphingolipid ceramide-1-phosphate (C1P) impact essential biological processes. C1P levels are spatiotemporally regulated by ceramide-1-phosphate transfer protein (CPTP), which efficiently shuttles C1P between organelle membranes. Yet, how CPTP rapidly extracts and inserts C1P into a membrane remains unknown. Here, we devise a multiscale simulation approach to elucidate biophysical details of CPTP-mediated C1P transport. We find that CPTP binds a membrane poised to extract and insert C1P and that membrane binding promotes conformational changes in CPTP that facilitate C1P uptake and release. By significantly disrupting a lipid's local hydrophobic environment in the membrane, CPTP lowers the activation free energy barrier for passive C1P desorption and enhances C1P extraction from the membrane. Upon uptake of C1P, further conformational changes may aid membrane unbinding in a manner reminiscent of the electrostatic switching mechanism used by other lipid transfer proteins. Insertion of C1P into an acceptor membrane, eased by a decrease in membrane order by CPTP, restarts the transfer cycle. Most notably, we provide molecular evidence for CPTP's ability to catalyze C1P extraction by breaking hydrophobic C1P-membrane contacts with compensatory hydrophobic lipid-protein contacts. Our work, thus, provides biophysical insights into how CPTP efficiently traffics C1P between membranes to maintain sphingolipid homeostasis and, additionally, presents a simulation method aptly suited for uncovering the catalytic mechanisms of other lipid transfer proteins.

Mechanistic Insights into Passive Membrane Permeability of Drug-like Molecules from a Weighted Ensemble of Trajectories

Passive permeability of a drug-like molecule is a critical property assayed early in a drug discovery campaign that informs a medicinal chemist how well a compound can traverse biological membranes, such as gastrointestinal epithelial or restrictive organ barriers, so it can perform a specific therapeutic function. However, the challenge that remains is the development of a method, experimental or computational, which can both determine the permeation rate and provide mechanistic insights into the transport process to help with the rational design of any given molecule. Typically, one of the following three methods are used to measure the membrane permeability: (1) experimental permeation assays acting on either artificial or natural membranes; (2) quantitative structure–permeability relationship models that rely on experimental values of permeability or related pharmacokinetic properties of a range of molecules to infer those for new molecules; and (3) estimation of permeability from the Smoluchowski equation, where free energy and diffusion profiles along the membrane normal are taken as input from large-scale molecular dynamics simulations. While all these methods provide estimates of permeation coefficients, they provide very little information for guiding rational drug design. In this study, we employ a highly parallelizable weighted ensemble (WE) path sampling strategy, empowered by cloud computing techniques, to generate unbiased permeation pathways and permeability coefficients for a set of drug-like molecules across a neat 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine membrane bilayer. Our WE method predicts permeability coefficients that compare well to experimental values from an MDCK-LE cell line and PAMPA assays for a set of drug-like amines of varying size, shape, and flexibility. Our method also yields a series of continuous permeation pathways weighted and ranked by their associated probabilities. Taken together, the ensemble of reactive permeation pathways, along with the estimate of the permeability coefficient, provides a clearer picture of the microscopic underpinnings of small-molecule membrane permeation.

Explaining reaction coordinates of alanine dipeptide isomerization obtained from deep neural networks using Explainable Artificial Intelligence (XAI)

A method for obtaining appropriate reaction coordinates is required to identify transition states distinguishing product and reactant in complex molecular systems. Recently, abundant research has been devoted to obtaining reaction coordinates using artificial neural networks from deep learning literature, where many collective variables are typically utilized in the input layer. However, it is difficult to explain the details of which collective variables contribute to the predicted reaction coordinates owing to the complexity of the nonlinear functions in deep neural networks. To overcome this limitation, we used Explainable Artificial Intelligence (XAI) methods of the Local Interpretable Model-agnostic Explanation (LIME) and the game theory-based framework known as Shapley Additive exPlanations (SHAP). We demonstrated that XAI enables us to obtain the degree of contribution of each collective variable to reaction coordinates that is determined by nonlinear regressions with deep learning for the committor of the alanine dipeptide isomerization in vacuum. In particular, both LIME and SHAP provide important features to the predicted reaction coordinates, which are characterized by appropriate dihedral angles consistent with those previously reported from the committor test analysis. The present study offers an AI-aided framework to explain the appropriate reaction coordinates, which acquires considerable significance when the number of degrees of freedom increases.

Membrane hydrophobicity determines the activation free energy of passive lipid transport

The collective behavior of lipids with diverse chemical and physical features determines a membrane's thermodynamic properties. Yet, the influence of lipid physicochemical properties on lipid dynamics, in particular interbilayer transport, remains underexplored. Here, we systematically investigate how the activation free energy of passive lipid transport depends on lipid chemistry and membrane phase. Through all-atom molecular dynamics simulations of 11 chemically distinct glycerophospholipids, we determine how lipid acyl chain length, unsaturation, and headgroup influence the free energy barriers for two elementary steps of lipid transport: lipid desorption, which is rate limiting, and lipid insertion into a membrane. Consistent with previous experimental measurements, we find that lipids with longer, saturated acyl chains have increased activation free energies compared to lipids with shorter, unsaturated chains. Lipids with different headgroups exhibit a range of activation free energies; however, no clear trend based solely on chemical structure can be identified, mirroring difficulties in the interpretation of previous experimental results. Compared to liquid-crystalline phase membranes, gel phase membranes exhibit substantially increased free energy barriers. Overall, we find that the activation free energy depends on a lipid's local hydrophobic environment in a membrane and that the free energy barrier for lipid insertion depends on a membrane's interfacial hydrophobicity. Both of these properties can be altered through changes in lipid acyl chain length, lipid headgroup, and membrane phase. Thus, the rate of lipid transport can be tuned through subtle changes in local membrane composition and order, suggesting an unappreciated role for nanoscale membrane domains in regulating cellular lipid dynamics.

Lipophilicity of Cationic Ligands Promotes Irreversible Adsorption of Nanoparticles to Lipid Bilayers

A mechanistic understanding of the influence of the surface properties of engineered nanomaterials on their interactions with cells is essential for designing materials for applications such as bioimaging and drug delivery as well as for assessing nanomaterial safety. Ligand-coated gold nanoparticles have been widely investigated because their highly tunable surface properties enable investigations into the effect of ligand functionalization on interactions with biological systems. Lipophilic ligands have been linked to adverse biological outcomes through membrane disruption, but the relationship between ligand lipophilicity and membrane interactions is not well understood. Here, we use a library of cationic ligands coated on 2 nm gold nanoparticles to probe the impact of ligand end group lipophilicity on interactions with supported phosphatidylcholine lipid bilayers as a model for cytoplasmic membranes. Nanoparticle adsorption to and desorption from the model membranes were investigated by quartz crystal microbalance with dissipation monitoring. We find that nanoparticle adsorption to model membranes increases with ligand lipophilicity. The effects of ligand structure on gold nanoparticle attachment were further analyzed using atomistic molecular dynamics simulations, which showed that the increase in ligand lipophilicity promotes ligand intercalation into the lipid bilayer. Together, the experimental and simulation results could be described by a two-state model that accounts for the initial attachment and subsequent conversion to a quasi-irreversibly bound state. We find that only nanoparticles coated with the most lipophilic ligands in our nanoparticle library undergo conversion to the quasi-irreversible state. We propose that the initial attachment is governed by interaction between the ligands and phospholipid tail groups, whereas conversion into the quasi-irreversibly bound state reflects ligand intercalation between phospholipid tail groups and eventual lipid extraction from the bilayer. The systematic variation of ligand lipophilicity enabled us to demonstrate that the lipophilicity of cationic ligands correlates with nanoparticle-bilayer adsorption and suggested that changing the nonpolar ligand R group promotes a mechanism of ligand intercalation into the bilayer associated with irreversible adsorption.

Partitioning of Aβ Peptide Fragments into Blood-Brain Barrier Mimetic Bilayer

We used all-atom replica-exchange umbrella sampling molecular dynamics simulations to investigate the partitioning of the charged tetrapeptide KLVF and its neutral apolar counterpart VVIA into the blood-brain barrier (BBB)-mimetic bilayer. Our findings allowed us to reconstruct the partitioning mechanism for these two Aβ peptide fragments. Despite dissimilar sequences, their permeation shares significant common features. Computations of free energies and permeabilities show that partitioning of both peptides is highly unfavorable, ruling out passive transport. The peptides experience multiple rotational transitions within the bilayer and typically cause considerable lipid disorder and bilayer thinning. Near the bilayer midplane, they lose almost entirely their solvation shell and the interactions with the lipid headgroups. The peptides cause complex reorganization within the proximal bilayer region. Upon insertion, they induce striking cholesterol influx reversed by its depletion and the influx of DMPC when the peptides reach the midplane. The differences in partitioning mechanisms are due to the much higher polarity of KLVF peptide, the permeation of which is more unfavorable and which exclusively assumes vertical orientations within the bilayer. In contrast, VVIA positions itself flat between the leaflets, causing minor disorder and even thickening of the BBB-mimetic bilayer. Due to the high density of the cholesterol-rich BBB bilayer, the unfavorable work associated with the peptide insertion provides a significant, but not dominant, contribution to the partition free energy, which is still governed by dehydration and loss of peptide-headgroup interactions. Comparison with experiments indicates that KLVF and VVIA permeation is similar to that of proline tetrapeptide, mannitol, or cimetidine, all of which exhibit no passive transport.