key: cord-0976578-mxnx4mmu
authors: Wang, Li; Bao, Bo-Bo; Song, Guo-Qing; Chen, Cheng; Zhang, Xu-Meng; Lu, Wei; Wang, Zefang; Cai, Yan; Li, Shuang; Fu, Sheng; Song, Fu-Hang; Yang, Haitao; Wang, Jian-Guo
title: Discovery of unsymmetrical aromatic disulfides as novel inhibitors of SARS-CoV main protease: Chemical synthesis, biological evaluation, molecular docking and 3D-QSAR study
date: 2017-09-08
journal: Eur J Med Chem
DOI: 10.1016/j.ejmech.2017.05.045
sha: 8c91f04711723df4b384552acd81f3d05569d36d
doc_id: 976578
cord_uid: mxnx4mmu

The worldwide outbreak of severe acute respiratory syndrome (SARS) in 2003 had caused a high rate of mortality. Main protease (M(pro)) of SARS-associated coronavirus (SARS-CoV) is an important target to discover pharmaceutical compounds for the therapy of this life-threatening disease. During the course of screening new anti-SARS agents, we have identified that a series of unsymmetrical aromatic disulfides inhibited SARS-CoV M(pro) significantly for the first time. Herein, 40 novel unsymmetrical aromatic disulfides were synthesized chemically and their biological activities were evaluated in vitro against SARS-CoV M(pro). These novel compounds displayed excellent IC(50) data in the range of 0.516–5.954 μM. Preliminary studies indicated that these disulfides are reversible and mpetitive inhibitors. A possible binding mode was generated via molecular docking simulation and a comparative field analysis (CoMFA) model was constructed to understand the structure-activity relationships. The present research therefore has provided some meaningful guidance to design and identify anti-SARS drugs with totally new chemical structures.

Severe acute respiratory syndrome (SARS) is a highly infective respiratory disease caused by SARS coronavirus (SARS-CoV). Its sudden emergence and rapid outbreak during 2002e2003 had resulted in~800 deaths among >8000 reported individual cases worldwide [1] . Although the SARS epidemic had been under control for years, reemergence of this threatening illness is still a possible risk and potentially new strains of SARS can be more dangerous than the previous ones. A number of important targets have been recognized to take part in the biological events critical to SARS-CoV replication, among which a papain-like protease (PL pro ) and a chymotrypsin-like protease (3CL pro ) are of significant importance to design anti-SARS inhibitors [2] . The 3CL pro , also known as the main protease (M pro ), has attracted much attention, which could be revealed from numerous publications about novel inhibitor discovery. Crystal structures of SARS-CoV M pro , either free enzyme alone or in complex with an inhibitor, had been determined to facilitate the structural and functional investigation of this protease [3, 4] . The active site of SARS-CoV M pro contains Cys145 and His41 to constitute a catalytic dyad, in which cysteine functions as the common nucleophile in the proteolytic process.

Biological active inhibitors against SARS-CoV M pro have been reported mainly from two different approaches: one is screening large library to identify new active compounds using highthroughput technique, the other is novel inhibitor design based on the substrate structure or active site properties rationally [5] . The inhibitory activities of these compounds were then validated by in vitro protease assays. In most cases the kinetic study indicated that the inhibitor is involved in an irreversible process by forming a covalent bond with Cys145, while in some other cases the inhibition is actually a reversible behavior. The reported SARS-CoV M pro inhibitors covered a variety of different chemical scaffolds, which contain peptidomimetic compounds, 3-quinoline carboxylic acid derivatives, thiophene-2-carboxylate derivatives, zinc-conjugated compounds, cinanserin, calmodulin, keto-glutamine analogues, anilide, bifunctional boronic acid compounds, isatin derivatives, etacrynic acid derivatives, serine derivatives, trifluoromethyl ketones, acetamides, pyrazolone and quercertins [5e12] . It is a pity that research on drugs and vaccines towards SARS or SARS-like coronavirus has not brought any candidate for clinical use. Hence there still exists an urgent need to discover and identify new SARS-CoV M pro agents, especially those compounds from totally new chemical families, to develop effective therapy against this fatal viral infection.

Disulfide bonds play essential roles for bioactive proteins to keep correct folding [13] . There are a few cases that simple disulfides such as diallyl disulfide and dimethyl disulfide exhibit hypochlorous acid scavenging activity and tyrosinase inhibitory activity ( Fig. 1A) [14, 15] . The unsymmetrical disulfide compounds are useful tools in the research of dynamic combinatorial chemistry [16] . These compounds have also been reported to display a variety of biological activities. For examples, Turos et al. reported that some unsymmetrical aryl-alkyl disulfides were inhibitors of methicillinresistant Staphylococcus aureus and Bacillus anthracis (Fig. 1B ) [17] , while Khosla and co-workers published some unsymmetrical disulfides that could selectively inhibit extracellular thioredoxin (Fig. 1C ) [18] . Yoon et al. showed that some unsymmetrical disulfide compounds were inhibitors of Mycobacterium tuberculosis and Haemophilus influenzae, by interfering with acetohydroxyacid synthase (AHAS), a key enzyme in the biosynthesis pathway of branched chain amino-acids (Fig. 1D) [19, 20] . In the past of our research, we found that some unsymmetrical aromatic disulfides could inhibit plant AHAS and was useful for herbicide research (Fig. 1E) [21, 22] .

There are indeed a few reports that some aromatic disulfides exhibit antiviral activity. The virucidal activity of NSC4492 was due to targeting of arenavirus RNA synthesis (Fig. 1F ) [23] , while 5,5'dithiobis-2-nitrobenzoic acid (DTNB) had antiviral propierties against T-tropic human immunodeficiency virus type 1 (HIV-1) (Fig. 1G ) [24] . An antiviral disulfide NSC20625 compound could block interaction between arenavirus Z protein and cellular promyelocytic leukemia protein (Fig. 1H ) [25] . However, the antiviral activities against arenavirus might be different from that against coronavirus. There is no evidence that the reported antiviral activities of aromatic disulfides have any direct relationships with SARS-CoV M pro inhibition.

In an effort to discover novel inhibitors of SARS-CoV M pro , we have synthesized a series of novel unsymmetrical aromatic disulfides and evaluated their biological activities in this study. The target compounds could inhibit main protease of SARS-CoV remarkably, with the best IC 50 value of 0.516 mM. Subsequent enzymatic kinetics study indicated that the aromatic disulfides acted as reversible and non-competitive inhibitors. Therefore we have demonstrated that unsymmetrical disulfide compounds with aromatic rings are novel inhibitors of SARS-CoV M pro from a totally new chemical family, which will provide helpful information for further drug discovery.

The target unsymmetrical disulfides were synthesized by the reaction of various substituted 2-mercapto- [1, 3, 4] oxadiazole, substituted 2-mercapto-thiazole, substituted 2-mercapto-1Himidazole or substituted 2-mercapto-pyrimidine with substituted arenesulfenyl chloride in ethyl ether under very mild condition as reported previously [21, 22] . It is a quite straightforward nucleophilic substitution, in which the thiol group in the mercapto compound serves as a nucleophilic reagent and attacks the sulfur atom in the arenesulfenyl chloride. Most yields for the reactions were satisfactory, showing that this is a simple and easy procedure to produce unsymmetrical aromatic disulfides, superior to the synthetic route developed by Hahn et al. [26] . In Hahn's paper, moist tetrahydrofuran gave better yields than dried benzene as the solvent for the reaction, so the authors added 5e10 M equivalent of water to the reaction mixture and observed satisfactory yields. Thus, the presence of water was thought to be preferable for the reaction by Hahn et al. [26] . In our synthesis experiment, only ethyl ether was used as the solvent and no additional water was added to the reaction, nevertheless, we observed very high yields for most of the target compounds.

The eS-Sbond had been confirmed in our previous paper [21] . Compound 3-4, 3-6, 3-7, 3-8, 3-9, 3-10 and 3-11 were further acetylated from corresponding parent compounds that had been published by us [21] . The molecular structures of the compounds are listed in Table 1 . The title compounds were fully characterized by 1 H NMR, 13 C NMR and HRMS (Part 1 of the supplementary data for the original figures).

All the synthesized unsymmetrical aromatic disulfides were subjected to the in vitro assay of SARS-CoV M pro . The results are also illustrated in Table 1 , expressed by IC 50 values. It could be seen that the target compounds exhibited encouraging biological potency, with excellent IC 50 values ranging from 0.516 mM to 5.954 mM (The inhibition curves of all target compounds can be found in part 2 of the supplementary data). This was surprising due to the fact that not any research group had ever identified the disulfide compounds as inhibitors of SARS-CoV M pro , not to mention such strong inhibition.

Inhibition type of the disulfides was determined by means of enzymatic kinetic study, for which 3-31 and 3-39 were used. From Fig. 2 it can be seen that in the plot of enzyme concentration versus reaction velocity, the lines represent different inhibitor concentrations intersect at a same point, suggesting that the inhibition is actually a reversible action. We then measured the enzymatic velocity of SARS-CoV M pro versus substrate concentrations in the presence of either 3-31 or 3-39 (Part 3 of the supplementary data). The lines displayed in reciprocal plots intersect at a same point, indicating that both inhibitors serve as a non-competitive inhibitor with a < 1 (Part 3 of the supplementary data) [27] . On this basis, the kinetic parameters (aK i ¼ 0.20 mM, K i ¼ 0.24 mM) of 3-31 were determined [28] (Fig. 3) , which clearly proved that, the noncompetitive inhibitor 3-31 is characterized by smaller equilibrium-binding constant compared to some known inhibitors such as N3 (K i ¼ 9.0 mM) and N9 (K i ¼ 6.7 mM) [3] .

Since there is a cysteine in the active site of SARS-CoV M pro (Cys145), which plays an essential role for the biological activity of this protease, it is possible that the disulfide compound reacts with Cys145 to form a new eS-Sbond and results in a loss of enzyme activity. It is known that, if a disulfide reacts with another thiol to give a new disulfide, the thiols that are parts of the old disulfide can also react directly with this thiol to form the same new disulfide [29] . We tested the biological activities of different aryl thiols derived from our disulfides, and no inhibition of SARS-CoV M pro could be detected for any of them even at very high concentration. Accordingly it seemed unlikely that Cys145 formed a eS-Sbond by reacting with the target disulfide compounds. Another method to rule out this possibility is to determine the change of the molecular weight of the protein before and after inhibition [30, 31] . If the disulfide compound reacts with Cys145, the molecular weight would have a shift after SARS-CoV M pro is inhibited, and this is approximately the mass of half moiety of the unsymmetrical disulfide. Bearing this in mind, we measured the molecular weight of protease before and after inhibition by three disulfide compounds with significant structure difference (3-8, 3-31 and 3-39, data shown in part 4 of the supplementary data). However, no such assumed shift in the molecular weight was observed to support this idea.

For another possibility, Khosla et al. had reported selective inhibition of extracellular thioredoxin by unsymmetrical disulfides [18] , in which two cysteine residues in a close distance form an intramolecular disulfide bond. After careful analysis of the SARS-CoV M pro crystal structure (pdb entry 2AMD) [3] , no other cysteine residue was found to be in a nearby space of Cys145. In fact not any two cysteine residues are in a reasonable distance to form possible intramolecular disulfide bond. Thus we also denied this probable inhibitory mechanism.

Since the unsymmetrical aromatic disulfides do not react with the residues of the SARS-CoV M pro , it means that these compounds act as intact molecules when inhibiting this protease. Therefore, in silico molecular docking technique was utilized to predict possible binding modes of the disulfide compounds with SARS-CoV M pro . In our previous study, we had docked a small library of isatin compounds to the active site of AHAS and probable binding modes were predicted by FlexX [32e37]. Here similar strategy was adopted to carry out database docking. After investigation of the resulting docked conformations, nineteen compounds were found to overlay one another quite well (Part 5 of the supplementary data), while all the other twenty-one compounds were in unreasonable binding space, and they failed to overlay well with one another. Thus, the docked conformations of the nineteen compounds were thought to be the possible binding conformations in this study. Compound 3-31 was chosen to depict the binding mode of the disulfide inhibitors. Fig. 4 is a two-dimensional illustration of the interactions between the inhibitor and the surrounding residues of SARS-CoV M pro drawn by LIGPLOT [38] . The location of 3-31 has some overlap with the inhibitor N9, the inhibitor in the original pdb file. The compound binds with SARS-CoV M pro via multiple hydrogen bonding contacts and hydrophobic contacts. Phe140, Leu141, His163, Met165, Glu166 and His172 form hydrophobic interactions with the small molecule, while Asn142, Gly143 and Cys145 form intermolecular H-bond with the inhibitor. The predicted binding mode therefore provides a useful clue to understand the possible molecular basis of these inhibitors.

Comparative field analysis (CoMFA) is a tool to generate 3D contour models to quantitatively analyze the structure-activity relationships of bioactive compounds by steric and electrostatic contributions [39e41]. On the basis of the docked conformation, the molecules in the database were aligned to construct a CoMFA model. Compounds 3-8, 3-23 and 3-40 were excluded from the database because they were statistical outliers in the training set, that is, the inclusion of any of these molecules did not yield a satisfactory leave-one-out q 2 . The training set without the outliers gave a leave-one-out q 2 of 0.681 when the optimum components was 6, and the non-crossvalidated r 2 was 0.916, with a standard error of estimate of 0.088 and F values of 37.968. The steric and electrostatic contributions were 43.6% and 56.4%, respectively.

Compound 3-31 was used to illustrate the steric and electrostatic contour maps, together with the neighboring residues in the docked binding pocket (Fig. 5) . For the steric contour map, a bulky Table 1 The novel unsymmetrical aromatic disulfide compounds and their SARS-CoV M pro inhibitory activities (IC 50 group is favorable for better inhibition in the green contour region and such a group is likely to decrease the activity in the yellow contour space. The green maps are mostly in three bulks: one formed by Ser1and Leu141, one formed by Asn142 and Gly143 and the last one formed by Ser144 and Cys145; whereas the yellow maps are just located in a space nearby Gly143. For the electrostatic contour map, in the blue contour region, an increase in the positive charge will lead to an increase of activity, yet in the red contour region, negative charge is favorable to enhance the activity. The blue map is in a region surrounded by Leu27, Cly143, Ser144 and Cys145, whereas the red maps are in three cavities: one formed by Leu141 and Asn142, one formed by Ser144, Cys145 and Met165 and the last one formed by Asn142 and Gly143. The 3D CoMFA maps have afforded important structural features of the unsymmetrical aromatic disulfides from steric and electrostatic views, which is valuable for further design and discovery of more potent inhibitors.

The lack of effective anti-SARS agents makes it a possible danger when SARS breaks out sometime in the future, numerous people will be killed again. Therefore it is still an urgent demand to discover novel anti-SARS inhibitors to combat this deadly disease. SARS CoV M pro is an important target for the design of therapeutically useful drugs. In the present study, in an effort to develop non-peptidic anti-SARS inhibitors, we have identified for the first time, that some unsymmetrical aromatic disulfides are excellent inhibitors of SARS CoV M pro , the mechanism of which seems distinct as they are reversible and mpetitive inhibitors. This suggests that the unsymmetrical disulfides are promising lead compounds identification and development of a new family of biologically active anti-SARS agents. A possible binding model of the disulfide inhibitor was built by molecular docking, and a CoMFA model was constructed subsequently to point out the structural features of these novel inhibitors of SARS CoV M pro . Based on this information, further structural modifications are ongoing for better pharmaceutical compounds. The Lipinski rules will also be utilized to help to develop compounds with a final in vivo activity [42] . We are also trying to co-crystallize the protease and the best inhibitor, to gain insight into a real binding mode and explain the molecular basis of these compounds.

Various heterocyclic aromatic thiols 1 were commercial procured from 5A Pharmatech (China), Apichemical (China), Aldrich and ACES pharma, which were all >95% purity grade. All solvents and liquid reagents were dried in advance using standard methods and distilled before use. Substituted arenesulfenyl chlorides 2 were synthesized as described in our previous publications [21, 22] . Synthetic methods for compounds 3′ had been reported and these parent compounds for 3-4, 3-6. 3-7, 3-8, 3-9, 3-10 and 3-11 had also been fully characterized before [21] . Melting points were determined using an X-4 melting apparatus and were uncorrected. 1 H NMR spectra and 13 C NMR were obtained using a 400 MHz Varian Mercury Plus 400 spectrometer. The chemical shift values (d) for the NMR spectra were reported as parts per million (ppm), using deuterated chloroform (CDCl 3 ) or dimethyl sulfoxide (DMSO-d 6 ) as the solvent and tetramethylsilane (TMS) as an internal reference standard. Mass spectra were recorded on a Thermo Finnigan LCQ Advantage LC/mass detector instrument.

Heterocyclic aromatic thiols (1, 5 mmol) was added to a solution of freshly prepared arenesulfenyl chlorides (2, 5 mmol) in 25 mL of anhydrous ethyl ether at room temperature. The mixture was then stirred for 5 h at the same temperature, after that the solvent was removed under reduced pressure. Products 3 (3′ for 3-4, 3-6. 3-7, 3-8, 3-9, 3-10 and 3-11) were purified by column chromatography in 75e95% yields. Unsymmetrical aromatic disulfide (3′, 5 mmol) was added to 10 mL of acetic anhydride. The reaction mixture was stirred for 0.5 h at 60 C and then 200 mL water was added to the mixture. The pH value was adjusted by sodium bicarbonate to 7e8. The product was extracted by ethyl acetate and the solvent was removed under reduced pressure. The target compounds were finally purified by 

The expression and purification of SARS CoV M pro was described by Rao et al. [43] . Basically, the sequence of SARS-CoV M pro cloned into the pGEX-6P-1 vector was transformed into E. coli BL21 (DE3) cells. The GST fusion protein, GST-SARS-CoV M pro , was purified by GST-glutathione affinity chromatography, cleaved with PreScission protease, and the recombinant SARS-CoV M pro was further purified by using anion-exchange chromatography. Eventually purified protein was of high purity (>95%) as judged by SDSPAGE analysis and the concentration is 0.5 mM, and the buffer contains 50 mM Tris-HCl, pH 7.3 and 1 mM EDTA. The fluorogenic substrate with consensus sequence of CoV M pro , MCA-AVLQSGFR-Lys(Dnp)-Lys-NH2 (95% purity), was synthesized in Shanghai Biological Engineering Company. The substrate was dissolved in DMSO in 0.8 mM liquid storage for use.

The inhibition assay was similar to Yang's procedure [11] . The SARS CoV M pro inhibition assays were conducted by fluorescence resonance energy transfer (FRET). The natural substrate amino acid sequence (AVLQSGFRKK) of SARS-CoV M pro started with the MCA fluorescent group and connected the Dnp fluorescence quenching group with penultimate K. The tested compounds were dissolved by sterilized DMSO and diluted to various concentrations. The settled concentrations of proteins, compounds and substrate were preheated at 37 C and oscillated. The excitation/emission light was 320/405 nm, and the test was carried out every 5 s for 200 times. Drawing curves, the maximum value of the negative control curve slope is V 0 , and the largest compound curve slope is V1. The inhibition ratio can be defined (1ÀV 1 /V 0 ). And the IC 50 value was calculated by equation (1) using CraphPad Prism5: 50 (1)

V 0 shows the initial rate of the reaction without inhibitor, V means the initial rate of the reaction with the inhibitor at various concentrations and [I] indicates the concentration of the inhibitor.

The determination of the inhibitor as a valent inhibitor employs the above methods as well, albeit with two modifications. Firstly, the inhibitor concentration was set to 2 or 4 mM. And for each inhibitor concentration, we measured the enzymatic activity of M pro whose concentration spans 0e2 mM. Secondly, M pro and inhibitor were first incubated for 20 min to ensure a thorough 'M pro -inhibitor' reaction and then the inhibition assay was initiated by adding substrate and characterized by fluorescence monitoring.

The further characterization of the inhibitor as a noncompetitive inhibitor employs the methods described in earlier work [28] . Basically, the enzymatic velocity of SARS-CoV M pro versus substrate concentrations with presence of inhibitors is depicted by equation (2) [27] , where K i is the dissociation constant for the SARS-CoV M pro complexed with inhibitor 3-31; factor a reflects the effect of inhibitor 3-31 on the affinity of the enzyme for its substrate; V max and K m represent the maximum velocity and Michaelis-Menten constant, respectively.

The values of V max and K m at different inhibitor concentrations were apparent V max and K m , called V app max andK app m , respectively. According to equation (2), V app max andK app m can be calculated by equation (3).

The kinetic parameters of V Mass spectra were recorded on Waters Xevo G2-XS Q-TOF mass spectrometry. Mass spectra were acquired in positive ion mode using a capillary voltage of 3 kV, a sampling cone voltage of 40 V and a source offset voltage of 80 V. The cone gas flow was set up to 50 L/h and desolvation gas flow was 800 L/h. Desolvation temperature and source temperature were set to 400 C and 100 C, respectively. The mass of intact protein was obtained by deconvolution of the raw data using MaxEnt1 tool. The samples were prepared at the similar condition with the IC 50 determination, except that the concentration of the inhibitors was 20 times of the concentration of SARS M pro .

Chemical structures of the compounds were built within Sybyl 7.3 (Tripos Inc., St Louis, MO). All the molecules were assigned Gasteiger-Hückel charges and minimized by the Tripos force field when convergence reached 0.001 kcal/mol/Å. Molecular docking of the unsymmetrical aromatic disulfides to the active site of SARS-CoV M pro was performed by FlexX. The crystal structure of SARS-CoV M pro in complex with inhibitor (pdb code 2AMD) was retrieved from the pdb databank. All water molecules were removed, and hydrogen atoms were added in the standard geometry. Any amino acid residue within 6.5 Å of the location of the original inhibitor N9 was considered to be in the binding pocket. Cscore calculation was enabled and set to serial mode. Database docking and subsequent scoring procedures were performed using the default parameters in the program.

For CoMFA, The molecules were superimposed using 3-31 from the molecular docking result as the template. All the parameters were used the default value within CoMFA module and the column filtering was set to 2.0 kcal/mol. The "leave-one-out" (LOO) cross validation method was applied to determine the optimum number of partial least squares (PLS) components. The non-cross validated method was used to derive the final model to explain the quantitative structure-activity relationship in a three dimensional manner.

-nitrophenyl)disulfanyl)thiazole (3-24) Yield 92%; m.p.: 61e63 C

33 (s, 1H, CH), 2.33 (s, 3H, CH 3 ); 13 C NMR (101 MHz

-methylthiazol-2-yl)disulfanyl)benzoate (3-25) Yield 91%; m.p.: 103e105 C; yellow solid; 1 H NMR (400 MHz, DMSO-d 6 ), d 8.09e8.03 (m, 1H, ArH), 8.01 (d, J ¼ 7.9 Hz, 1H, ArH), 7.74 (t, J ¼ 7.7 Hz, 1H, ArH), 7.45 (t, J ¼ 7.5 Hz, 1H, ArH)

-methyl-1,3,4-oxadiazol-2-yl)disulfanyl) benzoate (3-26) Yield 93%; m.p.: 80e82 C

78 (s, 1H, ArH), 7.47 (d, J ¼ 6.3 Hz, 1H, ArH), 3.89 (s, 3H, OCH 3 ), 2.49 (s, 3H, CH 3 ); 13 C NMR (101 MHz, DMSO-d 6 )

-methyl-1,3,4-oxadiazol-2-yl)disulfanyl)benzoate (3-27) Yield 88%; m.p.: 73e74 C

50 (s, 3H,CH 3 ), 1.34 (t

-28) Yield 92%; m.p.: 91e93 C; yellow solid; 1 H NMR (400 MHz, DMSO-d 6 ), d 8.34 (d, J ¼ 8.2 Hz, 1H, ArH), 8.30 (d

4-oxadiazol-2-yl)disulfanyl)benzoate (3-29) Yield 87%; m.p.: 99e101 C

DMSO-d 6 ), d 9.37 (s, 1H, CH), 8.14 (d, J ¼ 8.2, Hz, 1H, ArH), 8.05 (d, J ¼ 7.8, Hz, 1H, ArH), 7.87 (t, J ¼ 7.8, Hz, 1H, ArH), 7.48 (t, J ¼ 7.7, Hz, 1H, ArH)

HRMS(MALDI) m/z

45 (t, J ¼ 7.2 Hz, 1H, ArH), 3.90 (s, 3H, OCH 3 ), 2.07 (s 3H, CH 3 ); 13 C NMR (101 MHz, DMSO-d 6 )

12 (d, J ¼ 25.8 Hz, 1H, ArH), 2.35 (s, 6H, CH 3 ); 13 C NMR (101 MHz, DMSO-d 6 )

HRMS(MALDI) m/z: calculated for C 12 H 11

Yield 84%; m.p.: 76e79 C; yellow solid; 1 H NMR (400 MHz, DMSO-d 6 ), d 7.53 (dd, J ¼ 15.6, 7.2 Hz, 4H, ArH), 7.15 (s, 1H, ArH), 2.38 (s, 6H, CH 3 ); 13 C NMR (101 MHz

Yield 81%; m.p.: 61e63 C; yellow solid

Yield 84%; m.p.: 72e75 C; yellow solid; 1 H NMR (400 MHz, DMSO-d 6 )

-nitrophenyl)disulfanyl)pyrimidine (3-37) Yield 86%; m.p.: 130e132 C

17 (s, 1H, ArH)

-chlorophenyl)disulfanyl)pyrimidine (3-38) Yield 80%; m.p.: 79e80 C; yellow solid; 1 H NMR (400 MHz, DMSO-d 6 ), d 8.66e8.47 (m, 2H), 7.47 (t, J ¼ 9.0 Hz, 2H), 7.29e7.16 (m, 2H)

-bromophenyl)disulfanyl)pyrimidine (3-39) Yield 83%; m.p.: 75e78 C

World Health Organization, Summary of Probable SARS Cases with Onset of Illness from 1

Substrate specificity profiling and identification of a new class of inhibitor for the major protease of the SARS coronavirus

Design of wide-spectrum inhibitors targeting coronavirus main proteases

Coronavirus main proteinase (3CLpro) structure: basis for design of anti-SARS drugs

Isatin compounds as non-covalent SARS coronavirus 3C-like protease inhibitors

Binding interaction of quercetin-3-b-galactoside and its synthetic derivatives with SARS-CoV 3CL pro : structure-activity relationship studies reveal salient pharmacophore features

Synthesis and evaluation of pyrazolone compounds as SARS-coronavirus 3C-like protease inhibitors

Design, synthesis, and evaluation of trifluoromethyl ketones as inhibitors of SARS-CoV 3CL protease

Discovery, synthesis, and structure-based optimization of a series of N-(tert-butyl)-2-(N-arylamido)-2-(pyridin-3-yl) acetamides (ML188) as potent non-covalent small molecule inhibitors of the severe acute respiratory syndrome coronavirus (SARS-CoV) 3CL protease

Development of potent dipeptide-type SARS-CoV 3CL protease inhibitors with novel P3 scaffolds: design, synthesis, biological evaluation, and docking studies

Synthesis, modification and docking studies of 5-sulfonyl isatin derivatives as SARS-CoV 3C-like protease inhibitors

Design, synthesis, and biological evaluation of novel dipeptide-type SARS-CoV 3CL protease inhibitors: structure-activity relationship study

Disulfide bond formation and its impact on the biological activity and stability of recombinant therapeutic proteins produced by Escherichia coli expression system

Hypochlorous acid scavenging activities of thioallyl compounds from garlic

Effects of selected organosulfur compounds on melanin formation

Dynamic molecular networks: from synthetic receptors to self-replicators

Unsymmetric aryl-alkyl disulfide growth inhibitors of methicillin-resistant Staphylococcus aureus and Bacillus anthracis

Selective inhibition of extracellular thioredoxin by asymmetric disulfides

Characterization of acetohydroxyacid synthase from Mycobacterium tuberculosis and the identification of its new inhibitor from the screening of a chemical library

Identification of the catalytic subunit of acetohydroxyacid synthase in Haemophilus influenzae and its potent inhibitors

Synthesis, crystal structure, in vitro acetohydroxyacid synthase inhibition, in vivo herbicidal activity, and 3D-QSAR of new asymmetric aryl disulfides

Synthesis and biological evaluation of nonsymmetrical aromatic disulfides as novel inhibitors of acetohydroxyacid synthase

Targeting of arenavirus RNA synthesis by a carboxamide-derivatized aromatic disulfide with virucidal activity

Antiviral propierties of 5, 5'-dithiobis-2-nitrobenzoic acid and bacitracin against T-tropic human immunodeficiency virus type 1

An antiviral disulfide compound blocks interaction between arenavirus Z protein and cellular promyelocytic leukemia protein

A traceless, one-pot preparation of unsymmetric disulfides from symmetric disulfides through a repeated process involving sulfenic acid and thiosulfinate intermediates

Reversible inhibitors

Crystal structure of feline infectious peritonitis virus main protease in complex with synergetic dual inhibitors

Dynamic combinatorial chemistry

Discovery of a novel covalent non-blactam inhibitor of the metallo-b-lactamase NDM-1

Design, synthesis, and evaluation of inhibitors for severe acute respiratory syndrome 3C-like protease based on phthalhydrazide ketones or heteroaromatic esters

A fast flexible docking method using an incremental construction algorithm

Chemical synthesis, in vitro acetohydroxyacid synthase (AHAS) inhibition,herbicidal activity, and computational studies of isatin derivatives

Comparative evaluation of several docking tools for docking small molecule ligands to DC-SIGN

A structure guided QSAR: a rapid and accurate technique to predict IC 50 : a case study

Evaluation of various inverse docking schemes in multiple targets identification

How to improve docking accuracy of AutoDock4.2: a case study using different electrostatic potentials

LIGPLOT: a program to generate schematic diagrams of protein-ligand interactions

Comparative molecular field analysis (CoMFA). 1. Effect of shape on binding of steroids to carrier proteins

Bioactive conformation analysis of anthranilic diamide insecticides: DFT-based potential energy surface scanning and 3D-QSAR investigations

N-Fluorinated phenyl-N 0 -pyrimidyl urea derivatives: synthesis, biological evaluation and 3D-QSAR study

Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings

The crystal structures of severe acute respiratory syndrome virus main protease and its complex with an inhibitor

This work was financially supported by the Natural Science Foundation of China (No. 21272128 and 21672114), the "111" Project of Ministry of Education of China (No. B06005) and the National Basic Research Program of China (No. 2013CB734004). We appreciate Professor Chuanzheng Zhou in Nankai University for his kind assistance for the mass spectrometry of SARS CoV main protease.

Supplementary data related to this article can be found at http:// dx.doi.org/10.1016/j.ejmech.2017.05.045.