The additivity of substrate fragments in enzyme-ligand binding

The genetic landscape of a metabolic interaction

While much prior work has explored the constraints on protein sequence and evolution induced by physical protein-protein interactions, the sequence-level constraints emerging from non-binding functional interactions in metabolism remain unclear. To quantify how variation in the activity of one enzyme constrains the biochemical parameters and sequence of another, we focus on dihydrofolate reductase (DHFR) and thymidylate synthase (TYMS), a pair of enzymes catalyzing consecutive reactions in folate metabolism. We use deep mutational scanning to quantify the growth rate effect of 2696 DHFR single mutations in 3 TYMS backgrounds under conditions selected to emphasize biochemical epistasis. Our data are well-described by a relatively simple enzyme velocity to growth rate model that quantifies how metabolic context tunes enzyme mutational tolerance. Together our results reveal the structural distribution of epistasis in a metabolic enzyme and establish a foundation for the design of multi-enzyme systems.

The Genetic Landscape of a Metabolic Interaction

Enzyme abundance, catalytic activity, and ultimately sequence are all shaped by the need of growing cells to maintain metabolic flux while minimizing accumulation of deleterious intermediates. While much prior work has explored the constraints on protein sequence and evolution induced by physical protein-protein interactions, the sequence-level constraints emerging from non-binding functional interactions in metabolism remain unclear. To quantify how variation in the activity of one enzyme constrains the biochemical parameters and sequence of another, we focused on dihydrofolate reductase (DHFR) and thymidylate synthase (TYMS), a pair of enzymes catalyzing consecutive reactions in folate metabolism. We used deep mutational scanning to quantify the growth rate effect of 2,696 DHFR single mutations in 3 TYMS backgrounds under conditions selected to emphasize biochemical epistasis. Our data are well-described by a relatively simple enzyme velocity to growth rate model that quantifies how metabolic context tunes enzyme mutational tolerance. Together our results reveal the structural distribution of epistasis in a metabolic enzyme and establish a foundation for the design of multi-enzyme systems.

A defined structural unit enables de novo design of small-molecule-binding proteins

The de novo design of proteins that bind highly functionalized small molecules represents a great challenge. To enable computational design of binders, we developed a unit of protein structure-a van der Mer (vdM)-that maps the backbone of each amino acid to statistically preferred positions of interacting chemical groups. Using vdMs, we designed six de novo proteins to bind the drug apixaban; two bound with low and submicromolar affinity. X-ray crystallography and mutagenesis confirmed a structure with a precisely designed cavity that forms favorable interactions in the drug-protein complex. vdMs may enable design of functional proteins for applications in sensing, medicine, and catalysis.

Inter-Active Site Communication Mediated by the Dimer Interface β-Sheet in the Half-the-Sites Enzyme, Thymidylate Synthase

Thymidylate synthase (TS) is a dimeric enzyme conserved in all life forms that exhibits the allosteric feature of half-the-sites activity. Neither the reason for nor the mechanism of this phenomenon is understood. We used a combined nuclear magnetic resonance (NMR) and molecular dynamics approach to study a stable intermediate preceding hydride transfer, which is the rate-limiting and half-the-sites step. In NMR titrations with ligands leading to this intermediate, we measured chemical shifts of the apoenzyme (lig₀), the saturated holoenzyme (lig₂), and the typically elusive singly bound (lig₁) states. Approximately 40 amides showed quartet patterns providing direct NMR evidence of coupling between the active site and probes >30 Å away in the distal subunit. Quartet peak patterns have symmetrical character, indicating reciprocity in communicating the first and second binding events to the distal protomer. Quartets include key catalytic residues and map to the dimer interface β-sheet, which also represents the shortest path between the two active sites. Simulations corroborate the coupling observed in solution in that there is excellent overlap between quartet residues and main-chain atoms having intersubunit cross-correlated motions. Simulations identify five hot spot residues, three of which lie at the kink in the unique β-bulge abutting the active sites on either end of the sheet. Interstrand cross-correlated motions become more organized and pronounced as the enzyme progresses from lig₀ to lig₁ and ultimately lig₂. Coupling in the apparently symmetrical complex has implications for half-the-sites reactivity and potentially resolves the paradox of inequivalent TS active sites despite the vast majority of X-ray structures appearing to be symmetrical.

Structural Comparison of Enterococcus faecalis and Human Thymidylate Synthase Complexes with the Substrate dUMP and Its Analogue FdUMP Provides Hints about Enzyme Conformational Variabilities

Thymidylate synthase (TS) is an enzyme of paramount importance as it provides the only de novo source of deoxy-thymidine monophosphate (dTMP). dTMP, essential for DNA synthesis, is produced by the TS-catalyzed reductive methylation of 2′-deoxyuridine-5′-monophosphate (dUMP) using N⁵,N¹⁰-methylenetetrahydrofolate (mTHF) as a cofactor. TS is ubiquitous and a validated drug target. TS enzymes from different organisms differ in sequence and structure, but are all obligate homodimers. The structural and mechanistic differences between the human and bacterial enzymes are exploitable to obtain selective inhibitors of bacterial TSs that can enrich the currently available therapeutic tools against bacterial infections. Enterococcus faecalis is a pathogen fully dependent on TS for dTMP synthesis. In this study, we present four new crystal structures of Enterococcus faecalis and human TSs in complex with either the substrate dUMP or the inhibitor FdUMP. The results provide new clues about the half-site reactivity of Enterococcus faecalis TS and the mechanisms underlying the conformational changes occurring in the two enzymes. We also identify relevant differences in cofactor and inhibitor binding between Enterococcus faecalis and human TS that can guide the design of selective inhibitors against bacterial TSs.

Experimental and Computational Studies Delineate the Role of Asparagine 177 in Hydride Transfer for E. coli Thymidylate Synthase

Thymidylate synthase (TSase), an enzyme responsible for the de novo biosynthesis of 2'-deoxythymidine 5'-monophosphate (thymidylate, dTMP) necessary for DNA synthesis, has been a drug target for decades. TSase is a highly conserved enzyme across species ranging from very primitive organisms to mammals. Among the many conserved active site residues, an asparagine (N177, using Escherichia coli residues numbering) appears to make direct hydrogen bonds with both the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate and its pyrimidine ring's N3. Recent studies have reassessed the TSase catalytic mechanism, focusing on the degree of negative charge accumulation at the O4 carbonyl of the substrate during two critical H-transfers - a proton abstraction and a hydride transfer. To obtain insights into the role of this conserved N177 on the hydride transfer, we examined its aspartic acid (D) and serine (S) mutants - each of which is expected to alter hydrogen bonding and charge stabilization around the C4=O4 carbonyl of the 2'-deoxyuridine 5'-monophosphate (uridylate, dUMP) substrate. Steady-state kinetics, substrate binding order studies and temperature-dependency analysis of intrinsic KIEs for the hydride transfer step of the TSase catalytic cycle suggest the active site of N177D is not precisely organized for that step. A smaller disruption was observed for N177S, which could be rationalized by partial compensation by water molecules and rearrangement of other residues toward preparation of the system for the hydride transfer under study. These experimental findings are qualitatively mirrored by QM/MM computational simulations, thereby shedding light on the sequence and synchronicity of steps in the TSase-catalyzed reaction. This information could potentially inform the design of mechanism-based drugs targeting this enzyme.

Thermodynamic and NMR Assessment of Ligand Cooperativity and Intersubunit Communication in Symmetric Dimers: Application to Thymidylate Synthase

Thymidylate synthase (TS) is a homodimeric enzyme with evidence for negative regulation of one protomer while the other protomer acts on substrate, so called half-the-sites reactivity. The mechanisms by which multisubunit allosteric proteins communicate between protomers is not well understood, and the simplicity of dimeric systems has advantages for observing conformational and dynamic processes that functionally connect distance-separated active sites. This review considers progress in overcoming the inherent challenges of accurate thermodynamic and atomic-resolution characterization of interprotomer communication mechanisms in symmetric protein dimers, with TS used as an example. Isothermal titration calorimetry (ITC) is used to measure ligand binding cooperativity, even in cases where the two binding enthalpies are similar, and NMR spectroscopy is used to detect site-specific changes occurring in the two protomers. The NMR approach makes use of mixed-labeled dimers, enabling protomer-specific detection of signals in the singly ligated state. The rich informational content of the NMR signals from the singly ligated state, relative to the apo and saturated states, requires new considerations that do not arise in simple cases of 1:1 protein-ligand interactions.

A Single Mutation Traps a Half-Sites Reactive Enzyme in Midstream, Explaining Asymmetry in Hydride Transfer

In Escherichia coli thymidylate synthase (EcTS), rate-determining hydride transfer from the cofactor 5,10-methylene-5,6,7,8-tetrahydrofolate to the intermediate 5-methylene-2'-deoxyuridine 5'-monophosphate occurs by hydrogen tunneling, requiring precise alignment of reactants and a closed binding cavity, sealed by the C-terminal carboxyl group. Mutations that destabilize the closed conformation of the binding cavity allow small molecules such as β-mercaptoethanol (β-ME) to enter the active site and compete with hydride for addition to the 5-methylene group of the intermediate. The C-terminal deletion mutant of EcTS produced the β-ME adduct in proportions that varied dramatically with cofactor concentration, from 50% at low cofactor concentrations to 0% at saturating cofactor conditions, suggesting communication between active sites. We report the 2.4 Å X-ray structure of the C-terminal deletion mutant of E. coli TS in complex with a substrate and a cofactor analogue, CB3717. The structure is asymmetric, with reactants aligned in a manner consistent with hydride transfer in only one active site. In the second site, CB3717 has shifted to a site where the normal cofactor would be unlikely to form 5-methylene-2'-deoxyuridine 5'-monophosphate, consistent with no formation of the β-ME adduct. The structure shows how the binding of the cofactor at one site triggers hydride transfer and borrows needed stabilization from substrate binding at the second site. It indicates pathways through the dimer interface that contribute to allostery relevant to half-sites reactivity.

New insights into the molecular characteristics behind the function of Renilla luciferase

Renilla Luciferase (RLuc) is a blue light emitter protein which can be applied as a valuable tool in medical diagnosis. But due to lack of the crystal structure of RLuc-ligand complex, the functional motions and catalytic mechanism of this enzyme remain largely unknown. In the present study, the active site properties and the ligand-receptor interactions of the native RLuc and its red-shifted light emitting variant (Super RLuc 8) were investigated using molecular docking approach, molecular dynamics (MD) analysis, and MM-PBSA method. The detailed analysis of the main clusters led to identifying a lid-like structure and its functional motions. Furthermore, an induced-fit mechanism is proposed where ligand-binding induces conformational changes of the active site. Our findings give an insight into the deeper understanding of RLuc conformational changes during binding steps and ligand-receptor pattern. Moreover, our work broaden our understanding of how active site geometry is adjusted to support the catalytic activity and red-shifted light emission in Super RLuc 8.

Crystal structures of nematode (parasitic T. spiralis and free living C. elegans), compared to mammalian, thymidylate synthases (TS). Molecular docking and molecular dynamics simulations in search for nematode-specific inhibitors of TS

Three crystal structures are presented of nematode thymidylate synthases (TS), including Caenorhabditis elegans (Ce) enzyme without ligands and its ternary complex with dUMP and Raltitrexed, and binary complex of Trichinella spiralis (Ts) enzyme with dUMP. In search of differences potentially relevant for the development of species-specific inhibitors of the nematode enzyme, a comparison was made of the present Ce and Ts enzyme structures, as well as binary complex of Ce enzyme with dUMP, with the corresponding mammalian (human, mouse and rat) enzyme crystal structures. To complement the comparison, tCONCOORD computations were performed to evaluate dynamic behaviors of mammalian and nematode TS structures. Finally, comparative molecular docking combined with molecular dynamics and free energy of binding calculations were carried out to search for ligands showing selective affinity to T. spiralis TS. Despite an overall strong similarity in structure and dynamics of nematode vs mammalian TSs, a pool of ligands demonstrating predictively a strong and selective binding to TsTS has been delimited. These compounds, the E63 family, locate in the dimerization interface of TsTS where they exert species-specific interactions with certain non-conserved residues, including hydrogen bonds with Thr174 and hydrophobic contacts with Phe192, Cys191 and Tyr152. The E63 family of ligands opens the possibility of future development of selective inhibitors of TsTS and effective agents against trichinellosis.

Classification of ligand molecules in PDB with graph match-based structural superposition

The fast heuristic graph match algorithm for small molecules, COMPLIG, was improved by adding a structural superposition process to verify the atom-atom matching. The modified method was used to classify the small molecule ligands in the Protein Data Bank (PDB) by their three-dimensional structures, and 16,660 types of ligands in the PDB were classified into 7561 clusters. In contrast, a classification by a previous method (without structure superposition) generated 3371 clusters from the same ligand set. The characteristic feature in the current classification system is the increased number of singleton clusters, which contained only one ligand molecule in a cluster. Inspections of the singletons in the current classification system but not in the previous one implied that the major factors for the isolation were differences in chirality, cyclic conformations, separation of substructures, and bond length. Comparisons between current and previous classification systems revealed that the superposition-based classification was effective in clustering functionally related ligands, such as drugs targeted to specific biological processes, owing to the strictness of the atom-atom matching.

Structural basis of the substrate preference towards CMP for a thymidylate synthase MilA involved in mildiomycin biosynthesis

Modified pyrimidine monophosphates such as methyl dCMP (mdCMP), hydroxymethyl dUMP (hmdUMP) and hmdCMP in some phages are synthesized by a large group of enzymes termed as thymidylate synthases (TS). Thymidylate is a nucleotide required for DNA synthesis and thus TS is an important drug target. In the biosynthetic pathway of the nucleoside fungicide mildiomycin isolated from Streptomyces rimofaciens ZJU5119, a cytidylate (CMP) hydroxymethylase, MilA, catalyzes the conversion of CMP into 5′-hydroxymethyl CMP (hmCMP) with an efficiency (k_cat/K_M) of 5-fold faster than for deoxycytidylate (dCMP). MilA is thus the first enzyme of the TS superfamily preferring CMP to dCMP. Here, we determined the crystal structures of MilA and its complexes with various substrates including CMP, dCMP and hmCMP. Comparing these structures to those of dCMP hydroxymethylase (CH) from T4 phage and TS from Escherichia coli revealed that two residues in the active site of CH and TS, a serine and an arginine, are respectively replaced by an alanine and a lysine, Ala176 and Lys133, in MilA. Mutation of A176S/K133R of MilA resulted in a reversal of substrate preference from CMP to dCMP. This is the first study reporting the evolution of the conserved TS in substrate selection from DNA metabolism to secondary nucleoside biosynthesis.

When Does Chemical Elaboration Induce a Ligand To Change Its Binding Mode?

Traditional hit-to-lead optimization assumes that upon elaboration of chemical structure, the ligand retains its binding mode relative to the receptor. Here, we build a large-scale collection of related ligand pairs solved in complex with the same protein partner: we find that for 41 of 297 pairs (14%), the binding mode changes upon elaboration of the smaller ligand. While certain ligand physiochemical properties predispose changes in binding mode, particularly those properties that define fragments, simple structure-based modeling proves far more effective for identifying substitutions that alter the binding mode. Some ligand pairs change binding mode because the added substituent would irreconcilably conflict with the receptor in the original pose, whereas others change because the added substituent enables new, stronger interactions that are available only in a different pose. Scaffolds that can engage their target using alternate poses may enable productive structure-based optimization along multiple divergent pathways.

Widespread Perturbation of Function, Structure, and Dynamics by a Conservative Single-Atom Substitution in Thymidylate Synthase

Thymidylate synthase (TSase) is responsible for synthesizing the sole de novo source of dTMP in all organisms. TSase is a drug target, and as such, it has been well studied in terms of both structure and reaction mechanism. Cysteine 146 in Escherichia coli TSase is universally conserved because it serves as the nucleophile in the enzyme mechanism. Here we use the C146S mutation to probe the role of the sulfur atom in early events in the catalytic cycle beyond serving as the nucleophile. Surprisingly, the single-atom substitution severely decreases substrate binding affinity, and the unfavorable ΔΔG°_bind is comprised of roughly equal enthalpic and entropic components at 25 °C. Chemical shifts in the free and dUMP-bound states show the mutation causes perturbations throughout TSase, including regions important for complex stability, in agreement with a less favorable enthalpy change. We measured the nuclear magnetic resonance methyl symmetry axis order parameter (S²_axis), a proxy for conformational entropy, for TSase at all vertices of the dUMP binding/C146S mutation thermodynamic cycle and found that the calculated TΔΔS°_conf is similar in sign and magnitude to the calorimetric TΔΔS°. Further, we ascribed minor resonances in wild-type-dUMP spectra to a state with a covalent bond between Sγ of C146 and C6 of dUMP and find S²_axis values are unaffected by covalent bond formation, indicating this reaction step is neutral with respect to ΔS°_conf. Lastly, the C146S mutation allowed us to measure cofactor analog binding by isothermal titration calorimetry without the confounding heat signature of covalent bond formation. Raltitrexed binds free and singly bound TSase with similar affinities, yet the two binding events have different enthalpy changes, providing further evidence of communication between the two active sites.

Substrate Fragmentation for the Design of M. tuberculosis CYP121 Inhibitors

The cyclo-dipeptide substrates of the essential M. tuberculosis (Mtb) enzyme CYP121 were deconstructed into their component fragments and screened against the enzyme. A number of hits were identified, one of which exhibited an unexpected inhibitor-like binding mode. The inhibitory pharmacophore was elucidated, and fragment binding affinity was rapidly improved by synthetic elaboration guided by the structures of CYP121 substrates. The resulting inhibitors have low micromolar affinity, good predicted physicochemical properties and selectivity for CYP121 over other Mtb P450s. Spectroscopic characterisation of the inhibitors' binding mode provides insight into the effect of weak nitrogen-donor ligands on the P450 heme, an improved understanding of factors governing CYP121-ligand recognition and speculation into the biological role of the enzyme for Mtb.

Non-native ligands define the active site of Pennisetum glaucum (L.) R. Br dehydroascorbate reductase

Dehydroascorbate reductase (DHAR), a member of the glutathione-S-transferase (GST) family, reduces dehydroascorbate (DHA) to ascorbate (AsA; Vitamin-C) in a glutathione (GSH)-dependent manner and in doing so, replenishes the critical AsA pool of the cell. To understand the enzyme mechanism in detail, we determined the crystal structure of a plant DHAR from Pennisetum glaucum (PgDHAR) using Iodide-Single Anomalous Dispersion (SAD) and Molecular replacement methods, in two different space groups. Here, we show PgDHAR in complex with two non-native ligands, viz. an acetate bound at the G-site, which resembles the γ-carboxyl moiety of GSH, and a glycerol at the H-site, which shares the backbone of AsA. We also show that, in the absence of bound native substrates, these non-native ligands help define the critical 'hook points' in the DHAR enzyme active site. Further, our data suggest that these non-native ligands can act as the logical bootstrapping points for iterative design of inhibitors/analogs for DHARs.

Fragment-Based Approaches to the Development of Mycobacterium tuberculosis CYP121 Inhibitors

The essential enzyme CYP121 is a target for drug development against antibiotic resistant strains of Mycobacterium tuberculosis. A triazol-1-yl phenol fragment 1 was identified to bind to CYP121 using a cascade of biophysical assays. Synthetic merging and optimization of 1 produced a 100-fold improvement in binding affinity, yielding lead compound 2 (KD = 15 μM). Deconstruction of 2 into its component retrofragments allowed the group efficiency of structural motifs to be assessed, the identification of more LE scaffolds for optimization and highlighted binding affinity hotspots. Structure-guided addition of a metal-binding pharmacophore onto LE retrofragment scaffolds produced low nanomolar (KD = 15 nM) CYP121 ligands. Elaboration of these compounds to target binding hotspots in the distal active site afforded compounds with excellent selectivity against human drug-metabolizing P450s. Analysis of the factors governing ligand potency and selectivity using X-ray crystallography, UV-vis spectroscopy, and native mass spectrometry provides insight for subsequent drug development.

Bacterial Thymidylate Synthase Binds Two Molecules of Substrate and Cofactor without Cooperativity

Thymidylate synthase (TSase) is a clinically important enzyme because it catalyzes synthesis of the sole de novo source of deoxy-thymidylate. Without this enzyme, cells die a "thymineless death" since they are starved of a crucial DNA synthesis precursor. As a drug target, TSase is well studied in terms of its structure and reaction mechanism. An interesting mechanistic feature of dimeric TSase is that it is "half-the-sites reactive", which is a form of negative cooperativity. Yet, the basis for this is not well-understood. Some experiments point to cooperativity at the binding steps of the reaction cycle as being responsible for the phenomenon, but the literature contains conflicting reports. Here we use ITC and NMR to resolve these inconsistencies. This first detailed thermodynamic dissection of multisite binding of dUMP to E. coli TSase shows the nucleotide binds to the free and singly bound forms of the enzyme with nearly equal affinity over a broad range of temperatures and in multiple buffers. While small but significant differences in ΔC°P for the two binding events show that the active sites are not formally equivalent, there is little-to-no allostery at the level of ΔG°bind. In addition NMR titration data reveal that there is minor intersubunit cooperativity in formation of a ternary complex with the mechanism based inhibitor, 5F-dUMP, and cofactor. Taken together, the data show that functional communication between subunits is minimal for both binding steps of the reaction coordinate.

Intermolecular epistasis shaped the function and evolution of an ancient transcription factor and its DNA binding sites

Complexes of specifically interacting molecules, such as transcription factor proteins (TFs) and the DNA response elements (REs) they recognize, control most biological processes, but little is known concerning the functional and evolutionary effects of epistatic interactions across molecular interfaces. We experimentally characterized all combinations of genotypes in the joint protein-DNA sequence space defined by an historical transition in TF-RE specificity that occurred some 500 million years ago in the DNA-binding domain of an ancient steroid hormone receptor. We found that rampant epistasis within and between the two molecules was essential to specific TF-RE recognition and to the evolution of a novel TF-RE complex with unique derived specificity. Permissive and restrictive epistatic mutations across the TF-RE interface opened and closed potential evolutionary paths accessible by the other, making the evolution of each molecule contingent on its partner's history and allowing a molecular complex with novel specificity to evolve.

Substrate deconstruction and the nonadditivity of enzyme recognition

Predicting substrates for enzymes of unknown function is a major postgenomic challenge. Substrate discovery, like inhibitor discovery, is constrained by our ability to explore chemotypes; it would be expanded by orders of magnitude if reactive sites could be probed with fragments rather than fully elaborated substrates, as is done for inhibitor discovery. To explore the feasibility of this approach, substrates of six enzymes from three different superfamilies were deconstructed into 41 overlapping fragments that were tested for activity or binding. Surprisingly, even those fragments containing the key reactive group had little activity, and most fragments did not bind measurably, until they captured most of the substrate features. Removing a single atom from a recognized substrate could often reduce catalytic recognition by 6 log-orders. To explore recognition at atomic resolution, the structures of three fragment complexes of the β-lactamase substrate cephalothin were determined by X-ray crystallography. Substrate discovery may be difficult to reduce to the fragment level, with implications for function discovery and for the tolerance of enzymes to metabolite promiscuity. Pragmatically, this study supports the development of libraries of fully elaborated metabolites as probes for enzyme function, which currently do not exist.

Biofragments: an approach towards predicting protein function using biologically related fragments and its application to Mycobacterium tuberculosis CYP126

We present a novel fragment-based approach that tackles some of the challenges for chemical biology of predicting protein function. The general approach, which we have termed biofragments, comprises two key stages. First, a biologically relevant fragment library (biofragment library) can be designed and constructed from known sets of substrate-like ligands for a protein class of interest. Second, the library can be screened for binding to a novel putative ligand-binding protein from the same or similar class, and the characterization of hits provides insight into the basis of ligand recognition, selectivity, and function at the substrate level. As a proof-of-concept, we applied the biofragments approach to the functionally uncharacterized Mycobacterium tuberculosis (Mtb) cytochrome P450 isoform, CYP126. This led to the development of a tailored CYP biofragment library with notable 3D characteristics and a significantly higher screening hit rate (14%) than standard drug-like fragment libraries screened previously against Mtb CYP121 and 125 (4% and 1%, respectively). Biofragment hits were identified that make both substrate-like type-I and inhibitor-like type-II interactions with CYP126. A chemical-fingerprint-based substrate model was built from the hits and used to search a virtual TB metabolome, which led to the discovery that CYP126 has a strong preference for the recognition of aromatics and substrate-like type-I binding of chlorophenol moieties within the active site near the heme. Future catalytic analyses will be focused on assessing CYP126 for potential substrate oxidative dehalogenation.

Dissecting fragment-based lead discovery at the von Hippel-Lindau protein:hypoxia inducible factor 1α protein-protein interface

Fragment screening is widely used to identify attractive starting points for drug design. However, its potential and limitations to assess the tractability of often challenging protein:protein interfaces have been underexplored. Here, we address this question by means of a systematic deconstruction of lead-like inhibitors of the pVHL:HIF-1α interaction into their component fragments. Using biophysical techniques commonly employed for screening, we could only detect binding of fragments that violate the Rule of Three, are more complex than those typically screened against classical druggable targets, and occupy two adjacent binding subsites at the interface rather than just one. Analyses based on ligand and group lipophilicity efficiency of anchored fragments were applied to dissect the individual subsites and probe for binding hot spots. The implications of our findings for targeting protein interfaces by fragment-based approaches are discussed.

Fragment-based screening by protein crystallography: successes and pitfalls

Fragment-based drug discovery (FBDD) concerns the screening of low-molecular weight compounds against macromolecular targets of clinical relevance. These compounds act as starting points for the development of drugs. FBDD has evolved and grown in popularity over the past 15 years. In this paper, the rationale and technology behind the use of X-ray crystallography in fragment based screening (FBS) will be described, including fragment library design and use of synchrotron radiation and robotics for high-throughput X-ray data collection. Some recent uses of crystallography in FBS will be described in detail, including interrogation of the drug targets β-secretase, phenylethanolamine N-methyltransferase, phosphodiesterase 4A and Hsp90. These examples provide illustrations of projects where crystallography is straightforward or difficult, and where other screening methods can help overcome the limitations of crystallography necessitated by diffraction quality.

Structural insights into the molecular basis of the ligand promiscuity

Selectivity is a key factor in drug development. In this paper, we questioned the Protein Data Bank to better understand the reasons for the promiscuity of bioactive compounds. We assembled a data set of >1000 pairs of three-dimensional structures of complexes between a "drug-like" ligand (as its physicochemical properties overlap that of approved drugs) and two distinct "druggable" protein targets (as their binding sites are likely to accommodate "drug-like" ligands). Studying the similarity between the ligand-binding sites in the different targets revealed that the lack of selectivity of a ligand can be due (i) to the fact that Nature has created the same binding pocket in different proteins, which do not necessarily have otherwise sequence or fold similarity, or (ii) to specific characteristics of the ligand itself. In particular, we demonstrated that many ligands can adapt to different protein environments by changing their conformation, by using different chemical moieties to anchor to different targets, or by adopting unusual extreme binding modes (e.g., only apolar contact between the ligand and the protein, even though polar groups are present on the ligand or at the protein surface). Lastly, we provided new elements in support to the recent studies which suggest that the promiscuity of a ligand might be inferred from its molecular complexity.

Experimental and computational active site mapping as a starting point to fragment-based lead discovery

Small highly soluble probe molecules such as aniline, urea, N-methylurea, 2-bromoacetate, 1,2-propanediol, nitrous oxide, benzamidine, and phenol were soaked into crystals of various proteins to map their binding pockets and to detect hot spots of binding with respect to hydrophobic and hydrophilic properties. The selected probe molecules were first tested at the zinc protease thermolysin. They were then applied to a wider range of proteins such as protein kinase A, D-xylose isomerase, 4-diphosphocytidyl-2C-methyl-D-erythritol synthase, endothiapepsin, and secreted aspartic protease 2. The crystal structures obtained clearly show that the probe molecules populate the protein binding pockets in an ordered fashion. The thus characterized, experimentally observed hot spots of binding were subjected to computational active site mapping using HotspotsX. This approach uses knowledge-based pair potentials to detect favorable binding positions for various atom types. Good agreement between the in silico hot spot predictions and the experimentally observed positions of the polar hydrogen bond forming functional groups and hydrophobic portions was obtained. Finally, we compared the observed poses of the small-molecule probes with those of much larger structurally related ligands. They coincide remarkably well with the larger ligands, considering their spatial orientation and the experienced interaction patterns. This observation confirms the fundamental hypothesis of fragment-based lead discovery: that binding poses, even of very small molecular probes, do not significantly deviate or move once a ligand is grown further into the binding site. This underscores the fact that these probes populate given hot spots and can be regarded as relevant seeds for further design.

Computational medicinal chemistry in fragment-based drug discovery: what, how and when

The use of fragment-based drug discovery (FBDD) has increased in the last decade due to the encouraging results obtained to date. In this scenario, computational approaches, together with experimental information, play an important role to guide and speed up the process. By default, FBDD is generally considered as a constructive approach. However, such additive behavior is not always present, therefore, simple fragment maturation will not always deliver the expected results. In this review, computational approaches utilized in FBDD are reported together with real case studies, where applicability domains are exemplified, in order to analyze them, and then, maximize their performance and reliability. Thus, a proper use of these computational tools can minimize misleading conclusions, keeping the credit on FBDD strategy, as well as achieve higher impact in the drug-discovery process. FBDD goes one step beyond a simple constructive approach. A broad set of computational tools: docking, R group quantitative structure–activity relationship, fragmentation tools, fragments management tools, patents analysis and fragment-hopping, for example, can be utilized in FBDD, providing a clear positive impact if they are utilized in the proper scenario – what, how and when. An initial assessment of additive/non-additive behavior is a critical point to define the most convenient approach for fragments elaboration.

Mechanism of DNA methylation: the double role of DNA as a substrate and as a cofactor

Methylation of cytosine residues in the DNA is one of the most important epigenetic marks central to the control of differential expression of genes. We perform quantum mechanical calculations to investigate the catalytic mechanism of the bacterial HhaI DNA methyltransferase. We find that the enzyme nucleophile, Cys81, can attack C6 of cytosine only after it is deprotonated by the DNA phosphate group, a reaction facilitated by a bridging water molecule. This finding, which indicates that the DNA acts as both the substrate and the cofactor, can explain the total loss of activity observed in an analogous enzyme, thymidylate synthase, when the phosphate group of the substrate was removed. Furthermore, our results displaying the inability of the phosphate group to deprotonate the side chain of serine is in agreement with the total, or the large extent of, inactivity observed for the C81S mutant. In contrast to results from previous calculations, we find that the active site conserved residues, Glu119, Arg163, and Arg165, are crucial for catalysis. In addition, the enzyme-DNA adduct formation and the methyl transfer from the cofactor S-adenosyl-L-methionine are not concerted but proceed via stepwise mechanism. In many of the different steps of this methylation reaction, the transfer of a proton is found to be necessary. To render these processes possible, we find that several water molecules, found in the crystal structure, play an important role, acting as a bridge between the donating and accepting proton groups.

Fragment-based deconstruction of Bcl-xL inhibitors

Fragment-based drug design consists of screening low-molecular-weight compounds in order to identify low-affinity ligands that are then modified or linked to yield potent inhibitors. The method thus attempts to build bioactive molecules in a modular way and relies on the hypothesis that the fragment binding mode will be conserved upon elaboration of the active molecule. If the inverse process is considered, do the fragments resulting from the deconstruction of high-affinity inhibitors recapitulate their binding mode in the large molecule? Few studies deal with this issue. Here, we report the analysis of 22 fragments resulting from the dissection of 9 inhibitors of the antiapoptotic protein Bcl-x(L). To determine if the fragments retained affinity toward the protein and identify their binding site, ligand-observed and protein-observed NMR experiments were used. The analysis of the fragments behavior illustrates the complexity of low-affinity protein-ligand interactions involved in the fragment-based construction of bioactive molecules.

Explaining an unusually fast parasitic enzyme: folate tail-binding residues dictate substrate positioning and catalysis in Cryptosporidium hominis thymidylate synthase

The essential enzyme TS-DHFR from Cryptosporidium hominis undergoes an unusually rapid rate of catalysis at the conserved TS domain, facilitated by two nonconserved residues, Ala287 and Ser290, in the folate tail-binding region. Mutation of these two residues to their conserved counterparts drastically affects multiple steps of the TS catalytic cycle. We have determined the crystal structures of all three mutants (A287F, S290G, and A287F/S290G) in complex with active site ligands dUMP and CB3717. The structural data show two effects of the mutations: an increased distance between the ligands in the active site and increased flexibility of the folate ligand in the partially open enzyme state that precedes conformational change to the active catalytic state. The latter effect is able to be rescued by the mutants containing the A287F mutation. In addition, the conserved water network of TS is altered in each of the mutants. The structural results point to a role of the folate tail-binding residues in closely positioning ChTS ligands and restricting ligand flexibility in the partially open state to allow for a rapid transition to the active closed state and enhanced rate of catalysis. These results provide an explanation on how folate tail-binding residues at one end of the active site affect long-range interactions throughout the TS active site and validate these residues as targets for species-specific drug design.

A molecular modeling study of the interaction of 2'-fluoro-substituted analogues of dUMP/FdUMP with thymidylate synthase

Molecular dynamics simulations and free energy calculations are presented, exploring previously described experimentally studied interactions of a series of 2'-fluoro-substituted dUMP/FdUMP analogues with thymidylate synthase (TS). The results show the inhibitory behaviors of 2'-F-ara-UMP, 2',2''-diF-dUMP and 2',5-diF-ara-UMP to be dependent upon the binding positions and orientations adopted by the molecules of these compounds in the active site of TS. The binding mode of 2',5-diF-ara-UMP suggests a novel role of the active site residue Trp 80, stabilizing through hydrophobic stacking the binding position of the pyrimidine ring in 2',5-diF-ara-UMP.

Role of an invariant lysine residue in folate binding on Escherichia coli thymidylate synthase: calorimetric and crystallographic analysis of the K48Q mutant

Thymidylate synthase (TS) catalyzes the reductive methylation of deoxyuridine monophosphate (dUMP) using methylene tetrahydrofolate (CH(2)THF) as cofactor, the glutamate tail of which forms a water-mediated hydrogen bond with an invariant lysine residue of this enzyme. To understand the role of this interaction, we studied the K48Q mutant of Escherichia coli TS using structural and biophysical methods. The k(cat) of the K48Q mutant was 430-fold lower than wild-type TS in activity, while the K(m) for the (R)-stereoisomer of CH(2)THF was 300 microM, about 30-fold larger than K(m) from the wild-type TS. Affinity constants were determined using isothermal titration calorimetry, which showed that binding was reduced by one order of magnitude for folate-like TS inhibitors, such as propargyl-dideazafolate (PDDF) or compounds that distort the TS active site like BW1843U89 (U89). The crystal structure of the K48Q-dUMP complex revealed that dUMP binding is not impaired in the mutant, and that U89 in a ternary complex of K48Q-nucleotide-U89 was bound in the active site with subtle differences relative to comparable wild-type complexes. PDDF failed to form ternary complexes with K48Q and dUMP. Thermodynamic data correlated with the structural determinations, since PDDF binding was dominated by enthalpic effects while U89 had an important entropic component. In conclusion, K48 is critical for catalysis since it leads to a productive CH(2)THF binding, while mutation at this residue does not affect much the binding of inhibitors that do not make contact with this group.

Fragment-based approaches to enzyme inhibition

Fragment-based approaches have provided a new paradigm for small-molecule drug discovery. The methodology is complementary to high-throughput screening approaches, starting from fragments of low molecular complexity and high ligand efficiency, and building up to more potent inhibitors. The approach, which depends heavily on a number of biophysical techniques, is now being taken up by more groups in both industry and academia. This article describes key aspects of the process and highlights recent developments and applications.

The effect of 5-substitution in the pyrimidine ring of dUMP on the interaction with thymidylate synthase: molecular modeling and QSAR

Thymidylate synthase (TS) is a target enzyme for a number of anticancer agents including the 5-fluorouracil metabolite, FdUMP. The present paper reports on molecular modeling studies of the effect of substitution at C(5) position in the pyrimidine ring of the TS substrate, dUMP, on the binding affinity for the enzyme. The results of molecular dynamics simulations show that the binding of C(5) analogues of dUMP to TS in the binary complexes does not undergo changes, unless a substituent with a large steric effect, such as the propyl group, is involved. On the other hand, apparent differences in the binding of the TS cofactor, resulting from varying substitution at dUMP C(5), are observed in the modeled structures of the ternary complexes of TS. These binding characteristics are supplemented with a classical QSAR model quantifying the relation between the affinity for TS and the substituent electronic and steric effects of C(5) analogues of dUMP. Based on the findings from the present work, the perspectives for finding promising new C(5) analogues of dUMP as potential agents targeted against TS are discussed.

Deconstructing fragment-based inhibitor discovery

Fragment-based screens test multiple low-molecular weight molecules for binding to a target. Fragments often bind with low affinities but typically have better ligand efficiencies (DeltaG(bind)/heavy atom count) than traditional screening hits. This efficiency, combined with accompanying atomic-resolution structures, has made fragments popular starting points for drug discovery programs. Fragment-based design adopts a constructive strategy: affinity is enhanced either by cycles of functional-group addition or by joining two independent fragments together. The final inhibitor is expected to adopt the same geometry as the original fragment hit. Here we consider whether the inverse, deconstructive logic also applies--can one always parse a higher-affinity inhibitor into fragments that recapitulate the binding geometry of the larger molecule? Cocrystal structures of fragments deconstructed from a known beta-lactamase inhibitor suggest that this is not always the case.

Probing hot spots at protein-ligand binding sites: a fragment-based approach using biophysical methods

Mapping interactions at protein-ligand binding sites is an important aspect of understanding many biological reactions and a key part of drug design. In this paper, we have used a fragment-based approach to probe "hot spots" at the cofactor-binding site of a model dehydrogenase, Escherichia coli ketopantoate reductase. Our strategy involved the breaking down of NADPH (Kd = 300 nM) into smaller fragments and the biophysical characterization of their binding using WaterLOGSY NMR spectroscopy, isothermal titration calorimetry (ITC), and inhibition studies. The weak binding affinities of fragments were measured by direct ITC titrations under low c value conditions. The 2'-phosphate and the reduced nicotinamide groups were found to contribute a large part of the binding energy. A combination of ITC and site-directed mutagenesis enabled us to locate the fragments at separate hot spots on opposite ends of the cofactor-binding site. This study has identified structural determinants for cofactor recognition that represent a blueprint for future inhibitor design.

The role of protein dynamics in thymidylate synthase catalysis: variants of conserved 2'-deoxyuridine 5'-monophosphate (dUMP)-binding Tyr-261

The enzyme thymidylate synthase (TS) catalyzes the reductive methylation of 2'-deoxyuridine 5'-monophosphate (dUMP) to 2'-deoxythymidine 5'-monophosphate. Using kinetic and X-ray crystallography experiments, we have examined the role of the highly conserved Tyr-261 in the catalytic mechanism of TS. While Tyr-261 is distant from the site of methyl transfer, mutants at this position show a marked decrease in enzymatic activity. Given that Tyr-261 forms a hydrogen bond with the dUMP 3'-O, we hypothesized that this interaction would be important for substrate binding, orientation, and specificity. Our results, surprisingly, show that Tyr-261 contributes little to these features of the mechanism of TS. However, the residue is part of the structural core of closed ternary complexes of TS, and conservation of the size and shape of the Tyr side chain is essential for maintaining wild-type values of kcat/Km. Moderate increases in Km values for both the substrate and cofactor upon mutation of Tyr-261 arise mainly from destabilization of the active conformation of a loop containing a dUMP-binding arginine. Besides binding dUMP, this loop has a key role in stabilizing the closed conformation of the enzyme and in shielding the active site from the bulk solvent during catalysis. Changes to atomic vibrations in crystals of a ternary complex of Escherichia coli Tyr261Trp are associated with a greater than 2000-fold drop in kcat/Km. These results underline the important contribution of dynamics to catalysis in TS.

Crystal structures of thymidylate synthase mutant R166Q: structural basis for the nearly complete loss of catalytic activity

Thymidylate synthase (TS) catalyzes the folate-dependent methylation of deoxyuridine monophosphate (dUMP) to form thymidine monophosphate (dTMP). We have investigated the role of invariant arginine 166, one of four arginines that contact the dUMP phosphate, using site-directed mutagenesis, X-ray crystallography, and TS from Escherichia coli. The R166Q mutant was crystallized in the presence of dUMP and a structure determined to 2.9 A resolution, but neither the ligand nor the sulfate from the crystallization buffer was found in the active site. A second structure determined with crystals prepared in the presence of dUMP and the antifolate 10-propargyl-5,8-dideazafolate revealed that the inhibitor was bound in an extended, nonproductive conformation, partially occupying the nucleotide-binding site. A sulfate ion, rather than dUMP, was found in the nucleotide phosphate-binding site. Previous studies have shown that the substitution at three of the four arginines of the dUMP phosphate-binding site is permissive; however; for Arg166, all the mutations lead to a near-inactive mutant. The present structures of TS R166Q reveal that the phosphate-binding site is largely intact, but with a substantially reduced affinity for phosphate, despite the presence of the three remaining arginines. The position of Cys146, which initiates catalysis, is shifted in the mutant and resides in a position that interferes with the binding of the dUMP pyrimidine moiety.

In silico fragment-based discovery of DPP-IV S1 pocket binders

Dipeptidyl peptidase IV is a clinically validated target for type-2 diabetes and belongs to a family of peptidases with a quite unique post-proline cleavage specificity. Known inhibitors contain a limited number of molecular anchors occupying the small prototypical S1 pocket. A virtual screening approach for such S1-binding fragments was carried out using FlexX docking to evaluate its potential to confirm known and find novel compounds. Several low molecular weight inhibitors exhibiting activities in the micromolar range could be identified as starting points for structure-based design.

Calorimetric studies of ligand binding in R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) is a novel bacterial protein that possesses 222 symmetry and a single active site pore. Although the 222 symmetry implies that four symmetry-related binding sites must exist for each substrate as well as for each cofactor, various studies indicate only two molecules bind. Three possible combinations include two dihydrofolate molecules, two NADPH molecules, or one substrate plus one cofactor. The latter is the productive ternary complex. To explore the role of various ligand substituents during binding, numerous analogues, inhibitors, and fragments of NADPH and/or folate were used in both isothermal titration calorimetry (ITC) and K(i) studies. Not surprisingly, as the length of the molecule is shortened, affinity is lost, indicating that ligand connectivity is important in binding. The observed enthalpy change in ITC measurements arises from all components involved in the binding process, including proton uptake. As a buffer dependence for binding of folate was observed, this likely correlates with perturbation of the bound N3 pK(a), such that a neutral pteridine ring is preferred for pairwise interaction with the protein. Of interest, there is no enthalpic signal for binding of folate fragments such as dihydrobiopterin where the p-aminobenzoylglutamate tail has been removed, pointing to the tail as providing most of the enthalpic signal. For binding of NADPH and its analogues, the nicotinamide carboxamide is quite important. Differences between binary (binding of two identical ligands) and ternary complex formation are observed, indicating interligand pairing preferences. For example, while aminopterin and methotrexate both form binary complexes, albeit weakly, neither readily forms ternary complexes with the cofactor. These observations suggest a role for the O4 atom of folate in a pairing preference with NADPH, which ultimately facilitates catalysis.

Structures of Leishmania major pteridine reductase complexes reveal the active site features important for ligand binding and to guide inhibitor design

Pteridine reductase (PTR1) is an NADPH-dependent short-chain reductase found in parasitic trypanosomatid protozoans. The enzyme participates in the salvage of pterins and represents a target for the development of improved therapies for infections caused by these parasites. A series of crystallographic analyses of Leishmania major PTR1 are reported. Structures of the enzyme in a binary complex with the cofactor NADPH, and ternary complexes with cofactor and biopterin, 5,6-dihydrobiopterin, and 5,6,7,8-tetrahydrobiopterin reveal that PTR1 does not undergo any major conformational changes to accomplish binding and processing of substrates, and confirm that these molecules bind in a single orientation at the catalytic center suitable for two distinct reductions. Ternary complexes with cofactor and CB3717 and trimethoprim (TOP), potent inhibitors of thymidylate synthase and dihydrofolate reductase, respectively, have been characterized. The structure with CB3717 reveals that the quinazoline moiety binds in similar fashion to the pterin substrates/products and dominates interactions with the enzyme. In the complex with TOP, steric restrictions enforced on the trimethoxyphenyl substituent prevent the 2,4-diaminopyrimidine moiety from adopting the pterin mode of binding observed in dihydrofolate reductase, and explain the inhibition properties of a range of pyrimidine derivates. The molecular detail provided by these complex structures identifies the important interactions necessary to assist the structure-based development of novel enzyme inhibitors of potential therapeutic value.

Coincidence of cleavage sites of intron endonuclease I-TevI and critical sequences of the host thymidylate synthase gene

To maximize spread of their host intron or intein, many homing endonucleases recognize nucleotides that code for important and conserved amino acid residues of the target gene. Here, we examine the cleavage requirements for I-TevI, which binds a stretch of thymidylate synthase (TS) DNA that codes for functionally critical residues in the TS active site. Using an in vitro selection scheme, we identified two base-pairs in the I-TevI cleavage site region as important for cleavage efficiency. These were confirmed by comparison of I-TevI cleavage efficiencies on mutant and on wild-type substrates. We also showed that nicking of the bottom strand by I-TevI is not affected by mutation of residues surrounding the bottom-strand cleavage site, unlike other homing endonucleases. One of these two base-pairs is universally conserved in all TS sequences, and is identical with a previously identified cleavage determinant of I-BmoI, a related GIY-YIG endonuclease that binds a homologous stretch of TS-encoding DNA. The other base-pair is conserved only in a subset of TS genes that includes the I-TevI, but not the I-BmoI, target sequence. Both the I-TevI and I-BmoI cleavage site requirements correspond to functionally critical residues involved in an extensive hydrogen bond network within the TS active site. Remarkably, these cleavage requirements correlate with TS phylogeny in bacteria, suggesting that each endonuclease has individually adapted to efficiently cleave distinct TS substrates.

Relative free energies of binding to thymidylate synthase of 2- and/or 4-thio and/or 5-fluoro analogues of dUMP

Free energy perturbation calculations have been applied to evaluate the relative free energies of binding of 2'-deoxyuridine-5'-monophosphate (dUMP) and its 2- and/or 4-thio and/or 5-fluoro analogues to the wild-type E. coli thymidylate synthase (ecTS). The results accurately reproduce experimentally measured differences in the free energy of binding of dUMP versus 5-fluoro-dUMP to thymidylate synthase. They indicate that preferred binding of dUMP compared to 5-fluoro-dUMP in the binary complex is equally related to (i) more favorable electrostatic interactions of the dUMP molecule in the enzyme active site, and (ii) its less favorable solvation in the aqueous solution. The relative free energies of binding in the binary complex show moderate and qualitatively indistinguishable discrimination among the studied fluorinated and non-fluorinated 2- and/or 4-thio analogues of dUMP. The binding free energies of monothio analogues of dUMP and 5-fluoro-dUMP correspond quite well with experimentally measured activities of these nucleotides in the thymidylate synthase reaction. On the other hand, the binding free energies of both dithio analogues, 2,4-dithio-dUMP and 2,4-dithio-FdUMP, show lack of such correlation. The latter suggests that very low activities of the dithio analogues of dUMP and 5-fluoro-dUMP may relate more to the covalent reaction of these nucleotides within the ternary complex with TS and 5,10-methylenetetrahydrofolate, than to their pre-covalent binding. We speculate that a lack of substrate activity of 2,4-dithio-dUMP is related to the high aromaticity of its pyrimidine ring that prevents the Michael addition of the active site cysteine thiol to the pyrimidine C6 atom. A stronger affinity of the fluorinated analogues of dUMP to thymidylate synthase, compared to the non-fluorinated congeners, results from the fluorine substituent producing a local strain in the C6 region in the pyrimidine ring, thus sensitizing C6 to the Michael addition of the cysteine thiol.

Tethering: fragment-based drug discovery

The genomics revolution has provided a deluge of new targets for drug discovery. To facilitate the drug discovery process, many researchers are turning to fragment-based approaches to find lead molecules more efficiently. One such method, Tethering1, allows for the identification of small-molecule fragments that bind to specific regions of a protein target. These fragments can then be elaborated, combined with other molecules, or combined with one another to provide high-affinity drug leads. In this review we describe the background and theory behind Tethering and discuss its use in identifying novel inhibitors for protein targets including interleukin-2 (IL-2), thymidylate synthase (TS), protein tyrosine phosphatase 1B (PTP-1B), and caspases.

A minimalist approach to fragment-based ligand design using common rings and linkers: Application to kinase inhibitors

We present a novel method for stepwise scaffold assembly that integrates fragment-by-fragment ligand design approaches with high-throughput virtual library screening (COREGEN). As an extension of our earlier studies of common features present in drug molecules, we investigate the hypothesis that most pharmaceutically interesting ligands can be expressed in terms of the ring-linker frameworks that comprise them. Analysis of 119 published kinase inhibitors from at least 18 different targets illustrates that a basis set of 4 rings and 8 linkers is sufficient to describe approximately 90% of ring and linker occurrences, respectively. A similar result was derived from a larger set of approximately 40,000 kinase inhibitors from curated patents. A method for ring-linker-based assembly of scaffold libraries that uses experimental information to guide the placement of anchor fragments is validated using a set of reported kinase inhibitors of Bcr-Abl, Cdk2, and Src. In every case, the predominant structural motif of reported ligand cores is reproduced and variations are suggested. To underscore generality of this approach, a novel scaffold for a cyclooxygenase-2 (COX-2) selective ligand is proposed.

Conformational change in substrate binding, catalysis and product release: an open and shut case?

The role of conformational change in substrate binding, catalysis and product release is reviewed for 11 enzymes, for which crystal structures are available for the apo, substrate- and product-bound states. The extent of global conformational changes is measured, and the movements of the functional regions involved in catalysis and ligand binding are compared to the rest of the structure. We find that most of these enzymes undergo relatively small amounts of conformational change and particularly small changes in catalytic residue geometry, usually less than 1 A. In some enzymes there is significant movement of the binding residues, usually on surface loops.

Molecular recognition: the fragment approach in lead generation

The successful practice of medicinal chemistry is crucially dependent on the principles of molecular recognition: the first and "fundamental" requirement for a drug is to bind to its target; specificity, or at least selectivity, of binding is also a must. Subsequent optimization steps to develop a lead compound into a drug are a complex mixture of processes that are not yet fully understood or predictable. Fortunately, criteria exist to discard leads that would be intractable for optimization. The concepts of non-lead-likeness and lead-likeness, in respect to drug-likeness and non-drug-likeness, have prompted a rich discussion in the recent medicinal chemistry literature. The fragment approach is an emerging philosophy in the process of lead compound discovery. The basic interactions responsible for binding affinity are defined from the "protein interactions world" and key structural fragments are combined according to the criteria of three-dimensional diversity to find new leads. New techniques in screening are used for the detection of the weaker interactions of fragments with their targets that might be undetectable in classical biological assays.