Predicting doxorubicin drug delivery by single-walled carbon nanotube through cell membrane in the absence and presence of nicotine molecules: a molecular dynamics simulation study

Abstract

In order to study the interaction of the anticancer agent Doxorubicin with the single-walled carbon nanotubes with different diameters as drug delivery systems, the molecular dynamics (MD) simulations have been used. Also, for design and development of intracellular Doxorubicin drug delivery systems, a series of steered MD simulations are applied to explore the possibility of encapsulated Doxorubicin–carbon nanotube penetration through a lipid bilayer in presence and absence of Nicotine molecules at different pulling rates. Our simulation results showed that in spite of the adsorption of drug molecules on the outer sidewall of the nanotubes, the spontaneous localization of one Doxorubicin molecule into the cavity of the nanovectors with larger diameters is observed. It is found that the presence of Nicotine molecules in extracellular medium increases the required force for pulling nanotube-encapsulated drug as well as the required time for penetration process, especially at higher velocity. Also, the entering process of the Nicotine molecules into the carbon nanotube causes that the encapsulated drug molecule is fully released in the hydrophobic phase of the lipid bilayer.

Communicated by Ramaswamy H. Sarma     

Structural Effects and Translocation of Doxorubicin in a DPPC/Chol Bilayer: The Role of Cholesterol

We use molecular dynamics simulations to characterize the influence of cholesterol (Chol) on the interaction between the anticancer drug doxorubicin (DOX) and a dipalmitoyl phosphatidylcholine/Chol lipid bilayer. We calculate the potential of mean force, which gives us an estimate of the free energy barrier for DOX translocation across the membrane. We find free energy barriers of 23.1 ± 3.1 k_BT, 36.8 ± 5.1 k_BT, and 54.5 ± 4.7 k_BT for systems composed of 0%, 15%, and 30% Chol, respectively. Our predictions agree with Arrhenius activation energies from experiments using phospholipid membranes, including 20 k_BT for 0% Chol and 37.2 k_BT for 20% Chol. The location of the free energy barrier for translocation across the bilayer is dependent on composition. As Chol concentration increases, this barrier changes from the release of DOX into the water to flip-flop over the membrane center. The drug greatly affects local membrane structure by attracting dipalmitoyl phosphatidylcholine headgroups, curving the membrane, and allowing water penetration. Despite its hydrophobicity, DOX facilitates water transport via its polar groups.     

The interactions of bilirubin with model and biological membranes

The partitioning of bilirubin between albumin and model and biological membranes and the differential partitioning of bilirubin between membranes with different lipid and protein compositions were measured. Partition coefficients were independent of the concentration of bilirubin in membranes up to at least 7 mol of bilirubin/mol of phospholipid. The avidity of albumin for bilirubin was greater than that of membranes, but the avidity of the latter for bilirubin depended on the composition of the membrane. Bilirubin partitioned preferentially into model membranes comprised of microsomal lipids greater than dioleoylphosphatidylcholine = plasma membrane lipids much greater than egg phosphatidylcholine = dimyristoylphosphatidylcholine. Partitioning into membranes was increased if these contained proteins, but the effect of proteins could not be attributed to specific binding to sites on proteins, as reflected by the temperature independence of partition coefficients. Differential partitioning of bilirubin into different membranes of pure lipids also was independent of temperature. Differences in the bulk phase fluidity of membranes does not appear to account for the preferential partitioning of bilirubin into some membranes. It appears that bilirubin partitions into elements of free volume of differing sizes in membranes with variable lipid compositions and that the size of these elements can be increased by adding proteins to membranes.     

Combined effects of avasimibe immunotherapy,doxorubicin chemotherapy,and metal–organic frameworks nanoparticles on breast cancer

CD8⁺ T cells play a vital role in cancer immunotherapy and can be shaped by metabolism. Avasimibe is an acyl coenzyme A-cholesterol acyltransferase (ACAT) inhibitor, which has been clinically verified safe in other phase Ⅲ clinical trials. It can potentiate the killing function of CD8⁺ T cells by modulating cholesterol metabolism. Doxorubicin (DOX) is an anticancer drug widely used in many cancers to induce tumor cell apoptosis. Unfortunately, DOX also can induce toxic and side effects in many organs, compromising its usage and efficacy. Herein, we report the combinational usage of avasimibe and a safe pH sensitive nano-drug delivery system composing of DOX and metal–organic frameworks nanoparticles (MNPs). Our findings demonstrated that DOX–MNPs treatment inhibited tumor growth with good safety profile and avasimibe treatment combined DOX–MNPs treatment exhibited a better efficacy than monotherapies in 4T1 breast cancer therapy.     

Membrane Partitioning: ��Classical�� and ��Nonclassical�� Hydrophobic Effects

The free energy of transfer of nonpolar solutes from water to lipid bilayers is often dominated by a large negative enthalpy rather than the large positive entropy expected from the hydrophobic effect. This common observation has led to the idea that membrane partitioning is driven by the "nonclassical" hydrophobic effect. We examined this phenomenon by characterizing the partitioning of the well-studied peptide melittin using isothermal titration calorimetry (ITC) and circular dichroism (CD). We studied the temperature dependence of the entropic (-TΔS) and enthalpic (ΔH) components of free energy (ΔG) of partitioning of melittin into lipid membranes made of various mixtures of zwitterionic and anionic lipids. We found significant variations of the entropic and enthalpic components with temperature, lipid composition and vesicle size but only small changes in ΔG (entropy-enthalpy compensation). The heat capacity associated with partitioning had a large negative value of about -0.5 kcal mol(-1) K(-1). This hallmark of the hydrophobic effect was found to be independent of lipid composition. The measured heat capacity values were used to calculate the hydrophobic-effect free energy ΔG (hΦ), which we found to dominate melittin partitioning regardless of lipid composition. In the case of anionic membranes, additional free energy comes from coulombic attraction, which is characterized by a small effective peptide charge due to the lack of additivity of hydrophobic and electrostatic interactions in membrane interfaces [Ladokhin and White J Mol Biol 309:543-552, 2001]. Our results suggest that there is no need for a special effect-the nonclassical hydrophobic effect-to describe partitioning into lipid bilayers.     

Determining Peptide Partitioning Properties via Computer Simulation

The transfer of polypeptide segments into lipid bilayers to form transmembrane helices represents the crucial first step in cellular membrane protein folding and assembly. This process is driven by complex and poorly understood atomic interactions of peptides with the lipid bilayer environment. The lack of suitable experimental techniques that can resolve these processes both at atomic resolution and nanosecond timescales has spurred the development of computational techniques. In this review, we summarize the significant progress achieved in the last few years in elucidating the partitioning of peptides into lipid bilayer membranes using atomic detail molecular dynamics simulations. Indeed, partitioning simulations can now provide a wealth of structural and dynamic information. Furthermore, we show that peptide-induced bilayer distortions, insertion pathways, transfer free energies, and kinetic insertion barriers are now accurate enough to complement experiments. Further advances in simulation methods and force field parameter accuracy promise to turn molecular dynamics simulations into a powerful tool for investigating a wide range of membrane active peptide phenomena.     

Improved therapeutic efficiency of photothermal treatment and nursing care in prostate cancer by DOX loaded PEG coated Cu@Se nano-hybrid vesicle

Cancer is one of the major life threatening diseases, with higher mortality rate and morbidity. It is always a challenge for effective drug delivery and release of drug in specific tumor sites. Therefore to identify the synergistic effect of chemotherapeutic drug and photo thermal agent on tumor area, Doxorubicin (DOX) acts as anticancer agent but it has low aqueous solubility so its clinical application is limited. The present study developed doxorubicin (DOX) were designed to be with the poly ethylene glycol (PEG) functionalized copper and selenium (Cu-Se) nanoparticles (PEG@Cu-Se+DOX) and it is efficiently synthesized and enhance its aqueous formulation and improve the prostate cancer (DU145 and LNCaP) activity. The characteristics like mono dispersity, size stability and constant spectral of as-synthesized nanoparticles are comparably excellent than DOX alone. Also the enhanced cellular uptake and in vitro cytoxicicty suggests these nanoparticles selectively killing prostate cancer. In this present study explained that PEG@Cu-Se+DOX as a safe and hopeful strategy for chemotherapeutics of photothermal therapy and deserve for further clinical evaluations.     

Protein Partitioning into Ordered Membrane Domains: Insights from Simulations

Cellular membranes are laterally organized into domains of distinct structures and compositions by the differential interaction affinities between various membrane lipids and proteins. A prominent example of such structures are lipid rafts, which are ordered, tightly packed domains that have been widely implicated in cellular processes. The functionality of raft domains is driven by their selective recruitment of specific membrane proteins to regulate their interactions and functions; however, there have been few general insights into the factors that determine the partitioning of membrane proteins between coexisting liquid domains. In this work, we used extensive coarse-grained and atomistic molecular dynamics simulations, potential of mean force calculations, and conceptual models to describe the partitioning dynamics and energetics of a model transmembrane domain from the linker of activation of T cells. We find that partitioning between domains is determined by an interplay between protein-lipid interactions and differential lipid packing between raft and nonraft domains. Specifically, we show that partitioning into ordered domains is promoted by preferential interactions between peptides and ordered lipids, mediated in large part by modification of the peptides by saturated fatty acids (i.e., palmitoylation). Ordered phase affinity is also promoted by elastic effects, specifically hydrophobic matching between the membrane and the peptide. Conversely, ordered domain partitioning is disfavored by the tight molecular packing of the lipids therein. The balance of these dominant drivers determines partitioning. In the case of the wild-type linker of activation of T cells transmembrane domain, these factors combine to yield enrichment of the peptide at L_o/L_d interfaces. These results define some of the general principles governing protein partitioning between coexisting membrane domains and potentially explain previous disparities among experiments and simulations across model systems.     

Interaction of Doxorubicin and Dipalmitoylphosphatidylcholine Liposomes

The interaction between doxorubicin (DOX), an anthracycline antibiotic frequently used in chemotherapy, and zwitterionic dipalmitoylphosphatidylcholine (DPPC) was investigated using Fourier transform infrared (FTIR) spectroscopy, differential scanning calorimetry (DSC), and rheological measurements. FTIR results showed that DOX shifted the wavenumber of the PO₂ ⁻ band for pure DPPC to a higher wavenumber. This may have been because of the strong interactions between the NH₃ ⁺ group in DOX and the phosphate (PO₂ ⁻) group in the polar head of DPPC. The main transition temperature of DPPC liposomes was slightly shifted to a lower temperature for DPPC liposome-encapsulated DOX. This suggested that DOX had a significant effect on the acyl chains in the DPPC bilayers, and that its presence decreased the transition cooperativity of lipid acyl chains. There was also the appearance of an additional transition peak at nearly 136°C for the DPPC/DOX sample. These interactions between DOX and DPPC phospholipid would cause a decrease in the DPPC liposomes plastic viscosity and increase membrane fluidity. A better understanding of the interactions between DOX and lipid bilayers could help in the design and development of improved liposomal drug delivery systems.     

Lipid Clustering Correlates with Membrane Curvature as Revealed by Molecular Simulations of Complex Lipid Bilayers

Cell membranes are complex multicomponent systems, which are highly heterogeneous in the lipid distribution and composition. To date, most molecular simulations have focussed on relatively simple lipid compositions, helping to inform our understanding of in vitro experimental studies. Here we describe on simulations of complex asymmetric plasma membrane model, which contains seven different lipids species including the glycolipid GM3 in the outer leaflet and the anionic lipid, phosphatidylinositol 4,5-bisphophate (PIP₂), in the inner leaflet. Plasma membrane models consisting of 1500 lipids and resembling the in vivo composition were constructed and simulations were run for 5 µs. In these simulations the most striking feature was the formation of nano-clusters of GM3 within the outer leaflet. In simulations of protein interactions within a plasma membrane model, GM3, PIP₂, and cholesterol all formed favorable interactions with the model α-helical protein. A larger scale simulation of a model plasma membrane containing 6000 lipid molecules revealed correlations between curvature of the bilayer surface and clustering of lipid molecules. In particular, the concave (when viewed from the extracellular side) regions of the bilayer surface were locally enriched in GM3. In summary, these simulations explore the nanoscale dynamics of model bilayers which mimic the in vivo lipid composition of mammalian plasma membranes, revealing emergent nanoscale membrane organization which may be coupled both to fluctuations in local membrane geometry and to interactions with proteins.     

Partitioning of 2,6-Bis(1H-Benzimidazol-2-yl)pyridine fluorophore into a phospholipid bilayer: complementary use of fluorescence quenching studies and molecular dynamics simulations

Successful use of fluorescence sensing in elucidating the biophysical properties of lipid membranes requires knowledge of the distribution and location of an emitting molecule in the bilayer. We report here that 2,6-bis(1H-benzimidazol-2-yl)pyridine (BBP), which is almost non-fluorescent in aqueous solutions, reveals a strong emission enhancement in a hydrophobic environment of a phospholipid bilayer, making it interesting for fluorescence probing of water content in a lipid membrane. Comparing the fluorescence behavior of BBP in a wide variety of solvents with those in phospholipid vesicles, we suggest that the hydrogen bonding interactions between a BBP fluorophore and water molecules play a crucial role in the observed “light switch effect”. Therefore, the loss of water-induced fluorescence quenching inside a membrane are thought to be due to deep penetration of BBP into the hydrophobic, water-free region of a bilayer. Characterized by strong quenching by transition metal ions in solution, BBP also demonstrated significant shielding from the action of the quencher in the presence of phospholipid vesicles. We used the increase in fluorescence intensity, measured upon titration of probe molecules with lipid vesicles, to estimate the partition constant and the Gibbs free energy (ΔG) of transfer of BBP from aqueous buffer into a membrane. Partitioning BBP revealed strongly favorable ΔG, which depends only slightly on the lipid composition of a bilayer, varying in a range from − 6.5 to − 7.0 kcal/mol. To elucidate the binding interactions of the probe with a membrane on the molecular level, a distribution and favorable location of BBP in a POPC bilayer were modeled via atomistic molecular dynamics (MD) simulations using two different approaches: (i) free, diffusion-driven partitioning of the probe molecules into a bilayer and (ii) constrained umbrella sampling of a penetration profile of the dye molecule across a bilayer. Both of these MD approaches agreed with regard to the preferred location of a BBP fluorophore within the interfacial region of a bilayer, located between the hydrocarbon acyl tails and the initial portion of the lipid headgroups. MD simulations also revealed restricted permeability of water molecules into this region of a POPC bilayer, determining the strong fluorescence enhancement observed experimentally for the membrane-partitioned form of BBP.     

Radio Frequency-Activated Nanoliposomes for Controlled Combination Drug Delivery

This work was conducted in order to design, characterize, and evaluate stable liposomes containing the hydrophobic drug raloxifene HCl (RAL) and hydrophilic doxycycline HCl (DOX), two potentially synergistic agents for treating osteoporosis and other bone lesions, in conjunction with a radio frequency-induced, hydrophobic magnetic nanoparticle-dependent triggering mechanism for drug release. Both drugs were successfully incorporated into liposomes by lipid film hydration, although combination drug loading compromised liposome stability. Liposome stability was improved by reducing the drug load and by including Pluronics® (PL) in the formulations. DOX did not appear to interact with the phospholipid membranes comprising the liposomes, and its release was maximized in the presence of radio frequency (RF) heating. In contrast, differential scanning calorimetry (DSC) and phosphorus-31 nuclear magnetic resonance (³¹P-NMR) analysis revealed that RAL developed strong interactions with the phospholipid membranes, most notably with lipid phosphate head groups, resulting in significant changes in membrane thermodynamics. Likewise, RAL release from liposomes was minimal, even in the presence of RF heating. These studies may offer useful insights into the design and optimization of multidrug containing liposomes. The effects of RAL on liposome characteristics and drug release performance underscore the importance of appropriate physical-chemical analysis in order to identify and characterize drug-lipid interactions that may profoundly affect liposome properties and performance early in the formulation development process.KEY WORDS: controlled release, drug combination, liposomes, nanoparticles     

Protective role of CoQ10 or L-carnitine on the integrity of the myocardium in doxorubicin induced toxicity

Doxorubicin (DOX) is a chemotherapeutic agent used for treatment of different cancers and its clinical usage is hindered by the oxidative injury-related cardiotoxicity. This work aims to declare if the harmful effects of DOX on heart can be alleviated with the use of Coenzyme Q10 (CoQ10) or L-carnitine. The study was performed on seventy two female Wistar albino rats divided into six groups, 12 animals each: Control group; DOX group (10 mg/kg); CoQ10 group (200 mg/kg); L-carnitine group (100 mg/kg); DOX + CoQ10 group; DOX + L-carnitine group. CoQ10 and L-carnitine treatment orally started 5 days before a single dose of 10 mg/kg DOX that injected intraperitoneally (IP) then the treatment continued for 10 days. At the end of the study, serum biochemical parameters of cardiac damage, oxidative stress indices, and histopathological changes were investigated. CoQ10 or L-carnitine showed a noticeable effects in improving cardiac functions evidenced reducing serum enzymes as serum interleukin-1 beta (IL-1 β), tumor necrosis factor alpha (TNF-α), leptin, lactate dehydrogenase (LDH), Cardiotrophin-1, Troponin-I and Troponin-T. Also, alleviate oxidative stress, decrease of cardiac Malondialdehyde (MDA), Nitric oxide (NO) and restoring cardiac reduced glutathione levels to normal levels. Both corrected the cardiac alterations histologically and ultrastructurally. With a visible improvements in α-SMA, vimentin and eNOS immunohistochemical markers. CoQ10 or L-carnitine supplementation improves the functional and structural integrity of the myocardium.     

Preclinical Studies with Doxorubicin Encapsulated in Polyethyleneglycol-Coated Liposomes

Abstract

Doxorubicin (DOX) has been encapsulated with high efficiency in the water phase of small-sized lipid vesicles. Plasma-induced drug leakage from these vesicles is minimal when hydrogenated phosphatidylcholine is present as the main component. A prolonged circulation time of liposome-encapsulated DOX is observed in animal models when a small fraction of polyethyleneglycol-derivatized phospholipid (PEG) is present in the liposome bilayer. Using these PEG-coated liposomes, we found that the concentration of DOX in tumor implants of the mouse M-109 carcinoma is significantly enhanced by liposome delivery. The antitumor activity of liposome-encapsulated DOX in a lung metastases model of the M-109 carcinoma is superior to that of free DOX. The minimal lethal dose of DOX to tumor-free mice was substantially increased by encapsulation in PEG-coated liposomes, indicating that toxicity is reduced. We also found that the vesicant of DOX after intradermal injection is prevented by liposome encapsulation. These preclinical observations, suggesting that encapsulation of DOX in PEG-coated liposomes may lead to a significant improvement of the therapeutic index of DOX, have led to the initiation of clinical trials in cancer patients.     

Effects of curvature and composition on α-synuclein binding to lipid vesicles

Parkinson's disease is characterized by the presence of intracellular aggregates composed primarily of the neuronal protein α-synuclein (αS). Interactions between αS and various cellular membranes are thought to be important to its native function as well as relevant to its role in disease. We use fluorescence correlation spectroscopy to investigate binding of αS to lipid vesicles as a function of the lipid composition and membrane curvature. We determine how these parameters affect the molar partition coefficient of αS, providing a quantitative measure of the binding energy, and calculate the number of lipids required to bind a single protein. Specific anionic lipids have a large effect on the free energy of binding. Lipid chain saturation influences the binding interaction to a lesser extent, with larger partition coefficients measured for gel-phase vesicles than for fluid-phase vesicles, even in the absence of anionic lipid components. Although we observe variability in the binding of the mutant proteins, differences in the free energies of partitioning are less dramatic than with varied lipid compositions. Vesicle curvature has a strong effect on the binding affinity, with a >15-fold increase in affinity for small unilamellar vesicles over large unilamellar vesicles, suggesting that αS may be a curvature-sensing protein. Our findings provide insight into how physical properties of the membrane may modulate interactions of αS with cellular membranes.     

Molecular dynamics simulation of cation–phospholipid clustering in phospholipid bilayers: Possible role in stalk formation during membrane fusion

In this study, we performed all-atom long-timescale molecular dynamics simulations of phospholipid bilayers incorporating three different proportions of negatively charged lipids in the presence of K⁺, Mg^2 +, and Ca^2 + ions to systemically determine how membrane properties are affected by cations and lipid compositions. Our simulations revealed that the binding affinity of Ca^2 + ions with lipids is significantly stronger than that of K⁺ and Mg^2 + ions, regardless of the composition of the lipid bilayer. The binding of Ca^2 + ions to the lipids resulted in bilayers having smaller lateral areas, greater thicknesses, greater order, and slower rotation of their lipid head groups, relative to those of corresponding K⁺- and Mg^2 +-containing systems. The Ca^2 + ions bind preferentially to the phosphate groups of the lipids. The complexes formed between the cations and the lipids further assembled to form various multiple-cation-centered clusters in the presence of anionic lipids and at higher ionic strength—most notably for Ca^2 +. The formation of cation–lipid complexes and clusters dehydrated and neutralized the anionic lipids, creating a more-hydrophobic environment suitable for membrane aggregation. We propose that the formation of Ca^2 +–phospholipid clusters across apposed lipid bilayers can work as a “cation glue” to adhere apposed membranes together, providing an adequate configuration for stalk formation during membrane fusion.     

Influence of lipid rafts on CD1d presentation by dendritic cells

Our main objective was to analyze the role of lipid rafts in the activation of Vα-14^? and Vα-14⁺ T hybridomas by dendritic cells. We showed that activation of Vα-14⁺ hybridomas by dendritic cells or other CD1d-expressing cells was altered by disruption of lipid rafts with the cholesterol chelator MβCD. However, CD1d presentation to autoreactive Vα-14^? anti-CD1d hybridomas which do not require the endocytic pathway was not altered. Using partitioning of membrane fractions with Brij98 at 37°C, we confirmed that CD1d was enriched in subcellular fractions corresponding to lipid rafts and we describe that α-GalCer enhanced CD1d amount in the low density detergent insoluble fraction. We conclude that the membrane environment of CD1d can influence antigen presentation mainly when the endocytic pathway is required. Flow cytometry analysis can provide additional information on lipid rafts in plasma membranes and allows a dynamics follow-up of lipid rafts partitioning. Using this method, we showed that CD1d plasma membrane expression was sensitive to low concentrations of detergent. This may suggest either that CD1d is associated with lipid rafts mainly in intracellular membranes or that its association with the lipid rafts in the plasma membrane is weak.     

The influence of cholesterol on interactions and dynamics of ibuprofen in a lipid bilayer

In this work, molecular dynamics (MD) simulations with atomistic details were performed to examine the influence of the cholesterol on the interactions and the partitioning of the hydrophobic drug ibuprofen in a fully hydrated 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) bilayer. Analysis of MD simulations indicated that ibuprofen molecules prefer to be located in the hydrophobic acyl chain region of DMPC/cholesterol bilayers. This distribution decreases the lateral motion of lipid molecules. The presence of ibuprofen molecules in the bilayers with 0 and 25 mol% cholesterol increases the ordering of hydrocarbon tails of lipids whereas for the bilayers with 50 mol% cholesterol, ibuprofen molecules perturb the flexible chains of DMPC lipids which leads to the reduction of the acyl chain order parameter. The potential of the mean force (PMF) method was used to calculate the free energy profile for the transferring of an ibuprofen molecule from the bulk water into the DMPC/cholesterol membranes. The PMF studies indicated that the presence of 50 mol% cholesterol in the bilayers increases the free energy barrier and slows down the permeation of the ibuprofen drug across the DMPC bilayer. This can be due to the condensing and ordering effects of the cholesterol on the bilayer.     

Alternative Mechanisms for the Interaction of the Cell-Penetrating Peptides Penetratin and the TAT Peptide with Lipid Bilayers

Cell-penetrating peptides (CPPs) have recently attracted much interest due to their apparent ability to penetrate cell membranes in an energy-independent manner. Here molecular-dynamics simulation techniques were used to study the interaction of two CPPs: penetratin and the TAT peptide with 1,2-Dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) phospolipid bilayers shed light on alternative mechanisms by which these peptides might cross biological membranes. In contrast to previous simulation studies of charged peptides interacting with lipid bilayers, no spontaneous formation of transmembrane pores was observed. Instead, the simulations suggest that the peptides may enter the cell by micropinocytosis, whereby the peptides induce curvature in the membrane, ultimately leading to the formation of small vesicles within the cell that encapsulate the peptides. Specifically, multiple peptides were observed to induce large deformations in the lipid bilayer that persisted throughout the timescale of the simulations (hundreds of nanoseconds). Pore formation could be induced in simulations in which an external potential was used to pull a single penetratin or TAT peptide into the membrane. With the use of umbrella-sampling techniques, the free energy of inserting a single penetratin peptide into a DPPC bilayer was estimated to be ∼75 kJmol⁻¹, which suggests that the spontaneous penetration of single peptides would require a timescale of at least seconds to minutes. This work also illustrates the extent to which the results of such simulations can depend on the initial conditions, the extent of equilibration, the size of the system, and the conditions under which the simulations are performed. The implications of this with respect to the current systems and to simulations of membrane-peptide interactions in general are discussed.     

Interactions between model membranes and lignin-related compounds studied by immobilized liposome chromatography

In order to elucidate the modes of interaction between lignin precursors and membranes, we have studied the influence of temperature, lipid composition and buffer composition on the partitioning of monolignol and dilignol model substances into phospholipid bilayers. The partitioning was determined by immobilized liposome chromatography, which is an established method for studies of pharmaceutical drugs but a new approach in studies of lignin synthesis. The temperature dependence of the retention and the effect of a high ammonium sulfate concentration in the mobile phase demonstrated that the interaction involved both hydrophobic effects and polar interactions. There was also a good correlation between the partitioning and the estimated hydrophobicity, in terms of octanol/water partitioning. The partitioning behavior of the model substances suggests that passive diffusion over the cell membrane is a possible transport route for lignin precursors. This conclusion is strengthened by comparison of the present results with the partitioning of pharmaceutical drugs that are known to pass cell membranes by diffusion.