Crystallographic studies of two alcohol dehydrogenase-bound analogues of thiazole-4-carboxamide adenine dinucleotide (TAD), the active anabolite of the antitumor agent tiazofurin

From ribavirin to NAD analogues and back to ribavirin in search for anticancer agents

Ribavirin, a broad-spectrum antiviral agent is used in the clinic alone or in combination with other antivirals and/or interferons. Numerous structural analogues of ribavirin have been developed, among them tiazofurin, which is inactive against viruses but is a potent anticancer drug. Tiazofurin was found to inhibit nicotinamide adenine dinucleotide (NAD)-dependent inosine monophosphate dehydrogenase (IMPDH) after metabolic conversion into tiazofurin adenine dinucleotide (TAD), which binds well but could not serve as IMPDH cofactor. TAD showed high selectivity against human IMPDH vs. other cellular dehydrogenases. Mycophenolic acid (MPA) was even more specific, binding at the cofactor-binding domain of IMPDH. Ribavirin adenine dinucleotide, however, did not show any significant inhibition at the enzymatic level. We synthesized numerous NAD analogues in which natural nicotinamide riboside was replaced by tiazofurin, MPA moiety, or benzamide riboside, and the adenosine moiety as well as the pyrophosphate linker were broadly modified. Some of these compounds were found to be low nanomolar inhibitors of the enzyme and sub-micromolar inhibitors of cancer cell line proliferation. The best were as potent as tyrosine kinase inhibitor gleevec heralded as a ‘magic bullet’ against chronic myelogenous leukemia. In recent years, ribavirin was rediscovered as a potential anticancer agent against number of tumors including leukemia. It was clearly established that its antitumor activity is related to the inhibition of an oncogene, the eukaryotic translation initiation factor (eIF4E).

Stereoelectronic effects dictate molecular conformation and biological function of heterocyclic amides

Heterocycles adjacent to amides can have important influences on molecular conformation due to stereoelectronic effects exerted by the heteroatom. This was shown for imidazole- and thiazole-amides by comparing low energy conformations (ab initio MP2 and DFT calculations), charge distribution, dipole moments, and known crystal structures which support a general principle. Switching a heteroatom from nitrogen to sulfur altered the amide conformation, producing different three-dimensional electrostatic surfaces. Differences were attributed to different dipole and orbital alignments and spectacularly translated into opposing agonist vs antagonist functions in modulating a G-protein coupled receptor for inflammatory protein complement C3a on human macrophages. Influences of the heteroatom were confirmed by locking the amide conformation using fused bicyclic rings. These findings show that stereoelectronic effects of heterocycles modulate molecular conformation and can impart strikingly different biological properties.

Interacting with the biomolecular solvent accessible surface via a haptic feedback device

Background

From the 1950s computer based renderings of molecules have been produced to aid researchers in their understanding of biomolecular structure and function. A major consideration for any molecular graphics software is the ability to visualise the three dimensional structure of the molecule. Traditionally, this was accomplished via stereoscopic pairs of images and later realised with three dimensional display technologies. Using a haptic feedback device in combination with molecular graphics has the potential to enhance three dimensional visualisation. Although haptic feedback devices have been used to feel the interaction forces during molecular docking they have not been used explicitly as an aid to visualisation.

Results

A haptic rendering application for biomolecular visualisation has been developed that allows the user to gain three-dimensional awareness of the shape of a biomolecule. By using a water molecule as the probe, modelled as an oxygen atom having hard-sphere interactions with the biomolecule, the process of exploration has the further benefit of being able to determine regions on the molecular surface that are accessible to the solvent. This gives insight into how awkward it is for a water molecule to gain access to or escape from channels and cavities, indicating possible entropic bottlenecks. In the case of liver alcohol dehydrogenase bound to the inhibitor SAD, it was found that there is a channel just wide enough for a single water molecule to pass through. Placing the probe coincident with crystallographic water molecules suggests that they are sometimes located within small pockets that provide a sterically stable environment irrespective of hydrogen bonding considerations.

Conclusion

By using the software, named HaptiMol ISAS (available from ), one can explore the accessible surface of biomolecules using a three-dimensional input device to gain insights into the shape and water accessibility of the biomolecular surface that cannot be so easily attained using conventional molecular graphics software.

Role of tryptophan 95 in substrate specificity and structural stability of Sulfolobus solfataricus alcohol dehydrogenase

A mutant of the thermostable NAD(+)-dependent (S)-stereospecific alcohol dehydrogenase from Sulfolobus solfataricus (SsADH) which has a single substitution, Trp95Leu, located at the substrate binding pocket, was fully characterized to ascertain the role of Trp95 in discriminating between chiral secondary alcohols suggested by the wild-type SsADH crystallographic structure. The Trp95Leu mutant displays no apparent activity with short-chain primary and secondary alcohols and poor activity with aromatic substrates and coenzyme. Moreover, the Trp --> Leu substitution affects the structural stability of the archaeal ADH, decreasing its thermal stability without relevant changes in secondary structure. The double mutant Trp95Leu/Asn249Tyr was also purified to assist in crystallographic analysis. This mutant exhibits higher activity but decreased affinity toward aliphatic alcohols, aldehydes as well as NAD(+) and NADH compared to the wild-type enzyme. The crystal structure of the Trp95Leu/Asn249Tyr mutant apo form, determined at 2.0 A resolution, reveals a large local rearrangement of the substrate site with dramatic consequences. The Leu95 side-chain conformation points away from the catalytic metal center and the widening of the substrate site is partially counteracted by a concomitant change of Trp117 side chain conformation. Structural changes at the active site are consistent with the reduced activity on substrates and decreased coenzyme binding.

Conformational changes and catalysis by alcohol dehydrogenase

As shown by X-ray crystallography, horse liver alcohol dehydrogenase undergoes a global conformational change upon binding of NAD(+) or NADH, involving a rotation of the catalytic domain relative to the coenzyme binding domain and the closing up of the active site to produce a catalytically efficient enzyme. The conformational change requires a complete coenzyme and is affected by various chemical or mutational substitutions that can increase the catalytic turnover by altering the kinetics of the isomerization and rate of dissociation of coenzymes. The binding of NAD(+) is kinetically limited by a unimolecular isomerization (corresponding to the conformational change) that is controlled by deprotonation of the catalytic zinc-water to produce a negatively-charged zinc-hydroxide, which can attract the positively-charged nicotinamide ring. The deprotonation is facilitated by His-51 acting through a hydrogen-bonded network to relay the proton to solvent. Binding of NADH also involves a conformational change, but the rate is very fast. After the enzyme binds NAD(+) and closes up, the substrate displaces the hydroxide bound to the catalytic zinc; this exchange may involve a double displacement reaction where the carboxylate group of a glutamate residue first displaces the hydroxide (inverting the tetrahedral coordination of the zinc), and then the exogenous ligand displaces the glutamate. The resulting enzyme-NAD(+)-alcoholate complex is poised for hydrogen transfer, and small conformational fluctuations may bring the reactants together so that the hydride ion is transferred by quantum mechanical tunneling. In the process, the nicotinamide ring may become puckered, as seen in structures of complexes of the enzyme with NADH. The conformational changes of alcohol dehydrogenase demonstrate the importance of protein dynamics in catalysis.

Nucleoside-5'-phosphoimidazolides: reagents for facile synthesis of dinucleoside pyrophosphates

A facile method is presented for preparation of dinucleoside pyrophosphate derivatives based on reaction of a nucleoside 5'-monophosphate with carbonyldiimidazole followed by treatment of the resulting nucleoside 5'-phosphoimidazolide with a nucleoside 5'-phosphate. This method is suitable for preparation of pyrophosphates analogous to NAD, FAD, and related natural pyrophosphates. The resulting compounds are useful for mechanistic studies of enzymes that use natural pyrophosphates as co-factors or substrates, and in development of inhibitors that have potential applications as therapeutic agents.

Progress in the understanding of drug-receptor interactions, part 2: experimental and theoretical electrostatic moments and interaction energies of an angiotensin II receptor antagonist (C30H30N6(O)3S)

A combined experimental and theoretical charge density study of an angiotensin II receptor antagonist (1) is presented focusing on electrostatic properties such as atomic charges, molecular electric moments up to the fourth rank and energies of the intermolecular interactions, to gain an insight into the physical nature of the drug-receptor interaction. Electrostatic properties were derived from both the experimental electron density (multipole refinement of X-ray data collected at T=17 K) and the ab initio wavefunction (single molecule and fully periodic calculations at the DFT level). The relevance of SO and SN intramolecular interactions on the activity of 1 is highlighted by using both the crystal and gas-phase geometries and their electrostatic nature is documented by means of QTAIM atomic charges. The derived electrostatic properties are consistent with a nearly spherical electron density distribution, characterised by an intermingling of electropositive and -negative zones rather than by a unique electrophilic region opposed to a nucleophilic area. This makes the first molecular moment scarcely significant and ill-determined, whereas the second moment is large, significant and highly reliable. A comparison between experimental and theoretical components of the third electric moment shows a few discrepancies, whereas the agreement for the fourth electric moment is excellent. The most favourable intermolecular bond is show to be an NHN hydrogen bond with an energy of about 50 kJ mol(-1). Key pharmacophoric features responsible for attractive electrostatic interactions include CHX hydrogen bonds. It is shown that methyl and methylene groups, known to be essential for the biological activity of the drug, provide a significant energetic contribution to the total binding energy. Dispersive interactions are important at the thiophene and at both the phenyl fragments. The experimental estimates of the electrostatic contribution to the intermolecular interaction energies of six molecular pairs, obtained by a new model proposed by Spackman, predict the correct relative electrostatic energies with no exceptions.

Amino acid composition of parallel helix-helix interfaces

Amino acids at helix-helix parallel interfaces influence arrangement of helices and interhelical angles. Parallel interfaces in 79 proteins were considered. Location of amino acids at the positions analogous to a and d in GCN4 leucine zipper nomenclature shows that certain combinations of amino acids characteristic for parallel packing occur more often than could be expected by chance. Repeating sequence combinations occur at a and d positions of parallel helix-helix interfaces with similar values of interhelical angles not only in homologous proteins but also within the same protein and in nonhomologous proteins. Within each group of observed combinations correlation exists between the size of amino acid and magnitude of the interhelical angle.

Supramolecular Synthesis of Selenazoles Using Selenourea in Water in the Presence of β-Cyclodextrin under Atmospheric Pressure

Selenazoles were synthesized from alpha-bromo ketones and selenourea in the presence of beta-cyclodextrin in water at 50 degrees C under atmospheric pressure.

Progress in the understanding of drug-receptor interactions, Part 1: experimental charge-density study of an angiotensin II receptor antagonist (C30H30N6O3S) at T = 17 K

An experimental study of the electron-density distribution rho(r) in an angiotensin II receptor antagonist 1 has been made on the basis of single-crystal X-ray diffraction data collected at a low temperature. The crystal structure of 1 consists of infinite ribbons in which molecules are connected by an N-H...N hydrogen bond and several interactions of the C-H...O, C-H...N, and C-H...S type. The molecular conformation, characterized by the syn orientation of a tetrazole and a pyrimidinone ring with respect to a phenyl spacer group, is stabilized by two short SO and SN intramolecular contacts between a substituted thiophene fragment and the other two heterocycles of 1. The electrostatic nature of these interactions is documented. Furthermore, the Laplacian of rho(r) in the plane defined by the sulfur, oxygen, and nitrogen atoms involved in these interactions shows their strongly directional character as the regions of charge concentration on the valence shell of the nitrogen and oxygen atoms directly face the regions of charge depletion on the valence shell of the sulfur atom. All the chemical bonds and the relevant intra- and intermolecular interactions of 1 have been quantitatively described by the topological analysis of rho(r). Simple relationships between the bond path lengths (R(b)) and the values of rho at the bond critical points (rho(bcp)) have been obtained for the 28 C-C bonds, the seven N-C bonds, and the four O-C bonds. For the first two classes of bonds the relationship is in the form of a straight line, whose parameters, for the C-C bonds, agree, within experimental uncertainty, with those previously derived in our laboratory from a 19 K X-ray diffraction study of crystals of a different compound. Maps of the molecular electrostatic potential phi(r) derived from the experimental charge density display features that are important for the drug-receptor recognition of 1.

Molecular dynamics simulations of NAD+-induced domain closure in horse liver alcohol dehydrogenase

Horse liver alcohol dehydrogenase is a homodimer, the protomer having a coenzyme-binding domain and a catalytic domain. Using all available x-ray structures and 50 ns of molecular dynamics simulations, we investigated the mechanism of NAD+-induced domain closure. When the well-known loop at the domain interface was modeled to its conformation in the closed structure, the NAD+-induced domain closure from the open structure could be simulated with remarkable accuracy. Native interactions in the closed structure between Arg369, Arg47, His51, Ala317, Phe319, and NAD+ were seen to form at different stages during domain closure. Removal of the Arg369 side-chain charge resulted in the loss of the tendency to close, verifying that specific interactions do help drive the domains closed. Further simulations and a careful analysis of x-ray structures suggest that the loop prevents domain closure in the absence of NAD+, and a cooperative mechanism operates between the subunits for domain closure. This cooperative mechanism explains the role of the loop as a block to closure because in the absence of NAD+ it would prevent the occurrence of an unliganded closed subunit when the other subunit closes on NAD+. Simulations that started with one subunit open and one closed supported this.

Synthesis, conformational analysis, and biological activity of new analogues of thiazole-4-carboxamide adenine dinucleotide (TAD) as IMP dehydrogenase inhibitors

Thiazole-4-carboxamide adenine dinucleotide (TAD) analogues T-2'-MeAD (1) and T-3'-MeAD (2) containing, respectively, a methyl group at the ribose 2'-C-, and 3'-C-position of the adenosine moiety, were prepared as potential selective human inosine monophosphate dehydrogenase (IMPDH) type II inhibitors. The synthesis of heterodinucleotides was carried out by CDI-catalyzed coupling reaction of unprotected 2'-C-methyl- or 3'-C-methyl-adenosine 5'-monophosphate with 2',3'-O-isopropylidene-tiazofurin 5'-monophosphate, and then deisopropylidenation. Biological evaluation of dinucleotides 1 and 2 as inhibitors of recombinant human IMPDH type I and type II resulted in a good activity. Inhibition of both isoenzymes by T-2'-MeAD and T-3'-MeAD was noncompetitive with respect to NAD substrate. Binding of T-3'-MeAD was comparable to that of parent compound TAD, while T-2'-MeAD proved to be a weaker inhibitor. However, no significant difference was found in inhibition of the IMPDH isoenzymes. T-2'-MeAD and T-3'-MeAD were found to inhibit the growth of K562 cells (IC(50) 30.7 and 65.0muM, respectively).

The oxidized subunit B8 from human complex I adopts a thioredoxin fold

Subunit B8 from ubiquinone oxidoreductase (complex I) (CI-B8) is one of several nuclear-encoded supernumerary subunits that are not present in bacterial complex I. Its solution structure shows a thioredoxin fold with highest similarities to the human thioredoxin mutant C73S and thioredoxin 2 from Anabeana sp. Interestingly, these proteins contain active sites in the same area, where the disulfide bond of oxidized CI-B8 is located. The redox potential of this disulfide bond is -251.6 mV, comparing well to that of disulfides in other thioredoxin-like proteins. Analysis of the structure reveals a surface area that is exclusively composed of highly conserved residues and thus most likely a subunit interaction site within complex I.

Statistical and visual morph movie analysis of crystallographic mutant selection bias in protein mutation resource data

Structural studies of the effects of non-silent mutations on protein conformational change are an important key in deciphering the language that relates protein amino acid primary structure to tertiary structure. Elsewhere, we presented the Protein Mutant Resource (PMR) database, a set of online tools that systematically identified groups of related mutant structures in the Protein DataBank (PDB), accurately inferred mutant classifications in the Gene Ontology using an innovative, statistically rigorous data-mining algorithm with more general applicability, and illustrated the relationship of these mutant structures via an intuitive user interface. Here, we perform a comprehensive statistical analysis of the effect of PMR mutations on protein tertiary structure. We find that, although the PMR does contain spectacular examples of conformational change, in general there is a counter-intuitive inverse relationship between conformational change (measured as C-alpha displacement or RMS of the core structure) and the number of mutations in a structure. That is, point mutations by structural biologists present in the PDB contrast naturally evolved mutations. We compare the frequency of mutations in the PMR/PDB datasets against the accepted PAM250 natural amino acid mutation frequency to confirm these observations. We generated morph movies from PMR structure pairs using technology previously developed for the Macromolecular Motions Database (http://molmovdb.org), allowing bioinformaticians, geneticists, protein engineers, and rational drug designers to analyze visually the mechanisms of protein conformational change and distinguish between conformational change due to motions (e.g., ligand binding) and mutations. The PMR morph movies and statistics can be freely viewed from the PMR website, http://pmr.sdsc.edu.

Identification of Specific Interactions that Drive Ligand-induced Closure in Five Enzymes with Classic Domain Movements

In order to better understand ligand-induced closure in domain enzymes, open unliganded X-ray structures and closed liganded X-ray structures have been studied in five enzymes: adenylate kinase, aspartate aminotransferase, citrate synthase, liver alcohol dehydrogenase, and the catalytic subunit of cAMP-dependent protein kinase. A sequential model of ligand binding and domain closure was used to test the hypothesis that the ligand actively drives closure from an open conformation. The analysis supports the assumption that each enzyme has a dedicated binding domain to which the ligand binds first and a closing domain. In every case, a small number of residues are identified to interact with the ligand to initiate and drive domain closure. In all cases except adenylate kinase, the backbone of residues located in an interdomain-bending region (hinge site) is identified to interact with the ligand to aid in driving closure. In adenylate kinase, the side-chain of a residue located directly adjacent to a bending region drives closure. It is thought that by binding near a hinge site the ligand is able to get within interaction range of residues when the enzyme is in the open conformation. Interdomain bending regions not involved in inducing closure are involved in control, helping to determine the location of the hinge axis. Similarities have been discovered between aspartate aminotransferase and citrate synthase that only come to light in the context of their dynamical behaviour in response to binding their substrate. Similarity also exists between liver alcohol dehydrogenase and cAMP-dependent protein kinase whereby groups on NAD and ATP, respectively, mimic the backbone of a single amino acid residue in a process where a three residue segment located at the terminus of a beta-sheet, moves to form hydrogen bonds with the mimic that resemble those found in a parallel beta-sheet. This interaction helps to drive domain closure in a process that has analogy to protein folding.

Crystal structure of human type II inosine monophosphate dehydrogenase: implications for ligand binding and drug design

Inosine monophosphate dehydrogenase (IMPDH) controls a key metabolic step in the regulation of cell growth and differentiation. This step is the NAD-dependent oxidation of inosine 5' monophosphate (IMP) to xanthosine 5' monophosphate, the rate-limiting step in the synthesis of the guanine nucleotides. Two isoforms of IMPDH have been identified, one of which (type II) is significantly up- regulated in neoplastic and differentiating cells. As such, it has been identified as a major target in antitumor and immunosuppressive drug design. We present here the 2.9-A structure of a ternary complex of the human type II isoform of IMPDH. The complex contains the substrate analogue 6-chloropurine riboside 5'-monophosphate (6-Cl-IMP) and the NAD analogue selenazole-4-carboxamide adenine dinucleotide, the selenium derivative of the active metabolite of the antitumor drug tiazofurin. The enzyme forms a homotetramer, with the dinucleotide binding at the monomer-monomer interface. The 6 chloro-substituted purine base is dehalogenated, forming a covalent adduct at C6 with Cys-331. The dinucleotide selenazole base is stacked against the 6-Cl-IMP purine ring in an orientation consistent with the B-side stereochemistry of hydride transfer seen with NAD. The adenosine end of the ligand interacts with residues not conserved between the type I and type II isoforms, suggesting strategies for the design of isoform-specific agents.

Novel nicotinamide adenine dinucleotide analogues as potential anticancer agents: quest for specific inhibition of inosine monophosphate dehydrogenase

Synthetic nicotinamide adenine dinucleotide (NAD) analogues containing 5-beta-D-ribofuranosylnicotinamide (C-NAD), 6-beta-D-ribofuranosylpicolinamide (C-PAD), 3-beta-D-ribofuranosylbenzamide (BAD), and 2-beta-D-ribofuranosylthiazole-4-carboxamide (TAD) in place of the nicotinamide riboside moiety are described and evaluated as potential inhibitors of inosine monophosphate dehydrogenase (IMPDH). TAD and BAD showed potent inhibitory activity against the enzyme in the form of pyrophosphates, as well as metabolically stable methylene- and difluoromethylenebis(phosphonate)s. Fluorination at the C2' (ribo and arabino configuration) and C3' (ribo) of the adenosine moiety of TAD afforded analogues highly potent against IMPDH, but weakly active against alcohol dehydrogenase. With the exception of the methylenebis(phosphonate) analogue of TAD compounds containing a methylene bridge were poor inhibitors of growth of K562 cells. On the other hand, NAD analogues containing difluoromethylene linkage were highly effective in inhibition of K562 cell growth, as well as potent inducers of K562 cell differentiation. Such compounds, therefore, may be of potential therapeutic interest.

The protein conformation of Cd-substituted horse liver alcohol dehydrogenase and its metal-site coordination geometry in binary and ternary inhibitor complexes

The coordination geometry of the metal at the active site in Cd-substituted horse liver alcohol dehydrogenase (LADH) has been investigated for the binary complexes of LADH with imidazole, isobutyramide, decanoic acid and Cl-, and for the ternary complexes of LADH with NADH and imidazole, NADH and isobutyramide, NAD+ and decanoic acid and NAD+ and Cl-, by using the method of perturbed angular correlation of gamma-rays (PAC). The spectral results are consistent with a flexible structure around the metal for the binary complexes with inhibitors. For ternary complexes, however, a rigid structure is observed. An exception is the ternary complex between LADH, NADH and imidazole, in which the metal site is still flexible. Comparing with available structures determined by X-ray crystallography, we found a correlation between open structures and flexible metal sites, and between closed structures and rigid metal sites. This indicates that the PAC technique can be applied to distinguish the two conformations in solution. The spectral parameters, omega(o) and eta, of the experiments, except for the complexes with imidazole, fall into two groups: one with low omega(o) and one with high omega(o) (eta is relatively constant in all experiments). In this work it is clarified that the low omega(o) values are connected with the presence of a negatively charged solvent ligand. Using an angular-overlap approach to interpret the results, the low omega(o) values are found to be compatible with a coordination geometry where the S-Cd-S (Cys174 and Cys46 coordinate to the metal) angle is about 110 degrees as suggested in [Hemmingsen, L., Bauer, R., Danielsen, E., Bjerrum. M. J., Zeppezauer, M., Adolph, H. W., Formicka, G. & Cedergren-Zeppezauer, E. (1995) Biochemistry 34, 7145-7153], whereas high omega(o) values are compatible with an S-Cd-S angle of 130 degrees. The presence of a negatively charged metal ligand, therefore, might trigger the movement of the sulfur of Cys174. As it is believed that alcohols coordinate to the metal as alcoholate ions this could be important for catalysis.

Crystal structure of the catalytic domain of Pseudomonas exotoxin A complexed with a nicotinamide adenine dinucleotide analog: implications for the activation process and for ADP ribosylation

The catalytic, or third domain of Pseudomonas exotoxin A (PEIII) catalyzes the transfer of ADP ribose from nicotinamide adenine dinucleotide (NAD) to elongation factor-2 in eukaryotic cells, inhibiting protein synthesis. We have determined the structure of PEIII crystallized in the presence of NAD to define the site of binding and mechanism of activation. However, NAD undergoes a slow hydrolysis and the crystal structure revealed only the hydrolysis products, AMP and nicotinamide, bound to the enzyme. To better define the site of NAD binding, we have now crystallized PEIII in the presence of a less hydrolyzable NAD analog, beta-methylene-thiazole-4-carboxamide adenine dinucleotide (beta-TAD), and refined the complex structure at 2.3 angstroms resolution. There are two independent molecules of PEIII in the crystal, and the conformations of beta-TAD show some differences in the two binding sites. The beta-TAD attached to molecule 2 appears to have been hydrolyzed between the pyrophosphate and the nicotinamide ribose. However molecule 1 binds to an intact beta-TAD and has no crystal packing contacts in the vicinity of the binding site, so that the observed conformation and interaction with the PEIII most likely resembles that of NAD bound to PEIII in solution. We have compared this complex with the catalytic domains of diphtheria toxin, heat labile enterotoxin, and pertussis toxin, all three of which it closely resembles.

Crystallographic studies of isosteric NAD analogues bound to alcohol dehydrogenase: specificity and substrate binding in two ternary complexes

CNAD (5-beta-D-ribofuranosylnicotinamide adenine dinucleotide) is an isosteric C-glycosidic analogue of NAD(H) containing a neutral pyridine ring. CPAD (5-beta-D-ribofuranosylpicolinamide adenine dinucleotide) is a closely related pyridine-containing analogue with the pyridine nitrogen on the opposite side of the ring. CNAD is a potent and specific inhibitor of horse liver alcohol dehydrogenase (LADH), binding with a dissociation constant in the nanomolar range. CPAD binds LADH with an affinity comparable to that of NAD. Crystal structures of CNAD and CPAD bound to LADH are presented at 2.4 and 2.7 A, respectively. The two complexes are isomorphous, crystallizing in the triclinic system with cell dimensions different from those seen in previous ternary LADH complexes. Structures were solved using the molecular replacement method and refined to crystallographic R values of 18% (CNAD) and 17% (CPAD). Both inhibitors bind to the "closed" form of LADH in the normal cofactor-binding cleft. The conformation of LADH-bound CPAD closely mimics that of LADH-bound NAD(H). The data suggest that alcohol substrate binds directly to the catalytic zinc atom. In the CNAD complex, the pyridine nitrogen replaces alcohol as the fourth coordination ligand to the active site zinc atom, while all other polar interactions remain the same as those of bound NAD(H). The zinc-nitrogen ligand explains the high affinity of CNAD for LADH.