Search

Free energy perturbation (FEP) methodologies have become commonplace methods for modeling potency in hit-to-lead and lead optimization stages of drug discovery. The conformational states of the initial poses of compounds for FEP+ calculations are often set up by alignment to a cocrystal structure ligand, but it is not clear if this method provides the best result for all proteins or all ligands. Not only are ligand conformational states potential variables in modeling compound potency in FEP but also the selection of crystallographic water molecules for inclusion in the FEP input structures can impact FEP models. Here, we report the results of FEP calculations using FEP+ from Schrödinger and starting from maximum common substructure alignment and docked poses generated with an array of docking methodologies. As a benchmark data set, we use monoacylglycerol lipase (MAGL), an important clinical drug target in cancer malignancy, neurological diseases, and metabolic disorders, and a set of 17 MAGL inhibitors. We found a large variation among FEP+ correlations to experimental IC 50 values depending on the method used to generate the input pose and that the inclusion of ligand-based information in the docking process, with some methods, increases the correlation between FEP+ free energies and IC 50 values. Upon analysis of the initial poses, we found that the differences in FEP+ correlations stemmed from rotation around a tertiary amide bond as well as translation of the compound toward the more hydrophobic side of the MAGL pocket. FEP+ estimation improved across all pose modeling methods when hydrogen bond constraint information was added. However, simple maximum common substructure alignment in the presence of all crystallographic water molecules outperformed all other methods in correlation between estimated and experimental IC 50 values. Taken together, these findings suggest that pose selection and crystallographic water inclusion greatly impact how well FEP+ estimated IC 50 values align with experimental IC 50 values and that modelers should benchmark a few different pose generation methodologies and different water inclusion strategies for their hit-to-lead and lead optimization drug discovery projects.

Acid-sensing ion channels (ASICs) are proton-gated cation channels widely expressed in the nervous system. ASIC gating is modulated by divalent cations as well as small molecules; however, the molecular determinants of gating modulation by divalent cations are not well understood. Previously, we identified two small molecules that bind to ASIC1a at a novel site in the acidic pocket and modulate ASIC1 gating in a manner broadly resembling divalent cations, raising the possibility that these small molecules may help to illuminate the molecular determinants of gating modulation by divalent cations. Here, we examined how these two groups of modulators might interact as well as mutational effects on ASIC1a gating and its modulation by divalent cations. Our results indicate that binding of divalent cations to an acidic pocket site plays a key role in gating modulation of the channel., (© 2024. The Author(s).)

Activated factor XI (FXIa) inhibitors are promising novel anticoagulants with low bleeding risk compared with current anticoagulants. The discovery of potent FXIa inhibitors with good oral bioavailability has been challenging. Herein, we describe our discovery effort, utilizing nonclassical interactions to improve potency, cellular permeability, and oral bioavailability by enhancing the binding while reducing polar atoms. Beginning with literature-inspired pyridine N-oxide-based FXIa inhibitor 1 , the imidazole linker was first replaced with a pyrazole moiety to establish a polar C-H···water hydrogen-bonding interaction. Then, structure-based drug design was employed to modify lead molecule 2d in the P1' and P2' regions, with substituents interacting with key residues through various nonclassical interactions. As a result, a potent FXIa inhibitor 3f ( K i = 0.17 nM) was discovered. This compound demonstrated oral bioavailability in preclinical species (rat 36.4%, dog 80.5%, and monkey 43.0%) and displayed a dose-dependent antithrombotic effect in a rabbit arteriovenous shunt model of thrombosis.

Drug discovery building blocks available commercially or within an internal inventory cover a diverse range of chemical space and yet describe only a tiny fraction of all chemically feasible reagents. Vendors will eagerly provide tools to search the former; there is no straightforward method of mining the latter. We describe a procedure and use case in assembling chemical structures not available for purchase but that could likely be synthesized in one robust chemical transformation starting from readily available building blocks. Accessing this vast virtual chemical space dramatically increases our curated collection of reagents available for medicinal chemistry exploration and novel hit generation, almost tripling the number of those with 10 or fewer atoms., Competing Interests: The authors declare no competing financial interest., (© 2021 American Chemical Society.)

Background: Accurate prediction of absorption, distribution, metabolism and excretion (ADME) properties can facilitate the identification of promising drug candidates. Methodology & Results: The authors present the Janssen generic Target Product Profile (gTPP) model, which predicts 18 early ADME properties, employs a graph convolutional neural network algorithm and was trained on between 1000-10,000 internal data points per predicted parameter. gTPP demonstrated stronger predictive power than pretrained commercial ADME models and automatic model builders. Through a novel logging method, the authors report gTPP usage for more than 200 Janssen drug discovery scientists. Conclusion: The investigators successfully enabled the rapid and systematic implementation of predictive ML tools across a drug discovery pipeline in all therapeutic areas. This experience provides useful guidance for other large-scale AI/ML deployment efforts.

Acid-sensing ion channels (ASICs) are proton-gated cation channels critical for neuronal functions. Studies of ASIC1, a major ASIC isoform and proton sensor, have identified acidic pocket, an extracellular region enriched in acidic residues, as a key participant in channel gating. While binding to this region by the venom peptide psalmotoxin modulates channel gating, molecular and structural mechanisms of ASIC gating modulation by small molecules are poorly understood. Here, combining functional, crystallographic, computational and mutational approaches, we show that two structurally distinct small molecules potently and allosterically inhibit channel activation and desensitization by binding at the acidic pocket and stabilizing the closed state of rat/chicken ASIC1. Our work identifies a previously unidentified binding site, elucidates a molecular mechanism of small molecule modulation of ASIC gating, and demonstrates directly the structural basis of such modulation, providing mechanistic and structural insight into ASIC gating, modulation and therapeutic targeting.

G-protein coupled receptor kinase 2 (GRK2), which is upregulated in the failing heart, appears to play a critical role in heart failure (HF) progression in part because enhanced GRK2 activity promotes dysfunction of β-adrenergic signaling and myocyte death. An orally bioavailable GRK2 inhibitor could offer unique therapeutic outcomes that cannot be attained by current heart failure treatments that directly target GPCRs or angiotensin-converting enzyme. Herein, we describe the discovery of a potent, selective, and orally bioavailable GRK2 inhibitor, 8h, through high-throughput screening, hit-to-lead optimization, structure-based design, molecular modelling, synthesis, and biological evaluation. In the cellular target engagement assays, 8h enhances isoproterenol-mediated cyclic adenosine 3',5'-monophosphate (cAMP) production in HEK293 cells overexpressing GRK2. Compound 8h was further evaluated in a human stem cell-derived cardiomyocyte (HSC-CM) contractility assay and potentiated isoproterenol-induced beating rate in HSC-CMs., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Artificial intelligence and machine learning have demonstrated their potential role in predictive chemistry and synthetic planning of small molecules; there are at least a few reports of companies employing in silico synthetic planning into their overall approach to accessing target molecules. A data-driven synthesis planning program is one component being developed and evaluated by the Machine Learning for Pharmaceutical Discovery and Synthesis (MLPDS) consortium, comprising MIT and 13 chemical and pharmaceutical company members. Together, we wrote this perspective to share how we think predictive models can be integrated into medicinal chemistry synthesis workflows, how they are currently used within MLPDS member companies, and the outlook for this field.

A tandem one-pot reaction featuring a cross-coupling followed by an intramolecular oxetane ring opening by mild nucleophiles is reported. The overall transformation comprises a carbon-carbon bond formation along with a carbon-heteroatom bond construction providing diverse multicyclic ring systems with a pendant hydroxymethyl handle for further elaboration. This approach constitutes a convergent method for rapid access to various scaffolds. Furthermore, a comparison of computed low-energy conformers is presented to rationalize instances in which cyclization was not observed.

The topic of gender equality within the United States workforce is receiving a great deal of attention. The field of chemistry is no exception and is increasingly focused on taking steps to achieve gender diversity within the chemistry workforce. Over the past several years, many computational chemistry groups within large pharmaceutical companies have realized growth in the number of women, and here we discuss the key factors that we believe have played a role in attracting and retaining the authors of this review as computational chemists in pharma. Furthermore, we combine our professional experiences in the context of how computational methodology and technology have evolved over the past decades and how that evolution has facilitated the inclusion of more women into the field. Our hope is to be a part of a solution and provide insight that will allow the chemistry workforce to continue to make steps forward in attaining gender diversity in the workplace.

Type 2 diabetes mellitus (T2DM) is characterized by chronically elevated plasma glucose levels. The inhibition of glucagon-induced hepatic glucose output via antagonism of the glucagon receptor (GCGR) using a small-molecule antagonist is a promising mechanism for improving glycemic control in the diabetic state. The present work discloses the discovery of indazole-based β-alanine derivatives as potent GCGR antagonists through an efficient enantioselective synthesis and structure-activity relationship (SAR) exploration and optimization. Compounds within this class exhibited excellent pharmacokinetic properties in multiple preclinical species. In an acute dog glucagon challenge test, compound 13K significantly inhibited glucagon-mediated blood glucose increase when dosed orally at 10 mg/kg., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

A novel series of indazole/indole derivatives were discovered as glucagon receptor (GCGR) antagonists through scaffold hopping based on two literature leads: MK-0893 and LY-2409021. Further structure-activity relationship (SAR) exploration and optimization led to the discovery of multiple potent GCGR antagonists with excellent pharmacokinetic properties in mice and rats, including low systemic clearance, long elimination half-life, and good oral bioavailability. These potent GCGR antagonists could be used for potential treatment of type II diabetes., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

A novel series of 6-benzhydryl-4-amino-quinolin-2-ones was discovered as cannabinoid type 1 receptor (CB 1 R) inverse agonists based on the high-throughput screening hit, compound 1a. Structure-activity relationships were studied to improve in vitro/in vivo pharmacology and restrict distribution to the peripheral circulation. We adopted several strategies such as increasing topological polar surface area, incorporating discrete polyethylene glycol side chains, and targeting P-glycoprotein (P-gp) to minimize access to the brain. Compound 6a is a P-gp substrate and a potent and highly selective CB 1 R inverse agonist, demonstrating excellent in vivo metabolic stability and a low brain to plasma ratio. However, brain receptor occupancy studies showed that compound 6a may accumulate in brain with repeat dosing. This was evidenced by compound 6a inhibiting food intake and inducing weight loss in diet-induced obese mice. Thus, a strategy based on P-gp efflux may not be adequate for peripheral restriction of the disclosed quinolinone series.

Computer-aided drug discovery activities at Janssen are carried out by scientists in the Computational Chemistry group of the Discovery Sciences organization. This perspective gives an overview of the organizational and operational structure, the science, internal and external collaborations, and the impact of the group on Drug Discovery at Janssen.

Structure-activity relationship (SAR) studies on a highly potent series of arylamide FMS inhibitors were carried out with the aim of improving FMS kinase selectivity, particularly over KIT. Potent compound 17r (FMS IC50 0.7 nM, FMS cell IC50 6.1 nM) was discovered that had good PK properties and a greater than fivefold improvement in selectivity for FMS over KIT kinase in a cellular assay relative to the previously reported clinical candidate 4. This improved selectivity was manifested in vivo by no observed decrease in circulating reticulocytes, a measure of bone safety, at the highest studied dose. Compound 17r was highly active in a mouse pharmacodynamic model and demonstrated disease-modifying effects in a dose-dependent manner in a strep cell wall-induced arthritis model of rheumatoid arthritis in rats., (Copyright © 2013 Elsevier Ltd. All rights reserved.)

A small-molecule drug discovery effort can benefit from having several chemical series. Where multiple series are not available, it is often the goal of a project to find novel scaffolds. Structural studies of ligand/protein complexes provide important information on the interactions driving binding. By generalizing these, it is possible to find molecules lacking in similarity in their connectivity yet retaining the ability to interact with the same target protein. Our studies on inhibitors of the cFMS tyrosine kinase provide a dramatic example of three different chemical series that make the same key interactions with the target protein. Collectively, these structural data provide a striking example of the pharmacophore hypothesis at work. In addition, they should prompt one to employ a broad approach when attempting scaffold hopping or any search for a novel series. It is clear that molecules that bind with similar interactions to a target need not possess 2-dimensional molecular similarity.

A class of potent inhibitors of colony-stimulating factor-1 receptor (CSF-1R or FMS), as exemplified by 8 and 21, was optimized to improve pharmacokinetic and pharmacodynamic properties and potential toxicological liabilities. Early stage absorption, distribution, metabolism, and excretion assays were employed to ensure the incorporation of druglike properties resulting in the selection of several compounds with good activity in a pharmacodynamic screening assay in mice. Further investigation, utilizing the type II collagen-induced arthritis model in mice, culminated in the selection of anti-inflammatory development candidate JNJ-28312141 (23, FMS IC(50) = 0.69 nM, cell assay IC(50) = 2.6 nM). Compound 23 also demonstrated efficacy in rat adjuvant and streptococcal cell wall-induced models of arthritis and has entered phase I clinical trials.

A high-resolution structure of a ligand-bound, soluble form of human monoglyceride lipase (MGL) is presented. The structure highlights a novel conformation of the regulatory lid-domain present in the lipase family as well as the binding mode of a pharmaceutically relevant reversible inhibitor. Analysis of the structure lacking the inhibitor indicates that the closed conformation can accommodate the native substrate 2-arachidonoyl glycerol. A model is proposed in which MGL undergoes conformational and electrostatic changes during the catalytic cycle ultimately resulting in its dissociation from the membrane upon completion of the cycle. In addition, the study outlines a successful approach to transform membrane associated proteins, which tend to aggregate upon purification, into a monomeric and soluble form., (Copyright © 2011 The Protein Society.)

Fragment-based drug discovery has emerged over the past 15 years as an effective lead discovery paradigm that is complementary to traditional high-throughput screening. The starting point for fragment-based drug discovery is the identification of low-molecular weight, typically low-affinity compounds that bind to a target of interest. These fragments can then be elaborated by growing or linking to create compounds with high affinity and selectivity. A wide variety of techniques from the computational chemistry tool chest can be applied in a fragment-based project. The computational tools are equally useful in combination with experimental-binding determination or in a completely in silico design procedure. This chapter will outline these techniques, their utility, and their validation in the design of novel lead compounds., (Copyright © 2011 Elsevier Inc. All rights reserved.)

During efforts to improve the bioavailability of FMS kinase inhibitors 1 and 2, a series of saturated and aromatic 4-heterocycles of reduced basicity were prepared and evaluated in an attempt to also improve the cardiovascular safety profile over lead arylamide 1, which possessed ion channel activity. The resultant compounds retained excellent potency and exhibited diminished ion channel activity., (2010 Elsevier Ltd. All rights reserved.)

A series of pyrido[2,3-d]pyrimidin-5-ones has been synthesized and evaluated as inhibitors of the kinase domain of macrophage colony-stimulating factor-1 receptor (FMS). FMS inhibitors may be useful in treating rheumatoid arthritis and other chronic inflammatory diseases. Structure-based optimization of the lead amide analogue 10 led to hydroxamate analogue 37, which possessed excellent potency and an improved pharmacokinetic profile. During the chronic phase of streptococcal cell wall-induced arthritis in rats, compound 37 (10, 3, and 1 mg/kg) was highly effective at reversing established joint swelling. In an adjuvant-induced arthritis model in rats, 37 prevented joint swelling partially at 10 mg/kg. In this model, osteoclastogenesis and bone erosion were prevented by low doses (1 or 0.33 mg/kg) that had minimal impact on inflammation. These data underscore the potential of FMS inhibitors to prevent erosions and reduce symptoms in rheumatoid arthritis.

An anti-inflammatory 1,2,4-phenylenetriamine-containing series of FMS inhibitors with a potential to form reactive metabolites was transformed into a series with equivalent potency by incorporation of carbon-based replacement groups. Structure-based modeling provided the framework to efficiently effect this transformation and restore potencies to previous levels. This optimization removed a risk factor for potential idiosyncratic drug reactions.

A series of 8-substituted pyrazolo[1,5-a]-1,3,5-triazines were considered as a bioisosteric replacement for the 9-substituted adenine derivatives resulting in the discovery of 8-(2-methoxybenzyl)-4-(N-methylamino)-2-n-propylpyrazolo[1,5-a]-1,3,5-triazine (14d) and 2-trifluoromethyl-8-(2-methoxybenzyl)-4-(N-methylamino)pyrazolo[1,5-a]-1,3,5-triazine (14e) as a new structural class of potent phosphodiesterase type 4 inhibitors (IC(50)=13 nM and 11 nM, respectively) with high isoenzyme selectivity. An original tandem of reactions involving a palladium-mediated cross-coupling reaction (PMCCR) of the readily available 8-iodo-2-methyl-4-(N-methyl-N-phenylamino)pyrazolo[1,5-a]-1,3,5-triazine (11a) and arylboronic acids or alkynes followed by the displacement of the N-methyl-N-phenylamino group constitute the key steps in a novel synthetic approach developed herein. The treatment of 11a-c with n-BuLi and selected aldehydes represents an interesting alternative to the PMCCR for the synthesis of benzylic derivatives 14a-i. Preliminary biological testing has shown that compounds 14d and 14e strongly inhibit LPS-induced TNFalpha release from human mononuclear cells from healthy subjects. These two compounds were selected for further biological evaluation.

A series of pyrimidinopyridones has been designed, synthesized and shown to be potent and selective inhibitors of the FMS tyrosine kinase. Introduction of an amide substituent at the 6-position of the pyridone core resulted in a significant potency increase. Compound 24 effectively inhibited in vivo LPS-induced TNF in mice greater than 80%.

A series of 3,4,6-substituted 2-quinolones has been synthesized and evaluated as inhibitors of the kinase domain of macrophage colony-stimulating factor-1 receptor (FMS). The fully optimized compound, 4-(4-ethyl-phenyl)-3-(2-methyl-3H-imidazol-4-yl)-2-quinolone-6-carbonitrile 21b, has an IC(50) of 2.5 nM in an in vitro assay and 5.0 nM in a bone marrow-derived macrophage cellular assay. Inhibition of FMS signaling in vivo was also demonstrated in a mouse pharmacodynamic model.

Complement activation has been implicated in disease states such as hereditary angioedema, ischemia-reperfusion injury, acute respiratory distress syndrome, and acute transplant rejection. Even though the complement cascade provides several protein targets for potential therapeutic intervention only two complement inhibitors have been approved so far for clinical use including anti-C5 antibodies for the treatment of paroxysmal nocturnal hemoglobinuria and purified C1-esterase inhibitor replacement therapy for the control of hereditary angioedema flares. In the present study, optimization of potency and physicochemical properties of a series of thiophene amidine-based C1s inhibitors with potential utility as intravenous agents for the inhibition of the classical pathway of complement is described.

The optimization of the arylamide lead 2 resulted in identification of a highly potent series of 2,4-disubstituted arylamides. Compound 8 (FMS kinase IC(50)=0.0008 microM) served as a proof-of-concept candidate in a collagen-induced model of arthritis in mice.

The vitronectin receptor alpha(v)beta(3) has been identified as a promising potential target for the treatment of osteoporosis, diabetic retinopathy and cancer. We have recently reported 5-substituted indoles 3-[5-[2-(5,6,7,8-tetrahydro[1,8]naphthyridin-2-yl)ethoxy]indol-1-yl]-3-(3-pyridyl)propionic acid 3 and 3-[5-[2-(5,6,7,8-tetrahydro[1,8]naphthyridin-2-yl)ethoxy]indol-1-yl]-3-(3,4-methylenedioxyphenyl)propionic acid 4, as an original series of potent alpha(v)beta(3) antagonists with subnanomolar activity. Ligand-protein docking analyses have been performed to generate binding models of three different chemical classes of known alpha(v)beta(3) antagonists with alpha(v)beta(3). Results of this docking study suggested that indoles bearing the basic tetrahydronaphthyridine group at position 4 can easily adopt the correct binding conformation and should be as potent as our current 5-substituted indole leads 3 and 4. This hypothesis was nicely demonstrated by the synthesis of a series of 1,4-disubstituted indoles through a tandem of reactions involving: (i) the N-alkylation of indoles 15 and 22 with propargyl esters and cesium fluoride, and (ii) a Heck coupling reaction between 4-bromoindole and 7-vinyl-3,4-dihydro-2H-[1,8]naphthyridine-1-carboxylic acid tert-butyl ester 12, or (iii) a reductive amination involving the N-substituted-4-aminoindole 23 and the BOC-protected tetrahydro[1,8]naphthyridine aldehyde 13. Among the compounds assayed, 3-(3-pyridyl)-3-[4-[2-(5,6,7,8-tetrahydro[1,8]naphthyridin-2-yl)ethyl]indol-1-yl]propionic acid 21 showed the most promising activity on alpha(v)beta(3) (IC(50)=0.5 nM), and was found to have the same potency as our current leads 3 and 4, while maintaining selectivity over alpha(IIb)beta(IIIa). Moreover, based on the reasonable apparent permeability coefficient in an in vitro CACO-2 cell monolayer assay (P(app) apical/basolateral=2.2 x 10(-6)cm/s, P(app) basolateral/apical=2.5 x 10(-6)cm/s), compound 21 is expected to be absorbed through the intestine in human. Thus, 1,4-disubstituted indole 21 represents a new lead for this novel class of conformationally restricted alpha(v)beta(3) antagonists. Additionally, this study validates the pharmacophore model previously postulated and provides an improved basis for further structure-based drug design in the field of alpha(v)beta(3).

Virtual screening involves the mining of small molecule databases from various sources. The small molecule databases used in virtual screening are typically processed, from simple 2D representations, to maximize their information content and to optimize them for input to the particular virtual screening technology being used. Processing interprets or adds molecular information related to connectivity, stereochemistry, protonation, tautomers and conformation. For virtual screening with an automated docking protocol, a technique that relies on specific intermolecular atom-atom contacts for ranking molecules, it is expected that the pre-processing protocol can affect the results of the docking experiment. The possible effects of processing on docking results have not been extensively studied, and this topic has only recently emerged as a significant aspect of the docking-based virtual screening process. One recent report highlights significant effects of different processing procedures on docking enrichment, while another outlines a general ligand preparation strategy. Here we survey and comment on recent practice in the field.

MAPK-activated protein kinase-2 (MAPKAPK2) regulates the synthesis of tumor necrosis factor and other cytokines and is a potential drug target for inflammatory diseases. Five protein constructs were produced in 4-10mg quantities per liter of culture media using baculovirus-infected insect cells and characterized for kinase activity, thermal stability, and ligand-binding affinity. Compared to construct 1-370, removal of the C-terminal autoinhibitory peptide in 1-338 resulted in a destabilized but partially active nonphosphorylated enzyme; phosphorylation of 1-338 by p38alpha further increased activity 12-fold. A putative constitutively active mutant, 1-370/T222E/T334E, was 6.3-fold less active than phosphorylated 1-370. ThermoFluor, an equilibrium ligand-binding assay, was used to measure nucleotide analogue affinity for various constructs. Binding of phosphorylated nucleotides was Mg(2+)-dependent. Residues 1-40 were required for high-affinity binding of ADP, ATPgammaS, staurosporine, and K252a. A mutation M138A rendered 1-370 susceptible to p38-inhibitors SB-203580 and SB-202190 with IC50 values of 17.4 and 14.1 microM, respectively. Taken together, these studies provide information on the mechanism of ligand-binding to MAPKAPK2 that can be used in the search for selective small-molecule inhibitors.

Inhibiting the classical pathway of complement activation by attenuating the proteolytic activity of the serine protease C1s is a potential strategy for the therapeutic intervention in disease states such as hereditary angioedema, ischemia-reperfusion injury, and acute transplant rejection. A series of arylsulfonylthiophene-2-carboxamidine inhibitors of C1s were synthesized and evaluated for C1s inhibitory activity. The most potent compound had a Ki of 10nM and >1000-fold selectivity over uPA, tPA, FX(a), thrombin, and plasmin.

The discovery, SAR, and X-ray crystal structure of novel biarylaminoacyl-(S)-2-cyano-pyrrolidines and biarylaminoacylthiazolidines as potent inhibitors of dipeptidyl peptidase IV (DPP IV) are reported.

We describe the synthesis and structure/activity relationship of RGD mimetics that are potent inhibitors of the integrin alpha(v)beta3. Indol-1-yl propionic acids containing a variety of basic moieties at the 5-position, as well as substitutions alpha and beta to the carboxy terminus were synthesized and evaluated. Novel compounds with improved potency have been identified.

The protein product of an essential gene of unknown function from Streptococcus pneumoniae was expressed and purified for screening in the ThermoFluor affinity screening assay. This assay can detect ligand binding to proteins of unknown function. The recombinant protein was found to be in a dimeric, native-like folded state and to unfold cooperatively. ThermoFluor was used to screen the protein against a library of 3000 compounds that were specifically selected to provide information about possible biological functions. The results of this screen identified pyridoxal phosphate and pyridoxamine phosphate as equilibrium binding ligands (K(d) approximately 50 pM, K(d) approximately 2.5 microM, respectively), consistent with an enzymatic cofactor function. Several nucleotides and nucleotide sugars were also identified as ligands of this protein. Sequence comparison with two enzymes of known structure but relatively low overall sequence homology established that several key residues directly involved in pyridoxal phosphate binding were strictly conserved. Screening a collection of generic drugs and natural products identified the antifungal compound canescin A as an irreversible covalent modifier of the enzyme. Our investigation of this protein indicates that its probable biological role is that of a nucleoside diphospho-keto-sugar aminotransferase, although the preferred keto-sugar substrate remains unknown. These experiments demonstrate the utility of a generic affinity-based ligand binding technology in decrypting possible biological functions of a protein, an approach that is both independent of and complementary to existing genomic and proteomic technologies.

The performance of several commercially available docking programs is compared in the context of virtual screening. Five different protein targets are used, each with several known ligands. The simulated screening deck comprised 1000 molecules from a cleansed version of the MDL drug data report and 49 known ligands. For many of the known ligands, crystal structures of the relevant protein-ligand complexes were available. We attempted to run experiments with each docking method that were as similar as possible. For a given docking method, hit rates were improved versus what would be expected for random selection for most protein targets. However, the ability to prioritize known ligands on the basis of docking poses that resemble known crystal structures is both method- and target-dependent.

The state of the art of various computational aspects of docking-based virtual screening of database of small molecules is presented. The review encompasses the different search algorithms and the scoring functions used in docking methods and their applications to protein and nucleic acid drug targets. Recent progress made in the development and application of methods to include target flexibility are summarized. The fundamental issues and challenges involved in comparing various docking methods are discussed. Limitations of current technologies as well as future prospects are presented.

Activation of the classical pathway of complement has been implicated in disease states such as hereditary angioedema, ischemia-reperfusion injury and acute transplant rejection. The trypsin-like serine protease C1s represents a pivotal upstream point of control in the classical pathway of complement activation and is therefore likely to be a useful target in the therapeutic intervention of these disease states. A series of thiopheneamidine-based inhibitors of C1s has been optimized to give a 70 nM inhibitor that inhibits the classical pathway of complement activation in vitro.

The discovery of novel 5,7-disubstituted[1,6]naphthyridines as potent inhibitors of Spleen Tyrosine Kinase (SYK) is discussed. The SAR reveals the necessity for a 7-aryl group with preference towards para substitution and that this in combination with 5-aminoalkylamino substituents further improved the potency of the compounds. The initial SAR as well as a survey of the other positions is discussed in detail.

A study of the S1 binding of lead 5-methylthiothiophene amidine 3, an inhibitor of urokinase-type plasminogen activator, was undertaken by the introduction of a variety of substituents at the thiophene 5-position. The 5-alkyl substituted and unsubstituted thiophenes were prepared using organolithium chemistry. Heteroatom substituents were introduced at the 5-position using a novel displacement reaction of 5-methylsulfonylthiophenes and the corresponding oxygen or sulfur anions. Small alkyl group substitution at the 5-position provided inhibitors equipotent with but possessing improved solubility.

The serine protease urokinase (uPa) has been implicated in the progression of both breast and prostate cancer. Utilizing structure based design, the synthesis of a series of substituted 4-[2-amino-1,3-thiazolyl]-thiophene-2-carboxamidines is described. Further optimization of this series by substitution of the terminal amine yielded urokinase inhibitors with excellent activities.

The structure-based design and synthesis of lactam-constrained azapeptide inhibitors of human cathepsin K are described. Enhanced stability to proteolytic cleavage over acyclic analogues is discussed.

Using binding models which were based on the X-ray crystal structure of an amino acid-based active site-spanning inhibitor complexed with cathepsin K, Cbz-leucine mimics have been developed, leading ultimately to the design of a potent cathepsin K inhibitor free of amino acid components. These mimics, which consist of alpha-substituted biphenylacetyl groups in place of Cbz-leucine moieties, effectively mimic all aspects of the Cbz-leucine moieties which are important for inhibitor binding. The predicted directions of binding for the inhibitors were confirmed by mass spectral analysis of their complexes with cathepsin K, which gave results consistent with acylation of the enzyme and loss of the acylhydrazine portion of the inhibitor which binds on the S' side of the active site. The binding models were found to be very predictive of relative inhibitor potency as well as direction of inhibitor binding. These results strengthen the validity of a strategy involving iterative cycles of structure-based design and inhibitor synthesis and evaluation for the discovery of non-peptide inhibitors.

We tested the hypotheses that the carboxylate side chain of Asp147 of the mu opioid receptor interacts with the protonated nitrogen of naltrexone and morphine and that this interaction is important for pharmacological properties of the two compounds. Mutation of Asp147 to Ala or Asn substantially reduced the affinity of naltrexone and the affinity, potency and efficacy of morphine, while the Glu mutant had similar properties as the wildtype, indicating the significant role of the carboxylate group of Asp147 in receptor binding and activation. This role could be due to its direct interaction with ligands or involvement in interhelical interactions. The unprotonated analogs of naltrexone and morphine, cyclopropylcarbonyl noroxymorphone (CPCNOM) and N-formylnormorphine (NFNM), respectively, were used to discriminate between these mechanisms. CPCNOM was much less potent as an antagonist and had substantially lower affinity for the mu receptor than naltrexone. Similarly, NFNM was unable to activate the mu receptor and had much lower affinity than morphine. These results indicate the importance of the protonated nitrogen. Notably, the D147A and D147N mutations did not appreciably affect the binding affinities of CPCNOM and NFNM. In addition, the D147E mutant had similar affinities for CPCNOM and NFNM as the D147A and D147N mu receptors. Thus, the carboxylate group of Asp147 is not important for binding of the two unprotonated compounds. These results indicate that the carboxylate group of Asp147 of the mu receptor interacts directly with the protonated nitrogen of naltrexone and morphine and this interaction is important for binding and receptor activation.

Papain has been used as a surrogate enzyme in a drug design effort to obtain potent and selective inhibitors of cathepsin K, a new member of the papain superfamily of cysteine proteases that is selectively and highly expressed in osteoclasts and is implicated in bone resorption. Here we report the crystal structures of two papain-inhibitor complexes and the rational design of novel cathepsin K inhibitors. Unlike previously known crystal structures of papain-inhibitor complexes, our papain structures show ligand binding extending deep within the S'-subsites. The two inhibitor complexes, carbobenzyloxyleucinyl-leucinyl-leucinal and carbobenzyloxy-L-leucinyl-L-leucinyl methoxymethyl ketone, were refined to 2.2- and 2.5-A resolution with R-factors of 0.190 and 0. 217, respectively. The S'-subsite interactions with the inhibitors are dominated by an aromatic-aromatic stacking and an oxygen-aromatic ring edge interaction. The knowledge of S'-subsite interactions led to a design strategy for an inhibitor spanning both subsites and yielded a novel, symmetric inhibitor selective for cathepsin K. Simultaneous exploitation of both S- and S'-sites provides a general strategy for the design of cysteine protease inhibitors having high specificity to their target enzymes.

Peptidomimetic cathepsin K inhibitors have been designed using binding models which were based on the X-ray crystal structure of an amino acid-based, active site-spanning inhibitor complexed with cathepsin K. These inhibitors, which contain a benzyloxybenzoyl group in place of a Cbz-leucine moiety, maintained good inhibitory potency relative to the amino acid-based inhibitor, and the binding models were found to be very predictive of relative inhibitor potency. The binding mode of one of the inhibitors was confirmed by X-ray crystallography, and the crystallographically determined structure is in close qualitative agreement with the initial binding model. These results strengthen the validity of a strategy involving iterative cycles of structure-based design, inhibitor synthesis and evaluation, and crystallographic structure determination for the discovery of peptidomimetic inhibitors.

Potent and selective active-site-spanning inhibitors have been designed for cathepsin K, a cysteine protease unique to osteoclasts. They act by mechanisms that involve tight binding intermediates, potentially on a hydrolytic pathway. X-ray crystallographic, MS, NMR spectroscopic, and kinetic studies of the mechanisms of inhibition indicate that different intermediates or transition states are being represented that are dependent on the conditions of measurement and the specific groups flanking the carbonyl in the inhibitor. The species observed crystallographically are most consistent with tetrahedral intermediates that may be close approximations of those that occur during substrate hydrolysis. Initial kinetic studies suggest the possibility of irreversible and reversible active-site modification. Representative inhibitors have demonstrated antiresorptive activity both in vitro and in vivo and therefore are promising leads for therapeutic agents for the treatment of osteoporosis. Expansion of these inhibitor concepts can be envisioned for the many other cysteine proteases implicated for therapeutic intervention.

We previously demonstrated that [3H]beta-funaltrexamine ([3H]beta-FNA) labeled the rat mu opioid receptor expressed in Chinese hamster ovary cells with high specificity, and [3H]beta-FNA-labeled receptors migrated as one broad band with a mass of 80 kDa. In this study, we determined the region and then the amino acid residue of the mu receptor involved in the covalent binding of [3H]beta-FNA. [3H]beta-FNA-labeled receptors were solubilized and purified to approximately 10% purity by immunoaffinity chromatography with antibodies against a C-terminal domain peptide. The site of covalent bond formation was determined to be within Ala206-Met243 by CNBr cleavage of partially purified labeled mu receptors and determinations of sizes of labeled receptor fragments. The amino acid residue of beta-FNA covalent incorporation was then determined by site-directed mutagenesis studies within this region. Mutation of Lys233 to Ala, Arg, His, and Leu completely eliminated covalent binding of [3H]beta-FNA, although these mutants bound beta-FNA with high affinity. Mutations of other amino acid residues did not affect covalent binding of [3H]beta-FNA. These results indicate that [3H]beta-FNA binds covalently to Lys233. Since [3H]beta-FNA is a rigid molecule, the information will be very useful for molecular modeling of interaction between morphinans and the mu receptor.

We describe a structural validation of the use of presentation scaffolds for control and elucidation of bioactive conformations of peptides. The protein REI-RGD34--produced by inserting the sequence RIPRGDMP into the CDR1 loop region of the immunoglobulin VL domain REI--strongly inhibits fibrinogen binding to the integrins alpha IIb beta 3 and alpha V beta 3. In the X-ray crystal structure of their protein at 2.4 A resolution, the RGD-containing loop exhibits defined electron density that is consistent with models for the bioactive conformations of ligands of these receptors based on previous small-molecule studies. Furthermore, a search of a small-molecule database with conformational information derived from the structure of REI-RGD34 identified constrained peptides and peptidomimetics known to be antagonists of the platelet receptor alpha IIb beta 3.