- Research article
- Open Access
- Open Peer Review
Hydrocarbon molar water solubility predicts NMDA vs. GABAA receptor modulation
BMC Pharmacology and Toxicologyvolume 15, Article number: 62 (2014)
Many anesthetics modulate 3-transmembrane (such as NMDA) and 4-transmembrane (such as GABAA) receptors. Clinical and experimental anesthetics exhibiting receptor family specificity often have low water solubility. We hypothesized that the molar water solubility of a hydrocarbon could be used to predict receptor modulation in vitro.
GABAA (α1β2γ2s) or NMDA (NR1/NR2A) receptors were expressed in oocytes and studied using standard two-electrode voltage clamp techniques. Hydrocarbons from 14 different organic functional groups were studied at saturated concentrations, and compounds within each group differed only by the carbon number at the ω-position or within a saturated ring. An effect on GABAA or NMDA receptors was defined as a 10% or greater reversible current change from baseline that was statistically different from zero.
Hydrocarbon moieties potentiated GABAA and inhibited NMDA receptor currents with at least some members from each functional group modulating both receptor types. A water solubility cut-off for NMDA receptors occurred at 1.1 mM with a 95% CI = 0.45 to 2.8 mM. NMDA receptor cut-off effects were not well correlated with hydrocarbon chain length or molecular volume. No cut-off was observed for GABAA receptors within the solubility range of hydrocarbons studied.
Hydrocarbon modulation of NMDA receptor function exhibits a molar water solubility cut-off. Differences between unrelated receptor cut-off values suggest that the number, affinity, or efficacy of protein-hydrocarbon interactions at these sites likely differ.
Inhaled anesthetics interact with putative cell receptor targets in a manner uncharacteristic of most other pharmacologic agents. They exhibit immobilization efficacy in all animals, both vertebrates and invertebrates, and even in protozoa and plants. The inhaled agents range in diversity from single elements to diatomic molecules to complex hydrocarbons and share no conserved size, shape, or functional groups. Most agents also modulate multiple phylogenetically unrelated cell targets and generally inhibit excitatory channels or receptors and potentiate inhibitory ones, complicating the identification of a drug-receptor structural motif predictive of anesthetic molecular action .
Attempts to determine how disparate drugs act on so many unrelated receptors—and to define those targets essential to inhaled anesthetic actions—have thus far proved illusory. Consequently, in the present study, we posed the opposite question: Why are some anesthetics unable to modulate certain anesthetic-sensitive ion channels or receptors, even at supra-pharmacologic concentrations? To answer this question, we restricted consideration to two unrelated ion channels in particular: the N-methyl-d-aspartate (NMDA) receptor, a member of the 3-transmembrane (TM3) ion channel family, and the γ-aminobutyric acid type A (GABAA) receptor, a member of the 4-transmembrane (TM4) ion channel family. Strong evidence supports the role of both receptors in mediating various endpoints of general anesthesia, including immobility and amnesia [2–6]. Within the NMDA and GABAA channel proteins, hydrophilic or amphipathic cavities or pockets have been postulated near the solvent interface and within transmembrane segments of several subunits [7, 8]. Presumably, such cavities might conceivably contain solvent molecules or nothing at all, though it should be thermodynamically less favorable to maintain a vacuum within a pocket that is sufficiently large and polar enough to accommodate water.
The binding affinity of water within each pocket on an ion channel can be expressed as a standard dissociation constant, and the ability of a drug to displace water and either completely or partially fill the pocket can be expressed as a standard equilibrium constant. Since structures of phylogenetically unrelated channels are different, the hydrophilicity of their respective pockets should not be conserved. Therefore, dissociation constants are expected to be different for water bound within pockets of different channels. Consequently, the minimum concentration of drug necessary to displace water within a pocket capable of inducing a functional conformational change should be different between unrelated channels. If the drug cannot reach a concentration within the aqueous phase sufficient to displace water within the protein modulatory pocket, then that drug should exert no pharmacologic effect on the channel, even when delivered at a saturated aqueous phase concentration. The physical property that defines this maximum (saturated) concentration of a drug in the aqueous phase is the molar water solubility.
The effects of a number of conventional (modern and historic) and experimental anesthetics on NMDA and GABAA receptors have been studied previously. When plotted in order of their calculated molar water solubilities (Figure 1), there is an abrupt change in the ability of compounds to modulate NMDA receptors. Drugs with water solubility greater than about 1 mM all modulate NMDA and GABAA currents when administered in sufficiently large concentrations. However, drugs with lower molar water solubility modulate only GABAA or neither receptor type. These data suggest that molar water solubility may predict an NMDA receptor cut-off effect at values less than 1 mM and are consistent with non-specific interactions in which water displaced from modulatory pockets on the NMDA receptor contribute, at least in part, to molecular mechanisms of action. Compounds unable to reach approximately 1 mM aqueous concentration may thus be unable to competitively displace water within any of the NMDA critical modulatory sites.
We hypothesized that hydrocarbon molar water solubility below this critical cut-off value is associated with loss of NMDA receptor modulation, but not GABAA receptor modulation. As tests agents, we considered 14 different series of normal hydrocarbon chains or rings, where molecules within each series differed only by the length of the Ω chain or substitutent alkyl groups (Table 1). Because different functional groups differ in hydrophilicity, the length and volume of the hydrocarbon chain needed to achieve a given molar water solubility will likewise differ. Thus, if the hypothesis is correct, then observed cut-off responses should correlate with hydrocarbon molar water solubility rather than hydrocarbon size.
Although generally advocated to study anesthetic effects at clinically-relevant concentrations, we proposed to test this hypothesis at saturated aqueous drug concentrations for 3 reasons. First, because an aqueous binding site is postulated, it is important to study drug effects at an interfacial concentration that can be directly related to a bulk aqueous concentration. However, anesthetics do not distribute equally through the lipid bilayer. Halothane shows a preference for the phospholipid headgroup interface . Xenon atoms prefer regions at the lipid-water interface and the central region of the bilayer . The anesthetics cyclopropane, nitrous oxide, desflurane, isoflurane, and 1,1,2-trifluoroethane all preferentially concentrate at the interface between water and hexane, but the nonimmobilizer perfluoroethane does not exhibit a hydrophilic-hydrophobic interfacial maxima . Nonetheless, even without knowing the membrane distribution profile of an anesthetic, the interfacial concentration can be assumed maximal at a saturating aqueous phase concentration at equilibrium; thus drug responses are compared at their respective relative maximum interfacial concentrations.
Second, different anesthetic endpoints are achieved at different drug concentrations. Thus, a drug could exhibit relative receptor specificity; that is to say, a drug may act preferentially at one receptor to achieve one endpoint—such as amnesia—but act on additional receptors when administered at higher concentrations to achieve other endpoints—such as immobility. Failure to modulate a receptor when a drug is delivered at a saturated concentration implies a null receptor response at lower drug concentrations and for any therapeutic endpoint.
Third, the drug concentrations that produce different anesthetic endpoints—amnesia, unconsciousness and immobility—are unknown for many experimental compounds. Some anesthetic endpoints may not even be achievable with some compounds, such as the nonimmobilizers. Absolute receptor specificity means that there is also relative specificity for any drug effect produced. Accordingly, comparisons of receptor effects using saturated concentrations in vitro are not prepositioned upon knowledge of in vivo drug anesthetic effects.
Hence, this study aimed to test whether NMDA versus GABAA receptor modulation was correlated with calculated molar water solubility “cut-off” values for diverse series of hydrocarbon functional groups.
Oocyte collection and receptor expression
An ovary from tricaine-anesthetized Xenopus laevis frogs was surgically removed using a protocol approved by the Institutional Animal Care and Use Committee at the University of California, Davis (Protocol #12030). Following manual disruption of the ovarian lobule septae, the ovary was incubated in 0.2% Type I collagenase (Worthington Biochemical, Lakewood, NJ) to defolliculate oocytes which were washed and stored in fresh and filtered modified Barth’s solution composed of 88 mM NaCl, 1 mM KCl, 2.4 mM NaHCO3, 20 mM HEPES, 0.82 mM MgSO4, 0.33 mM Ca(NO3)2, 0.41 mM CaCl2, 5 mM sodium pyruvate, gentamycin, penicillin, streptomycin, and corrected to pH = 7.4. All salts and antibiotics were A.C.S. grade (Fisher Scientific, Pittsburgh, PA).
Clones used were provided as a gift from Dr. R.A. Harris (University of Texas, Austin) and were sequenced and compared to references in the National Center for Biotechnology Information database to confirm the identity of each gene. GABAA receptors were expressed using clones for the human GABAA α1 and the rat GABAA β2 and γ2s subunits in pCIS-II vectors. Approximately 0.25-1 ng total plasmid mixture containing either α1, β2, or γ2 genes in a respective ratio of 1:1:10 was injected intranuclearly through the oocyte animal pole and studied 2–4 days later. These plasmid ratios ensured incorporation of the γ subunit into expressed receptors, as confirmed via receptor potentiation to 10 μM chlordiazepoxide or insensitivity to 10 μM zinc chloride during co-application with GABA. In separate oocytes, glutamate receptors were expressed using rat NMDA NR1 clones in a pCDNA3 vector and rat NMDA NR2A clones in a Bluescript vector. RNA encoding each subunit, prepared using a commercial transcription kit (T7 mMessage mMachine, Ambion, Austin, TX), was mixed in a 1:1 ratio, and 1–10 ng of total RNA was injected into oocytes and studied 1–2 days later. Oocytes injected with similar volumes of water served as controls.
GABAA receptor electrophysiology studies
Oocytes were studied in a 250 μL linear-flow perfusion chamber with solutions administered by syringe pump at 1.5 ml/min with gastight glass syringes and Teflon tubing. Oocyte GABAA currents were studied using standard two-electrode voltage clamping techniques at a holding potential of -80 mV using a 250 μL channel linear-flow perfusion chamber with solutions administered by syringe pump at 1.5 mL/min.
Frog Ringer’s (FR) solution composed of 115 mM NaCl, 2.5 mM KCl, 1.8 mM CaCl2, and 10 mM HEPES prepared in 18.2 MΩ H2O and filtered and adjusted to pH = 7.4 was used to perfuse oocytes. Agonist solutions also contained 20-to-40 μM, equal to an EC10–20, of 4-aminobutanoic acid (FR-GABA) [30–32]. After FR perfusion for 5 min, oocytes were exposed to 30 sec FR-GABA followed by another 5 min FR washout; this was repeated until stable GABAA-elicited peaks were obtained. Next, FR containing a saturated solution of the study drug (Table 1)—or for gaseous study compounds a vapor pressure equal to 90% of barometric pressure with balance oxygen—was used to perfuse the oocyte chamber for 2 min followed by perfusion with a FR-GABA solution containing the identical drug concentration for 30 sec. FR was next perfused for 5 min to allow drug washout, and oocytes were finally perfused with FR-GABA for 30 sec to confirm return of currents to within 10% of the initial baseline response.
NMDA receptor electrophysiology studies
Methods for measurement of whole-cell NMDA receptor currents have been described [33, 34]. Briefly, baseline perfusion solutions were the same as for GABAA with the substitution of equimolar BaCl2 for calcium salts and the addition of 0.1 mM EGTA; this constituted barium frog Ringer’s solution (BaFR). Agonist solutions for NMDA studies also contained 0.1 mM glutamate (E) and 0.01 mM glycine (G) to constitute a BaFREG solution that produced a NMDA receptor current ≥ EC99.
The syringe pump and perfusion chamber apparatus as well as the clamp holding potential and baseline-agonist exposure protocols were identical to that described for the GABAA studies. The same test compounds, concentrations, and preparative methods were used in NMDA voltage clamp studies as in the GABAA voltage clamp studies (Table 1).
Response calculations and data analysis
Modulating drug responses were calculated as the percent of the control (baseline) peak as follows: where ID and IB were the peak currents measured during agonist + drug and agonist baseline perfusions, respectively. When present, direct receptor activation by a drug was similarly calculated as a percent of the agonist response. Average current responses for each drug and channel were described by mean ± SD. A lack of receptor response (cut-off) was defined as a <10% change from baseline current that was statistically indistinguishable from zero using a two-tailed Student t-test. Hence, drug responses ≥110% of the baseline peak showed potentiation of receptor function, and drug responses ≤90% of the baseline peak showed inhibition of receptor function.
The log10 of the calculated solubility (log10S) for compounds immediately below and above the cut-off for each hydrocarbon functional group were used to determine the receptor cut-off. For each hydrocarbon, there was a “grey area” of indeterminate solubility effect (Figure 2) between sequentially increasing hydrocarbon chain lengths. Mean solubility cut-offs were calculated as the average log10S for the least soluble compound that modulated receptor function and the most soluble neighboring compound for which no effect was observed. From this result, a 95% confidence interval for log10S was calculated for receptor solubility cut-offs.
Hydrocarbon effects on NMDA and GABAA receptors are summarized in Table 2, and sample recordings are presented in Figure 2. All of the compounds tested positively modulated GABAA receptor function, and a few of the 5-to-6 carbon compounds caused mild direct GABAA receptor activation, particularly the 1-fluoroalkanes and thiols. Mild direct receptor activation also occurred with dibutylether. With the exception of the aldehydes, alkynes, and cycloalkanes, GABAA receptor inhibition tended to decrease with increasing hydrocarbon chain length. No water solubility cut-off effect was observed for GABAA receptors for the compounds tested.
In contrast, NMDA receptors currents were decreased by the shorter hydrocarbons within each functional group (Table 1), but lengthening the hydrocarbon chain eventually produced a null response—a cut-off effect. No direct hydrocarbon effects on NMDA receptor function were detected in the absence of glutamate and glycine agonist.The cut-off effect for NMDA receptor current modulation was associated with a hydrocarbon water solubility of 1.1 mM with a 95% confidence interval between 0.45 mM and 2.8 mM (Figure 3). More soluble hydrocarbons consistently inhibited NMDA receptor currents when applied at saturated aqueous concentrations, and hydrocarbons below this range had no appreciable effect on NMDA receptor function. Moveover, during the course of the study, water solubility was sufficiently predictive of an NMDA receptor cut-off so as to require identify and testing of only single pair of compounds bracketing this critical solubility value, as occurred with the alkenes, amines, cyclic hydrocarbons, and sulfur-containing compounds.Increasing hydrocarbon chain length decreases water solubility, but also increases molecular size. However, when graphed as a function of either carbon number (Figure 4) or molecular volume (Figure 5), the observed NMDA receptor cut-off effects show no consistent pattern. For example, the n-alkanes, 1-alkenes, and 1-alkynes show progressive lengthening of the hydrocarbon chain cut-off, presumably as a result of the increasing aqueous solubility conferred by the double and triple carbon bonds, respectively. There was also tremendous variation in molecular size of compounds exhibiting NMDA receptor cut-off. Alkanes exhibited NMDA receptor cut-off between butane and pentane, respectively 4 and 5 carbons in length, whereas the primary amines exhibited cut-off between 1-octadecanamine and 1-eicosanamine, respectively 18 and 20 carbons in length. As expected, the molecular volume of these compounds associated with NMDA receptor cut-off is also quite different, with the primary amine being over 3 times larger than the alkane.
NMDA receptor modulation is associated with an approximate 1.1 mM water solubility cut-off (Figure 3). In contrast, GABAA receptors potentiated all study compounds; this may be because a GABAA cut-off occurs at a lower water solubility value. Sequentially increasing hydrocarbon length to find a receptor cut-off effect potentially introduces confounding factors of carbon number and molecular volume that could in turn be responsible for the cut-off effect [35–38]. However, an aggregate comparison of cut-off values for all hydrocarbon groups as a function of carbon number (Figure 4) or molecular volume (Figure 5) shows no discernible pattern, suggesting that these physical properties are unlikely the primary limiting factors for drug-receptor modulation.
The correlation between molar water solubility and the NMDA receptor cut-off suggests hydrocarbons compete with water for amphipathic binding pockets within anesthetic-sensitive ion channels. Most inhaled anesthetics exhibit low-affinity binding on receptors as evidenced by generally large median effective concentrations for anesthesia—in the 230–290 μM range for isoflurane and halothane —as compared to agents that exert narcotic effects via a singular or primary molecular targets. These specific interactions—exemplified by ketamine antagonism of NMDA receptors , etomidate agonism of GABAA receptors , dexmedetomidine agonism of α2-adrenoreceptors , and morphine agonism of μ-opioid receptors —typically require only a few μM or less of drug and are consistent with high affinity interactions to resulting in induced fit binding. Instead, volatile anesthetics likely bind to pre-existing pockets and surfaces on or within the protein . Amphipathic pockets likely contain water molecules; when these are displaced by amphipathic drugs, fewer strong hydrophilic interactions and more hydrophobic interactions are possible with amino acid side chains in the cavity. We propose such nonspecific binding causes a change in pocket shape and, in consequence, the larger three-dimensional protein structure that affects channel gating or conductance.
Hydrocarbons act as hydrogen bond donors—or in the case of electrophiles, as hydrogen bond acceptors—with amino acid residues on anesthetic-sensitive receptors, resulting in displacement of water molecules from these binding pockets and alteration of protein function [44–46]. These low energy anesthetic-protein interactions are postulated to be enthalpically favorable since the displaced water molecules should be better able to hydrogen bond with like molecules in the bulk solvent rather than with amino acids [44, 46]. Halothane and isoflurane both have been shown to bind in water accessible pockets formed between α-helices in δ-subunits of the nicotinic acetylcholine receptor , a member of the 4-transmembrane receptor superfamily that includes the GABAA receptor. Models of nicotinic acetylcholine receptors and GABAA receptors further suggest that endogenous agonist or anesthetic binding might increase water accumulation in hydrophilic pockets and increase the number and accessibility of hydrophilic sites that are important for channel gating [48, 49]. However, molecules that are insufficiently water soluble may not be able to displace enough water molecules at enough critical sites in order to modulate channel function.
NMDA receptor modulation by inhaled anesthetics such as isoflurane, xenon, and carbon dioxide occurs—at least in part—at hydrophilic agonist binding sites [33, 50]. However, data from the present study show that less hydrophilic drugs—those with lower molar water solubilities—are still able to modulate GABAA receptor currents even when NMDA receptor efficacy is lost. Since these receptors belong to different and phylogenetically distinct superfamilies, it seems likely that either the number of displaced water molecules required to effect modulation, the relative affinities of the hydrocarbon versus water molecule for a critical hydrophilic protein pocket, and/or the number of hydrophilic sites necessary for allosteric modulation should also be different between proteins. Simply stated, there is no reason to suppose that unrelated channels are likely to have conserved affinity constants for water within protein cavities capable of inducing an effect on ion conductance. As evidence, GABA receptors currents, but not those of NMDA receptors, can be modulated by compounds with molar aqueous solubilities much less than 1.1 mM. Presumably, the water dissociation constant in the NMDA receptor binding site is lower than that in the analogous GABAA receptor binding site.
The locations within NMDA receptor subunits responsible for anesthetic-mediated modulation of protein function remain unresolved. However, it is likely that volatile anesthetics bind at multiple extracellular and transmembrane interfacial sites and cavities, as has been observed with nicotinic acetylcholine receptors [47, 51], Gloeobacter ligand-gated ion channels , and voltage-gated sodium channels . Whether an anesthetic positively or negatively modulates an ion channel may be a function of competing drug interactions at different sites within the protein . The protein cavities must to some extent be hydrophilic so that they are occupied by water in the native state. However, if these sites are too hydrophilic, then the energy necessary for a low-affinity drug to displace the water molecules becomes too great, requiring either very high drug concentrations that alter cell effects due to changes in serum osmolality or highly polar or charged drugs which, by virtue of these properties, are impermeable to the cell membrane and therefore cannot access critical transmembrane modulatory sites.
Although the present study assessed the association between water solubility and drug efficacy on anesthetic-sensitive ion channels, a relationship between solubility and potency may exist as well. For example, site-directed mutagenesis of the Ala825 residue on the M4 domain of the NR2A subunit, a region at which alcohols bind and can negatively modulate NMDA receptor currents [55, 56], affects potency of the slightly water-soluble hypnotic agent tribromoethanol. In fact, the hydrophobicity of the substituted amino acid negatively correlates with tribromoethanol potency as an NMDA receptor inhibitor . Interestingly, except when replaced by the extremely hydrophobic tryptophan residue, which may change binding or access of water itself to the cavity, mutagenesis at this same site had no effect on the infinitely water-soluble anesthetic ethanol. Perhaps the small and highly polar ethanol molecule mimics many of the intramolecular forces, such as hydrogen bonding, of the water molecules it displaces within the binding site, thus extreme changes within the pocket are required to affect potency.
Water and lipid solubility also affect drug potency for wild-type receptors. Increasing chain length of straight-chain alcohols or diols increases hydrophobicity and is initially associated with increased inhibitory potency of NMDA receptors [58, 59]. Similarly, the magnitude of GABAA receptor positive modulation in the present study tended to increase as hydrocarbon chains lengthened within any functional group. This is consistent with the Meyer-Overton prediction of increased anesthetic potency as a function of increasing hydrophobicity [60, 61]. The probability of a hydrocarbon passively entering a lipid cell membrane is parabolically related to hydrocarbon hydrophobicity; more hydrophobic molecules—up to a point—may be able to more easily access transmembrane modulatory sites . However, drug solubility in the lipid membrane should not affect drug concentration in the amphipathic protein pocket, since the net energy involved in moving into and out of the lipid membrane are equal and opposite. The change in state—in this case, drug diffusion in a reversible path from perfusate to lipid membrane to protein pocket—is thermodynamically defined only by the initial state (free energy in perfusate) and final state (free energy of receptor binding), and the change in free energy is independent of the membrane path between these two states . However, increasingly hydrophobic molecules should differ more in their intermolecular interactions with surrounding amino acid side chains compared to the water molecules they displaced. Therefore, if they can successfully access this amphipathic pocket, increasingly hydrophobic molecules may be capable of producing larger conformational changes in the protein and greater modulation of protein function. However, as molecules become even more hydrophobic and water solubility falls below the cut-off value, there are simply insufficient molecules in the aqueous phase to successfully compete with water at hydrophilic modulation or transduction sites on a receptor to alter its function. There is progressive loss of modulatory efficacy at sites with higher cut-offs that reduces the maximum drug-effect magnitude. Finally, when the drug water solubility becomes such that it is insufficient in concentration to out-compete the water in the lowest-cut-off site, the drug effect is reduced to zero.
Likewise, in the whole animal, this plausibly explains why transitional compounds and nonimmobilizers predicted by the Meyer-Overton correlation to produce anesthesia either have lower than expected potency or lack anesthetic efficacy altogether. As with the NMDA cut-off hydrocarbons presented in the present study, transitional compounds and nonimmobilizers all share a common property of low aqueous solubility . Nonimmobilizers such as 1,2-dichlorohexafluorocyclobutane fail to depress GABAA-dependent pyramidal cells  or NMDA-dependent CA1 neurons  in the hippocampus, and likely lack these effects elsewhere in the central nervous system. With decreasing water solubility, there is differential loss of receptor effects—such as occurred with NMDA receptors versus GABAA receptors in the present study. The anesthetic cut-off effect in whole animal models correlates with agent water solubility, and might be explained by the loss of one or more anesthetic-receptor contributions to central nervous system depression. Conversely, receptor molar water solubility cut-off values may be used to define those ion channels that are not essential for volatile anesthetic potency. Inhaled agents likely act via low affinity interactions with multiple cell receptors and ion channels to decrease neuronal excitability in the brain and spinal cord, but a loss or inadequate contribution from certain targets—perhaps GABAA or glycine receptors—as water solubility decreases may render a drug a nonimmobilizer. Additionally, agents having a water solubility below the cut-off value for some anesthetic-sensitive receptors may also produce undesirable pharmacologic properties, such as seizures following the loss of GABAA receptor modulation . In contrast, NMDA receptors can contribute to immobilizing actions of conventional volatile anesthetics , but they are not as a general principle essential for inhaled anesthetic action since an agent like pentane does not modulate NMDA receptors—even at a saturated aqueous concentration (Table 2)—yet has a measurable minimum alveolar concentration [68, 69].
As shown in Figure 3, the different hydrocarbon series exhibit small variability about the 1.1 mM cut-off. Some variability is due simply to the lack of compounds of intermediate solubility within a functional group series. For example, pentanethiol inhibited NMDA receptors, whereas the 1-carbon longer hexanethiol did not (Table 2). This pre-cut-off thiol is nearly 3-times more soluble in water than its post-cut-off cognate; yet it is not possible to obtain a more narrowly defined cut-off delineation for 1-thiols. Even larger variability was observed with the dialkylbenzene series, to which 1 additional carbon was added to each 1- and 3-alkyl group. The solubility ratio between the NMDA antagonist 1,3-dimethylbenzene and its cut-off cognate 1,3-diethylbenzene is more than 18 (Table 2).
Variability about the molar water solubility NMDA receptor cut-off may also have arisen from the use of calculated, rather than measured, values for hydrocarbon molar water solubility. Aqueous solubility is difficult to measure accurately, particularly for poorly soluble substances. Calculated solubilities are more accurate for small uncharged compounds, but still can have an absolute error within 1 log unit . However, even predicted values for nonpolar n-alkanes may show large deviations from experimental data as the hydrocarbon chain length increases .
Furthermore, the molar solubility values used in the present study were calculated for pure water at 25°C and at pH = 7.0. These were not the conditions under which drug-receptor effects were studied. Ringer’s oocyte perfusates contained buffers and physiologic concentrations of sodium, potassium, and chloride resulting in a 250 mOsm solution. The solubility of haloether and haloalkane anesthetic vapors vary inversely with osmolarity , as do the water-to-saline solubility ratio of benzenes, amines, and ketones . The presence of salts could have caused overestimation of aqueous solubility for some compounds when using values calculated for pure water. Likewise, solubility is also temperature-dependent. Studies were conducted at 22°C; solubility of gases in water should be greater than values calculated at 25°C. In contrast, most solutes used in the present study have negative enthalpy for dissolution , so solubility should be decreased at the lower ambient temperature. The reverse should occur for exothermic solutions, as predicted by the Le Chatelier principle. As for hydronium ion concentration, the solubility of most study compounds is trivially affected at pH values between 7-to-8. However, hydrocarbons containing an amine group have pKa values that are closer to physiologic pH, and the calculated aqueous solubility of 1-eicosanamine and 1-octadecanamine (Table 1) decreases by about 66% as pH increases from 7 to 8. Calculated molar water solubilities for the amines in this study were probably modestly overestimated at a physiologic pH equal to 7.4.
Only water solubility reliably predicted the NMDA receptor cut-off. Yet, molecular size and shape still likely influence this effect to some lesser degree. Most of the hydrocarbons examined in the present study had functional groups located on the 1- or 2-carbon position. However, the ethers were all 1,1′-oxybisalkanes; each member of the ether consisted of symmetrical 1-carbon additions to alkyl groups on either side of the oxygen atom (Table 1). Consequently, this electron-rich oxygen atom allowing hydrogen bonding with water molecules or amino acid residues with strong partial positive charges lies buried in the middle of the ether. For hydrocarbons with equivalent molar water solubilities, it may be more difficult for dialkyl ether to form hydrogen bonds in hydrophilic receptor pockets compared to a long primary amine (Table 1) that might more easily insert its nucleophilic terminus into the anesthetic-binding pocket while the long hydrophobic carbon chain remains in the lipid membrane. This may explain why ethers in this study appear to exhibit an NMDA cut-off that is slightly greater than hydrocarbons with other functional groups. Perhaps if a methyl-alkyl ether series were used instead of a dialkyl ether series, the apparent molar water solubility cut-off for this group would have been lower. Nonetheless, the cut-off variability is sufficiently small to allow a priori predictions of low-affinity hydrocarbon modulation of NMDA receptors. It remains to be shown whether other anesthetic-sensitive ion channels exhibit distinct cut-off effects that may also be predicted by a single physical property: molar water solubility.
Cut-off responses for NMDA receptor inhibition by diverse hydrocarbons occurs when drug molar water solubility is less than approximately 1.1 mM. However, hydrocarbons having lower molar water solubilities are still able to potentiate GABAA receptors. This finding supports a hypothesis that volatile compounds, such as inhaled anesthetics, access one or more amphipathic “pockets” within a protein, such as the NMDA receptor, to displace resident water molecules and induce a conformational change. Volatile anesthetics and other compounds with insufficient molar water solubility can never achieve sufficient concentration at the amphipathic pocket site to displace water and modulate protein function—even when such compounds are administered at a saturated aqueous phase concentration.
Brosnan RJ: Inhaled anesthetics in horses. Vet Clin North Am Equine Pract. 2013, 29 (1): 69-87. 10.1016/j.cveq.2012.11.006.
Chau PL: New insights into the molecular mechanisms of general anaesthetics. Br J Pharmacol. 2010, 161 (2): 288-307. 10.1111/j.1476-5381.2010.00891.x.
Stabernack C, Sonner JM, Laster M, Zhang Y, Xing Y, Sharma M, Eger EI: Spinal N-methyl-d-aspartate receptors may contribute to the immobilizing action of isoflurane. Anesth Analg. 2003, 96 (1): 102-107.
Rau V, Oh I, Liao M, Bodarky C, Fanselow MS, Homanics GE, Sonner JM, Eger EI: Gamma-aminobutyric acid type A receptor beta3 subunit forebrain-specific knockout mice are resistant to the amnestic effect of isoflurane. Anesth Analg. 2011, 113 (3): 500-504.
Rau V, Iyer SV, Oh I, Chandra D, Harrison N, Eger EI, Fanselow MS, Homanics GE, Sonner JM: Gamma-aminobutyric acid type A receptor alpha 4 subunit knockout mice are resistant to the amnestic effect of isoflurane. Anesth Analg. 2009, 109 (6): 1816-1822. 10.1213/ANE.0b013e3181bf6ae6.
Dutton RC, Laster MJ, Xing Y, Sonner JM, Raines DE, Solt K, Eger EI: Do N-methyl-D-aspartate receptors mediate the capacity of inhaled anesthetics to suppress the temporal summation that contributes to minimum alveolar concentration?. Anesth Analg. 2006, 102 (5): 1412-1418. 10.1213/01.ane.0000205759.67123.76.
Mayer ML: Glutamate receptors at atomic resolution. Nature. 2006, 440 (7083): 456-462. 10.1038/nature04709.
Ernst M, Bruckner S, Boresch S, Sieghart W: Comparative models of GABAA receptor extracellular and transmembrane domains: important insights in pharmacology and function. Mol Pharmacol. 2005, 68 (5): 1291-1300. 10.1124/mol.105.015982.
Nistri A, Berti C: Potentiating action of midazolam on GABA-mediated responses and its antagonism by Ro 14–7437 in the frog spinal cord. Neurosci Lett. 1983, 39 (2): 199-204. 10.1016/0304-3940(83)90077-0.
Macdonald R, Barker JL: Benzodiazepines specifically modulate GABA-mediated postsynaptic inhibition in cultured mammalian neurones. Nature. 1978, 271 (5645): 563-564. 10.1038/271563a0.
Dildy-Mayfield JE, Mihic SJ, Liu Y, Deitrich RA, Harris RA: Actions of long chain alcohols on GABAA and glutamate receptors: relation to in vivo effects. Br J Pharmacol. 1996, 118 (2): 378-384. 10.1111/j.1476-5381.1996.tb15413.x.
Flood P, Krasowski MD: Intravenous anesthetics differentially modulate ligand-gated ion channels. Anesthesiology. 2000, 92 (5): 1418-1425. 10.1097/00000542-200005000-00033.
Kendig JJ, Kodde A, Gibbs LM, Ionescu P, Eger EI: Correlates of anesthetic properties in isolated spinal cord: cyclobutanes. Eur J Pharmacol. 1994, 264 (3): 427-436. 10.1016/0014-2999(94)00499-4.
Jenkins A, Franks NP, Lieb WR: Effects of temperature and volatile anesthetics on GABA(A) receptors. Anesthesiology. 1999, 90 (2): 484-491. 10.1097/00000542-199902000-00024.
Krasowski MD, Harrison NL: The actions of ether, alcohol and alkane general anaesthetics on GABAA and glycine receptors and the effects of TM2 and TM3 mutations. Br J Pharmacol. 2000, 129 (4): 731-743. 10.1038/sj.bjp.0703087.
Hollmann MW, Liu HT, Hoenemann CW, Liu WH, Durieux ME: Modulation of NMDA receptor function by ketamine and magnesium. Part II: interactions with volatile anesthetics. Anesth Analg. 2001, 92 (5): 1182-1191.
Nishikawa K, Harrison NL: The actions of sevoflurane and desflurane on the gamma-aminobutyric acid receptor type A: effects of TM2 mutations in the alpha and beta subunits. Anesthesiology. 2003, 99 (3): 678-684. 10.1097/00000542-200309000-00024.
Yamakura T, Sakimura K, Shimoji K, Mishina M: Effects of propofol on various AMPA-, kainate- and NMDA-selective glutamate receptor channels expressed in Xenopus oocytes. Neurosci Lett. 1995, 188 (3): 187-190. 10.1016/0304-3940(95)11431-U.
Hales TG, Lambert JJ: The actions of propofol on inhibitory amino acid receptors of bovine adrenomedullary chromaffin cells and rodent central neurones. Br J Pharmacol. 1991, 104 (3): 619-628. 10.1111/j.1476-5381.1991.tb12479.x.
Yamakura T, Chavez-Noriega LE, Harris RA: Subunit-dependent inhibition of human neuronal nicotinic acetylcholine receptors and other ligand-gated ion channels by dissociative anesthetics ketamine and dizocilpine. Anesthesiology. 2000, 92 (4): 1144-1153. 10.1097/00000542-200004000-00033.
Yamakura T, Harris RA: Effects of gaseous anesthetics nitrous oxide and xenon on ligand-gated ion channels. Comparison with isoflurane and ethanol. Anesthesiology. 2000, 93 (4): 1095-1101. 10.1097/00000542-200010000-00034.
Ogata J, Shiraishi M, Namba T, Smothers CT, Woodward JJ, Harris RA: Effects of anesthetics on mutant N-methyl-D-aspartate receptors expressed in Xenopus oocytes. J Pharmacol Exp Ther. 2006, 318 (1): 434-443. 10.1124/jpet.106.101691.
Martin DC, Plagenhoef M, Abraham J, Dennison RL, Aronstam RS: Volatile anesthetics and glutamate activation of N-methyl-D-aspartate receptors. Biochem Pharmacol. 1995, 49 (6): 809-817. 10.1016/0006-2952(94)00519-R.
Wong EH, Kemp JA, Priestley T, Knight AR, Woodruff GN, Iversen LL: The anticonvulsant MK-801 is a potent N-methyl-D-aspartate antagonist. Proc Natl Acad Sci U S A. 1986, 83 (18): 7104-7108. 10.1073/pnas.83.18.7104.
Hara K, Eger EI, Laster MJ, Harris RA: Nonhalogenated alkanes cyclopropane and butane affect neurotransmitter-gated ion channel and G-protein-coupled receptors: differential actions on GABAA and glycine receptors. Anesthesiology. 2002, 97 (6): 1512-1520. 10.1097/00000542-200212000-00025.
Krasowski MD, Jenkins A, Flood P, Kung AY, Hopfinger AJ, Harrison NL: General anesthetic potencies of a series of propofol analogs correlate with potency for potentiation of gamma-aminobutyric acid (GABA) current at the GABA(A) receptor but not with lipid solubility. J Pharmacol Exp Ther. 2001, 297 (1): 338-351.
Vemparala S, Saiz L, Eckenhoff RG, Klein ML: Partitioning of anesthetics into a lipid bilayer and their interaction with membrane-bound peptide bundles. Biophys J. 2006, 91 (8): 2815-2825. 10.1529/biophysj.106.085324.
Stimson LM, Vattulainen I, Rog T, Karttunen M: Exploring the effect of xenon on biomembranes. Cell Mol Biol Lett. 2005, 10 (4): 563-569.
Pohorille A, Wilson MA, New MH, Chipot C: Concentrations of anesthetics across the water-membrane interface; the Meyer-Overton hypothesis revisited. Toxicol Lett. 1998, 100–101: 421-430.
Brosnan R, Gong D, Cotten J, Keshavaprasad B, Yost CS, Eger EI, Sonner JM: Chirality in anesthesia II: stereoselective modulation of ion channel function by secondary alcohol enantiomers. Anesth Analg. 2006, 103 (1): 86-91. 10.1213/01.ane.0000221437.87338.af.
Yang L, Milutinovic PS, Brosnan RJ, Eger EI, Sonner JM: The plasticizer di(2-ethylhexyl) phthalate modulates gamma-aminobutyric acid type A and glycine receptor function. Anesth Analg. 2007, 105 (2): 393-396. 10.1213/01.ane.0000267336.37735.d7.
Yang L, Zhao J, Milutinovic PS, Brosnan RJ, Eger EI, Sonner JM: Anesthetic properties of the ketone bodies beta-hydroxybutyric acid and acetone. Anesth Analg. 2007, 105 (3): 673-679. 10.1213/01.ane.0000278127.68312.dc.
Brosnan RJ, Pham TL: Does anesthetic additivity imply a similar molecular mechanism of anesthetic action at N-methyl-D-aspartate receptors?. Anesth Analg. 2011, 112 (3): 568-573. 10.1213/ANE.0b013e3182080599.
Brosnan RJ, Pham TL: Carbon dioxide negatively modulates N-methyl-D-aspartate receptors. Br J Anaesth. 2008, 101 (5): 673-679. 10.1093/bja/aen266.
Eger EI, Halsey MJ, Harris RA, Koblin DD, Pohorille A, Sewell JC, Sonner JM, Trudell JR: Hypothesis: volatile anesthetics produce immobility by acting on two sites approximately five carbon atoms apart. Anesth Analg. 1999, 88 (6): 1395-1400.
Jenkins A, Greenblatt EP, Faulkner HJ, Bertaccini E, Light A, Lin A, Andreasen A, Viner A, Trudell JR, Harrison NL: Evidence for a common binding cavity for three general anesthetics within the GABAA receptor. J Neurosci. 2001, 21 (6): RC136-
Wick MJ, Mihic SJ, Ueno S, Mascia MP, Trudell JR, Brozowski SJ, Ye Q, Harrison NL, Harris RA: Mutations of gamma-aminobutyric acid and glycine receptors change alcohol cutoff: evidence for an alcohol receptor?. Proc Natl Acad Sci U S A. 1998, 95 (11): 6504-6509. 10.1073/pnas.95.11.6504.
Eger EI, Laster MJ: The effect of rigidity, shape, unsaturation, and length on the anesthetic potency of hydrocarbons. Anesth Analg. 2001, 92 (6): 1477-1482.
Franks NP, Lieb WR: Selective actions of volatile general anaesthetics at molecular and cellular levels. Br J Anaesth. 1993, 71 (1): 65-76. 10.1093/bja/71.1.65.
Tomlin SL, Jenkins A, Lieb WR, Franks NP: Stereoselective effects of etomidate optical isomers on gamma-aminobutyric acid type A receptors and animals. Anesthesiology. 1998, 88 (3): 708-717. 10.1097/00000542-199803000-00022.
Tonner PH, Scholz J, Koch C, Schulte am Esch J: The anesthetic effect of dexmedetomidine does not adhere to the Meyer-Overton rule but is reversed by hydrostatic pressure. Anesth Analg. 1997, 84 (3): 618-622.
Pert CB, Snyder SH: Properties of opiate-receptor binding in rat brain. Proc Natl Acad Sci U S A. 1973, 70 (8): 2243-2247. 10.1073/pnas.70.8.2243.
Trudell JR, Bertaccini E: Molecular modelling of specific and non-specific anaesthetic interactions. Br J Anaesth. 2002, 89 (1): 32-40. 10.1093/bja/aef157.
Bertaccini EJ, Trudell JR, Franks NP: The common chemical motifs within anesthetic binding sites. Anesth Analg. 2007, 104 (2): 318-324. 10.1213/01.ane.0000253029.67331.8d.
Abraham MH, Lieb WR, Franks NP: Role of hydrogen bonding in general anesthesia. J Pharm Sci. 1991, 80 (8): 719-724. 10.1002/jps.2600800802.
Streiff JH, Jones KA: Volatile anesthetic binding to proteins is influenced by solvent and aliphatic residues. J Chem Inf Model. 2008, 48 (10): 2066-2073. 10.1021/ci800206a.
Chiara DC, Dangott LJ, Eckenhoff RG, Cohen JB: Identification of nicotinic acetylcholine receptor amino acids photolabeled by the volatile anesthetic halothane. Biochemistry. 2003, 42 (46): 13457-13467. 10.1021/bi0351561.
Williams DB, Akabas MH: Gamma-aminobutyric acid increases the water accessibility of M3 membrane-spanning segment residues in gamma-aminobutyric acid type A receptors. Biophys J. 1999, 77 (5): 2563-2574. 10.1016/S0006-3495(99)77091-8.
Willenbring D, Xu Y, Tang P: The role of structured water in mediating general anesthetic action on alpha4beta2 nAChR. Phys Chem Chem Phys. 2010, 12 (35): 10263-10269. 10.1039/c003573d.
Dickinson R, Peterson BK, Banks P, Simillis C, Martin JC, Valenzuela CA, Maze M, Franks NP: Competitive inhibition at the glycine site of the N-methyl-D-aspartate receptor by the anesthetics xenon and isoflurane: evidence from molecular modeling and electrophysiology. Anesthesiology. 2007, 107 (5): 756-767. 10.1097/01.anes.0000287061.77674.71.
Brannigan G, LeBard DN, Henin J, Eckenhoff RG, Klein ML: Multiple binding sites for the general anesthetic isoflurane identified in the nicotinic acetylcholine receptor transmembrane domain. Proc Natl Acad Sci U S A. 2010, 107 (32): 14122-14127. 10.1073/pnas.1008534107.
Chen Q, Cheng MH, Xu Y, Tang P: Anesthetic binding in a pentameric ligand-gated ion channel: GLIC. Biophys J. 2010, 99 (6): 1801-1809. 10.1016/j.bpj.2010.07.023.
Raju SG, Barber AF, LeBard DN, Klein ML, Carnevale V: Exploring volatile general anesthetic binding to a closed membrane-bound bacterial voltage-gated sodium channel via computation. PLoS Comput Biol. 2013, 9 (6): e1003090-10.1371/journal.pcbi.1003090.
Bromstrup T, Howard RJ, Trudell JR, Harris RA, Lindahl E: Inhibition versus potentiation of ligand-gated ion channels can be altered by a single mutation that moves ligands between intra- and intersubunit sites. Structure. 2013, 21 (8): 1307-1316. 10.1016/j.str.2013.06.018.
Xu M, Smothers CT, Trudell J, Woodward JJ: Ethanol inhibition of constitutively open N-methyl-D-aspartate receptors. J Pharmacol Exp Ther. 2012, 340 (1): 218-226. 10.1124/jpet.111.187179.
Ren H, Zhao Y, Dwyer DS, Peoples RW: Interactions among positions in the third and fourth membrane-associated domains at the intersubunit interface of the N-methyl-D-aspartate receptor forming sites of alcohol action. J Biol Chem. 2012, 287 (33): 27302-27312. 10.1074/jbc.M111.338921.
Salous AK, Ren H, Lamb KA, Hu XQ, Lipsky RH, Peoples RW: Differential actions of ethanol and trichloroethanol at sites in the M3 and M4 domains of the NMDA receptor GluN2A (NR2A) subunit. Br J Pharmacol. 2009, 158 (5): 1395-1404. 10.1111/j.1476-5381.2009.00397.x.
Peoples RW, Ren H: Inhibition of N-methyl-D-aspartate receptors by straight-chain diols: implications for the mechanism of the alcohol cutoff effect. Mol Pharmacol. 2002, 61 (1): 169-176. 10.1124/mol.61.1.169.
Peoples RW, Weight FF: Cutoff in potency implicates alcohol inhibition of N-methyl-D-aspartate receptors in alcohol intoxication. Proc Natl Acad Sci U S A. 1995, 92 (7): 2825-2829. 10.1073/pnas.92.7.2825.
Mihic SJ, McQuilkin SJ, Eger EI, Ionescu P, Harris RA: Potentiation of gamma-aminobutyric acid type A receptor-mediated chloride currents by novel halogenated compounds correlates with their abilities to induce general anesthesia. Mol Pharmacol. 1994, 46 (5): 851-857.
Horishita T, Eger EI, Harris RA: The effects of volatile aromatic anesthetics on voltage-gated Na + channels expressed in Xenopus oocytes. Anesth Analg. 2008, 107 (5): 1579-1586. 10.1213/ane.0b013e318184b966.
Hansch C, Fujita T: p-σ-π Analysis. A Method for the Correlation of Biological Activity and Chemical Structure. J Am Chem Soc. 1964, 86 (8): 1616-1626. 10.1021/ja01062a035.
Fermi E: Thermodynamics. 1937, New York, NY: Prentice-Hall
Eger EI: Mechanisms of Inhaled Anesthetic Action. The Pharmacology of Inhaled Anesthetics. Volume 1. Edited by: Eger EI2nd. 2002, IL, USA: Baxter Healthcare Corporation, 33-42.
Perouansky M, Banks MI, Pearce RA: The differential effects of the nonimmobilizer 1,2-dichlorohexafluorocyclobutane (F6, 2 N) and isoflurane on extrasynaptic gamma-aminobutyric acid A receptors. Anesth Analg. 2005, 100 (6): 1667-1673. 10.1213/01.ANE.0000150942.68328.66.
Taylor DM, Eger EI, Bickler PE: Halothane, but not the nonimmobilizers perfluoropentane and 1,2-dichlorohexafluorocyclobutane, depresses synaptic transmission in hippocampal CA1 neurons in rats. Anesth Analg. 1999, 89 (4): 1040-1045.
Raines DE: Anesthetic and nonanesthetic halogenated volatile compounds have dissimilar activities on nicotinic acetylcholine receptor desensitization kinetics. Anesthesiology. 1996, 84 (3): 663-671. 10.1097/00000542-199603000-00022.
Liu J, Laster MJ, Taheri S, Eger EI, Koblin DD, Halsey MJ: Is there a cutoff in anesthetic potency for the normal alkanes?. Anesth Analg. 1993, 77 (1): 12-18.
Taheri S, Laster MJ, Liu J, Eger EI, Halsey MJ, Koblin DD: Anesthesia by n-alkanes not consistent with the Meyer-Overton hypothesis: determinations of the solubilities of alkanes in saline and various lipids. Anesth Analg. 1993, 77 (1): 7-11.
Delaney JS: Predicting aqueous solubility from structure. Drug Discov Today. 2005, 10 (4): 289-295. 10.1016/S1359-6446(04)03365-3.
Ferguson AL, Debenedetti PG, Panagiotopoulos AZ: Solubility and molecular conformations of n-alkane chains in water. J Phys Chem. 2009, 113 (18): 6405-6414. 10.1021/jp811229q.
Lerman J, Willis MM, Gregory GA, Eger EI: Osmolarity determines the solubility of anesthetics in aqueous solutions at 37 degrees C. Anesthesiology. 1983, 59 (6): 554-558. 10.1097/00000542-198312000-00013.
Long FA, McDevit WF: Activity coefficients of nonelectrolyte solutes in aqueous salt solutions. Chem Rev. 1952, 51 (1): 119-169. 10.1021/cr60158a004.
Abraham MH: Free energies, enthalpies, and entropies of solution of gaseous nonpolar nonelectrolytes in water and nonaqueous solvents The hydrophobic effect. J Am Chem Soc. 1982, 104 (8): 2085-2094. 10.1021/ja00372a001.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/2050-6511/15/62/prepub
This entirety of this work was funded by the National Institutes of Health National Institutes of General Medical Sciences (Grant GM092821-02).
The authors declare that they have no competing interests.
RB: Developed hypothesis, designed experiments, assisted with experiments, analyzed data, and wrote the manuscript. TP: Conducted experiments and reviewed the manuscript. Both authors read and approved the final manuscript.