Structure-based combinatorial library design: discovery of non-peptidic inhibitors of caspases 3 and 8.

Structure-based design of a combinatorial array was carried out in order to identify non-peptidic thiomethylketone inhibitors of caspases 3 and 8. Five compounds from the designed array were active against caspase 3, and two were active against caspase 8. A 2.5-A resolution co-crystal structure of caspase 3 and a thiomethylketone array member is reported. The structure-based design strategy has proved useful for identifying caspase inhibitors.

Potent and selective nonpeptide inhibitors of caspases 3 and 7 inhibit apoptosis and maintain cell functionality.

Caspases have been strongly implicated to play an essential role in apoptosis. A critical question regarding the role(s) of these proteases is whether selective inhibition of an effector caspase(s) will prevent cell death. We have identified potent and selective non-peptide inhibitors of the effector caspases 3 and 7. The inhibition of apoptosis and maintenance of cell functionality with a caspase 3/7-selective inhibitor is demonstrated for the first time, and suggests that targeting these two caspases alone is sufficient for blocking apoptosis. Furthermore, an x-ray co-crystal structure of the complex between recombinant human caspase 3 and an isatin sulfonamide inhibitor has been solved to 2.8-A resolution. In contrast to previously reported peptide-based caspase inhibitors, the isatin sulfonamides derive their selectivity for caspases 3 and 7 by interacting primarily with the S(2) subsite, and do not bind in the caspase primary aspartic acid binding pocket (S(1)). These inhibitors blocked apoptosis in murine bone marrow neutrophils and human chondrocytes. Furthermore, in camptothecin-induced chondrocyte apoptosis, cell functionality as measured by type II collagen promoter activity is maintained, an activity considered essential for cartilage homeostasis. These data suggest that inhibiting chondrocyte cell death with a caspase 3/7-selective inhibitor may provide a novel therapeutic approach for the prevention and treatment of osteoarthritis, or other disease states characterized by excessive apoptosis.

Nucleic acid base-pairing and N-methylacetamide self-association in chloroform: affinity and conformation.

A recently developed computational method, 'mining minima', is used to examine the hydrogen-bonding interactions of nucleic acid base-pairs and of the N-methylacetamide homodimer in chloroform. The mining minima algorithm aggressively samples molecular conformations, identifies the most important local minima, and computes their contributions to the overall free energy of the system. Here, the CHARMM 98 parameter set is used for the potential energy and the generalized Born/surface area solvent model is used to account for the influence of the solvent. Good agreement with experiment is obtained for the non-covalent binding affinities of a series of complexes. The computational approach used here is applicable to a range of molecular systems.

Thermodynamic linkage between the binding of protons and inhibitors to HIV-1 protease.

The aspartyl dyad of free HIV-1 protease has apparent pK(a)s of approximately 3 and approximately 6, but recent NMR studies indicate that the aspartyl dyad is fixed in the doubly protonated form over a wide pH range when cyclic urea inhibitors are bound, and in the monoprotonated form when the inhibitor KNI-272 is bound. We present computations and measurements related to these changes in protonation and to the thermodynamic linkage between protonation and inhibition. The Poisson-Boltzmann model of electrostatics is used to compute the apparent pK(a)s of the aspartyl dyad in the free enzyme and in complexes with four different inhibitors. The calculations are done with two parameter sets. One assigns epsilon = 4 to the solute interior and uses a detailed model of ionization; the other uses epsilon = 20 for the solute interior and a simplified representation of ionization. For the free enzyme, both parameter sets agree well with previously measured apparent pK(a)s of approximately 3 and approximately 6. However, the calculations with an internal dielectric constant of 4 reproduce the large pKa shifts upon binding of inhibitors, but the calculations with an internal dielectric constant of 20 do not. This observation has implications for the accurate calculation of pK(a)s in complex protein environments. Because binding of a cyclic urea inhibitor shifts the pK(a)s of the aspartyl dyad, changing the pH is expected to change its apparent binding affinity. However, we find experimentally that the affinity is independent of pH from 5.5 to 7.0. Possible explanations for this discrepancy are discussed.

A new class of models for computing receptor-ligand binding affinities.

Models for predicting the binding affinities of molecules in solution are either very detailed, making them computationally intensive and hard to test, or very simple, and thus less informative than one might wish. A new class of models that focus on the predominant states of the binding molecules promise to capture the essential physics of binding at modest computational cost.