Diffusion of 1,2-Dimethoxyethane and 1,2-Dimethoxypropane through Phosphatidycholine Bilayers: A Molecular Dynamics Study

Mechanism of the enhancing effect of glycyrrhizin on nifedipine penetration through a lipid membrane

The saponin glycyrrhizin from liquorice root shows the ability to enhance the therapeutic activity of other drugs when used as a drug delivery system. Due to its amphiphilic properties, glycyrrhizin can form self-associates (dimers, micelles) and supramolecular complexes with a wide range of hydrophobic drugs, which leads to an increase in their solubility, stability and bioavailability. That is why the mechanism of the biological activity of glycyrrhizin is of considerable interest and has been the subject of intensive physical and chemical research in the last decade. Two mechanisms have been proposed to explain the effect of glycyrrhizin on drug bioavailability, namely, the increase in drug solubility in water and enhancement of the membrane permeability. Interest in the membrane-modifying ability of glycyrrhizic acid (GA) is also growing at present due to its recently discovered antiviral activity against SARS-CoV-2 Bailly and Vergoten (2020) [1]. In the present study, the passive permeability of the DOPC lipid membrane for the calcium channel blocker nifedipine was elucidated by parallel artificial membrane permeability assay (PAMPA) and full atomistic molecular dynamics (MD) simulation with free energy calculations. PAMPA experiments show a remarkable increase in the amount of nifedipine (NF) permeated with glycyrrhizin compared to free NF. In previous studies, we have shown using MD techniques that glycyrrhizin molecules can integrate into the lipid bilayer. In this study, MD simulation demonstrates a significant decrease in the energy barrier of NF penetration through the lipid bilayer in the presence of glycyrrhizin both in the pure DOPC membrane and in the membrane with cholesterol. This effect can be explained by the formation of hydrogen bonds between NF and GA in the middle of the bilayer.

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

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

Glycyrrhizin-Assisted Transport of Praziquantel Anthelmintic Drug through the Lipid Membrane: An Experiment and MD Simulation

Praziquantel (PZQ) is one of the most widespread anthelmintic drugs. However, the frequent insufficient application of PZQ after oral administration is associated with its low solubility, penetration rate, and bioavailability. In the present study, the permeation of PZQ through a 1,2-dioleoyl- sn-glycero-3-phosphocholine (DOPC) membrane was investigated to probe glycyrrhizin-assisted transport. Glycyrrhizin (or glycyrrhizic acid, GA), a natural saponin, shows the ability to enhance the therapeutic activity of various drugs when it is used as a drug delivery system. However, the molecular mechanism of this effect is still under debate. In the present study, the transport rate was measured experimentally by a parallel artificial membrane permeation assay (PAMPA) and molecular dynamics (MD) simulation with DOPC lipid bilayers. The formation of the noncovalent supramolecular complex of PZQ with disodium salt of GA (Na2GA) in an aqueous solution was proved by the NMR relaxation technique. PAMPA experiments show a strong increase in the amount of the penetrating praziquantel molecules in comparison with a saturated aqueous solution of pure drug used as a control. MD simulation of PZQ penetration through the bilayer demonstrates an increase in permeability into the membrane in the presence of a glycyrrhizin molecule. A decrease in the free energy barrier in the middle of the lipid bilayer was obtained, associated with the hydrogen bond between PZQ and GA. Also, GA reduces the local bilayer surface resistance to penetration of PZQ by rearranging the surface lipid headgroups. This study clarifies the mechanism of increasing the drug's bioavailability in the presence of glycyrrhizin.

In situ examination of a charged amino acid-induced structural change in lipid bilayers by sum frequency generation vibrational spectroscopy

The interactions between amino acids (AAs) and membranes represent various short-range and long-range interactions for biological phenomena; however, they are still poorly understood. In this study, we used cationic lysine and arginine as AA models, and systematically investigated the interactions between charged AAs and lipid bilayers using sum frequency generation vibrational spectroscopy (SFG-VS) in situ and in real time. The AA-induced dynamic structural changes of the lipid bilayer were experimentally monitored using the spectral features of CD₂, CD₃, the lipid head phosphate, and carbonyl groups in real time. Time-dependent SFG changes in the structure of the lipid bilayer provide direct evidence for the different interactions of lysine and arginine with the membrane. It was found that the discrepancy between lysine and arginine in binding with the lipid bilayer is due to the nature of the terminal functional groups. Arginine exhibits a more drastic impact on the membrane than lysine. SFG responses of the acyl chains, phosphate groups, and carbonyl groups provide evidence that the interaction between AAs and the membrane most likely follows an electrostatics and hydrogen bond-induced defect model. This work presents an exemplary method for comprehensive investigations of interactions between membranes and other functionally significant substances.

Molecular view of the interaction of S-methyl methanethiosulfonate with DPPC bilayer

We present molecular dynamics (MD) simulation studies of the interaction of a chemo preventive and protective agent, S-methyl methanethiosulfonate (MMTS), with a model bilayer of 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC). We analyzed and compared its diffusion mechanisms with the related molecule dimethyl sulfoxide (DMSO). We obtained spatially resolved free energy profiles of MMTS partition into a DPPC bilayer in the liquid-crystalline phase through potential of mean force (PMF) calculations using an umbrella sampling technique. These profiles showed a minimum for MMTS close to the carbonyl region of DPPC. The location of MMTS molecules in the DPPC bilayer observed in the MD was confirmed by previous SERS studies. We decomposed PMF profiles into entropic and enthalpic contributions. These results showed that the driving force for the partitioning of MMTS into the upper region of DPPC is driven by a favorable entropy change while partitioning into the acyl chains is driven by enthalpy. On the other hand, the partition of DMSO into the membrane is not favored, and is driven by entropy instead of enthalpy. Free diffusion MD simulations using all atom and coarse grained (CG) models of DPPC in presence of MMTS were used to analyze the effect of DPPC-MMTS interaction. Density profiles showed that MMTS locates preferentially in the carbonyl region, as expected according to the PMF profile and the experimental evidence. MMTS presented two differential effects over the packing of DPPC hydrocarbonate chains at low or at high molar ratios. An ordering effect was observed when a CG MMTS model was used. Finally, free diffusion MD and PMF decomposition for DMSO were used for comparison.

Understanding the interaction of block copolymers with DMPC lipid bilayer using coarse-grained molecular dynamics simulations

In this paper, we present a computational model of the adsorption and percolation mechanism of poloxamers (poly(ethylene oxide) (PEO) and poly(propylene oxide) (PPO) triblock copolymers) across a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid bilayer. A coarse-grained model was used to cope with the long time scale of the percolation process. The simulations have provided details of the interaction mechanism of Pluronics with lipid bilayer. In particular, the results have shown that polymer chains containing a PPO block with a length comparable to the DMPC bilayer thickness, such as P85, tends to percolate across the lipid bilayer. On the contrary, Pluronics with a shorter PPO chain, such as L64 and F38, insert partially into the membrane with the PPO block part while the PEO blocks remain in water on one side of the lipid bilayer. The percolation of the polymers into the lipid tail groups reduces the membrane thickness and increases the area per lipid. These effects are more evident for P85 than L64 or F38. Our findings are qualitatively in good agreement with published small-angle X-ray scattering experiments that have evidenced a thinning effect of Pluronics on the lipid bilayer as well as the role of the length of the PPO block on the permeation process of the polymer through the lipid bilayer. Our theoretical results complement the experimental data with a detailed structural and dynamic model of poloxamers at the interface and inside the lipid bilayer.

Behavior of 4-hydroxynonenal in phospholipid membranes

Under conditions of oxidative stress, 4-hydroxy-2-nonenal (4-HNE) is commonly present in vivo. This highly reactive and cytotoxic compound is generated by oxidation of lipids in membranes and can be easily transferred from a membrane to both cytosol and the extracellular space. Employing time-dependent fluorescence shift (TDFS) method and molecular dynamics simulations, we found that 4-HNE is stabilized in the carbonyl region of a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. 4-HNE is thus able to react with cell membrane proteins and lipids. Stabilization in the membrane is, however, moderate and a transfer of 4-HNE to either extra- or intracellular space occurs on a microsecond time scale. These molecular-level details of 4-HNE behavior in the lipid membrane rationalize the experimentally observed reactivity of 4-HNE with proteins inside and outside the cell. Furthermore, these results support the view that 4-HNE may play an active role in cell signaling pathways.