key: cord-0970588-9ww58w3k authors: Bibi, Shamsa; Urrehman, Shafiq; Khalid, Laryeb; Yaseen, Muhammad; Khan, Abdul Quyyam; Jia, Ran title: Metal doped fullerene complexes as promising drug delivery materials against COVID-19 date: 2021-08-09 journal: Chem Zvesti DOI: 10.1007/s11696-021-01815-4 sha: d28553f80c02ea12c7137d3ac2de5f165642a5b6 doc_id: 970588 cord_uid: 9ww58w3k An outbreak of respiratory disorder caused by coronavirus has been named as coronavirus infection 2019 (COVID-19). To find a specific treatment against this disease researchers are at the frontline. To cure COVID-19, favipiravir (FPV) has been reported as an effective drug based on its high recovery rate. Among nanomaterials, fullerene C60 has achieved enormous attention as a drug delivery vehicle due to its good bioavailability and low toxicity. Hence, in this work, we have investigated the potential of metal-doped fullerene as a drug carrier, based on DFT calculations by using M06-2X functional and 6-31G(d) basis set in water media. In this research electronic parameters and adsorption energy of FPV on interaction with metal-doped (Cr, Fe, and Ni) fullerene is studied. The charge transfer between drug and doped fullerene has been studied through electrophilicity indexes. The structural and electronic properties are explored in terms of adsorption energy through frontier molecular orbital (FMO) and density of state (DOS). It is observed that doping of fullerene C60 with Cr, Fe, and Ni metals significantly enhances the drug delivery rate and provides numerous advantages including controlled drug release at specific target sites which minimize the generic collection in vivo and reduce the side effects. Thusly, it is suggested that our designed metal-doped complexes might be efficient candidates as drug delivery materials for COVID-19 infection. [Image: see text] SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11696-021-01815-4. In December 2019, the first case of coronavirus known as coronavirus infection 2019 (COVID-19) similar to severe acute respiratory syndrome (SARS) was reported in China and spread across the globe within two months (Zu et al. 2020) . It is a respiratory disorder and has become the major pathogen of emergency respiratory outbreaks (Harapan et al. 2020 ). According to the latest studies, the most common symptoms of coronavirus illness are fever, dry cough, fatigue, myalgia, and the less common symptoms are abdominal pain, headache, nausea, vomiting, and diarrhoea (Song et al. 2020) . Often infected patients of coronavirus lead to multiple organ failure due to severe respiratory disorder that results in a high mortality rate . The SARS was started about 18 years ago, genetically indistinguishable and having 79.6% uniqueness with COVID-19. This similarity can be advantageous for theoretical investigation (Hiscott et al. 2020) . However, to control the spread of this disease and to discover an effective antiviral drug, the healthcare workforces are in the research battle zone (Huang et al. 2020) . For the treatment of coronavirus different antiviral drugs, including Remdesivir, Favipiravir, Arbidol, and Chloroquine are under investigation (Dong et al. 2020; Wu et al. 2020) . Comparison between Remdesivir and Favipiravir (FPV) suggests that FPV can be most effective to cure coronavirus based on its high recovery rate and reduce the pyrexia and cough (Sreekanth Reddy and Lai 2021) . FPV is an RNA-dependent RNA polymerase (RdRp) inhibitor that is effective for the treatment of the Ebola virus and influenza (Sissoko et al. 2016) . Recently it has been highlighted that FPV as a prodrug was efficacious in reducing coronavirus infection (Udwadia et al. 2021) . So clinical researchers are required to check the effectiveness and safety of this antiviral drug against coronavirus. In pharmaceuticals, nanostructure-based drug delivery system has become most efficient to reduce side effects and improve the bioavailability of the drug. In comparison with conventional dosage, these drug delivery systems provide numerous advantages including controlled drug release at the specific effective site, minimize generic collection in vivo tissues, enhance cellular uptake, improve adsorption and reduce the harmful side effect of drugs (Slepička et al. 2013; . Various nanostructures have been established for drug delivery applications which are superior to microstructures because of the high surface to volume ratio. Among nanomaterial carbon-based nanomaterials including graphene, fullerene, and carbon nanotubes (CNT) have been widely considered because of their unique physical and chemical properties in nanotechnology as drug delivery vehicles (Raphey et al. 2019) . A drug delivery system is considered a major therapy, as it employs nanomaterials to deliver a specific amount of drug at the targeted site in a controlled way (Hazrati and Hadipour 2016a, b) . Nanostructures help to load and conjugate the drug molecule to the targeted cells so can be used as an ideal candidate for drug delivery. Among nanomaterials, fullerene and its derivatives have achieved enormous attention to use as drug delivery purposes due to their unique properties such as hollow cagelike structure, versatile physical, chemical, and biological properties, effective drug loading capability, and antioxidant capacity (Hazrati and Hadipour 2016a, b; Ergürhan et al. 2018 ). On the other hand, also having fewer biological side effects due to its hydrophilicity. As carbon is hydrophobic, so the major limitation of carbon-based nanomaterials for biological application is their hydrophobicity (Kazemzadeh and Mozafari 2019). To overcome this problem various methods such as surface functionalization with specific groups, encapsulation, and chemical modification through doping by using specific atoms are commonly used to show promising results and to enhance the hydrophilicity of nanostructures (Kumar and Raza 2017; Henna et al. 2020; Maleki et al. 2020) . It is also investigated that through the doping phenomenon adsorption properties and sensitivity can be improved. Magnetic and electronic properties can also be improved through doping and functionalization to get the desired results (Parlak et al. 2020) . Different theoretical and experimental studies have been done to investigate the interaction between different drugs and fullerene (Bashiri et al. 2017; Samanta and Das 2017; Parlak et al. 2020) . It has been concluded that doping of fullerene with metal is suitable to improve the drug delivery system through enhancing adsorption potential (Hazrati and Hadipour 2016a, b) . It also has been investigated from previous work that Ni and Cr doped fullerene significantly improves the adsorption of uracil, thymine and adenine (Rad and Aghaei 2018) . In this work, we aimed to check the ability of fullerene as a favourable candidate of drug delivery to adsorb FPV drugs through density functional theory (DFT) calculations. Recently a theoretical study was carried out to investigate adsorption between doped and undoped fullerene and FPV. In which nine different possible interactions were checked. It was observed that -N,-OH, and -OH interaction edges have the highest adsorption energy in the water phase . In another study, it was also investigated the adsorption between transition metal-doped C20 and FPV (Rad et al. 2020) . Notably, there is no work on the interaction between -N, -OH interaction edges of FPV and chromium (Cr), iron (Fe), and nickel (Ni) doped C60. Therefore, in the following research, we are going to design FPV…CrC59, FPV… FeC59, and FPV…NiC59 complexes, to investigate adsorption energy between FPV and iron (Fe), chromium (Cr), and nickel (Ni) doped C60 for the first time (Scheme 1). The experimental investigation on new material is also essential but they are time-consuming and costly and do not give enough information. Therefore, computational work has been employed to understand the nature of interaction and its mechanism (Hazrati and Hadipour 2016a, b; Novir and Aram 2020) . Some of the important chemical properties such as binding energy (E b ), energy gap (E g ), electrophilicity indexes (ω), chemical hardness (η), the electronegativity of a system (χ), and maximum electronic charge (ΔN max ), density of states (DOS) and molecular electrostatic potential maps (MEPS) have been calculated in this work. All calculations were carried out with Gaussian 09 program (Shi et al. 2010) . Initially, all the structures of complexes were built through the GaussView program 6.0 (Dikmen et al. 2018) . To get the most stable configuration, optimization of all structures was computed by density functional theory (DFT) method with M06-2X functional and 6-31G(d) basis set (Patel and Singh, 2018) . PyMOlyze 2.0 program was used to calculate the density of state (DOS). To define the relative stability regarding the nanostructures-drug interactions of complexes adsorption or binding energy has been performed through the following equation. In Eq. (1), E complex is the total energy of doped fullerene attached with drug, E C59, and E FPV is individual energies of optimized structures. The counterpoise (cp) correction method is used to eliminate the basis set superposition error (BSSE) (Boys and Bernardi 2006) . The effect of FPV adsorption on electronic parameters has been investigated through the energy gap (E g ). E g is the difference between the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO). To measure stabilization energy, the electrophilicity index (ω) is one of the focal parameters when a system gets an additional electronic charge. Electrophilicity index is a charge transfer descriptor and calculated by the following equation (Korotkikh et al. 2019 ). (1) Scheme 1 Scheme of designed metal-doped complexes, FPV…CrC59, FPV…FeC59, and FPV…NiC59 In the above equations, μ and Ƞ represent the chemical potential and chemical hardness, respectively. The (ΔN max ), maximum electronic charge accepted from the surrounding have been calculated through the following formula. where χ is the electronegativity of the system and ω is the electrophilic index. The amount of charge transfer between drug and doped fullerene can be defined through electrophilicity. In the above equation χ A is electronegativity of drug molecule and χ B is electronegativity of doped fullerene, if ECT < 0 drug will act as donor, and if ECT > 0 drug will act as acceptor. To calculate the electronic sensitivity and to correlate it with electrical conductivity and energy gap following expression is used; In the above equation, k represents the Boltsman constant and an exponential relationship present between electrical conductivity and Eg. Molecular electrostatic potential (ESP) has also been calculated by using the same basis set to investigate the nature of interaction and to study the charge distribution. Due to the absence of charge distribution, C60 exhibits zero dipole moment. Through doping one of the carbon atoms from the C60 cage-like structure has been replaced by Cr, Fe, and Ni, and the final complexes have been optimized with the same M06-2X/6-31G(d) method in water media. The optimized structures of FPV…Cr59, FPV…FeC59, and FPV…NiC59 is presented in Fig. 1 . According to the atomic charges, active sites of FeC59, CrC59, and NiC59 are Fe, Cr, and Ni, respectively, due to the distribution of positive charge around the metal atom. All the dopant atoms cause deformation in fullerene at the doped points. For numerical evaluation distance between different atoms of drug and metal-doped fullerene bond lengths is reported in Table 1 . Bond angles and dihedral angles also have been calculated to study structural changes. Table 1 indicates that doping of Fe, Cr, and Ni increases the distance of dopant atoms with other atoms, due to which deformation in a structure occurs and this deformation is more prominent for large size dopant atoms. The calculated bond lengths suggest that the interaction of FPV drug with metal-doped fullerene is considerable. It can be visualized that for each metal atom FPV had a unique configuration and orientation due to the hybridization of different orbital between them. Overall, it proves planner confirmation for all doped metal complexes. It can be visualized from the results given in Table1 that slight variation in bond length and dihedral angle occur after the substitution of one of the carbon atoms with a metal atom. Among metal atoms, the Ni atom is more closely bonded with surrounding carbon atoms with a bond length of 1.869, 1.868, and 1.833. Also, the drug molecule is most strongly bound onto NiC59 with a bond length of 2.071. Hence, the highest adsorption energy was observed for the Ni atom making it the most stable adsorbent for FPV drug than others. Furthermore, the torsion angle expresses the rotatable bond conformation around their axis of rotation. The negative values of the dihedral angle express the anticlockwise, while positive values express the clockwise rotations. This study shows the electron flow from the drug (donor) to metal-doped fullerenes (acceptor) with less steric hindrance and great convenience. However, these results verified that metal doping may affect the geometrical parameters, but any drastic effect to distort the geometry and stability of the compound was not observed. Through doping of electropositive metals, charge distribution occurs which results in increased dipole moment, electrophilicity, and more charge transfer from drug to fullerene which in turn increases the binding energy and decreases the E g . After doping maximum charges are transferred to the doped electropositive metal atoms. So, these metal atoms are the active site for the binding energy of the drug. The negative values of binding energy indicate that adsorption of FPV with doped fullerene is an exothermic and spontaneous process (Table 2) . Among all possible interaction edges of FPV to interact with doped fullerene, the most favourable edge is -N, -OH. On the other side, the dipole moment value increases after doping of metal, and charge density around the metal is maximized. Therefore, for fullerene, the metal atom is the most favourable site for energetic interaction. By comparison Cr, Fe and Ni, it can be visualized that, Nidoped fullerene showed the highest binding energy towards drug, while the binding energy of Cr and Fe doped fullerenes are also considerable. These resulting binding energy values reveal that from moving left to right in the first row of transition metals E g decreases and binding energy increases. Interestingly, it is also observed that a decrease in E g and increase in binding energy for metal-doped C60 fullerene is more than that of C20 (Rad et al. 2020 ). Among all investigating compounds, FPV…NiC59 complex exhibits the highest binding energy and the strongest interaction in the water phase. Hence, results prove that the interaction of a drug with doped fullerene influencing the magnetic and electronic properties and this change is more prominent for Ni atom. Thus, doping of electropositive metals results in increased dipole moment which in turn enhances charge transfer from drug to fullerene with an increase in binding energy. The molecular structure of FPV and doped fullerene complexes were fully optimized in the ground state. Figure 2 displays that HOMO and LUMO energy levels were different for all complexes due to which changes in the electronic structure of doped metallofullerene were also different. It can be seen that HOMOs shifted to lower negative and LUMOs shifted to higher negative after the adsorption of the drug to the doped fullerene (Table S2 ). Biological activities also increase with an increase in E HOMO or by a decrease in E LUMO . Therefore, changes in E g occur due to changes in HOMO and LUMO energy levels. The irregular trend of HOMO and LUMO energy levels is due to doping of metals and this effect varies from metal to metal. The density distribution of HOMO and LUMO orbitals was also affected after the adsorption of the FPV drug on metallofullerenes (Fig. 3) . For LUMO density distribution, the effect of drug adsorption was far obvious. In the case of FPV…CrC59, LUMO orbital was almost shifted to the drug, while in the case of Fe and Ni complex this shift was also noticeable. This shift in the distribution of orbitals density occurs due to the transfer of charge from drug to metal atom. Graphical presentation of HOMO and LUMO along with E g is provided in Fig. 2 , Table S2 . As HOMO energy levels of all the complexes increases and LUMO energy levels decrease, E g values for all complexes also decrease. From the obtained results, we can conclude that after the adsorption of FPV on the doped C60 E g value decreases for FPV…FeC59, FPV…CrC59 and FPV… NiC59 complexes but this decrease in value is higher for Ni-doped C60. This means that the change in E g for the configuration of Ni-doped fullerene interacting with FPV drug is more as compared to the Cr and Fe doped fullerene. Thus, we concluded that an electronically favourable interactions occur which change the FPV drug properties significantly. This analysis gives a comparison of HOMO and LUMO energy levels of designed complexes and each constituent of a complex. Partial density of state (PDOS) shows the orbital distribution of individual component complex and total density of state (TDOS) shows the overall distribution of all orbitals in the complex. The change in E g value in complexes after adsorption can be visualized. The DOS of each complex is given in Fig. 4 , which reflects the change in height and displacement in their HOMO and LUMO levels. It can be estimated from the high relative intensity of drug and metal-doped complex that high conjugation results in high relative intensity. Observation of a change in DOS curves reflects the sensitivity of each complex which can be used to make an FPV drug a good electrochemical sensor. Furthermore, mulliken analysis represents the %age contribution around the FMOs (Fig. 5 , Table S1 ). The %age contribution of C59, FPV drug and Ni for HOMO is 65.7%, 11.2% and 23.1%, respectively, whereas C59 72.8%, FPV 17.9% and Ni 9.3% for LUMO. This contribution is also evident from DOS of FPV…NiC59. The %age composition of each building component in FPV…NiC59 is different for HOMO and LUMO. Likewise, C59, FPV drug and doped metal atom showed different %age contribution for HOMO and LUMOs in FPV…FeC59, FPV…CrC59 complexes. By comparing the %age contribution of all doped metals it can be visualized that contribution of Ni is higher for LUMO which results in efficient electronic transition (Table S1 ). Electronic properties are vigorously correlated to the HOMO and LUMO, their energy values and distribution patterns. Therefore, our major focus is on a detailed investigation of E g and the distribution pattern of FMOs. To observe the electronic sensitivity and to correlate it with electrical conductivity and energy gap, values were calculated through Eq. 8. An exponential relationship is present between electrical conductivity and E g . Compounds with Table 2 that after adsorption of the drug to doped fullerene conductivity has been increased. The values of the electrophilic index and chemical reactivity were applied to estimate the reactivity and chemical stability of the compound. E g plays a significant role in determining chemical stability. The polarizability, chemical hardness or softness, and reactivity are also calculated through formulas. To measure the reactivity of new molecules chemical hardness (Ƞ) is an important factor because molecules can easily alter their electronic density if they are structurally soft and biological activities are also higher for softer molecules. Moreover, for softer compounds, electronic transitions are better and easier which are essential for a chemical reaction. The values given in Table 2 show that after adsorption of FPV drug to metal doped C60 chemical hardness of complex molecule decreases as compared to FPV drug. Hence, doping of Cr, Fe and Ni causes a decrease in hardness. Therefore, a considerable decrease in hardness occurs for FPV…CrC59, FPV… FeC59, and FPV…NiC59 complexes. Thus, the comparison of all complexes indicates that complex FPV…NiC59 has greater softness so biologically more active and has higher adsorption power. The electrophilicity (ω) and maximum electronic charge accepted from the environment (ΔN max ) are important factors for determining the tendency of a molecule to adsorb electrons. Whenever two molecules react one acts as a nucleophile and the others as an electrophile. The compounds with low electrophilicity can transfer electrons and the compounds with high electrophilicity can accept electrons. The electrophilicity of FPV drug after binding reduces which means that it has a lower capacity to accept electrons after adsorption. The amount of charge transfer between drug and metallofullerene can also be explained through electrophilicity called electrophilicity-based charge transfer (ECT) can be calculated through the equation. The value of ECT indicates ECT < 0, so the drug behaves as a donor and the charge is transferred from drug to metal atom. The dipole moment is another good parameter to measure the solubility of compounds in polar solvents. Therefore, compounds with higher dipole moments are highly soluble in water. According to results, dipole moment increases after adsorption of FPV drug with doped fullerene due to which polarity of complex increases (Table S3) . Hence an increase in dipole moment is a favourable property to deliver the drug in the biological system. Based on calculated parameters, we conclude that after adsorption of the drug to fullerene Eg decreases for FPV…CrC59, FPV…FeC59, and FPV… NiC59 complexes due to which charge transfer from drug to fullerene increases, electrophilicity decreases and adsorption energy increases. Hence, this analysis indicates that complex of FPV drug with Ni-doped fullerene having the lowest Eg value than those of Cr and Fe doped complexes exhibits the highest adsorption energy. Molecular electrostatic potential surfaces (MEPS) are an important parameter to explain the relationship of electrophile and nucleophile attachment points. MEPS have been calculated to understand the nature of nucleophilic and electrophilic attacking sites. MEPS drawn from the optimized structure are shown in Fig. 6 . In MEPS, red and blue colours indicate the nucleophilic and electrophilic regions, respectively (Bibi et al. 2019) . MEPS maps show that for nucleophilic attack positive region is more favourable and for an electrophile to attack the negative region is more favourable. Possible interaction sites of FPV and doped fullerene have been purposed through MEPS maps. By considering the previous research and considering the interaction edge of the drug, it can be estimated from MEPS of FPV and Cr, Fe, and Ni-doped fullerene that for a ligand to interact, Fe, Cr, and Ni metal atoms are considered as the most active site because of positive charge distribution around the metal. In water solvent, the strongest interaction occurs from the -N, -OH edge of the drug and this edge is energetically more favourable than others. One reason for the favourable geometry of this edge is that there is a possibility of two bond formation (Rad et al. 2021) . The binding energy of FPV interacted with metal doped fullerene is given in Table 2 . From binding energy, it can suggest that chemisorption occurs between drug and doped molecule. In chemisorption new chemical bond formation occurs between interacting molecules, therefore interatomic distance must be at a feasible length. Thus, it can be concluded that favourable interactions occur in FPV…CrC59, FPV…FeC59, and FPV… NiC59 complexes. So, these metal doped fullerenes might be an advantageous candidate for drug delivery. It can be visualized from binding energy and electrophilicity index values ( Table 2 ) that among all metal atoms Ni doping significantly enhances the tendency of fullerene to adsorb FPV drug. In this work, DFT calculations were carried out with M06-2X/6-31G(d) method to explore the interaction between FPV drug and three metal-doped fullerenes (CrC59, FeC59, NiC59) to investigate the final three complexes (FPV… CrC59, FPV…FeC59, and FPV…NiC59) as drug delivery materials for COVID-19 treatment. All the calculations were carried out using Gaussian09 in water media. Among various interaction edges of FPV drug, the -N, and -OH edges were used to study the binding energy and other electronic parameters. After the interaction of the drug to doped fullerene, HOMO shifted to lower negative and LUMO shifted to higher negative values, which resulted in a decrease Eg. The compounds with lower Eg easily transfer electrons across the conduction band which resulted in higher electronic properties. The most negative binding energy was observed for FPV…NiC59, it showed that there exists the strongest interaction between FPV drug and Ni-doped fullerene. The decrease in Eg proved that Cr, Fe, and Ni-doped complexes are electronically harmless and significantly alter the FPV drug properties. Among all complexes, FPV…NiC59 exhibited the highest value of conductivity. An increase in the electrophilic index (ω) of metal-doped fullerene complexes as compared to a pure drug indicated that drug after adsorption with doped fullerene showed more electrophilicity for the observed complexes. Similarly, values of electrophilic charge transfer (ECT) showed that doped C59 acts as electron acceptor and drug as the electron donor. After the adsorption of FPV drug with the metal-doped fullerene, dipole moment value was increased which is a favourable property for these complexes to be used as drug delivery vehicle in a biological system for the purpose of antiviral diseases. In addition, a decrease in Eg and increase in binding energy for metal doped C60 fullerene was more than that of previously studied C20 doped complexes. Hence, based on these results, we can conclude that Cr, Fe and particularly Ni doped fullerene can efficiently adsorb FPV drug and would be the best nano-drug delivery material in COVID-19 treatment due to improved electronic and enhanced adsorption properties. NMR determination of solvent dependent behavior and XRD structural properties of 4-carboxy phenylboronic acid: a DFT supported study Utility of extrinsic [60] fullerenes as work function type sensors for amphetamine drug detection: DFT studies Effect of different topological structures (D-π-D and D-π-A-π-D) on the optoelectronic properties of benzo[2,1-B:3,4-B]dithiophene based donor molecules toward organic solar cells The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors Discovering drugs to treat coronavirus disease 2019 (COVID-19) Conformational and electronic properties of hydroquinone adsorption on C60 fullerenes: doping atom, solvent and basis set effects Coronavirus disease 2019 (COVID-19): a literature review Adsorption behavior of 5-fluorouracil on pristine, B-, Si-, and Al-doped C60 fullerenes: A firstprinciples study A DFT study on the functionalization of C60 fullerene with 1, 2-benzoquinone Carbon nanostructures: the drug and the delivery system for brain disorders The global impact of the coronavirus pandemic COVID-19: a call for physical scientists and engineers Fullerene-based delivery systems The electronic properties of carbenes C60-fullerenes as drug delivery carriers for anticancer agents: promises and hurdles Efficacy and safety of integrated traditional chinese and western medicine for corona virus disease 2019 (COVID-19): a systematic review and meta-analysis Molecular insight into the smart functionalized TMC-Fullerene nanocarrier in the pH-responsive adsorption and release of anti-cancer drugs Quantum mechanical simulation of Chloroquine drug interaction with C60 fullerene for treatment of COVID-19 A density functional theory investigation on amantadine drug interaction with pristine Computational study on favipiravir adsorption onto undoped-and silicon-decorated C60 fullerenes Investigation of detection and adsorption properties of β-propiolactone with silicon and aluminum doped fullerene C60 using density functional theory Synthesis, structural characterization, DFT studies and in-vitro antidiabetic activity of new mixed ligand oxovanadium (IV) complex with tridentate Schiff base Potential of metal-fullerene hybrids as strong nanocarriers for cytosine and guanine nucleobases: a detailed DFT study DFT calculations towards the geometry optimization, electronic structure, infrared spectroscopy and UV-vis analyses of Favipiravir adsorption on the first-row transition metals doped fullerenes; a new strategy for COVID-19 therapy DFT calculations towards the geometry optimization, electronic structure, infrared spectroscopy and UV-vis analyses of Favipiravir adsorption on the first-row transition metals doped fullerenes; a new strategy for COVID-19 therapy Advanced biomedical applications of carbon nanotube Tackling COVID-19 using remdesivir and favipiravir as therapeutic options Noncovalent interaction assisted fullerene for the transportation of some brain anticancer drugs: a theoretical study Nanotechnology in drug delivery and tissue engineering: from discovery to applications Experimental treatment with favipiravir for Ebola virus disease (the JIKI Trial): a historically controlled, single-arm proof-of-concept trial in Guinea Surface characterization of plasma treated polymers for applications as biocompatible carriers SARS-CoV-2 induced diarrhoea as onset symptom in patient with COVID-19 Efficacy and safety of favipiravir, an oral RNA-dependent RNA polymerase inhibitor, in mild-to-moderate COVID-19: a randomized, comparative, open-label, multicenter, phase 3 clinical trial An update on current therapeutic drugs treating COVID-19 Coronavirus disease 2019 (COVID-19): a perspective from China The calculations were performed in the computational and theoretical chemistry laboratory at the department of chemistry, University of Agriculture Faisalabad, Faisalabad, Pakistan. The authors vastly acknowledge the Higher education commission (HEC) of Pakistan for financial support. The authors also acknowledge Institute of Theoretical Chemistry and College of Chemistry, Jilin University, Changchun, China for technical support.