The Influence of Cholesterol on the Properties and Permeability of Hypericin Derivatives in Lipid Membranes

Increasing the efficiency of hyperthermic intraperitoneal chemotherapy (HIPEC) by combination with a photosensitive drug in pediatric rhabdomyosarcoma in an animal model

BACKGROUND

Cytoreductive surgery (CRS) in combination with hyperthermic intraperitoneal chemotherapy (HIPEC) is an option in advanced peritoneal sarcomatosis. Nevertheless, CRS and HIPEC are not successful in all patients. An enhancement of HIPEC using photodynamic therapy (PDT) might be beneficial. Therefore, a combination of the photosensitizer hypericin (HYP) with HIPEC was evaluated in an animal model.

PROCEDURE

An established HIPEC animal model for rhabdomyosarcoma (NOD/LtSz-scid IL2Rγnullmice, n = 80) was used. All groups received HYP (100 μg/200 μl) intraperitoneally with and without cisplatin-based (30 or 60 mg/m² ) HIPEC (37°C or 42°C, for 60 minutes) (five groups, each n = 16). Peritoneal cancer index (PCI) was documented visually and by HYP-based photodynamic diagnosis (PDD). HYP-based PDT of the tumor was performed. Tissue samples were evaluated regarding proliferation (Ki-67) and apoptosis (TUNEL).

RESULTS

HYP uptake was detected even in smallest tumor nodes (<1 mm) with improved tumor detection during PDD (PCI with PDD vs. PCI without PDD: 8.5 vs. 7, p < .001***). Apoptotic effects after PDT without HIPEC were limited to the tumor surface, whereas PDT after HIPEC (60 mg/m² , 42°C) showed additional reduction of tumor proliferation in the top nine to 11 cell layers (50 μm).

CONCLUSION

HYP as fluorescent photosensitizer offers an intraoperative diagnostic advantage detecting intraperitoneal tumor dissemination. The combination of HYP and cisplatin-based HIPEC was feasible in vivo, showing enhanced effects on tumor proliferation and apoptosis induction across the tumor surface. Further studies combining HYP and HIPEC will follow to establish a clinical application.

Visible light-driven i-motif-based DNAzymes

DNA foldings provide variant possibilities to develop DNAzymes with remarkable catalytic performance. In spite of fruitful reports on G-quadruplex DNAzymes, four-stranded cytosine-rich i-motifs have not been explored as the potential skeletons of DNAzymes. In this work, we developed a visible light-driven DNAzyme based on human telomeric i-motifs using a natural photosensitizer of hypericin (Hyp) as the cofactor and dissolved oxygen as the oxidant source. The i-motif folding in acidic solution caused the distal thymine overhangs at the 3' and 5' ends to approach each other to provide a favorable binding site for Hyp via an interaction of fully complementary hydrogen bonding. However, the i-motifs without the distal overhangs or with the inappropriate overhang length and the base identity exhibited no binding with Hyp. The binding event converted Hyp from the fully dark state to the emissive state under visible light illumination. Subsequently, the excited Hyp had an opportunity to transfer energy to dissolved oxygen. Resultantly, singlet oxygen (¹O₂) was generated to initiate the substrate oxidation. The catalytic performance of the DNAzyme can be improved using a long-lived mediator. Our developed i-motif-based DNAzyme can be driven by almost the whole range of visible lights, suggesting broad applications in the photocatalytic fields, for example, as an alternative strategy in developing biodevices.

Photosensitizing Antivirals

Antiviral action of various photosensitizers is already summarized in several comprehensive reviews, and various mechanisms have been proposed for it. However, a critical consideration of the matter of the area is complicated, since the exact mechanisms are very difficult to explore and clarify, and most publications are of an empirical and "phenomenological" nature, reporting a dependence of the antiviral action on illumination, or a correlation of activity with the photophysical properties of the substances. Of particular interest is substance-assisted photogeneration of highly reactive singlet oxygen (1O2). The damaging action of 1O2 on the lipids of the viral envelope can probably lead to a loss of the ability of the lipid bilayer of enveloped viruses to fuse with the lipid membrane of the host cell. Thus, lipid bilayer-affine 1O2 photosensitizers have prospects as broad-spectrum antivirals against enveloped viruses. In this short review, we want to point out the main types of antiviral photosensitizers with potential affinity to the lipid bilayer and summarize the data on new compounds over the past three years. Further understanding of the data in the field will spur a targeted search for substances with antiviral activity against enveloped viruses among photosensitizers able to bind to the lipid membranes.

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.

Computational and Experimental Approaches to Investigate Lipid Nanoparticles as Drug and Gene Delivery Systems

Lipid nanoparticles (LNPs) have been widely applied in drug and gene delivery. More than twenty years ago, DoxilTM was the first LNPs-based drug approved by the US Food and Drug Administration (FDA). Since then, with decades of research and development, more and more LNP-based therapeutics have been used to treat diverse diseases, which often offer the benefits of reduced toxicity and/or enhanced efficacy compared to the active ingredients alone. Here, we provide a review of recent advances in the development of efficient and robust LNPs for drug/gene delivery. We emphasize the importance of rationally combining experimental and computational approaches, especially those providing multiscale structural and functional information of LNPs, to the design of novel and powerful LNP-based delivery systems.

Mucoadhesive In Situ Gelling Liquid Crystalline Precursor System to Improve the Vaginal Administration of Drugs

The vaginal mucosa is a very promising route for drug administration due to its high permeability and the possibility to bypass first pass metabolism; however, current vaginal dosage forms present low retention times due to their dilution in vaginal fluids, which hampers the efficacy of many pharmacological treatments. In order to overcome these problems, this study proposes to develop a mucoadhesive in situ gelling liquid crystalline precursor system composed of 30% of oleic acid and cholesterol (7:1), 40% of ethoxylated and propoxylated cetyl alcohol, and 30% of a dispersion of 16% Poloxamer 407. The effect of the dilution with simulated vaginal fluid (SVF) on this system was evaluated by polarized light microscopy (PLM), small-angle X-ray scattering (SAXS), rheological studies, texture profile analysis (TPA), mucoadhesion study, in vitro drug release test using hypericin (HYP) as drug model, and cytotoxicity assay. PLM and SAXS confirmed the formation of an isotropic system. After the addition of three different concentrations of SVF (30, 50, and 100%), the resultant formulations presented anisotropy and characteristics of viscous lamellar phases. Rheology shows that formulations with SVF behaved as a non-Newtonian fluid with suitable shear thinning for vaginal application. TPA and mucoadhesion assays indicated the formation of long-range ordered systems as the amount of SVF increases which may assist in the fixation of the formulation on the vaginal mucosa. The formulations were able to control about 75% of the released HYP demonstrating a sustained release profile. Finally, all formulations acted as safe vaginal drug delivery systems.

Probing the interactions of phthalocyanine-based photosensitizers with model phospholipid bilayer by molecular dynamics simulations

Phthalocyanines (Pc) have received considerable attention in the design of photosensitizers for photodynamic therapy (PDT). It is of great interest to design novel Pc-based photosensitizer with improved biological efficiency, which largely relies on the understanding of the interactions of Pc with cell membranes at a molecular level. Here, via all-atom molecular dynamics simulations, we explored the interaction mechanism between a model phospholipid bilayer and three zinc Pc (ZnPc) molecules with different hydrophobicity in nature. We find that the adsorption and insertion behaviors of ZnPc molecules in the model bilayer are different due to the differing hydrophobicity and interaction patterns with phospholipids. Moreover, our simulations demonstrate that the conjunction of a ZnPc skeleton with a cholesterol moiety may reduce the intrinsic molecular rotation of ZnPc in membranes, presumably leading to an increase of the generation efficiency of reactive oxygen species for PDT. The molecular insights obtained here are likely to help improve the rational design of novel photosensitizers with enhanced cellular uptake and photocytotoxic activity.

Rheological Droplet Interface Bilayers (rheo-DIBs): Probing the Unstirred Water Layer Effect on Membrane Permeability via Spinning Disk Induced Shear Stress

A new rheological droplet interface bilayer (rheo-DIB) device is presented as a tool to apply shear stress on biological lipid membranes. Despite their exciting potential for affecting high-throughput membrane translocation studies, permeability assays conducted using DIBs have neglected the effect of the unstirred water layer (UWL). However as demonstrated in this study, neglecting this phenomenon can cause significant underestimates in membrane permeability measurements which in turn limits their ability to predict key processes such as drug translocation rates across lipid membranes. With the use of the rheo-DIB chip, the effective bilayer permeability can be modulated by applying shear stress to the droplet interfaces, inducing flow parallel to the DIB membranes. By analysing the relation between the effective membrane permeability and the applied stress, both the intrinsic membrane permeability and UWL thickness can be determined for the first time using this model membrane approach, thereby unlocking the potential of DIBs for undertaking diffusion assays. The results are also validated with numerical simulations.

Formation of Large Hypericin Aggregates in Giant Unilamellar Vesicles-Experiments and Modeling

The incorporation of hypericin (Hyp) from aqueous solutions into giant unilamellar vesicle (GUV) membranes has been studied experimentally and by means of kinetic Monte Carlo modeling. The time evolution of Hyp fluorescence originating from Hyp monomers dissolved in the GUV membrane has been recorded by confocal microscopy and while trapping individual GUVs in optical tweezers. It was shown that after reaching a maximum, the fluorescence intensity gradually decreased toward longer times. Formation of oversized Hyp clusters has been observed on the GUV surface at prolonged time. A simplified kinetic Monte Carlo model is presented to follow the aggregation/dissociation processes of Hyp molecules in the membrane. The simulation results reproduced the basic experimental observations: the scaling of the characteristic fluorescence decay time with the vesicle diameter and the buildup of large Hyp clusters in the GUV membrane.

Physical insights into the blood-brain barrier translocation mechanisms

The number of individuals suffering from diseases of the central nervous system (CNS) is growing with an aging population. While candidate drugs for many of these diseases are available, most of these pharmaceutical agents cannot reach the brain rendering most of the drug therapies that target the CNS inefficient. The reason is the blood-brain barrier (BBB), a complex and dynamic interface that controls the influx and efflux of substances through a number of different translocation mechanisms. Here, we present these mechanisms providing, also, the necessary background related to the morphology and various characteristics of the BBB. Moreover, we discuss various numerical and simulation approaches used to study the BBB, and possible future directions based on multi-scale methods. We anticipate that this review will motivate multi-disciplinary research on the BBB aiming at the design of effective drug therapies.

Differential diagnosis of gallstones by using hypericin as a fluorescent optical imaging agent

AIM: To explore the feasibility of using hypericin as an optical imaging probe with affinity for cholesterol for differential fluorescent detection of human gallstones.

METHODS: Cholesterol, mixed and pigment stones from cholecystectomy patients were incubated with hypericin or solvent. After 72 h, the stones were analysed for fluorescence (365 nm) and treated with 2-propanol/dimethyl sulfoxide for high performance liquid chromatography (HPLC) analysis. Rats with virtual gallbladder containing human cholesterol, mixed or pigment gallstones (VGHG) received 5 mg/kg hypericin or solvent and VGHG rats with cholesterol stones were given different hypericin doses (5-15 mg/kg). Twelve hours later, the stones were analysed at 365 nm. Biliary excretion and metabolites of hypericin were assessed in common bile duct (CBD) cannulated rats for 9 h using fluorospectrometry, HPLC and matrix-assisted laser desorption/ionization-time-of-flight mass spectrometry (MALDI-TOF MS).

RESULTS: Homogeneous high fluorescence was seen on cholesterol stones either pre-incubated with hypericin or extracted from VGHG rats receiving hypericin. Mixed stones showed a dotted fluorescent pattern, whereas pigment and solvent-treated ones lacked fluorescence. HPLC showed 7.68, 6.65 and 0.08 × 10^-3 M of cholesterol in extracts from cholesterol, mixed, and pigment gallstones, respectively. Hypericin accounted for 2.0, 0.5 and 0.2 × 10^-6 M in that order. On cholesterol stones from VGHG rats receiving different hypericin doses, a positive correlation was observed between dose and fluorescence. In the bile from CBD-cannulated rats, fluorescence represented 20% of the injected dose with two peaks in 9 h. HPLC analysis revealed that hypericin conjugates reached 60% of the peak area. By MALDI-TOF MS, hypericin-glucuronide was detected.

CONCLUSION: This study proves the potential use of hypericin for differential fluorescent detection of human gallstones regarding their chemical composition.

In silico pharmacology: Drug membrane partitioning and crossing

Over the past decade, molecular dynamics (MD) simulations have become particularly powerful to rationalize drug insertion and partitioning in lipid bilayers. MD simulations efficiently support experimental evidences, with a comprehensive understanding of molecular interactions driving insertion and crossing. Prediction of drug partitioning is discussed with respect to drug families (anesthetics; β-blockers; non-steroidal anti-inflammatory drugs; antioxidants; antiviral drugs; antimicrobial peptides). To accurately evaluate passive permeation coefficients turned out to be a complex theoretical challenge; however the recent methodological developments based on biased MD simulations are particularly promising. Particular attention is paid to membrane composition (e.g., presence of cholesterol), which influences drug partitioning and permeation. Recent studies concerning in silico models of membrane proteins involved in drug transport (influx and efflux) are also reported here. These studies have allowed gaining insight in drug efflux by, e.g., ABC transporters at an atomic resolution, explicitly accounting for the mandatory forces induced by the surrounded lipid bilayer. Large-scale conformational changes were thoroughly analyzed.

Molecular Dynamics Studies of Liposomes as Carriers for Photosensitizing Drugs: Development, Validation, and Simulations with a Coarse-Grained Model

Liposomes are proposed as drug delivery systems and can in principle be designed so as to cohere with specific tissue types or local environments. However, little detail is known about the exact mechanisms for drug delivery and the distributions of drug molecules inside the lipid carrier. In the current work, a coarse-grained (CG) liposome model is developed, consisting of over 2500 lipids, with varying degrees of drug loading. For the drug molecule, we chose hypericin, a natural compound proposed for use in photodynamic therapy, for which a CG model was derived and benchmarked against corresponding atomistic membrane bilayer model simulations. Liposomes with 21-84 hypericin molecules were generated and subjected to 10 microsecond simulations. Distribution of the hypericins, their orientations within the lipid bilayer, and the potential of mean force for transferring a hypericin molecule from the interior aqueous "droplet" through the liposome bilayer are reported herein.

Molecular simulation of nonfacilitated membrane permeation

This is a review. Non-electrolytic compounds typically cross cell membranes by passive diffusion. The rate of permeation is dependent on the chemical properties of the solute and the composition of the lipid bilayer membrane. Predicting the permeability coefficient of a solute is important in pharmaceutical chemistry and toxicology. Molecular simulation has proven to be a valuable tool for modeling permeation of solutes through a lipid bilayer. In particular, the solubility-diffusion model has allowed for the quantitative calculation of permeability coefficients. The underlying theory and computational methods used to calculate membrane permeability are reviewed. We also discuss applications of these methods to examine the permeability of solutes and the effect of membrane composition on permeability. The application of coarse grain and polarizable models is discussed. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Sphingolipids contribute to acetic acid resistance in Zygosaccharomyces bailii

Lignocellulosic raw material plays a crucial role in the development of sustainable processes for the production of fuels and chemicals. Weak acids such as acetic acid and formic acid are troublesome inhibitors restricting efficient microbial conversion of the biomass to desired products. To improve our understanding of weak acid inhibition and to identify engineering strategies to reduce acetic acid toxicity, the highly acetic‐acid‐tolerant yeast Zygosaccharomyces bailii was studied. The impact of acetic acid membrane permeability on acetic acid tolerance in Z. bailii was investigated with particular focus on how the previously demonstrated high sphingolipid content in the plasma membrane influences acetic acid tolerance and membrane permeability. Through molecular dynamics simulations, we concluded that membranes with a high content of sphingolipids are thicker and more dense, increasing the free energy barrier for the permeation of acetic acid through the membrane. Z. bailii cultured with the drug myriocin, known to decrease cellular sphingo­lipid levels, exhibited significant growth inhibition in the presence of acetic acid, while growth in medium without acetic acid was unaffected by the myriocin addition. Furthermore, following an acetic acid pulse, the intracellular pH decreased more in myriocin‐treated cells than in control cells. This indicates a higher inflow rate of acetic acid and confirms that the reduction in growth of cells cultured with myriocin in the medium with acetic acid was due to an increase in membrane permeability, thereby demonstrating the importance of a high fraction of sphingolipids in the membrane of Z. bailii to facilitate acetic acid resistance; a property potentially transferable to desired production organisms suffering from weak acid stress. Biotechnol. Bioeng. 2016;113: 744–753. © 2015 The Authors. Biotechnology and Bioengineering Published by Wiley Periodicals, Inc.

A simple and transferable all-atom/coarse-grained hybrid model to study membrane processes

We present an efficient all-atom/coarse-grained hybrid model and apply it to membrane processes. This model is an extension of the all-atom/ELBA model applied previously to processes in water. Here, we improve the efficiency of the model by implementing a multiple-time step integrator that allows the atoms and the coarse-grained beads to be propagated at different timesteps. Furthermore, we fine-tune the interaction between the atoms and the coarse-grained beads by computing the potential of mean force of amino acid side chain analogs along the membrane normal and comparing to atomistic simulations. The model was independently validated on the calculation of small-molecule partition coefficients. Finally, we apply the model to membrane peptides. We studied the tilt angle of the Walp23 and Kalp23 helices in two different model membranes and the stability of the glycophorin A dimer. The model is efficient, accurate, and straightforward to use, as it does not require any extra interaction particles, layers of atomistic solvent molecules or tabulated potentials, thus offering a novel, simple approach to study membrane processes.

Cholesterol, sphingolipids, and glycolipids: what do we know about their role in raft-like membranes?

Lipids rafts are considered to be functional nanoscale membrane domains enriched in cholesterol and sphingolipids, characteristic in particular of the external leaflet of cell membranes. Lipids, together with membrane-associated proteins, are therefore considered to form nanoscale units with potential specific functions. Although the understanding of the structure of rafts in living cells is quite limited, the possible functions of rafts are widely discussed in the literature, highlighting their importance in cellular functions. In this review, we discuss the understanding of rafts that has emerged based on recent atomistic and coarse-grained molecular dynamics simulation studies on the key lipid raft components, which include cholesterol, sphingolipids, glycolipids, and the proteins interacting with these classes of lipids. The simulation results are compared to experiments when possible.

An overview of translational (radio)pharmaceutical research related to certain oncological and non-oncological applications

Translational medicine pursues the conversion of scientific discovery into human health improvement. It aims to establish strategies for diagnosis and treatment of diseases. Cancer treatment is difficult. Radio-pharmaceutical research has played an important role in multiple disciplines, particularly in translational oncology. Based on the natural phenomenon of necrosis avidity, OncoCiDia has emerged as a novel generic approach for treating solid malignancies. Under this systemic dual targeting strategy, a vascular disrupting agent first selectively causes massive tumor necrosis that is followed by iodine-131 labeled-hypericin (¹²³I-Hyp), a necrosis-avid compound that kills the residual cancer cells by crossfire effect of beta radiation. In this review, by emphasizing the potential clinical applicability of OncoCiDia, we summarize our research activities including optimization of radioiodinated hypericin Hyp preparations and recent studies on the biodistribution, dosimetry, pharmacokinetic and, chemical and radiochemical toxicities of the preparations. Myocardial infarction is a global health problem. Although cardiac scintigraphy using radioactive perfusion tracers is used in the assessment of myocardial viability, searching for diagnostic imaging agents with authentic necrosis avidity is pursued. Therefore, a comparative study on the biological profiles of the necrosis avid ¹²³I-Hyp and the commercially available ^99mTc-Sestamibi was conducted and the results are demonstrated. Cholelithiasis or gallstone disease may cause gallbladder inflammation, infection and other severe complications. While studying the mechanisms underlying the necrosis avidity of Hyp and derivatives, their naturally occurring fluorophore property was exploited for targeting cholesterol as a main component of gallstones. The usefulness of Hyp as an optical imaging agent for cholelithiasis was studied and the results are presented. Multiple uses of automatic contrast injectors may reduce costs and save resources. However, cross-contaminations with blood-borne pathogens of infectious diseases may occur. We developed a radioactive method for safety evaluation of a new replaceable patient-delivery system. By mimicking pathogens with a radiotracer, we assessed the feasibility of using the system repeatedly without septic risks. This overview is deemed to be interesting to those involved in the related fields for translational research.

Effects of Phospholipid Composition on the Transfer of a Small Cationic Peptide Across a Model Biological Membrane

The transfer of a lysine amino acid analogue across phospholipid membrane models was investigated using molecular-dynamics simulations. The evolution of the protonation state of this small peptide as a function of its position inside the membrane was studied by determining the local pKa by means of free-energy calculations. Permeability and mean-first-passage time were evaluated and showed that the transfer occurs on the submillisecond time scale. Comparative studies were conducted to evaluate changes in the pKa arising from differences in the phospholipid chemical structure. We compared, hence, the effect of an ether vs an ester linkage of the lipid headgroup as well as linear vs branched lipid tails. The study reveals that protonated lysine residues can be buried further inside an ether lipid membrane than an ester lipid membrane, while branched lipids are found to stabilize less the charged form compared to their unbranched lipid chain counterparts.

Improved clearance of radioiodinated hypericin as a targeted anticancer agent by using a duodenal drainage catheter in rats

We sought to reduce the radioactive intestinal waste after intravenous injection of necrosis avid iodine-131-labeled hypericin in dual-targeting anticancer radiotherapy and to study its pharmacokinetics in rats using a newly designed catheter. Iodine-123-labeled hypericin was prepared with iodogen as oxidant and characterized by high-performance liquid chromatography and mass spectrometry. After iodine-123-labeled hypericin administration, duodenal juice was collected via a catheter from groups of rats (n = 5) at intervals of 0-4, 4-8 or 20-24 h. The content was assessed by gamma-counting. The biodistribution and pharmacokinetics of iodine-123-labeled hypericin were investigated in rats without (n = 5) and with continuous catheterization (n = 5) for 9 h. After labeling, a high radiochemical yield was obtained with iodine-123-labeled hypericin (>95%), as confirmed by high-performance liquid chromatography and mass spectrometry. In the duodenal aspirate from animals with intermittent catheterization during 24 h, radioactivity accounted for 46% of the total with two peaks at 3 h and 8 h, suggesting enterohepatic circulation. Rats with 9 h of catheterization exhibited one peak representing 20% of the radioactivity. Major metabolites appeared to be conjugated iodine-123-labeled hypericin forms. In rats without and with catheter, iodine-123-labeled hypericin showed exponential elimination from plasma with no significant dehalogenation. Delayed iodine-123-labeled hypericin excretion, a higher maximum concentration (Cmax), larger area under concentration-time curve [AUC(0-∞)] and a longer mean residence time were observed in non-catheterized animals (P < 0.05). The catheterized group exhibited lower urinary excretion than non-catheterized group (P < 0.05). Rats with a catheter showed lower radioactivity (P = 0.01) in the small intestines than those without a catheter (1.82 ± 0.41 versus 18.95 ± 4.32 percentage of the injected dose). After iodine-123-labeled hypericin administration, the radioactivity excreted into bile was efficiently removed from the body via a duodenal catheter. Radiation overexposure due to the prolonged elimination of iodine-131-labeled hypericin can be prevented using this approach.

Quantitative analysis of molecular transport across liposomal bilayer by J-mediated 13C Overhauser dynamic nuclear polarization

We introduce a new NMR technique to dramatically enhance the solution-state (13)C NMR sensitivity and contrast at 0.35 T and at room temperature by actively transferring the spin polarization from Overhauser dynamic nuclear polarization (ODNP)-enhanced (1)H to (13)C nuclei through scalar (J) coupling, a method that we term J-mediated (13)C ODNP. We demonstrate the capability of this technique by quantifying the permeability of glycine across negatively charged liposomal bilayers composed of dipalmitoylphosphatidylcholine (DPPC) and dipalmitoylphosphatidylglycerol (DPPG). The permeability coefficient of glycine across this DPPC/DPPG bilayer is measured to be (1.8 ± 0.1) × 10(-11)m/s, in agreement with the literature value. We further observed that the presence of 20 mol % cholesterol within the DPPC/DPPG lipid membrane significantly retards the permeability of glycine by a factor of 4. These findings demonstrate that the high sensitivity and contrast of J-mediated (13)C ODNP affords the measurement of the permeation kinetics of small hydrophilic molecules across lipid bilayers, a quantity that is difficult to accurately measure with existing techniques.

Another Piece of the Membrane Puzzle: Extending Slipids Further

To be able to model complex biological membranes in a more realistic manner, the force field Slipids (Stockholm lipids) has been extended to include parameters for sphingomyelin (SM), phosphatidylglycerol (PG), phosphatidylserine (PS) lipids, and cholesterol. Since the parametrization scheme was faithful to the scheme used in previous editions of Slipids, all parameters are consistent and fully compatible. The results of careful validation of a number of key structural properties for one and two component lipid bilayers are in excellent agreement with experiments. Potentials of mean force for transferring water across binary mixtures of lipids and cholesterol were also computed in order to compare water permeability rates to experiments. In agreement with experimental and simulation studies, it was found that the permeability and partitioning of water is affected by cholesterol in lipid bilayers made of saturated lipids to the largest extent. With the extensions of Slipids presented here, it is now possible to study complex systems containing many different lipids and proteins in a fully atomistic resolution in the isothermic-isobaric (NPT) ensemble, which is the proper ensemble for membrane simulations.

Identifying the sarco(endo)plasmic reticulum Ca2+ ATPase (SERCA) as a potential target for hypericin--a theoretical study

The exact cellular target for the potent anti-cancer agent hypericin has not yet been determined; this thus encourages the application of computational chemistry tools to be employed in order to provide insights that can be employed in further drug development studies. In the present study computational docking and molecular dynamics simulations are applied to investigate possible interactions between hypericin and the Ca(2+) pump SERCA as proposed in the literature. Hypericin was found to bind strongly both in pockets within the transmembrane region and in the cytosolic region of the protein, although the two studied isoforms of SERCA differ slightly in their preferred binding sites. The calculated binding energies for hypericin in the four investigated sites were of the same magnitude as for thapsigargin (TG), the most potent SERCA inhibitor, or in the range between TG and di-tert-butylhydroquinone (BHQ), which is also known to possess inhibitory activity. The hydrophobic character of hypericin indicates that the molecule initially binds in the ER membrane from which it diffuses into the transmembrane region of the protein and to binding pockets therein. The transmembrane TG and BHQ binding pockets provide suitable locations for hypericin as they allow for favourable interactions with the lipid tails that surround these. High binding energies were noted for hypericin in these pockets and are expected to constitute highly possible binding sites due to their accessibility from the ER membrane. Hypericin most likely binds to both isoforms of SERCA and acts as an inhibitor or, under light irradiation, as a singlet oxygen generator that in turn degrades the protein or induces lipid peroxidation.

Testing physical models of passive membrane permeation

The biophysical basis of passive membrane permeability is well-understood, but most methods for predicting membrane permeability in the context of drug design are based on statistical relationships that indirectly capture the key physical aspects. Here, we investigate molecular mechanics-based models of passive membrane permeability and evaluate their performance against different types of experimental data, including parallel artificial membrane permeability assays (PAMPA), cell-based assays, in vivo measurements, and other in silico predictions. The experimental data sets we use in these tests are diverse, including peptidomimetics, congeneric series, and diverse FDA approved drugs. The physical models are not specifically trained for any of these data sets; rather, input parameters are based on standard molecular mechanics force fields, such as partial charges, and an implicit solvent model. A systematic approach is taken to analyze the contribution from each component in the physics-based permeability model. A primary factor in determining rates of passive membrane permeation is the conformation-dependent free energy of desolvating the molecule, and this measure alone provides good agreement with experimental permeability measurements in many cases. Other factors that improve agreement with experimental data include deionization and estimates of entropy losses of the ligand and the membrane, which lead to size-dependence of the permeation rate.

Large influence of cholesterol on solute partitioning into lipid membranes

Cholesterol plays an important role in maintaining the correct fluidity and rigidity of the plasma membrane of all animal cells, and hence, it is present in concentrations ranging from 20 to 50 mol %. Whereas the effect of cholesterol on such mechanical properties has been studied exhaustively over the last decades, the structural basis for cholesterol effects on membrane permeability is still unclear. Here we apply systematic molecular dynamics simulations to study the partitioning of solutes between water and membranes. We derive potentials of mean force for six different solutes permeating across 20 different lipid membranes containing one out of four types of phospholipids plus a cholesterol content varying from 0 to 50 mol %. Surprisingly, cholesterol decreases solute partitioning into the lipid tail region of the membranes much more strongly than expected from experiments on macroscopic membranes, suggesting that a laterally inhomogeneous cholesterol concentration and permeability may be required to explain experimental findings. The simulations indicate that the cost of breaking van der Waals interactions between the lipid tails of cholesterol-containing membranes account for the reduced partitioning rather than the surface area per phospholipid, which has been frequently suggested as a determinant for solute partitioning. The simulations further show that the partitioning is more sensitive to cholesterol (i) for larger solutes, (ii) in membranes with saturated as compared to membranes with unsaturated lipid tails, and (iii) in membranes with smaller lipid head groups.

Convergence of Free Energy Profile of Coumarin in Lipid Bilayer

Atomistic molecular dynamics (MD) simulations of druglike molecules embedded in lipid bilayers are of considerable interest as models for drug penetration and positioning in biological membranes. Here we analyze partitioning of coumarin in dioleoylphosphatidylcholine (DOPC) bilayer, based on both multiple, unbiased 3 μs MD simulations (total length) and free energy profiles along the bilayer normal calculated by biased MD simulations (∼7 μs in total). The convergences in time of free energy profiles calculated by both umbrella sampling and z-constraint techniques are thoroughly analyzed. Two sets of starting structures are also considered, one from unbiased MD simulation and the other from "pulling" coumarin along the bilayer normal. The structures obtained by pulling simulation contain water defects on the lipid bilayer surface, while those acquired from unbiased simulation have no membrane defects. The free energy profiles converge more rapidly when starting frames from unbiased simulations are used. In addition, z-constraint simulation leads to more rapid convergence than umbrella sampling, due to quicker relaxation of membrane defects. Furthermore, we show that the choice of RESP, PRODRG, or Mulliken charges considerably affects the resulting free energy profile of our model drug along the bilayer normal. We recommend using z-constraint biased MD simulations based on starting geometries acquired from unbiased MD simulations for efficient calculation of convergent free energy profiles of druglike molecules along bilayer normals. The calculation of free energy profile should start with an unbiased simulation, though the polar molecules might need a slow pulling afterward. Results obtained with the recommended simulation protocol agree well with available experimental data for two coumarin derivatives.