Continuum solvent model studies of the interactions of an anticonvulsant drug with a lipid bilayer

Polarity of Hydrated Phosphatidylcholine Headgroups

The headgroup (H) stratum (sometimes called the polar region) of membrane bilayers is a relevant yet poorly understood solvation phase for small molecules and macromolecules interacting with the membranes. Solvation of compounds in bilayer strata is characterized experimentally by wide- and small-angle X-ray scattering, neutron diffraction, and various NMR techniques. The quantification is tedious and only available for a limited set of small molecules. Our recently published model of liposome partitioning of small molecules shows that solvation of compounds in the H-stratum of fluid phosphatidylcholine (PC) bilayers correlates well with their solvation in hydrated diacetyl phosphatidylcholine (DAcPC), and solvation in the core (C) depends in a similar way on that in n-hexadecane. These two correlations became a basis for a model describing the location of compounds in the H- and C-strata and at the connecting interface as a nonlinear function of the fragment solvation characteristics of the compounds. In this study, refractivity of hydrated DAcPC phases with varying water contents was measured and polarity was determined using the steady-state fluorescence of indole and Nile Red. The results were compared with the published data obtained by other techniques for PC bilayers in liposomes or on solid supports. The demonstrated qualitative agreement, as well as the polarity and refractivity dependencies on the DAcPC concentration, supports the suitability of hydrated DAcPC as the H-stratum surrogate. Interestingly, depending on hydrations typical for the H-strata of fluid PC bilayers, the dielectric constant could decrease significantly from 31.0 to 7.3 for 16 and 8 water molecules per headgroup, respectively. Although additional experiments are needed for confirmation, this observation could help set proper dielectric constant magnitudes in continuum-based computational models of accumulation and crossing of the PC bilayers with varying hydration levels thanks to the temperature or the structure of fatty acid chains.

A Permeability Study of O2 and the Trace Amine p-Tyramine through Model Phosphatidylcholine Bilayers

We study here the permeability of the hydrophobic O2 molecule through a model DPPC bilayer at 323K and 350K, and of the trace amine p-tyramine through PC bilayers at 310K. The tyramine results are compared to previous experimental work at 298K. Nonequilibrium work methods were used in conjunction to simultaneously obtain both the potential of mean force (PMF) and the position dependent transmembrane diffusion coefficient, D(z), from the simulations. These in turn were used to calculate the permeability coefficient, P, through the inhomogeneous solubility-diffusion model. The results for O2 are consistent with previous simulations, and agree with experimentally measured P values for PC bilayers. A temperature dependence in the permeability of O2 through DPPC was obtained, with P decreasing at higher temperatures. Two relevant species of p-tyramine were simulated, from which the PMF and D(z) were calculated. The charged species had a large energetic barrier to crossing the bilayer of ~ 21 kcal/mol, while the uncharged, deprotonated species had a much lower barrier of ~ 7 kcal/mol. The effective in silico permeability for p-tyramine was calculated by applying three approximations, all of which gave nearly identical results (presented here as a function of the pKa). As the permeability value calculated from simulation was highly dependent on the pKa of the amine group, a further pKa study was performed that also varied the fraction of the uncharged and zwitterionic p-tyramine species. Using the experimental P value together with the simulated results, we were able to label the phenolic group as responsible for the pKa1 and the amine for the pKa2, that together represent all of the experimentally measured pKa values for p-tyramine. This agrees with older experimental results, in contrast to more recent work that has suggested there is a strong ambiguity in the pKa values.

Structure-based prediction of drug distribution across the headgroup and core strata of a phospholipid bilayer using surrogate phases

Solvation of drugs in the core (C) and headgroup (H) strata of phospholipid bilayers affects their physiological transport rates and accumulation. These characteristics, especially a complete drug distribution profile across the bilayer strata, are tedious to obtain experimentally, to the point that even simplified preferred locations are only available for a few dozen compounds. Recently, we showed that the partition coefficient (P) values in the system of hydrated diacetyl phosphatidylcholine (DAcPC) and n-hexadecane (C16), as surrogates of the H- and C-strata of the bilayer composed of the most abundant mammalian phospholipid, PC, agree well with the preferred bilayer location of compounds. High P values are typical for lipophiles accumulating in the core, and low P values are characteristic of cephalophiles preferring the headgroups. This simple pattern does not hold for most compounds, which usually have more even distribution and may also accumulate at the H/C interface. To model complete distribution, the correlates of solvation energies are needed for each drug state in the bilayer: (1) for the H-stratum it is the DAcPC/W P value, calculated as the ratio of the C16/W and C16/DAcPC (W for water) P values; (2) for the C-stratum, the C16/W P value; (3) for the H/C interface, the P values for all plausible molecular poses are characterized using the fragment DAcPC/W and C16/W solvation parameters for the parts of the molecule embedded in the H- and C-strata, respectively. The correlates, each scaled by two Collander coefficients, were used in a nonlinear, mass-balance based model of intrabilayer distribution, which was applied to the easily measurable overall P values of compounds in the DMPC (M = myristoyl) bilayers and monolayers as the dependent variables. The calibrated model for 107 neutral compounds explains 94% of experimental variance, achieves similar cross-validation levels, and agrees well with the nontrivial, experimentally determined bilayer locations for 27 compounds. The resulting structure-based prediction system for intrabilayer distribution will facilitate more realistic modeling of passive transport and drug interactions with those integral membrane proteins, which have the binding sites located in the bilayer, such as some enzymes, influx and efflux transporters, and receptors. If only overall bilayer accumulation is of interest, the 1-octanol/W P values suffice to model the studied set.

Structural determinants of drug partitioning in surrogates of phosphatidylcholine bilayer strata

The knowledge of drug concentrations in bilayer headgroups, core, and at the interface between them is a prerequisite for quantitative modeling of drug interactions with many membrane-bound transporters, metabolizing enzymes and receptors, which have the binding sites located in the bilayer. This knowledge also helps understand the rates of trans-bilayer transport because balanced interactions of drugs with the bilayer strata lead to high rates, while excessive affinities for any stratum cause a slowdown. Experimental determination of bilayer location is so tedious and costly that the data are only available for some fifty compounds. To extrapolate these valuable results to more compounds at a higher throughput, surrogate phases have been used to obtain correlates of the drug affinities for individual strata. We introduced a novel system, consisting of a diacetyl phosphatidylcholine (DAcPC) solution with the water content of the fluid bilayer as the headgroup surrogate and n-hexadecane (C16) representing the core. The C16/DAcPC partition coefficients were measured for 113 selected compounds, containing structural fragments that are frequently occurring in approved drugs. The data were deconvoluted into the ClogP-based fragment solvation characteristics and processed using a solvatochromic correlation. Increased H-bond donor ability and excess molar refractivity of compounds promote solvation in the DAcPC phase as compared to bulk water, contrary to H-bond acceptor ability, dipolarity/polarizability, and volume. The results show that aromates have more balanced distribution in bilayer strata, and thus faster trans-bilayer transport, than similar alkanes. This observation is in accordance with the frequent occurrence of aromatic rings in approved drugs and with the role of rigidity of drug molecules in promoting intestinal absorption. Bilayer locations, predicted using the C16/DAcPC system, are in excellent agreement with available experimental data, in contrast to other surrogate systems.

Structural determinants of drug partitioning in n-hexadecane/water system

Surrogate phases have been widely used as correlates for modeling transport and partitioning of drugs in biological systems, taking advantage of chemical similarity between the surrogate and the phospholipid bilayer as the elementary unit of biological phases, which is responsible for most of the transport and partitioning. Solvation in strata of the phospholipid bilayer is an important drug characteristic because it affects the rates of absorption and distribution, as well as the interactions with the membrane proteins having the binding sites located inside the bilayer. The bilayer core can be emulated by n-hexadecane (C16), and the headgroup stratum is often considered a hydrophilic phase because of the high water content. Therefore, we tested the hypothesis that the C16/water partition coefficients (P) can predict the bilayer locations of drugs and other small molecules better than other surrogate systems. Altogether 514 PC16/W values for nonionizable (458) and completely ionized (56) compounds were collected from the literature or measured, when necessary. With the intent to create a fragment-based prediction system, the PC16/W values were factorized into the fragment solvation parameters (f) and correction factors based on the ClogP fragmentation scheme. A script for the PC16/W prediction using the ClogP output is provided. To further expand the prediction system and reveal solvation differences, the fC16/W values were correlated with their more widely available counterparts for the 1-octanol/water system (O/W) using solvatochromic parameters. The analysis for 50 compounds with known bilayer location shows that the available and predicted PC16/W and PO/W values alone or the PC16/O values representing their ratio do not satisfactorily predict the preference for drug accumulation in bilayer strata. These observations indicate that the headgroups stratum, albeit well hydrated, does not have solvation characteristics similar to water and is also poorly described by the O/W partition characteristics.

Identifying an uptake mechanism for the antiepileptic and bipolar disorder treatment valproic acid using the simple biomedical model Dictyostelium

Valproic acid (VPA) is the most highly prescribed epilepsy treatment worldwide and is also used to prevent bipolar disorder and migraine. Surprisingly, very little is known about its mechanisms of cellular uptake. Here, we employ a range of cellular, molecular and genetic approaches to characterize VPA uptake using a simple biomedical model, Dictyostelium discoideum. We show that VPA is taken up against an electrochemical gradient in a dose-dependent manner. Transport is protein-mediated, dependent on pH and the proton gradient and shows strong substrate structure specificity. Using a genetic screen, we identified a protein homologous to a mammalian solute carrier family 4 (SLC4) bicarbonate transporter that we show is involved in VPA uptake. Pharmacological and genetic ablation of this protein reduces the uptake of VPA and partially protects against VPA-dependent developmental effects, and extracellular bicarbonate competes for VPA uptake in Dictyostelium. We further show that this uptake mechanism is likely to be conserved in both zebrafish (Danio rerio) and Xenopus laevis model systems. These results implicate, for the first time, an uptake mechanism for VPA through SLC4-catalysed activity.

Valproic acid alters mitochondrial cholesterol transport in Y1 adrenocortical cells

Several reports suggest putative interactions between valproic acid (VPA) treatment and the hypothalamus-pituitary-adrenal axis. Given that VPA alters mitochondrial functions, an action of this drug on a mitochondrial process such as steroid synthesis in adrenal cells should be expected. In order to disclose a putative action of VPA on the adrenocortical cell itself we evaluated VPA effects on regulatory steps of the acute stimulation of steroidogenesis in Y1 adrenocortical cells. This study demonstrates that VPA increases progesterone production in non-stimulated cells without inducing the levels of Steroidogenic Acute Regulatory (StAR) protein, which facilitates cholesterol transport. This result suggests that VPA increases mitochondrial cholesterol transport through a StAR-independent mechanism and is further supported by the fact that in isolated mitochondria VPA stimulates exogenous cholesterol metabolization to progesterone. VPA also reduces the cAMP-mediated increase of the StAR protein, mRNA levels, promoter activity and progesterone production. In summary, the present data show that VPA can alter steroid production in adrenal cells by a complex mechanism that mainly involves an action on cholesterol access to the inner mitochondrial membrane. The VPA-mediated increase of basal steroidogenesis could be linked to the increase of basal cortisolemia described in patients under VPA treatment.

High-throughput pKa screening and prediction amenable for ADME profiling

Recent technological advances have made it possible for several new pK(a) assays to be used in drug screening. In this review, a critical overview is provided of current new methodologies for high-throughput screening and prediction of pK(a). Typical applications of using pK(a )constants and charge state for absorption, distribution, metabolism and excretion (ADME) profiling and quantitative structure-activity relationship modelling complements the methodological comparisons and discussions. The experimental methods discussed include high-throughput screening of pK(a) by multiplexed capillary with ultraviolet absorbance detection on a 96-capillary format instrument, capillary electrophoresis and mass spectrometry (CEMS) based on sample pooling, determination of pK(a) by pH gradient high-performance liquid chromatography, and measurement of pK(a) by a mixed-buffer liner pH gradient system. Comparisons of the different experimental assays are made with emphasis on the newly developed CEMS method. The current status and recent progress in computational approaches to pK(a) prediction are also discussed. In particular, the accuracy limits of simple fragment-based approaches as well as quantum mechanical methods are addressed. Examples of pK(a) prediction from in-house drug candidates as well as commercially available drug molecules are shown and an outline is provided for how drug discovery companies can integrate experiments with computational approaches for increased applications for ADME profiling.

Membrane-bound ARF1 peptide: interpretation of neutron diffraction data by molecular dynamics simulation methods

Adenosine diphosphate ribosylation factor-1 (ARF1) is activated by cell membrane binding of a self-folding N-terminal domain. We have previously presented four possible conformations of the membrane bound, human ARF1 N-terminal peptide in planar lipid bilayers of DOPC and DOPG (7:3 molar ratio), determined from lamellar neutron diffraction and circular dichroism data. In this paper we analyse the four possible conformations by molecular dynamics simulations. The aim of these simulations was to use MD to distinguish which of the four possible membrane bound structures was the most likely. The most likely conformation was determined according to the following criteria: (a) location of label positions on the peptide in relation to the bilayer, (b) lowest mean square displacement from the initial structure, (c) lowest system energy, (d) most peptide-lipid headgroup hydrogen bonding, (e) analysis of phi/psi angles of the peptide. These findings demonstrate the application of molecular dynamics simulations to explore neutron diffraction data.

[Drugs for status epilepticus treatment]

The pharmacokinetics and pharmacodynamics of major antiepileptic agents are presented. The onset of action and the factors leading to extraction across the blood brain barrier are described as well as the mechanism and extent of metabolism, and the main interactions with other drugs. For each class, the dosing scheme and practical issues related to administration are described, based on evidence when available in the literature.

Time-dependent interactions of the two porphyrinic compounds chlorin e6 and mono-L-aspartyl-chlorin e6 with phospholipid vesicles probed by NMR spectroscopy

The distribution processes of chlorin e6 (CE) and monoaspartyl-chlorin e6 (MACE) between the outer and inner phospholipid monolayers of 1,2-dioleoyl-phosphatidylcholine (DOPC) vesicles were monitored by 1H NMR spectroscopy through analysis of chemical shifts and line widths of the DOPC vesicle resonances. Chlorin adsorption to the outer vesicle monolayer induced changes in the DOPC 1H NMR spectrum. Most pronounced was a split of the N-methyl choline resonance, allowing for separate analysis of inner and outer vesicle layers. Transbilayer distribution of the chlorin compounds was indicated by time-dependent characteristic spectral changes of the DOPC resonances. Kinetic parameters for the flip-flop processes, that is, half-lives and rate constants, were obtained from the experimental data points. In comparison to CE, MACE transbilayer movement was significantly reduced, with MACE remaining more or less attached to the outer membrane layer. The distribution coefficients for CE and MACE between the vesicular and aqueous phase were determined. Both CE and MACE exhibited a high affinity for the vesicular phase. For CE, a positive correlation was found between transfer rate and increasing molar ratio CE/DOPC. Enhanced membrane rigidity induced by increasing amounts of cholesterol into the model membrane was accompanied by a decrease of CE flip-flop rates across the membrane. The present study shows that the movement of porphyrins across membranes can efficiently be investigated by 1H NMR spectroscopy and that small changes in porphyrin structure can have large effects on membrane kinetics.

Computational model for predicting chemical substituent effects on passive drug permeability across parallel artificial membranes

Drug permeability is often a limiting step in drug action, requiring chemical optimization of a drug candidate to improve this property. Such optimization is typically performed in the context of a congeneric series, where substituents are varied to optimize the target property. Motivated by this need the present work examines the influence of chemical substituents on passive permeability (log P pass) across parallel artificial membranes (PAMPA) undertaken for three congeneric series of compounds; benzoic acids, pyridines and quinolines. PAMPA showed pyridine and quinoline to have high permeability and chemical substituents to typically reduce the permeability. On the contrary, benzoic acid showed poor permeability and chemical substituents typically increased the permeability. To quantitate these effects with respect to physical properties, models were built to explain and predict the permeability of these classes of compounds based on computed molecular descriptors. Models for the benzoic acid series in the ionized state indicated the solvent accessible surface area, cavity dispersion and the free energy of solvation in hexane as well as in water to dominate permeability. However, when the acid group is treated as neutral, the free energy of solvation in water, the fraction polar surface area, the polar surface area and difference in the free energy of solvation between hexane and water were important; these terms, among others, were also important for the neutral pyridine-quinoline series. Considering that the permeability of the benzoic acid series is about 2 orders of magnitude lower than the pyridines and quinolines and that a shift of approximately two pH units in the p K a of the acid group of benzoic acid will allow for the neutral species of the molecule to dominate under experimental conditions (pH = 6.5), these results suggest that the additional energy barrier associated with permeation of the benzoic acid series is associated with the need to protonate the acidic group, thereby forming the neutral species which may then cross the hydrophobic region of the membrane.

Valproate and copper accelerate TRH-like peptide synthesis in male rat pancreas and reproductive tissues

Treatment with valproate (Valp) facilitates the synthesis of TRH-like peptides (pGlu-X-Pro-NH(2)) in rat brain where "X" can be any amino acid residue. Because high levels of TRH-like peptides occur in the pancreas and pGlu-Glu-Pro-NH(2) (Glu-TRH) has been shown to be a fertilization promoting peptide, we hypothesized that these peptides mediate some of the metabolic and reproductive side effects of Valp. Male WKY rats were treated with Valp acutely (AC), chronically (CHR) or chronically followed by a 2 day withdrawal (WD). AC, CHR and WD treatments significantly altered TRH and/or TRH-like peptide levels in pancreas and reproductive tissues. Glu-TRH was the predominant TRH-like peptide in epididymis, consistent with its fertilization promoting activity. Glu-TRH levels in the epididymis increased 3-fold with AC Valp. Phe-TRH, the most abundant TRH-like peptide in the pancreas, increased 4-fold with AC Valp. Phe-TRH inhibits both basal and TRH-stimulated insulin release. Large dense core vesicles (LDCV's) contain a copper-dependent enzyme responsible for the post-translational processing of precursors of TRH and TRH-like peptides. Copper (500 microM) increased the in vitro C-terminal amidation of TRH-like peptides by 8- and 4-fold during 24 degrees C incubation of homogenates of pancreas and testis, respectively. Valp (7 microM) accelerated 3-fold the processing of TRH and TRH-like peptide precursors in pancreatic LDCV's incubated at 24 degrees C. We conclude that copper, an essential cofactor for TRH and TRH-like peptide biosynthesis that is chelated by Valp, mediates some of the metabolic and reproductive effects of Valp treatment via acceleration of intravesicular synthesis and altered release of these peptides.

Cooperative adsorption of proteins onto lipid membranes

The adsorption of proteins onto a lipid membrane depends on and thus reflects the energetics of the underlying substrate. This is particularly relevant for mixed membranes that contain lipid species with different affinities for the adsorbed proteins. In this case, there is an intricate interplay between lateral membrane organization and the number of adsorbed proteins. Most importantly, proteins often tend to enhance the propensity of the lipid mixture to form clusters, domains, or to macroscopically phase separate. Sigmoidal binding isotherms are the typical signature of the corresponding cooperativity in protein adsorption. We discuss the underlying thermodynamic basis, and compare various theoretical binding models for protein adsorption onto mixed membranes. We also present experimental data for the adsorption of the C2A protein motif and analyse to what extent these data reflect cooperative binding.

Physiologically based pharmacokinetic modelling 2: Predicting the tissue distribution of acids, very weak bases, neutrals and zwitterions

A key component of whole body physiologically based pharmacokinetic (WBPBPK) models is the tissue-to-plasma water partition coefficients (Kpu's). The predictability of Kpu values using mechanistically derived equations has been investigated for 7 very weak bases, 20 acids, 4 neutral drugs and 8 zwitterions in rat adipose, bone, brain, gut, heart, kidney, liver, lung, muscle, pancreas, skin, spleen and thymus. These equations incorporate expressions for dissolution in tissue water and, partitioning into neutral lipids and neutral phospholipids. Additionally, associations with acidic phospholipids were incorporated for zwitterions with a highly basic functionality, or extracellular proteins for the other compound classes. The affinity for these cellular constituents was determined from blood cell data or plasma protein binding, respectively. These equations assume drugs are passively distributed and that processes are nonsaturating. Resultant Kpu predictions were more accurate when compared to published equations, with 84% as opposed to 61% of the predicted values agreeing with experimental values to within a factor of 3. This improvement was largely due to the incorporation of distribution processes related to drug ionisation, an issue that is not addressed in earlier equations. Such advancements in parameter prediction will assist WBPBPK modelling, where time, cost and labour requirements greatly deter its application.

Sodium valproate inhibits glucose transport and exacerbates Glut1-deficiency in vitro

Anticonvulsant sodium valproate interferes with brain glucose metabolism. The mechanism underlying such metabolic disturbance is unclear. We tested the hypothesis that sodium valproate interferes with cellular glucose transport with a focus on Glut1 since glucose transport across the blood-brain barrier relies on this transporter. Cell types enriched with Glut1 expression including human erythrocytes, human skin fibroblasts, and rat astrocytes were used to study the effects of sodium valproate on glucose transport. Sodium valproate significantly inhibited Glut1 activity in normal and Glut1-deficient erythrocytes by 20%-30%, causing a corresponding reduction of Vmax of glucose transport. Similarly, in primary astrocytes as well as in normal and Glut1-deficient fibroblasts, sodium valproate inhibited glucose transport by 20%-40% (P < 0.05), accompanied by an up to 60% downregulation of GLUT1 mRNA expression (P < 0.05). In conclusion, sodium valproate inhibits glucose transport and exacerbates Glut1 deficiency in vitro. Our findings imply the importance of prudent use of sodium valproate for patients with compromised Glut1 function.

Thermodynamic and electrostatic properties of ternary Oxytricha nova TEBP-DNA complex

Telomeres constitute the nucleoprotein ends of eukaryotic chromosomes which are essential for their proper function. Telomere end binding protein (TEBP) from Oxytricha nova was among the first telomeric proteins, which were well characterized biologically. TEBP consists of two protein subunits (alpha, beta) and forms a ternary complex with single stranded telomeric DNA containing tandem repeats TTTTGGGG. This work presents the characterization of the thermodynamic and electrostatic properties of this complex by computational chemistry methods (continuum Poisson-Boltzmann and solvent accessible surface calculations). Our calculations give a new insight into molecular properties of studied system. Based on the thermodynamic analysis we provide a rationale for the experimental observation that alpha and ssDNA forms a binary complex and the beta subunit joins alpha:ssDNA complex only after the latter is formed. Calculations of distribution of the molecular electrostatic potential for protein subunits alone and for all possible binary complexes revealed the important role of the "guiding funnel" potential generated by alpha:ssDNA complex. This potential may help the beta subunit to dock to the already formed alpha:DNA intermediate in highly steric and electrostatic favorable manner. Our pK(a) calculations of TEBP are able to explain the experimental mobility shifts of the complex in electrophoretic non-denaturating gels.

Free diffusion of steroid hormones across biomembranes: a simplex search with implicit solvent model calculations

Steroid hormones such as progesterone, testosterone, and estradiol are derived from cholesterol, a major constituent of biomembranes. Although the hormones might be expected to associate with the bilayer in a fashion similar to that of cholesterol, their biological action in regulating transcription of target genes involves transbilayer transfer by free diffusion, which is not observed for cholesterol. We used a novel combination of a continuum-solvent model and the downhill simplex search method for the calculation of the free energy of interaction of these hormones with lipid membranes, and compared these values to that of cholesterol-membrane interaction. The hormones were represented in atomic detail and the membrane as a structureless hydrophobic slab embedded in implicit water. A deep free-energy minimum of approximately -15 kcal/mol was obtained for cholesterol at its most favorable location in the membrane, whereas the most favorable locations for the hormones were associated with shallower minima of -5.0 kcal/mol or higher. The free-energy difference, which is predominantly due to the substitution of cholesterol's hydrophobic tail with polar groups, explains the different manner in which cholesterol and the hormones interact with the membrane. Further calculations were conducted to estimate the rate of transfer of the hormones from the aqueous phase into hexane, and from hexane back into the aqueous phase. The calculated rates agreed reasonably well with measurements in closely related systems. Based on these calculations, we suggest putative pathways for the free diffusion of the hormones across biomembranes. Overall, the calculations imply that the hormones may rapidly cross biomembrane barriers. Implications for gastrointestinal absorption and transfer across the blood-brain barrier and for therapeutic uses are discussed.

Combined quantitative and mechanistic study of drug-membrane interactions using a novel 2H NMR approach

Several analytical methods are available for determining the partition coefficients of drug compounds in model phospholipid membranes, but such methods provide little information at the molecular level about how the membrane affinity of drugs relates to their interactions with the lipid molecules. A new (2)H nuclear magnetic resonance (NMR) approach has been developed here that quantifies the affinity of (2)H-labeled small molecules for different phospholipid membranes and, simultaneously, provides information on the mechanism of the drug-membrane interaction. In the example given, (2)H NMR analysis of a weakly basic ion pump inhibitor found that the drug partitioned preferentially into membranes of predominantly unsaturated or short-chain phospholipids. The (2)H NMR analysis also suggested that the membrane specificity of the drug was directly correlated to the ability of its phenyl moiety to penetrate into the interior of the lipid bilayer. The (2)H NMR approach could be of value in guiding medicinal chemistry toward or away from structures promoting interactions with specific types of biological membranes.

Interactions of the M2delta segment of the acetylcholine receptor with lipid bilayers: a continuum-solvent model study

M2delta, one of the transmembrane segments of the nicotinic acetylcholine receptor, is a 23-amino-acid peptide, frequently used as a model for peptide-membrane interactions. In this and the companion article we describe studies of M2delta-membrane interactions, using two different computational approaches. In the present work, we used continuum-solvent model calculations to investigate key thermodynamic aspects of its interactions with lipid bilayers. M2delta was represented in atomic detail and the bilayer was represented as a hydrophobic slab embedded in a structureless aqueous phase. Our calculations show that the transmembrane orientation is the most favorable orientation of the peptide in the bilayer, in good agreement with both experimental and computational data. Moreover, our calculations produced the free energy of association of M2delta with the lipid bilayer, which, to our knowledge, has not been reported to date. The calculations included 10 structures of M2delta, determined by nuclear magnetic resonance in dodecylphosphocholine micelles. All the structures were found to be stable inside the lipid bilayer, although their water-to-membrane transfer free energies differed by as much as 12 kT. Although most of the structures were roughly linear, a single structure had a kink in its central region. Interestingly, this structure was found to be the most stable inside the lipid bilayer, in agreement with molecular dynamics simulations of the peptide and with the recently determined structure of the intact receptor. Our analysis showed that the kink reduced the polarity of the peptide in its central region by allowing the electrostatic masking of the Gln13 side chain in that area. Our calculations also showed a tendency for the membrane to deform in response to peptide insertion, as has been previously found for the membrane-active peptides alamethicin and gramicidin. The results are compared to Monte Carlo simulations of the peptide-membrane system, as presented in the accompanying article.

Permeation across hydrated DPPC lipid bilayers: simulation of the titrable amphiphilic drug valproic acid

Valproic acid is a short branched fatty acid used as an anticonvulsant drug whose therapeutic action has been proposed to arise from membrane-disordering properties. Static and kinetic properties of valproic acid interacting with fully hydrated dipalmitoyl phosphatidylcholine lipid bilayers are studied using molecular-dynamics simulations. We calculate spatially resolved free energy profiles and local diffusion coefficients using the distance between the bilayer and valproic acid respective centers-of-mass along the bilayer normal as reaction coordinate. To investigate the pH dependence, we calculate profiles for the neutral valproic acid as well as its water-soluble anionic conjugate base valproate. The local diffusion constants for valproate/valproic acid along the bilayer normal are found to be approximately 10(-6) to 10(-5) cm2 s(-1). Assuming protonation of valproic acid upon association with--or insertion into--the lipid bilayer, we calculate the permeation coefficient to be approximately 2.0 10(-3) cm s(-1), consistent with recent experimental estimates of fast fatty acid transport. The ability of the lipid bilayer to sustain local defects such as water intrusions stresses the importance of going beyond mean field and taking into account correlation effects in theoretical descriptions of bilayer translocation processes.

Effective energy function for proteins in lipid membranes

A simple extension of the EEF1 energy function to heterogeneous membrane-aqueous media is proposed. The extension consists of (a) development of solvation parameters for a nonpolar phase using experimental data for the transfer of amino acid side-chains from water to cyclohexane, (b) introduction of a heterogeneous membrane-aqueous system by making the reference solvation free energy of each atom dependent on the vertical coordinate, (c) a modification of the distance-dependent dielectric model to account for reduced screening of electrostatic interactions in the membrane, and (d) an adjustment of the EEF1 aqueous model in light of recent calculations of the potential of mean force between amino acid side-chains in water. The electrostatic model is adjusted to match experimental observations for polyalanine, polyleucine, and the glycophorin A dimer. The resulting energy function (IMM1) reproduces the preference of Trp and Tyr for the membrane interface, gives reasonable energies of insertion into or adsorption onto a membrane, and allows stable 1-ns MD simulations of the glycophorin A dimer. We find that the lowest-energy orientation of melittin in bilayers varies, depending on the thickness of the hydrocarbon layer.

Thermodynamic Study of Small Hydrophobic Ions at the Water–Lipid Interface

The thermodynamics of binding of two small hydrophobic ions such as norharman and tryptophan to neutral and negatively charged small unilamellar vesicles was investigated at pH 7.4 using fluorescence spectroscopy. Vesicles were formed at room temperature from dimyristoyl phosphatidylcholine (DMPC) or DMPC/dimyristoylphosphatidic acid and DMPC/dimyristoylphosphatidylglycerol. The changes in fluorescence properties were used to obtain association isotherms at variable membrane surface negative charge and at different ionic strengths. The binding of both ions was found to be quantitatively enhanced as the percentage of negative phospholipid increases in the membrane. Also, a decrease in ion binding was found to occur as the concentration of monovalent salt was increased (0.045-0.345 M). If electrostatic effects were ignored, the experimental data showed biphasic behavior in Scatchard plots. When electrostatic effects were taken into account by means of the Gouy-Chapman theory, the same data yielded linear Scatchard plots that were described by a simple partition equilibrium of the hydrophobic ion into the lipid-water interface. We demonstrate that the effective interfacial charge, nu, of the ion is a determinant factor to obtain a unique value of the intrinsic (hydrophobic) binding constant independently of the surface charge density of the lipid membrane.

Stability of an ion channel in lipid bilayers: implicit solvent model calculations with gramicidin

Gramicidin is a helical peptide, 15 residues in length, which dimerizes to form ion-conducting channels in lipid bilayers. Here we report calculations of its free energy of transfer from the aqueous phase into bilayers of different widths. The electrostatic and nonpolar contributions to the desolvation free energy were calculated using implicit solvent models, in which gramicidin was described in atomic detail and the hydrocarbon region of the membrane was described as a slab of hydrophobic medium embedded in water. The free energy penalties from the lipid perturbation and membrane deformation effects, and the entropy loss associated with gramicidin immobilization in the bilayer, were estimated from a statistical thermodynamic model of the bilayer. The calculations were carried out using two classes of experimentally observed conformations: a head-to-head dimer of two single-stranded (SS) beta-helices and a double-stranded (DS) intertwined double helix. The calculations showed that gramicidin is likely to partition into the bilayer in all of these conformations. However, the SS conformation was found to be significantly more stable than the DS in the bilayer, in agreement with most of the experimental data. We tested numerous transmembrane and surface orientations of gramicidin in bilayers of various widths. Our calculations indicate that the most favorable orientation is transmembrane, which is indeed to be expected from a channel-forming peptide. The calculations demonstrate that gramicidin insertion into the membrane is likely to involve a significant deformation of the bilayer to match the hydrophobic width of the peptide (22 A), again in good agreement with experimental data. Interestingly, deformation of the bilayer was induced by all of the gramicidin conformations.