Stefan Paula

Permeation of protons, potassium ions, and small polar molecules through phospholipid bilayers as a function of membrane thickness.

Two mechanisms have been proposed to account for solute permeation of lipid bilayers. Partitioning into the hydrophobic phase of the bilayer, followed by diffusion, is accepted by many for the permeation of water and other small neutral solutes, but transient pores have also been proposed to account for both water and ionic solute permeation. These two mechanisms make distinctively different predictions about the permeability coefficient as a function of bilayer thickness. Whereas the solubility-diffusion mechanism predicts only a modest variation related to bilayer thickness, the pore model predicts an exponential relationship. To test these models, we measured the permeability of phospholipid bilayers to protons, potassium ions, water, urea, and glycerol. Bilayers were prepared as liposomes, and thickness was varied systematically by using unsaturated lipids with chain lengths ranging from 14 to 24 carbon atoms. The permeability coefficient of water and neutral polar solutes displayed a modest dependence on bilayer thickness, with an approximately linear fivefold decrease as the carbon number varied from 14 to 24 atoms. In contrast, the permeability to protons and potassium ions decreased sharply by two orders of magnitude between 14 and 18 carbon atoms, and leveled off, when the chain length was further extended to 24 carbon atoms. The results for water and the neutral permeating solutes are best explained by the solubility-diffusion mechanism. The results for protons and potassium ions in shorter-chain lipids are consistent with the transient pore model, but better fit the theoretical line predicted by the solubility-diffusion model at longer chain lengths.

TWO MECHANISMS OF PERMEATION OF SMALL NEUTRAL MOLECULES AND HYDRATED IONS ACROSS PHOSPHOLIPID BILAYERS

Abstract The Gibbs free energy of dipole or iron permeation of lipid bilayers is calculated as the sum of all electrostatic, solvophobic and specific interactions. Partitioning models are consistent with dipole permeation and some features of ionic permeation, particularly if the solvophobic energy is taken into account. Ionic and dipole permeabilities are extremely sensitive to the ionic/dipolar radius. Despite this sensitivity, calculations of the permeability can be carried out for typical monovalent cations, and provide reasonable estimates, but only for hydrated species. An alternative mechanism proposed for ionic permeation involves the occurrence of transient pore-like defects in lipid bilayers which permit ions to bypass the Born energy barrier. The two alternative hypothesis, partitioning vs. transient pores, can be tested by measuring the ionic and dipolar permeation through bilayers of varying thickness. Experimental observations for both potassium and proton permeability are consistent with the transient pore mechanism for shorter chain lipids, but tend towards the theoretical line for partitioning models for longer chain lipids. Results for small neutral solutes are best explained by the solubility-diffusion mechanism. The proposed method of calculation of the Gibbs free energy of ion or dipole membrane transfer and the liquid membrane permittivity can be effectively used not only in describing the biophysical properties of membranes, but also in extraction processes, pharmaceutical applications and liquid membrane separations.

Permeation of Halide Anions through Phospholipid Bilayers Occurs by the Solubility-Diffusion Mechanism

Two alternative mechanisms are frequently used to describe ionic permeation of lipid bilayers. In the first, ions partition into the hydrophobic phase and then diffuse across (the solubility-diffusion mechanism). The second mechanism assumes that ions traverse the bilayer through transient hydrophilic defects caused by thermal fluctuations (the pore mechanism). The theoretical predictions made by both models were tested for halide anions by measuring the permeability coefficients for chloride, bromide, and iodide as a function of bilayer thickness, ionic radius, and sign of charge. To vary the bilayer thickness systematically, liposomes were prepared from monounsaturated phosphatidylcholines (PC) with chain lengths between 16 and 24 carbon atoms. The fluorescent dye MQAE (N-(ethoxycarbonylmethyl)-6-methoxyquinolinium bromide) served as an indicator for halide concentration inside the liposomes and was used to follow the kinetics of halide flux across the bilayer membranes. The observed permeability coefficients ranged from 10(-9) to 10(-7) cm/s and increased as the bilayer thickness was reduced. Bromide was found to permeate approximately six times faster than chloride through bilayers of identical thickness, and iodide permeated three to four times faster than bromide. The dependence of the halide permeability coefficients on bilayer thickness and on ionic size were consistent with permeation of hydrated ions by a solubility-diffusion mechanism rather than through transient pores. Halide permeation therefore differs from that of a monovalent cation such as potassium, which has been accounted for by a combination of the two mechanisms depending on bilayer thickness.

Development of a new class of aromatase inhibitors: Design, synthesis and inhibitory activity of 3-phenylchroman-4-one (isoflavanone) derivatives

Aromatase (CYP19) catalyzes the aromatization reaction of androgen substrates to estrogens, the last and rate-limiting step in estrogen biosynthesis. Inhibition of aromatase is a new and promising approach to treat hormone-dependent breast cancer. We present here the design and development of isoflavanone derivatives as potential aromatase inhibitors. Structural modifications were performed on the A and B rings of isoflavanones via microwave-assisted, gold-catalyzed annulation reactions of hydroxyaldehydes and alkynes. The in vitro aromatase inhibition of these compounds was determined by fluorescence-based assays utilizing recombinant human aromatase (baculovirus/insect cell-expressed). The compounds 3-(4-phenoxyphenyl)chroman-4-one (1h), 6-methoxy-3-phenylchroman-4-one (2a) and 3-(pyridin-3-yl)chroman-4-one (3b) exhibited potent inhibitory effects against aromatase with IC(50) values of 2.4 μM, 0.26 μM and 5.8 μM, respectively. Docking simulations were employed to investigate crucial enzyme/inhibitor interactions such as hydrophobic interactions, hydrogen bonding and heme iron coordination. This report provides useful information on aromatase inhibition and serves as a starting point for the development of new flavonoid aromatase inhibitors.

Proton Translocation by Bacteriorhodopsin in the Absence of Substantial Conformational Changes

Unlike wild-type bacteriorhodopsin (BR), the BR triple mutant D96G/F171C/F219L has been shown to undergo only minor structural rearrangements during its photocycle. Nonetheless, the mutant is capable of transporting protons at a rate of 125(+/-40) H+/BR per minute under light-saturating conditions. Light adaptation of the triple mutants retinal proceeds in a pH-dependent manner up to a maximum of 63% all-trans. These two findings imply that the transport activity of the triple mutant comprises 66% of the wild-type activity. Time-resolved spectroscopy reveals that the identity and sequence of intermediates in the photocycle of the triple mutant in the all-trans configuration correspond to that of wild-type BR. The only differences relate to a slower rise and decay of the M and O intermediates, and a significant spectral contribution from a 13-cis component. No indication for accumulation of the N intermediate is found under a variety of conditions that normally favor the formation of this species in wild-type BR. The Fourier transform infrared (FTIR) spectrum of the M intermediate in the triple mutant resembles that of wild type. Minor changes in the amide I region during the photocycle suggest that only small movements of the protein backbone occur. Electron microscopy reveals large differences in conformation between the unilluminated state of the mutant protein and wild-type but no light-induced changes in time-resolved measurements. Evidently, proton transport by the triple mutant does not require the major conformational rearrangements that occur on the same time-scale with wild-type. Thus, we conclude that large conformational changes observed in the photocycle of the wild-type and many BR mutants are not a prerequisite for the change in accessibility of the Schiff base nitrogen atom that must occur during vectorial catalysis to allow proton transport.

Molecular determinants of thapsigargin binding by SERCA Ca2+-ATPase: A computational docking study

Thapsigargin (TG) is a potent and commonly used inhibitor of the ion transport activity of sarco/endoplasmic reticulum Ca2+‐ATPases (SERCA). Based on the recently published crystal structures of rabbit muscle SERCA1a in the Ca2+/E1 (E1) and TG/E2 (E2) conformations, we performed computational docking studies to characterize the molecular interactions that govern binding of TG and TG‐analogs by the enzyme. Using the program GOLD (genetic optimization for ligand docking) in combination with the scoring function ChemScore, TG was docked into the binding site of the E1 and E2 conformations of SERCA1a. The docking results revealed a consensus ligand‐binding mode consistent with the crystal structure and showed that hydrophobic interactions are the primary driving force of TG binding by SERCA. Moreover, it was shown that the conformational changes accompanying the E2 to E1 transition in the enzyme likely displace TG from its favored orientation in the binding site, thereby substantially reducing its binding affinity. This finding illustrates on the molecular level how TG may exert its inhibitory effect in binding tightly to the E2 form and preventing it from converting into its E1 form, a requirement for catalytic function. We also docked 9 TG analogs into the E2 conformation of the enzyme. Eight of the analogs adopted a binding mode very similar to that of TG, whereas one compound preferred a different orientation in the binding site. Analysis of the predicted binding affinities showed a good correlation with the experimentally observed inhibitory potencies of the analogs. Docking was also performed with several modeled mutants of SERCA1a, whose phenylalanine residue in position 256 (Phe256) had been modified. The experimentally observed declines in TG sensitivity in most of the Phe256 mutants was qualitatively accounted for and appears, at least in part, be due to a slightly altered TG‐binding mode. Proteins 2004.

Structure-based virtual screening for novel inhibitors of the sarco/endoplasmic reticulum calcium ATPase and their experimental evaluation

A public compound library with 260,000 compounds was screened virtually by computational docking for novel inhibitors of the transmembrane enzyme sarco/endoplasmic reticulum calcium ATPase (SERCA). Docking was performed with the program GOLD in conjunction with a high resolution X-ray crystal structure of SERCA. Compounds that were predicted to be active were tested in bioassays. Nineteen novel compounds were discovered that were capable of inhibiting the ATP hydrolysis activity of SERCA at concentrations below 50 microM. Crucial enzyme/inhibitor interactions were identified by analyzing the docking-predicted binding poses of active compounds. Like other SERCA inhibitors, the newly discovered compounds are of considerable medicinal interest because of their potential for cancer chemotherapy.

Molecular determinants of sarco/endoplasmic reticulum calcium ATPase inhibition by hydroquinone‐based compounds

The ion transport activity of the sarco/endoplasmic reticulum calcium ATPase (SERCA) is specifically and potently inhibited by the small molecule 2,5‐di‐tert‐butylhydroquinone (BHQ). In this study, we investigated the relative importance of the nature and position of BHQs four substituents for enzyme inhibition by employing a combination of experimental and computational techniques. The inhibitory potencies of 21 commercially available or synthesized BHQ derivatives were determined in ATPase activity assays, and 11 compounds were found to be active. Maximum inhibitory potency was observed in compounds with two para hydroxyl groups, whereas BHQ analogues with only one hydroxyl group were still active, albeit with a reduced potency. The results also demonstrated that two alkyl groups were an absolute requirement for activity, with the most potent compounds having 2,5‐substituents with four or five carbon atoms at each position. Using the program GOLD in conjunction with the ChemScore scoring function, the structures of the BHQ analogues were docked into the crystal structure of SERCA mimicking the enzymes E2 conformation. Analysis of the docking results indicated that inhibitor binding to SERCA was primarily mediated by a hydrogen bond between a hydroxyl group and Asp‐59 and by hydrophobic interactions involving the bulky inhibitor alkyl groups. Attempts to dock BHQ into crystal structures corresponding to the E1 conformation of the enzyme failed, because the conformational changes accompanying the E2/E1 transition severely restricted the size of the binding site, suggesting that BHQ stabilizes the enzyme in its E2 form. The potential role of Glu309 in enzyme inhibition is discussed in the context of the computational results. The docking scores correlated reasonably well with the measured inhibitory potencies and allowed the distinction between active and inactive compounds, which is a key requirement for future virtual screening of large compound databases for novel SERCA inhibitors. Proteins 2008.

WATER TRANSPORT BY THE BACTERIAL CHANNEL ALPHA -HEMOLYSIN

This study is an investigation of the ability of the bacterial channel alpha-hemolysin to facilitate water permeation across biological membranes. alpha-Hemolysin channels were incorporated into rabbit erythrocyte ghosts at varying concentrations, and water permeation was induced by mixing the ghosts with hypertonic sucrose solutions. The resulting volume decrease of the ghosts was followed by time-resolved optical absorption at pH 5, 6, and 7. The average single-channel permeability coefficient of alpha-hemolysin for water ranged between 1.3x10-12 cm/s and 1.5x10-12 cm/s, depending on pH. The slightly increased single-channel permeability coefficient at lower pH-values was attributed to an increase in the effective pore size. The activation energy of water transport through the channel was low (Ea=5.4 kcal/mol), suggesting that the properties of water inside the alpha-hemolysin channel resemble those of bulk water. This conclusion was supported by calculations based on macroscopic hydrodynamic laws of laminar water flow. Using the known three-dimensional structure of the channel, the calculations accurately predicted the rate of water flow through the channel. The latter finding also indicated that water permeation data can provide a good estimate of the pore size for large channels.

Novel phenolic inhibitors of the sarco/endoplasmic reticulum calcium ATPase: identification and characterization by quantitative structure–activity relationship modeling and virtual screening

Abstract Inhibitors of the sarco/endoplasmic reticulum calcium ATPase (SERCA) are valuable research tools and hold promise as a new generation of anti-prostate cancer agents. Based on previously determined potencies of phenolic SERCA inhibitors, we created quantitative structure–activity relationship (QSAR) models using three independent development strategies. The obtained QSAR models facilitated virtual screens of several commercial compound collections for novel inhibitors. Sixteen compounds were subsequently evaluated in SERCA activity inhibition assays and 11 showed detectable potencies in the micro- to millimolar range. The experimental results were then incorporated into a comprehensive master QSAR model, whose physical interpretation by partial least squares analysis revealed that properly positioned substituents at the central phenyl ring capable of forming hydrogen bonds and of undergoing hydrophobic interactions were prerequisites for effective SERCA inhibition. The established SAR was in good agreement with findings from previous structural studies, even though it was obtained independently using standard QSAR methodologies.