Symmetric fluoro-substituted diol-based HIV protease inhibitors. Ortho-fluorinated and meta-fluorinated P1/P1'-benzyloxy side groups significantly improve the antiviral activity and preserve binding efficacy

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.

Single atom changes in newly synthesized HIV protease inhibitors reveal structural basis for extreme affinity, high genetic barrier, and adaptation to the HIV protease plasticity

HIV-1 protease inhibitors (PIs), such as darunavir (DRV), are the key component of antiretroviral therapy. However, HIV-1 often acquires resistance to PIs. Here, seven novel PIs were synthesized, by introducing single atom changes such as an exchange of a sulfur to an oxygen, scission of a single bond in P2′-cyclopropylaminobenzothiazole (or -oxazole), and/or P1-benzene ring with fluorine scan of mono- or bis-fluorine atoms around DRV’s scaffold. X-ray structural analyses of the PIs complexed with wild-type Protease (PR_WT) and highly-multi-PI-resistance-associated PR_DRV^R_P51 revealed that the PIs better adapt to structural plasticity in PR with resistance-associated amino acid substitutions by formation of optimal sulfur bond and adaptation of cyclopropyl ring in the S2′-subsite. Furthermore, these PIs displayed increased cell permeability and extreme anti-HIV-1 potency compared to DRV. Our work provides the basis for developing novel PIs with high potency against PI-resistant HIV-1 variants with a high genetic barrier.

Novel radial distribution function approach in the study of point mutations: the HIV-1 protease case study

Background: Mutations are one of the engines of evolution. Under constant stress pressure, mutations can lead to the emergence of unwanted, drug-resistant entities. Methodology: The radial distribution function weighted by the number of valence shell electrons is used to design quantitative structure-activity relationship (QSAR) model relating descriptors with the inhibition constant for a series of wild-type HIV-1 protease inhibitor complexes. The residuals of complexes with mutant HIV-1 protease were correlated with the energy of the highest occupied molecular orbitals of the residues introduced to enzyme via point mutations. Conclusion: Successful identification of residues Ile3, Asp25, Val32 and Ile50 as the one whose substitution influences the inhibition constant the most, demonstrates the potential of the proposed methodology for the study of the effects of point mutations.

Effective estimation of the inhibitor affinity of HIV-1 protease via a modified LIE approach

The inhibition of the Human Immunodeficiency Virus Type 1 Protease (HIV-1 PR) can prevent the synthesis of new viruses. Computer-aided drug design (CADD) would enhance the discovery of new therapies, through which the estimation of ligand-binding affinity is critical to predict the most efficient inhibitor. A time-consuming binding free energy method would reduce the usefulness of CADD. The modified linear interaction energy (LIE) approach emerges as an appropriate protocol that performs this task. In particular, the polar interaction free energy, which is obtained via numerically resolving the linear Poisson-Boltzmann equation, plays as an important role in driving the binding mechanism of the HIV-1 PR + inhibitor complex. The electrostatic interaction energy contributes to the attraction between two molecules, but the vdW interaction acts as a repulsive factor between the ligand and the HIV-1 PR. Moreover, the ligands were found to adopt a very strong hydrophobic interaction with the HIV-1 PR. Furthermore, the results obtained corroborate the high accuracy and precision of computational studies with a large correlation coefficient value R = 0.83 and a small RMSE δ RMSE = 1.25 kcal mol-1. This method is less time-consuming than the other end-point methods, such as the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) and free energy perturbation (FEP) approaches. Overall, the modified LIE approach would provide ligand-binding affinity with HIV-1 PR accurately, precisely, and rapidly, resulting in a more efficient design of new inhibitors.

Performance of radial distribution function-based descriptors in the chemoinformatic studies of HIV-1 protease

Aim: This letter investigates the role of radial distribution function-based descriptors for in silico design of new drugs. Methodology: The multiple linear regression models for HIV-1 protease and its complexes with a series of inhibitors were constructed. A detailed analysis of major atomic contributions to the radial distribution function descriptor weighted by the number of valence shell electrons identified residues Arg8, Asp29 and residues of the catalytic triad as crucial for the correlation with the inhibition constant, together with residues Asp30 and Ile50, whose mutations are known to cause an emergence of drug resistant variants. Conclusion: This study demonstrates an easy and fast assessment of the activity of potential drugs and the derivation of structural information of their complexes with the receptor or enzyme.

DockRMSD: an open-source tool for atom mapping and RMSD calculation of symmetric molecules through graph isomorphism

Comparison of ligand poses generated by protein–ligand docking programs has often been carried out with the assumption of direct atomic correspondence between ligand structures. However, this correspondence is not necessarily chemically relevant for symmetric molecules and can lead to an artificial inflation of ligand pose distance metrics, particularly those that depend on receptor superposition (rather than ligand superposition), such as docking root mean square deviation (RMSD). Several of the commonly-used RMSD calculation algorithms that correct for molecular symmetry do not take into account the bonding structure of molecules and can therefore result in non-physical atomic mapping. Here, we present DockRMSD, a docking pose distance calculator that converts the symmetry correction to a graph isomorphism searching problem, in which the optimal atomic mapping and RMSD calculation are performed by an exhaustive and fast matching search of all isomorphisms of the ligand structure graph. We show through evaluation of docking poses generated by AutoDock Vina on the CSAR Hi-Q set that DockRMSD is capable of deterministically identifying the minimum symmetry-corrected RMSD and is able to do so without significant loss of computational efficiency compared to other methods. The open-source DockRMSD program can be conveniently integrated with various docking pipelines to assist with accurate atomic mapping and RMSD calculations, which can therefore help improve docking performance, especially for ligand molecules with complicated structural symmetry.

Effect of polarization on HIV-1protease and fluoro-substituted inhibitors binding energies by large scale molecular dynamics simulations

Molecular dynamics simulations in explicit water are carried out to study the binding of six inhibitors to HIV-1 protease (PR) for up to 700 ns using the standard AMBER force field and polarized protein-specific charge (PPC). PPC is derived from quantum mechanical calculation for protein in solution and therefore it includes electronic polarization effect. Our results show that in all six systems, the bridging water W301 drifts away from the binding pocket in AMBER simulation. However, it is very stable in all six complexes systems using PPC. Especially, intra-protease, protease-inhibitor hydrogen bonds are dynamic stabilized in MD simulation. The computed binding free energies of six complexes have a significantly linear correlation with those experiment values and the correlation coefficient is found to be 0.91 in PPC simulation. However, the result from AMBER simulation shows a weaker correlation with the correlation coefficient of −0.51 due to the lack of polarization effect. Detailed binding interactions of W301, inhibitors with PR are further analyzed and discussed. The present study provides important information to quantitative understanding the interaction mechanism of PR-inhibitor and PR-W301 and these data also emphasizes the importance of both the electronic polarization and the bridging water molecule in predicting precisely binding affinities.

Conformer generation with OMEGA: learning from the data set and the analysis of failures

We recently published a high quality validation set for testing conformer generators, consisting of structures from both the PDB and the CSD (Hawkins, P. C. D. et al. J. Chem. Inf. Model. 2010, 50, 572.), and tested the performance of our conformer generator, OMEGA, on these sets. In the present publication, we focus on understanding the suitability of those data sets for validation and identifying and learning from OMEGA's failures. We compare, for the first time we are aware of, the coverage of the applicable property spaces between the validation data sets we used and the parent compound sets to determine if our data sets adequately sample these property spaces. We also introduce the concept of torsion fingerprinting and compare this method of dissimilation to the more traditional graph-centric diversification methods we used in our previous publication. To improve our ability to programmatically identify cases where the crystallographic conformation is not well reproduced computationally, we introduce a new metric to compare conformations, RMSTanimoto. This new metric is used alongside those from our previous publication to efficiently identify reproduction failures. We find RMSTanimoto to be particularly effective in identifying failures for the smallest molecules in our data sets. Analysis of the nature of these failures, particularly those for the CSD, sheds further light on the issue of strain in crystallographic structures. Some of the residual failure cases not resolved by simple changes in OMEGA's defaults present significant challenges to conformer generation engines like OMEGA and are a source of new avenues to further improve their performance, while others illustrate the pitfalls of validating against crystallographic ligand conformations, particularly those from the PDB.

Analysis of structural water and CH···π interactions in HIV-1 protease and PTP1B complexes using a hydrogen bond prediction tool, HBPredicT

A hydrogen bond prediction tool HBPredicT is developed for detecting structural water molecules and CH···π interactions in PDB files of protein-ligand complexes. The program adds the missing hydrogen atoms to the protein, ligands, and oxygen atoms of water molecules and subsequently all the hydrogen bonds in the complex are located using specific geometrical criteria. Hydrogen bonds are classified into various types based on (i) donor and acceptor atoms, and interactions such as (ii) protein-protein, (iii) protein-ligand, (iv) protein-water, (v) ligand-water, (vi) water-water, and (vii) protein-water-ligand. Using the information in category (vii), the water molecules which form hydrogen bonds with the ligand and the protein simultaneously-the structural water-is identified and retrieved along with the associated ligand and protein residues. For CH···π interactions, the relevant portions of the corresponding structures are also extracted in the output. The application potential of this program is tested using 19 HIV-1 protease and 11 PTP1B inhibitor complexes. All the systems showed presence of structural water molecules and in several cases, the CH···π interaction between ligand and protein are detected. A rare occurrence of CH···π interactions emanating from both faces of a phenyl ring of the inhibitor is identified in HIV-1 protease 1D4L.

A new structural theme in C2-symmetric HIV-1 protease inhibitors: ortho-substituted P1/P1' side chains

In this report, the rapid syntheses of 24 novel C2-symmetric HIV-1 protease inhibitors are described. Two ortho-iodobenzyloxy containing C-terminal duplicated inhibitors served as starting materials for microwave-enhanced palladium(0)-catalyzed carbon-carbon bond forming reactions (Suzuki, Sonogashira, Heck, and Negishi). Highly potent inhibitors equipped with ortho-functionalized P1/P1' side chains as the structural theme were identified. Computational efforts were applied to study the binding mode of this class of inhibitors and to establish structure-activity relationships. The overall orientation of the inhibitors in the active site was reproduced by docking which suggested three possible conformations of the P1/P1' groups of which two seem more plausible.

Crystallographic studies on two bioisosteric analogues, N-acetyl-beta-D-glucopyranosylamine and N-trifluoroacetyl-beta-D-glucopyranosylamine, potent inhibitors of muscle glycogen phosphorylase

Structure-based inhibitor design has led to the discovery of a number of potent inhibitors of glycogen phosphorylase b (GPb), N-acyl derivatives of beta-D-glucopyranosylamine, that bind at the catalytic site of the enzyme. The first good inhibitor in this class of compounds, N-acetyl-beta-D-glucopyranosylamine (NAG) (K(i) = 32 microM), has been previously characterized by biochemical, biological and crystallographic experiments at 2.3 angstroms resolution. Bioisosteric replacement of the acetyl group by trifluoroacetyl group resulted in an inhibitor, N-trifluoroacetyl-beta-D-glucopyranosylamine (NFAG), with a K(i) = 75 microM. To elucidate the structural basis of its reduced potency, we determined the ligand structure in complex with GPb at 1.8 angstroms resolution. To compare the binding mode of N-trifluoroacetyl derivative with that of the lead molecule, we also determined the structure of GPb-NAG complex at a higher resolution (1.9 angstroms). NFAG can be accommodated in the catalytic site of T-state GPb at approximately the same position as that of NAG and stabilize the T-state conformation of the 280 s loop by making several favourable contacts to Asn284 of this loop. The difference observed in the K(i) values of the two analogues can be interpreted in terms of subtle conformational changes of protein residues and shifts of water molecules in the vicinity of the catalytic site, variations in van der Waals interaction, and desolvation effects.