Antibody modeling assessment

A blinded study to assess the state of the art in three-dimensional structure modeling of the variable region (Fv) of antibodies was conducted. Nine unpublished high-resolution x-ray Fab crystal structures covering a wide range of antigen-binding site conformations were used as benchmark to compare Fv models generated by four structure prediction methodologies. The methodologies included two homology modeling strategies independently developed by CCG (Chemical Computer Group) and Accerlys Inc, and two fully automated antibody modeling servers: PIGS (Prediction of ImmunoGlobulin Structure), based on the canonical structure model, and Rosetta Antibody Modeling, based on homology modeling and Rosetta structure prediction methodology. The benchmark structure sequences were submitted to Accelrys and CCG and a set of models for each of the nine antibody structures were generated. PIGS and Rosetta models were obtained using the default parameters of the servers. In most cases, we found good agreement between the models and x-ray structures. The average rmsd (root mean square deviation) values calculated over the backbone atoms between the models and structures were fairly consistent, around 1.2 Å. Average rmsd values of the framework and hypervariable loops with canonical structures (L1, L2, L3, H1, and H2) were close to 1.0 Å. H3 prediction yielded rmsd values around 3.0 Å for most of the models. Quality assessment of the models and the relative strengths and weaknesses of the methods are discussed. We hope this initiative will serve as a model of scientific partnership and look forward to future antibody modeling assessments.

Design, synthesis and biological evaluation of a library of thiocarbazates and their activity as cysteine protease inhibitors

Recently, we identified a novel class of potent cathepsin L inhibitors, characterized by a thiocarbazate warhead. Given the potential of these compounds to inhibit other cysteine proteases, we designed and synthesized a library of thiocarbazates containing diversity elements at three positions. Biological characterization of this library for activity against a panel of proteases indicated a significant preference for members of the papain family of cysteine proteases over serine, metallo-, and certain classes of cysteine proteases, such as caspases. Several potent inhibitors of cathepsin L and S were identified. The SAR data were employed in docking studies in an effort to understand the structural elements required for cathepsin S inhibition. This study provides the basis for the design of highly potent and selective inhibitors of the papain family of cysteine proteases.

Molecular docking of cathepsin L inhibitors in the binding site of papain

The papain/CLIK-148 coordinate system was employed as a model to study the interactions of a nonpeptide thiocarbazate inhibitor of cathepsin L ( 1). This small molecule inhibitor, a thiol ester containing a diacyl hydrazine functionality and one stereogenic center, was most active as the S-enantiomer, with an IC 50 of 56 nM; the R-enantiomer ( 2) displayed only weak activity (33 microM). Correspondingly, molecular docking studies with Extra Precision Glide revealed a correlation between score and biological activity for the two thiocarbazate enantiomers when a structural water was preserved. The molecular interactions between 1 and papain were very similar to the interactions observed for CLIK-148 ( 3a and 3b) with papain, especially with regard to the hydrogen-bonding and lipophilic interactions of the ligands with conserved residues in the catalytic binding site. Subsequent docking of virtual compounds in the binding site led to the identification of a more potent inhibitor ( 5), with an IC 50 of 7.0 nM. These docking studies revealed that favorable energy scores and correspondingly favorable biological activities could be realized when the virtual compound design included occupation of the S2, S3, and S1' subsites by hydrophobic and aromatic functionalities of the ligand, and at least three hydrogen bonding contacts between the ligand and the conserved binding site residues of the protein.

Kinetic characterization and molecular docking of a novel, potent, and selective slow-binding inhibitor of human cathepsin L

A novel small molecule thiocarbazate (PubChem SID 26681509), a potent inhibitor of human cathepsin L (EC 3.4.22.15) with an IC(50) of 56 nM, was developed after a 57,821-compound screen of the National Institutes of Health Molecular Libraries Small Molecule Repository. After a 4-h preincubation with cathepsin L, this compound became even more potent, demonstrating an IC(50) of 1.0 nM. The thiocarbazate was determined to be a slow-binding and slowly reversible competitive inhibitor. Through a transient kinetic analysis for single-step reversibility, inhibition rate constants were k(on) = 24,000 M(-1)s(-1) and k(off) = 2.2 x 10(-5) s(-1) (K(i) = 0.89 nM). Molecular docking studies were undertaken using the experimentally derived X-ray crystal structure of papain/CLIK-148 (1cvz. pdb). These studies revealed critical hydrogen bonding patterns of the thiocarbazate with key active site residues in papain. The thiocarbazate displayed 7- to 151-fold greater selectivity toward cathepsin L than papain and cathepsins B, K, V, and S with no activity against cathepsin G. The inhibitor demonstrated a lack of toxicity in human aortic endothelial cells and zebrafish. In addition, the thiocarbazate inhibited in vitro propagation of malaria parasite Plasmodium falciparum with an IC(50) of 15.4 microM and inhibited Leishmania major with an IC(50) of 12.5 microM.

Identification and characterization of 3-substituted pyrazolyl esters as alternate substrates for cathepsin B: the confounding effects of DTT and cysteine in biological assays

Substituted pyrazole esters were identified as hits in a high throughput screen (HTS) of the NIH Molecular Libraries Small Molecule Repository (MLSMR) to identify inhibitors of the enzyme cathepsin B. Members of this class, along with functional group analogs, were synthesized in an effort to define the structural requirements for activity. Analog characterization was hampered by the need to include a reducing agent such as dithiothreitol (DTT) or cysteine in the assay, highlighting the caution required in interpreting biological data gathered in the presence of such nucleophiles. Despite the confounding effects of DTT and cysteine, our studies demonstrate that the pyrazole 1 acts as alternate substrate for cathepsin B, rather than as an inhibitor.

Metastin (KiSS-1) mimetics identified from peptide structure-activity relationship-derived pharmacophores and directed small molecule database screening

Metastin, also known as KiSS-1, the cognate ligand for the metastin receptor GPR54, is a peptide known to dramatically reduce metastasis in experimental models. Despite this, there is no reported structure for metastin nor any small molecule modulators of metastin function that could be used either clinically or experimentally. Here we report the NMR solution structure of a 13-residue metastin peptide in a membrane-like environment (SDS micelles) and find it to have a relatively stable helix conformation from residues 7 to 13. In assays for metastin receptor binding and calcium flux with receptor-transfected HEK-293 cells, we demonstrate through alanine scanning and amino acid substitutions that the peptide C-terminus shows helix periodicity in an NMR structural model and that Phe9, Arg12, and Phe13 are crucial to the activity of the peptide. These three residues lie on one face of the helix and define a pharmacophore site for metastin. We used these pharmacophore features in small molecule database searches to identify hits with submicromolar affinity for the metastin receptor. We also show here that molecules mimicking key elements of this pharmacophore site bind to the metastin receptor and act as full agonists, albeit with reduced potency compared to that of metastin itself. Together this structure-activity approach may yield pharmacologically useful compounds relevant in defining and modulating metastin receptor function.

Indole-glucosides as novel sodium glucose co-transporter 2 (SGLT2) inhibitors. Part 2

A series of indole-O-glucosides and C-glucosides was synthesized and evaluated in SGLT1 and SGLT2 cell-based functional assays. Compounds 2a and 2o were identified as potent SGLT2 inhibitors and screened in ZDF rats.

Heteroaryl-O-glucosides as novel sodium glucose co-transporter 2 inhibitors. Part 1

A series of benzo-fused heteroaryl-O-glucosides was synthesized and evaluated in SGLT1 and 2 cell-based functional assays. Indole-O-glucoside 10a and benzimidazole-O-glucoside 18 exhibited potent in vitro SGLT2 inhibitory activity.

Synthesis and evaluation of 3-anilino-quinoxalinones as glycogen phosphorylase inhibitors

A series of 3-anilino-quinoxalinones has been identified as a new class of glycogen phosphorylase inhibitors. The lead compound 1 was identified through high throughput screening as well as through pharmacophore-based electronic screening. Modifications were made to the scaffold of 1 to produce novel analogues, some of which are 25 times more potent than the lead compound.

The nociceptin pharmacophore site for opioid receptor binding derived from the NMR structure and bioactivity relationships

Nociceptin, a 17 amino acid opioid-like peptide that has an inhibitory effect on synaptic transmission in the nervous system, is involved in learning, memory, attention, and emotion and is also implicated in the perception of pain and visual, auditory, and olfactory functions. In this study, we investigated the NMR solution structure of nociceptin in membrane-like environments (trifluoroethanol and SDS micelles) and found it to have a relatively stable helix conformation from residues 4-17 with functionally important N-terminal residues being folded aperidoically on top of the helix. In functional assays for receptor binding and calcium flux, alanine-scanning variants of nociceptin indicated that functionally important residues generally followed helix periodicity, consistent with the NMR structural model. Structure-activity relationships allowed identification of pharmacophore sites that were used in small molecule data base searches, affording hits with demonstrated nociceptin receptor binding affinities.

Nonpeptidic, noncovalent inhibitors of the cysteine protease cathepsin S

The first nonpeptidic, noncovalent inhibitors of the cysteine protease cathepsin S (CatS) are described. Electronic database searching using the program DOCK generated a screening set of potential CatS inhibitors from which two lead structures were identified as promising starting points for a drug discovery effort. Lead optimization afforded potent (IC(50) < 50 nM) and selective inhibitors of CatS demonstrating cellular activity and reversibility of enzyme inhibition.

Synthesis and biological evaluation of novel macrocyclic bis-7-azaindolylmaleimides as potent and highly selective glycogen synthase kinase-3 beta (GSK-3 beta) inhibitors

Palladium catalyzed cross-coupling reactions were used to synthesize two key intermediates 3 and 5 that resulted in the synthesis of novel series of macrocyclic bis-7-azaindolylmaleimides. Among the three series of macrocycles, the oxygen atom and thiophene containing linkers yielded molecules with higher inhibitory potency at GSK-3 beta (K(i)=0.011-0.079 microM) while the nitrogen atom containing linkers yielded molecules with lower potency (K(i)=0.150->1 microM). Compound 33 and 36 displayed 1-2 orders of magnitude selectivity at GSK-3 beta against CDK2, PKC beta II, Rsk3 and little or no inhibitions to the other 62 protein kinases. Compound 46 was at least 100-fold more selective towards GSK-3 beta than PKC beta II, and it had little or no activity against a panel of 65 protein kinases, almost behaved as a GSK-3 beta 'specific inhibitor'. All three compounds showed good potency in GS assay. Molecular docking studies were conducted in an attempt to rationalize the GSK-3 beta selectivity of azaindolylmaleimides. The high selectivity, inhibitory potency and cellular activities of these non-crown-ether typed molecules may provide them as a valuable pharmacological tools in elucidating the complex roles of GSK-3 beta in cell signaling pathways and the potential usage for the treatment of elevated level of GSK-3 beta involved diseases.

Synthesis and discovery of macrocyclic polyoxygenated bis-7-azaindolylmaleimides as a novel series of potent and highly selective glycogen synthase kinase-3beta inhibitors

Attempts to design the macrocyclic maleimides as selective protein kinase C gamma inhibitors led to the unexpected discovery of a novel series of potent and highly selective glycogen synthase kinase-3beta (GSK-3beta) inhibitors. Palladium-catalyzed cross-coupling reactions were used to synthesize the key intermediates 17 and 22 that resulted in the synthesis of novel macrocycles. All three macrocyclic series (bisindolyl-, mixed 7-azaindoleindolyl-, and bis-7-azaindolylmaleimides) were found to have submicromolar inhibitory potency at GSK-3beta with various degrees of selectivity toward other protein kinases. To gain the inhibitory potency at GSK-3beta, the ring sizes of these macrocycles may play a major role. To achieve the selectivity at GSK-3beta, the additional nitrogen atoms in the indole rings may contribute to a significant degree. Overall, the bis-7-azaindolylmaleimides 28 and 29 exhibited little or no inhibitions to a panel of 50 protein kinases. Compound 29 almost behaved as a GSK-3beta specific inhibitor. Both 28 and 29 displayed good potency in GS cell-based assay. Molecular docking studies were conducted in an attempt to rationalize the GSK-3beta selectivity of azaindolylmaleimides.

Structure-based combinatorial library design: methodologies and applications

Rational design of small focused libraries that are biased toward specific therapeutic targets is currently at the forefront of combinatorial library design. Various structure-based design strategies can be implemented in focused library design when the 3D structure of the target is available through X-ray or NMR determination. This review discusses the major methods and programs specifically developed for the purpose of designing combinatorial libraries under the constraint of the binding site of a biological target, with emphasis on their advantages and disadvantages. Examples of the successful application of these methodologies are highlighted, demonstrating their performances within the practical drug discovery process.

Discovery of the first non-peptide antagonist of the motilin receptor

A first-in-class non-peptide antagonist of the motilin receptor was identified through electronic screening of our corporate database against a 3D pharmacophore. The pharmacophore was developed from the motilin 22 residue endogenous peptide using NMR structural data, principles of peptide folding, and peptide structure activity relationships. The NMR data supported helical content within the peptide, and both the hydrophobic staple and N-capping box motifs were identified in the motilin sequence. The conformational features of these motifs were imposed on the peptide structure, providing a constrained conformer as a starting point for database searching. A trisubstituted cyclopentene lead was identified directly from the electronic search. Compounds in this series inhibit the binding of 125I-motilin to human antral smooth muscle membrane and antagonize motilin-induced intracellular calcium mobilization in cells expressing the human motilin receptor. A potent compound developed through optimization, RWJ 68023, is active in binding and cell-based functional assays and is also effective in inhibiting motilin-induced contractility in segments of rabbit duodenum. This orally active compound is currently undergoing clinical evaluation for the treatment of gastrointestinal disorders associated with altered motility.

NOE-derived conformation of GRGDSP cell adhesion recognition site in the presence of SDS micelles and integrin receptor GPIIB/IIIA

The tripeptide RGD is well known for its role in integrin receptor-mediated cell-cell surface adhesion. Here, NMR and transferred NOE studies have been done with the fibrinogen/fibronectin-derived hexapeptide GRGDSP in the presence of sodium dodecyl sulfate (SDS) and purified platelet glycoprotein integrin receptor GPIIb/IIIa. In the presence of SDS and absence of receptor, GRGDSP gives NOE-based distance geometry-generated structures characteristic of two "nested' beta-turns centered at RG and GD. In the presence of integrin GPIIb/IIIa, GRGDSP resonances are chemically shifted and broadened consistent with a dynamic equilibrium between free and receptor "bound' peptide. NOEs characteristic of the nested beta-turns are either absent or weaker indicating a significant conformational change in GRGDSP in the receptor bound state. GRGDSP appears to bind the receptor in a more extended backbone conformation which positions aspartic acid and arginine residues spatially close for potential electrostatic interactions.

RGD induces conformational transition in purified platelet integrin GPIIb/IIIa-SDS system yielding multiple binding states for fibrinogen gamma-chain C-terminal peptide

Fibrinogen gamma-chain C-terminal peptide HHLG-GAKQAGDV (gamma 12) and alpha-chain peptide GRGDSP are known to inhibit fibrinogen-mediated platelet cell aggregation via competitive interactions with platelet integrin receptor GPIIb/IIIa. NMR studies of gamma 12 in the presence of purified GPIIb/IIIa in SDS/water solution have demonstrated the presence of two gamma 12 binding states, one of which is eliminated by GRGDSP (RGD) up to a RGD: gamma 12 ratio of 2:1. RGD: gamma 12 ratios greater than 2:1 produce multiple sets of gamma 12 NMR signals in TOCSY spectra. At a ratio of 4:1, two to four such resonance sets can be resolved for A405, Q407, A408, G409, D410 and V411 spin systems. The number of multiple resonances remains unchanged at ratios of 6:1 and 8:1. Addition of gamma 12 to reverse the ratio to 8:8 (1:1) has no apparent effect on the RGD-induced distribution. Results suggest that RGD irreversibly induces a conformational transition(s) in GPIIb/IIIa to produce multiple gamma 12 binding sites on the receptor.

Design and evaluation of nonpeptide fibrinogen gamma-chain based GPIIb/IIIa antagonists

Two series of nonpeptide turn mimetics were designed by analysis of the solution NMR structure of the 385-411 sequence of the gamma-chain of fibrinogen. These compounds, based on the KQAGD (Lys-Gln-Ala-Gly-Asp, 406-410) sequence, were synthesized and studied in vitro. The most interesting compound from our study, RWJ 50042 (25), exhibits potent inhibition of fibrinogen binding to GPIIb/IIIa (IC50 = 0.009 microM), as well as thrombin- or collagen-induced platelet aggregation (IC50 = 0.76, 0.14 microM). Since the 400-411 sequence is required for gamma-chain bioactivity and is a unique recognition sequence among ligands for integrins, vis-a-vis other RGD (Arg-Gly-Asp)-presenting proteins, these turn mimetics may represent a new, selective approach to antagonism of the fibrinogen receptor.