Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.     

Resolution of discordant HIV-1 protease resistance rankings using molecular dynamics simulations

The emergence of drug resistance is a major challenge for the effective treatment of HIV. In this article, we explore the application of atomistic molecular dynamics simulations to quantify the level of resistance of a patient-derived HIV-1 protease sequence to the inhibitor lopinavir. A comparative drug ranking methodology was developed to compare drug resistance rankings produced by the Stanford HIVdb, ANRS, and RegaDB clinical decision support systems. The methodology was used to identify a patient sequence for which the three rival online tools produced differing resistance rankings. Mutations at only three positions ( L10I , A71IV, and L90M ) influenced the resistance level assigned to the sequence. We use ensemble molecular dynamics simulations to elucidate the origin of these discrepancies and the mechanism of resistance. By simulating not only the full patient sequences but also systems containing the constituent mutations, we gain insight into why resistance estimates vary and the interactions between the various mutations. In the same way, we also gain valuable knowledge of the mechanistic causes of resistance. In particular, we identify changes in the relative conformation of the two beta sheets that form the protease dimer interface which suggest an explanation of the relative frequency of different amino acids observed in patients at residue 71.     

Structural and electronic properties of new fullerene derivatives and their possible application as HIV-1 protease inhibitors

Density functional theory (DFT) calculations have been carried out at the hybrid Becke 3-Lee–Yang–Parr; B3LYP/3-21G** level of theory to study two series of hydroxy-chalca-acetic acid-(4-pyrrolidin-1-yl-phenyl) ester [C₆₀-C₂H₄N-(4-XCOCH₂OH)C₆H₄] and hydroxy-chalcoacetic acid-[2-(2-hydroxy-acetylchalcanyl)-4-pyrrolidin-1-yl-phenyl] ester[C₆₀-C₂H₄N-(3,4-XCOCH₂OH)C₆H₄]. The X atom is O, S or Se for the two series. The vibrational spectra, physical, chemical, thermodynamics and Quantitative Structure Activity Relationship (QSAR) properties of the studied molecules are calculated and discussed. We have evaluated these molecules as HIV-1 protease inhibitors based on the hydrogenation interaction between the hydroxymethylcarbonyl (HMC) groups and the two aspartic acid of the HIV-1 protease active site. Results show that some of the investigated fullerene-based derivatives can be considered promising as HIV-1 protease inhibitors.     

Synthesis of a small library of non-symmetric cyclic sulfamide HIV-1 protease inhibitors

A set of 11 non-symmetric cyclic sulfamide HIV-1 protease inhibitors were synthesized and evaluated. The use of a key microwave-assisted silver(I) oxide mediated selective mono N-benzylation reaction enabled fast and straightforward synthesis. The K_i values of the new inhibitors ranged between 0.28 μM and >20 μM.     

Accurate prediction of relative binding affinities of a series of HIV-1 protease inhibitors using semi-empirical quantum mechanical charge

A major challenge in computer-aided drug design is the accurate estimation of ligand binding affinity. Here, a new approach that combines the adaptive steered molecular dynamics (ASMD) and partial atomic charges calculated by semi-empirical quantum mechanics (SQMPC), namely ASMD-SQMPC, is suggested to predict the ligand binding affinities, with 24 HIV-1 protease inhibitors as testing examples. In the ASMD-SQMPC, the relative binding free energy (ΔG) is reflected by the average maximum potential of mean force (<PMF>_max) between bound and unbound states. The correlation coefficient (R²) between the <PMF>_max and experimentally determined ΔG is 0.86, showing a significant improvement compared with the conventional ASMD (R² = 0.52). Therefore, this study provides an efficient approach to predict the relative ΔG and reveals the significance of precise partial atomic charges in the theoretical simulations.     

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.     

Rapid synthesis and in situ screening of potent HIV-1 protease dimerization inhibitors

A library of dimerization inhibitors of HIV-1 protease is described based on crosslinked interfacial peptides. The 54 component library was designed to contain two modifications to the starting structure, one each in the Northern and Southern fragments. A rapid synthesis and in situ screening method in microtiter plates was developed to facilitate the generation and evaluation of the library members. More than 90% of the doubly modified agents were more potent than their respective singly mutated parent compounds, and five of the most potent dimerization inhibitors of HIV-1 protease described to date were identified. The free energy of binding for the combined two modifications was generally found to be additive, demonstrating the predictive value of earlier libraries.     

Receptor-independent 4D-QSAR analysis of a set of norstatine derived inhibitors of HIV-1 protease

A training set of 27 norstatine derived inhibitors of HIV-1 protease, based on the 3(S)-amino-2(S)-hydroxyl-4-phenylbutanoic acid core (AHPBA), for which the -log IC(50) values were measured, was used to construct 4D-QSAR models. Five unique RI-4D-QSAR models, from two different alignments, were identified (q(2) = 0.86-0.95). These five models were used to map the atom type morphology of the lining of the inhibitor binding site at the HIV protease receptor site as well as predict the inhibition potencies of seven test set compounds for model validation. The five models, overall, predict the -log IC(50) activity values for the test set compounds in a manner consistent with their q(2) values. The models also correctly identify the hydrophobic nature of the HIV protease receptor site, and inferences are made as to further structural modifications to improve the potency of the AHPBA inhibitors of HIV protease. The finding of five unique, and nearly statistically equivalent, RI-4D-QSAR models for the training set demonstrates that there can be more than one way to fit structure-activity data even within a given QSAR methodology. This set of unique, equally good individual models is referred to as the manifold model.     

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.     

Atypical protonation states in the active site of HIV-1 protease: a computational study

The HIV protease (HIVP) is a prominent example for successful structure-based drug design. Besides its pharmaceutical impact, it is a well-studied system for which, as experimentally evidenced, protonation changes in the active site occur upon ligand binding. Therefore, it serves as an ideal candidate for a case study using our newly developed partial charge model, which was optimized toward the application of Poisson-Boltzmann based pK(a) calculations. The charge model suggests reliably experimentally determined protonation states in the active site of HIVP. Furthermore, we perform pKa calculations for two HIVP complexes with novel types of inhibitors developed and synthesized in our group. For these complexes, no experimental knowledge about the protonation states is given. For one of the compounds, containing a central pyrrolidine ring, the calculations predict that both catalytic aspartates should be deprotonated upon ligand binding.     

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.     

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.     

Fast and selective synthesis of novel cyclic sulfamide HIV-1 protease inhibitors under controlled microwave heating

A novel and highly selective silver-promoted monobenzylation method was developed to promote synthesis of nonsymmetrical sulfamide-based HIV-1 inhibitors. Microwave-accelerated palladium-catalyzed N-amide arylation- and aminocarbonylation reactions were employed for rapid and reliable compound generation. With this class of inhibitory agents, six active inhibitors were identified, the most potent inhibitor possessing a K_i-value of 20 nM.     

Design and synthesis of stereochemically defined novel spirocyclic P2-ligands for HIV-1 protease inhibitors

The synthesis of a series of stereochemically defined spirocyclic compounds and their use as novel P2-ligands for HIV-1 protease inhibitors are described. The bicyclic core of the ligands was synthesized by an efficient nBu 3SnH-promoted radical cyclization of a 1,6-enyne followed by oxidative cleavage. Structure-based design, synthesis of ligands, and biological evaluations of the resulting inhibitors are reported.