key: cord-0944067-ampxugzf
authors: Im, Isak; Lee, Eui Seung; Choi, Soo Jeong; Lee, Ju-Yeon; Kim, Yong-Chul
title: Structure–activity relationships of heteroaromatic esters as human rhinovirus 3C protease inhibitors
date: 2009-07-01
journal: Bioorg Med Chem Lett
DOI: 10.1016/j.bmcl.2009.04.114
sha: a99d23691f8e8353280d1da87cedcabd40df768a
doc_id: 944067
cord_uid: ampxugzf

Human rhinovirus 3C protease (HRV 3C(pro)) is known to be a promising target for development of therapeutic agents against the common cold because of the importance of the protease in viral replication as well as its expression in a large number of serotypes. To explore non-peptidic inhibitors of HRV 3C(pro), a series of novel heteroaromatic esters was synthesized and evaluated for inhibitory activity against HRV 3C(pro), to determine the structure–activity relationships. The most potent inhibitor, 7, with a 5-bromopyridinyl group, had an IC(50) value of 80 nM. In addition, the binding mode of a novel analog, 19, with the 4-hydroxyquinolinone moiety, was explored by molecular docking, suggesting a new interaction in the S1 pocket.

Human rhinovirus 3C protease (HRV 3C pro ) is known to be a promising target for development of therapeutic agents against the common cold because of the importance of the protease in viral replication as well as its expression in a large number of serotypes. To explore non-peptidic inhibitors of HRV 3C pro , a series of novel heteroaromatic esters was synthesized and evaluated for inhibitory activity against HRV 3C pro , to determine the structure-activity relationships. The most potent inhibitor, 7, with a 5-bromopyridinyl group, had an IC 50 value of 80 nM. In addition, the binding mode of a novel analog, 19, with the 4hydroxyquinolinone moiety, was explored by molecular docking, suggesting a new interaction in the S1 pocket.

Ó 2009 Elsevier Ltd. All rights reserved.

Picornaviruses include important human pathogens such as human rhinovirus (HRV), enterovirus (EV), coxsackievirus (CV), poliovirus (PV), and hepatitis A virus (HAV). HRVs are the single major cause of the common cold in populations of all ages. 1,2 Although rhinovirus infection is self-limiting, complications still occur in patients with asthma, congestive heart failure, bronchiectasis, and cystic fibrosis. 1 To date, no agent has been approved by the FDA for treatment of rhinovirus infection. In addition, the large number (>100) of known HRV serotypes makes vaccine development difficult. 3 Clinical treatments are directed toward relief of the most prominent symptoms of each clinical syndrome. Previous work on anti-rhinoviral therapeutic agents focused principally on capsid-binding anti-picornaviral compounds 4 and HRV 3C protease (3C pro ) inhibitors. 5, 6 HRV 3C pro is essential for virus replication, specifically catalyzing the cleavage of glutamine-glycine peptide bonds in a 230 kDa viral polyprotein produced by cellular translation of the viral RNA genome. HRV 3C pro belongs to the cysteine protease family, showing structural similarity to trypsin proteases, but sharing minimal homology with common mammalian enzymes. 7, 8 Because of the importance of the enzyme in viral replication, and the high conservation of enzyme active site components in all known HRV serotypes, potent and selective HRV 3C pro inhibitors could be potential broad-spectrum anti-rhinoviral therapeutic agents. [9] [10] [11] As part of an effort to develop HRV 3C pro inhibitors, peptidomimetic approaches have been widely employed based on active site structure and cleavage specificity. To date, one peptidomimetic analog, rupintrivir (AG7088) (1), shown in Figure 1 , is the only HRV 3C Pro inhibitor that has advanced to the stage of human clinical trials. Phase II clinical trials of rupintrivir showed that intranasal delivery reduced the mean total daily symptom score by 33%. 12 Recently, several efforts to develop orally bioavailable anti-rhinoviral agents have been reported; these include a rupintrivir analog, 2, 13 with an improved pharmacokinetic profile, and the non-peptidic inhibitors 3 and 4. [14] [15] [16] The crystal structure of severe acute respiratory syndrome (SARS) coronavirus 3C-like protease (3CL pro ) 17 revealed a high similarity of active site architecture to that of HRV 3C pro , suggesting that picornavirus 3C pro and SARS 3CL pro inhibitors should cross-react. In fact, SARS 3CL pro inhibitors were developed by modifying rupintrivir 18 and common inhibitors with IC 50 values in the micromolar range were recently identified for both picornavirus 3C pro and coronavirus 3CL pro using high-throughput screening. 19, 20 The pyridinyl thiophene ester 5 was discovered by screening a small molecule library (50,000 molecules) and has significant inhibitory activity against both HAV 3C pro and SARS 3CL pro (IC 50 = 0.5 lM). 19, 21 Further work on a series of heteroaromatic esters resulted in the discovery of the most potent ester 6, inhibiting SARS 3CL pro with an IC 50 of 30 nM and an antiviral activity (EC 50 value) of 6.9 lM. 22 Another analog, 7, also showed potent inhibition of HAV 3C pro with an IC 50 of 53 nM. 23 Studies on the inhibitory and molecular docking mechanisms of heteroaromatic esters have shown that the 5-halopyridine moiety is crucial in binding to the S1 pocket, because both 3C pro and 3CL pro show high specificity for Gln as the P1 residue. Herein, we report on the structure-activity relationships of heteroaromatic esters with respect to inhibitory activity against HRV 3C pro . Thirty-one heteroaromatic esters, including 15 newly designed derivatives, were synthesized and screened against HRV 3C pro using a fluorescent peptide substrate. We also performed molecular docking studies to understand the potential inhibitor binding modes to the active site of HRV 3C pro .

A general synthetic method for heteroaromatic esters is described in Scheme 1. A series of compounds was synthesized by esterification of various carboxylic acids and hydroxy-heteroaromatic groups. Initial replacement of the pyridine moiety with various heteroaromatic rings was performed with retention of the 2-furoyl group. Commercially available hydroxy-heteroaromatic rings were esterified with 2-furoyl chloride in pyridine and CH 2 Cl 2 (7-18). Compound 19 was prepared from 4-hydroxyquinolinone (39) in two steps (Scheme 2). The structure and inhibitory activities of furan-2-carboxylate analogs (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) are shown in Table  1 . As a second series of analogs, with 5-halopyridyl esters, we used 5-chloro (or -bromo)-3-pyridinol for esterification with diverse carboxylic acid moieties. Compounds 20-29 were prepared with commercially available carboxylic acids. Compound 30 was synthesized by arylation of the -N-H group of the imidazole ring of 28, in a reaction with 4-trifluoromethylphenyl boronic acid followed by pyridine, with use of cupric acetate and 4-Å molecular sieves in CH 2 Cl 2 . 24 The synthesis of the 2-substituted oxazole-carboxylic acids of compounds 31-35 is outlined in Scheme 3. Briefly, acylation of 2-amino diethyl malonate (40a-e) using various acid chlorides followed by dehydrative cyclization employing trifluoroacetic anhydride or trichloroacetyl chloride under microwave conditions at 160°C for 5-10 min afforded the oxazoles 41a-e, 25 which were subjected to basic hydrolysis to yield the corresponding oxazole-carboxylic acids 42a-e. The structures and activities of the 5-chloro/bromopyridyl ester analogs are shown in Table 2 .

A FRET-based assay system was used to measure the enzyme activity of HRV16 3C pro , expressed from cDNA 26 , using the fluorogenic substrate NMN-ALFQGPPVK-DNP. 27 As shown in Table 1 , the 5-bromopyridine moiety of compound 7, a representative analog of the heteroaromatic inhibitors, was replaced by various heteroaromatic groups. The 3-pyridyl furonate, 8, showed weak inhibition of HRV 3C pro (8% inhibition at 1 lM), whereas the 2-and 4-pyridinyl furonates, 9 and 10 showed no inhibitory effects. Chloro-and bromo-substituents at the 5 position of the pyridinyl ring significantly increased inhibitory effects (compounds 7, 11), with compound 7 showing the most potent activity with an IC 50 value of 80 nM. 5-Halopyridine esters functioned well as HRV 3C pro inhibitors as a result of S1 pocket similarities with both HAV 3C pro and SARS 3CL pro .

Molecular docking of compound 7 into the active site of X-ray crystal structure of HRV 3C pro (PDB: 1CQQ) was performed by a CHARMm force field based docking tool, CDOCKER of Discovery Studio 2.0 (Accelrys Inc., San Diego, CA). In a rigid receptor structure, random orientations of ligand conformations were subjected to simulated annealing molecular dynamics. Final minimized poses of docked ligand were sorted by CHARMm energy (interaction energy plus ligand strain). Mass spectrometry and HPLC studies have already revealed the inhibition mechanism of heteroaromatic esters acting in the active site of SARS 3CL pro,20,22 Such esters initially bind competitively and strongly to the active site, being recognized by S1 pocket. The carbonyl carbon is then attacked by the nucleophilic Cys residue, and become covalently connected with the protease prior to eventual release. Thus in our study, possible prehydrolysis-binding complex models of 7 was generated (Fig. 2) . The lone pair electron of nitrogen moiety of 5-bromopyridine in 7 occupies the S1 binding pocket, forming a hydrogen bond with N e2 of His161 (3.1 Å). The bromo group occupies S1 pocket near the protein surface and, notably, a hydrogen bond with the -NH of Gly166 (3.0 Å) was observed, explaining the importance of the presence of a halogen group on the pyridine ring. The ester carbonyl oxygen forms a hydrogen bond with the -NH of Gly164 (2.5 Å). This predicted binding model is consistent with the binding mode of halopyridinyl esters to SARS 3CL pro and HAV 3C pro .

To explore the possibility of new binding interaction, new motives for the substitution of 2-position of the pyridine ring were designed for pointing toward N e2 of His161 in S1 pocket. Therefore, 2-substituted pyridinyl compounds (12) (13) (14) (15) (16) were investigated. Amine (14) and amide (16) substituted compounds showed moderate activity at 1 lM in the absence of chlorine at the 5 position.

The new S1 pocket binding moiety was also evaluated using para, meta benzamide (17, 18) and 4-quinolinone (19) , instead of 5-halopyridine. Compounds 17-19 may possibly mimic Gln, the substrate cleavage moiety. While 17 showed weak inhibitiory effect (5% inhibition at 1 lM, 21% inhibition at 10 lM), it was notable that compound 19 showed a potent inhibitory effect at 1 lM (98% inhibition), with an IC 50 value of 200 nM.

To gain molecular insight into the binding mode of the 4-quinolinone moiety, compound 19 was also docked in the HRV 3C pro active site (Fig. 3) . The oxygen of quinolinone forms a hydrogen bond with N e2 of His161 (2.8 Å) and the benzene ring of quinolinone occupies the outer space of the S1 pocket as does the bromo group of 5-bromopyridine, but there is no hydrogen bond with the -NH of Gly166. The hydrogen bond of the ester carbonyl oxygen with the -NH of Gly164 (2.0 Å) was also conserved as well. Although there is different spatial position of the 4-quinolinone carbonyl oxygen of 19 from the pyridyl nitrogen of 7, the requirement of hydrogen bond angle between the carbonyl oxygen and N e2 of His161 was satisfied with the ideal hydrogen bond at 30-60°to the O@C axis within 30°of the carbonyl plane, while the sp 2 lone pair of pyridyl nitrogen of 7 is lying toward N e2 of His161 forming optimal hydrogen bond angle (180°). This precise hydrogen bond with His161 seems to be the key initial binding unit which mimics the Gln moiety (the P1 residue), which resulted in the dramatic inhibitory activity changes depending on the substituents R 1 . Interaction energy of compound 7$19 with 3C pro was calculated by CDOCKER program. The most potent compounds 7 and 19 showed À24.5 and À26.3 kcal/mol, respectively. Comparison of the interaction energy of the ortho (9, À21.9 kcal/mol), meta (8, À22.2 kcal/mol), and para (10, À20.0 kcal/mol) position of the pyridine nitrogen showed a parallel results with the biological activities, suggesting that the meta position could prefer the orientation for the hydrogen bond. Compound 12-16 with additional hydrogen bond acceptor at 2-position of the pyridine ring might disturb the optimal hydrogen bond of the pyridyl nitrogen showing the less interaction energies (À19.9 kcal/mol to À23.5 kcal/mol). Although weak inhibitors of benzamide analogs, 17-18 could form hydrogen bonds with His161 (À24.2 kcal/mol for 17, À22.6 kcal/ mol for 18), the distance between ester carbonyl carbon and the nucleophilic -SH group of Cys might be changed unfavorably resulting in weak or no inhibitory activity.

To search for effective moieties other than the 2-furoyl group, a series of 5-halo-pyridinyl esters from various carboxylic acids was synthesized and tested. This R 2 carboxylic acids were expected to provide site specificity at S2 hydrophobic pocket and affect the covalently connected binding mode at the active site. Most compounds showed moderate-to-good inhibitory effects at 1 lM except for 29 and 31 (Table 2) . Compounds with thiophen-2carbonyl (20) , benzoyl (21) , phenylpropanoyl groups (36), and cinnamoyl (37) showed lower activities than did the 2-furoyl analogs (7 and 11) . Substitution of the 5 position of the furan ring with aromatic groups allowed retention of good activity (22) (23) (24) (25) . The steric effect of the additional aromatic groups could stabilize the postreaction state by p-stacking interaction with His40 22 rather than tight binding to S2 pocket. The 2-naphthoyl (26), 1-naphthoyl (27) , and imidazole (28) groups were useful building blocks, showing potent inhibitory activities (IC 50 of 290 nM for 28). However, arylation of the imidazole ring of 28 showed twofold decrease in activity (30), which could be caused by unfavorable constraint compared to furan ring.

In further efforts to replace the furoyl ring with other heterocyclic carboxylate moieties, isoxazole and oxazole groups were investigated. In the case of 3-methylisooxazole derivative, 29, the replaced position of furan oxygen by carbon atom resulted in the loss of activity significantly. However, oxazole derivatives (31-35) demonstrated a broad range of inhibitory activities depending on substitutions at the 2 position of the oxazole group. A cinnamyloxazole analog, 34, showed the highest activity among these compounds with 87% inhibition at 1 lM and an IC 50 value of 690 nM. Lower electron density of oxazol ring may result in weaker binding affinity than furan or imidazol moiety, but additional hydrophobic phenyl group in a proper position connected to 2oxazolic position with two carbon chain (34) significantly enhanced the inhibitory activity compared to the compounds with shorter chains or bulky aromatic groups (31-33, 35).

In conclusion, 31 heteroaromatic esters were synthesized and screened as non-peptidic inhibitors against HRV 3C pro . Compound 7, which was one of the most potent inhibitors in an earlier series, with activity against both SARS 3CL pro and HAV 3C pro , was also found to be the most potent in the current work (IC 50 of 80 nM). Substitution of 5-halopyridine with various heteroaromatic rings resulted in the discovery of the 4-quinolinone group as an alternative key structure. Further optimization of 4-quinolinone ester analogs as Gln skeleton mimics is underway, with the aim of enhancing inhibitory activities.

Fields Virology

Principles and Practices of Infectious Diseases

Viral Infections of Humans

Fields Virology

The HRV16 3C pro coding region was amplified by PCR, subcloned into the pTYB12 expression vector

A protease stock solution was maintained in 20 mM HEPES/ NaOH, 100 mM NaCl, 1 mM EDTA, and 1 mM dithiothreitol

Assays were performed at 30°C in 96-well microplates in reaction volumes of 100 lL with 50 mM Tris (pH 7.6), 1 mM EDTA, 50 lM substrate

Inhibitors and enzymes were incubated for 10 min in reaction buffer and reactions were initiated by addition of FRET-substrate. Fluorescence values were monitored at 340 nm (excitation) and 440 nm (emission)

IC 50 values of six potent inhibitors were determined using various concentrations with three independent experiments