1. Introduction
Lactoferrin (Lf) is a member of the transferrin, that binds and transfers iron in the blood and it can also remove excess iron from the body, which is a most abundant glycoprotein in human and ruminant milk resources [
1]. The previous studies of lactoferrin disclosed a wide display with different functions including anticancer activity, anti-inflammatory, and cognitive function enhancement in patients with Alzheimer’s. Lactoferrin has 1–4 glycans with single-chain polypeptides of about 80,000 Da. Lactoferrin is present in large scales in milk with multifunctional glycoprotein including N-glycans, which are active with functional groups and depending on the species, and it makes a significant contribution to the host that defines the system. [
2,
3,
4]. In addition, lactoferrin carries many important biological functions, including N-glycans bonding to iron or others, being bioactive in cell explosion and diversity, as an anti-parasitic protein. Lactoferrin is known to generate host protective responses against Mycobacterium tuberculosis, anti-bacterial and anti-viral. These functions diverge from lactoferrin considerable attention as the primary nutritional contribution to iron-binding by the role of glycosylation[
5,
6]. The molecular structure of human Lactoferin and amino acids were discovered in two globular lobes of ~700 amino acids stabilized by disulfide bonds, which are linked by a flexible alpha helix- amino and carboxy, N-lobe and C-lobe [
7].
Ionic liquids containing N-active group, that’s make the critical role of carbonyl groups of glycan binding in many biological processes very easily. The chain of saccharides-glycan moieties in lactoferrin is likely to contribute significantly to the N-ionic liquids roles by carbonyl of saccharides. Despite the high amino group of ionic liquid sequence homology in different with excellent results, which exhibits a unique N-glycosylation for heterogeneity of the biological properties and Lactoferrin is chosen as a good example source of N-glycans [
6,
8,
9].
Exploring and identifying the new characterization of novel ionic liquid delivers reacting with oligosaccharides of glycoproteins. Several interesting studies are encouraging for discovers the new application of ionic liquids as potential drugs antimicrobial including anti- coronavirus disease 2019 [
10]. Studying the activity of glycoproteins is the assessment of the contributions of individual glycans to the observed bioactivities. This work examines how the study of N-link glycosylation in Lactoferrin which reacted with ionic liquid MIE-NH
2 increases the understanding of ionic liquid functionality [
2,
11,
12]. Since 2019, researchers have revealed various structural and non-structural targets of SARS-CoV-2 that have been utilized in drug repurposing in the fight against COVID-19. In this study, some of such targets that have been investigated well were selected for the computational study. There are clinically available drugs that target main protease (Mpro), RNA-dependent RNA polymerase (RdRp), Transmembrane protease serine 2 (TMPRSS2), and Papain-like protease (PLpro) [
13,
14,
15]. Drugs like favipiravir and remdesivir, which have been used for the treatment of COVID-19, are examples for RdRp inhibitors. Lopinavir and ritonavir are PLpro and M inhibitors whereas nafamostat and camostat act on TMPRSS2 [
16]. These drugs (IL-Lfs) have small molecular structures with various scaffolds. The probability of the synthesized molecule to inhibit any of these targets was investigated to find out its potency as a SARS-CoV-2 inhibitor.
3. Results
The glycoproteins in bovine lactoferrin (BLf) were chosen as good substrates source to prepare N-glycans with high structure including oligosaccharides [
28]. The high proportion of glycosylation verifies the methodology of MIE-NH
2 following reductive amination, it was used for labelling of N-glycans [
17]. According to our previous studies, ionic liquids drive lactoferrin iron free (IL-Lf-iron free)[
18]. There was found proximity 42 types of N-glycans with diverse potential sites N-glycosylation in bovine Lactoferin (BLf), the different N-glycans with all structures sites [
29,
30]. The lactoferrin IL-Lf molecules modified by using ionic liquid, there are 14 different MIE-NH
2 derivative lactoferrin-N-glycans were deduced according to UPLC profile and MS spectrum see in
Figure 1. The corresponding structures of Lactoferrin-MIE-NH
2 were assigned as shown in
Figure 1 and
Figure 2. The results of the detection in figure 1 suggested the possible structures of compounds were modified by IL-MIE-NH
2 and this result was confirmed by MALDI ToF analysis. The m/z values of structures either with mono-charge or di-charge were calculated related to the signals of MIE-NH
2 linked to N-glycans was observed. In the extracted ion chromatogram of the products of N-glycans linked to MIE-NH
2 from Lactoferrin by HPLC two peaks exhibited the same m/z value of new products 1716.50 which assigned and identified with theoretical m/z = 1716.70 [m]
+). In this case, we suggested the new product is MIE-NH
2 linked monofucosylated monogalactosylated bi antennary complex N-glycan isoforms [
17,
31]. For example, from LC-MS analysis, it was found the peak of 13.7 min was assigned as MIE-NH
2-FA2G1 and the peak at 14.6 min was derived as MIE-NH
2-A2G1F. This work demonstrated the catalytic mechanism of the derivatization of Lactoferrin-N-glycans with ionic liquid MIE-NH
2 following the reductive amination. The free aldehyde realized in the acidic medium and by reducing ligand as sodium cyanoborohydride, which possesses significant converted the carbonyl to an imine by the NH
2 group of MIE-NH
2 (
Scheme 1,b)[
18].
Recent scientific research focused on the progress in protein-based nanomedicine, albumin-paclitaxel as nanoparticles have been introduced in novel therapeutics and used for the treatment of cancer and viral infections [
10,
29]. However, specific drug targeting of SARS-CoV-2 is almost challenging and absent until now, premature drug release and supports the poor pharmaceutical stores for resistance COVID-19 and its mutations. Therefore, some studies with alternative protein-based nanomedicines have opening the eyes to the use ionic liquids for extend and developing a novel of small molecules form glycoproteins. Regarding to this challenge, lactoferrin (Lf-iron free) offers a promising bioactive well as potentials therapeutic and drug nano carrier. In this work, we focused on the major pharmacological actions of modified glycoproteins form lactoferrin with ionic liquids to produce new molecules including antiviral, anti-cancer, and/ or improve immunology.
To enhance the efficacy of glycoproteins as potential drugs anti SARS-CoV-2 it was functionalization of N-glycans with an emphasis of lactoferrin. Besides this technique wide application, it’s depended on the recent advances of ionic liquids-Lf-based small molecules as efficient platforms for delivering novel drugs anti-viral drugs, particularly for treating the COVID-19 infections.
The binding potential of the new IL-Lf molecule to four targets (Mpro, RdRp, TMPRSS2, PLpro) was investigated via molecular docking. The resulting interaction was compared to a standard drug, remdesivir, which is one of the drugs that has been used in the fight against COVID-19 and approved by the FDA for this indication. The crystal structure of RdRp utilized in this study has remdesivir inside it. A previous study suggested that remdesivir exhibited its activity by binding to RdRp [
23]. As a result, we investigated the binding potential of remdesivir to RdRp first in order to validate the docking process. The docking investigation demonstrated that remdesivir interacted to the enzyme with ten conventional hydrogen bonds Urd7 (2), Ade8, Gua10, Cyt11, Cyt12, Gua13, Asn496, Asn497 (2)) and one other interaction (Ade10). The compound interacted with the enzyme very well with -8.9 kcal/mol binding energy. The ligand had interactions mostly on the nucleotide residues (
Table 1,
Figure 3). Similarly, in the previous crystallographic analysis, remdesivir had interactions with the nucleotides to exert remdesivir-induced RdRp stalling [
23]. The high level of binding observed for remdesivir leads us to assume that the docking protocol would give reliable interaction of the compounds with the targets. Thereafter, molecular docking of the modified molecule and remdesivir to the targets was pursued.
The derivatized molecule had a good level of interaction with the main protease. It interacted with the enzyme through nine conventional hydrogen bonds (Thr26, Asn142, Ala145, His163, His164, Glu166, Val186, Arg188 (2)) and one other interaction types (Thr25, Met49, Gln189) (
Figure 3, Table1). It interacted with the enzyme stronger than remdesivir as it formed four more conventional hydrogen bonds. The binding energy of remdesivir was slightly lower than that of the molecule. Therefore, the two compounds are expected to have a similar affinity to the enzyme with a slightly higher affinity for remdesivir. A previous crystallographic study revealed that ligands had interactions with the enzyme at His41, Ala145, His163, His164, Met165, Glu166, and Gln189 [
22]. In the computational study, the interactions at all of these residues were observed (
Figure 3,
Table 1). In this regard, the computational study gave similar interaction points with the experimental study. The modified compound interacted with RdRp very well. It formed ten conventional hydrogen bonds (Urd7, Ade8 (2), Cyt9(2), Cyt12, Gua13, Asn497(3)). Together with this, its interaction was weaker than that of remdesivir as the later had one more carbon hydrogen bond interaction (
Figure 3,
Table 1). In addition to this, remdesivir had lower binding energy than the molecule that implicated a better affinity for it. Therefore, remdesivir interacted with RdRp stronger and had a higher affinity towards the enzyme. The interactions detected were mostly with the nucleotide residues of the enzyme for both of them as observed in a previous study[
23].
The derivatized IL-Lf molecule had strong interaction with TMPRSS2 with seven conventional hydrogen bonds (Glu299, Lys300, Gly439, Ser441, Ser460, Gly462, Gly464) and three other interactions (His296(2), Gln438) (
Figure 3,
Table 1). The molecule had better interaction with TMPRSS2 in relative to remdesivir as it formed more hydrogen bonds. However, the binding affinity of remdesivir is expected to be slightly higher than the molecule as it had lower binding energy. An experimental study reported that TMPRSS2 had interactions with ligands at Ile256, His296, Asp345, Asp435, Ser436, Gly439, and Ser441 residues [
24]. In this computational study, the interactions at His296, Ser436, Gly439, and Ser441 were observed. There was the similarity between the experimental and computational studies as more than half of the residues were common for both methods. In addition to this, the computational study gave strong interaction to TMPRSS2 for the two compounds. The molecule interacted with PLpro through four conventional hydrogen bonds (Leu162, Glu167, Gln269 (2)) and three other interactions (Glu161, Glu167, Gln269) (
Figure 3,
Table 1). The interaction was good but weaker than the interaction of remdesivir with the enzyme. The binding affinity of remdesivir was also better as it had lower binding energy [
33]. Therefore, remdesivir is expected to have better interaction with PLpro. Furthermore, the synthesized molecule had the weakest interaction with this enzyme. A previous crystallographic study reported that various ligands had interactions with Asp164, Arg166, Glu167, Tyr264, Tyr268, and Gln269 [
23,
34]. In this study, the synthesized small-molecule ionic liquid-Lf had interactions with Glu167 and Gln269 residues. The molecule had some level of interaction similarity with the previous study. On the other hand, remdesivir didn’t have any common interaction residue with the experimental study. This has implicated that it could interact with the enzyme but at a different binding region.
The computational study revealed that the synthesized molecule interacted with the targets very well generally. The strength of the interaction was different from each other with a decreasing order of RdRp, Mpro, TMPRSS2, and PLpro. The molecule had a stronger interaction than remdesivir with Mpro and TMPRSS2. Together with this, remdesivir had slightly lower binding energy than the molecule with the four target structures. Overall, the computational study demonstrated that the designed IL-Lf molecule had high binding potential toward the targets, especially towards Mpro and TMPRSS2.