Biomedical subjects

Irwin D Kuntz

Publications and source records attributed to Irwin D Kuntz.

At least 19 recordsLinked to original sources

PubMed2006

Development and validation of a modular, extensible docking program: DOCK 5.

We report on the development and validation of a new version of DOCK. The algorithm has been rewritten in a modular format, which allows for easy implementation of new scoring functions, sampling methods and analysis tools. We validated the sampling algorithm with a test set of 114 protein-ligand complexes. Using an optimized parameter set, we are able to reproduce the crystal ligand pose to within 2 A of the crystal structure for 79% of the test cases using our rigid ligand docking algorithm with an average run time of 1 min per complex and for 72% of the test cases using our flexible ligand docking algorithm with an average run time of 5 min per complex. Finally, we perform an analysis of the docking failures in the test set and determine that the sampling algorithm is generally sufficient for the binding pose prediction problem for up to 7 rotatable bonds; i.e. 99% of the rigid ligand docking cases and 95% of the flexible ligand docking cases are sampled successfully. We point out that success rates could be improved through more advanced modeling of the receptor prior to docking and through improvement of the force field parameters, particularly for structures containing metal-based cofactors.

Algorithms↗

PubMed2006

Rational design of a nonpeptide general chemical scaffold for reversible inhibition of PDZ domain interactions.

Novel small molecules were designed to specifically target the ligand-binding pocket of a PDZ domain. Iterative molecular docking and modeling allowed the design of an indole scaffold 10a as a reversible inhibitor of ligand binding. The 10a scaffold inhibited the interaction between MAGI-3 and PTEN and showed cellular activities that are consistent with the inhibition of NHERF-1 function.

Adaptor Proteins, Signal Transducing↗

PubMed2005

Apolipoprotein (apo) E4 enhances amyloid beta peptide production in cultured neuronal cells: apoE structure as a potential therapeutic target.

Apolipoprotein (apo) E4 is a major risk factor for Alzheimer's disease, and many studies have suggested that apoE has isoform-specific effects on the deposition or clearance of amyloid beta (Abeta) peptides. We examined the effects of apoE isoforms on the processing of amyloid precursor protein (APP) and on Abeta production in rat neuroblastoma B103 cells stably transfected with human wild-type APP695 (B103-APP). Lipid-poor apoE4 increased Abeta production in B103-APP cells to a greater extent than lipid-poor apoE3 (60% vs. 30%) due to more pronounced stimulation of APP recycling by apoE4 than apoE3. The difference in Abeta production was abolished by preincubating the cells with the receptor-associated protein (25 nM), which blocks the low-density lipoprotein receptor-related protein (LRP) pathway, or by reducing LRP expression by small interference RNA. The differences were also attenuated by replacing Arg-61 with threonine in apoE4 or pretreating apoE4 with small molecules, both of which abolish apoE4 intramolecular domain interaction. Thus, apoE4 appears to modulate APP processing and Abeta production through both the LRP pathway and domain interaction. These findings provide insights into why apoE4 is associated with increased risk for Alzheimer's disease and may represent a potential target for drug development.

Amyloid beta-Peptides↗

PubMed2005

Evaluating the high-throughput screening computations.

The judges evaluated the submissions for the McMaster University High-Throughput Data-Mining and Docking Competition based on 3 criteria: identification of active compounds, percent enrichment, and overview of the competition. Using these metrics, 4 of the participating groups found meaningful enrichment, and 3 groups made perceptive comments about the general nature of the competition.

Computational Biology↗

PubMed2005

Hierarchical database screenings for HIV-1 reverse transcriptase using a pharmacophore model, rigid docking, solvation docking, and MM-PB/SA.

In this work, an efficient strategy was presented to search drug leads for human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) using hierarchical database screenings, which included a pharmacophore model, multiple-conformation rigid docking, solvation docking, and molecular mechanics-Poisson-Boltzmann/surface area (MM-PB/SA) sequentially. Encouraging results were achieved in searching a refined available chemical directory (ACD) database: the enrichment factor after the first three filters was estimated to be 25-fold; the hit rate for all the four filters was predicted to be 41% in a control test using 37 known HIV-1 non-nucleoside reverse transcriptase inhibitors; 10 out of 30 promising solvation-docking hits had MM-PB/SA binding free energies better than -6.8 kcal/mol and the best one, HIT15, had -17.0 kcal/mol. In conclusion, the hierarchical multiple-filter database searching strategy is an attractive strategy in drug lead exploration.

Allosteric Site↗

PubMed2005

An information theoretic approach to macromolecular modeling: I. Sequence alignments.

We are interested in applying the principles of information theory to structural biology calculations. In this article, we explore the information content of an important computational procedure: sequence alignment. Using a reference state developed from exhaustive sequences, we measure alignment statistics and evaluate gap penalties based on first-principle considerations and gap distributions. We show that there are different gap penalties for different alphabet sizes and that the gap penalties can depend on the length of the sequences being aligned. In a companion article, we examine the information content of molecular force fields.

Algorithms↗

PubMed2005

An information theoretic approach to macromolecular modeling: II. Force fields.

In this article, we explore the information content of molecular force-field calculations. We make use of exhaustive lattice models of molecular conformations and reduced alphabet sequences to determine the relative resolving power of pairwise interaction-based force fields. We find that sequence-specific interactions that operate over longer distances offer greater amounts of information than nearest-neighbor or non-sequence-specific interactions. In a companion article in this issue, we explored the information content of sequence alignment procedures and the calculation of gap penalties. Both articles have implications for protein and nucleic-acid computations.

Computational Biology↗

PubMed2004

Photochemically knocking out glutamate receptors in vivo.

AMPA (alpha-amino-3-hydroxy-5-methyl-4-isooxazole) receptors, a major subtype of ionotropic glutamate receptors (iGluRs), mediate the majority of the fast communication between neurons, and the activity-dependent trafficking of AMPA receptors at synapses plays a role in mammalian learning and memory. Here we describe the design, synthesis, and evaluation of a photoreactive AMPA receptor antagonist that provides a means of "knocking out" AMPA receptors present on the surface of cells. The antagonist, 6-azido-7-nitro-1,4-dihydroquinoxaline-2,3-dione (ANQX), was designed by introducing a photoreactive azido group onto a quinoxalinedione inhibitor scaffold. Computational docking of ANQX to the AMPA receptor ligand-binding core predicted efficient binding to AMPA receptors. Glutamate-evoked currents were reversibly blocked at micromolar ANQX concentrations prior to photolysis and irreversibly blocked following photolysis. ANQX provides a means of directly evaluating the trafficking of native AMPA receptors with unparalleled spatiotemporal resolution.

Animals↗

PubMed2004

Computer-assisted design of selective imidazole inhibitors for cytochrome p450 enzymes.

A modified version of the DOCK program has been used to predict inhibitors for cytochrome P450cam and its L244A mutant. A library of azole compounds was designed in silico and screened for binding to wild-type P450cam. Lead compounds were synthesized and found to inhibit wild-type P450cam. To test our approach to designing ligands that discriminate between closely related sites, the azole library was DOCKed into both the active sites of wild-type P450cam and its L244A mutant. The L244A active site is predicted to be slightly larger than that of wild-type P450cam. Ligands predicted to have a high affinity for the mutant alone were synthesized and assayed with the recombinant enzymes. All of the compounds showed inhibition of the L244A enzyme (IC(50) = 6-40 microM), and the compounds that were predicted to be too large to bind to the wild-type showed poor inhibition (IC(50) > or = 1 mM). The binding mode was shown to be similar to that predicted by our modified version of DOCK by spectroscopic analysis. A discrepancy between the IC(50) values and spectroscopic K(s) values indicates that the spectroscopic binding constants do not accurately estimate inhibitory activity. This study, the first report of computer-assisted ligand (drug) design for P450 enzymes in which the coordination bond between imidazole and the heme is explicitly considered in structural modeling, opens a promising design avenue because azole compounds are widely used as P450 enzyme inhibitors and drugs.

Computer-Aided Design↗

PubMed2004

A molecular basis for the selectivity of thiadiazole urea inhibitors with stromelysin-1 and gelatinase-A from generalized born molecular dynamics simulations.

Matrix metalloproteinases (MMPs) represent a potentially important class of therapeutic targets for the treatment of diseases such as cancer. Selective inhibition of MMPs will be required given the high sequence identity across the family and the discovery that individual MMPs also regulate the natural angiogenesis inhibitor angiostatin. In this study, we have used computational methods to model the selectivity for six thiadiazole urea inhibitors with stromelysin-1 and gelatinase-A, two homologous MMPs that have been implicated in breast cancer. From continuum Generalized Born molecular dynamics (GB-MD) and MM-GBSA analysis, we estimated ligand free energies of binding using 200 snapshots obtained from a short 40 ps simulation of the relevant protein-ligand complex. The MM-GBSA free energies, computed from the continuum GB-MD trajectories, show strong correlation with the experimental affinities (r(2) = 0.74); prior studies have employed explicit water MD simulations. Including estimates for changes in solute entropy in the binding calculations slightly diminishes the overall correlation with experiment (r2 = 0.71). Notably, in every case, the simulation results correctly predict that a given ligand will bind selectively to stromelysin-1 over gelatinase-A which is gratifying given the high degree of structural homology between the two proteins. The increased selectivity for stromelysin-1 appears to be driven by (1) increased favorable van der Waals interactions, (2) increased favorable Coulombic interactions, and (3) decreased unfavorable total electrostatic energies (Coulombic plus desolvation).

Amino Acid Sequence↗

PubMed2004

Calculation of ligand-nucleic acid binding free energies with the generalized-born model in DOCK.

The calculation of ligand-nucleic acid binding free energies is investigated by including solvation effects computed with the generalized-Born model. Modifications of the solvation module in DOCK, including introduction of all-atom parameters and revision of coefficients in front of different terms, are shown to improve calculations involving nucleic acids. This computing scheme is capable of calculating binding energies, with reasonable accuracy, for a wide variety of DNA-ligand complexes, RNA-ligand complexes, and even for the formation of double-stranded DNA. This implementation of GB/SA is also shown to be capable of discriminating strong ligands from poor ligands for a series of RNA aptamers without sacrificing the high efficiency of the previous implementation. These results validate this approach to screening large databases against nucleic acid targets.

DNA↗

PubMed2004

Inhibition of Mycobacterium tuberculosis AhpD, an element of the peroxiredoxin defense against oxidative stress.

The resistance of Mycobacterium tuberculosis to isoniazid (INH) is largely linked to suppression of a catalase-peroxidase enzyme (KatG) that activates INH. In the absence of KatG, antioxidant protection is provided by enhanced expression of the peroxiredoxin AhpC, which is itself reduced by AhpD, a protein with low alkylhydroperoxidase activity of its own. Inhibition of AhpD might therefore impair the antioxidant protection afforded by AhpC and make KatG-negative strains more sensitive to oxidative stress. We report here that the 3(E),17-dioxime of testosterone is a potent competitive AhpD inhibitor, with a K(i) of 50 +/- 2 nM. The inhibitor is stereospecific, in that the 3(E) but not 3(Z) isomer is active. Computational studies provide support for a proposed AhpD substrate binding site. However, the inhibitor does not completely suppress the in vitro activity of AhpC/AhpD, because a low titer of AhpD suffices to maintain AhpC activity. This finding, and the low solubility of the inhibitor, explains its inability to suppress the growth of INH-resistant M. tuberculosis in infected mouse lungs.

Animals↗

PubMed2004

Distributions in protein conformation space: implications for structure prediction and entropy.

By considering how polymer structures are distributed in conformation space, we show that it is possible to quantify the difficulty of structural prediction and to provide a measure of progress for prediction calculations. The critical issue is the probability that a conformation is found within a specified distance of another conformer. We address this question by constructing a cumulative distribution function (CDF) for the average probability from observations about its limiting behavior at small displacements and numerical simulations of polyalanine chains. We can use the CDF to estimate the likelihood that a structure prediction is better than random chance. For example, the chance of randomly predicting the native backbone structure of a 150-amino-acid protein to low resolution, say within 6 A, is 10(-14). A high-resolution structural prediction, say to 2 A, is immensely more difficult (10(-57)). With additional assumptions, the CDF yields the conformational entropy of protein folding from native-state coordinate variance. Or, using values of the conformational entropy change on folding, we can estimate the native state's conformational span. For example, for a 150-mer protein, equilibrium alpha-carbon displacements in the native ensemble would be 0.3-0.5 A based on T Delta S of 1.42 kcal/(mol residue).

Algorithms↗

PubMed2003

Investigation of neuraminidase-substrate recognition using molecular dynamics and free energy calculations.

Development of the new generation of therapeutics against the influenza viral coat protein neuraminidase is a response to the continuing threat of influenza epidemics. A variety of structurally similar compounds have been reported that vary greatly in their ability to inhibit neuraminidase, a critical enzyme that cleaves sialic acid and promotes virion release. To determine how neuraminidase exhibits this wide range of affinities with structurally similar compounds, molecular dynamic simulations, coupled with free energy calculations, were used to determine the binding components of a series of neuraminidase inhibitors. Using four cocrystal structures of neuraminidase-inhibitor complexes, we examined the structural and energetic components of ligand potency and selectivity. An in-depth energetic analysis, including internal energy, entropy, and nonbonded interactions, reveals that potency of ligand binding is governed by nonpolar contacts. Electrostatic components generally oppose binding, although two of the best inhibitors use electrostatic interactions to orient the ligand. This investigation suggests that the enhanced selectivity and potency of the better ligands may arise from an improved positioning of their ligand atoms in the active site due to polar and hydrophobic functionalities. Simulations that included crystal water molecules in the active site indicate that the more potent ligands make less use of water-mediated interactions.

Alanine↗

PubMed2003

A selective irreversible inhibitor targeting a PDZ protein interaction domain.

Irreversible inhibitors of proteases have proven themselves useful tools for determining which proteases are active under given conditions in tissues or cells and for studying the functional role that a protease plays in physiological processes. The application of such techniques to the study of the activity and function of protein-protein interactions has been hindered by the lack of guiding principles for the mechanistic design of irreversible inhibitors targeting the "active site" of a protein interaction. We report herein the first example of a mechanism-based irreversible inhibitor of a protein interaction that has been specifically targeted to one member of the PDZ family of protein interaction domains: the second PDZ domain of the membrane-associated guanylate kinase MAGI3. This inhibitor was designed using rationally directed computational evaluation to take advantage of a conserved histidine in the PDZ domain by introducing an ionizable group that will be held in close proximity to that nucleophile during binding. The novel compound exhibits all of the characteristics of an irreversible inhibitor of the interaction of the tumor suppressor PTEN with MAGI3 in in vitro models. In cells, the inhibitor can be shown to release PTEN from sequestration by MAGI3 and consequently upregulate the PKB signaling pathway.

Drug Design↗

PubMed2003

Molecular recognition and docking algorithms.

Molecular docking is an invaluable tool in modern drug discovery. This review focuses on methodological developments relevant to the field of molecular docking. The forces important in molecular recognition are reviewed and followed by a discussion of how different scoring functions account for these forces. More recent applications of computational chemistry tools involve library design and database screening. Last, we summarize several critical methodological issues that must be addressed in future developments.

Algorithms↗

PubMed2003

Free energy calculations for theophylline binding to an RNA aptamer: Comparison of MM-PBSA and thermodynamic integration methods.

We have applied the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method (J. Srinivasan, T. E. Cheatham, P. Cieplak, P. A. Kollman, and D. A. Case, Journal of the American Chemical Society, 1998, Vol. 120, pp. 9401-9409) to study the interaction of an RNA aptamer with theophylline and its analogs. The MM-PBSA free energy analysis provides a reasonable absolute binding free energy for the RNA aptamer-theophylline complex formation. Energetic analysis reveals that the van der Waals interaction and the nonpolar contribution to solvation provide the basis for the favorable absolute free energy of complex. This trend is similar to other protein-ligand interactions studied previously. The MM-PBSA method also ranks the relative binding energies of five theophylline analogs approximately correctly, but not as well as the more conventional thermodynamic integration calculations, which were carried out to convert theophylline into its analogs. The comparison of MM-PBSA with TI suggests that the MM-PBSA method has some difficulties with the first-solvation-shell energetics.

Binding Sites↗

PubMed2003

Direct hydroxide attack is a plausible mechanism for amidase antibody 43C9.

Direct hydroxide attack on the scissile carbonyl of the substrate has been suggested as a likely mechanism for esterase antibodies elicited by phosphonate haptens, which mimic the transition states for the alkaline hydrolysis of esters.1 The unique amidase activity of esterase antibody 43C9 has been attributed to nucleophilic attack by an active-site histidine residue.2 Yet, the active site of 43C9 is strikingly similar to those of other esterase antibodies, particularly 17E8. We have carried out quantum mechanical calculations, molecular dynamics simulations, and free energy calculations to assess the mechanism involving direct hydroxide attack for 43C9. Results support this mechanism and suggest that the mechanism is plausible for other antiphosphonate antibodies that catalyze the hydrolysis of (p-nitro)phenyl esters.

Amidohydrolases↗