Structural modification of antineoplastic drug carmofur designed to the inhibition of SARS-CoV-2 main protease: A theoretical investigation.

A coherent account of the reaction mechanistic details, structural modifications, and inhibition potentials of antineoplastic drug carmofur and its modified analogs to inhibition of SARS-CoV-2 main protease (Mpro) is reported. The survey is performed by integrating the density functional based tight binding (DFTB3) with density functional theory (DFT) calculations. The inhibition process commences with nucleophilic attack from the sulfur atom on the carbonyl group, yielding a C-S bond formation, followed by a bond formation of the H-O9 by 2.07 Å, which results in a transition state contains a ring of six atoms. We found that although the direct addition of sulfhydryl group hydrogen to the N3 position is likely to happen, the proper position of the hydrogen to O9 decreases its accessibility. The thermodynamic stability of the complex was calculated to be highly sensitive to the substituent on the N11 position. Compounds with CH2NH2 and CH2F at N11 positions of carmofur revealed high thermodynamic stability to complexation with Mpro but induced no change in substrate-binding pocket comparable to carmofur. Replacing the N11 of carmofur with carbon (C-carmofur) was effective in terms of complexation stability at CH2CH2CH2F and CH2CH2CH2OH substitutions and occupation of S1 subsite by these structures in addition to the S2 subsite. Based on the resulted data, increasing the length of the carbon chain at introduced substitutions in N-carmofur almost decreases the complexation stability while in C-carmofur the trend is reversed. Throughout these information outputs, it was suggested that compounds d, e, i', and k' might be novel and more efficacious drug candidates instead of carmofur. We believe that our characterization of mechanistic details and structural modification on Mpro/carmofur complex will significantly intensify researchers' understanding of this system, and consequently help them to take advantage of results into practice and design various valuable derivatives for inhibition of SARS-CoV-2 main protease.

Introduction

COVID-19 as a highly contagious pathogenic has rapidly spread worldwide in December 2019 from Wuhan city in Hubei Province of China [1]. While the number of reported confirmed cases worldwide continues to rise, efforts are underway to overcome this challenge. So far, no known effective drug is available to treat or alleviate the disease symptoms [2]. The development of new generation of antiviral drugs is urgently needed. Since, the process of new drug design is highly time‐consuming, risky, and costly, drug repurposing can be used as an alternative strategy, determining the new indications of existing drugs for another disease [3]. Reducing the risk of adverse side effects, drug interactions, and drug development time and expenditure are advantages of this method [4]. Considering the matter of time, the computational approaches for drug repurposing provide the best possible chance of selecting the most effective drug among the broad list of approved drugs for the life‐threatening emergency condition of COVID‐19 [5].

Carmofur, a derivative of 5-fluorouracil, is an antineoplastic drug containing an electrophilic carbonyl reactive group that inhibits human acid ceramidase (AC) through covalent modification of its catalytic cysteine [6]. Raised AC levels are linked to several malignancies, including breast [7], prostate [8], colorectal [9], melanoma [10], and brain [11]. The ability of carmofur to inhibit AC led to its anticancer activity. Since the 1980s, carmofur has been used to treat colorectal cancer [12]. It has also been shown clinical benefits against breast [13], gastric [14], and bladder cancers [15].

Similar to SARS and MERS, the coronavirus genome encodes four non-structural proteins including Spike (S‐protein), RNA-dependent RNA polymerase (RdRp), the main protease (MPro), and papain-like protease (PLpro), which recognized as promising targets for developing drugs and treatment against the recent coronavirus epidemic [16], [17], [18], [19]. A recent study shows that carmofur inhibits the SARS-CoV-2 main protease (Mpro) [20]. The X-ray crystal structure of the Mpro/carmofur complex shows that the carbonyl reactive group of carmofur is covalently bound to the Sγ atom of Cys145, while its fatty acid tail is inserted into the hydrophobic S2 subsite. Based on these observations, the authors suggested that the sulfhydryl group of Cys145 attacks the electrophilic carbonyl group of carmofur, resulting in covalent modification of Cys145 and release of the 5-FU moiety. Considering this, the present study is designed to elucidate the detailed molecular mechanism of inhibition of MPro by carmofur along with structural modification of the drug to inhibit SARS-CoV-2 main protease in a more effective affinity. We hope that the obtained results will help in the repurposing of already approved drugs to combat the recent dangerous coronavirus epidemic.

Model and computational details

The three-dimensional structure of the SARS-CoV-2 main protease (PDB ID: 7BUY) was downloaded from RCSB PDB database (http://www.rcsb.org/pdb). Before calculation, the protein was edited by cutting the Arg40 to Met49 residues, Leu141 to Cys145 residues, and His163 to Glu166 residues from the crystal structure. The starting configuration consisted of the crystal structure of the Mpro/carmofur complex, which has been fully optimized using the density functional based tight binding (DFTB3) method [21] by third-order parametrization for organic and biological systems (3OB)-3–1 parameter [22], [23], [24]. Vibrational frequencies were computed to obtain the Gibbs free energy for all the Mpro/carmofur analog complexes. The DFTB + program [25] was used for this part of the calculations. B3LYP hybrid functional (20% HF exchange) with Grimme’s DFT-D3 dispersion correction [26] were used for all geometry optimization related to transition state structures. The basis set was set to 6–311 + G(d,p) [27] for all atoms during the calculations. The Berny algorithm was used to transition states optimization. Vibrational frequencies are computed by using the same method/basis set as the geometries optimized. The first-order saddle point (transition state) was confirmed by only one imaginary frequency on computed vibrational frequencies. Intrinsic reaction coordinate (IRC) calculation was computed to confirmed transition structure connection to minima [28], [29]. SMD/Water variation was considered to consider the solvent effects on all the calculations. All the optimization and frequency calculations in this part of the study were executed with the ORCA program package [30]. In the present study, the free energy for each compound in solution was computed through the following formula:where ΔG1atm→1M = 1.89 kcal/mol is the change of free-energy for the contraction of 1 mol of an ideal gas from 1 atm to the 1 M solution-phase standard state. The atom numbering used throughout the study is illustrated in Fig. 1 .

Fig. 1

The molecular structure of carmofur with atom numbering.

Results and discussion

The main purpose of the present study is to modify the carmofur molecular structure to develop novel carmofur analogs that inhibit SARS-CoV-2 main protease and therefore cause the virus to malfunction. The crystal structure of Mpro in complex with carmofur shows that the hydrophobic sequence tail of carmofur is inserted into the S2 subsite, while the carbonyl moiety (CO) is linked to the Sγ atom of Cys145 through a 1.8 Å covalent bond [20]. Moreover, a hydrogen bond between the O12 of the inhibitor and the NH2 hydrogen of Gly143 further stabilized the complex. The carmofur only occupies the S2 subsite of the Mpro [20]. Considering that the carbonyl group of carmofur is easy to be attacked by sulfhydryl group [31], we investigated the mechanism of Mpro/carmofur complexation according to Scheme 1 .

Scheme 1

The different reaction pathways for inhibition of Mpro by carmofur.

The reaction pathway leading to products Mpro/carmofur complex and 5-fluoro-2-hydroxypyrimidin-4(3H)-one, overcomes an energy barrier of as high as 45.24 kcal/mol through the transition state shown in Fig. 2 . The transition state contains a ring of six atoms, in which the sulfhydryl group acts as an entering group while 5-fluoro-2-hydroxypyrimidin-4(3H)-one leaves the carbonyl group. Nucleophilic attack from the sulfur atom on the carbonyl group leads to the bond formation of C10-S by 2.24 Å. Moreover, the process involves the changing of the C10-N3 bond from 1.47 to 2.04 Å and a bond formation of the SH-O9 by 2.07 Å. With the formation of the SH-O9 bond, the C10-N3 bond length increases from 1.47 to 2.04 Å, and hence the C4-N3 bond length decreases from 1.38 to 1.35 Å. In addition to the mentioned bond formations, the transition state is stabilized by two hydrogen bonds of H19-O9 (2.10 Å) and the hydrogen of the NH2 group at Cys145 and O9 (2.34 Å). This indicates that the formation of 5-fluoropyrimidine-2,4(1H,3H)-dione is unlikely through the direct addition of sulfhydryl group hydrogen to the N3 position as proposed by Jin, Z. et al [20], while it can be formed in the next steps by a tautomerization reaction with an energy barrier of as high as 9.08 kcal/mol [32] (see Scheme 1).

Fig. 2

Key parameters for the transition state structure located for Mpro/carmofur complexation using the B3LYP-D3BJ/6–311 + G(d,p) level of calculations. The interatomic distances are given in Å.

To apply the structural modification found in this study to future drug discovery research, especially for drug candidates containing carboxamides, this work aimed to stabilized the Mpro/inhibitor complexation along with the inhibitor structural flexibility to block different sites on Mpro. Because the CONH2 moiety in Asn142 can increase the stability of the complex by creating a hydrogen bond (see Fig. 3 ), we first modified carmofur by introducing substitutions at the N11 position of this inhibitor (see Fig. 4 , A).

Fig. 3

The crystal structure of Mpro/carmofur complex labeled Cys145 and Asn142.

Fig. 4

Modification of carmofur analogs by introducing substitutions at the N11 position (N-Carmofur-(CH2)nX) and or by altering the N11 to carbon and introducing substitutions at the new position of C11 (C-Carmofur-(CH2)nX).

Although increasing the length of the substitutions brings the hydrogen donor closer to the oxygen lone pair electrons of the CONH2 moiety, but according to results in Table 1 , the most complexation stability is seen at n = 1. While introducing NH2 and OH into the N-Carmofur (g and k) exhibited a modest increase in relative Gibbs’s free energy of the complex, introducing NH2 and F into the N11 position at n = 1 (N-Carmofur-(CH2)2NH2 and N-Carmofur-(CH2)3OH) made the complex highly thermodynamically stable (see d and e in Table 1). The d and e structures with −12.73 and −10.12 kcal/mol respectively are both stable than the Mpro/carmofur complex. The d inhibitor is stabilized by both a 2.04 Å hydrogen bond between the NH2 hydrogen of introduced substitution and the oxygen lone pair electrons of CONH2 moiety in Asn142, and also hydrophobic interactions (see Fig. 5 ). In the same manner, e structure also forms a hydrogen bond of 2.30 Å between the CH2F hydrogen of introduced substitution and the oxygen lone pair electrons of CONH2 moiety. Hydrogen bonding in e structure can occur due to the high electronegativity of the fluoride and so leaving the carbon relatively electron-poor. These data indicate that compounds d and e might inhibit Mpro by a relatively high negative binding energy. Besides, the results show that H at n = 1, F at n = 2, and H, NH2, and F at n = 3 are not preferable for Mpro inhibition (see b, I, j, l, and m in Table 1).

Table 1

The relative free energies for carmofur (black row) and its analogs in complex with Mpro. The relative free energies are given in kcal/mol.

n/X	N-Carmofur-(CH2)nX		C-Carmofur-(CH2)nX
n = 0, X = H	a	−7.17	a′	−2.60
n = 1, X = H	b	−6.12	b′	−6.40
n = 1, X = OH	c	−7.01	c′	−2.28
n = 1, X = NH2	d	−12.73	d′	−5.75
n = 1, X = F	e	−10.12	e′	−3.63
n = 2, X = H	f	−9.46	f′	−7.50
n = 2, X = OH	g	−7.58	g′	−7.51
n = 2, X = NH2	h	−8.53	h′	−5.64
n = 2, X = F	i	−4.93	i′	−13.15
n = 3, X = H	j	−5.12	j′	−5.98
n = 3, X = OH	k	−8.57	k′	−11.77
n = 3, X = NH2	l	−5.60	l′	−5.99
n = 3, X = F	m	−5.78	m′	−4.64

Fig. 5

Overall structure of Mpro in complex with d, e, i′, and k′ structures.

Next, we modified carmofur by altering the N13 to carbon and introducing substitutions at the new position of C11 (see Fig. 4, B). Replacing nitrogen with carbon without any other substitution modifications drastically reduces the Mpro/carmofur complexation stability (see a′ at Table 1). Also, unlike the previous case, all of the n = 1 compounds are highly unstable than the Mpro/carmofur complex, a. Nevertheless, compounds i′ and k′ are more thermodynamically stable than a and d. The i′ analog form a hydrogen bond of 2.34 Å between the CH2F hydrogen of its modified substitution and the oxygen lone pair electrons of CONH2 moiety. In addition to this hydrogen bond, the i′ inhibitor is also stabilized by hydrophobic interactions of modified substitution moiety (CH2CH2F) with Met165 of S1 subsite (see Fig. 5). Altogether, this structure with complexation energy of −13.15 kcal/mol is −5.98 kcal/mol more stable than a. On the other hand k′ by form a strong hydrogen bond of 1.95 Å between the OH hydrogen of its modified substitution and the oxygen lone pair electrons of CONH2 moiety and also hydrophobic interactions with Met165, His163 and His164 of S1 subsite is −4.60 kcal/mol more stable than a. Therefore, our results show that the mechanism of C-carmofur modifications is different from that of the N-carmofur compounds. In addition to the greater thermodynamic stability of i′ and k′ than N-Carmofur modifications, the largest conformational differences occur in the substrate-binding pocket, d and e analogs only occupy the S2 subsite, whereas C-Carmofur-(CH2)2F and C-Carmofur-(CH2)3OH occupy S1 subsite in addition to the S2 subsite. These findings indicate the structural elaboration potential of carmofur and will be helpful for the design of more potent derivatives against the Mpro.

Conclusions

In addition to the characterization of mechanistic details and specify the exact geometry of transition state involved in Mpro/carmofur complexation, we suggested several carmofur analogs and evaluated their complexation stability and intermolecular interactions with M pro. Compounds d, e, i′, and k′ were highly stable to Mpro and their Gibs free energy was much stronger than that of carmofur. Although these four compounds interact relatively identically with Mpro, compounds i′ and k′ occupy S1 subsite in addition to the S2 subsite. We propose that all d, e, i′, and k′ compounds might serve as drug candidates that could be developed for use as efficacious Mpro inhibitors drugs with much stronger binding energy than carmofur.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

CRediT authorship contribution statement

Niloofar Hemati: Project administration. Saba Hadidi: Conceptualization, Formal analysis, Investigation, Resources, Software, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing. Farshad Shiri: Formal analysis, Investigation, Methodology. Mohammad Hosein Farzaei: Supervision, Project administration.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.