Molecular and structural basis for the roles of hepatitis C virus polymerase NS5B amino acids 15, 223, and 321 in viral replication and drug resistance.

Resistance to the 2'-F-2'-C-methylguanosine monophosphate nucleotide hepatitis C virus (HCV) inhibitors PSI-352938 and PSI-353661 was associated with a combination of amino acid changes (changes of S to G at position 15 [S15G], C223H, and V321I) within the genotype 2a nonstructural protein 5B (NS5B), an RNA-dependent RNA polymerase. To understand the role of these residues in viral replication, we examined the effects of single and multiple point mutations on replication fitness and inhibition by a series of nucleotide analog inhibitors. An acidic residue at position 15 reduced replicon fitness, consistent with its proximity to the RNA template. A change of the residue at position 223 to an acidic or large residue reduced replicon fitness, consistent with its proposed proximity to the incoming nucleoside triphosphate (NTP). A change of the residue at position 321 to a charged residue was not tolerated, consistent with its position within a hydrophobic cavity. This triple resistance mutation was specific to both genotype 2a virus and 2'-F-2'-C-methylguanosine inhibitors. A crystal structure of the NS5B S15G/C223H/V321I mutant of the JFH-1 isolate exhibited rearrangement to a conformation potentially consistent with short primer-template RNA binding, which could suggest a mechanism of resistance accomplished through a change in the NS5B conformation, which was better tolerated by genotype 2a virus than by viruses of other genotypes.

Structure of hepatitis C virus polymerase in complex with primer-template RNA.

The replication of the hepatitis C viral (HCV) genome is accomplished by the NS5B RNA-dependent RNA polymerase (RdRp), for which mechanistic understanding and structure-guided drug design efforts have been hampered by its propensity to crystallize in a closed, polymerization-incompetent state. The removal of an autoinhibitory ß-hairpin loop from genotype 2a HCV NS5B increases de novo RNA synthesis by >100-fold, promotes RNA binding, and facilitated the determination of the first crystallographic structures of HCV polymerase in complex with RNA primer-template pairs. These crystal structures demonstrate the structural realignment required for primer-template recognition and elongation, provide new insights into HCV RNA synthesis at the molecular level, and may prove useful in the structure-based design of novel antiviral compounds. Additionally, our approach for obtaining the RNA primer-template-bound structure of HCV polymerase may be generally applicable to solving RNA-bound complexes for other viral RdRps that contain similar regulatory ß-hairpin loops, including bovine viral diarrhea virus, dengue virus, and West Nile virus.

Inhibition of hepatitis C virus NS5A by fluoro-olefin based γ-turn mimetics.

The HCV non-structural protein NS5A has been established as a viable target for the development of direct acting antiviral therapy. From computational modeling studies strong intra-molecular hydrogen bonds were found to be a common structural moiety within known NS5A inhibitors that have low pico-molar replicon potency. Efforts to reproduce these γ-turn-like substructures provided a novel NS5A inhibitor based on a fluoro-olefin isostere. This fluoro-olefin containing inhibitor exhibited picomolar activity (EC(50)=79 pM) against HCV genotype 1b replicon without measurable cytotoxicity. This level of activity is comparable to the natural peptide-based inhibitors currently under clinic evaluation, and demonstrates that a peptidomimetic approach can serve as a useful strategy to produce potent and structurally unique inhibitors of HCV NS5A.

Phenylpropenamide derivatives: anti-hepatitis B virus activity of the Z isomer, SAR and the search for novel analogs.

Phenylpropenamides have been reported to be a class of non-nucleoside inhibitors of the hepatitis B virus (HBV). This class of compounds was explored with the objective of developing potent anti-HBV agents, with a novel mechanism of action, that could be combined with nucleos(t)ide analogs currently used to treat HBV infection. To accomplish this objective a series of substituted arylpropenamide derivatives were prepared and the E and Z geometrical isomers were separated. The structural identity of each of the E and Z isomers was determined by single crystal X-ray crystallography. Contrary to previous reports, the activity of this class of molecules resides in the Z isomer. Further structure-activity relationship studies around the active Z isomer identified compounds that displayed potent antiviral activity against HBV with EC(90) value of approximately 0.5 µM in vitro. Attempts to develop ring constrained analogs did not lead to active HBV inhibitors.

Discovery of 5-aryloxy-2,4-thiazolidinediones as potent GPR40 agonists.

Systematic structure-activity relationship (SAR) studies of a screening lead led to the discovery of a series of thiazolidinediones (TZDs) as potent GPR40 agonists. Among them, compound C demonstrated an acute mechanism-based glucose-lowering in an intraperitoneal glucose tolerance test (IPGTT) in lean mice, while no effects were observed in GPR40 knock-out mice.

Molecular mechanism of action of pharmacoperone rescue of misrouted GPCR mutants: the GnRH receptor.

The human GnRH receptor (hGnRHR), a G protein-coupled receptor, is a useful model for studying pharmacological chaperones (pharmacoperones), drugs that rescue misfolded and misrouted protein mutants and restore them to function. This technique forms the basis of a therapeutic approach of rescuing mutants associated with human disease and restoring them to function. The present study relies on computational modeling, followed by site-directed mutagenesis, assessment of ligand binding, effector activation, and confocal microscopy. Our results show that two different chemical classes of pharmacoperones act to stabilize hGnRHR mutants by bridging residues D(98) and K(121). This ligand-mediated bridge serves as a surrogate for a naturally occurring and highly conserved salt bridge (E(90)-K(121)) that stabilizes the relation between transmembranes 2 and 3, which is required for passage of the receptor through the cellular quality control system and to the plasma membrane. Our model was used to reveal important pharmacophoric features, and then identify a novel chemical ligand, which was able to rescue a D(98) mutant of the hGnRHR that could not be rescued as effectively by previously known pharmacoperones.

What's all the FLAP about?: 5-lipoxygenase-activating protein inhibitors for inflammatory diseases.

Leukotrienes have physiological roles in innate immune responses and pathological roles in inflammatory diseases, such as asthma, allergic rhinitis and atherosclerosis. Anti-leukotriene therapy has proven benefits in the treatment of respiratory disease, either through the inhibition of leukotriene synthesis or the selective antagonism of leukotriene receptors. The first committed step in the synthesis of leukotrienes is the oxidation of arachidonic acid (AA) by 5-lipoxygenase (5-LO), and the integral membrane protein 5-lipoxygenase-activating protein (FLAP) is an essential partner of 5-LO for this process. FLAP was molecularly identified via a photoaffinity probe and an affinity gel based on MK-886, a selective leukotriene inhibitor that has no activity against broken-cell preparations of 5-LO. Several FLAP inhibitors showed efficacy in early clinical trials in asthma but were not developed commercially for unpublished reasons. Recently, the FLAP (ALOX5AP) gene has been linked to risk for myocardial infarction, stroke and restenosis, reigniting pharmaceutical interest in this target. In addition, the recent determination of the crystal structure of inhibitor-bound FLAP offers exciting potential for novel FLAP inhibitor design.

The differential interactions of peroxisome proliferator-activated receptor gamma ligands with Tyr473 is a physical basis for their unique biological activities.

Despite their proven antidiabetic efficacy, widespread use of peroxisome proliferator-activated receptor (PPAR)gamma agonists has been limited by adverse cardiovascular effects. To overcome this shortcoming, selective PPARgamma modulators (SPPARgammaMs) have been identified that have antidiabetic efficacy comparable with full agonists with improved tolerability in preclinical species. The results of structural studies support the proposition that SPPARgammaMs interact with PPARgamma differently from full agonists, thereby providing a physical basis for their novel activities. Herein, we describe a novel PPARgamma ligand, SPPARgammaM2. This compound was a partial agonist in a cell-based transcriptional activity assay, with diminished adipogenic activity and an attenuated gene signature in cultured human adipocytes. X-ray cocrystallography studies demonstrated that, unlike rosiglitazone, SPPARgammaM2 did not interact with the Tyr473 residue located within helix 12 of the ligand binding domain (LBD). Instead, SPPARgammaM2 was found to bind to and activate human PPARgamma in which the Tyr473 residue had been mutated to alanine (hPPARgammaY473A), with potencies similar to those observed with the wild-type receptor (hPPARgammaWT). In additional studies, we found that the intrinsic binding and functional potencies of structurally distinct SPPARgammaMs were not diminished by the Y473A mutation, whereas those of various thiazolidinedione (TZD) and non-TZD PPARgamma full agonists were reduced in a correlative manner. These results directly demonstrate the important role of Tyr473 in mediating the interaction of full agonists but not SPPARgammaMs with the PPARgamma LBD, thereby providing a precise molecular determinant for their differing pharmacologies.

Amino acid substitution of arginine 80 in 17beta-hydroxysteroid dehydrogenase type 3 and its effect on NADPH cofactor binding and oxidation/reduction kinetics.

17beta-Hydroxysteroid dehydrogenase type 3 (17beta-HSD-3) is a member of the short-chain dehydrogenase/reductase (SDR) family and is essential for the reductive conversion of inactive C(19)-steroid, androstenedione, to the biologically active androgen, testosterone, which plays a central role in the development of the male phenotype. Mutations that inactivate this enzyme give rise to a rare form of male pseudohermaphroditism, referred to as 17beta-HSD-3 deficiency. One such mutation is the replacement of arginine at position 80 with glutamine, compromising enzyme activity by increasing the cofactor binding constant 60-fold. In the absence of a 17beta-HSD-3 crystal structure, we have grafted its amino acid sequence for the NADPH binding site on the X-ray crystal structures of glutathione reductase (Protein Data Bank code 1gra) and 17beta-HSD type 1 (Protein Data Bank codes 1fdv and 1fdu) where we find the trunk of the arginine 80 side chain forms part of the hydrophobic pocket for the purine ring of adenosine while its guanidinium moiety interacts with the 2'-phosphate to both stabilize cofactor binding and neutralize its intrinsic negative charge through two hydrogen bonds. To qualitatively assess the role arginine 80 plays in both selecting and stabilizing NADPH binding, it was replaced with each amino acid and the mutant enzymes subjected to enzymatic analysis. There are only seven enzymes exhibiting any measurable enzymatic activity with arginine approximately lysine>leucine>glutamine>methionine>tyrosine>isoleucine. With an aspartic acid at position 58 in 17beta-HSD-3 occupying the equivalent space in the cofactor binding pocket as arginine 224 in glutathione reductase or serine 12 in 17beta-HSD-1, there was an expectation that some of the mutants might use NADH as a cofactor. In no case was NADH found to substitute for NADPH.

Web enabling technology for the design, enumeration, optimization and tracking of compound libraries.

Motivated by the need to augment Merck's in-house small molecule collection, web-based tools for designing, enumerating, optimizing and tracking compound libraries have been developed. The path leading to the current version of this Virtual Library Tool Kit (VLTK) is discussed in context of the (then) available commercial offerings and the constraints and requirements imposed by the end users. Though the effort was initiated to simplify the tasks of designing novel, drug-like and diverse compound libraries containing between 2K-10K unique entities, it has also evolved into a powerful tool for outsourcing syntheses as well as lead identification and optimization. The web tool includes components that select reagents, analyze synthons, identify backup reagents, enumerate libraries, calculate properties, optimize libraries and finally track the synthesized compounds through biological assays. In addition to accommodating project specific designs and virtual 3D library scanning, the application includes tools for parallel synthesis, laboratory automation and compound registration.

Structure-activity relationship of linear tetrapeptides Tic-DPhe-Arg-Trp-NH2 at the human melanocortin-4 receptor and effects on feeding behaviors in rat.

The melanocortin subtype-4 receptor (MC4R) has been implicated in the control of feeding behavior and body weight regulation. A series of tetrapeptides, based on Tic-DPhe-Arg-Trp-NH2-a mimic of the putative message sequence "His-Phe-Arg-Trp" and modified at the DPhe position, were prepared and pharmacologically characterized for potency and selectivity. Substitution of His with Tic gave peptides with significant increases in selectivity. The effects of the substitution pattern of DPhe were investigated and it has significant influences on potency and the level of the maximum cAMP accumulation. Intracerebroventricular administration of peptide 10 induced significant inhibition of cumulative overnight food intake and feeding duration in rats.

Modeling aided design of potent glycogen phosphorylase inhibitors.

Molecular modeling has been used to assist in the development of a novel series of potent glycogen phosphorylase inhibitors based on a phenyl diacid lead, compound 1. In the absence of suitable competitive binding assays, compound 1 was predicted to bind at the AMP allosteric site based on superposition onto known inhibitors which bind at different sites in the enzyme and analyses of the surrounding protein environment associated with these distinct sites. Possible docking modes of compound 1 at the AMP allosteric site were further explored using the crystal structure of rabbit muscle glycogen phosphorylase complexed with a Bayer diacid compound W1807 (PDB entry 3AMV). Compound 1 was predicted to interact with positively charged arginines at the AMP allosteric site in the docking model. Characterization of the binding pocket by a grid-based surface calculation of the docking model revealed a large unfilled hydrophobic region near the central phenyl ring, suggesting that compounds with larger hydrophobic groups in this region would improve binding. A series of naphthyl diacid compounds were designed and synthesized to access this hydrophobic cleft, and showed significantly improved potency.

Distinct properties and advantages of a novel peroxisome proliferator-activated protein [gamma] selective modulator.

Antidiabetic thiazolidinediones (TZDs) and non-TZD compounds have been shown to serve as agonists of the peroxisome proliferator-activated receptor gamma (PPARgamma). Here, we report the identification and characterization of a novel non-TZD selective PPARgamma modulator (nTZDpa). nTZDpa bound potently to PPARgamma with high selectivity vs. PPARalpha or PPARdelta. In cell-based assays for transcriptional activation, nTZDpa served as a selective, potent PPARgamma partial agonist and was able to antagonize the activity of PPARgamma full agonists. nTZDpa also displayed partial agonist effects when its ability to promote adipogenesis in 3T3-L1 cells was evaluated. Assessment of protein conformation using protease protection or solution nuclear magnetic resonance spectroscopy methods showed that nTZDpa produced altered PPARgamma conformational stability vs. full agonists, thereby establishing a physical basis for its observed partial agonism. DNA microarray analysis of RNA from 3T3-L1 adipocytes treated with nTZDpa or several structurally diverse PPARgamma full agonists demonstrated qualitative differences in the affected gene expression profile for nTZDpa. Chronic treatment of fat-fed, C57BL/6J mice with nTZDpa or a TZD full agonist ameliorated hyperglycemia and hyperinsulinemia. However, unlike the TZD, nTZDpa caused reductions in weight gain and adipose depot size. Feed efficiency was also substantially diminished. Unlike TZDs, nTZDpa did not cause cardiac hypertrophy in mice. When a panel of PPARgamma target genes was examined in white adipose tissue, nTZDpa produced a different in vivo expression pattern vs. the full agonist. These findings establish that novel selective PPARgamma modulators can produce altered receptor conformational stability leading to distinctive gene expression profiles, reduced adipogenic cellular effects, and potentially improved in vivo biological responses. Such compounds may lead to preferred therapies for diabetes, obesity, or metabolic syndrome.

Estrogen receptor ligands. II. Discovery of benzoxathiins as potent, selective estrogen receptor alpha modulators.

The discovery and synthesis of dihydrobenzoxathiins as potent, ERalpha subtype selective ligands are described. The most active analogue, 4-D, was found to be 50-fold selective in a competitive binding assay and 100-fold selective in a transactivation assay in HEK-293 cells. The alpha selectivity was postulated to lie in the interaction of the sulfur atom of the benzoxathiin ring with the two discriminating residues in the binding pocket of the receptor isoforms.

Design and pharmacology of N-[(3R)-1,2,3,4-tetrahydroisoquinolinium- 3-ylcarbonyl]-(1R)-1-(4-chlorobenzyl)- 2-[4-cyclohexyl-4-(1H-1,2,4-triazol- 1-ylmethyl)piperidin-1-yl]-2-oxoethylamine (1), a potent, selective, melanocortin subtype-4 receptor agonist.

Synthetic and natural peptides that act as nonselective melanocortin receptor agonists have been found to be anorexigenic and to stimulate erectile activity. We report the design and development of 1, a potent, selective (1184-fold vs MC3R, 350-fold vs MC5R), small-molecule agonist of the MC4 receptor. Pharmacological testing confirms the food intake lowering effects of MC4R agonism and suggests another role for the receptor in the stimulation of erectile activity.

Selective glucocorticoid receptor nonsteroidal ligands completely antagonize the dexamethasone mediated induction of enzymes involved in gluconeogenesis and glutamine metabolism.

Glucocorticoids (GCs) are vital multi-faceted hormones with recognized effects on carbohydrate, protein and lipid metabolism. Previous studies with the steroid antagonist, RU486 have underscored the essential role of GCs in the regulation of these metabolic pathways. This article describes the discovery and characterization of novel GRalpha selective nonsteroidal antagonists (NSGCAs). NSGCAs 2 and 3 are spirocyclic dihydropyridine derivatives that selectively bind the GRalpha with IC(50s) of 2 and 1.5 nM, respectively. Importantly, these compounds are full antagonists of the induction by dexamethasone (Dex) of marker genes for glucose and glutamine metabolism; the tyrosine amino transferase (TAT) and glutamine synthetase (GS) enzymes, respectively. In contrast, GC-dependent transcriptional repression of the collagenase 1 (MMP-1) enzyme, an established GRalpha responsive proinflammatory gene; is poorly antagonized by these compounds. These NSGCAs might have useful applications as tools in metabolic research and drug discovery.

A simple method for visualizing the differences between related receptor sites.

Pastor and Cruciani [J. Med. Chem. 38 (1995) 4637] and Kastenholz et al. [J. Med. Chem. 43 (2000) 3033] pioneered methods for comparing related receptors, with the ultimate goal of designing selective ligands. Such methods start with a reasonable superposition of high-resolution three-dimensional (3D) structures of the receptors. Next, molecular field maps are calculated for each receptor. Then the maps are analyzed to determine which map features are correlated with a particular subset of receptors. We present a method FLOGTV, based on the trend vector paradigm [J. Chem. Inf. Comput. Sci. 25 (1985) 64] to perform the analysis. This is mathematically simpler than the GRID/CPCA method of Kastenholz et al. and allows for the simultaneous comparison of many receptor structures. Also, the trend vector paradigm provides a method of selecting isopotential contours that are well above "noise". We demonstrate the method on four examples: HIV proteases versus two-domain acid proteases, thrombin versus trypsin and factor Xa, bacterial dihydrofolate reductases (DHFRs) versus vertebrate DHFRs, and P38 versus ERK protein kinases.

A simple method for visualizing the differences between related receptor sites.

Pastor and Cruciani [J. Med. Chem. 38 (1995) 4637] and Kastenholz et al. [J. Med. Chem. 43 (2000) 3033] pioneered methods for comparing related receptors, with the ultimate goal of designing selective ligands. Such methods start with a reasonable superposition of high-resolution three-dimensional (3D) structures of the receptors. Next, molecular field maps are calculated for each receptor. Then the maps are analyzed to determine which map features are correlated with a particular subset of receptors. We present a method FLOGTV, based on the trend vector paradigm [J. Chem. Inf. Comput. Sci. 25 (1985) 64] to perform the analysis. This is mathematically simpler than the GRID/CPCA method of Kastenholz et al. and allows for the simultaneous comparison of many receptor structures. Also, the trend vector paradigm provides a method of selecting isopotential contours that are well above "noise". We demonstrate the method on four examples: HIV proteases versus two-domain acid proteases, thrombin versus trypsin and factor Xa, bacterial dihydrofolate reductases (DHFRs) versus vertebrate DHFRs, and P38 versus ERK protein kinases.

Use of 2'-spirocyclic ethers in HCV nucleoside design.

Conformationally restricted 2'-spironucleosides and their prodrugs were synthesized as potential anti-HCV agents. Although the replicon activity of the new agents containing pyrimidine bases was modest, the triphosphate of a 2'-oxetane cytidine analogue demonstrated potent intrinsic biochemical activity against the NS5B polymerase, with IC50 = 8.48 µM. Activity against NS5B bearing the S282T mutation was reduced. Phosphoramidate prodrugs of a 2'-oxetane 2-amino-6-O-methyl-purine nucleoside demonstrated potent anti-HCV activity in vitro, and the corresponding triphosphate retained similar potent activity against both wild-type and S282T HCV NS5B polymerase.

Discovery of a novel class of potent HCV NS4B inhibitors: SAR studies on piperazinone derivatives.

HTS screening identified compound 2a (piperazinone derivative) as a low micromolar HCV genotype 1 (GT-1) inhibitor. Resistance mapping studies suggested that this piperazinone chemotype targets the HCV nonstructural protein NS4B. Extensive SAR studies were performed around 2a and the amide function and the C-3/C-6 cis stereochemistry of the piperazinone core were essential for HCV activity. A 10-fold increase in GT-1 potency was observed when the chiral phenylcyclopropyl amide side chain of 2a was replaced with p-fluorophenylisoxazole-carbonyl moiety (67). Replacing the C-6 nonpolar hydrophobic moiety of 67 with a phenyl moiety (95) did not diminish the GT-1 potency. A heterocyclic thiophene moiety (103) and an isoxazole moiety (108) were incorporated as isosteric replacements for the C-6 phenyl moiety (95), resulting in significant improvement in GT-1b and 1a potency. However, the piperazonone class of compounds lacks GT-2 activity and, consequently, were not pursued further into development.