Fluorogenic Chemical Probes for Wash-free Imaging of Cell Membrane Damage in Ferroptosis, Necrosis, and Axon Injury.

Selectively labeling cells with damaged membranes is needed not only for identifying dead cells in culture, but also for imaging membrane barrier dysfunction in pathologies in vivo. Most membrane permeability stains are permanently colored or fluorescent dyes that need washing to remove their non-uptaken extracellular background and reach good image contrast. Others are DNA-binding environment-dependent fluorophores, which lack design modularity, have potential toxicity, and can only detect permeabilization of cell volumes containing a nucleus (i.e., cannot delineate damaged volumes in vivo nor image non-nucleated cell types or compartments). Here, we develop modular fluorogenic probes that reveal the whole cytosolic volume of damaged cells, with near-zero background fluorescence so that no washing is needed. We identify a specific disulfonated fluorogenic probe type that only enters cells with damaged membranes, then is enzymatically activated and marks them. The esterase probe MDG1 is a reliable tool to reveal live cells that have been permeabilized by biological, biochemical, or physical membrane damage, and it can be used in multicolor microscopy. We confirm the modularity of this approach by also adapting it for improved hydrolytic stability, as the redox probe MDG2. We conclude by showing the unique performance of MDG probes in revealing axonal membrane damage (which DNA fluorogens cannot achieve) and in discriminating damage on a cell-by-cell basis in embryos in vivo. The MDG design thus provides powerful modular tools for wash-free in vivo imaging of membrane damage, and indicates how designs may be adapted for selective delivery of drug cargoes to these damaged cells: offering an outlook from selective diagnosis toward therapy of membrane-compromised cells in disease.

MemCross: Accelerated Weight Histogram method to assess membrane permeability.

Passive permeation events across biological membranes are determining steps in the pharmacokinetics of xenobiotics. To reach an accurate and rapid prediction of membrane permeation coefficients of drugs is a complex challenge, which can efficiently support drug discovery. Such predictions are indeed highly valuable as they may guide the selection of potential leads with optimum bioavailabilities prior to synthesis. Theoretical models exist to predict these coefficients. Many of them are based on molecular dynamics (MD) simulations, which allow calculation of permeation coefficients through the evaluation of both the potential of mean force (PMF) and the diffusivity profiles. However, these simulations still require intensive computational efforts, and novel methodologies should be developed and benchmarked. Free energy perturbation (FEP) method was recently shown to estimate PMF with a significantly reduced computational cost compared to the adaptive biasing force method. This benchmarking was achieved with small molecules, namely short-chain alcohols. Here, we show that to estimate the PMF of bulkier, drug-like xenobiotics, conformational sampling is a critical issue. To reach a sufficient sampling with FEP calculations requires a relatively long time-scale, which can lower the benefits related to the computational gain. In the present work, the Accelerated Weight Histogram (AWH) method was employed for the first time in all-atom membrane models. The AWH-based protocol, named MemCross, appears affordable to estimate PMF profiles of a series of drug-like xenobiotics, compared to other enhanced sampling methods. The continuous exploration of the crossing pathway by MemCross also allows modeling subdiffusion by computing fractional diffusivity profiles. The method is also versatile as its input parameters are largely insensitive to the molecule properties. It also ensures a detailed description of the molecule orientations along the permeation pathway, picturing all intermolecular interactions at an atomic resolution. Here, MemCross was applied on a series of 12 xenobiotics, including four weak acids, and a coherent structure-activity relationship was established.

MRP4 is responsible for the efflux transport of mycophenolic acid ß-d glucuronide (MPAG) from hepatocytes to blood.

Mycophenolic acid (MPA) has become a cornerstone of immunosuppressive therapy, in particular for transplant patients. In the gastrointestinal tract, the liver and the kidney, MPA is mainly metabolized into phenyl-ß-d glucuronide (MPAG). Knowledge about the interactions between MPA/MPAG and membrane transporters is still fragmented. The aim of the present study was to explore these interactions with the basolateral hepatic MRP4 transporter. The inhibition of the MRP4-driven transport by various drugs which can be concomitantly prescribed was also evaluated. In vitro experiments using vesicles overexpressing MRP4 showed an ATP-dependent transport of MPAG driven by MRP4 (Michaelis-Menten constant of 233.9 ± 32.8 µM). MPA was not effluxed by MRP4. MRP4-mediated transport of MPAG was inhibited (from -43% to -84%) by ibuprofen, cefazolin, cefotaxime and micafungin. An in silico approach based on molecular docking and molecular dynamics simulations rationalized the mode of binding of MPAG to MRP4. The presence of the glucuronide moiety in MPAG was highlighted as key, being prone to make electrostatic and H-bond interactions with specific residues of the MRP4 protein chamber. This explains why MPAG is a substrate of MRP4 whereas MPA is not.

Acetylated lignin nanoparticles as a possible vehicle for photosensitizing molecules.

Lignins are underused and abundant bio-sourced polymers with various potential applications. An attractive one is the development of nanoparticles for bioactive compound delivery. Here, we optimized the synthesis of hydrodispersible nanoparticles of acetylated lignin by comparing different lignin sources, degrees of acetylation and preparation methods. The formation of acetylated lignin nanoparticles in various solvents was probed by both experiments and, for the first time, a molecular dynamics simulation. We showed that dialysis is more suitable to obtain these nanoparticles than anti-solvent addition. The encapsulation of hydrophobic photosensitizing porphyrin in these nanoparticles was also demonstrated and rationalized at the molecular level, together with experiments, docking and molecular dynamics simulations. As acetylated lignin has been demonstrated to exhibit photosensitizing activity, the encapsulation of bioactive compounds in lignin nanoparticles opens the doors to a broad range of potential applications.

Lignin Modification Supported by DFT-Based Theoretical Study as a Way to Produce Competitive Natural Antioxidants.

The valorization of lignins as renewable aromatic feedstock is of utmost importance in terms of the use of sustainable resources. This study provides a deductive approach towards market-oriented lignin-derived antioxidants by ascertaining the direct effect of different structural features of lignin on the reactivity of its phenolic OH groups in the radical scavenging reactions. The antioxidant activity of a series of compounds, modeling lignin structural units, was experimentally characterized and rationalized, using thermodynamic descriptors. The calculated O-H bond dissociation enthalpies (BDE) of characteristic lignin subunits were used to predict the modification pathways of technical lignins. The last ones were isolated by soda delignification from different biomass sources and their oligomeric fractions were studied as a raw material for modification and production of optimized antioxidants. These were characterized in terms of chemical structure, molecular weight distribution, content of the functional groups, and the antioxidant activity. The developed approach for the targeted modification of lignins allowed the products competitive with two commercial synthetic phenolic antioxidants in both free radical scavenging and stabilization of thermooxidative destruction of polyurethane films.

Lipid bilayer position and orientation of novel carprofens, modulators of γ-secretase in Alzheimer's disease.

γ-Secretase is an integral membrane protein complex and is involved in the cleavage of the amyloid precursor protein APP to produce amyloid-ß peptides. Amyloid-ß peptides are considered causative agents for Alzheimer's disease and drugs targeted at γ-secretase are investigated as therapeutic treatments. We synthesized new carprofen derivatives, which showed γ-secretase modulating activity and determined their precise position, orientation, and dynamics in lipid membranes by combining neutron diffraction, solid-state NMR spectroscopy, and molecular dynamics simulations. Our data indicate that the carprofen derivatives are inserted into the membrane interface, where the exact position and orientation depends on the lipid phase. This knowledge will help to understand the docking of carprofen derivatives to γ-secretase and in the design of new potent drugs. The approach presented here promises to serve as a general guideline how drug/target interactions in membranes can be analyzed in a comprehensive manner.

Combining (Non)linear Optical and Fluorescence Analysis of DiD To Enhance Lipid Phase Recognition.

The widespread interest in phase recognition of lipid membranes has led to the use of different optical techniques to enable differentiation of healthy and not fully functional cells. In this work, we show how the combination of different (non)linear optical methods such as one-photon absorption (OPA), two-photon absorption (TPA), and second harmonic generation (SHG) as well as the study of the fluorescence decay time leads to an enhanced screening of membrane phases using a fluorescent 1,1'-dioctadecyl-3,3,3',3'-tetramethylindocarbocyanine (DiD) probe. In the current study we consider the pure liquid disordered phases of DOPC (dioleoyl- sn-glycero-3-phosphocholine, room temperature) and DPPC (1,2-dipalmitoyl- sn-glycero-3-phosphocholine, 323 K), the solid gel phase of DPPC (298 K), and the liquid ordered phase of a 2:1 binary mixture of sphingomyelin and cholesterol. By means of extensive hybrid quantum mechanics-molecular mechanics calculations and based upon the (non)linear absorption of the embedded probes, it is found that DiD can be used to identify the lipid bilayer phase. The joint TPA and SHG as well as fluorescence analyses qualifies DiD as a versatile probe for phase recognition. In particular, the SHG data obtained by means of hyper-Rayleigh scattering and by electric field induced second harmonic generation reveal differences in polarization of the probe in the different environments. The TPA results finally confirm the particular location of the probe in between the polar headgroup region of the 2:1 SM:Chol mixture in the liquid ordered phase.

Multidrug resistance-associated protein 4 (MRP4) controls ganciclovir intracellular accumulation and contributes to ganciclovir-induced neutropenia in renal transplant patients.

Ganciclovir (GCV) is the cornerstone of cytomegalovirus prevention and treatment in transplant patients. It is associated with problematic adverse hematological effects in this population of immunosuppressed patients, which may lead to dose reduction thus favoring resistance. GCV crosses the membranes of cells, is activated by phosphorylation, and then stops the replication of viral DNA. Its intracellular accumulation might favor host DNA polymerase inhibition, hence toxicity. Following this hypothesis, we investigated the association between a selected panel of membrane transporter polymorphisms and the evolution of neutrophil counts in n=174 renal transplant recipients. An independent population of n=96 renal transplants served as a replication and experiments using HEK293T-transfected cells were performed to validate the clinical findings. In both cohorts, we found a variant in ABCC4 (rs11568658) associated with decreased neutrophil counts following valganciclovir (GCV prodrug) administration (exploratory cohort: ß±SD=-0.68±0.28, p=0.029; replication cohort: ß±SD=-0.84±0.29, p=0.0078). MRP4-expressing cells showed decreased GCV accumulation as compared to negative control cells (transfected with an empty vector) (-61%; p<0.0001). The efflux process was almost abolished in cells expressing MRP4 rs11568658 variant protein. Molecular dynamic simulations of GCV membrane crossing showed a preferred location of the drug just beneath the polar head group region, which supports its interaction with efflux transporters.

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.

A Blue-Light-Emitting BODIPY Probe for Lipid Membranes.

Here we describe a new BODIPY-based membrane probe (1) that provides an alternative to dialkylcarbocyanine dyes, such as DiI-C18, that can be excited in the blue spectral region. Compound 1 has unbranched octadecyl chains at the 3,5-positions and a meso-amino function. In organic solvents, the absorption and emission maxima of 1 are determined mainly by solvent acidity and dipolarity. The fluorescence quantum yield is high and reaches 0.93 in 2-propanol. The fluorescence decays are well fitted with a single-exponential in pure solvents and in small and giant unilamellar vesicles (GUV) with a lifetime of ca. 4 ns. Probe 1 partitions in the same lipid phase as DiI-C18(5) for lipid mixtures containing sphingomyelin and for binary mixtures of dipalmitoylphosphatidylcholine (DPPC) and dioleoylphosphatidylcholine (DOPC). The lipid phase has no effect on the fluorescence lifetime but influences the fluorescence anisotropy. The translational diffusion coefficients of 1 in GUVs and OLN-93 cells are of the same order as those reported for DiI-C18. The directions of the absorption and emission transition dipole moments of 1 are calculated to be parallel. This is reflected in the high steady-state fluorescence anisotropy of 1 in high ordered lipid phases. Molecular dynamic simulations of 1 in a model of the DOPC bilayer indicate that the average angle of the transition moments with respect to membrane normal is ca. 70°, which is comparable with the value reported for DiI-C18.

Lipocarbazole, an efficient lipid peroxidation inhibitor anchored in the membrane.

Lipid peroxidation is a major deleterious effect caused by oxidative stress. It is involved in various diseases such as atherosclerosis, rheumatoid arthritis and neurodegenerative diseases. In order to inhibit lipid peroxidation, antioxidants must efficiently scavenge free radicals and penetrate inside biological membranes. Lipocarbazole has recently been shown to be a powerful antioxidant in solution. Here, we show its powerful capacity as lipid peroxidation inhibitor. Its mechanism of action is rationalized based on molecular dynamics simulations on a biomembrane model, quantum calculations and experimental evaluation. The role of the lipocarbazole side chain is particularly highlighted as a critical chemical feature responsible for its antioxidant activity.

Synergism of antioxidant action of vitamins E, C and quercetin is related to formation of molecular associations in biomembranes.

Vitamins E, C and polyphenols (flavonoids and non-flavonoids) are major natural antioxidants capable of preventing damage generated by oxidative stress. Here we show the capacity of these antioxidants to form non-covalent association within lipid bilayers close to the membrane/cytosol interface. Antioxidant regeneration is significantly enhanced in these complexes.

Benchmarking of Force Fields for Molecule-Membrane Interactions.

Studies of drug-membrane interactions witness an ever-growing interest, as penetration, accumulation, and positioning of drugs play a crucial role in drug delivery and metabolism in human body. Molecular dynamics simulations complement nicely experimental measurements and provide us with new insight into drug-membrane interactions; however, the quality of the theoretical data dramatically depends on the quality of the force field used. We calculated the free energy profiles of 11 molecules through a model dimyristoylphosphatidylcholine (DMPC) membrane bilayer using five force fields, namely Berger, Slipids, CHARMM36, GAFFlipids, and GROMOS 43A1-S3. For the sake of comparison, we also employed the semicontinuous tool COSMOmic. High correlation was observed between theoretical and experimental partition coefficients (log K). Partition coefficients calculated by all-atomic force fields (Slipids, CHARMM36, and GAFFlipids) and COSMOmic differed by less than 0.75 log units from the experiment and Slipids emerged as the best performing force field. This work provides the following recommendations (i) for a global, systematic and high throughput thermodynamic evaluations (e.g., log K) of drugs COSMOmic is a tool of choice due to low computational costs; (ii) for studies of the hydrophilic molecules CHARMM36 should be considered; and (iii) for studies of more complex systems, taking into account all pros and cons, Slipids is the force field of choice.

Synthesis, molecular structure, computational study and in vitro anticancer activity of dinuclear thiolato-bridged pentamethylcyclopentadienyl Rh(III) and Ir(III) complexes.

Neutral dinuclear dithiolato-bridged pentamethylcyclopentadienyl Rh(III) complexes of the type (C5Me5)2Rh2(µ-SR)2Cl2 (R = CH2Ph, 1; R = CH2CH2Ph, 2) and cationic dinuclear trithiolato-bridged pentamethylcyclopentadienyl Rh(III) and Ir(III) complexes of the type [(C5Me5)2M2(µ-SR)3](+) (M = Rh, R = CH2Ph, 3; M = Rh, R = CH2CH2Ph, 5; M = Rh, R = CH2C6H4-p-(t)Bu, 7: M = Ir, R = CH2Ph, 4; M = Ir, R = CH2CH2Ph, 6; M = Ir, R = CH2C6H4-p-(t)Bu, 8) have been synthesized from the chloro-bridged pentamethylcyclopentadienyl Rh(III) and Ir(III) dimers (C5Me5)2M2(µ-Cl)2Cl2 by reaction with the corresponding thiol derivative (RSH). Complexes 3-8 were isolated as chloride salts. All complexes were obtained in good yield and were fully characterized by spectroscopic methods. The molecular structures of the neutral complexes (1 and 2) show interesting features: the two rhodium atoms are bridged by two thiolato ligands with no metal-metal bonds and the pentamethylcyclopentadienyl and chlorido ligands are oriented syn to each other, an uncommon conformation for such dinuclear complexes. These structural features were rationalized using DFT calculations. Additionally, the antiproliferative activity of the complexes was evaluated against the cancerous A2780 (cisplatin sensitive) and A2780cisR (cisplatin resistant) human ovarian cell lines and on the noncancerous HEK293 human embryonic kidney cells. All complexes were found to be active and the cationic iridium complexes , and are particularly cytotoxic, with IC50 values in the nanomolar range (IC50 < 0.1 µM). The catalytic activity of the complexes for the oxidation of glutathione (GSH) to GSSG was evaluated by NMR spectroscopy.

Lipid bilayer membrane affinity rationalizes inhibition of lipid peroxidation by a natural lignan antioxidant.

Lipid peroxidation is a degenerative oxidative process that modifies the structure of membranes, influencing their biological functions. Lignans, natural polyphenolic antioxidants widely distributed in plants, can prevent this membrane damage by free-radical scavenging. Here, we rationalize the difference in lipid peroxidation inhibition activity of argenteane, a natural dilignan isolated from wild nutmeg, and 3,3'-dimethoxy-1,1'-biphenyl-2,2'-diol, which represents the central part of argenteane responsible for its antioxidant activity. Although both compounds have the same capacity to scavenge free radicals, argenteane is a more active inhibitor of lipid peroxidation. We show that both compounds penetrate into DOPC and PLPC lipid bilayers and adopt similar positions and orientations, which therefore does not explain the difference in their lipid peroxidation inhibition activity. However, free energy profiles indicate that argenteane has a significantly higher affinity to the lipid bilayer, and thus a higher effective concentration to scavenge radicals formed during lipid peroxidation. This finding explains the higher activity of argenteane to inhibit lipid peroxidation.