Evolution of anti-HIV drug candidates. Part 2: Diaryltriazine (DATA) analogues.

A synthesis program directed toward improving the stability of imidoyl thiourea based non-nucleoside reverse transcriptase inhibitors (NNRTIs) led to the discovery of diaryltriazines (DATAs), a new class of potent NNRTIs. The synthesis and anti-HIV structure-activity relationship (SAR) studies of a series of DATA derivatives are described.

Classification and identification of proteins by means of common and specific amino acid n-tuples in unaligned sequences.

Unaligned amino acid sequences can be characterized by their composition of amino acid n-tuples (i.e. doublets, triplets, quadruplets, etc.). In this study we investigated the performance of two statistics, termed commonality and specificity, that are derived from n-tuple counts using a set of G-protein coupled receptor (GPCR) sequences. The commonality of a tuple is defined as its relative occurrence in the sequences that belong to a given GPCR subtype. The specificity of a tuple is derived from its relative occurrence in the sequences of a given GPCR subtype and from its relative non-occurrence in the sequences that do not belong to this subtype. A graphical presentation, termed 'polygram', is described for the visualization of common and specific tuples. The method can be applied to the classification of unknown GPCR sequences. It can also be applied to the identification of fragments of GPCRs, such as may occur in chimeric receptors. The method is generally applicable to other protein families and other types of coding.

The alpha and omega of G-protein coupled receptors: a novel method for classification. Part 2. Bin classification.

A novel way of classification of G-protein coupled receptors is presented that is only based on receptor sequence information by counting of amino acid residues. It involves the number of amino acid residues between the Asn residue in TM1 and the residue Cys in the loop between TM4 and TM5, the number of residues between the latter Cys residue and Pro residue in TM6, and the number of residues between the latter Pro and the last amino acid residue (called omega) in the sequence. The classification of 131 sequences, covering biogenic amine, opioid and somatostatin receptors, is visualized by means of a diagram which is referred to as a bin map. Each bin in the diagram encloses all the sequences that belong to one and only one receptor type or subtype. This so-called bin classification was obtained by means of the genetic algorithm methodology, which offers new opportunities for classifying proteins.

Crystal structures of 8-Cl and 9-Cl TIBO complexed with wild-type HIV-1 RT and 8-Cl TIBO complexed with the Tyr181Cys HIV-1 RT drug-resistant mutant.

Human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) is an important target for chemotherapeutic agents used in the treatment of AIDS; the TIBO compounds are potent non-nucleoside inhibitors of HIV-1 RT (NNRTIs). Crystal structures of HIV-1 RT complexed with 8-Cl TIBO (R86183, IC50 = 4.6 nM) and 9-Cl TIBO (R82913, IC50 = 33 nM) have been determined at 3.0 A resolution. Mutant HIV-1 RT, containing Cys in place of Tyr at position 181 (Tyrl81Cys), is highly resistant to many NNRTIs and HIV-1 variants containing this mutation have been selected in both cell culture and clinical trials. We also report the crystal structure of Tyrl81Cys HIV-1 RT in complex with 8-Cl TIBO (IC50 = 130 nM) determined at 3.2 A resolution. Averaging of the electron density maps computed for different HIV-1 RT/NNRTI complexes and from diffraction datasets obtained using a synchrotron source from frozen (-165 degrees C) and cooled (-10 degrees C) crystals of the same complex was employed to improve the quality of electron density maps and to reduce model bias. The overall locations and conformations of the bound inhibitors in the complexes containing wild-type HIV-1 RT and the two TIBO inhibitors are very similar, as are the overall shapes and volumes of the non-nucleoside inhibitor-binding pocket (NNIBP). The major differences between the two wild-type HIV-1 RT/TIBO complexes occur in the vicinity of the TIBO chlorine substituents and involve the polypeptide segments around the beta5-beta6 connecting loop (residues 95 to 105) and the beta13-beta14 hairpin (residues 235 and 236). In all known structures of HIV-1 RT/NNRTI complexes, including these two, the position of the beta12-beta13 hairpin or the "primer grip" is significantly displaced relative to the position in the structure of HIV-1 RT complexed with a double-stranded DNA and in unliganded HIV-1 RT structures. Since the primer grip helps to position the template-primer, this displacement suggests that binding of NNRTIs would affect the relative positions of the primer terminus and the polymerase active site. This could explain biochemical data showing that NNRTI binding to HIV-1 RT reduces efficiency of the chemical step of DNA polymerization, but does not prevent binding of either dNTPs or DNA. When the structure of the Tyr181Cys mutant HIV-1 RT in complex with 8-Cl TIBO is compared with the corresponding structure containing wild-type HIV-1 RT, the overall conformations of Tyr181Cys and wild-type HIV-1 RT and of the 8-Cl TIBO inhibitors are very similar. Some positional changes in the polypeptide backbone of the beta6-beta10-beta9 sheet containing residue 181 are observed when the Tyr181Cys and wild-type complexes are compared, particularlty near residue Val179 of beta9. In the p51 subunit, the Cys181 side-chain is oriented in a similar direction to the Tyr181 side-chain in the wild-type complex. However, the electron density corresponding to the sulfur of the Cys181 side-chain in the p66 subunit is very weak, indicating that the thiol group is disordered, presumably because there is no significant interaction with either 8-Cl TIBO or nearby amino acid residues. In the mutant complex, there are slight rearrangements of the side-chains of other amino acid residues in the NNIBP and of the flexible dimethylallyl group of 8-Cl TIBO; these conformational changes could potentially compensate for the interactions that were lost when the relatively large tyrosine at position 181 was replaced by a less bulky cysteine residue. In the corresponding wild-type complex, Tyr181 iin the p66 subunit has significant interactions with the bound inhibitor and the position of the Tyr181 side-chain is well defined in both subunits. Apparently the Tyr181 --> Cys mutation eliminates favorable contacts of the aromatic ring of the tyrosine and the bou

The alpha and omega of G-protein coupled receptors: a novel method for classification. Part 1.

In this paper, a novel method is presented for classification of 113 neurotransmitter and opioid receptors that is based on the counting of amino acid residues. With the use of sequence alignment, ten conserved key residues were identified: alpha (the first Met residue of the sequence), Asn in the tentative first transmembrane domain (TM1), Asp (TM2), Cys (top of TM3), Arg (bottom of TM3), Trp (TM4), Cys (in the loop between TM4 and TM5), Pro (TM5), Pro (TM6), Pro (TM7) and omega (the last residue of the sequence). The number of residues between these key residues is unique for each receptor or receptor subtype and is used for classification. The number of residues between two key residues defines a segment. The sum of the segments before or after a key residue is defined as a partition. In total, 73 possible classification schemes were found using two or three segments, partitions or a combination of segments and partitions. The surprising and striking results is that each of the sequences examined can be characterized by a code consisting of two or three figures. Each figure represents a number of amino acid residues.

Molecular modeling studies of HIV-1 reverse transcriptase nonnucleoside inhibitors: total energy of complexation as a predictor of drug placement and activity.

Computer modeling studies have been carried out on three nonnucleoside inhibitors complexed with human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT), using crystal coordinate data from a subset of the protein surrounding the binding pocket region. Results from the minimizations of solvated complexes of 2-cyclopropyl-4-methyl-5,11-dihydro-5H-dipyrido[3,2-b :2',3'-e][1,4] diazepin-6-one (nevirapine), alpha-anilino-2, 6-dibromophenylacetamide (alpha-APA), and 8-chloro-tetrahydro-imidazo(4,5,1-jk)(1,4)-benzodiazepin-2(1H)-thi one (TIBO) show that all three inhibitors maintain a very similar conformational shape, roughly overlay each other in the binding pocket, and appear to function as pi-electron donors to aromatic side-chain residues surrounding the pocket. However, side-chain residues adapt to each bound inhibitor in a highly specific manner, closing down around the surface of the drug to make tight van der Waals contacts. Consequently, the results from the calculated minimizations reveal that only when the inhibitors are modeled in a site constructed from coordinate data obtained from their particular RT complex can the calculated binding energies be relied upon to predict the correct orientation of the drug in the pocket. In the correct site, these binding energies correlate with EC50 values determined for all three inhibitors in our laboratory. Analysis of the components of the binding energy reveals that, for all three inhibitors, solvation of the drug is endothermic, but solvation of the protein is exothermic, and the sum favors complex formation. In general, the protein is energetically more stable and the drug less stable in their complexes as compared to the reactant conformations. For all three inhibitors, interaction with the protein in the complex is highly favorable. Interactions of the inhibitors with individual residues correlate with crystallographic and site-specific mutational data. pi-Stacking interactions are important in binding and correlate with drug HOMO RHF/6-31G* energies. Modeling results are discussed with respect to the mechanism of complex formation and the design of nonnucleoside inhibitors that will be more effective against mutants of HIV-1 RT that are resistant to the currently available drugs.

A molecular model for the interaction between vorozole and other non-steroidal inhibitors and human cytochrome P450 19 (P450 aromatase).

In a previous study (Vanden Bossche et al., Breast Cancer Res. Treat. 30 (1994) 43) the interaction between (+)-S-vorozole and the I-helix of cytochrome P450 19 (P450 aromatase) has been reported. In the present study we extended the "I-helix model" by incorporating the C-terminus of P450 aromatase. The crystal structures of P450 101 (P450 cam), 102 (P450 BM-3) and 108 (P450 terp) reveal that the C-terminus is structurally conserved and forms part of their respective substrate binding pocket. Furthermore, the present study is extended to the interaction between P450 aromatase and its natural substrate androstenedione and the non-steroidal inhibitors (-)-R-vorozole, (-)-S-fadrozole, R-liarozole and (-)-R-aminoglutethimide. It is found that (+)-S-vorozole, (-)-S-fadrozole and R-liarozole bind in a comparable way to P450 aromatase and interact with both the I-helix (Glu302 and Asp309) and C-terminus (Ser478 and His480). The weak activity of (-)-R-aminoglutethimide might be attributed to a lack of interaction with the C-terminus.

Structure of HIV-1 RT/TIBO R 86183 complex reveals similarity in the binding of diverse nonnucleoside inhibitors.

We report the structure of HIV-1 reverse transcriptase (RT) complexed with the nonnucleoside inhibitor TIBO R 86183 at 3.0 A resolution. Comparing this structure with those of complexes of HIV-1 RT/alpha-APA R 95845 and HIV-1 RT/nevirapine provides a basis for understanding the nature of nonnucleoside inhibitor binding, the structure of the binding site and the interactions between the bound inhibitors and surrounding amino acid residues as well as for understanding mechanisms of inhibition by and resistance to nonnucleoside inhibitors. All three inhibitors considered assume a similar butterfly-like shape and bind to HIV-1 RT in a very similar way. Important differences occur in the conformation of amino acid residues that form the binding pocket.

P450 inhibitors of use in medical treatment: focus on mechanisms of action.

A number of cytochrome P450s are targets for compounds that are clinically used or under clinical evaluation for treatment of patients with mycotic infections, such as dermatophytosis, superficial and systemic candidiasis, cryptococcosis and aspergillosis, with skin diseases, such as psoriasis or ichthyosis, and other retinoid-sensitive malignancies, e.g., neuro-ectodermal glioma. Some of the P450 inhibitors are candidates for the treatment of hirsutism or prostate cancer, others are potent inhibitors of the P450 isomerase involved in the synthesis of thromboxane A2, a potent platelet aggregation inducer and vasoconstrictor.

Aromatase inhibitors--mechanisms for non-steroidal inhibitors.

The conversion of androgens to estrogens occurs in a variety of cells and tissues, such as ovarian granulosa and testicular cells, placenta, adipose tissue, and various sites of the brain. The extragonadal synthesis of estrogens has great pathophysiological importance. Estrogens produced by, for example, adipose tissue have a role in the pathogenesis of certain forms of breast cancer and endometrial adenocarcinoma. The biosynthesis of estrogens is catalyzed by the aromatase, an enzyme localized in the endoplasmic reticulum that consists of two components: a cytochrome P450 (P450 Arom, P450 19 product of the CYP19 gene) and the NADPH cytochrome P450 reductase. The alignment of the amino acid sequences of human P450 19 with other mammalian P450s shows little sequence similarity, which indicates not only that P450 19 is a unique form of the P450 superfamily but also that the aromatase may be a good target for the development of selective P450 inhibitors. Aminoglutethimide (AG) is the pioneer drug of the reversible competitive nonsteroidal aromatase inhibitors. Since AG is a nonspecific aromatase inhibitor and presents some problems with tolerability, a number of structural analogues have been synthesized. For example, rogletimide is slightly less potent than AG but has the advantage of not inhibiting the cholesterol side-chain cleavage and is devoid of sedative action. Elongation of the ethyl substituent of AG and rogletimide leads to an increase in aromatase inhibition. Further studies led to the discovery of a new generation of much more potent aromatase inhibitors. An example is fadrozole. However, although fadrozole is a poor inhibitor of the cholesterol side-chain cleavage, it suppresses aldosterone release by ACTH-stimulated human adrenocortical cells. More selective aromatase inhibitors are the triazole derivatives. Examples are CGS 20267, CGS 47645, R 76 713, and ICI D1033. R 76 713's aromatase inhibitory effect is largely due to its (+)-S-enantiomer, vorozole. Computer modeling studies of the interaction of vorozole with part of the "I-helix" of P450 19 suggest that the chlorine-substituted phenyl ring of vorozole interacts with the gamma-carbonyl group of Glu-302. Thr-310, which corresponds to the highly conserved Thr-252 in P450 101, interacts with vorozole's triazole ring, and the 1-methyl-benzotriazole moiety binds near Asp-309.

Genomic cloning, heterologous expression and pharmacological characterization of a human histamine H1 receptor.

A human histamine H1 receptor gene lacking introns was isolated by screening a human genomic library with a bovine histamine H1 receptor probe. The deduced protein of 487 amino acids showed characteristic properties of G-protein-coupled receptors. The coding region was subcloned into the expression vector pSVL (Pharmacia), and the resulting construct transfected into COS-7 cells. Binding studies with [3H]pyrilamine on membranes from transfected cells revealed saturable specific binding with a KD of 1.2 nM and a Bmax of 3400 fmol/mg protein. Binding affinities of histamine and known histamine antagonists were similar to those for histamine H1 receptors in guinea-pig cerebellum. In transfected COS-7 cells, histamine induced inositol phosphate formation, that was inhibitable by pyrilamine.

Novel computational model for the interaction of dopamine with the D2 receptor.

With the use of molecular modelling and computational chemistry, two models were obtained showing the binding of dopamine to part of the D2 receptor. These models consist of a molecular complex containing dopamine, the transmembrane domains (TM) 3, 4 and 5, a small part of the extracellular loop just before TM3 and a larger part of the loop connecting TM4 and TM5. The models are further characterized by the presence of two disulphide bonds: one between TM3 and TM4 and a second linking the Cys residue before TM3 with the Cys residue in the loop between TM4 and TM5. In the first model, the beta-adrenergic receptor is used as a reference with the catechol moieties of dopamine interacting with two Ser residues in TM5, and the Asp residue in TM3 interacting with the protonated nitrogen of dopamine. In the second model, dopamine is positioned between the TM without interaction assumptions. In both models the TM positions are close to the arrangement of the alpha-helix ordering as found in bacteriorhodopsin. Energy minimization was accomplished using these two starting structures. The two minimized complexes show different molecular interactions that make good chemical sense and are in agreement with reported X-ray data of different proteins. The disulphide bond between TM3 and 4 is apparently unique to the D2 and D3 receptors. The first minimized complex revealed that not Ser194, as in the beta-receptor, but Ser 193 interacted with the meta-hydroxyl of dopamine, although a proper hydrogen bond with Ser194 was initially present in the starting structure. The para-hydroxyl group interacted with Ser197.(ABSTRACT TRUNCATED AT 250 WORDS)

The sequence homologies of cytochromes P-450 and active-site geometries.

The amino acid sequence alignment of 16 cytochrome P-450 proteins representative of the major families is reported. The sequence matching process has been carried out on the basis of maximum homology by residue type, retention of secondary structure and minimization of deletions/insertions except where additional loop regions exist. From the starting point of known reported sequence homology matching from the literature, a realignment on the basis of conserved residues involved in both structure and function gives rise to a self-consistent set of sequences which correlates with known mechanistic and structural data. Once fitted, these archetypal sequences form a straightforward template for the alignment of all P-450 subfamilies. Computer modelling of the active-site regions constructed from homology with the bacterial form of the enzyme (P-450CAM) evinces the correct substrate specificity. Furthermore, the construction of the macromolecular assembly of components of the cytochrome P-450 system on the microsomal endoplasmic reticular membrane is presented from the evidence of site-directed mutagenesis, analysis by molecular probes, X-ray crystallography and molecular modelling.

R 76713 and enantiomers: selective, nonsteroidal inhibitors of the cytochrome P450-dependent oestrogen synthesis.

The triazole derivative, R 76713 and its enantiomers R 83839(-) and R 83842(+) are effective inhibitors of the aromatization of androstenedione. For human placental microsomes, the (+) enantiomer (R 83824) is about 1.9- and 32-times more active than the racemate (IC50 2.6 nM) and the (-) enantiomer, respectively. R 83842 is about 30- and 1029-times more active than 4-hydroxyandrostene-3,17-dione and aminoglutethimide. This potency might originate from its high affinity for the microsomal cytochrome P450 (P450). Indeed, R 83842, compared to R 76713 and R 83839, forms a more stable P450-drug complex. Difference spectral measurements indicate that the triazole nitrogen N-4 coordinates to the haem iron. The reversed type 1 spectral changes suggest that R 76713 is able to displace the substrate from its binding place and the stable complex formed in particular with the (+) enantiomer suggests that its N-1-substituent occupies a lipophilic region of the apoprotein moiety. Kinetic analysis implies that there is a competitive part in the inhibition of the human placental aromatase by R 76713. The Ki values for R 76713, R 83842 and R 83839 are 1.3 nM, 0.7 nM and 18 nM, respectively. These results are indicative of stereospecificity for binding. Up to 10 microM, R 76713 and its enantiomers have no statistically significant effect on the regio- and stereoselective oxidations of testosterone in male rat liver microsomes. All three compounds have no effect on the P450-dependent cholesterol synthesis, cholesterol side-chain cleavage and 7 alpha-hydroxylation and 21-hydroxylase. At 10 microM, R 76713 has a slight effect on the bovine adrenal 11 beta-hydroxylase. This effect originates mainly from R 83839, the less potent aromatase inhibitor. On the other hand, the inhibition of the 17,20-lyase of rat testis observed at concentrations greater than or equal to 0.5 microM, originates rather from R 83842. However, 50% inhibition is only achieved at 1.8 microM R 83842, i.e. at a concentration about 1300-times higher than that needed to reach 50% inhibition of the human placental aromatase.

CGEMA and VGAP: a Colour Graphics Editor for Multiple Alignment using a variable GAP penalty. Application to the muscarinic acetylcholine receptor.

Today, more than 40 protein amino acid (AA) sequences of membrane receptors coupled to guanine nucleotide binding proteins (G-proteins) are available. For those working in the field of medicinal chemistry, these sequences present a new type of information that should be taken into consideration. To make maximal use of sequence data it is essential to be able to compare different protein sequences in a similar way to that used for small molecules. A prerequisite, however, is the availability of a processing environment that enables one to handle sequences in an easy way, both by hand and by computer. In order to meet these ends, the package CGEMA (Colour Graphics Editor for Multiple Alignment) was developed in our laboratory. The programme uses a user-definable colour coding for the different AAs. Sequences can be aligned by hand or by computer, using VGAP, and both approaches can be combined. VGAP is a novel in-house written alignment programme with a variable gap penalty that also handles consecutive alignments using one sequence as a probe. In addition, secondary structure prediction tools are available. From the 20 protein sequences, available for the muscarinic acetylcholine receptor, 13 different sequences were selected, covering the subtypes m1 to m5. By comparing the sequences, two major groups are revealed that correspond to those found by considering the transducing system coupled to the various receptor subtypes. Different parts of the protein sequences are identified as characterizing the subtype and binding the ligands, respectively.

Two groups of rhinoviruses revealed by a panel of antiviral compounds present sequence divergence and differential pathogenicity.

A variety of chemically different compounds inhibit the replication of several serotypes of rhinoviruses (common-cold viruses). We noticed that one of these antiviral compounds, WIN 51711, had an antiviral spectrum clearly distinctive from a consensus spectrum or other capsid-binding compounds, although all of them were shown to share the same binding site. A systematic evaluation of all known rhinovirus capsid-binding compounds against all serotyped rhinoviruses was therefore initiated. Multivariate analysis of the results revealed the existence of two groups of rhinoviruses, which we will call antiviral groups A and B. The differential sensitivity of members of these groups to antiviral compounds suggests the existence of a dimorphic binding site. The antiviral groups turned out to be a reflection of a divergence of rhinovirus serotypes on a much broader level. Similarities in antiviral spectra were highly correlated with sequence similarities, not only of amino acids lining the antiviral compound-binding-site, but also of amino acids of the whole VP1 protein. Furthermore, analysis of epidemiological data indicated that group B rhinoviruses produced more than twice as many clinical infections per serotype than group A rhinoviruses did. Rhinoviruses belonging to the minor receptor group were without exception all computed to lie in the same region of antiviral group B.

Biochemical approaches to selective antifungal activity. Focus on azole antifungals.

Azole antifungals (e.g. the imidazoles: miconazole, clotrimazole, bifonazole, imazalil, ketoconazole, and the triazoles: diniconazole, triadimenol, propiconazole, fluconazole and itraconazole) inhibit in fungal cells the 14 alpha-demethylation of lanosterol or 24-methylenedihydrolanosterol. The consequent inhibition of ergosterol synthesis originates from binding of the unsubstituted nitrogen (N-3 or N-4) of their imidazole or triazole moiety to the heme iron and from binding of their N-1 substituent to the apoprotein of a cytochrome P-450 (P-450(14)DM) of the endoplasmic reticulum. Great differences in both potency and selectivity are found between the different azole antifungals. For example, after 16h of growth of Candida albicans in medium supplemented with [14C]-acetate and increasing concentrations of itraconazole, 100% inhibition of ergosterol synthesis is achieved at 3 x 10(-8) M. Complete inhibition of this synthesis by fluconazole is obtained at 10(-5) M only. The agrochemical imidazole derivative, imazalil, shows high selectivity, it has almost 80 and 98 times more affinity for the Candida P-450(s) than for those of the piglet testes microsomes and bovine adrenal mitochondria, respectively. However, the topically active imidazole antifungal, bifonazole, has the highest affinity for P-450(s) of the testicular microsomes. The triazole antifungal itraconazole inhibits at 10(-5) M the P-450-dependent aromatase by 17.9, whereas 50% inhibition of this enzyme is obtained at about 7.5 x 10(-6)M of the bistriazole derivative fluconazole. The overall results show that both the affinity for the fungal P-450(14)DM and the selectivity are determined by the nitrogen heterocycle and the hydrophobic N-1 substituent of the azole antifungals. The latter has certainly a greater impact. The presence of a triazole and a long hypdrophobic nonligating portion form the basis for itraconazole's potency and selectivity.