Dataset Information

Interaction of HIV-1 aspartic protease with its inhibitor, by molecular dynamics and ab initio fragment molecular orbital method.

ABSTRACT: For the three complex crystal structures of HIV-1 aspartic protease (an enzyme of AIDS) with its inhibitor in the Protein Data Bank, molecular dynamics of the generalized Born surface area and the ab initio fragment molecular orbital of an ABINIT-MP calculation was performed to obtain the binding free energy, the molecular orbital energy, the interaction energy of residues with an inhibitor and the charge transfer at the active site. The inhibitors are five symmetric cyclic ureas, of which three were modelled, and an asymmetric dipeptide. The interaction energy of the inhibitor at the active sites of aspartic acid is as great as 50 kcal mol(-1), coinciding with a tetrahedral transition state. For the inhibitor with a higher affinity, charge was transferred to the inhibitor from the active site. The difference in symmetry of the inhibitor was not evident. Binding free energy corresponds to the experimental value of the binding constant, while molecular orbital energy does not always, which is considered to be an entropy effect.

SUBMITTER: Koyano K

PROVIDER: S-EPMC2394819 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BACKGROUND:The reliable and robust estimation of ligand binding affinity continues to be a challenge in drug design. Many current methods rely on molecular mechanics (MM) calculations which do not fully explain complex molecular interactions. Full quantum mechanical (QM) computation of the electronic state of protein-ligand complexes has recently become possible by the latest advances in the development of linear-scaling QM methods such as the ab initio fragment molecular orbital (FMO) method. This approximate molecular orbital method is sufficiently fast that it can be incorporated into the development cycle during structure-based drug design for the reliable estimation of ligand binding affinity. Additionally, the FMO method can be combined with approximations for entropy and solvation to make it applicable for binding affinity prediction for a broad range of target and chemotypes. RESULTS:We applied this method to examine the binding affinity for a series of published cyclin-dependent kinase 2 (CDK2) inhibitors. We calculated the binding affinity for 28 CDK2 inhibitors using the ab initio FMO method based on a number of X-ray crystal structures. The sum of the pair interaction energies (PIE) was calculated and used to explain the gas-phase enthalpic contribution to binding. The correlation of the ligand potencies to the protein-ligand interaction energies gained from FMO was examined and was seen to give a good correlation which outperformed three MM force field based scoring functions used to appoximate the free energy of binding. Although the FMO calculation allows for the enthalpic component of binding interactions to be understood at the quantum level, as it is an in vacuo single point calculation, the entropic component and solvation terms are neglected. For this reason a more accurate and predictive estimate for binding free energy was desired. Therefore, additional terms used to describe the protein-ligand interactions were then calculated to improve the correlation of the FMO derived values to experimental free energies of binding. These terms were used to account for the polar and non-polar solvation of the molecule estimated by the Poisson-Boltzmann equation and the solvent accessible surface area (SASA), respectively, as well as a correction term for ligand entropy. A quantitative structure-activity relationship (QSAR) model obtained by Partial Least Squares projection to latent structures (PLS) analysis of the ligand potencies and the calculated terms showed a strong correlation (r2 = 0.939, q2 = 0.896) for the 14 molecule test set which had a Pearson rank order correlation of 0.97. A training set of a further 14 molecules was well predicted (r2 = 0.842), and could be used to obtain meaningful estimations of the binding free energy. CONCLUSIONS:Our results show that binding energies calculated with the FMO method correlate well with published data. Analysis of the terms used to derive the FMO energies adds greater understanding to the binding interactions than can be gained by MM methods. Combining this information with additional terms and creating a scaled model to describe the data results in more accurate predictions of ligand potencies than the absolute values obtained by FMO alone.

Dataset Information

Interaction of HIV-1 aspartic protease with its inhibitor, by molecular dynamics and ab initio fragment molecular orbital method.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure