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.

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

Partitioning of bioactive molecules, including drugs, into cell membranes may produce indiscriminate changes in membrane protein function. As a guide to safe drug development, it therefore becomes important to be able to predict the bilayer-perturbing potency of hydrophobic/amphiphilic drugs candidates. Toward this end, we exploited gramicidin channels as molecular force probes and developed in silico and in vitro assays to measure drugs' bilayer-modifying potency. We examined eight drug-like molecules that were found to enhance or suppress gramicidin channel function in a thick 1,2-dierucoyl-sn-glycero-3-phosphocholine (DC22:1PC) but not in thin 1,2-dioleoyl-sn-glycero-3-phosphocholine (DC18:1PC) lipid bilayer. The mechanism underlying this difference was attributable to the changes in gramicidin dimerization free energy by drug-induced perturbations of lipid bilayer physical properties and bilayer-gramicidin interactions. The combined in silico and in vitro approaches, which allow for predicting the perturbing effects of drug candidates on membrane protein function, have implications for preclinical drug safety assessment.

Molecular Mechanism for Gramicidin Dimerization and Dissociation in Bilayers of Different Thickness.

Membrane protein functions can be altered by subtle changes in the host lipid bilayer physical properties. Gramicidin channels have emerged as a powerful system for elucidating the underlying mechanisms of membrane protein function regulation through changes in bilayer properties, which are reflected in the thermodynamic equilibrium distribution between nonconducting gramicidin monomers and conducting bilayer-spanning dimers. To improve our understanding of how subtle changes in bilayer thickness alter the gramicidin monomer and dimer distributions, we performed extensive atomistic molecular dynamics simulations and fluorescence-quenching experiments on gramicidin A (gA). The free-energy calculations predicted a nonlinear coupling between the bilayer thickness and channel formation. The energetic barrier inhibiting gA channel formation was sharply increased in the thickest bilayer (1,2-dierucoyl-sn-glycero-3-phosphocholine). This prediction was corroborated by experimental results on gramicidin channel activity in bilayers of different thickness. To further explore the mechanism of channel formation, we performed extensive unbiased molecular dynamics simulations, which allowed us to observe spontaneous gA dimer formation in lipid bilayers. The simulations revealed structural rearrangements in the gA subunits and changes in lipid packing, as well as water reorganization, that occur during the dimerization process. Together, the simulations and experiments provide new, to our knowledge, insights into the process and mechanism of gramicidin channel formation, as a prototypical example of the bilayer regulation of membrane protein function.

Gramicidin Increases Lipid Flip-Flop in Symmetric and Asymmetric Lipid Vesicles.

Unlike most transmembrane proteins, phospholipids can migrate from one leaflet of the membrane to the other. Because this spontaneous lipid translocation (flip-flop) tends to be very slow, cells facilitate the process with enzymes that catalyze the transmembrane movement and thereby regulate the transbilayer lipid distribution. Nonenzymatic membrane-spanning proteins with unrelated primary functions have also been found to accelerate lipid flip-flop in a nonspecific manner and by various hypothesized mechanisms. Using deuterated phospholipids, we examined the acceleration of flip-flop by gramicidin channels, which have well-defined structures and known functions, features that make them ideal candidates for probing the protein-membrane interactions underlying lipid flip-flop. To study compositionally and isotopically asymmetric proteoliposomes containing gramicidin, we expanded a recently developed protocol for the preparation and characterization of lipid-only asymmetric vesicles. Channel incorporation, conformation, and function were examined with small angle x-ray scattering, circular dichroism, and a stopped-flow spectrofluorometric assay, respectively. As a measure of lipid scrambling, we used differential scanning calorimetry to monitor the effect of gramicidin on the melting transition temperatures of the two bilayer leaflets. The two calorimetric peaks of the individual leaflets merged into a single peak over time, suggestive of scrambling, and the effect of the channel on the transbilayer lipid distribution in both symmetric 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine and asymmetric 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine/1,2-dimyristoyl-sn-glycero-3-phosphocholine vesicles was quantified from proton NMR measurements. Our results show that gramicidin increases lipid flip-flop in a complex, concentration-dependent manner. To determine the molecular mechanism of the process, we used molecular dynamics simulations and further computational analysis of the trajectories to estimate the extent of membrane deformation. Together, the experimental and computational approaches were found to constitute an effective means for studying the effects of transmembrane proteins on lipid distribution in both symmetric and asymmetric model membranes.

Antidepressants are modifiers of lipid bilayer properties.

The two major classes of antidepressants, tricyclic antidepressants (TCAs) and selective serotonin reuptake inhibitors (SSRIs), inhibit neurotransmitter reuptake at synapses. They also have off-target effects on proteins other than neurotransmitter transporters, which may contribute to both desired changes in brain function and the development of side effects. Many proteins modulated by antidepressants are bilayer spanning and coupled to the bilayer through hydrophobic interactions such that the conformational changes underlying their function will perturb the surrounding lipid bilayer, with an energetic cost (ΔG def) that varies with changes in bilayer properties. Here, we test whether changes in ΔG def caused by amphiphilic antidepressants partitioning into the bilayer are sufficient to alter membrane protein function. Using gramicidin A (gA) channels to probe whether TCAs and SSRIs alter the bilayer contribution to the free energy difference for the gramicidin monomer⇔dimer equilibrium (representing a well-defined conformational transition), we find that antidepressants alter gA channel activity with varying potency and no stereospecificity but with different effects on bilayer elasticity and intrinsic curvature. Measuring the antidepressant partition coefficients using isothermal titration calorimetry (ITC) or cLogP shows that the bilayer-modifying potency is predicted quite well by the ITC-determined partition coefficients, and channel activity is doubled at an antidepressant/lipid mole ratio of 0.02-0.07. These results suggest a mechanism by which antidepressants could alter the function of diverse membrane proteins by partitioning into cell membranes and thereby altering the bilayer contribution to the energetics of membrane protein conformational changes.