In Vitro Activation of Cytochrome P450 46A1 (CYP46A1) by Efavirenz-Related Compounds.

Cytochrome P450 46A1 (CYP46A1) is a central nervous system-specific enzyme, which catalyzes cholesterol 24-hydroxylation. Currently CYP46A1 is being evaluated in a clinical trial for activation by small doses of the anti-HIV drug efavirenz. Eight efavirenz-related compounds were investigated for CYP46A1 activation in vitro, induction of a CYP46A1 spectral response, spectral Kd values, interaction with the P450 allosteric sites, and a model of binding to the enzyme active site. We gained insight into structure-activity relationships of efavirenz for CYP46A1 activation and found that the investigated efavirenz primary metabolites are stronger and better activators of CYP46A1 than efavirenz. We also established that CYP46A1 is activated by racemates and that a conformational-selection mechanism is operative in CYP46A1. The results suggest structural modifications of efavirenz to further increase CYP46A1 activation without inhibition at high compound concentrations. It is possible that not only efavirenz but its metabolites activate CYP46A1 in vivo.

Introduction

CYP46A1 is an endoplasmic reticulum cytochrome P450 enzyme, which is expressed in the central nervous system, normally in the soma and dendrites of multiple neuron types.[2−4] n class="Gene">CYP46A1 converts cholesterol to 24-hydroxycholesterol (24HC), an oxysterol that rapidly diffuses across the blood–brain barrier to the systemic circulation and is delivered to the liver for further degradation to bile acids.[2,5,6] Initially CYP46A1 was recognized for its key roles in brain cholesterol elimination and turnover.[2,7,8] CYP46A1 was later discovered to be linked to other brain processes including memory and learning.[9−12] CYP46A1 is now a focus for intensive investigations as a potential target for treatment of different brain disorders with both enzyme inhibition and activation having therapeutic potential.[13] Animal studies suggest that increased CYP46A1 activity can ameliorate the manifestations of Alzheimer’s and Huntington’s diseases.[14−18] Conversely, a decreased CYP46A1 activity is hypothesized to reduce seizure frequency and is now being tested in two clinical trials in subjects with Dravet, Lennox-Gastaut, and Duplication 15q syndromes as well as CDKL5 deficiency disorder.[12,13] In addition to the brain, CYP46A1 is of importance for the retina and retinal blood vessels.[19]

This laboratory discovered that n class="Gene">CYP46A1 could be activated pharmacologically in mice by (S)-efavirenz (EFV), an anti-HIV drug.[20−22] Remarkably, the CYP46A1 activating EFV dose was very low (0.1 (mg/day)/kg body weight), at least 80 times lower than that (600 mg/day) given to HIV-positive individuals to keep their viral load low. When administered orally to normal (C57BL/6J) mice or a model of Alzheimer’s disease (5XFAD mice), this low EFV dose not only activated CYP46A1 in the brain but also elicited a compensatory increase in cholesterol biosynthesis and hence turnover.[20−22] As a result, the steady-state levels of cholesterol in the brain remained unchanged in C57BL/6J mice and were normalized in 5XFAD mice. In the latter, EFV treatment also led to behavioral improvements with the effect on the total brain amyloid β load being the treatment paradigm-specific.[21,22] However, when the EFV dose was increased to 0.22 (mg/day)/kg body weight, CYP46A1 activation was no longer observed and the higher 0.32 (mg/da)/kg body dose even inhibited the enzyme.[20] Thus, the therapeutic window for CYP46A1 activation by EFV in mice appears to be quite narrow.

In vitro studies utilizing purified recombinant CYP46A1 revealed that the concentn class="Species">ration dependence curves of EFV for the CYP46A1-mediated cholesterol 24-hydroxylation were bell-shaped with the enzyme activation at low compound concentrations and inhibition at high concentartions.[20] Subsequent mechanistic investigations suggested an allosteric site for EFV on CYP46A1,[20] which was then mapped experimentally to the cytosolic (proximal) surface of the P450 molecule.[23] This site was close to the redox partner binding site but away from the CYP46A1 active site, a banana-shaped tunnel extending from the lipid-embedded distal CYP46A1 surface to the interior of the protein molecule.[24] We proposed that at low concentrations EFV binds to the allosteric site and activates CYP46A1, whereas at high concentrations the drug likely binds to both the allosteric and active sites and inhibits CYP46A1 because of the competition with cholesterol for the active site.[20] EFV was also found to increase the rate of the fast CYP46A1 reduction and the total amount of reduced P450.[1] Furthermore, the analysis of the CYP46A1 crystal structure[24] suggested that EFV could displace the water molecule in the allosteric site and thereby alter the hydrogen bond network that connects the allosteric and active sites and determines in part the hydration of the active site.[23]

Currently, we are involved in an ongoing clinical trial of EFV in n class="Species">patients with mild cognitive impairment due to Alzheimer’s disease (NCT03706885, ClinicalTrials.gov), which tests the 50 and 200 mg/day EFV doses. In parallel, we have a goal of generating a compound with a broader therapeutic window for CYP46A1 activation than EFV and preferably more potent to lower a possibility of CYP46A1 inhibition and hence side effects while perhaps enhancing the potentially beneficial drug effects. Therefore we initiated the current study to gain insights into the role of some of the EFV functionalities. We tested two types of EFV-related compounds for the effects on purified CYP46A1 in vitro: some of the primary EFV metabolites produced during hepatic drug clearance[4,25−28] and EFV analogs with changes in the oxazinone and cyclopropyl rings.

Results

CYP46A1 Activation by (S)-EFV and Related Compounds

In mammals, (S)-EFV is mainly oxidized to n class="Chemical">(S)-8-hydroxyEFV (1), with the metabolism to (S)-7-hydroxyEFV (2) occurring to a lesser extent (Figure a). Both of these metabolites can be further hydroxylated to yield (S)-8,14-dihydroxyEFV (3) and (S)-7,8-dihydroxyEFV (4).[4,25−28] (S)-EFV also has secondary metabolites produced as a result of glucuronidation and/or sulfation of the EFV oxidation products.[25,29] Of the primary EFV metabolites, only (S)-1 is commercially available along with racemic (rac)-1, rac-2, and rac-3. Hence, all these compounds were tested in the enzyme assay using purified recombinant CYP46A1. Similar to CYP46A1 activation by (S)-EFV, the concentration dependence curves for (S)-1, rac-1, and rac-2 were bell-shaped, whereas that for rac-3 was essentially hyperbolic (Figure b). All of the tested metabolites had a maximal CYP46A1 activation higher than that with (S)-EFV (5.8-6.2-fold vs 4.1-fold, Table ), which was observed at the same 20 μM metabolite concentration. (S)-1 and rac-1 had almost identical concentration dependence curves and along with rac-2 did not inhibit CYP46A1 at the 60–100 μM concentrations that were inhibitory for (S)-EFV. Thus, EFV hydroxylations at C-7, -8, and -14 render the drug properties that we sought when the study was initiated, namely, an increase in CYP46A1 activation at lower metabolite concentrations and a decrease in enzyme inhibition at higher compound concentrations. We then evaluated commercially available EFV analogs (Figure c): rac-5, which has a C atom instead of N-1; (S)-6, lacking the keto group; rac-7, which has a methyl group instead of the keto group; rac-8, which has a C atom instead of O-3; rac-9, which has an ethyl group instead of a cyclopropyl ring. All of these compounds, except rac-5, had bell-shaped concentration dependence curves, and CYP46A1 activation was either similar (rac-8) or lower than that with (S)-EFV(rac-5, (S)-6, rac-7, and rac-9). P450 inhibition was decreased for all compounds at their higher concentrations, and rac-5 as well as (S)-6 did not inhibit the enzyme at the concentrations tested. Thus, lack of O-3 (rac-8) did not significantly impair CYP46A1 activation by the derivative but decreased CYP46A1 inhibition at higher compound concentrations. Similarly, CYP46A1 inhibition was decreased by the replacement or elimination of EFV functionalities at positions 1 (rac-5), 2 ((S)-6), and 14 (rac-9), but the enzyme activation was less than that by (S)-EFV.

Figure 1

EFV metabolites and analogs. (a, c) Chemical structures of EFV, its metabolites, and analogs studied in the present work. Structural divergence from EFV is highlighted in either red font or circles. EFV metabolism by the cytochrome P450 enzymes is also shown and taken with permission from ref (4). (b, d) Dependence of CYP46A1 activity in vitro on the concentration of (S)-EFV or an EFV-related compound. CYP46A1 activity (the rate of cholesterol 24-hydroxylation) is presented as nanomoles of 24-hydroxycholesterol (24HC) formed per nmole of CYP46A1 per min. The range of the tested compound concentrations is based on our previous study,[20] in which 100 μM (S)-EFV almost completely inhibited CYP46A1. The results are the mean ± SD of the measurements from the three independent experiments.

Table 1

Summary of the CYP46A1 Responses to the Studied EFV-Related Compounds

					interacts via the allosteric site for
compd	substrate-free CYP46A1a	cholesterol-bound CYP46A1a	maximal activation (fold)	maximal activation in the presence of l-Glu (fold)	EFV	Glu
(S)-EFV	type 2; H	reverse type 1; S	4.1	4.9	+	–
(S)-1	type 2; H	reverse type 1; S	6.2	7.4	+	+
rac-2	type 2; H	reverse type 1; S	5.8	6.8	+	+
rac-5	type 1; S	no spectral response	1.9	2.8	+	+
rac-3	type 2; H	reverse type 1; H	6.1	4.8	–	+
(S)-6	type 1; H	no spectral response	2.5	abolished	–	+
rac-7	type 1; H	no spectral response	2.6	abolished	–	+
rac-8	type 1; H	reverse type 1; S	3.9	4.6	+	–
rac-9	type I; H	reverse type 1; S	3.1	3.7	+	–

Presented is the spectral response type and fit: H, hyperbolic fit; S, sigmoidal fit.

EFV metabolites and analogs. (a, c) Chemical structures of n class="Chemical">EFV, its metabolites, and analogs studied in the present work. Structural divergence from EFV is highlighted in either red font or circles. EFV metabolism by the cytochrome P450 enzymes is also shown and taken with permission from ref (4). (b, d) Dependence of CYP46A1 activity in vitro on the concentration of (S)-EFV or an EFV-related compound. CYP46A1 activity (the rate of cholesterol 24-hydroxylation) is presented as nanomoles of 24-hydroxycholesterol (24HC) formed per nmole of CYP46A1 per min. The range of the tested compound concentrations is based on our previous study,[20] in which 100 μM (S)-EFV almost completely inhibited CYP46A1. The results are the mean ± SD of the measurements from the three independent experiments.

Spectral Responses in CYP46A1 Induced by (S)-EFV and Related Compounds

CYP46A1 is a hemoprotein and, as with most P450s, has a n class="Chemical">water molecule coordinating its heme iron in the sixth axial position when in the resting, substrate-free state.[24] When this water molecule is displaced and the heme iron is left uncoordinated, the CYP46A1 difference spectrum has a peak at ∼390 nm and a trough at ∼422 nm,[20] a so-called type I spectral response, elicited in many cases by the P450 substrates.[30−32] However, when the heme iron becomes coordinated with a nitrogen atom, a trough around ∼410 nm and a peak at ∼421–435 nm appear in the P450 difference spectrum,[33] a so-called type 2 spectral response, often produced by P450 inhibitors.[30−32] A third type of the spectral response for the P450 enzymes has also been described, reverse type I, which represents a mirror reflection of a type I difference spectrum (a trough at 380–390 nm and peak at 415–420 nm).[32] This response is believed to reflect the formation of an Fe–O bond, where the oxygen donor atom is from either a water molecule or the added compound.[32,34] Previously we reported that titrations with (S)-EFV induced a type 2 spectral response in substrate-free CYP46A1 (Figure ), which could be fit (ΔA) to a hyperbolic curve.[20] These data suggested that in the absense of cholesterol, EFV binds to the CYP46A1 active site and coordinates the heme iron with its N-1 atom. Similarly, type 2 spectral responses that fit to a hyperbolic curve were observed in substrate-free CYP46A1 upon the addition of (S)-1, rac-1, rac-2, and rac-3, an indication that the hydroxyl groups at C-7, C-8, and C-8,14 do not prevent EFV derivatives from binding to the CYP46A1 active site in an orientation that enables the formation of the N–Fe coordinate bond (Figure ). Yet the apparent Kd values of the hydroxylated EFV metabolites were 2–5 times higher than that of (S)-EFV, suggesting that the affinity of these compounds for CYP46A1 was decreased (Figure b). Conversely, rac-5, (S)-6, rac-7, rac-8, and rac-9 elicited a type 1 spectral response in substrate-free CYP46A1 consistent with absence of N-1 (rac-5) as well as the hydrogen bond-forming functionalities ((S)-6, rac-7, rac-8) that may help EFV bind in the CYP46A1 active site in the orientation suitable for the ligation of the heme iron. In addition, testing of rac-9 revealed the importance of the EFV cyclopropyl ring for the drug coordination of the heme iron. As with type 2 ligands, the binding of type 1 ligands can be fit to a hyperbolic curve except in the case of rac-5, whose curve was sigmoidal and may underlie the sigmoidal-like curve in the concentration dependence enzyme assay (Figure d).

Figure 2

Spectral responses induced in substrate-free CYP46A1 by EFV-related compounds. Data for (S)-EFV are shown in the center. ΔA is the amplitude of spectral changes in the P450 difference spectra (insets). The results are presented as the mean ± SD of the measurements in the three spectral titrations.

Spectral responses induced in substrate-free n class="Gene">CYP46A1 by EFV-related compounds. Data for (S)-EFV are shown in the center. ΔA is the amplitude of spectral changes in the P450 difference spectra (insets). The results are presented as the mean ± SD of the measurements in the three spectral titrations.

The allosteric response of CYP46A1 was only revealed when n class="Chemical">EFV was added to the P450 in the presence of the saturating (20 μM) cholesterol, i.e., substrate-bound enzyme.[20] The spectral titration curve became sigmoidal, and the spectral response was of a reverse type 1 nature (Figure ). On the basis of this and other data, we suggested EFV binding to an allosteric site and interpreted this response as an altered hydration of the CYP46A1 active site with the subsequent formation of the Fe–O coordinate bond.[20,23] Similarly, sigmoidal curves and a reverse type 1-like spectral response were observed in CYP46A1 upon addition of (S)-1, rac-1, rac-2, rac-8, and rac-9, consistent with comparable (≤25%) or higher CYP46A1 activation by these compounds at lower concentrations and their bell-shape concentration dependence curves (Figure b,d). rac-3 also elicited a reverse type 1-like spectral response in cholesterol-bound CYP46A1, but its spectral binding was of a hyperbolic nature. Since this metabolite essentially did not inhibit CYP46A1 in the concentration dependence enzyme assay (Figure b), we suggest that this response could reflect preferential EFV binding only to the allosteric site. Lastly, rac-5, (S)-6, and rac-7 did not seem to elicit a quantifiable spectral response in cholesterol-bound CYP46A1 while still activating the P450 (Figure d), although to a lower extent than EFV (up to 2.6-fold vs 4.1-fold). These data suggest that compound interactions with the allosteric site(s) could be spectrally silent and that this spectral silence could be due to absence of the hydrogen-forming N-1 (rac-5) as well as the 2-keto group ((S)-6 and rac-7) in the tested derivatives. CYP46A1 spectral titrations with different EFV-related compounds are summarized in Table .

Figure 3

Spectral responses induced in cholesterol-bound CYP46A1 by EFV-related compounds. Data for (S)-EFV are shown in the center. ΔA is the amplitude of spectral changes in the P450 difference spectra (insets). The results are presented as the mean ± SD of the measurements in the three spectral titrations.

Spectral responses induced in cholesterol-bound n class="Gene">CYP46A1 by EFV-related compounds. Data for (S)-EFV are shown in the center. ΔA is the amplitude of spectral changes in the P450 difference spectra (insets). The results are presented as the mean ± SD of the measurements in the three spectral titrations.

Mapping of the Allosteric Site(s) for Compound Binding

We first conducted in silico docking of the (S)- and (R)-isomers of 3, the best n class="Gene">CYP46A1 activator in vitro. Autodock, VINA, and Ledock were used, and none of the docking programs could fit these isomers to the allosteric site for (S)-EFV. The three programs did provide a consensus model for (S)-3 and (R)-3 binding to the so-called glutamate (Glu)-binding site, a region on the CYP46A1 proximal surface adjacent to the site EFV-binding (Figure a–c).[1] Glu activates CYP46A1 in vitro either alone or synergistically with (S)-EFV,[1] and we tested how CYP46A1 co-incubations with l-Glu and rac-3 or other compounds affected the enzyme activity under the conditions of the maximal P450 activation (Figure d). l-Glu did not synergistically activate CYP46A1 in the co-incubations with rac-3, (S)-6, or rac-7, with the P450 activation becoming either lower (rac-3) or even abolished ((S)-6 or rac-7) relative to the P450 activity in the presence of either activator. Conversely, co-incubations with l-Glu and (S)-1, rac-2, rac-5, rac-8, or rac-9 synergistically activated CYP46A1 in vitro. These data suggested that rac-3 and (S)-6 or rac-7 probably bind mainly to the allosteric site for Glu and compete with the neurotransmitter for this site, consistent with the computational predictions for (S)-3 and (R)-3. In contrast, other EFV-related compounds likely bind to both EFV- and Glu-binding sites. Accordingly we then tested co-incubations of (S)-EFV with rac-3, (S)-6, or rac-7 and found a synergistic CYP46A1 activation (Figure d). Thus we obtained support for the interaction of rac-3, (S)-6, and rac-7 mainly with the allosteric site for Glu.

Figure 4

Mapping of the CYP46A1 allosteric site for rac-3. (a–c) Computational models of (S)-EFV (black), l-Glu (yellow), (S)-3 (dark green), and (R)-3 (light green) binding to the CYP46A1 surface. Docking experiments were carried out on PDB code 2Q9G. Models for (S)-EFV and l-Glu are taken from ref (1), and those for (S)-3 and (R)-3 were generated in the present work. The heme group is in red. The nitrogen, oxygen, chlorine, and fluorine atoms are in blue, red, light green, and cyan, respectively. (d) Effect of l-Glu (100 μM), (S)-EFV (20 μM), or EFV analog (20 μM) on the maximal extent of CYP46A1 activation. CYP46A1 activity is presented as nanomoles of 24-hydroxycholesterol (24HC) formed per nmole of CYP46A1 per min. Control incubations with no activator, l-Glu only, and (S)-EFV only are shown as a white, yellow, and black bar, respectively. Incubations with an EFV-related compound only are shown as solid colored bars, and co-incubations with l-Glu and (S)-EFV are shown as dashed (yellow and black dashes, respectively) colored bars. The results are presented as the mean ± SD of the measurements from the three independent experiments.

Mapping of the CYP46A1 allosteric site for n class="Gene">rac-3. (a–c) Computational models of (S)-EFV (black), l-Glu (yellow), (S)-3 (dark green), and (R)-3 (light green) binding to the CYP46A1 surface. Docking experiments were carried out on PDB code 2Q9G. Models for (S)-EFV and l-Glu are taken from ref (1), and those for (S)-3 and (R)-3 were generated in the present work. The heme group is in red. The nitrogen, oxygen, chlorine, and fluorine atoms are in blue, red, light green, and cyan, respectively. (d) Effect of l-Glu (100 μM), (S)-EFV (20 μM), or EFV analog (20 μM) on the maximal extent of CYP46A1 activation. CYP46A1 activity is presented as nanomoles of 24-hydroxycholesterol (24HC) formed per nmole of CYP46A1 per min. Control incubations with no activator, l-Glu only, and (S)-EFV only are shown as a white, yellow, and black bar, respectively. Incubations with an EFV-related compound only are shown as solid colored bars, and co-incubations with l-Glu and (S)-EFV are shown as dashed (yellow and black dashes, respectively) colored bars. The results are presented as the mean ± SD of the measurements from the three independent experiments.

To further map compound binding, we used the CYP46A1 mutants genen class="Species">rated for this study as well as previously[1,23] and assessed how the replacement of the amino acid residues specific for EFV binding site (R415A and F405A) and Glu binding site (K358R/L/A, F416A, and Y427F/A) as well as the loop region forming the common border between the two sites (K422A and R424A) affect maximal CYP46A1 activation by different EFV metabolites and derivatives. The P450 activation by rac-3, (S)-6, or rac-7 was not significantly affected by the smaller sized replacements within the EFV-binding site (R415A and F405A) but was largely reduced or abolished by the smaller sized replacements within the Glu-binding site (K358A, F416A, and Y427A) (Figure a). The effect of the alanine substitutions of the border-forming Lys-422 and Arg-424 was differential, and activation pattern of the CYP46A1 mutants by rac-8 and rac-9 was different with only EFV-site but not the Glu-site smaller sized mutations mainly affecting the P450 activation, evidence for the preferential EFV-site binding (Figure b). Lastly, the activation of all of the smaller sized CYP46A1 mutants was altered by (S)-1, rac-2, and rac-5, suggesting that these compounds bind to both EFV- and Glu-binding sites. Thus, depending on the structure, EFV-related compounds could bind either to both or mainly one of the allosteric sites and activate CYP46A1 in vitro with the extent of the P450 activation being compound-specific but not the allosteric site-specific (Table ).

Figure 5

Mapping of the CYP46A1 allosteric site(s) for binding of EFV-related compounds. (a, b) Maximal activity of wild type (WT) CYP46A1 and the P450 mutants in the absence (white bars) or presence of different EFV-related compounds (20 μM, colored bars). CYP46A1 activity is presented as nanomoles of 24-hydroxycholesterol (24HC) formed per nmole of CYP46A1 per min. The results are presented as the mean ± SD of the measurements from the three independent experiments. (c–e). Some of the amino acid residues forming EFV binding site (c) and Glu binding site (d and e) as well as the border between the sites (c–e). The positions of (S)-EFV, (S)-3, and (R)-3 within these sites are also shown. Docking experiments were carried out on PDB code 2Q9G.

Mapping of the CYP46A1 allosteric site(s) for binding of n class="Chemical">EFV-related compounds. (a, b) Maximal activity of wild type (WT) CYP46A1 and the P450 mutants in the absence (white bars) or presence of different EFV-related compounds (20 μM, colored bars). CYP46A1 activity is presented as nanomoles of 24-hydroxycholesterol (24HC) formed per nmole of CYP46A1 per min. The results are presented as the mean ± SD of the measurements from the three independent experiments. (c–e). Some of the amino acid residues forming EFV binding site (c) and Glu binding site (d and e) as well as the border between the sites (c–e). The positions of (S)-EFV, (S)-3, and (R)-3 within these sites are also shown. Docking experiments were carried out on PDB code 2Q9G.

Studies of a Model of (S)-EFV and rac-3 Binding to the CYP46A1 Active Site

These studies were carried out to gain insight into why there was no decrease in CYP46A1 activation and ultimately enzyme inhibition by n class="Gene">rac-3 at higher metabolite concentrations as in the case with (S)-EFV (Figure b). Indeed, both (S)-EFV and rac-3 bind to the active site of substrate-free CYP46A1 and likely act as inhibitors by coordinating the P450 heme iron with the secondary amine group (Figure ). Further, only a 2-fold difference in the apparent spectral Kd of rac-3 and (S)-EFV for substrate-free CYP46A1 is not large enough to predict that there would be no competition between rac-3 and cholesterol for the CYP46A1 active site as likely in the case with (S)-EFV. We first determined the binding rates of (S)-EFV and rac-3 for substrate-free CYP46A1 at different compound concentrations as diagnostic for a model of compound binding: conformation-selection vs induced-fit[35,36] with the former discovered recently to be operative in several human P450 enzymes.[37−40] Increasing compound concentrations increased the magnitude of the spectral changes with both (S)-EFV and rac-3 (Figure a–d) but clearly decreased the rates of binding (kobserved) (Figure ). Notably, for rac-3, these rates were always lower (at least ∼7-fold) than those for (S)-EFV, and rac-3 showed essentially no spectral binding at 15–30 μM, the concentrations that still led to a measurable (S)-EFV binding. Similarly, when (S)-EFV and rac-3 were added to cholesterol-bound CYP46A1 (Figure e,f), the absorbance change for the type 2 spectral response was still observed for (S)-EFV at 50 μM but not rac-3 at 2 μM. A decrease in binding rates with increasing compound concentration indicates a conformation-selection model of binding.[35,36] Collectively, the data obtained suggested that different CYP46A1 conformations are necessary for (S)-EFV and rac-3 binding to the P450 active site and that much less enzyme molecules exist in solution in the conformation favorable for rac-3, thus preventing it from efficient competition with cholesterol for the CYP46A1 active site.

Figure 6

Binding of (S)-EFV and rac-3 to CYP46A1. The P450 (4 μM) was mixed with varying concentrations of a compound, and spectra were collected and analyzed. (a, b) Mixing CYP46A1 with low (2 μM) and (c, d) high (50 μM or 60 μM) compound concentrations. (e, f) Mixing of 4 μM P450 46A1·4 μM cholesterol complex with 50 μM or 2 μM compound. The fits shown are single exponentials. At least four individual traces were collected and averaged.

Figure 7

Dependence of rates of binding to CYP46A1 as a function of compound concentration. Traces of Δ(A437 – A416) (e.g., Figure ) were fit to single exponential plots. The error bars show the SD of the OLIS software fits. At least four separate reactions were averaged for each data point.

Binding of (S)-EFV and n class="Gene">rac-3 to CYP46A1. The P450 (4 μM) was mixed with varying concentrations of a compound, and spectra were collected and analyzed. (a, b) Mixing CYP46A1 with low (2 μM) and (c, d) high (50 μM or 60 μM) compound concentrations. (e, f) Mixing of 4 μM P450 46A1·4 μM cholesterol complex with 50 μM or 2 μM compound. The fits shown are single exponentials. At least four individual traces were collected and averaged.

Dependence of rates of binding to n class="Gene">CYP46A1 as a function of compound concentration. Traces of Δ(A437 – A416) (e.g., Figure ) were fit to single exponential plots. The error bars show the SD of the OLIS software fits. At least four separate reactions were averaged for each data point.

Next, we determined the apparent koff n class="Species">rates for cholesterol, (S)-EFV and rac-3 (Table , Figure S1). Release of cholesterol (koff) was very slow, 2 orders of magnitude slower than that of (S)- EFV and rac-3, perhaps a reason why the apparent Kd of cholesterol for CYP46A1 was also much lower than those of (S)-EFV and rac-3 (9 nM vs 5 μM). Accordingly the calculated kon rate of cholesterol was 1.7-fold lower than that for (S)-EFV but 1.9-fold higher than that for rac-3, another advantage of (S)-EFV relative to rac-3 in terms of binding to the CYP46A1 active site. Thus, not only the predominant enzyme conformation but also the steady state kon puts (S)-EFV at advantage,

Table 2

Compound Binding to CYP46A1a

compd	Kd for substrate-free CYP46A1 (M)	calculated konb (M–1 s–1)	experimental koff (s–1)
cholesterol	9 × 10–9 c	1.3 × 106	0.012 ± 0.002
(S)-EFV	5 × 10–6 d ± 0.4 × 10–6	2.2 × 106	11.2 ± 0.9
rac-3	5 × 10–6 e ± 0.3 × 10–6	0.7 × 106	3.6 ± 0.2

Experimental data are presented as the mean ± SD of the measurements from the three independent experiments.

Values calculated from Kd and koff assuming .

This value comes from an apparent Kd of 80 nM [20] after the correction for HPCD binding. See Experimental Section for details.

This is the Kd value for titrations at 24 °C [20] because the stopped-flow experiments were conducted at 23 °C.

This is the Kd value for titrations at 24 °C determined in the present work in addition to the Kd value at 18 °C (Figure ) as compared to rac-3, for the CYP46A1 active site binding and hence inhibition.

This is the Kd value for titrations at 24 °C determined in the present work in addition to the Kd value at 18 °C (Figure ) as compared to n class="Gene">rac-3, for the CYP46A1 active site binding and hence inhibition.

Discussion and Conclusion

The present work provides insight into the structure–activity relationships of (S)-EFV for n class="Gene">CYP46A1 activation. Of the four functionalities tested (the N-1, 2-keto group, O-3, and cyclopropyl ring, Figure a), only one (the O-3) does not seem to be important for CYP46A1 activation as indicated by the concentration dependence curve for rac-8 (Figure d). This curve was similar to that of (S)-EFV, consistent with compound binding to the same allosteric site as EFF (Figure b). Yet interaction with the CYP46A1 the active site was different as indicted by a type 1 spectral response in substrate-free CYP46A1 (Figure ). In addition, our studies revealed that introducing the hydroxyl group at C-7, C-8, or both C-8 and C-14 of EFV increased the maximal extent of CYP46A1 activation and eliminated the undesired CYP46A1 inhibition at higher compound concentrations (Figure b). All three hydroxylated compounds ((S)-1, rac-1, rac-2, and rac-3) appeared to bind similarly to the CYP46A1 active site (i.e., elicited a type 2 spectral response in substrate-free CYP46A1, Figure ) and seemed to interact with either the allosteric site for Glu (rac-3) or both EFV and Glu ((S)-1 and rac-2). As with rac-8, a reverse type 1 spectral response induced in cholesterol-bound CYP46A1 was a common feature of the hydroxylated compounds. However, only in the case of rac-3 did this spectral response fit to a hyperbolic curve. It is thus possible that a reverse type 1 spectral response in cholesterol-bound CYP46A1 and a hyperbolic spectral fit are the features of a compound that does not inhibit the P450 at higher concentrations while activating CYP46A1 to the extent comparable with that of (S)-EFV (Table ). Further studies are required to address this possibility.

Notably the apparent Kd values for substn class="Species">rate-free CYP46A1 (Figure ) did not appear to be indicative of CYP46A1 inhibition in the concentration dependence curves, as was demonstrated by rac-1, rac-2, and rac-3 (Figure b). To clarify the reason for a weak predictive power of an apparent Kd value, we capitalized on a recent discovery in the P450 field showing that the P450-compound interactions are complex and that conformational-selection is a dominant feature of many of these interactions for human P450s 2C6, 2D6, 3A4, 4A11, 17A1, and 21A2.[37−40] By conducting the stopped-flow experiments, we ascertained that (S)-EFV and rac-3 bind to substrate-free CYP46A1 in a mechanism involving conformational-selection, which is apparently more favorable to (S)-EFV than rac-3, a bulkier compound requiring more space in the P450 active site. In addition, a 3-fold lower kon of rac-3 than that of (S)-EFV when cholesterol release from the CYP46A1 active site (the koff) is extremely slow (Table ) could be another reason for caution when considering the apparent Kd values.

An important finding of the present work was that CYP46A1 could be activated by n class="Gene">racemic mixtures, which is in contrast to HIV reverse transcriptase, which is inhibited only by the (S)-isomer of EFV.[41] CYP46A1 activation and response to (S)-1 and rac-1 were similar but not identical (Figures b, 2, 3), and molecular docking of (S)-3 and (R)-3 suggested the same binding site but a different binding mode (Figure b,c). Thus, at least in some cases, using a racemate as compared to an expensive (S)-isomer could be more advantageous (besides lower compound cost) for activating CYP46A1 in vivo. First, we could still see the benefits of the enzyme activation while possibly minimizing the negative, non-CYP46A1 related effects[42−44] suggested by several cell culture studies, in which (S)-8-hydroxyEFV (as well as (S)-7-hydroxyEFV but to a lesser extent and not always) was shown to induce cell damage or death.[45−47] However, these cell culture findings are not supported by a study in humans showing that neither plasma nor cerebrospinal fluid levels of (S)-8-hydroxyEFV are associated with central nervous system toxicity.[29] Also, the metabolite concentrations in the cell culture studies were designed to mimic the 600 mg/day EFV dose, which is >80-fold higher (∼8.5 (mg/kg of body weight)/day for a 70 kg person) than that (0.1 (mg/kg body weight)/day) given to mice. Second, CYP2B6 metabolizes (S)-EFV 10 times more efficiently in vitro than (R)-EFV and does not generate the (R) isomers of 7-hydroxyEFV (2), 7,8-dihydroxyEFV (4), and 8,14-dihydroxyEFV (3) from (R)-EFV.[48] Therefore the plasma half-life of racemic EFV or a hydroxylated EFV metabolite could be longer and hence requires a less frequent drug administration to activate CYP46A1. Lastly, potential toxic effects of the R-isomers of EFV metabolites have not been studied and may not be so detrimental as those suggested for the (S)-isomers.[45−47] An in vivo investigation with racemic EFV and/or 8,14-dihydroxyEFV is required to evaluate potential benefits of the racemate use indicated by the current work in vitro.

Of immediate clinical significance is our finding that EFV hydroxylations at positions 7, 8, and 14, which represent initial steps in n class="Chemical">EFV clearance by the liver,[4,25−28] enhance CYP46A1 activation at lower metabolite concentrations without inhibiting the enzyme at higher concentrations (Figure b). This finding raises a possibility that not only EFV but also its metabolites activate CYP46A1 in vivo. Indeed, all these compounds are present in the plasma of EFV-treated subjects: EFV at 2170–3350 ng/mL; 8-hydroxyEFV (1) at 315–2160 ng/mL; 7-hydroxyEFV (2) at 8.8–225 ng/mL; 8,14-dihydroxyEFV (3) at 10.2 ng/mL.[4,49−54] Also, the total polar surface areas of 38.3 Å2 and 78.8 Å2 for EFV and 8,14-dihydroxyEFV, respectively, and the predicted log P values of 4.53 and 3.17, respectively (calculated using the Molinspiration property calculation service), are within the desirable range for good CNS availability (<90–120 and 2–5).[55] Consistent with these properties, EFV (13.9–19.5 ng/mL), 8-hydroxyEFV (1, 3.1–5.3 ng/mL), 7-hydroxyEFV (2, 0.63 ng/mL), and 8,14-dihydroxyEFV (3, 0.38–0.44 ng/mL) were detected in the cerebrospinal fluid, an established marker for free compound levels in the brain.[4,29,49−53] The detected compound levels in the cerebrospinal fluid were however low, probably because in the blood, EFV and EFV metabolites are ≥99.6% protein (albumin)-bound;[4,50,51] hence only small drug amounts likely cross the blood–brain barrier.[44] Then once in the brain, EFV seems to readily accumulate in the tissue due to high propensity for protein binding as indicated by the experiments showing that EFV brain levels could exceed those in the plasma 4.6 times within 1 h of an intraperitoneal administration.[56] In addition, EFV metabolites in the brain could be produced locally because EFV-metabolizing CYP2B6, CYP1A1, and CYP2C8 are not only expressed in the liver but also in the brain.[57−59] Thus, EFV metabolites could be present in the brain and activate CYP46A1 along with EFV, a possibility that would explain in part why only a small EFV dose (0.1 (mg /kg body weight)/day) is sufficient to activate CYP46A1 in mouse brain.

To conclude, structure–activity relationships of EFV for n class="Gene">CYP46A1 activation were investigated and revealed that hydroxylations at C-7, -8, and -14 increase CYP46A1 activation at lower compound concentrations and decrease enzyme inhibition at higher compound concentrations. EFV N-1, 2-keto group, and the cyclopropyl ring were found to contribute to CYP46A1 activation, whereas the drug O-3 did not seem to be important. All of the investigated EFV metabolite and analogs decreased CYP46A1 inhibition at higher concentrations, possibly because their binding to the P450 active site involves the conformational selection mechanism. We also mapped the allosteric sites of compound binding on the CYP46A1 surface. The data obtained suggest the future direction of our work. This includes testing the following compounds, first in vitro and then, if justified, in vivo as potentially better CYP46A1 activators than (S)-EFV: (S)- and rac-7,8-dihydroxyEFV as well as (S)-8,14-dihydroxyEFV. Efforts are in progress for the generation of these compounds.

Experimental Section

Materials

(S)-EFV (n class="Chemical">(S)-6-chloro-4-(cyclopropylethynyl)-4-(trifluoromethyl)-1,4-dihydro-2H-benzo[d][1,3]oxazin-2-one) and all EFV metabolites, as well as analogs, were purchased from Toronto Research Chemicals (Toronto, ON, Canada) and were >95% pure according to the manufacturer’s certificates of analysis (Figures S2–S11). 1H NMR and 19F NMR spectra of these compounds conformed to their structures. l-Glutamate was from MilliporeSigma (St. Louis, MO), and a 2 mM stock solution was prepared in water. Cholesterol was obtained from Steraloids (Newport, RI), and 24-hydroxy-[25,26,26,26,27,27,27-2H7]-cholesterol was from Medical Isotopes (Pelham, NH). Cholesterol was added from a 1 mM stock in 4.5%, w/v, aqueous 2-hydroxypropyl-β-cyclodextrin (HPCD), and deuterated 24-hydroxycholesterol was added from a 0.1 mM stock in methanol. Human truncated Δ(2–50)CYP46A1 with a four-histidine tag on the C terminus and rat cytochrome P450 oxidoreductase were expressed in Escherichia coli and purified as described.[60,61] The F405A Δ(2–50)CYP46A1 mutant was generated by using an in vitro QuikChange site-directed mutagenesis kit (Stratagene, San Diego, CA) according to the instructions. The correct generation of the desired mutation and absence of undesired mutations were confirmed by nucleotide sequencing of the entire CYP46A1 coding region as well as by the restriction analysis.

Enzyme Assays

Incubations for the concentration dependence curves were carried out in 1 mL of 50 mM n class="Chemical">KPi buffer (pH 7.2) containing 100 mM NaCl, 40 μg/mL l-α-1,2-dilauroyl-sn-glycero-3-phosphocholine, 0.5 μM purified Δ(2-50)CYP46A1, 1.0 μM cytochrome P450 oxidoreductase, 40 μM cholesterol, varying concentrations of a test compound (0–100 μM added from 0.5 to 5 mM methanol stocks to keep the 2% final methanol concentration), 2 units of catalase, and an NADPH-regenerating system (1 mM NADPH, 10 mM glucose 6-phosphate, and 2 units of glucose 6-phosphate dehydrogenase). Enzymatic reactions proceeded for 30 min at 37 °C, followed by sterol extraction with 5 mL of dichloromethane containing 1 nmol of 24-hydroxy-[25,26,26,26,27,27,27-2H7]-cholesterol, which served as an internal standard. Sterol extracts were then processed and analyzed by gas chromatography–mass spectrometry as described.[62] Incubations to determine for the maximal extent of CYP46A1 activation and the activation of the CY46A1 mutants were carried out similarly, except a fixed concentration (20 μM) of a test compound (added from a 5 mM stock in methanol) was used. Co-incubations with l-Glu were with 100 μM neurotransmitter.

Spectral Titrations

Spectral titrations were carried out at 18 °C as described,[63] i.e., in 1 mL of 50 mM n class="Chemical">KPi buffer (pH 7.2) containing 100 mM NaCl and 0.4–0.5 μM Δ(2-50)CYP46A1, either substrate free or substrate-bound, i.e., in the presence of 20 μM cholesterol. Stock solutions of a 0.5–5 mM test compound in methanol or 0.5 mM cholesterol in 4.5% (w/v) aqueous HPCD were used. Data were fit by the GraphPad Prism software to either of the following equations:or an allosteric sigmoidal curve, in which [E] is the enzyme concentration, ΔA is the spectral response at different compound concentrations [L], and ΔAmax is the maximal amplitude of the spectral response.

Computational Dockings

The structures of (S)-3 and n class="Chemical">(R)-3 were optimized with density functional theory calculations at the B3LYP/6-31+G* level of theory using Gaussian 16.[64] These architectures were then used as the input for docking experiments, using Autodock 4.2.6,[65] Autodock VINA,[66] and Ledock.[67] Input files were generated using Autodock Tools 1.5.6[68] and converted using OpenBabel 2.3.1.[69] Docking experiments were carried out on PDB code 2Q9G,[24] with a search box centered on the l-Glu binding pocket. The search box was similarly sized for all docking approaches (Autodock, 74 × 74 × 74, grid spacing 0.375 Å; Autodock VINA, 28 Å × 28 Å × 28 Å; Ledock, 28 Å × 28 Å × 28 Å). Visualization of docked compounds on the protein surface was carried out using PyMol.

Stopped-Flow Kinetics

Measurements were made in an OLIS RSM-1000 stopped-flow spectrophotometer (On-Line Instrument Systems, Bogart, GA) in the rapid-scanning mode, using a 20 mm × 4 mm cell, 16 mm × 0.2 mm Scandisk, 1.24 mm slits (8 nm bandpass), and 600 lines/500 nm gratings, open class="Species">rating at 23 °C. When the collection period was ≤4 s, data were collected at 1000 scans/s; when the collection time was >4 s, 62 scans/s were collected in the signal average mode. The wavelength range used was 330–570 nm. Measurements typically involved mixing one syringe containing 4 μM CYP46A1 in 50 mM KPi (pH 7.2), supplemented 100 mM NaCl, and another syringe containing a varying compound concentration in an equal volume of the same buffer. When CYP46A1 was premixed with equimolar cholesterol, the latter was added from a 2 mM stock in 4.5% (w/v) 2-hydroxypropyl-β-cyclodextrin. The data were saved as OLIS (.ols), and Excel files and converted to ΔAmax – ΔAmin plots. The resulting Excel files were corrected to ΔA > 0, saved as txt files, and imported into KinTek Explorer software version 8.0 (KinTek, Snowshoe, PA).[70] Analyses utilized an Apple iMac OSX 10.13.6 system.

The apparent rates of release (koff) of n class="Chemical">cholesterol, (S)-EFV, and rac-3 from CYP46A1 were estimated by mixing a prebound complex with compound under study (4 μM CYP46A1 and 4 μM compound) with an excess of the inhibitor clotrimazole (40 μM), which binds tightly and has a different spectral complex.[33,71] An increase in the Δ(A437-A416) absorbance was measured, and data from 8 to 12 plots were averaged. KinTek Explorer was used for data fitting. Since , the kon rates were then calculated from the compound Kd values for substrate-free CYP46A1 (Figure ). Cholesterol binds to both HPCD and CYP46A1; therefore the true Kd of cholesterol for CYP46A1 was calculated (9 nM) using KinTek Explorer and an apparent Kd of cholesterol for CYP46A1 (0.08 μM [20]), the Kd of cholesterol for HPCD in 50 mM KPi buffer (pH 7.2) containing 100 mM NaCl at 24 °C (540 μM, determined in the present work as described[72]), and the HPCD concentration (3.1 mM).