Marc René De Jonge

Evolution of anti-HIV drug candidates. Part 3: diarylpyrimidine (DAPY) analogues

The synthesis and anti-HIV-1 activity of a series of diarylpyrimidines (DAPYs) are described. Several members of this novel class of non-nucleoside reverse transcriptase inhibitors (NNRTIs) are extremely potent against both wild-type and a panel of clinically significant single- and double-mutant strains of HIV-1.

Evolution of anti-HIV drug candidates. Part 2: Diaryltriazine (DATA) analogues.

A synthesis program directed toward improving the stability of imidoyl thiourea based non-nucleoside reverse transcriptase inhibitors (NNRTIs) led to the discovery of diaryltriazines (DATAs), a new class of potent NNRTIs. The synthesis and anti-HIV structure-activity relationship (SAR) studies of a series of DATA derivatives are described.

A computational model of the inhibition of Mycobacterium tuberculosis ATPase by a new drug candidate R207910

Diarylquinolines (DARQs) are a new class of potent inhibitors of the ATPase of Mycobacterium tuberculosis. We have created a homology model of a binding site for this class of compounds located on the contact area of the a‐subunit (gene atpB) and c‐subunits (gene atpE) of Mycobacterium tuberculosis ATPase. The binding pocket that was identified from the analysis of the homology model is formed by 4 helices of three c‐subunits and 2 helices of the a‐subunit. The lead compound of the DARQ series, R207910, was docked into the pocket using a simulated annealing, multiple conformer, docking algorithm. Different stereoisomers were treated separately. The best docking pose for each stereoisomer was optimized by molecular dynamics simulation on the 5300 atoms of the binding region and ligand. The interaction energies in the computed complexes enable us to rank the different stereoisomers in order of interaction strength with the ATPase binding pockets. We propose that the activity of R207910 against Mycobacterium tuberculosis is based on interference of the compound with the escapement geometry of the proton transfer chain. Upon binding the compound mimicks the conserved Arg‐186 residue of the a‐subunit and interacts in its place with the conserved acidic residue Glu‐61 of the c‐subunit. This mode of action is corroborated by the good agreement between the computed interaction energies and the observed pattern of stereo‐specificity in the model of the binding region. Proteins 2007.

On the detection of multiple-binding modes of ligands to proteins, from biological, structural, and modeling data

There are several indications that a given compound or a set of related compounds can bind in different modes to a specific binding site of a protein. This is especially evident from X-ray crystallographic structures of ligand-protein complexes. The availability of multiple binding modes of a ligand in a binding site may present an advantage in drug design when simultaneously optimizing several criteria. In the case of the design of anti-HIV compounds we observed that the more active compounds that are also resilient against mutation of the non-nucleoside binding site of HIV1-reverse transcriptase make use of more binding modes than the less active and resilient compounds.

Correlations between Factors Determining the Pharmacokinetics and Antiviral Activity of HIV-1 Non-Nucleoside Reverse Transcriptase Inhibitors of the Diaryltriazine and Diarylpyrimidine Classes of Compounds

AbstractObjective: To investigate the important factors that determine the bioavailability and the antiviral activity of the diaryltriazine (DATA) and diarylpyrimidine (DAPY) non-nucleoside reverse transcriptase inhibitors (NNRTIs) of HIV-1 in animal species and humans using cell-based assays, physicochemical and comput ed parameters. Methods: This naturalistic study included 15 parameters ranging from molecular mechanics calculations to phase I clinical trials. The calculated parameters were solvent-accessible surface area (SASA), polar surface area and Gibbs free energy of solvation. Physicochemical parameters comprised lipophilicity (octanol/water partition coefficient [cLogP]), ionisation constant (pKa), solubility and aggregate radius. Cell-based assays included human colonic adenocarcinoma cell (Caco-2) permeability (transepithelial transport), drug metabolism and antiviral activity (negative logarithm of the molar effective concentration inhibiting viral replica tion by 50% [pEC50]). Exposure was tested in rats, dogs and human volunteers. Results: Of the 15 parameters, eight correlated consistently among one another. Exposure (area under the plasma concentration-time curve [AUC]) in humans correlated positively with that in rats (r = 1.00), with transepithelial transport (r = 0.83), lipophilicity (r = 0.60), ionisability (r = 0.89), hydrodynamic radius of aggregates (r = 0.66) and with antiviral activity (r = 0.61). Exposure in humans was also seen to correlate negatively with SASA (r = −0.89). No consistent correlation was found between exposure in dogs and the eight parameters. Of the 14 DATA/DAPY molecules, 11 form aggregates with radii between 34 and 100nm. Conclusions: We observed correlations between exposure in humans with expo sure in rats, transepithelial transport (Caco-2 cells), ionisability, lipophilicity, aggregate radius and SASA in the class of DATA/DAPY NNRTI compounds. The lipophilic DATA/DAPY compounds form aggregates. It can be assumed that absorption in the intestinal tract and endocytosis in infected cells of these lipophilic compounds are governed by the common phenomenon of aggregate formation. As the lymphatic system offers a pathway for intestinal uptake of aggregates, this may offer a therapeutic advantage in the treatment of HIV-1 infection. Although it was not the objective of the study, we found that the rat was a better in vivo model than the dog for the prediction of systemic exposure in this particular set of compounds.

An ant algorithm for the conformational analysis of flexible molecules.

Originally, the ant system was developed for optimization in discrete search spaces such as the traveling salesman problem. We detail our adaptation of the algorithm to optimization in the continuous search space of conformational analysis. The parameters of the algorithm were tuned using a simple test molecule, undecane, and a drug molecule, imatinib. The algorithm is further tested on four more drug or drug‐like molecules, on vitamin A and on alanine tetrapeptide.

Ab Initio potential grid based docking: From High Performance Computing to In Silico Screening

We present a new and completely parallel method for protein ligand docking. The potential of the docking target structure is obtained directly from the electron density derived through an ab initio computation. A large subregion of the crystal structure of Isocitrate Lyase, was selected as docking target. To allow the full ab initio treatment of this region special care was taken to assign optimal basis functions. The electrostatic potential is tested by docking a small charged molecule (succinate) into the binding site. The ab initio grid yields a superior result by producing the best binding orientation and position, and by recognizing it as the best. In contrast the same docking procedure, but using a classical point‐charge based potential, produces a number of additional incorrect binding poses, and does not recognize the correct pose as the best solution.