The additivity of substrate fragments in enzyme-ligand binding

Information decay in molecular docking screens against holo, apo, and modeled conformations of enzymes

Molecular docking uses the three-dimensional structure of a receptor to screen a small molecule database for potential ligands. The dependence of docking screens on the conformation of the binding site remains an open question. To evaluate the information loss that occurs as the active site conformation becomes less defined, a small molecule database was docked against the holo (ligand bound), apo, and homology modeled structures of 10 different enzyme binding sites. The holo and apo representations were crystallographic structures taken from the Protein Data Bank (PDB), and the homology-modeled structures were taken from the publicly available resource ModBase. The database docked was the MDL Drug Data Report (MDDR), a functionally annotated database of 95000 small molecules that contained at least 35 ligands for each of the 10 systems. In all sites, at least 99% of the molecules in the MDDR were treated as nonbinding decoys. For each system, the holo, apo, and modeled structures were used to screen the MDDR, and the ability of each structure to enrich the known ligands for that system over random selection was evaluated. The best overall enrichment was produced by the holo structure in seven systems, the apo structure in two systems, and the modeled structure in one system. These results suggest that the performance of the docking calculation is affected by the particular representation of the receptor used in the screen, and that the holo structure is the one most likely to yield the best discrimination between known ligands and decoy molecules, but important exceptions to this rule also emerge from this study. Although each of the holo, apo, and modeled conformations led to enrichment of known ligands in all systems, the enrichment did not always rise to a level judged to be sufficient to justify the effort of a docking screen. Using a 20-fold enrichment of known ligands over random selection as a rough guideline for what might be enough to justify a docking screen, the holo conformation of the enzyme met this criterion in eight of 10 sites, whereas the apo conformation met this criterion in only two sites and the modeled conformation in three.

The consequences of translational and rotational entropy lost by small molecules on binding to proteins

When a small molecule binds to a protein, it loses a significant amount of rigid body translational and rotational entropy. Estimates of the associated energy barrier vary widely in the literature yet accurate estimates are important in the interpretation of results from fragment-based drug discovery techniques. This paper describes an analysis that allows the estimation of the rigid body entropy barrier from the increase in binding affinities that results when two fragments of known affinity and known binding mode are joined together. The paper reviews the relatively rare number of examples where good quality data is available. From the analysis of this data, we estimate that the barrier to binding, due to the loss of rigid-body entropy, is 15-20 kJ/mol, i.e. around 3 orders of magnitude in affinity at 298 K. This large barrier explains why it is comparatively rare to observe multiple fragments binding to non-overlapping adjacent sites in enzymes. The barrier is also consistent with medicinal chemistry experience where small changes in the critical binding regions of ligands are often poorly tolerated by enzymes.

Tryptophan 80 and leucine 143 are critical for the hydride transfer step of thymidylate synthase by controlling active site access

Mutant forms of thymidylate synthase (TS) with substitutions at the conserved active site residue, Trp 80, are deficient in the hydride transfer step of the TS reaction. These mutants produce a beta-mercaptoethanol (beta-ME) adduct of the 2'-deoxyuridine-5'-monophosphate (dUMP) exocyclic methylene intermediate. Trp 80 has been proposed to assist hydride transfer by stabilizing a 5,6,7,8-tetrahydrofolate (THF) radical cation intermediate [Barrett, J. E., Lucero, C. M., and Schultz, P. G. (1999) J. Am. Chem. Soc. 121, 7965-7966.] formed after THF changes its binding from the cofactor pocket to a putative alternate site. To understand the molecular basis of hydride transfer deficiency in a mutant in which Trp 80 was changed to Gly, we determined the X-ray structures of this mutant Escherichia coli TS complexed with dUMP and the folate analogue 10-propargyl-5,8-dideazafolate (CB3717) and of the wild-type enzyme complexed with dUMP and THF. The mutant enzyme has a cavity in the active site continuous with bulk solvent. This cavity, sealed from bulk solvent in wild-type TS by Leu 143, would allow nucleophilic attack of beta-ME on the dUMP C5 exocyclic methylene. The structure of the wild-type enzyme/dUMP/THF complex shows that THF is bound in the cofactor binding pocket and is well positioned to transfer hydride to the dUMP exocyclic methylene. Together, these results suggest that THF does not reorient during hydride transfer and indicate that the role of Trp 80 may be to orient Leu 143 to shield the active site from bulk solvent and to optimally position the cofactor for hydride transfer.

High-throughput crystallography for lead discovery in drug design

Knowledge of the three-dimensional structures of protein targets now emerging from genomic data has the potential to accelerate drug discovery greatly. X-ray crystallography is the most widely used technique for protein structure determination, but technical challenges and time constraints have traditionally limited its use primarily to lead optimization. Here, we describe how significant advances in process automation and informatics have aided the development of high-throughput X-ray crystallography, and discuss the use of this technique for structure-based lead discovery.

Site-directed ligand discovery

We report a strategy (called "tethering") to discover low molecular weight ligands ( approximately 250 Da) that bind weakly to targeted sites on proteins through an intermediary disulfide tether. A native or engineered cysteine in a protein is allowed to react reversibly with a small library of disulfide-containing molecules ( approximately 1,200 compounds) at concentrations typically used in drug screening (10 to 200 microM). The cysteine-captured ligands, which are readily identified by MS, are among the most stable complexes, even though in the absence of the covalent tether the ligands may bind very weakly. This method was applied to generate a potent inhibitor for thymidylate synthase, an essential enzyme in pyrimidine metabolism with therapeutic applications in cancer and infectious diseases. The affinity of the untethered ligand (K(i) approximately 1 mM) was improved 3,000-fold by synthesis of a small set of analogs with the aid of crystallographic structures of the tethered complex. Such site-directed ligand discovery allows one to nucleate drug design from a spatially targeted lead fragment.

Catalytic cysteine of thymidylate synthase is activated upon substrate binding

The role of Ser 167 of Escherichia coli thymidylate synthase (TS) in catalysis has been characterized by kinetic and crystallographic studies. Position 167 variants including S167A, S167N, S167D, S167C, S167G, S167L, S167T, and S167V were generated by site-directed mutagenesis. Only S167A, S167G, S167T, and S167C complemented the growth of thymidine auxotrophs of E. coli in medium lacking thymidine. Steady-state kinetic analysis revealed that mutant enzymes exhibited k(cat) values 1.1-95-fold lower than that of the wild-type enzyme. Relative to wild-type TS, K(m) values of the mutant enzymes for 2'-deoxyuridylate (dUMP) were 5-90 times higher, while K(m) values for 5,10-methylenetetrahydrofolate (CH(2)H(4)folate) were 1.5-16-fold higher. The rate of dehalogenation of 5-bromo-2'-deoxyuridine 5'-monophosphate (BrdUMP), a reaction catalyzed by TS that does not require CH(2)H(4)folate as cosubstrate, by mutant TSs was analyzed and showed that only S167A and S167G catalyzed the dehalogenation reaction and values of k(cat)/K(m) for the mutant enzymes were decreased by 10- and 3000-fold, respectively. Analysis of pre-steady-state kinetics of ternary complex formation revealed that the productive binding of CH(2)H(4)folate is weaker to mutant TSs than to the wild-type enzyme. Chemical transformation constants (k(chem)) for the mutant enzymes were lower by 1.1-6.0-fold relative to the wild-type enzyme. S167A, S167T, and S167C crystallized in the I2(1)3 space group and scattered X-rays to either 1.7 A (S167A and S167T) or 2.6 A (S167C). The high-resolution data sets were refined to a R(crys) of 19.9%. In the crystals some cysteine residues were derivatized with 2-mercaptoethanol to form S,S-(2-hydroxyethyl)thiocysteine. The pattern of derivatization indicates that in the absence of bound substrate the catalytic cysteine is not more reactive than other cysteines. It is proposed that the catalytic cysteine is activated by substrate binding by a proton-transfer mechanism in which the phosphate group of the nucleotide neutralizes the charge of Arg 126', facilitating the transfer of a proton from the catalytic cysteine to a His 207-Asp 205 diad via a system of ordered water molecules.

Energetic contributions of four arginines to phosphate-binding in thymidylate synthase are more than additive and depend on optimization of "effective charge balance"

In thymidylate synthase, four conserved arginines provide two hydrogen bonds each to the oxygens of the phosphate group of the substrate, 2'-deoxyuridine-5'-monophosphate. Of these, R23, R178, and R179 are far removed from the site of methyl transfer and contribute to catalysis solely through binding and orientation of ligands. These arginines can be substituted by other residues, while still retaining more than 1% activity of the wild-type enzyme. We compared the kinetics and determined the crystal structures of dUMP complexes of three of the most active, uncharged single mutants of these arginines, R23I, R178T, and R179T, and of double mutants (R23I, R179T) and (R178T, R179T). The dramatically higher K(m) for R178T compared to the other two single mutants arises from the effects of R178 substitution on the orientation of dUMP; 10-15-fold increases in for R23I and R178T reflect the role of these residues in stabilizing the closed conformation of TS in ternary complexes. The free energy for productive dUMP binding, DeltaG(S), increases by at least 1 kcal/mol for each mutant, even when dUMP orientation and mobility in the crystal structure is the same as in wild-type enzyme. Thus, the four arginines do not contribute excess positive charge to the PO(4)(-2) binding site; rather, they ideally complement the charge and geometry of the phosphate moiety. More-than-additive increases in DeltaG(S) seen in the double mutants are consistent with quadratic increases in DeltaG(S) predicted for deviations from ideal electrostatic interactions and may also reflect cooperative binding of the arginines to the phosphate oxygens.

Inactivity of N229A thymidylate synthase due to water-mediated effects: isolating a late stage in methyl transfer

Mutation of thymidylate synthase N229(177) to alanine results in an essentially inactive enzyme, yet it leads to formation of a stable ternary complex. The kinetics of N229(177)A show that kcat for Escherichia coli is reduced by 200-fold while the Km for dUMP is increased 200-fold and the Km for folate increased by tenfold versus the wild-type enzyme. The crystal structures of N229(177)A in complex with dUMP and CB3717, and in complex with dUMP alone are determined at 2.4 A, and 2.5 A resolution. These structures identify the covalently bound ternary complex and show how N229(177)A traps an intermediate, and so becomes inactive in a later step of the reaction. Since the smaller alanine side-chain at N229(177)A does not directly sterically impair binding of ligands, the structures implicate, and place quantitative limits on the involvement of the structured water network in the active site of thymidylate synthase in both catalysis and in determining the binding affinity for dUMP (in contrast, the N229(177)V mutation in Lactobacillus casei has minimal effect on activity).

D221 in thymidylate synthase controls conformation change, and thereby opening of the imidazolidine

In thymidylate synthase (TS), the invariant residue Asp-221 provides the only side chain that hydrogen bonds to the pterin ring of the cofactor, 5,10-methylene-5,6,7,8-tetrahydrofolate. All mutants of D221 except cysteine abolish activity. We have determined the crystal structures of two ternary complexes of the Escherichia coli mutant D221N. In a complex with dUMP and the antifolate 10-propargyl-5,8-dideazafolate (CB3717), dUMP is covalently bound to the active site cysteine, as usual. CB3717, which has no imidazolidine ring, is also bound in the usual productive orientation, but is less ordered than in wild-type complexes. The side chain of Asn-221 still hydrogen bonds to N3 of the quinazoline ring of CB3717, which must be in the enol form. In contrast, the structure of D221N with 5-fluoro-dUMP and 5,10-methylene-5,6,7, 8-tetrahydrofolate shows the cofactor bound in two partially occupied, nonproductive binding sites. In both binding modes, the cofactor has a closed imidazolidine ring and adopts the solution conformation of the unbound cofactor. In one of the binding sites, the pterin ring is turned around such that Asn-221 hydrogen bonds to the unprotonated N1 instead of the protonated N3 of the cofactor. This orientation blocks the conformational change required for forming covalent ternary complexes. Taken together, the two crystal structures suggest that the hydrogen bond between the side chain of Asp-221 and N3 of the cofactor is most critical during the early steps of cofactor binding, where it enforces the correct orientation of the pterin ring. Proper orientation of the cofactor appears to be a prerequisite for opening the imidazolidine ring prior to formation of the covalent steady-state intermediate in catalysis.