Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, February 2004, p. 568-574, Vol. 48, No. 2
0066-4804/04/$08.00+0 DOI: 10.1128/AAC.48.2.568-574.2004
Copyright © 2004, American Society for Microbiology. All Rights Reserved.
-Sterol Demethylases from Aspergillus fumigatus and Candida albicans Provide Insights into Posaconazole Binding
Departments of Structural Chemistry,1 Antimicrobial Therapy and Molecular Genetics, Schering-Plough Research Institute, Kenilworth, New Jersey 070332
Received 27 August 2003/ Returned for modification 3 October 2003/ Accepted 24 October 2003
|
|
|---|
-demethylase enzyme (CYP51) is the target of azole antifungals. Azoles block ergosterol synthesis, and thereby fungal growth, by binding in the active-site cavity of the enzyme and ligating the iron atom of the heme cofactor through a nitrogen atom of the azole. Mutations in and around the CYP51 active site have resulted in azole resistance. In this work, homology models of the CYP51 enzymes from Aspergillus fumigatus and Candida albicans were constructed based on the X-ray crystal structure of CYP51 from Mycobacterium tuberculosis. Using these models, binding modes for voriconazole (VOR), fluconazole (FLZ), itraconazole (ITZ), and posaconazole (POS) were predicted from docking calculations. Previous work had demonstrated that mutations in the vicinity of the heme cofactor had a greater impact on the binding of FLZ and VOR than on the binding of POS and ITZ. Our modeling data suggest that the long side chains of POS and ITZ occupy a specific channel within CYP51 and that this additional interaction, which is not available to VOR and FLZ, serves to stabilize the binding of these azoles to the mutated CYP51 proteins. The model also predicts that mutations that were previously shown to specifically impact POS susceptibility in A. fumigatus and C. albicans act by interfering with the binding of the long side chain. |
|
|---|
-demethylase (CYP51, also called Erg11p). Ergosterol depletion, coupled with the accumulation of methylated sterol precursors, has been proposed to affect both membrane integrity and the function of some membrane-bound proteins, including chitin synthase (27). The net result is an inhibition of fungal growth. Resistance to azoles is a concern, particularly during the long-term treatment of AIDS patients with oropharyngeal candidiasis (25). The two most-prevalent causes of resistance are mutations in the target site, resulting in reduced azole binding to CYP51, and decreased intracellular drug accumulation, resulting from increased expression of efflux pump genes (recently reviewed in reference 24). To combat azole resistance, and to extend the spectrum of treatable pathogens, more-potent azoles have been developed. One such agent is posaconazole (POS) (SCH56592), a broad-spectrum triazole in phase III trials. Unlike fluconazole (FLZ) and voriconazole (VOR), POS is not effluxed by the pumps encoded by FLU1 and MDR1. However, POS is a substrate for the ATP-dependent pumps encoded by CDR1 and CDR2. In addition, POS appears to be relatively insensitive to amino acid substitutions in CYP51 (D. Sanglard, F. Ischer, and J. Bille, Abstr. 42nd Intersci. Conf. Antimicrob. Agents Chemother., abstr. M-221, 2002). Finally, POS is active against Aspergillus spp. (17).
The CYP51 proteins belong to the cytochrome P450 superfamily. Unlike the soluble bacterial P450s, all the fungal CYP51 proteins characterized to date are integral membrane proteins, making structural and biophysical characterization more challenging. The publication of the X-ray crystal structure of CYP51 from Mycobacterium tuberculosis (MT-CYP51), which has >25% sequence identity to most fungal CYP51s, provided a new opportunity to study the fungal enzymes (18, 19). In this work, we used the MT-CYP51 structure to construct homology models for the CYP51 proteins from Aspergillus fumigatus (AF-CYP51A) and Candida albicans (CA-CYP51). To gain insight into how the various azoles bind these proteins, we combined docking calculations with information gained from studying substitutions that resulted in azole resistance. The resultant binding models help explain why particular CYP51 mutations have more-profound affects on the binding of some azoles than others. They may also provide insights into the design of azoles with improved affinity for their target site.
|
|
|---|
Docking of ligands. The torsion angles of several residue side chains around the azole binding cavity in the CYP51 models were manually adjusted and in some cases further minimized using the Tripos force field, Amber charges, and a distance-dependent dielectric model. The three-dimensional structures of all azole ligands were constructed using Sybyl (version 6.9) accessed via CONCORD (Tripos Inc., St. Louis, Mo.). Gasteiger-Huckel charges (Sybyl; Tripos Inc.) were assigned to the azole ligand atoms and heme.
In the MT-CYP51/FLZ X-ray structure, FLZ binds with its triazole ring perpendicular to the heme porphyrin plane and with N-4 of the triazole coordinated to the heme iron. The other azoles were constrained so that their triazole rings overlapped the FLZ triazole. Since VOR is very similar to FLZ, its initial binding conformation was derived by overlap with FLZ in MT-CYP51/FLZ and then overlap of the AF-CYP51A and CA-CYP51 models with MT-CYP51. The program GOLD (version 2.0; CCDC, Cambridge, United Kingdom), an automated docking program (8), was used to dock POS and ITZ into the binding site in AF-CYP51A and CA-CYP51. The active-site radius was set at 15 Å. For the docking calculations, all residues of the protein and heme were fixed. while the ligands were flexible. A covalent-bond distance of 2 Å from N-4 of the triazole ring to the heme iron was set as a constraint. Twenty-one genetic algorithm runs were performed.
|
|
|---|
root mean square deviation between the template and the model was 0.5 Å over 442 pairs. The largest insertion, a 15-residue loop (residues 427 to 441), was located far from the active site and was not modeled in detail. Instead, an approximate loop model from MODELLER was used.
![]() View larger version (69K): [in a new window] |
FIG. 1. Alignment of the amino acid sequences of MT-CYP51 with CA-CYP51 and AF-CYP51A. The N-terminal membrane-spanning domains of CA-CYP51 (residues 1 to 40) and AF-CYP51A (residues 1 to 29) were deleted prior to alignment. The secondary structure elements are labeled, and residues involved in binding to POS are highlighted.
|
![]() View larger version (45K): [in a new window] |
FIG. 2. Overall structure of the homology model of AF-CYP51A. Helices and strands are shown in brown and blue. The heme cofactor is highlighted in green. Two potential ligand access channels are indicated.
|
Identification of potential azole binding sites. The heme cofactor was extracted from the MT-CYP51 structure and merged into the AF-CYP51A and CA-CYP51 models. To accommodate the heme, Phe134 in AF-CYP51A and Phe145 in CA-CYP51 were both changed to Ala so that only the Cß side chain atom was retained. Both residues are located in the middle of helix C with their side chains pointing directly at the heme. For the corresponding residue in the MT-CYP51 structure (Leu100) only the Cß atom coordinates of the side chain were reported, as all the other side chain atoms were disordered. The heme structure and surrounding side chain conformations were energy minimized to avoid steric conflicts. Two CA-CYP51 models, based on MT-CYP51, have been recently published (4, 11). In both models, a two-residue insertion that we placed before helix C (Fig. 1) was placed in the middle of the helix. Helix C is in close proximity to the heme and the ligand binding site. In the MT-CYP51/FLZ X-ray structure, helix C exhibited high levels of B-factors along its entire length, indicating thermal motion or multiple conformations. In addition, some side chain atoms of helix C were disordered. Experimental uncertainties in this region resulted in imprecise models; nevertheless, the helix C conformation of MT-CYP51 was retained in our models.
Through an analysis of all known bacterial and mammalian P450 X-ray structures (11, 18, 28, 29), two putative ligand entry channels were identified that could accommodate the long side chains of POS and ITZ (Fig. 2 and 3). Channel 1 is parallel to the plane of the heme and is formed by the bent helix I, the open BC loop, and helices B' and G. This channel is apparent in both MT-CYP51 and rabbit P450 2C5 X-ray structures (18, 28). Channel 2, based on bacterial P450BM3 structures (10), is perpendicular to the plane of the heme and is located between the ß-sheet and helical domains. The surface entrance of channel 2 in our model is an area surrounded by the FG loop, helix A', and the ß turn connecting strands ß4-1 and ß4-2. We explored both channels in AF-CYP51A and CA-CYP51 for binding of the long POS side chain. Only the docking results which placed POS in channel 2 are consistent with the reported resistance mutations (see below). Although channel 2 is apparent in the MT-CYP51 X-ray structure, it is not open to the surface. In the model of AF-CYP51A, the conformations of the side chains of residues Tyr107, Phe119, and Leu494 were modified to open an entrance to the surface. In addition, the side chain of Pro216, located in the FG loop, was deleted (specifically, side chain atoms Cß, C
, and C
were deleted) to leave enough space for the distal portion of the long side chain of POS. This modification served as an approximation for a possible backbone rearrangement of the FG loop. The latter is plausible since the FG loop has large B factors in the X-ray structure of MT-CYP51 indicating thermal motion or multiple conformations (18). Similar changes for the corresponding side chains in the CA-CYP51 model, including Pro230, were made.
![]() View larger version (22K): [in a new window] |
FIG. 3. Structures of POS (A), ITZ (B), VOR (C), and FLZ (D).
|
![]() ![]() View larger version (58K): [in a new window] |
FIG. 4. Azole binding sites in the CYP51 proteins. POS (in stereoview) (A) and VOR (B) binding sites in AF-CYP51A; POS (C) and VOR (D) binding sites in CA-CYP51. For clarity, only selected binding site residues are shown. The heme cofactor is highlighted in green. The azoles are shown as ball-and-stick models with their carbon atoms in yellow. Similarly, for the binding site residues the carbon, oxygen, nitrogen, and fluorine atoms are shown in grey, red, blue, and light blue, respectively.
|
VOR was docked into AF-CYP51A by overlapping with the X-ray structure of FLZ bound to MT-CYP51. The minimized binding conformation, highlighting active-site residues within 5 Å of VOR, is shown in Fig. 4B. Since VOR is a relatively compact molecule, binding does not involve residues from helix A' or the FG loop. Instead, the primary binding determinants are hydrophobic interactions between VOR and residues Phe115, Met292, Ala293, Ile364, Ile367, and Leu494. In addition, the Tyr107 side chain displays an aromatic interaction with the pyrimidine ring of VOR. It is noteworthy that VOR and POS appear to have similar interactions with the protein near the heme.
Similar docking strategies were applied to perform docking calculations on POS, ITZ, and VOR in the CA-CYP51 binding site (Fig. 4C and D). Residues that interact with POS are Ala61, Ala62, Tyr64, Gly65, Leu87, Leu88, Tyr118, Leu121, Thr122, Phe126, Phe145, Phe228, Thr229, Pro230, Ile231, Phe233, Met306, Gly307, His310, Thr311, Leu376, Ser378, Ile379, Phe380, Ser506, Ser507, and Met508. Residues that interact with VOR are Tyr118, Thr122, Phe126, Phe145, Phe228, Met306, Gly307, His310, Thr311, Leu376, Ile379, and Met508.
Locations of substitutions in CYP51 resulting in azole resistance. (i) A. fumigatus. Resistance to POS and ITZ in A. fumigatus results from substitutions at Gly54 in AF-CYP51A (3, 13, 15). Gly54 is located on helix A' next to the triazole-3-one five-membered ring in the long side chain of POS (Fig. 5). This part of the POS molecule fits snugly inside the pocket formed by helix A', the FG loop, and the ß turn between ß-strands ß4-1 and ß4-2. Replacing Gly54 with amino acids that have side chains would create van der Waals (VDW) conflicts between the side chains and POS. The spectrum of substitutions at Gly54 and the resultant changes in MIC appear to support this prediction (Table 1). Replacing Gly54 with either Glu or Arg resulted in a 30-fold increase in the POS MIC. Introduction of Trp, with its large rigid side chain, resulted in a >250-fold MIC increase (13). ITZ occupies the same binding site as POS; it was therefore not surprising that the substitutions at Gly54 also resulted in large increases in the ITZ MIC.
![]() View larger version (30K): [in a new window] |
FIG. 5. Mapping of substitutions in AF-CYP51A that cause azole resistance. The heme cofactor is highlighted in green. POS (A) and VOR (B) are shown as ball-and-stick models with their carbon atoms in yellow. The atoms of the residues at the sites of substitutions are depicted using space filling models.
|
|
View this table: [in a new window] |
TABLE 1. Substitutions in AF-CYP51A resulting in azole resistance in A. fumigatus isolates
|
(ii) C. albicans. An analysis of two series of C. albicans isolates, from two different patients, revealed that stepwise reductions in azole susceptibility were accompanied by the appearance of specific substitutions in CYP51. The first isolate from patient I had two substitutions in CA-CYP51, Tyr257His and Gly464Ser (Fig. 6). Tyr257 is far from the azole binding site and its effect on azole binding is not immediately obvious (Table 2). Gly464 corresponds to Gly448 in AF-CYP51A and is therefore predicted to disturb the heme environment. As seen for the analogous substitution in A. fumigatus, Gly464Ser resulted in a larger decrease in susceptibility to VOR and FLZ than to POS. The second isolate from patient I acquired a third substitution, Gly307Ser. Since Gly307 is located on helix I and contacts both the triazole and phenyl rings, the substitution would be predicted to negatively impact binding of all azoles through steric conflicts between the polar Ser side chain with the triazole and phenyl rings. In practice, the substitution impacted binding of VOR more than POS. The final isolate from patient I acquired a fourth substitution, Ala61Val. Like Gly54 in AF-CYP51A, Ala61 is located on helix A'. Consequently, the additional increase in the MIC of POS most likely resulted from VDW conflicts between the large Val side chain and the end of the POS side chain.
![]() View larger version (42K): [in a new window] |
FIG. 6. Mapping of substitutions in CA-CYP51 that cause azole resistance. The heme cofactor is highlighted in green. The azoles are shown as ball-and-stick models with their carbon atoms in yellow. The atoms of the residues at the sites of substitutions are depicted using space filling models.
|
|
View this table: [in a new window] |
TABLE 2. Substitutions in CA-CYP51 resulting in azole resistance in clinical C. albicans isolates
|
|
|
|---|
-sterol demethylase from M. tuberculosis, coupled with the identification of substitutions in fungal CYP51 enzymes that confer resistance to specific structural classes of azoles, has provided an opportunity to construct homology models that are potentially superior to those based on other bacterial P450 structures (1, 6, 7, 9, 18, 26). This study provides a model for the binding of azoles with extended side chains such as POZ and ITZ. The key finding from our work is that the long side chains of POS and ITZ appear to be contained within channel 2 and consequently make extensive hydrophobic contacts along their entire lengths. This channel is evident in the X-ray structure of MT-CYP51 as a cavity of 2,600 Å3 (18). The long-chain azoles span channel 2, extending from the heme to the surface of the protein, ending in a region surrounded by helix A', the ß-turn connecting ß4-1 and ß4-2, and the FG loop (compare Fig. 2 and 4). The model elegantly explains how a substitution, at Gly54 in the A' helix of AF-CYP51A, confers resistance to POS and ITZ by perturbing the binding of the long side chain in channel 2. Similarly, substitutions in either the A' helix (at residue Ala61) or the FG loop (at residue Pro230) of CA-CYP51, in combination with substitutions near the heme site, specifically impact resistance to POS by the same mechanism. However, unlike Gly54 in AF-CYP51A, substitutions in CA-CYP51 at Ala61 or Pro230 alone did not confer reduced susceptibility to POS when the corresponding CYP51 allele was expressed in Saccharomyces cerevisiae (10a; A. Chau et al., unpublished data).
From the model, it seems probable that the added interactions furnished by the long side chain of POS and ITZ would result in inherently tighter binding affinities compared with azoles lacking such substituents. Although precise biophysical measurements are lacking to prove this hypothesis, the resistance data appear to support this contention. Specifically, mutations near the heme site, which result in significant levels of resistance to FLZ and VOR in C. albicans and A. fumigatus, have far less impact on the susceptibility of the organisms to POS and ITZ (this study; Sanglard et al., 42nd ICAAC; Manavathu et al, 41st ICAAC). We propose that the lack of impact of these mutations on POS and ITZ results from inherently tighter affinity and/or reflects compensatory adjustments for binding made along channel 2 by these drugs. This hypothesis is in direct contrast to the work of Ji et al. (7). These workers constructed a CA-CYP51 model based on the bacterial P450BM3 X-ray structure and evaluated the binding of ITZ. Although they also placed the long side chain of ITZ in a channel, the side chain reached the surface of the protein at a different location (surrounded by helix A, ß1-1, ß1-2, and the FG loop); their alignment for the A' helix also differed from the one described here. These are both features that in the present model help explain how particular mutations specifically impact the binding of POS and ITZ. Based on their model, Ji et al. (7) asserted that "the side chains of the inhibitors were not the determinants of activity" but merely adjusted physicochemical or pharmacokinetic attributes. Our model for POS binding, along with the supporting resistance data, argues for a much more significant role of the side chains in determining activity. Furthermore, our model may explain differential effects of substituents on the phenyl side chain close to the heme. The MIC data suggested that the Glu and Arg substitutions of A. fumigatus residue Gly54 impact the binding of ITZ more than POS. This may be because the larger substituents on the phenyl ring of ITZ, in a sterically restricted pocket, exacerbate the effects of the substitutions at Gly54.
Two other CA-CYP51 models based on MT-CYP51 were recently published (4, 11). Neither model addressed the binding of long-chain azoles. In addition, both publications used sequence alignments that differ from ours, which resulted in differences in both the pockets close to the heme, as well as in the remote binding determinants for POS and ITZ.
The present study provides the first description of a potential binding mode for azoles in CYP51A from A. fumigatus and may explain how specific mutations in AF-CYP51A confer resistance to one structural group of azoles (e.g., the compact azoles FLZ and VOR) but not to another (e.g., extended azoles ITZ and POS) and vice versa. The model also provides a working framework for the rational design of other azoles with extended substituents that could similarly exploit channel 2 for added binding interactions.
|
|
|---|
-sterol demethylase gene cyp51A contributes to itraconazole resistance in Aspergillus fumigatus. Antimicrob. Agents Chemother. 47:1120-1124.
-demthylase of Candida albicans and its interaction with azole antifungals. J. Med. Chem. 43:2493-2505.[CrossRef][Medline]
-demethylase. Antimicrob. Agents Chemother. 47:577-581.
-sterol demethylase (CYP51) from Mycobacterium tuberculosis in complex with azole inhibitors. Proc. Natl. Acad. Sci. USA 98:3068-3073.
-sterol demethylase from comparative analysis of amino acid sequences and X-ray structure of Mycobacterium tuberculosis CYP51. J. Inorg. Biochem. 87:227-235.[CrossRef][Medline]
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»