Novel non-active site inhibitor of Cryptosporidium hominis TS-DHFR identified by a virtual screen.

The essential enzyme thymidylate synthase-dihydrofolate reductase (TS-DHFR) is a validated drug target for many pathogens, but has been elusive in Cryptosporidium hominis, as active site inhibitors of the enzymes from related parasitic protozoa show decreased potency and lack of species specificity over the human enzymes. As a rational approach to discover novel inhibitors, we conducted a virtual screen of a non-active site pocket in the DHFR linker region. From this screen, we have identified and characterized a noncompetitive inhibitor, flavin mononucleotide (FMN), with micromolar potency that is selective for ChTS-DHFR versus the human enzymes. These results describe a novel allosteric pocket amenable to inhibitor targeting, and a lead compound with which to move towards potent, selective inhibitors of ChTS-DHFR.

Analyses of activity for factor Xa inhibitors based on Monte Carlo simulations.

Monte Carlo/Extended Linear Response (MC/ELR) simulations have been conducted on 60 inhibitors of human factor Xa to determine the important interactions associated with their activity. A variety of physicochemical descriptors were configurationally averaged during the course of the simulations of each inhibitor bound to factor Xa and free in water. A regression equation was then derived; it reproduces the experimental inhibition data with a correlation coefficient, r(2), of 0.74, an rms error of 0.67 kcal/mol, and an average unsigned error of 0.60 kcal/mol using only two physically reasonable descriptors. The two factors that emerged as important in determining inhibitory potential are (1) favorable van der Waals interactions between protein and ligand and (2) direct hydrogen bonding between the inhibitor and protein. The conclusions were supported with structural analyses and results of MC/free energy perturbation (FEP) calculations.

Structural and energetic analyses of the effects of the K103N mutation of HIV-1 reverse transcriptase on efavirenz analogues.

The effect of the K103N mutation of HIV-1 reverse transcriptase (RT) on the activity of efavirenz analogues was studied via Monte Carlo/free energy perturbation calculations. The relative fold resistance energies indicate that efavirenz binds to K103N RT in a manner similar to the wild-type enzyme. The improved performance of the quinazolinones against the mutant enzyme is attributed to formation of a more optimal hydrogen-bonding network with bridging water molecules between the ligands and Glu138.

Prediction of activity for nonnucleoside inhibitors with HIV-1 reverse transcriptase based on Monte Carlo simulations.

Results of Monte Carlo (MC) simulations for more than 200 nonnucleoside inhibitors of HIV-1 reverse transcriptase (NNRTIs) representing eight diverse chemotypes have been correlated with their anti-HIV activities in an effort to establish simulation protocols and methods that can be used in the development of more effective drugs. Each inhibitor was modeled in a complex with the protein and by itself in water, and potentially useful descriptors of binding affinity were collected during the MC simulations. A viable regression equation was obtained for each data set using an extended linear response approach, which yielded r(2) values between 0.54 and 0.85 and an average unsigned error of only 0.50 kcal/mol. The most common descriptors confirm that a good geometrical match between the inhibitor and the protein is important and that the net loss of hydrogen bonds with the inhibitor upon binding is unfavorable. Other physically reasonable descriptors of binding are needed on a chemotype case-by-case basis. By including descriptors in common from the individual fits, combination regressions that include multiple data sets were also developed. This procedure led to a refined "master" regression for 210 NNRTIs with an r(2) of 0.60 and a cross-validated q(2) of 0.55. The computed activities show an rms error of 0.86 kcal/mol in comparison with experiment and an average unsigned error of 0.69 kcal/mol. Encouraging results were obtained for the predictions of 27 NNRTIs, representing a new chemotype not included in the development of the regression model. Predictions for this test set using the master regression yielded a q(2) value of 0.51 and an average unsigned error of 0.67 kcal/mol. Finally, additional regression analysis reveals that use of ligand-only descriptors leads to models with much diminished predictive ability.

Effects of Arg90 Neutralization on the Enzyme-Catalyzed Rearrangement of Chorismate to Prephenate.

Chorismate mutase (CM) is an enzyme that catalyzes the Claisen rearrangement of chorismate to prephenate. In a recent effort to understand the basis for catalysis by CM, Kienhöfer and co-workers (J. Am. Chem. Soc. 2003, 125, 3206-3207) reported results on the mutation of Arg90 in Bacillus subtilis CM (BsCM) to citrulline (Cit), an isosteric but neutral arginine analogue. An ca. 10(4)-fold decrease in kcat or 5.9 kcal/mol increase in the free-energy barrier (ΔG(‡)) for the overall catalysis was observed upon mutation. In this work, attention is turned to determining the key factors that contribute to the reduced catalytic efficiency of Arg90Cit BsCM. Using a combined QM/MM Monte Carlo/Free-Energy Perturbation method, a ΔΔG(‡) value of 3.3 kcal/mol is obtained. The higher free-energy barrier for the mutant is exclusively related to inferior stabilization of the TS, particularly one of its carboxylate groups, by neutral Cit. In addition, the reaction becomes 2.0 kcal/mol more exergonic. As BsCM is limited by product release, this step contributes to the remainder of the 10(4)-fold decrease in the rate constant in going from Arg90 to Cit.

Macrophomate synthase: QM/MM simulations address the Diels-Alder versus Michael-Aldol reaction mechanism.

Macrophomate synthase (MPS) of the phytopathogenic fungus Macrophoma commelinae catalyzes the transformation of 2-pyrone derivatives to the corresponding benzoate analogues. The transformation proceeds through three separate chemical reactions, including decarboxylation of oxalacetate to produce pyruvate enolate, two C-C bond formations between 2-pyrone and pyruvate enolate that form a bicyclic intermediate, and final decarboxylation with concomitant dehydration. Although some evidence suggests that the second step of the reaction catalyzed by MPS is a Diels-Alder reaction, definite proof that the C-C bond formations occur via a concerted mechanism is still required. An alternative route for formation of the C-C bonds is a stepwise Michael-aldol reaction. In this work, mixed quantum and molecular mechanics (QM/MM) combined with Monte Carlo simulations and free-energy perturbation (FEP) calculations were performed to investigate the relative stabilities of the transition states (TS) for both reaction mechanisms. The key results are that the Diels-Alder TS model is 17.7 and 12.1 kcal/mol less stable than the Michael and aldol TSs in the stepwise route, respectively. Therefore, this work indicates that the Michael-aldol mechanism is the route used by MPS to catalyze the second step of the overall transformation, and that the enzyme is not a natural Diels-Alderase, as claimed by Ose and co-workers (Nature 2003, 422, 185-189; Acta Crystallogr. 2004, D60, 1187-1197). A modified link-atom treatment for the bonds at the QM/MM interface is also presented.

Validation of a model for the complex of HIV-1 reverse transcriptase with nonnucleoside inhibitor TMC125.

The structure for the complex of nonnucleoside inhibitor TMC125 and HIV-1 reverse transcriptase has been determined and validated through computation of resistance profiles using Monte Carlo/free-energy perturbation calculations. The good quantitative agreement between the computed and experimental anti-HIV activities for TMC125, nevirapine, and efavirenz with wild-type RT and four common mutants (L100I, K103N, Y181C, and Y188L) confirms the correctness of the predicted structure and provides insights into the improved potency of this novel NNRTI. The blue shading in the figure indicates basic residues.

Accuracy of free energies of hydration using CM1 and CM3 atomic charges.

Absolute free energies of hydration (DeltaGhyd) have been computed for 25 diverse organic molecules using partial atomic charges derived from AM1 and PM3 wave functions via the CM1 and CM3 procedures of Cramer, Truhlar, and coworkers. Comparisons are made with results using charges fit to the electrostatic potential surface (EPS) from ab initio 6-31G* wave functions and from the OPLS-AA force field. OPLS Lennard-Jones parameters for the organic molecules were used together with the TIP4P water model in Monte Carlo simulations with free energy perturbation theory. Absolute free energies of hydration were computed for OPLS united-atom and all-atom methane by annihilating the solutes in water and in the gas phase, and absolute DeltaGhyd values for all other molecules were computed via transformation to one of these references. Optimal charge scaling factors were determined by minimizing the unsigned average error between experimental and calculated hydration free energies. The PM3-based charge models do not lead to lower average errors than obtained with the EPS charges for the subset of 13 molecules in the original study. However, improvement is obtained by scaling the CM1A partial charges by 1.14 and the CM3A charges by 1.15, which leads to average errors of 1.0 and 1.1 kcal/mol for the full set of 25 molecules. The scaled CM1A charges also yield the best results for the hydration of amides including the E/Z free-energy difference for N-methylacetamide in water.

Activity predictions for efavirenz analogues with the K103N mutant of HIV reverse transcriptase.

Monte Carlo-extended linear response (MC/ELR) calculations are used to examine the binding of efavirenz analogues with the K103N mutant of HIV-1 reverse transcriptase (HIVRT). A regression equation previously reported for the wild type (WT) enzyme is shown to predict 47 experimental activities for the K103N mutant with a q(2)=0.55 and avg error of only 0.46 kcal/mol. Further analysis identifies the key features for binding to the K103N mutant: ligand flexibility, burial of hydrophobic surface area, and protein-ligand van der Waals interactions.