Crystal structure of a novel synthetic inhibitor of HIV-1 protease

The crystal structure of a novel non-peptidic HIV-1 protease inhibitor derived by simple solid-state dimerization of 4-aryl-1,4-dihydropyridines, reveals a strained central cage and the conformation of its phenyl, benzyl, and hydroxymethylene substituents. The polycyclic cage includes two nearly flat cyclobutane rings and four fused piperidine rings in boat conformations. The cage geometry reveals two unexpected features, namely marked distortions of the valence angles in every second piperidine and a shortening of one of the cyclobutane bonds. The molecule displays exact centrosymmetry, but the central cage and the hydroxymethylene substituents also approximate the C₂-symmetry of the target enzyme. The two independent hydroxyl groups are involved in intermolecular hydrogen bonding, one as a donor, the other as an acceptor. The disposition of the hydroxyl groups in the molecular framework is compatible with the dual role of the inhibitor in the active-site cavity of HIV-1 protease, whereby one OH group is hydrogen-bonded to the catalytic aspartates, whereas another one provides an interface to the locked flaps of the enzyme.     

Computational design of novel cyclic urea as HIV-1 protease inhibitor

A series of novel cyclic urea molecules 5,6-dihydroxy-1,3-diazepane-2,4,7-trione as HIV-1 protease inhibitors were designed using computational techniques. The designed molecules were compared with the known cyclic urea molecules by performing docking studies, calculating their ADME (Absorption, Distribution, Metabolism, and Excretion) properties and protein ligand interaction energy. These novel molecules were designed by substituting the P ₁/P′ ₁ positions (4 ^th and 7 ^th position of 1, 3-diazepan-2-one) with double bonded oxygens. This reduces the molecular weight and increases the bioavailability, indicating better ADME properties. The docking studies showed good binding affinity towards HIV-1 protease. The biological activity of these inhibitors were predicted by a model equation generated by the regression analysis between biological activity (log 1/K ⁱ ) of known inhibitors and their protein ligand interaction energy. The synthetic studies are in progress.      

Conformational analysis of TMC114, a novel HIV-1 protease inhibitor

TMC114, a potent novel HIV-1 protease inhibitor, remains active against a broad spectrum of mutant viruses. In order to bind to a variety of mutants, the compound needs to make strong, preferably backbone, interactions and have enough conformational flexibility to adapt to the changing geometry of the active site. The conformational analysis of TMC114 in the gas phase yielded 43 conformers in which five types of intramolecular H-bond interactions could be observed. All 43 conformers were subject to both rigid and flexible ligand docking in the wild-type and a triple mutant (L63P/V82T/I84V) of HIV-1 protease. The largest binding energy was calculated for the conformations that are close to the conformation observed in the X-ray complexes of TMC114 and HIV-1 protease.     

Heterocyclic HIV-1 protease inhibitors

[formula: see text] A series of simple heterocyclic HIV-1 protease inhibitors were developed on the basis of size, shape, and electronic complementarity to the active site of the enzyme. The C2-symmetric heterocycles do not contain a transition-state isostere nor are they active site directed irreversible inhibitors; thus, they represent the success of a new design strategy. The first generation heterocycles inhibit the protease in the micromolar range, whereas control compounds show no bioactivity at the same concentrations.     

Binding free energy contributions of interfacial waters in HIV-1 protease/inhibitor complexes

Water molecules are commonly observed in crystal structures of protein-ligand complexes where they mediate protein-ligand binding. It is of considerable theoretical and practical importance to determine quantitatively the individual free energy contributions of these interfacial water molecules to protein-ligand binding and to elucidate factors that influence them. The double-decoupling free energy molecular dynamics simulation method has been used to calculate the binding free energy contribution for each of the four interfacial water molecules observed in the crystal structure of HIV-1 protease complexed with KNI-272, a potent inhibitor. While two of these water molecules contribute significantly to the binding free energy, the other two have close to zero contribution. It was further observed that the protonation states of two catalytic aspartate residues, Asp25 and Asp125, strongly influence the free energy contribution of a conserved water molecule Wat301 and that different inhibitors significantly influence the free energy contribution of Wat301. Our results have important implications on our understanding of the role of interfacial water molecules in protein-ligand binding and to structure-based drug design aimed at incorporating these interfacial water molecules into ligands.     

HPLC quantification of the HIV-1 protease inhibitor saquinavir in brain and testis of mice

A rapid, reliable HPLC method with UV detection (240 nm) was developed and validated for quantitation of saquinavir in mice brain and testis. Saquinavir and the internal standard were isolated from homogenized tissue matrices using liquid-liquid extraction procedure and were then analyzed using an isocratic mobile phase by reversed-phase liquid chromatography. The lower limit of quantification was 50 ng/g for both brain and testis. A linear dynamic range of 50-5000 ng/g for both brain and testis was established. This HPLC method was validated with between-batch precision of 0.5-4.4 and 1.5-5.5% for brain and testis, respectively. The between-batch accuracy was 94.7-105.9% and 97.5-105.0% for brain and testis, respectively. The present method was applied for tissue distribution studies of the novel drug delivery systems of saquinavir in mice.     

3D-QSAR and molecular modeling of HIV-1 integrase inhibitors

Three-dimensional quantitative structure-activity relationship (3D QSAR) methods were applied on a series of inhibitors of HIV-1 integrase with respect to their inhibition of 3-processing and 3-end joining steps in vitro.The training set consisted of 27 compounds belonging to the class of thiazolothiazepines. The predictive ability of each model was evaluated using test set I consisting of four thiazolothiazepines and test set II comprised of seven compounds belonging to an entirely different structural class of coumarins. Maximum Common Substructure (MCS) based method was used to align the molecules and this was compared with other known methods of alignment. Two methods of 3D QSAR: comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) were analyzed in terms of their predictive abilities. CoMSIA produced significantly better results for all correlations. The results indicate a strong correlation between the inhibitory activity of these compounds and the steric and electrostatic fields around them. CoMSIA models with considerable internal as well as external predictive ability were obtained. A poor correlation obtained with hydrophobic field indicates that the binding of thiazolothiazepines to HIV-1 integrase is mainly enthalpic in nature. Further the most active compound of the series was docked into the active site using the crystal structure of integrase. The binding site was formed by the amino acid residues 64-67, 116, 148, 151-152, 155-156, and 159. The comparison of coefficient contour maps with the steric and electrostatic properties of the receptor shows high level of compatibility.     

Use of 3D QSAR to investigate the mode of binding of pyrazinones to HIV-1 RT

Abstract  

Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) based on the docked conformation were performed for 24 pyrazinone derivatives. All compounds were docked into the wild-type HIV-1 RT binding pocket and the lowest-energy docked configurations were used to construct the 3D QSAR models. The CoMFA and CoMSIA models enable good prediction of inhibition by the pyrazinones, with r_\textcv² r_{\text{cv}}^{2}  = 0.703 and 0.735. Results obtained from CoMFA and CoMSIA based on the docking conformation of the pyrazinones are, therefore, powerful means of elucidating the mode of binding of pyrazinones and suggesting the design of new potent NNRTIs.     

Use of 3D QSAR to investigate the mode of binding of pyrazinones to HIV-1 RT

Abstract  Comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) based on the docked conformation were performed for 24 pyrazinone derivatives. All compounds were docked into the wild-type HIV-1 RT binding pocket and the lowest-energy docked configurations were used to construct the 3D QSAR models. The CoMFA and CoMSIA models enable good prediction of inhibition by the pyrazinones, with  = 0.703 and 0.735. Results obtained from CoMFA and CoMSIA based on the docking conformation of the pyrazinones are, therefore, powerful means of elucidating the mode of binding of pyrazinones and suggesting the design of new potent NNRTIs. Graphical Abstract        

Comparison of azacyclic urea A-98881 as HIV-1 protease inhibitor with cage dimeric N-benzyl 4-(4-methoxyphenyl)-1,4- dihydropyridine as representative of a novel class of HIV-1 protease inhibitors: A molecular modeling study

The functional groups of cage dimeric N-alkyl substituted 3,5-bis(hydroxymethyl)-4-(4-methoxyphenyl)-1,4-dihydropyridines are similar to those of cyclic and azacyclic ureas that are potent inhibitors of HIV-1 protease of the dihydroxyethylene- and hydroxyethylene type, respectively. In the following study the conformity of common functional groups is investigated concerning their orientation in space as well as in the enzyme HIV-1 protease. Starting from X-ray crystal data of the centrosymmetric cage dimeric N-benzyl derivative with ester groups, the derivative with hydroxymethylene groups was built and a systematic conformational search was performed for the conformationally important torsion angles considering electrostatic and van der Waals interactions. From the huge number of conformations those comprising centrosymmetrical and C2-symmetrical energy minima were selected and minimized. The three remaining conformers were fitted to the azacyclic urea A-98881 selected from the HIV-1 protease enzyme- inhibitor complex using the centroids of the corresponding aromatic residues and additionally by the field fit option of the Advanced CoMFA module of SYBYL. Interestingly, the energetically most favourable one, which, additionally, possesses C2-symmetry like the active site cavity of HIV-1 protease, showed the best fit. Comparing the electrostatic potential (EP) of the latter with the EP of A-98881 the aromatic residues show excellent accordance. Slight differences in the extent of the EP were found in the areas of the hydroxymethylene groups of the cage dimer and the single hydroxy group as well as the urea carbonyl group of A- 98881, respectively. In order to compare the binding possibilities to the enzyme HIV-1 protease for the cage dimer and A-98881, their interaction fields with certain probes (CH3 for alkyl, NHamide, and carbonyl, O– of COO–), representing the decisive functional groups of the active site, have been calculated using GRID and projected into the enzyme placing the structures according to the position of A-98881 in the enzyme- inhibitor complex. The strongest calculated fields of the O– probe were found near Asp 25 for both structures. Another respective conformity consists in the overlap of the fields for the NHamide probe near Ile 50 and 50 for the investigated cage dimer and A-98881.     

Stereoselective synthesis of 3-aminoindan-1-ones and subsequent incorporation into HIV-1 protease inhibitors

A new method for the stereoselective synthesis of 3-aminoindan-1-ones from triflates of salicylic sulfinyl imines and ethylene glycol vinyl ether has been developed. The reaction sequence starts with a regioselective Heck reaction followed by stereoselective Lewis acid mediated annulation. Acidic cleavage of the sulfinamides produced pure (R)- and (S)-3-aminoindan-1-ones, which were successfully isolated and incorporated into active HIV-1 protease inhibitors.     

Absolute free energies of binding of peptide analogs to the HIV-1 protease from molecular dynamics simulations

The constants of binding of five peptide analogs to the active site of the HIV-1 aspartic-protease are calculated based on a novel sampling scheme that is efficient and does not introduce any approximations in addition to the energy function used to describe the system. The results agree with experiments. The squared correlation coefficient of the calculated vs. the measured values is 0.79. The sampling scheme consists of a series of molecular dynamics integrations with biases. The biases are selected based on an estimate of the probability density function of the system in a way to explore the conformational space and to reduce the statistical error in the calculated binding constants. The molecular dynamics integrations are done with a vacuum potential using a short cutoff scheme. To estimate the probability density of the simulated system, the results of the molecular dynamics integrations are combined using an extension of the weighted histogram analysis method (C. Bartels, Chem. Phys. Letters 331 (2000) 446-454). The probability density of the solvated ligand-protein system is obtained by applying a correction for the use of the short cutoffs in the simulations and by taking into account solvation with an electrostatic term and a hydrophobic term. The electrostatic part of the solvation is determined by finite difference Poisson-Boltzmann calculations; the hydrophobic part of the solvation is set proportional to the solvent accessible surface area. Setting the hydrophobic surface tension parameter equal to 8 mol(-1) K(-1) A(-2), absolute binding constants are in the muM to nM range. This is in agreement with experiments. The standard errors determined from eight repeated binding constant determinations are a factor of 14 to 411. A single determination of a binding constant is done with 499700 steps of molecular dynamics integration and 4500 finite difference Poisson-Boltzmann calculations. The simulations can be analyzed with respect to conformational changes of the active site of the HIV-1 protease or the ligands upon binding and provide information that complements experiments and can be used in the drug development process.