Assessing the Perturbing Effects of Drugs on Lipid Bilayers Using Gramicidin Channel-Based In Silico and In Vitro Assays

Nonspecific membrane bilayer perturbations by ivermectin underlie SARS-CoV-2 in vitro activity

Since it was proposed as a potential host-directed antiviral agent for SARS-CoV-2, the antiparasitic drug ivermectin has been investigated thoroughly in clinical trials, which have provided insufficient support for its clinical efficacy. To examine the potential for ivermectin to be repurposed as an antiviral agent, we therefore undertook a series of preclinical studies. Consistent with early reports, ivermectin decreased SARS-CoV-2 viral burden in in vitro models at low micromolar concentrations, five- to ten-fold higher than the reported toxic clinical concentration. At similar concentrations, ivermectin also decreased cell viability and increased biomarkers of cytotoxicity and apoptosis. Further mechanistic and profiling studies revealed that ivermectin nonspecifically perturbs membrane bilayers at the same concentrations where it decreases the SARS-CoV-2 viral burden, resulting in nonspecific modulation of membrane-based targets such as G-protein coupled receptors and ion channels. These results suggest that a primary molecular mechanism for the in vitro antiviral activity of ivermectin may be nonspecific membrane perturbation, indicating that ivermectin is unlikely to be translatable into a safe and effective antiviral agent. These results and experimental workflow provide a useful paradigm for performing preclinical studies on (pandemic-related) drug repurposing candidates.

Membrane manipulation by free fatty acids improves microbial plant polyphenol synthesis

Microbial synthesis of nutraceutically and pharmaceutically interesting plant polyphenols represents a more environmentally friendly alternative to chemical synthesis or plant extraction. However, most polyphenols are cytotoxic for microorganisms as they are believed to negatively affect cell integrity and transport processes. To increase the production performance of engineered cell factories, strategies have to be developed to mitigate these detrimental effects. Here, we examine the accumulation of the stilbenoid resveratrol in the cell membrane and cell wall during its production using Corynebacterium glutamicum and uncover the membrane rigidifying effect of this stilbenoid experimentally and with molecular dynamics simulations. A screen of free fatty acid supplements identifies palmitelaidic acid and linoleic acid as suitable additives to attenuate resveratrol's cytotoxic effects resulting in a three-fold higher product titer. This cost-effective approach to counteract membrane-damaging effects of product accumulation is transferable to the microbial production of other polyphenols and may represent an engineering target for other membrane-active bioproducts.

Activity modulation of the Escherichia coli F1FO ATP synthase by a designed antimicrobial peptide via cardiolipin sequestering

Most antimicrobial peptides (AMPs) exert their microbicidal activity through membrane permeabilization. The designed AMP EcDBS1R4 has a cryptic mechanism of action involving the membrane hyperpolarization of Escherichia coli, suggesting that EcDBS1R4 may hinder processes involved in membrane potential dissipation. We show that EcDBS1R4 can sequester cardiolipin, a phospholipid that interacts with several respiratory complexes of E. coli. Among these, F₁F_O ATP synthase uses membrane potential to fuel ATP synthesis. We found that EcDBS1R4 can modulate the activity of ATP synthase upon partition to membranes containing cardiolipin. Molecular dynamics simulations suggest that EcDBS1R4 alters the membrane environment of the transmembrane F_O motor, impairing cardiolipin interactions with the cytoplasmic face of the peripheral stalk that binds the catalytic F₁ domain to the F_O domain. The proposed mechanism of action, targeting membrane protein function through lipid reorganization may open new venues of research on the mode of action and design of other AMPs.

Screening for bilayer-active and likely cytotoxic molecules reveals bilayer-mediated regulation of cell function

A perennial problem encountered when using small molecules (drugs) to manipulate cell or protein function is to assess whether observed changes in function result from specific interactions with a desired target or from less specific off-target mechanisms. This is important in laboratory research as well as in drug development, where the goal is to identify molecules that are unlikely to be successful therapeutics early in the process, thereby avoiding costly mistakes. We pursued this challenge from the perspective that many bioactive molecules (drugs) are amphiphiles that alter lipid bilayer elastic properties, which may cause indiscriminate changes in membrane protein (and cell) function and, in turn, cytotoxicity. Such drug-induced changes in bilayer properties can be quantified as changes in the monomer↔dimer equilibrium for bilayer-spanning gramicidin channels. Using this approach, we tested whether molecules in the Pathogen Box (a library of 400 drugs and drug-like molecules with confirmed activity against tropical diseases released by Medicines for Malaria Venture to encourage the development of therapies for neglected tropical diseases) are bilayer modifiers. 32% of the molecules in the Pathogen Box were bilayer modifiers, defined as molecules that at 10 µM shifted the monomer↔dimer equilibrium toward the conducting dimers by at least 50%. Correlation analysis of the molecules' reported HepG2 cell cytotoxicity to bilayer-modifying potency, quantified as the shift in the gramicidin monomer↔dimer equilibrium, revealed that molecules producing <25% change in the equilibrium had significantly lower probability of being cytotoxic than molecules producing >50% change. Neither cytotoxicity nor bilayer-modifying potency (quantified as the shift in the gramicidin monomer↔dimer equilibrium) was well predicted by conventional physico-chemical descriptors (hydrophobicity, polar surface area, etc.). We conclude that drug-induced changes in lipid bilayer properties are robust predictors of the likelihood of membrane-mediated off-target effects, including cytotoxicity.

Regulation of Gramicidin Channel Function Solely by Changes in Lipid Intrinsic Curvature

Membrane protein function is regulated by the lipid bilayer composition. In many cases the changes in function correlate with changes in the lipid intrinsic curvature (c 0), and c 0 is considered a determinant of protein function. Yet, water-soluble amphiphiles that cause either negative or positive changes in curvature have similar effects on membrane protein function, showing that changes in lipid bilayer properties other than c 0 are important-and may be dominant. To further investigate the mechanisms underlying the bilayer regulation of protein function, we examined how maneuvers that alter phospholipid head groups effective "size"-and thereby c 0-alter gramicidin (gA) channel function. Using dioleoylphospholipids and planar bilayers, we varied the head groups' physical volume and the electrostatic repulsion among head groups (and thus their effective size). When 1,2-dioleyol-sn-glycero-3-phosphocholine (DOPC), was replaced by 1,2-dioleyol-sn-glycero-3-phosphoethanolamine (DOPE) with a smaller head group (causing a more negative c 0), the channel lifetime (τ) is decreased. When the pH of the solution bathing a 1,2-dioleyol-sn-glycero-3-phosphoserine (DOPS) bilayer is decreased from 7 to 3 (causing decreased head group repulsion and a more negative c 0), τ is decreased. When some DOPS head groups are replaced by zwitterionic head groups, τ is similarly decreased. These effects do not depend on the sign of the change in surface charge. In DOPE:DOPC (3:1) bilayers, pH changes from 5→9 to 5→0 (both increasing head group electrostatic repulsion, thereby causing a less negative c 0) both increase τ. Nor do the effects depend on the use of planar, hydrocarbon-containing bilayers, as similar changes were observed in hydrocarbon-free lipid vesicles. Altering the interactions among phospholipid head groups may alter also other bilayer properties such as thickness or elastic moduli. Such changes could be excluded using capacitance measurements and single channel measurements on gA channels of different lengths. We conclude that changes gA channel function caused by changes in head group effective size can be predicted from the expected changes in c 0.

Peptide-induced membrane elastic deformations decelerate gramicidin dimer-monomer equilibration

Gramicidin A (gA) is a hydrophobic pentadecapeptide readily incorporating into a planar bilayer lipid membrane (BLM), thereby inducing a large macroscopic current across the BLM. This current results from ion-channel formation due to head-to-head transbilayer dimerization of gA monomers with rapidly established monomer-dimer equilibrium. Any disturbance of the equilibrium, e.g., by sensitized photoinactivation of a portion of gA monomers, causes relaxation toward a new equilibrium state. According to previous studies, the characteristic relaxation time of the gA-mediated electric current decreases as the current increases upon elevating the gA concentration in the membrane. Here, we report data on the current relaxation kinetics for gA analogs with N-terminal valine replaced by glycine or tyrosine. Surprisingly, the relaxation time increased rather than decreased upon elevation of the total membrane conductance induced by these gA analogs, thus contradicting the classical kinetic scheme. We developed a general theoretical model that accounts for lateral interaction of monomers and dimers mediated by membrane elastic deformations. The modified kinetic scheme of the gramicidin dimerization predicts the reverse dependence of the relaxation time on membrane conductance for gA analogs, with a decreased dimerization constant that is in a good agreement with our experimental data. The equilibration process may be also modulated by incorporation of other peptides ("impurities") into the lipid membrane.

Amphiphiles capsaicin and triton X-100 regulate the chemotherapy drug colchicine's membrane adsorption and ion pore formation potency

Chemotherapy drugs (CDs), e.g. colchicine derivative thiocolchicoside (TCC) and taxol, have been found to physically bind with lipid bilayer membrane and induce ion pores. Amphiphiles capsaicin (Cpsn) and triton X-100 (TX100) are known to regulate lipid bilayer physical properties by altering bilayer elasticity and lipid monolayer curvature. Both CDs and amphiphiles are predicted to physically accommodate alongside lipids in membrane to exert their membrane effects. The effects of their binary accommodation in the lipid membrane are yet to be known. Firstly, we have performed experimental studies to inspect whether membrane adsorption of CDs (colchicine or TCC) gets regulated due to any membrane effects of Cpsn or TX100. We find that the aqueous phase presence of these amphiphiles, known to reduce the membrane stiffness, works towards enhancing the membrane adsorption of CDs. Our recently patented technology 'direct detection method' helps address the membrane adsorption mechanisms. Secondly, in electrophysiology records, we measured the amphiphile effects on the potency of ion channel induction due to CDs. We find that amphiphiles increase the CD induced channel induction potency. Specifically, the membrane conductance, apparently due to the ion channel induction by the TCC, increases substantially due to the Cpsn or TX100 induced alterations of the bilayer physical properties. Thus we may conclude that the binary presence of CDs and amphiphiles in lipid membrane may influence considerably in CD's membrane adsorption, as well as the membrane effects, such as ion pore formation.

The Antimicrobial Peptide Gramicidin S Enhances Membrane Adsorption and Ion Pore Formation Potency of Chemotherapy Drugs in Lipid Bilayers

We recently published two novel findings where we found the chemotherapy drugs (CDs) thiocolchicoside (TCC) and taxol to induce toroidal type ion pores and the antimicrobial peptide gramicidin S (GS) to induce transient defects in model membranes. Both CD pores and GS defects were induced under the influence of an applied transmembrane potential (≈100 mV), which was inspected using the electrophysiology record of membrane currents (ERMCs). In this article, I address the regulation of the membrane adsorption and pore formation of CDs due to GS-induced possible alterations of lipid bilayer physical properties. In ERMCs, low micromolar (≥1 μM) GS concentrations in the aqueous phase were found to cause an induction of defects in lipid bilayers, but nanomolar (nM) concentration GS did nothing. For the binary presence of CDs and GS in the membrane-bathing aqueous phase, the TCC pore formation potency is found to increase considerably due to nM concentration GS in buffer. This novel result resembles our recently reported finding that due to the binary aqueous presence of two AMPs (gramicidin A or alamethicin and GS), the pore or defect-forming potency of either AMP increases considerably. To reveal the underlying molecular mechanisms, the influence of GS (0-400 nM) on the quantitative liposome (membrane) adsorption of CD molecules, colchicine and TCC, was tested. I used the recently patented direct detection method, which helps detect the membrane active agents directly at the membrane in the mole fraction relative to its concentrations in aqueous phase. We find that GS, at concentrations known to do nothing to the lipid bilayer electrical barrier properties in ERMCs, increases the membrane adsorption (membrane uptake) of CDs considerably. This phenomenological finding along with the GS effects on CD-induced membrane conductance increase helps predict an important conclusion. The binary presence of AMPs alongside CDs in the lipid membrane vicinity may work toward enhancing the physical adsorption and pore formation potency of CDs in lipid bilayers. This may help understand why CDs cause considerable cytotoxicity.

Atomistic Characterization of Gramicidin Channel Formation

We investigated gramicidin A (gA) subunit dimerization in lipid bilayers using microsecond-long replica-exchange umbrella sampling simulations, millisecond-long unbiased molecular dynamics simulations, and machine learning. Our simulations led to a dimer structure that is indistinguishable from the experimentally determined gA channel structures, with the two gA subunits joined by six hydrogen bonds (6HB). The simulations also uncovered two additional dimer structures, with different gA-gA stacking orientations that were stabilized by four or two hydrogen bonds (4HB or 2HB). When examining the temporal evolution of the dimerization, we found that two bilayer-inserted gA subunits can form the 6HB dimer directly, with no discernible intermediate states, as well as through paths that involve the 2HB and 4HB dimers.