The structure of HIV-1 reverse transcriptase complexed with 9-chloro-TIBO: lessons for inhibitor design

Synthesis of 2-(aminocarbonylmethylthio)-1H-imidazoles as novel Capravirine analogues

Different analogues of Capravirine (AG-1549) or S-1153 were prepared by synthesis of 2-(5-benzyl-4-isopropyl-1-methyl-2,3-dihydro-1H-imidazol-2-ylthio)acetamide (3a-c), ethyl [5-benzyl-1-(ethoxymethyl)-4-ethyl-1H-imidazol-2-ylthio]acetate (10), 2-[5-alkyl-4-substituted 1-(pyridin-4-ylmethyl)-1H-imidazol-2-ylthio]acetamides (12a-f), and 2-[5-benzyl-1-(benzyloxymethyl)-4-isopropyl-1H-imidazol-2-ylthio]acetamides (14a-l) from their corresponding amino acids through a sequence of reactions: Dakin-West reaction, hydrolysis, condensation with thiocyanate derivatives, alkylation with 2-iodoacetamide and ethyl chloroacetate, and coupling with 4-pyridylmethyl chloride, ethoxymethyl chloride and benzyloxymethyl chloride. All the synthesized compounds were screened for their activity against HIV-1 (wild type) and some of them (especially Capravirine like structures) were found active.

Crystallography and the design of anti-AIDS drugs: conformational flexibility and positional adaptability are important in the design of non-nucleoside HIV-1 reverse transcriptase inhibitors

Drug resistance is a key cause of failure for treatment of HIV infection. The efficacy of non-nucleoside reverse transcriptase inhibiting (NNRTI) drugs is impaired by rapid emergence of drug-resistance mutations. A multidisciplinary effort led to the discovery of the potent NNRTIs dapivirine and etravirine, both of which are diarylpyrimidine (DAPY) derivatives. Systematic structural and molecular modeling studies of HIV-1 reverse transcriptase (RT)/NNRTI complexes revealed different modes of inhibitor binding, and some of the DAPY inhibitors can bind to RT in different conformations. The torsional flexibility ("wiggling") of the inhibitors can generate numerous conformational variants and the compactness of the inhibitors permits significant repositioning and reorientation (translation and rotation) within the pocket ("jiggling"). Such adaptations appear to be critical for the ability of the diarylpyrimidine NNRTIs to retain their potency against a wide range of drug-resistant HIV-1 RTs. Exploitation of inhibitor conformational flexibility (such as torsional flexibility about strategically located chemical bonds) can be a powerful element of drug design, especially for the design of drugs that will be effective against rapidly mutating targets (which is a collection of related targets).

Quantitative Structure–activity Relationship Analysis of Pyridinone HIV-1 Reverse Transcriptase Inhibitors using the k Nearest Neighbor Method and QSAR-based Database Mining

We have developed quantitative structure-activity relationship (QSAR) models for 44 non-nucleoside HIV-1 reverse transcriptase inhibitors (NNRTIs) of the pyridinone derivative type. The k nearest neighbor (kNN) variable selection approach was used. This method utilizes multiple descriptors such as molecular connectivity indices, which are derived from two-dimensional molecular topology. The modeling process entailed extensive validation including the randomization of the target property (Y-randomization) test and the division of the dataset into multiple training and test sets to establish the external predictive power of the training set models. QSAR models with high internal and external accuracy were generated, with leave-one-out cross-validated R2 (q2) values ranging between 0.5 and 0.8 for the training sets and R2 values exceeding 0.6 for the test sets. The best models with the highest internal and external predictive power were used to search the National Cancer Institute database. Derivatives of the pyrazolo[3,4-d]pyrimidine and phenothiazine type were identified as promising novel NNRTIs leads. Several candidates were docked into the binding pocket of nevirapine with the AutoDock (version 3.0) software. Docking results suggested that these types of compounds could be binding in the NNRTI binding site in a similar mode to a known non-nucleoside inhibitor nevirapine.

Conformational changes in HIV-1 reverse transcriptase induced by nonnucleoside reverse transcriptase inhibitor binding

Nonnucleoside reverse transcriptase inhibitors (NNRTI) are a group of small hydrophobic compounds with diverse structures that specifically inhibit HIV-1 reverse transcriptase (RT). NNRTIs interact with HIV-1 RT by binding to a single site on the p66 subunit of the p66/p51 heterodimeric enzyme, termed the NNRTI-binding pocket (NNRTI-BP). This binding interaction results in both short-range and long-range distortions of RT structure. In this article, we review the structural, computational and experimental evidence of the NNRTI-induced conformational changes in HIV-1 RT and relate them to the mechanism by which these compounds inhibit HIV-1 reverse transcription.

Peptide models XLV: Conformational properties of N-formyl-L-methioninamide and its relevance to methionine in proteins

The conformational space of the most biologically significant backbone folds of a suitable methionine peptide model was explored by density functional computational method. Using a medium [6-31G(d)] and a larger basis set [6-311++G(2d,2p)], the systematic exploration of low-energy backbone structures restricted for the "L-region" in the Ramachandran map of N-formyl-L-methioninamide results in conformers corresponding to the building units of an extended backbone structure (betaL), an inverse gamma-turn (gammaL), or a right-handed helical structure (alphaL). However, no poly-proline II type (epsilonL) fold was found, indicating that this conformer has no intrinsic stability, and highlighting the effect of molecular environment in stabilizing this backbone structure. This is in agreement with the abundance of the epsilonL-type backbone conformation of methionine found in proteins. Stability properties (DeltaE) and distinct backbone-side-chain interactions support the idea that specific intramolecular contacts are operative in the selection of the lowest energy conformers. Apart from the number of different folds, all stable conformers are within a 10 kcal x mol(-1) energy range, indicating the highly flexible behavior of methionine. This conformational feature can be important in supporting catalytic processes, facilitating protein folding and dimerization via metal ion binding. In both of the biological examples discussed (HIV-1 reverse transcriptase and PcoC copper-resistant protein), the conformational properties of Met residues were found to be of key importance. Spatial proximity to other types of residues or the same type of residue seems to be crucial for the structural integrity of a protein, whether Met is buried or exposed.

Flexible docking of pyridinone derivatives into the non-nucleoside inhibitor binding site of HIV-1 reverse transcriptase

Potent non-nucleoside reverse transcriptase inhibitors (NNRTIs) of the pyridinone derivative type were docked into nine NNRTIs binding pockets of HIV-1 reverse transcriptase (RT) structures. The docking results indicate that pyridinone analogues adopt a butterfly conformation and share the same binding mode as the crystal inhibitors in the pocket geometries of nevirapine, 1051U91, 9-Cl-TIBO, Cl-alpha-APA, efavirenz, UC-781, and S-1153. The results are in agreement with the data concerning mutational and structure-activity relationships available for pyridinone analogues and aid in the understanding, at the molecular level, of the biological response of published hybrid pyridinone molecules. Strategies to design further pyridinone derivatives active against RT containing mutations are discussed.

Structure of HIV-1 reverse transcriptase bound to an inhibitor active against mutant reverse transcriptases resistant to other nonnucleoside inhibitors

We have determined the crystal structure of the HIV type 1 reverse transcriptase complexed with CP-94,707, a new nonnucleoside reverse transcriptase inhibitor (NNRTI), to 2.8-A resolution. In addition to inhibiting the wild-type enzyme, this compound inhibits mutant enzymes that are resistant to inhibition by nevirapine, efavirenz, and delaviridine. In contrast to other NNRTI complexes where tyrosines 181 and 188 are pointing toward the enzyme active site, the binding pocket in this complex has the tyrosines pointing the opposite direction, as in the unliganded protein structure, to accommodate CP-94,707. This conformation of the pocket has not been observed previously in NNRTI complexes and substantially alters the shape and surface features that are available for interactions with the inhibitor. One ring of CP-94,707 makes extensive stacking interactions with tryptophan 229, one of the few residues in the NNRTI-binding pocket that cannot readily mutate to give rise to drug resistance. In this conformation of the pocket, mutations of tyrosines 181 and 188 are less likely to disrupt inhibitor binding. Modeling the asparagine mutation of lysine 103 shows that a hydrogen bond between it and tyrosine 188 could form as readily in the CP-94,707 complex as it does in the apoenzyme structure, providing an explanation for the activity of this inhibitor against this clinically important mutant.

Relative free energy of binding and binding mode calculations of HIV-1 RT inhibitors based on dock-MM-PB/GS

Tetrahydroimidazo-[4,5,l-jk][1,4]-benzodiazepin-2-(1H)-one (TIBO) derivatives are important nonnucleoside human immunodeficiency virus-1 reverse transcriptase inhibitors (NNRTI). Several TIBO derivatives have shown high potency to inhibit reverse transcriptase (RT) and one (Tivirapine) has entered into clinical trials. The free energy of binding (FEB) is a numerical way to express the binding affinity of a ligand to its receptor and has been applied in screening candidates in rational drug design. In this work, the FEB of 42 TIBOs in RT was studied. Relative FEB is expressed in the form of a linear combination of vdW, electrostatic, solvation, and nonpolar solvation energy terms. The predicted FEB activity of the TIBOs studied has a good correlation (r(2) = 0.8680, q(2) = 0.8298) with respect to the experimental activity (pIC(50)). Based on the data reported here, the Finite Difference Poisson Boltzmann with a Gaussian Smooth Dielectric Constant Function method (PB/GS) solvation energy term is very important in predicting the binding affinity of TIBOs in RT. In summary, the Dock-Molecular Mechanics (MM)-PB/GS method is a promising technique in predicting ligand/receptor binding affinity and it can be used to screen relatively large sets of molecules in a reasonable amount of computer time.

Effect of stereo and regiochemistry towards wild and multidrug resistant HIV-1 virus: viral potency of chiral PETT derivatives

Chiral derivatives of several substituted halopyridyl and thiazolyl PETT compounds were synthesized as non-nucleoside inhibitors of the reverse transcriptase (RT) enzyme of the human immunodeficiency virus (HIV-1). Molecular modeling studies indicated that because of the asymmetric geometry of the non-nucleoside inhibitors (NNRTI) binding pocket, the "R" stereoisomers would fit the NNRTI binding pocket of the HIV-1 RT much better than the corresponding "S" stereoisomers, as reflected by their 10(4)-fold lower K(i) values. The "R" stereoisomers of several PETT derivatives inhibited the recombinant RT in vitro with lower IC(50) values than their enantiomers. The active compounds were further evaluated for their ability to inhibit HIV-1 replication in human peripheral blood mononuclear cells (PBMCs). All the "R" isomers again showed potent anti-HIV activity and inhibited the replication of the HIV-1 strains HTLV(IIIB) in PBMCs at nanomolar concentrations whereas their enantiomers were less potent. The lead compounds for the respective groups were further tested against A17 (NNRTI-resistant, Y181C mutant RT), and A17Var (NNI-resistant Y181C +/- K103N mutant RT) as well as multidrug resistant viral strains. The results indicated that the lead compounds were several logs more potent than the standard NNRTI drug nevirapine. Structure-activity relationship among the derivatives showed preference of pyridyl unit with halo substitutions primarily at 5-position demonstrating the importance of both the stereochemistry as well as regiochemistry. Our data provides experimental evidence that the stereochemistry and the regiochemistry of non-nucleoside inhibitors can profoundly affect their anti-HIV activity.

Towards a MIP-based alignment and docking in computer-aided drug design

Structural alignment of ligands in their biological conformation is a crucial step in the building of pharmacophoric models in structure-based drug design. In addition, docking algorithms are limited in some cases by the quality of the scoring functions and the limited flexibility of the environment that the different programs allow. On the other hand, GRID molecular interaction potentials (MIPs) have been used for a long time in 3D-QSAR studies. However, in most of these studies the alignment of the molecules is performed on the basis of geometrical or physico-chemical criteria that differ from the MIPs used in the partial least squares statistical analysis. We have previously developed a method to use the same scoring function for the molecular alignment and for 3D-QSAR studies. This methodology, based on the use of GRID potentials, consists in the weighted averaging of similarities of the relevant MIPs of the molecules to be aligned. Here we present a method to obtain the weights for the different GRID probes in the average based on the structural information on protein-ligand complexes for relevant systems. The method, implemented in MIPSIM, is shown to yield good accuracy in the prediction of the alignments for two systems: a set of three inhibitors of dihydrofolate reductase and a set of fifteen non-nucleoside HIV-1 reverse transcriptase inhibitors (NNRTIs). The smooth GRID potentials are shown to capture the flexible character of the active site, as opposed to traditional docking scoring energy functions.

Nonnucleoside inhibitor binding affects the interactions of the fingers subdomain of human immunodeficiency virus type 1 reverse transcriptase with DNA

Site-directed photoaffinity cross-linking experiments were performed by using human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) mutants with unique cysteine residues at several positions (i.e., positions 65, 67, 70, and 74) in the fingers subdomain of the p66 subunit. Since neither the introduction of the unique cysteine residues into the fingers nor the modification of the SH groups of these residues with photoaffinity cross-linking reagents caused a significant decrease in the enzymatic activities of RT, we were able to use this system to measure distances between specific positions in the fingers domain of RT and double-stranded DNA. HIV-1 RT is quite flexible. There are conformational changes associated with binding of the normal substrates and nonnucleoside RT inhibitors (NNRTIs). Cross-linking was used to monitor intramolecular movements associated with binding of an NNRTI either in the presence or in the absence of an incoming deoxynucleoside triphosphate (dNTP). Binding an incoming dNTP at the polymerase active site decreased the efficiency of cross-linking but caused only modest changes in the preferred positions of cross-linking. This finding suggests that the fingers of p66 are closer to an extended template in the "open" configuration of the enzyme with the fingers away from the active site than in the closed configuration with the fingers in direct contact with the incoming dNTP. NNRTI binding caused increased cross-linking in experiments with diazirine reagents (especially with a diazirine reagent with a longer linker) and a moderate shift in the preferred sites of interaction with the template. Cross-linking occurred closer to the polymerase active site for RTs modified at positions 70 and 74. The effects of NNRTI binding were more pronounced in the absence of a bound dNTP; pretreatment of HIV-1 RT with an NNRTI reduced the effect of dNTP binding. These observations can be explained if the binding of NNRTI causes a decrease in the flexibility in the fingers subdomain of RT-NNRTI complex and a decrease in the distance from the fingers to the template extension.

The structure of the RNA-dependent RNA polymerase from bovine viral diarrhea virus establishes the role of GTP in de novo initiation

The bovine viral diarrhea virus (BVDV) RNA-dependent RNA polymerase can initiate RNA replication by a de novo mechanism without a primer. The structure of BVDV polymerase, determined to 2.9-A resolution, contains a unique N-terminal domain, in addition to the fingers, palm, and thumb domains common to other polymerases. The structure of BVDV polymerase complexed with GTP, which is required for de novo (primer-independent) initiation, shows that GTP binds adjacent to the initiation NTP, suggesting that the GTP mimics a vestigial RNA product. Comparison of five monomers in two different crystal forms showed conformational changes in the fingertip region and in the thumb domain that may help to translocate the RNA template and product strands during elongation. The putative binding sites of previously reported BVDV inhibitors are also discussed.

Crystal Structures of HIV-1 Reverse Transcriptases Mutated at Codons 100, 106 and 108 and Mechanisms of Resistance to Non-nucleoside Inhibitors

Leu100Ile, Val106Ala and Val108Ile are mutations in HIV-1 reverse transcriptase (RT) that are observed in the clinic and give rise to resistance to certain non-nucleoside inhibitors (NNRTIs) including the first-generation drug nevirapine. In order to investigate structural mechanisms of resistance for different NNRTI classes we have determined six crystal structures of mutant RT-inhibitor complexes. Val108 does not have direct contact with nevirapine in wild-type RT and in the RT(Val108Ile) complex the biggest change observed is at the distally positioned Tyr181 which is > 8 A from the mutation site. Thus in contrast to most NNRTI resistance mutations RT(Val108Ile) appears to act via an indirect mechanism which in this case is through alterations of the ring stacking interactions of the drug particularly with Tyr181. Shifts in side-chain and inhibitor positions compared to wild-type RT are observed in complexes of nevirapine and the second-generation NNRTI UC-781 with RT(Leu100Ile) and RT(Val106Ala), leading to perturbations in inhibitor contacts with Tyr181 and Tyr188. Such perturbations are likely to be a factor contributing to the greater loss of binding for nevirapine compared to UC-781 as, in the former case, a larger proportion of binding energy is derived from aromatic ring stacking of the inhibitor with the tyrosine side-chains. The differing resistance profiles of first and second generation NNRTIs for other drug resistance mutations in RT may also be in part due to this indirect mechanism.

A pharmacophore docking algorithm and its application to the cross-docking of 18 HIV-NNRTI's in their binding pockets

The docking of small molecules into the binding site of a target protein is an important but difficult step in structure-based drug design. The performance of a docking algorithm is usually evaluated by re-docking ligands into their native binding sites. We have explored the cross-docking of 18 HIV-NNRTIs (non-nucleoside inhibitors of HIV reverse transcriptase) of which the ligand-protein structure has been determined: each of the 18 ligands was docked into each of the 18 binding sites. The docking algorithms studied are an energy-based simulated annealing algorithm and a novel pharmacophore docking algorithm. It turns out that the energy-based docking of the ligands into non-native pockets is far less successful than the docking into their native pockets. The results can be improved by using explicit pharmacophore information, and by docking a ligand into a panel of protein structures and selecting the ligand-protein combination with the lowest interaction energy as the final result.

In silico structure-based design of a potent, mutation resilient, small peptide inhibitor of HIV-1 reverse transcriptase

A crucial step in the replication of HIV-1 is the conversion of its single-stranded RNA to double-stranded DNA, which is catalyzed by the virally encoded reverse transcriptase (RT). The latter is therefore a key target for the development of anti-HIV drugs. Currently approved anti-RT drugs fall into two main classes: (i). nucleoside analog inhibitors which are incorporated into the primer strand in their metabolically activated triphosphate forms, causing termination of DNA synthesis due to their 3'-deoxy configuration and (ii). the non-nucleoside inhibitors (NNIs), which are generally specific for HIV-1 RT and bind at an allosteric site approximately 10 A from the active site causing a displacement of the catalytic aspartate residues. The so-called "first generation" NNI drugs are generally susceptible to the effects of single-point mutations within RT, while more recent "second generation" NNIs, such as efavirenz, the carboxanilide UC-781 and certain quinoxalines demonstrate much greater resilience to mutations in RT. The crystal structures of the complexes of wild type and mutant RTs with first and second generation NNIs have shown that, for an inhibitor to be potent as well as mutation resilient, it should (i). make hydrogen bonds with the main chain of RT, (ii). have a large number of interactions with RT and (iii). have the ability to rearrange and adapt to a mutated NNI pocket. Based on the crystal structures of the complexes of wild type RT and Tyr188Cys mutant of RT with UC-781, we have designed a small peptide inhibitor. Docking results on this peptide using AutoDock3.0 and SYBYL 6.8.1 indicate that the peptide has a potency comparable to that of UC-781 with a retention of activity against the Tyr188Cys mutant RT. The proposed, small peptide is seen to possess all the desirable features of a potent and mutation resilient inhibitor and is hence a potential lead compound.

Steered molecular dynamics simulation on the binding of NNRTI to HIV-1 RT

HIV-1 reverse transcriptase (RT) is the primary target for anti-AIDS chemotherapy. Nonnucleoside RT inhibitors (NNRTIs) are very potent and most promising anti-AIDS drugs that specifically inhibit HIV-1 RT. The binding and unbinding processes of alpha-APA, an NNRTI, have been studied using nanosecond conventional molecular dynamics and steered molecular dynamics simulations. The simulation results show that the unbinding process of alpha-APA consists of three phases based on the position of alpha-APA in relation to the entrance of the binding pocket. When alpha-APA is bound in the binding pocket, the hydrophobic interactions between HIV-1 RT and alpha-APA dominate the binding; however, the hydrophilic interactions (both direct and water-bridged hydrogen bonds) also contribute to the stabilizing forces. Whereas Tyr-181 makes significant hydrophobic interactions with alpha-APA, Tyr-188 forms a strong hydrogen bond with the acylamino group (N14) of alpha-APA. These two residues have very flexible side chains and appear to act as two "flexible clamps" discouraging alpha-APA to dissociate from the binding pocket. At the pocket entrance, two relatively inflexible residues, Val-179 and Leu-100, gauge the openness of the entrance and form the bottleneck of the inhibitor-unbinding pathway. Two special water molecules at the pocket entrance appear to play important roles in inhibitor recognition of binding and unbinding. These water molecules form water bridges between the polar groups of the inhibitor and the residues around the entrance, and between the polar groups of the inhibitor themselves. The water-bridged interactions not only induce the inhibitor to adopt an energetically favorable conformation so the inhibitor can pass through the pocket entrance, but also stabilize the binding of the inhibitor in the pocket to prevent the inhibitor's dissociation. The complementary steered molecular dynamics and conventional molecular dynamics simulation results strongly support the hypothesis that NNRTIs inhibit HIV-1 RT polymerization activity by enlarging the DNA-binding cleft and restricting the flexibility and mobility of the p66 thumb subdomain that are believed to be essential during DNA translocation and polymerization.

Structural requirements for potent anti-human immunodeficiency virus (HIV) and sperm-immobilizing activities of cyclohexenyl thiourea and urea non-nucleoside inhibitors of HIV-1 reverse transcriptase

The current pandemic of sexually transmitted human immunodeficiency virus (HIV)/acquired immunodeficiency syndrome (AIDS) has created an urgent need for a new type of microbicide, one that is both a spermicide and a virucide. In a systematic effort to identify a non-detergent-type antiviral spermicide, we have rationally designed and synthesized a series of cyclohexenyl thiourea (CHET) nonnucleoside inhibitors (NNIs) of HIV-1 reverse transcriptase (RT) with sperm-immobilizing activity (SIA). To gain further insight into the structural requirements for the optimal activity of these dual-function NNIs, we compared the effects of thiazolyl, benzothiazolyl, and pyridyl ring substitutions and functionalization with electron-donating and electron-withdrawing groups as well as the importance of thiourea and urea moieties of 15 heterocyclic ring-substituted NNIs. RT activity and p24 antigen production in HIV-infected peripheral blood mononuclear cells were used as markers of viral replication. Computer-assisted sperm analysis was used for evaluating SIA of CHET compounds. The rabbit model was used for evaluation of in vivo mucosal toxicity and contraceptive activity of the lead NNIs. Three CHET-NNIs with a bromo, chloro, or methyl substitution at the 5 position of the pyridyl ring exhibited potent anti-HIV activity at nanomolar concentrations (IC(50) = 3-5 nM) and SIA at micromolar concentrations (EC(50) = 45-96 micro M). The dual-function CHET-NNIs were potent inhibitors of drug-resistant HIV-1 strains with genotypic and phenotypic NNI resistance. Upon substitution of the sulfur atom of the thiourea moiety with an oxygen atom, the most striking difference noted was a 38-fold reduction in time required for 50% sperm immobilization (T(1/2)). A quantitative structure-activity relationship (QSAR) analysis was used in deriving regression equations between 20 physicochemical properties and SIA of NNIs. QSAR analysis showed that the T(1/2) values positively correlated with values for molecular refractivity (r = 0.88), hydrophobicity (r = 0.72), atomic polarizability (r = 0.70), and principal moment of inertia (r = 0.63) of spermicidal NNIs. A stepwise multiple regression model to describe the relationship of T(1/2) values with these four regressors provided excellent predictability (r = 0.93). Exposure of semen to thiourea/urea NNIs either alone or in combination at the time of artificial insemination led to marked or complete inhibition of pregnancy in rabbits as assessed by the number of embryo implants versus corpora lutea on Day 8 of pregnancy. Repeated intravaginal application of a gel-microemulsion with and without 0.5%, 1%, and 2% CHET-NNI or its urea analog either alone or in combination did not induce mucosal toxicity. We hypothesize that the gain of spermicidal function by CHET-NNIs is due to their metabolic oxidation to urea analogs by sperm. Three reaction pathways are discussed. The extremely rapid SIA of the urea analog as well as the broad-spectrum anti-HIV activity of spermicidal CHET-NNIs together with their lack of mucosal toxicity and the marked ability to reduce in vivo fertility is particularly useful for the clinical development of a dual-function spermicidal microbicide. The cyclohexenyl pyridyl NNIs, especially N-[2-(1-cyclohexenyl)ethyl] N'-[2-(5-bromopyridyl)]-thiourea in combination with the urea analog, show unique clinical potential as anti-HIV spermicides aimed at curbing the sexual transmission of multidrug-resistant HIV-1 while providing effective fertility control for women.

Destabilization of the HIV-1 reverse transcriptase dimer upon interaction with N-acyl hydrazone inhibitors

N-(4-tert-butylbenzoyl)-2-hydroxy-1-naphthaldehyde hydrazone (BBNH) inhibits both the DNA polymerase and ribonuclease H (RNase H) activities of the human immunodeficiency virus type 1 reverse transcriptase. In this study, we show that BBNH binding impacts on the stability of the human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) heterodimer. The Gibbs free energy of dimer dissociation of HIV-1 RT is decreased in the presence of increasing concentrations of BBNH, resulting in a loss in stability of 3.8 kcal mol(-1). To evaluate whether this observed phenomenon was mediated by BBNH binding to one or more sites in RT, we synthesized a variety of BBNH analogs and identified (4-t-butylbenzoyl)-2-hydroxy-1-salicylyl hydrazone (BBSH) and (4,N,N-dimethylaminobenzoyl)-2-hydroxy-1-naphthyl hydrazone as specific inhibitors of RT DNA polymerase or RT RNase H activity, respectively. Interestingly, only BBSH provided significant destabilization of the HIV-1 RT dimer. The identification of these specific inhibitors, in combination with other biochemical data, suggests a model in which two molecules of BBNH bind per RT heterodimer. In this regard, only the binding of hydrazone molecules in the DNA polymerase domain activity elicits the observed destabilization of the HIV-1 RT heterodimer.

A quantum-dynamics study of the prototropic tautomerism of guanine and its contribution to spontaneous point mutations in Escherichia coli

High-level quantum-chemical and quantum-dynamics calculations are reported on the tautomerization equilibrium and rate constants of guanine and its complexes with one and two water molecules. The results are used to estimate the fraction of guanine present in the cell during DNA synthesis as the unwanted tautomer that forms an irregular base pair with thymine, thus giving rise to a spontaneous GC --> AT point mutation. A comparison of the estimated mutation frequency with the observed frequency in Escherichia coli is used to analyze two proposed mechanisms, differing in the extent of equilibration reached in the tautomerization reaction. In the absence of water, the equilibrium concentration of tautomeric forms is relatively large, but the barrier to their formation is high. If water is present, tautomeric forms are less favored, but water molecules may serve as efficient proton conduits causing rapid tautomerization. It is tentatively concluded that the mechanism in which a high tautomerization barrier keeps the tautomeric transformation far from a state of equilibrium is more likely than a mechanism in which water and/or polymerases produce a low equilibrium concentration of the tautomeric forms.

Conformational analysis of nevirapine, a non-nucleoside HIV-1 reverse transcriptase inhibitor, based on quantum mechanical calculations

The structure and the conformational behavior of the HIV-1 reverse transcriptase inhibitor, 11-cyclopropyl-5,11dihydro-4-methyl-6H-dipyrido[3,2-b2',3'-e][1,4]diazepin-6-one (nevirapine), is investigated by semiempirical (MNDO, AMI and PM3) method, ab initio at the HF/3-21G and HF/6-31G** levels and density functional theory at the B3LYP/6-31G** level. The fully optimized structure and rotational potential of the nitrogen and carbon bond in the cyclopropyl ring were examined in detail. A similar geometrical minimum is obtained from all methods which shows an almost identical structure to the geometry of the molecule in the complex structure with HIV-1 reverse transcriptase. To get some information on the structure in solution, NMR chemical shift calculations were also performed by a density functional theory at the B3LYP/6-31G** level, using GIAO approximation. The calculated 1H-NMR and 13C-NMR spectra for the energy minimum geometry agree well with the experimental results, which indicated that the geometry of nevirapine in solution is very similar to that of the molecule in the inhibition complex. Furthermore, the obtained results are compared to the conformational studies of other non-nucleoside reverse transcriptase inhibitors and reveal a common agreement of the non-nucleoside reverse transcriptase inhibitors. The specific butterfly-like shape and conformational flexibility within the side chain of the non-nucleoside reverse transcriptase inhibitors play an important role inducing conformational change of HIV-1 reverse transcriptase structure and are essential for the association at the inhibition pocket.

Cross-linking of the fingers subdomain of human immunodeficiency virus type 1 reverse transcriptase to template-primer

Cross-linking experiments were performed with human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) mutants with unique cysteine residues at several positions (positions 65, 67, 70, and 74) in the fingers subdomain of the p66 subunit. Two approaches were used--photoaffinity cross-linking and disulfide chemical cross-linking (using an oligonucleotide that contained an N(2)-modified dG with a reactive thiol group). In the former case, cross-linking can occur to any nucleotide in either DNA strand, and in the latter case, a specific cross-link is produced between the template and the enzyme. Neither the introduction of the unique cysteine residues into the fingers nor the modification of these residues with photocross-linking reagents caused a significant decrease in the enzymatic activities of RT. We were able to use this model system to investigate interactions between specific points on the fingers domain of RT and double-stranded DNA (dsDNA). Photoaffinity cross-linking of the template to the modified RTs with Cys residues in positions 65, 67, 70, and 74 of the fingers domain of the p66 subunit was relatively efficient. Azide-modified Cys residues produced 10 to 25% cross-linking, whereas diazirine modified residues produced 5 to 8% cross-linking. Disulfide cross-linking yields were up to 90%. All of the modified RTs preferentially photocross-linked to the 5' extended template strand of the dsDNA template-primer substrate. The preferred sites of interactions were on the extended template, 5 to 7 bases beyond the polymerase active site. HIV-1 RT is quite flexible. There are conformational changes associated with substrate binding. Cross-linking was used to detect intramolecular movements associated with binding of the incoming deoxynucleoside triphosphate (dNTP). Binding an incoming dNTP at the polymerase active site decreases the efficiency of cross-linking, but causes only modest changes in the preferred positions of cross-linking. This suggests that the interactions between the fingers of p66 and the extended template involve the "open" configuration of the enzyme with the fingers away from the active site rather than the closed configuration with the fingers in direct contact with the incoming dNTP. This experimental approach can be used to measure distances between any site on the surface of the protein and an interacting molecule.

Structural mechanisms of drug resistance for mutations at codons 181 and 188 in HIV-1 reverse transcriptase and the improved resilience of second generation non-nucleoside inhibitors

Mutations at either Tyr181 or Tyr188 within HIV-1 reverse transcriptase (RT) give high level resistance to many first generation non-nucleoside inhibitors (NNRTIs) such as the anti-AIDS drug nevirapine. By comparison second generation inhibitors, for instance the drug efavirenz, show much greater resilience to these mutations. In order to understand the structural basis for these differences we have determined a series of seven crystal structures of mutant RTs in complexes with first and second generation NNRTIs as well as one example of an unliganded mutant RT. These are Tyr181Cys RT (TNK-651) to 2.4 A, Tyr181Cys RT (efavirenz) to 2.6 A, Tyr181Cys RT (nevirapine) to 3.0 A, Tyr181Cys RT (PETT-2) to 3.0 A, Tyr188Cys RT (nevirapine) to 2.6 A, Tyr188Cys RT (UC-781) to 2.6 A and Tyr188Cys RT (unliganded) to 2.8 A resolution. In the two previously published structures of HIV-1 reverse transcriptase with mutations at 181 or 188 no side-chain electron density was observed within the p66 subunit (which contains the inhibitor binding pocket) for the mutated residues. In contrast the mutated side-chains can be seen in the NNRTI pocket for all seven structures reported here, eliminating the possibility that disordering contributes to the mechanism of resistance. In the case of the second generation compounds efavirenz with Tyr181Cys RT and UC-781 with Tyr188Cys RT there are only small rearrangements of either inhibitor within the binding site compared to wild-type RT and also for the first generation compounds TNK-651, PETT-2 and nevirapine with Tyr181Cys RT. For nevirapine with the Tyr188Cys RT there is however a more substantial movement of the drug molecule. We conclude that protein conformational changes and rearrangements of drug molecules within the mutated sites are not general features of these particular inhibitor/mutant combinations. The main contribution to drug resistance for Tyr181Cys and Tyr188Cys RT mutations is the loss of aromatic ring stacking interactions for first generation compounds, providing a simple explanation for the resilience of second generation NNRTIs, as such interactions make much less significant contribution to their binding.

The Lys103Asn mutation of HIV-1 RT: a novel mechanism of drug resistance

Inhibitors of human immunodeficiency virus (HIV) reverse transcriptase (RT) are widely used in the treatment of HIV infection. Loviride (an alpha-APA derivative) and HBY 097 (a quinoxaline derivative) are two potent non-nucleoside RT inhibitors (NNRTIs) that have been used in human clinical trials. A major problem for existing anti-retroviral therapy is the emergence of drug-resistant mutants with reduced susceptibility to the inhibitors. Amino acid residue 103 in the p66 subunit of HIV-1 RT is located near a putative entrance to a hydrophobic pocket that binds NNRTIs. Substitution of asparagine for lysine at position 103 of HIV-1 RT is associated with the development of resistance to NNRTIs; this mutation contributes to clinical failure of treatments employing NNRTIs. We have determined the structures of the unliganded form of the Lys103Asn mutant HIV-1 RT and in complexes with loviride and HBY 097. The structures of wild-type and Lys103Asn mutant HIV-1 RT in complexes with NNRTIs are quite similar overall as well as in the vicinity of the bound NNRTIs. Comparison of unliganded wild-type and Lys103Asn mutant HIV-1 RT structures reveals a network of hydrogen bonds in the Lys103Asn mutant that is not present in the wild-type enzyme. Hydrogen bonds in the unliganded Lys103Asn mutant but not in wild-type HIV-1 RT are observed between (1) the side-chains of Asn103 and Tyr188 and (2) well-ordered water molecules in the pocket and nearby pocket residues. The structural differences between unliganded wild-type and Lys103Asn mutant HIV-1 RT may correspond to stabilization of the closed-pocket form of the enzyme, which could interfere with the ability of inhibitors to bind to the enzyme. These results are consistent with kinetic data indicating that NNRTIs bind more slowly to Lys103Asn mutant than to wild-type HIV-1 RT. This novel drug-resistance mechanism explains the broad cross-resistance of Lys103Asn mutant HIV-1 RT to different classes of NNRTIs. Design of NNRTIs that make favorable interactions with the Asn103 side-chain should be relatively effective against the Lys103Asn drug-resistant mutant.

Anti-HIV activity of aromatic and heterocyclic thiazolyl thiourea compounds

Several thiazolyl thiourea derivatives were designed and synthesized as non-nucleoside inhibitors (NNRTI) of HIV-1 reverse transcriptase. Six lead compounds were identified that showed subnanomolar IC50 values for the inhibition of HIV replication, were minimally toxic to human peripheral blood mononuclear cells (PBMC) with CC50 values ranging from 28 to >100 microM, and showed remarkably high selectivity indices ranging from 28,000 to >100,000. The most promising compound was N-[1-(1-furoylmethyl)]-N'-[2-(thiazolyl)]thiourea (compound 6), which showed potency against two NNRTI-resistant HIV-1 isolates (A17 and A17 variant) at nanomolar to low micromolar concentrations, exhibited much greater potency against both wild-type as well as NNRTI-resistant HIV-1 than nevirapine, delavirdine, HI-443, and HI-244, was minimally toxic to PBMC, and had a selectivity index of > 100,000. The potency and minimal cytotoxicity of these aromatic/heterocyclic thiourea compounds suggest that they may be potentially useful as anti-AIDS drugs.

Structure-based drug design of non-nucleoside inhibitors for wild-type and drug-resistant HIV reverse transcriptase

The generation of anti-HIV agents using structure-based drug design methods has yielded a number of promising non-nucleoside inhibitors (NNIs) of HIV reverse transcriptase (RT). Recent successes in identifying potent NNIs are reviewed with an emphasis on the recent trend of utilizing a computer model of HIV RT to identify space in the NNI binding pocket that can be exploited by carefully chosen functional groups predicted to interact favorably with binding pocket residues. The NNI binding pocket model was used to design potent NNIs against both wild-type RT and drug-resistant RT mutants. Molecular modeling and score functions were used to analyze how drug-resistant mutations would change the RT binding pocket shape, volume, and chemical make-up, and how these changes could affect inhibitor binding. Modeling studies revealed that for an NNI of HIV RT to be active against RT mutants such as the especially problematic Y181C RT mutant, the following features are required: (a) the inhibitor should be highly potent against wild-type RT and therefore capable of tolerating a considerable activity loss against RT mutants (i.e. a picomolar-level inhibitor against wild-type RT may still be effective against RT mutants at nanomolar concentrations), (b) the inhibitor should maximize the occupancy in the Wing 2 region of the NNI binding site of RT, and (c) the inhibitor should contain functional groups that provide favorable chemical interactions with Wing 2 residues of wild-type as well as mutant RT. Our rationally designed NNI compounds HI-236, HI-240, HI-244, HI-253, HI-443, and HI-445 combine these three features and outperform other anti-HIV agents examined.

Structural basis for the resilience of efavirenz (DMP-266) to drug resistance mutations in HIV-1 reverse transcriptase

BACKGROUND

Efavirenz is a second-generation non-nucleoside inhibitor of HIV-1 reverse transcriptase (RT) that has recently been approved for use against HIV-1 infection. Compared with first-generation drugs such as nevirapine, efavirenz shows greater resilience to drug resistance mutations within HIV-1 RT. In order to understand the basis for this resilience at the molecular level and to help the design of further-improved anti-AIDS drugs, we have determined crystal structures of efavirenz and nevirapine with wild-type RT and the clinically important K103N mutant.

RESULTS

The relatively compact efavirenz molecule binds, as expected, within the non-nucleoside inhibitor binding pocket of RT. There are significant rearrangements of the drug binding site within the mutant RT compared with the wild-type enzyme. These changes, which lead to the repositioning of the inhibitor, are not seen in the interaction with the first-generation drug nevirapine.

CONCLUSIONS

The repositioning of efavirenz within the drug binding pocket of the mutant RT, together with conformational rearrangements in the protein, could represent a general mechanism whereby certain second-generation non-nucleoside inhibitors are able to reduce the effect of drug-resistance mutations on binding potency.

Piperidinylethyl, phenoxyethyl and fluoroethyl bromopyridyl thiourea compounds with potent anti-HIV activity

Derivatives of piperidinylethyl, phenoxyethyl and fluoroethyl bromopyridyl thioureas were designed and synthesized as non-nucleoside reverse transcriptase inhibitors (NNRTIs) of HIV-1 reverse transcriptase (RT). The anti-HIV activity of these compounds was examined by determining their ability to inhibit the replication of the HIV-1 strain HTLV(IIIB) in human peripheral blood mononuclear cells. The unsubstituted parent pyridyl thiourea compound N-[2-(1-piperidine)ethyl]-N'-[2-(pyridyl)] thiourea (1) exhibited no anti-HIV activity, even at 100 microM. However, the thiourea derivatives that contain a bromo- or chloro-substituted pyridyl group, compounds 2 and 5, inhibited HIV-1 replication at nanomolar concentrations. The addition of a methyl group onto the piperidine ring significantly altered the potency of these compounds; while methyl substitution at the 3-position of the piperidine ring reduced the activity, methyl substitution at the 2-position enhanced the anti-HIV activity. The IC50 value of the lead piperidinyl compound, N-[2-(2-methylpiperidinylethyl)]-N'-[2-(5-bromopyridyl)] thiourea (4) was <0.001 microM. All three phenoxyethyl derivatives, including the unsubstituted parent phenoxyethyl pyridyl thiourea compound N-[2-(phenoxy)ethyl]-N'-[2-(pyridyl)]thiourea (8) and the bromo-/chloro-substituted phenoxyethyl halopyridyl thiourea compounds N-[2-(phenoxy)ethyl]-N'-[2-(5-chloropyridyl)]thiourea (9) and N-[2-(phenoxy)ethyl]-N'-[2-(5-bromopyridyl)]thiourea (10) exhibited potent anti-HIV activity with nanomolar IC values. The corresponding fluoroethyl halopyridyl thiourea compounds beta-fluoro[2-phenethyl]-N'[2-(5-chloropyridyl)]thiourea (11) and beta-fluoro[2-phenethyl]-N'[2-(5-bromopyridyl)]thiourea (12) inhibited HIV-1 replication in PBMC with subnanomolar IC50 values and selectivity indices >30000. Compared to the corresponding phenoxyethyl thiourea compounds 9 and 10, these compounds were >4-5-fold more active as anti-HIV agents. Notably, the lead fluorothiourea compounds 11 and 12 were both substantially more active against the NNRTI-resistant HIV strains RT-MDR (V106A) and A17 (Y181C) than nevirapine or delavirdine. Taken together, our results provide additional experimental evidence that the structural features of the 'linker unit' between the pyridyl and phenyl moieties and changes in the phenyl group of PETT-related thiourea compounds significantly affects their biological activity as NNRTIs of HIV-1 RT.

Temperature-dependent equilibrium between the open and closed conformation of the p66 subunit of HIV-1 reverse transcriptase revealed by site-directed spin labelling

X-ray crystallographic studies of human immunodeficiency virus type 1 reverse transcriptase complexed with or without substrates or inhibitors show that the heterodimeric enzyme adopts distinct conformations that differ in the orientation of the so-called thumb subdomain in the large subunit. Site-directed spin labelling of mutated residue positions W24C and K287C is applied here to determine the distances between the fingers and thumb subdomains of liganded and unliganded RT in solution. The inter-spin distances of a DNA/DNA and a pseudoknot RNA complexed reverse transcriptase in solution was found to agree with the respective crystal data of the open and closed conformations. For the unliganded reverse transcriptase a temperature-dependent equilibrium between these two states was observed. The fraction of the closed conformation decreased from 95% at 313 K to 65% at 273 K. The spectral separation between the two structures was facilitated by the use of a perdeuterated ([15)N]nitroxide methane-thiosulfonate spin label.

Novel thiourea compounds as dual-function microbicides

Sexually active women represent the fastest growing human immunodeficiency virus (HIV)/acquired immunodeficiency syndrome risk group. In an effort to develop a vaginal microbicidal contraceptive potentially capable of preventing HIV transmission as well as providing fertility control, we previously reported the synthesis of novel nonnucleoside inhibitors (NNIs) of HIV-1 reverse transcriptase with sperm-immobilizing activity (SIA). To gain further insight into the structure-function relationship controlling these two properties of NNIs, we have rationally designed and synthesized 30 novel thiourea compounds and examined them for dual-function, anti-HIV and spermicidal activity. Twelve of the 30 thiourea compounds exhibited potent anti-HIV activity in the nanomolar range (IC(50) = <1-9 nM). Nine of the 30 thiourea derivatives exhibited both anti-HIV and spermicidal activity. Among the phenyl ring-containing thioureas, the 2-fluoro (HI-240) -substituted and 2-chloro (HI-253) -substituted derivatives exhibited potent anti-HIV activity (IC(50) = <1 nM) with SIA (EC(50) = 70 microM and 147 microM). Among the alicyclic ring-containing thioureas, the 5-bromo (HI-346) and 5-chloro (HI-445) functionalized cyclohexenyl ring-substituted thioureas were the most potent dual-function spermicides (EC(50) = 42 and 57 microM), with anti-HIV activity at nanomolar range (IC(50) = 3 nM). Unlike nonoxynol-9 (N-9), none of the potent dual-function thiourea compounds were cytotoxic to normal human vaginal, ectocervical, and endocervical epithelial cells at spermicidal concentrations. We conclude that as potent anti-HIV agents with SIA and reduced cytotoxicity when compared with N-9, the phenyl-substituted and cyclohexenyl-substituted thiourea derivatives, especially compounds HI-253 (N-[2-(2-chlorophenethyl)]-N'-[2-(5-bromopyridyl)-thiourea), HI-346 (N-[2-(5-bromopyridinyl)]-N'-[2-(1-cyclohexenyl)ethyl-thiourea), and HI-445 (N-[2-(5-chloropyridinyl)]-N'-[2-(1-cyclohexenyl)ethyl-thiourea) show unique clinical potential to become the active ingredients of a vaginal contraceptive for women who are at high risk for acquiring HIV by heterosexual vaginal transmission.

Protein-nucleic acid interactions and DNA conformation in a complex of human immunodeficiency virus type 1 reverse transcriptase with a double-stranded DNA template-primer

The conformation of the DNA and the interactions of the nucleic acid with the protein in a complex of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) and 19-mer/18-mer double-stranded DNA template-primer (dsDNA) are described. The structure of this HIV-1 RT complex with dsDNA serves as a useful paradigm for studying aspects of nucleotide polymerases such as catalysis, fidelity, drug inhibition, and drug resistance. The bound dsDNA has a bend of approximately 41 degrees at the junction of an A-form region (first five base pairs near the polymerase active site) and a B-form region (the last nine base pairs toward the RNase H active site). The 41 degrees bend occurs smoothly over the four base pairs between the A-form portion and the B-form portion in the vicinity of helices alpha H and alpha I of the p66 thumb subdomain. The interactions between the dsDNA and protein primarily involve the sugar-phosphate backbone of the nucleic acid and structural elements of the palm, thumb, and RNase H of p66, and are not sequence specific. Amino acid residues from the polymerase active site region, including amino acid residues of the conserved Tyr-Met-Asp-Asp (YMDD) motif and the "primer grip," interact with 3'-terminal nucleotides of the primer strand and are involved in positioning the primer terminal nucleotide and its 3'-OH group at the polymerase active site. Amino acid residues of the "template grip" have close contacts with the template strand and aid in positioning the template strand near the polymerase active site. Helix alpha H of the p66 thumb is partly inserted into the minor groove of the dsDNA and helix alpha I is directly adjacent to the backbone of the template strand. Amino acid residues of beta 1', alpha A', alpha B', and the loop containing His539 of the RNase H domain interact with the primer strand of the dsDNA.

Development of computational and graphical tools for analysis of movement and flexibility in large molecules

We developed a computer program for the calculation and display of the difference distance matrices (DDMs) of macromolecules that has the ability to compare multiple structures simultaneously. To demonstrate its use, a data set of atoms for superimposition of the HIV-1 reverse transcriptase enzyme was defined using the coordinates for the 21 available crystal structures of this enzyme and its complexes. The DDM technique for superimposition data set generation allows selection of atoms that are invariant in all structures, is free from user bias, and represents the most accurate and precise method of producing such subsets. Comparison of this technique was made against other published methods of generating superimposition data sets, and it was found that significant differences in magnitude and trends of atom movements are observed depending on which data set was used.

Binding of the second generation non-nucleoside inhibitor S-1153 to HIV-1 reverse transcriptase involves extensive main chain hydrogen bonding

S-1153 (AG1549) is perhaps the most promising non-nucleoside inhibitor of HIV-1 reverse transcriptase currently under development as a potential anti-AIDS drug, because it has a favorable profile of resilience to many drug resistance mutations. We have determined the crystal structure of S-1153 in a complex with HIV-1 reverse transcriptase. The complex possesses some novel features, including an extensive network of hydrogen bonds involving the main chain of residues 101, 103, and 236 of the p66 reverse transcriptase subunit. Such interactions are unlikely to be disrupted by side chain mutations. The reverse transcriptase/S-1153 complex suggests different ways in which resilience to mutations in the non-nucleoside inhibitors of reverse transcriptase binding site can be achieved.

The HIV-1 reverse transcription (RT) process as target for RT inhibitors

Since the Human Immunodeficiency Virus Type 1 (HIV-1) was identified as the etiologic agent of the Acquired Immune Deficiency Syndrome (AIDS), the HIV-1 reverse transcriptase (RT) has been the subject of intensive study. The reverse transcription entails the transition of the single-stranded viral RNA into double-stranded proviral DNA, which is then integrated into the host chromosome. Therefore, the HIV-1 reverse transcriptase plays a pivotal role in the life cycle of the virus and is consequently an interesting target for anti-HIV drug therapy. In the first section, we describe the complex process of reverse transcription and the different activities involved in this process. We then highlight the structure-function relationship of the HIV-1 reverse transcriptase, which is of great importance for a better understanding of resistance development, a major problem in anti-AIDS therapies. Finally, we summarize the mechanisms of HIV resistance toward various RT inhibitors and the implications thereof for the current anti-HIV drug therapies.

Phenylethylthiazolylthiourea (PETT) non-nucleoside inhibitors of HIV-1 and HIV-2 reverse transcriptases. Structural and biochemical analyses

Most non-nucleoside reverse transcriptase (RT) inhibitors are specific for HIV-1 RT and demonstrate minimal inhibition of HIV-2 RT. However, we report that members of the phenylethylthiazolylthiourea (PETT) series of non-nucleoside reverse transcriptase inhibitors showing high potency against HIV-1 RT have varying abilities to inhibit HIV-2 RT. Thus, PETT-1 inhibits HIV-1 RT with an IC(50) of 6 nM but shows only weak inhibition of HIV-2 RT, whereas PETT-2 retains similar potency against HIV-1 RT (IC(50) of 5 nM) and also inhibits HIV-2 RT (IC(50) of 2.2 microM). X-ray crystallographic structure determinations of PETT-1 and PETT-2 in complexes with HIV-1 RT reveal the compounds bind in an overall similar conformation albeit with some differences in their interactions with the protein. To investigate whether PETT-2 could be acting at a different site on HIV-2 RT (e.g. the dNTP or template primer binding site), we compared modes of inhibition for PETT-2 against HIV-1 and HIV-2 RT. PETT-2 was a noncompetitive inhibitor with respect to the dGTP substrate for both HIV-1 and HIV-2 RTs. PETT-2 was also a noncompetitive inhibitor with respect to a poly(rC).(dG) template primer for HIV-2 RT. These results are consistent with PETT-2 binding in corresponding pockets in both HIV-1 and HIV-2 RT with amino acid sequence differences in HIV-2 RT affecting the binding of PETT-2 compared with PETT-1.

Quantitative structure-activity relationships and comparative molecular field analysis of TIBO derivatised HIV-1 reverse transcriptase inhibitors

Quantitative structure-activity relationships (QSAR) and Comparative Molecular Field Analysis (CoMFA) have been applied in order to explain the structural requirements of HIV-1 reverse transcriptase (HIV-1 RT) inhibitory activity of TIBO derivatives on the MT-4 cells. The best QSAR model is satisfactory in both statistical significance and predictive ability. The derived structural descriptors indicate the importance of electronic contributions toward the HIV-1 RT inhibition of this class of compounds. However, it could not reveal any hydrophobic influence because of high collinearity between C2 and log P variables. In order to cope with steric interaction in the correlation, 3D-QSAR was performed using CoMFA. The obtained CoMFA model shows high predictive ability, rcv2 = 0.771, and clearly demonstrates its potential in the steric feature of the molecules through contour maps, explaining a majority (81.8%) of the variance in the data. Consequently, these results can be useful in identifying the structural requirements of TIBO derivatives and helpful for better understanding the HIV-1 RT inhibition. Eventually, they provide a beneficial basis to design new and more potent inhibitors of HIV-1 RT.

Non-nucleoside reverse transcriptase inhibitors: the NNRTI boom

Non-nucleoside reverse transcriptase inhibitors (NNRTIs) are promising drugs for the treatment of HIV when used in combination with other anti-HIV drugs such as nucleoside reverse transcriptase (RT) inhibitors and protease inhibitors. The first generation of NNRTIs have, however, suffered from the rapid development of resistance. This review discusses the properties of the FDA-approved NNRTI drugs and focuses on the recent efforts being made to produce second generation inhibitors that circumvent this resistance problem.

Structure-based design of non-nucleoside reverse transcriptase inhibitors of drug-resistant human immunodeficiency virus

A computer model of reverse transcriptase (RT) from human immunodeficiency virus type 1 (HIV-1) was used to design thiourea compounds that were predicted to inhibit RT. The RT model was used to approximate how changes in binding pocket shape, volume and chemical properties resulting from residue mutations would affect inhibitor binding. Our lead compound, N-[2-(2,5-dimethoxyphenylethyl)]-N'-[2-(5-bromopyridyl)]-thi ourea (HI-236) was tested against clinically observed non-nucleoside inhibitor (NNI)-resistant mutated strains of HIV. HI-236 was more potent than trovirdine, MKC-442 and zidovudine against the drug-sensitive HIV-1 strain IIIB, 50-100 times more effective than delavirdine or nevirapine and twice as effective as our recently reported lead compound N-[2-(2-fluorophenethyl)]-N'-[2-(5-bromopyridyl)]-thiourea (HI-240) against the NNI-resistant Y181C mutant HIV-1 strain A17. HI-236 was highly effective against the multidrug-resistant HIV-1 strain RT-MDR containing multiple mutations involving the RT residues 74V, 41L, 106A and 215Y. In general, thiourea compounds such as HI-236 and HI-240 showed better inhibition of drug-resistant strains of HIV-1 than thioalkylbenzyl-pyrimidine compounds such as HI-280 and HI-281. The improved activity of thioureas against RT mutants is consistent with a structural analysis of the NNI binding pocket model of RT. The activity of HI-236 against RT-MDR was superior to that of other anti-HIV agents tested, in the following order, from high to low activity; HI-236 (IC50 5 nM), HI-240 (IC50 6 nM), trovirdine (IC50 20 nM), zidovudine (IC50 150 nM), MKC-442 (IC50 300 nM), delavirdine (IC50 400 nM) and nevirapine (IC50 5 microM).

Rational design of N-[2-(2,5-dimethoxyphenylethyl)]-N'-[2-(5-bromopyridyl)]-thiourea (HI-236) as a potent non-nucleoside inhibitor of drug-resistant human immunodeficiency virus

The novel thiourea compound N-[2-(2,5-dimethoxyphenylethyl)]-N'-[2-(5-bromopyridyl)]-thi ourea (HI-236) targeting the non-nucleoside inhibitor (NNI) binding pocket of HIV-1 reverse transcriptase (RT) was rationally designed using a computer model of the NNI binding pocket. The NNI binding pocket model takes into consideration changes in binding pocket size, shape, and changes in residue character that result from clinically-observed NNI resistance-associated mutations of HIV RT. RT assays revealed that HI-236 was not only more potent than trovirdine, MKC-442, and AZT against the drug-sensitive HIV-1 strain HTLV(IIIB), it was also 50-100 times more effective than delavirdine or nevirapine and twice as effective as our recently reported lead compound N-[2-(2-fluorophenethyl)]-N'-[2-(5-bromopyridyl)]-thiourea (HI-240) against the NNI-resistant Y181C mutant HIV-1 strain A17. Most importantly, HI-236 was highly effective against the multidrug-resistant HIV-1 strain RT-MDR with multiple mutations involving the RT residues 74V, 41L, 106A, and 215Y. The activity of HI-236 against RT-MDR was superior to that of other anti-HIV agents tested, which are listed in the following order: HI-236 (IC50: 5 nM) > HI-240 (IC50: 6 nM) > trovirdine (IC50: 20 nM) > AZT (IC50: 150 nM) > MKC-442 (IC50: 300 nM) > delavirdine (IC50: 400 nM) > nevirapine (IC50: 5 microM).

Novel derivatives of phenethyl-5-bromopyridylthiourea and dihydroalkoxybenzyloxopyrimidine are dual-function spermicides with potent anti-human immunodeficiency virus activity

Sexually active women represent the fastest growing HIV/AIDS (human immunodeficiency virus/acquired immunodeficiency syndrome) risk group. In an effort to develop a vaginal microbicidal contraceptive potentially capable of preventing HIV transmission as well as providing fertility control, we have synthesized novel non-nucleoside inhibitors (NNIs) of HIV-1 reverse transcriptase (RT) and examined them for dual-function anti-HIV and spermicidal activity. Structure-based drug design by use of a computer docking procedure for the NNI binding pocket generated from nine RT-NNI crystal structures led to the synthesis of three novel NNIs: N-[2-(2, 5-dimethoxyphenethyl)]-N'-[2-(5-bromopyridyl)]-thiourea (D-PBT); N-[2-(2-fluorophenethyl)]-N'-[2-(5-bromopyridyl)]-thiourea (F-PBT); and 5-isopropyl-2-[(methylthiomethyl)thio]-6-(benzyl)-pyrimidin-4-(1H)-on e (S-DABO). The anti-HIV activity of these NNIs was compared with that of trovirdine and virucidal/spermicide, nonoxynol-9 (N-9), by measuring viral RT activity and p24 antigen production as markers of viral replication using HTLVIIIB-infected human peripheral blood mononuclear cells (PBMCs). The effects on sperm motion kinematics and sperm membrane integrity were examined by computer-assisted sperm analysis and by confocal laser scanning microscopy (CLSM), respectively. The growth-inhibitory effects of NNI versus N-9 against normal human ectocervical and endocervical epithelial cells were tested using the MTT (3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyl tetrazolium bromide) assay. All three NNIs were potent inhibitors of purified recombinant HIV RT and abrogated HIV replication in PBMCs at nanomolar concentrations (IC50 < 1 nM) when compared with N-9 or trovirdine (IC50 values of 2.2 microM and 0.007 microM, respectively). Two NNIs, F-PBT and S-DABO, also exhibited concentration- and time-dependent spermicidal activity. The drug concentration required to inhibit sperm motility by 50% (EC50 values) for the lead compound F-PBT versus N-9 was 147 microM and 81 microM, respectively. Sperm-immobilizing activity induced by F-PBT and S-DABO was rapid (t1/2 = 7-13 min) and irreversible. Unlike that of N-9, spermicidal activity of F-PBT and S-DABO was not accompanied by loss of acrosomal membrane as detected by fluorescent-lectin binding assay and CLSM. Whereas N-9 was cytotoxic to normal human ectocervical and endocervical cells at spermicidal doses, both F-PBT and S-DABO were selectively spermicidal. We conclude that as potent anti-HIV agents with spermicidal activity and reduced cytotoxicity, F-PBT and S-DABO show unique clinical potential to become the active ingredients of a vaginal contraceptive for women who are at high risk for acquiring HIV by heterosexual vaginal transmission.

Docking experiments in the flexible non-nucleoside inhibitor binding pocket of HIV-1 reverse transcriptase

Docking experiments were undertaken using a number of published crystal structures of HIV-1 reverse transcriptase complexes with various non-nucleoside inhibitors. The docking method was validated by successfully docking each ligand, in the conformation found in the crystal structure of the complex with the enzyme, back into its binding pocket in the right orientation and position. Each ligand was then subjected to conformational searching and a database of unique low-energy conformations of all ligands established. Docking this database into each of the reverse transcriptase binding pockets showed that all inhibitors could be fitted into each different pocket, without alteration of the pocket geometry. This contradicts findings from earlier docking investigations and implies that the conformation of the binding pocket in each different complex is conserved sufficiently to allow particular uniform ligand binding modes. The inhibitor conformations selected by this docking process are mostly the same as the one the ligand adopts in its own pocket and the selected conformations and orientations exhibit an impressive degree of similarity in the arrangement of their steric and electronic features. A correlation has also been observed between inhibitor flexibility and tightness of fit into the pockets with the more flexible inhibitors achieving a tighter fit and thus fewer favourable orientations upon docking.

Telomerases

Telomerases are RNA-dependent polymerases that catalyse the synthesis of the telomeric DNA at the tips of eukaryotic chromosomes. The recent identification of the catalytic subunit of telomerases from several different species suggests that the core of the telomerase is conserved. The proposed sequence and structural homology between the telomerase catalytic subunit and reverse transcriptases, together with a wealth of genetic and biochemical information, has led to significant advances in our understanding of the mechanism by which telomerases synthesise telomeric DNA.

Perspectives of non-nucleoside reverse transcriptase inhibitors (NNRTIs) in the therapy of HIV-1 infection

Non-nucleoside reverse transcriptase inhibitors (NNRTIs) have, in addition to the nucleoside reverse transcriptase inhibitors (NRTIs) and protease inhibitors (PIs), gained a definitive place in the treatment of HIV-1 infections. Starting from the HEPT and TIBO derivatives, more than thirty structurally different classes of compounds have been identified as NNRTIs, that is compounds that are specifically inhibitory to HIV-1 replication and targeted at the HIV-1 reverse transcriptase (RT). Two NNRTIs (nevirapine and delavirdine) have been formally licensed for clinical use and several others are (or have been) in preclinical and/or clinical development [tivirapine (TIBO R-86183), loviride (alpha-APA R89439), thiocarboxanilide UC-781, HEPT derivative MKC-442, quinoxaline HBY 097, DMP 266 (efavirenz), PETT derivatives (trovirdine, PETT-4, PETT-5) and the dichlorophenylthio(pyridyl)imidazole derivative S-1153]. The NNRTIs interact with a specific 'pocket' site of HIV-1 RT that is closely associated with, but distinct from, the NRTI binding site. NNRTIs are notorious for rapidly eliciting resistance due to mutations of the amino acids surrounding the NNRTI-binding site. However, the emergence of resistant HIV strains can be circumvented if the NNRTIs, preferably in combination with other anti-HIV agents, are used from the start at sufficiently high concentrations. In vitro, this procedure has been shown to 'knock-out' virus replication and to prevent resistance from arising. In vivo, various triple-drug combinations containing NNRTIs, NRTIs and/or PIs may result in an effective viral suppression and ensuing immune recovery. However, this so-called HAART (highly active antiretroviral therapy) may also fail, and this necessitates the design of new and more effective drugs and drug cocktails.

Mutational analysis of Tyr-318 within the non-nucleoside reverse transcriptase inhibitor binding pocket of human immunodeficiency virus type I reverse transcriptase

The highly conserved Tyr-318 is part of the non-nucleoside reverse transcriptase inhibitor (NNRTI)-specific lipophilic pocket of human immunodeficiency virus type I reverse transcriptase (RT) and makes contact within 4 A with the NNRTIs in all reported RT/NNRTI complexes. Using site-directed mutagenesis, six mutant RTs were constructed bearing the mutations Y318H, Y318K, Y318L, Y318C, Y318W, and Y318F. We found that only the Y318W and Y318F mutant RTs retained substantial RT activity, whereas the catalytic activities of the Y318K, Y318C, Y318H, and Y318L RT mutants were less than 5% of the wild-type activity. The Y318F mutant RT retained substantial sensitivity to the majority of NNRTIs tested, whereas the Y318W mutant RT showed varying degrees of resistance to NNRTIs. Subunit-specific site-directed mutagenesis revealed that there was no difference in the catalytic activity or resistance/sensitivity spectrum toward NNRTIs regardless of whether the Tyr-318 mutation was introduced in both subunits or only in the p66 subunit of RT. Recombinant viruses harboring the Y318F or Y318W mutation in the RT showed a similar resistance/sensitivity pattern to NNRTIs as their corresponding 318 mutant recombinant RTs. Our findings stress a functional or structural role for Tyr-318 in wild-type RT and argue for the design of novel NNRTIs that interact more closely with this amino acid in the NNRTI-specific pocket of human immunodeficiency virus type I RT.

Structure and functional implications of the polymerase active site region in a complex of HIV-1 RT with a double-stranded DNA template-primer and an antibody Fab fragment at 2.8 A resolution

The structure of human immunodeficiency virus type 1 (HIV-1) reverse transcriptase (RT) complexed with a 19-mer/18-mer double-stranded DNA template-primer (dsDNA) and the Fab fragment of monoclonal antibody 28 (Fab28) has been refined at 2.8 A resolution. The structures of the polymerase active site and neighboring regions are described in detail and a number of novel insights into mechanisms of polymerase catalysis and drug inhibition are presented. The three catalytically essential amino acid residues (Asp110, Asp185, and Asp186) are located close to the 3' terminus of the primer strand. Observation of a hydrogen bond between the 3'-OH of the primer terminus and the side-chain of Asp185 suggests that the carboxylate of Asp185 could act as a general base in initiating the nucleophilic attack during polymerization. Nearly all of the close protein-DNA interactions involve atoms of the sugar-phosphate backbone of the nucleic acid. However, the phenoxyl side-chain of Tyr183, which is part of the conserved YMDD motif, has hydrogen-bonding interactions with nucleotide bases of the second duplex base-pair and is predicted to have at least one hydrogen bond with all Watson-Crick base-pairs at this position. Comparison of the structure of the active site region in the HIV-1 RT/dsDNA complex with all other HIV-1 RT structures suggests that template-primer binding is accompanied by significant conformational changes of the YMDD motif that may be relevant for mechanisms of both polymerization and inhibition by non-nucleoside inhibitors. Interactions of the "primer grip" (the beta12-beta13 hairpin) with the 3' terminus of the primer strand primarily involve the main-chain atoms of Met230 and Gly231 and the primer terminal phosphate. Alternative positions of the primer grip observed in different HIV-1 RT structures may be related to conformational changes that normally occur during DNA polymerization and translocation. In the vicinity of the polymerase active site, there are a number of aromatic residues that are involved in energetically favorable pi-pi interactions and may be involved in the transitions between different stages of the catalytic process. The protein structural elements primarily responsible for precise positioning of the template-primer (including the primer grip, template grip, and helices alphaH and alphaI of the p66 thumb) can be thought of functioning as a "translocation track" that guides the relative movement of nucleic acid and protein during polymerization.