» Articles » PMID: 8164256

Calculation of Relative Differences in the Binding Free Energies of HIV1 Protease Inhibitors: a Thermodynamic Cycle Perturbation Approach

Overview
Journal J Med Chem
Specialty Chemistry
Date 1994 Apr 15
PMID 8164256
Citations 8
Authors
Affiliations
Soon will be listed here.
Abstract

An iterative computer-assisted drug design (CADD) method that combines molecular mechanics, dynamics, thermodynamic cycle perturbation (TCP) calculations, molecular design, synthesis, and biochemical testing of peptidomimetic inhibitors and crystallographic structure determination of the protein-inhibitor complexes has been successfully applied to the design of novel inhibitors for the HIV1 protease. The first "designer" compound in this series (I) was designed by replacing the C-terminal Val-Val methyl ester of a known hydroxyethylene inhibitor with a diphenhydramine amide derivative in which two phenyl groups fill the p2' and p3' side-chain binding pockets in the HIV1 protease. Subsequent testing showed modest inhibition (Ki = 1.67 microM). Concurrently, molecular mechanics calculations on designed analogs indicated the feasibility of replacement of a phenyl ring with an indole ring (II). Synthesis and biochemical testing resulted in better inhibition potency for II. X-ray crystal structure determination of HIV1 protease complexed with I and II provided structural information for subsequent design and TCP calculations. A TCP protocol was established and validated for the mutation of I-->II. TCP results showed a net gain of 2.1 (+/- 0.9) kcal/mol in replacing II with I, which agreed with experimental result within an error margin of 0.8 kcal/mol. TCP calculations for six other mutations (I-->III, II-->III, IV, V, VI, and VII) were performed prior to synthesis and testing. These results allowed for the prioritization of design ideas for synthesis. In all cases where experimental results are available, TCP calculations showed good agreement. These results demonstrate that the TCP approach can be used with medicinal chemistry and crystallography for screening the proposed derivatives of a lead compound prior to synthesis, thus potentially reducing the time for the discovery of new drugs.

Citing Articles

Identification of HIV inhibitors guided by free energy perturbation calculations.

Acevedo O, Ambrose Z, Flaherty P, Aamer H, Jain P, Sambasivarao S Curr Pharm Des. 2012; 18(9):1199-216.

PMID: 22316150 PMC: 3605731. DOI: 10.2174/138161212799436421.


Relative energies of binding for antibody-carbohydrate-antigen complexes computed from free-energy simulations.

Pathiaseril A, Woods R J Am Chem Soc. 2007; 122(2):331-8.

PMID: 17211491 PMC: 1764640. DOI: 10.1021/ja9914994.


Improved convergence of binding affinities with free energy perturbation: application to nonpeptide ligands with pp60src SH2 domain.

Price D, Jorgensen W J Comput Aided Mol Des. 2001; 15(8):681-95.

PMID: 11718474 DOI: 10.1023/a:1012266200343.


Protein ligand docking based on empirical method for binding affinity estimation.

Tao P, Lai L J Comput Aided Mol Des. 2001; 15(5):429-46.

PMID: 11394737 DOI: 10.1023/a:1011188704521.


Very empirical treatment of solvation and entropy: a force field derived from log Po/w.

Kellogg G, Burnett J, Abraham D J Comput Aided Mol Des. 2001; 15(4):381-93.

PMID: 11349819 DOI: 10.1023/a:1011136228678.