Crystal structures of HIV-1 protease-inhibitor complexes

Exploration of imatinib and nilotinib-derived templates as the P2-Ligand for HIV-1 protease inhibitors: Design, synthesis, protein X-ray structural studies, and biological evaluation

Structure-based design, synthesis, X-ray structural studies, and biological evaluation of a new series of potent HIV-1 protease inhibitors are described. These inhibitors contain various pyridyl-pyrimidine, aryl thiazole or alkylthiazole derivatives as the P2 ligands in combination with darunavir-like hydroxyethylamine sulfonamide isosteres. These heterocyclic ligands are inherent to kinase inhibitor drugs, such as nilotinib and imatinib. These ligands are designed to make hydrogen bonding interactions with the backbone atoms in the S2 subsite of HIV-1 protease. Various benzoic acid derivatives have been synthesized and incorporation of these ligands provided potent inhibitors that exhibited subnanomolar level protease inhibitory activity and low nanomolar level antiviral activity. Two high resolution X-ray structures of inhibitor-bound HIV-1 protease were determined. These structures provided important ligand-binding site interactions for further optimization of this class of protease inhibitors.

1,3-Diazepine: A privileged scaffold in medicinal chemistry

Privileged structures have been widely used as effective templates for drug discovery. While benzo-1,4-diazepine constitutes the first historical example of such a structure, the 1,3 analogue is just as rich in terms of applications in medicinal chemistry. The 1,3-diazepine moiety is present in numerous biological active compounds including natural products, and is used to design compounds displaying a large range of biological activities. It is present in the clinically used anticancer compound pentostatin, in several recent FDA approved β-lactamase inhibitors (e.g., avibactam) and also in coformycin, a natural product known as a ring-expanded purine analogue displaying antiviral and anticancer activities. Several other 1,3-diazepine containing compounds have entered into clinical trials. This heterocyclic structure has been and is still widely used in medicinal chemistry to design enzyme inhibitors, GPCR ligands, and so forth. This review endeavours to highlight the main use of the 1,3-diazepine scaffold and its derivatives, and their applications in medicinal chemistry, drug design, and therapy. We will focus more particularly on the development of enzyme inhibitors incorporating this scaffold, with a strong emphasis on the molecular interactions involved in the inhibition mechanism.

Protease Inhibitors for the Treatment of HIV/AIDS: Recent Advances and Future Challenges

Acquired Immunodeficiency Syndrome (AIDS) is a chronic disease characterized by multiple life-threatening illnesses caused by a retro-virus, Human Immunodeficiency Virus (HIV). HIV infection slowly destroys the immune system and increases the risk of various other infections and diseases. Although, there is no immediate cure for HIV infection/AIDS, several drugs targeting various cruxes of HIV infection are used to slow down the progress of the disease and to boost the immune system. One of the key therapeutic strategies is Highly Active Antiretroviral Therapy (HAART) or ' AIDS cocktail' in a general sense, which is a customized combination of anti-retroviral drugs designed to combat the HIV infection. Since HAART's inception in 1995, this treatment was found to be effective in improving the life expectancy of HIV patients over two decades. Among various classes of HAART treatment regimen, Protease Inhibitors (PIs) are known to be widely used as a major component and found to be effective in treating HIV infection/AIDS. For the past several years, a variety of protease inhibitors have been reported. This review outlines the drug design strategies of PIs, chemical and pharmacological characteristics of some mechanism-based inhibitors, summarizes the recent developments in small molecule based drug discovery with HIV protease as a drug target. Further discussed are the pharmacology, PI drug resistance on HIV PR, adverse effects of HIV PIs and challenges/impediments in the successful application of HIV PIs as an important class of drugs in HAART regimen for the effective treatment of AIDS.

Identifying chemicals with potential therapy of HIV based on protein-protein and protein-chemical interaction network

Acquired immune deficiency syndrome (AIDS) is a severe infectious disease that causes a large number of deaths every year. Traditional anti-AIDS drugs directly targeting the HIV-1 encoded enzymes including reverse transcriptase (RT), protease (PR) and integrase (IN) usually suffer from drug resistance after a period of treatment and serious side effects. In recent years, the emergence of numerous useful information of protein-protein interactions (PPI) in the HIV life cycle and related inhibitors makes PPI a new way for antiviral drug intervention. In this study, we identified 26 core human proteins involved in PPI between HIV-1 and host, that have great potential for HIV therapy. In addition, 280 chemicals that interact with three HIV drugs targeting human proteins can also interact with these 26 core proteins. All these indicate that our method as presented in this paper is quite promising. The method may become a useful tool, or at least plays a complementary role to the existing method, for identifying novel anti-HIV drugs.

Accounting for ligand conformational restriction in calculations of protein-ligand binding affinities

The conformation adopted by a ligand on binding to a receptor may differ from its lowest-energy conformation in solution. In addition, the bound ligand is more conformationally restricted, which is associated with a configurational entropy loss. The free energy change due to these effects is often neglected or treated crudely in current models for predicting binding affinity. We present a method for estimating this contribution, based on perturbation theory using the quasi-harmonic model of Karplus and Kushick as a reference system. The consistency of the method is checked for small model systems. Subsequently we use the method, along with an estimate for the enthalpic contribution due to ligand-receptor interactions, to calculate relative binding affinities. The AMBER force field and generalized Born implicit solvent model is used. Binding affinities were estimated for a test set of 233 protein-ligand complexes for which crystal structures and measured binding affinities are available. In most cases, the ligand conformation in the bound state was significantly different from the most favorable conformation in solution. In general, the correlation between measured and calculated ligand binding affinities including the free energy change due to ligand conformational change is comparable to or slightly better than that obtained by using an empirically-trained docking score. Both entropic and enthalpic contributions to this free energy change are significant.

Possible allosteric interactions of monoindazole-substituted P2 cyclic urea analogues with wild-type and mutant HIV-1 protease

Our ongoing efforts to understand the difference in the binding pattern of HIV-1 protease inhibitor (HIVPI) with the wild-type and mutant HIV-1 protease (HIVPR) and to provide mechanistic insight are continued further. We report here the results of a recent quantitative structure-activity relationship (QSAR) study on monoindazole-substituted P2 analogues of cyclic urea HIVPIs. The QSAR models revealed an inverted parabolic relationship between biological activity and calculated molar refractivity (CMR). That is, biological activity first decreases with increase in CMR and at a certain minimum point (inversion point) it suddenly changes and increases with further increase in CMR. CMR is a measure of volume-dependent-polarizability and is an indication of the polar interactions between ligand and receptor. The results seem to be best rationalized by larger molecules inducing a change in a receptor unit that allows for a new mode of interaction. Similar QSAR models were also observed for the biological activity of these molecules tested against a panel of mutant viruses including mutant strains with single amino acid substitution (I84V), double amino acid substitutions (I84V/V82F), and multiple amino acid changes corresponding to mutations observed in clinical isolates of patients treated with Ritonavir((R)). Interestingly the inversion points for these mutant strains were found larger than for wild-type. The subtle but significant difference in the inversion point indicates change in the shape and size of the binding pocket. Earlier QSAR studies have shown that the correlation of biological activity with an inverted parabola is an indicative of the 'allosteric interaction' of the ligands with the receptor. This report presents a detail analysis of these observations.

A new sequence representation as applied in better specificity elucidation for human immunodeficiency virus type 1 protease

Factor analysis scales of generalized amino acid information (FASGAI) involving hydrophobicity, alpha and turn propensities, bulky properties, compositional characteristics, local flexibility, and electronic properties were derived from 516 property parameters of 20-coded amino acids, and was then employed to represent sequence structures of 746 peptides with 8 amino acid residues. Cleavage site prediction models for human immunodeficiency virus type 1 protease by linear discriminant analysis and support vector machine with radial basis function kernel were constructed to identify if they could be cleaved or not, and were further utilized to investigate the cleavage specificity. These diversified properties, including the bulky properties, secondary conformation characteristics, electronic properties, and hydrophobicity at the first, the second, the fourth, the fifth, and the sixth residue, are possibly important factors in determining HIV PR cleavage or not. Particularly, maximal positive and negative influences result from the bulky properties of different sites. Further results from analysis of variance also likely reflect that the HIV PR recognizes diversified key properties of various sites in the octameric sequences. Satisfactory results show that FASGAI can not only be used to represent sequence structures of various functional peptides, but alsoprovide a potential feasible measure for exploring relationship between protein motif sequences and their functions.

De novo ligand design to an ensemble of protein structures

We describe a combinatorial method for de novo ligand design to an ensemble of receptor structures. Receptor conformations, protonation states, and structural water molecules are considered consistently within the framework of de novo ligand design. The method relies on Monte Carlo optimization to search the space of ligand structures, conformations, and rigid-body movements as well as receptor models. The method is applied to an ensemble of HIV protease and human collagenase receptor models. Ligand structures generated de novo exhibit the correct hydrogen-bonding pattern in the core of the active site, with hydrophobic groups extending into the receptor S1 and S1' pocket space. Furthermore, it is shown that known ligands are recovered in the correct binding mode and in the native, most tightly binding receptor model.

Comparing the conformational behavior of a series of diastereomeric cyclic urea HIV-1 inhibitors using the low mode:monte carlo conformational search method

The conformational flexibility of a series of diastereomeric cyclic urea HIV-1 protease inhibitors has been examined using the Low Mode:Monte Carlo conformational search method. Force fields were validated by a comparison of the energetic ordering of the minimum energy structures on the AMBER/GBSA(water), OPLSAA/GBSA(water) and HF/6-311G/SCRF(water) surfaces. The energetic ordering of the minima on the OPLSAA /GBSA(water) surface was in better agreement with the quantum calculations than the ordering on the AMBER/GBSA(water) surface. An ensemble of low energy structures was generated using OPLSAA/GBSA(water) and used to compare the molecular shape and flexibility of each diastereomer to the experimentally determined binding affinities and crystal structures of closely related systems. The results indicate that diastereomeric solution-phase energetic stability, conformational rigidity and ability to adopt a chair conformation correlate strongly with experimental binding affinities. Rigid body docking suggests that all of the diastereomers adopt solution-phase conformations suitable for alignment with the HIV-1 protease; however, these results indicate that the binding affinities are dependent upon subtle differences in the P1/P1' and P2/P2' substituent orientations.

X-ray crystallographic structure of ABT-378 (lopinavir) bound to HIV-1 protease

The crystal structure of ABT-378 (lopinavir), bound to the active site of HIV-1 protease is described. A comparison with crystal structures of ritonavir, A-78791, and BILA-2450 shows some analogous features with previous reported compounds. A cyclic urea unit in the P(2) position of ABT-378 is novel and makes two bidentate hydrogen bonds with Asp 29 of HIV-1 protease. In addition, a previously unreported shift in the Gly 48 carbonyl position is observed. A discussion of the structural features responsible for its high potency against wild-type HIV protease is given along with an analysis of the effect of active site mutations on potency in in vitro assays.

Design, synthesis and QSAR studies on N-aryl heteroarylisopropanolamines, a new class of non-peptidic HIV-1 protease inhibitors

A series of N-aryl heteroarylisopropanolamines in which an indole or a 3-arylpyrrole moiety was linked to an aryl group through an isopropanolamine linker, were designed and synthesized as potential anti-HIV-1-PR agents. Series was tested for their ability in blocking PR activity. As a rule, indole derivatives of class 1 exhibited more potency than pyrrole analogues of class 2 while tert-butylamide substituents increased anti-PR potency. In fact, bis tert-butylamide 1e showed the highest activity with IC(50)=25 microM. Even if not very potent, a simple class of anti-PR agents, with a facile synthetic pathway was discovered. QSAR studies on isopropanolamines 1 and 2 were performed in comparison with diarylbutanols, a new class of non peptidic anti-PR agents, recently discovered by Agouron Pharmaceuticals. QSAR and CoMFA models based on 30 diarylbutanols used as a training set were developed. The obtained models were used to investigate the binding mode of the newly synthesized derivatives 1 and 2. The results of this study suggest that N-aryl heteroarylisopropanolamines bind to the PR active site similarly to the diarylbutanols of Agouron.

Rational approach to AIDS drug design through structural biology

The discovery and development of more than a dozen drugs in the past 15 years for the treatment of AIDS offer an excellent example of progress in the field of rational drug design. At this time, the principal targets are reverse transcriptase and protease, enzymes encoded by the human immunodeficiency virus. The introduction of protease inhibitors, in particular, has drastically decreased the mortality and morbidity associated with AIDS. This review presents the methods used to develop such drugs and discusses the remaining problems, such as the rapid emergence of drug resistance.

Site-specific molecular design and its relevance to pharmacogenomics and chemical biology

The emergence of the new discipline of pharmacogenomics reflects the growing convergence of chemical and genomic space. The massive information-driven growth in both computational chemistry and structural biology is leading to unprecedented opportunities in both chemical and biological design. In this paper we relate current opinion in structural biology to recent developments in computational drug design. Sequence information now permits protein structure prediction and, together with experimental protein structure determination, a complete database of ligand-binding sites and protein-protein interactions can be assembled. When aligned with site exploration and virtual screening, this information provides a foundation for structure-based pharmacogenomics. In association with chemical genomics, structure-based design will allow major new insights into a compound's biological and pharmaceutical properties.

Quantitative structure-activity relationship of some HIV-1 protease inhibitors: a Fujita-Ban type analysis

A Fujita-Ban type analysis has been made on a few series of HIV-1 (human immunodeficiency virus of type 1) protease inhibitors and the activity contributions of various substituents obtained. From these activity contributions, a compound is predicted that may have better activity than ritonavir, presently prescribed for the treatment of patients suffering from HIV-1. A few other compounds are also suggested.

CoMFA and CoMSIA 3D QSAR and docking studies on conformationally-restrained cinnamoyl HIV-1 integrase inhibitors: exploration of a binding mode at the active site

Anti-HIV (human immunodeficiency virus) drug discovery has been increasingly focusing on HIV integrase (IN) as a potential therapeutic target. This enzyme is required for the integration of reverse transcribed proviral DNA into the host cell's genome and is essential for the propagation of the HIV life cycle. Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) three-dimensional quantitative structure-activity relationship (3D QSAR) studies and docking simulations were conducted on a series of potent conformationally restrained cinnamoyl inhibitors of HIV-1 IN (Artico; et al. J. Med. Chem. 1998, 41, 3948-3960). Predictive 3D QSAR models were established using SYBYL multifit molecular alignment rule, which had conventional r(2) and cross-validated coeffiecient (q(2)) values up to 0.981 and 0.721 for CoMFA and 0.975 and 0.804 for CoMSIA, respectively. These models were validated by an external test set (Burke; et al. J. Med. Chem. 1995, 38, 4171-4178). CoMFA and CoMSIA 3D QSAR models were also derived using a molecular alignment obtained by docking the compounds into the active site of HIV IN. These latter models were comparable to multifit-derived models in terms of relative descriptor field contributions and the partial least squares (PLS) contour maps. The CoMSIA 3D QSAR models performed better than the CoMFA models. The superior performance of CoMSIA was attributed to the large contribution of hydrogen-bonding interactions to the inhibitory activity differences among the compounds. This was supported by FlexX binding energy scores that correlated well with the inhibitory activity differences between hydroxylated compounds and their corresponding methoxy or deoxy counterparts. The CoMFA and CoMSIA PLS contour maps and MOLCAD-generated active site electrostatic, lipophilicity, and hydrogen-bonding potential surface maps, as well as the docking results, were integrated to propose a binding mode for the cinnamoyl inhibitors at the active site of HIV-1 IN.

A quantitative structure-activity relationship study on some HIV-1 protease inhibitors using molecular connectivity index

A quantitative structure-activity relationship (QSAR) study has been made on two different series of tetrahydropyrimidinones acting as HIV-1 protease inhibitors. A structural parameter, the first order valence molecular connectivity index ((1)chi(v)), has been used to account for the variation in the activity. The protease inhibition activity as well as the antiviral potency of the compounds are found to be significantly correlated with (1)chi(v) of P(2)/P(2') substituents attached to the two nitrogens N1 and N3, suggesting that substituents containing less electronegative and more saturated atoms, meaning thereby the less polar or more hydrophobic substituents, will be more advantageous. Further, if P(2) and P(2') are dissimilar, the former is found to be more effective than the latter. This difference is attributed to a conformational change in the enzyme that may be more favorable to P(2) binding than to P(2') binding.

Quantitative structure-activity relationship studies on cyclic cyanoguanidines acting as HIV-1 protease inhibitors

A quantitative structure-activity relationship study has been performed on some cyclic cyanoguanidines that inhibit the enzyme HIV-1 protease (HIV-1-PR) and exhibit antiviral potency, and the results have been compared with those of cyclic urea derivatives. Both the enzyme inhibition activity and antiviral potency in cyclic cyanoguanidines as well as in cyclic urea derivatives are found to be primarily governed by hydrophobic property of substituents attached to nitrogen (P2/P2') and further enhanced by OH or NH2 group, if any, present in the substituents. However, aromatic substituents are found to be unfavourable to both the activities of cyclic cyanoguanidines but not to any activity of cyclic urea derivatives. Cyclic urea derivatives are indicated to be more potent than cyclic cyanoguanidines. A model for the interaction of cyclic cyanoguanidines with the receptor is proposed.

Structure based prediction of binding affinity of human immunodeficiency virus-1 protease inhibitors

A series of computations were performed to derive a strategy for the prediction of binding affinities of non-peptidic human immunodeficiency virus-1 (HIV-1) protease inhibitors. This paper describes the development of a 3D quantitative structure-activity relationship (3D-QSAR) methodology by using receptor information of HIV-1 protease. The docking and molecular dynamics simulations were performed on a model ligand/enzyme complex to optimize the variables involved in the generation of ligand/enzyme models. The protonation scheme of the active site aspartic acid residues of HIV-1 protease was derived from a computational study. The active site aspartate is monoprotonated with a proton placed on the OD1 atom of the ASP B25. This protocol of docking and molecular dynamics (MD) simulation was then used to derive the ligand-enzyme complexes of the molecules used in the present study. The molecular mechanics interaction descriptors were calculated from these ligand/enzyme models. A partial least squares (PLS) method was used to derive a linear correlation between the interaction descriptors and the biological activity. A good correlation was observed when the change in the energy of the ligand upon complex formation and the electrostatic contributions to the solvation energy of the ligand were included in the QSAR analysis. A highest cross-validated q2 value of 0.649 was observed. This model had a conventional r2 of 0.725, and when this model was used to predict the activity of the external test set, it produced a r2pred of 0.761. The total interaction energy was partitioned into interactions in different subsites and interactions with each of the amino acid residues of the enzyme. The PLS analysis using these descriptors helped to identify the important interactions which can be exploited for the design of HIV-1 protease inhibitors.

Structure-based drug design approaches for predicting binding affinities of HIV1 protease inhibitors

Computational assessment of the binding affinity of enzyme inhibitors prior to synthesis is an important component of computer-assisted drug design (CADD) paradigms. The free energy perturbation (FEP) methodology is the most accurate means of estimating relative binding affinities between two inhibitors. However, due to its complexity and computation-intensive nature, practical applications are restricted to analysis of structurally-related inhibitors. Accordingly, there is a need for methods that enable rapid assessment of large number of structurally-unrelated molecules in a suitably accurate manner. In this review, the FEP method is compared with regression-based methods that employ multivariate models to assess the advantages of each in the estimation of relative binding affinities of inhibitors to an enzyme. Semiquantitative predictions of relative binding free energies of human immunodeficiency virus 1 (HIV1) protease inhibitors are also presented and compared with the corresponding FEP results. The results indicate that the regression-based methods and the FEP method are useful in the semi-quantitative and quantitative assessment of relative binding affinities of enzyme inhibitors, respectively, prior to synthesis.

Inhibitors of the C(2)-symmetric HIV-1 protease: nonsymmetric binding of a symmetric cyclic sulfamide with ketoxime groups in the P2/P2' side chains

Symmetric cyclic sulfamides, substituted in the P2/P2' position with functional groups foreseen to bind preferentially to the S2/S2' subsites of HIV-1 protease, have been prepared. Despite efforts to promote a symmetric binding, the sulfamides seemed prone to bind nonsymmetrically, as deduced from X-ray crystal structure analysis of one of the most potent inhibitors, possessing ketoxime groups in the P2/P2' side chains. Ab initio calculations suggested that the nonsymmetric conformation of the cyclic sulfamide scaffold had lower energy than the corresponding symmetric, cyclic urea-like conformation.

Quantitative structure-activity relationships of some HIV-protease inhibitors

A quantitative structure-activity relationship (QSAR) study has been made on different series of cycloalkylpyranones acting as human-immunodeficiency-virus type 1 (HIV-1) protease inhibitors. The results suggest that the enzyme binding affinity of the compounds would be favoured by a cyclooctyl ring, a 3-cyclopropylphenylmethyl substituent at the pyranone ring, and a 4-CN-2-pyridine-, an N-Me-imidazole-, or a 3- or 4-CN-phenyl-sulfonamide group at the meta position of the phenyl ring of the 3-substituent.

X-ray structure and conformational dynamics of the HIV-1 protease in complex with the inhibitor SDZ283-910: agreement of time-resolved spectroscopy and molecular dynamics simulations

Based on the X-ray structure of the human immunodeficiency virus type-1 (HIV-1) protease in complex with the statine-derived inhibitor SDZ283-910, a 542 ps molecular dynamics trajectory was computed. For comparison with the 805 ps trajectory obtained for the uncomplexed enzyme, the theoretical fluorescence anisotropy decay of the unliganded protease and the inhibitor complex was calculated from the trajectories of the Trp6A/Trp6B and Trp42A/Trp42B transition dipole moments. This enabled us to directly compare the simulated data with the experimental picosecond time-resolved fluorescence data. Fitting both experimental and simulated data to the Kohlrausch-Williams-Watts (KWW) function exp(-t/tauk)beta revealed a very good agreement for the uncomplexed protease as well as for the SDZ283-910 complex. Binding of the inhibitor induced a faster decay of both the experimental and the computed protease fluorescence anisotropy decay. By this integrative approach, the atomic detail of inhibitor-induced changes in the conformational dynamics of the HIV-1 protease was experimentally verified and will be used for further inhibitor optimisation.

Three-dimensional quantitative structure-activity relationship study on cyclic urea derivatives as HIV-1 protease inhibitors: application of comparative molecular field analysis

Three-dimensional quantitative structure-activity relationship (3D-QSAR) models have been developed using comparative molecular field analysis (CoMFA) on a large data set (118 compounds) of diverse cyclic urea derivatives as protease inhibitors against the human immunodeficiency virus type 1 (HIV-1). X-ray crystal structures of HIV-1 protease bound with this class of inhibitors were used to derive the most probable bioactive conformations of the inhibitors. The enzyme active site was used as a constraint to limit the number of possible conformations that are sterically accessible. The test sets have been created keeping in mind structural diversity as well as the uniform simple statistical criteria (mean, standard deviation, high and low values) of the protease inhibitory activities of the molecules compared to the training sets. Multiple predictive models have been developed with the training sets (93 compounds in each set) and validated with the corresponding test sets (25 compounds in each set). All the models yielded high predictive correlation coefficients (q2 from 0.699 to 0.727), substantially high fitted correlation coefficients (r2 from 0.965 to 0.973), and reasonably low standard errors of estimates (S from 0. 239 to 0.265). The steric and electrostatic effects have approximately equal contributions, 45% and 55% (approximately), respectively, toward explaining protease inhibitory activities. This analysis yielded models with significant information on steric and electrostatic interactions clearly discerned by the respective coefficient contour plots when overlapped on the X-ray structure of the HIV-1 protease. The HINT CoMFA study revealed significant contribution of hydrophobicity toward protease inhibitory activity. The 3D visualization technique utilizing these contour plots as well as the receptor site geometry may significantly improve our understanding of the inhibitor-protease (HIV-1) interactions and help in designing compounds with improved activity.

Tricyclic ureas: a new class of HIV-1 protease inhibitors

A new class of tricyclic ureas containing a conformationally constrained proline was designed with the aid of molecular modeling. Efficient stereoselective intermolecular pinacol coupling represented the highlight of the synthesis. These rigid cyclic ureas are active towards HIV-1 protease, with 9 being the most potent compound (Ki = 9 nM) despite interacting with only three side chain binding pockets of HIV protease.

A quantitative structure-activity relationship study on some sulfolanes and arylthiomethanes acting as HIV-1 protease inhibitors

A quantitative structure-activity relationship (QSAR) study on some sulfolanes and arylthiomethanes acting as human immunodeficiency virus-1 (HIV-1) protease inhibitors reveals that in the case of sulfolanes an octahydropyrindene ring and a five-membered 3(S)-sulfolane ring with a hydrophobic 2-substituent (cis to 3-substituent) will be crucial for the inhibition activity. The binding of a sulfolane, which is a nonpeptidic molecule, with the enzyme is shown to partly mimic the binding of a peptidic inhibitor. The 2-substituent is found to have strong hydrophobic interaction with the receptor. Similarly, in the case of arylthiomethanes, one of the substituents of the methane is found to have strong hydrophobic interaction with the enzyme, while the aryl substituent (4-hydroxy-6-phenyl-2-oxo-2H-pyran-3-yl) is assumed to be involved in the hydrogen bondings.

Geometrically and conformationally restrained cinnamoyl compounds as inhibitors of HIV-1 integrase: synthesis, biological evaluation, and molecular modeling

Various cinnammoyl-based structures were synthesized and tested in enzyme assays as inhibitors of the HIV-1 integrase (IN). The majority of compounds were designed as geometrically or conformationally constrained analogues of caffeic acid phenethyl ester (CAPE) and were characterized by a syn disposition of the carbonyl group with respect to the vinylic double bond. Since the cinnamoyl moiety present in flavones such as quercetin (inactive on HIV-1-infected cells) is frozen in an anti arrangement, it was hoped that fixing our compounds in a syn disposition could favor anti-HIV-1 activity in cell-based assays. Geometrical and conformational properties of the designed compounds were taken into account through analysis of X-ray structures available from the Cambridge Structural Database. The polyhydroxylated analogues were prepared by reacting 3,4-bis(tetrahydropyran-2-yloxy)benzaldehyde with various compounds having active methylene groups such as 2-propanone, cyclopentanone, cyclohexanone, 1,3-diacetylbenzene, 2, 4-dihydroxyacetophenone, 2,3-dihydro-1-indanone, 2,3-dihydro-1, 3-indandione, and others. While active against both 3'-processing and strand-transfer reactions, the new compounds, curcumin included, failed to inhibit the HIV-1 multiplication in acutely infected MT-4 cells. Nevertheless, they specifically inhibited the enzymatic reactions associated with IN, being totally inactive against other viral (HIV-1 reverse transcriptase) and cellular (RNA polymerase II) nucleic acid-processing enzymes. On the other hand, title compounds were endowed with remarkable antiproliferative activity, whose potency correlated neither with the presence of catechols (possible source of reactive quinones) nor with inhibition of topoisomerases. The SARs developed for our compounds led to novel findings concerning the molecular determinants of IN inhibitory activity within the class of cinnamoyl-based structures. We hypothesize that these compounds bind to IN featuring the cinnamoyl residue C=C-C=O in a syn disposition, differently from flavone derivatives characterized by an anti arrangement about the same fragment. Certain inhibitors, lacking one of the two pharmacophoric catechol hydroxyls, retain moderate potency thanks to nonpharmacophoric fragments (i.e., a m-methoxy group in curcumin) which favorably interact with an "accessory" region of IN. This region is supposed to be located adjacent to the binding site accommodating the pharmacophoric dihydroxycinnamoyl moiety. Disruption of coplanarity in the inhibitor structure abolishes activity owing to poor shape complementarity with the target or an exceedingly high strain energy of the coplanar conformation.

Quantitative structure-activity relationship studies on cyclic urea-based HIV protease inhibitors

A quantitative structure-activity relationship (QSAR) study is described on some cyclic ureas that inhibit the enzyme HIV-1 protease (HIV-1-PR) and exhibit antiviral potency. Both the enzyme inhibition activity and the antiviral potency were found to be primarily governed by the hydrophobic property of the substituents at the nitrogens (N2/N2') of the urea. Adjacent to the nitrogens, the C1/C1'-substituents are, however, found to affect the activity (inhibition) by their molecular size. The essential binding of the ureas with the receptor is, however, through multiple hydrogen bonding, where the substituents, too, can participate in such binding if they are capable of doing so. A schematic diagram of the overall interaction of the inhibitors with the receptor is presented.

Drug discovery: past, present and future

New drug discovery from early on involved a trial-and-error approach on naturally derived materials and substances until the end of the nineteenth century. The first half of the twentieth century witnessed systematic pharmacological evaluations of both natural and synthetic compounds. However, most new drugs until the 1970s were discovered by serendipity. With the exponential development of molecular biology on one hand and computer technology on the other, it became possible from 1980 onwards to place drug discovery on a rational basis. Cloning of genes has led to the development of methodologies for specific receptor-directed and enzyme-directed drug discoveries. Advances in recombinant DNA and transgenic technologies have enabled the production of human hormonal and other endogenous biomolecules as new drugs. As we understand more about the co-ordinating and regulating powers of the cerebral cortex during the next century, especially of the frontal lobe, man may be able to use bio-feedback training to voluntarily regulate the release of neurotransmitters, hormones, and other molecules involved in the regulation of various physiological processes in health as well as in disease.

Inhibitors of HIV-1 protease: a major success of structure-assisted drug design

Retroviral protease (PR) from the human immunodeficiency virus type 1 (HIV-1) was identified over a decade ago as a potential target for structure-based drug design. This effort was very successful. Four drugs are already approved, and others are undergoing clinical trials. The techniques utilized in this remarkable example of structure-assisted drug design included crystallography, NMR, computational studies, and advanced chemical synthesis. The development of these drugs is discussed in detail. Other approaches to designing HIV-1 PR inhibitors, based on the concepts of symmetry and on the replacement of a water molecule that had been found tetrahedrally coordinated between the enzyme and the inhibitors, are also discussed. The emergence of drug-induced mutations of HIV-1 PR leads to rapid loss of potency of the existing drugs and to the need to continue the development process. The structural basis of drug resistance and the ways of overcoming this phenomenon are mentioned.

Time-resolved fluorescence anisotropy of HIV-1 protease inhibitor complexes correlates with inhibitory activity

The tryptophan time-resolved fluorescence intensity and anisotropy of the HIV-1 protease dimer is shown to be a quick and efficient method for the conformational characterization of protease inhibitor complexes. Four fluorescence lifetimes were needed to adequately describe the fluorescence decay of the two tryptophan residues, W6 and W42, per protease monomer. As a result of the wavelength dependence of the respective amplitudes, the 2.06 ns and the 4.46 ns decay constants were suggested to be the intrinsic fluorescence lifetimes of the more solvent-exposed W6 and the less exposed W42 residues, respectively. Analysis of the fluorescence anisotropy decay yielded a short correlation time of 250 ps corresponding to local chromophore motions, and a long correlation time of 12.96 ns resulting from overall rotation of the protease enzyme. Fluorescence lifetimes and rotational correlation times changed when inhibitors of the HIV-1 protease were added. The effects of 11 different inhibitors including statine-derived, hydroxyethylamine-derived, and 2 symmetrical inhibitors on the protease fluorescence dynamics were investigated. Inhibitor binding is shown to induce an increase of the mean fluorescence lifetime taumean, an increase of the short rotational correlation time phi1, as well as a decrease of the long rotational correlation time phi2. The mean rotational correlation time phimean was identified as the global dynamic parameter for a given molecular complex, which correlates with the inhibitor dissociation constant Ki, and therefore with the activity of the inhibitor.

Structure of HOE/BAY 793 complexed to human immunodeficiency virus (HIV-1) protease in two different crystal forms--structure/function relationship and influence of crystal packing

Human immunodeficiency virus 1 (HIV-1) protease is a prime target in the search for drugs to combat the AIDS virus. The enzyme functions as a C2-symmetric dimer, cleaving the gag and gag-pol viral polyproteins at distinct sites. The possession of a twofold axis passing through the active site, has led to the design of C2-symmetrical inhibitors in the form of substrate-based transition-state analogs. One of the most active compounds of this class of inhibitors is HOE/BAY 793, which contains a vicinal diol central unit [Budt, K.-H., Hansen, J., Knolle, J., Meichsner, C., Paessens, A., Ruppert, D. & Stowasser, B. & Winkler, I. (1990) European Patent application EP0428,849; Budt, K.-H., Hansen, J., Knolle, J., Meichsner, C., Ruppert, D., Paessens, A. & Stowasser B. (1993) IXth International Conference on AIDS; Budt, K.-H., Peyman, A., Hansen, J., Knolle, J., Meichsner, C., Paessens, A., Ruppert, D. & Stowasser, B. (1995) Bioorg. Med. Chem. 3, 559-571.] The structure of this inhibitor bound to HIV-1 protease, in two different crystal forms, has been solved at 0.24-nm resolution using X-ray crystallography. In both forms, the details of the inhibitor-protease interactions revealed an overall asymmetric binding mode, especially between the central diol unit and the active-site aspartates. The main binding interactions comprise several specific H-bonds and hydrophobic contacts, which rationalize many of the characteristics of the structure/activity relationship in the class of vicinal diol inhibitors. In a general analysis of the mobility of the flap regions, which cover the active site and participate directly in binding, using our structures and the HIV protease models present in the Brookhaven databank, we found that in most structures the flexibility of the flaps is limited by local crystal contacts. However, in one of the structures presented here, no significant crystal contacts to the flap regions were present, and as a result the flexibility of the inhibitor bound flaps increased significantly. This suggests that the mobility and conformational flexibility of the flap residues are important in the functioning of HIV-1 protease, and must be considered in the future design of drugs against HIV protease and in structure-based drug design in general.

Potent HIV protease inhibitors containing a novel (hydroxyethyl)amide isostere

A series of HIV protease inhibitors containing a novel (hydroxyethyl)amidosuccinoyl core has been synthesized. These peptidomimetic structures inhibit viral protease activity at low nanomolar concentrations (IC50 < 10 nM for HIV-1 protease). The inhibition constant (Ki) for inhibitor 19 was determined to be 7.5 pM against HIV-1 and 1.2 nM against HIV-2 proteases, respectively. Several compounds (19-24) inhibited HIV-1 replication in cell culture assays with 50% effective concentrations (EC50) = 3.7-35 nM. This series of inhibitors was found to exhibit poor bioavailability (< 10%) in the rat, following oral administration. The synthesis and biological properties of these compounds are discussed. In addition, an X-ray structure of one of these inhibitors (23) in complex with HIV-2 protease provides insight into the binding mode of this novel class of HIV protease inhibitors.

HIV-1 protease inhibitors containing a novel C2 symmetrical hydroxyalkylgem-diamino core structure

Two series of peptidomimetics containing a novel C2 symmetrical hydroxyalkylgem-diamino core structure were prepared, from amino acid starting materials, and evaluated as inhibitors of HIV-1 protease (HIV-1 Pr). 1,1-Diamino-3-hydroxypropane (gHse) derivatives showed weak inhibitory potency (IC50 > 10 microM). In the 1,1-diamino-2-hydroxyethane (gSer) series, a compound containing P1/P1' benzyl and P2/P2'Fmoc substituents, displayed a significant HIV-1 Pr inhibition (IC50 = 440 nM).