key: cord-0781264-y21gjj5v authors: Pari, Azar Asgari; Yousefi, Mohammad title: Interactions between favipiravir and a BNC cage towards drug delivery applications date: 2021-09-08 journal: Struct Chem DOI: 10.1007/s11224-021-01833-8 sha: 46da86c92fabfce09827a1ef25f7538619131dfd doc_id: 781264 cord_uid: y21gjj5v Electronic structure analysis of bimolecular formation of favipiravir (Fav) and a representative model of boron-nitrogen-carbon (BNC) cage was performed in this work for providing more insightful information regarding the drug delivery purposes by the importance of Fav drug for medication of COVID-19. To achieve the purpose of this work, density functional theory (DFT) calculations were carried out to obtain the stabilized structures and corresponding molecular and atomic scale descriptors. Six models of BNC-Fav complexes were obtained reading the participation of different atomic positions of Fav to interactions with the BNC cage surface. The results yielded BNC-Fav2 at the highest strength and BNC-Fav4 at the lowest strength of bimolecular formations. Molecular orbital–related features and atomic scale quadrupole coping constants all revealed that BNC-Fav2 complex could be proposed for employing in drug delivery process by managing the loaded Fav contribution to future interactions. By the end of year 2019, the world was shocked by an attack of a wide-spreading type of coronavirus causing the coronavirus disease (COVID- 19) with harmful impacts on human health system and life worldwide [1] [2] [3] . To this date, the problem of treating COVID-19 has not been solved yet and the deathful impacts are still seen for the people [4] [5] [6] . Regarding the importance of this issue, considerable efforts of researchers have been dedicated to find a way of prevention or treatment of COVID-19, but no certain solution has been found up to now [7] [8] [9] . Therefore, performing further research works on the topic of COVID-19 is still a must to be done [10] [11] [12] . Among the investigated examinations, favipiravir (Fav) has been seen a useful drug for medication of coronavirus-infected patients to save their life or increase health level [13] [14] [15] . To this point, several attempts have been done to show advantages and disadvantages of medication of patients by Fav, in which there is still no certain achievement about its prescription [16] [17] [18] . Therefore, several researchers have tried to focus on improving the features of Fav for more efficient prescription and medication [19] [20] [21] . Employing carriers has been also proposed for targeted delivery of Fav to the infected tissues, in which such carriers should be carefully investigated for appropriate application [22] [23] [24] . Nanostructures have been always expected to play import role of carriers in drug delivery processes, in which they have been also assumed to work for delivery of Fav [25] [26] [27] . However, the works are still under improvement regarding achievements of researchers of various fields of study [28] [29] [30] . Several attempts have been always dedicated to develop procedures of identifications of components living systems [31, 32] . Among the methods of dealing with such problems, computer-based methods based on electronic structure concepts could help to provide insightful information for going forward to reach a proposed solution for such a complicated problem [33] [34] [35] . Indeed, molecular and atomic scale descriptors of matters could be obtained by computer-based works to investigate the characteristic features of matters at the smallest scales [36] [37] [38] . In this work, the electronic structure features of Fav were analyzed in combination with a representative nanostructure as an estimated carrier of Fav for drug delivery purposes. To approach this goal, quantum calculations were performed to stabilize each of the singular and bimolecular models to show the impacts of nanostructure on the characteristic features of loaded Fav. The representative nanostructure was a fullerene-like cage structure with boron, nitrogen, and carbon atomic counterparts with the stoichiometry of C 8 B 6 N 6 called a BNC cage. It is noted that cage structures have been seen as interesting for investigations by earlier works [39] [40] [41] . Singular models of Fav and BNC are shown in Figs. 1 and 2. After obtaining optimized singular structures, their combination for BNC-Fav bimolecular models (Fig. 3) was investigated by performing additional optimization calculations and evaluating the corresponding electronic structure features. The obtained molecular and atomic scale descriptors are summarized in Tables 1 and 2 The major materials of this work were Fav (C 5 H 4 FN 3 O 2 ) and BNC (C 8 B 6 N 6 ) models, in which their bimolecular formations were investigated by generating further materials ( Figs. 1-3) . First, the singular models were optimized to obtain the stabilized strictures. Next, the already singular models were combined with each other through various positioning faces resulting in six bimolecular models, from BNC-Fav1 up to BNC-Fav6 as shown in Fig. 3 . Hence, the stoichiometry of bimolecular models was identical but their relaxing configurations were different. To stabile the strictures, optimization calculations were performed using the WB97XD/6-31+G* level of density functional theory (DFT) as implemented in the Gaussian program [42] . By doing such calculations, the models were prepared for running further analyses of electronic structure features. To this aim, values of total energy (E tot ), interaction energy (E int ), and the basis set superposition error (BSSE) were evaluated for the model systems. Furthermore, molecular orbital features including energies of the highest occupied and the lowest unoccupied molecular orbitals (HOMO and LUMO), energy gap (E gap ), and chemical hardness and softness (η and σ) were evaluated for the optimized models. Additionally, values of dipole moment (D m ) were evaluated for the investigated models to show electric charge distribution orientation. All the obtained descriptors for the optimized models are summarized in Table 1 . In addition to quantitative descriptors, qualitative representations of optimized models are exhibited in Figs. 1-3, and distribution patterns for HOMO and LUMO are exhibited in Fig. 4 in addition to the visualized electrostatic potential (ESP) surfaces. As a result of performing such computations, the required features of molecular scale were provided for analyzing electronic structure features of the models. To analyze atomic features of Fav, values of quadrupole coupling constants (C q ) were evaluated for hydrogen, nitrogen, and oxygen atoms of Fav to show the impacts of bimolecular formation with BNC cage on the original features of Fav at the atomic scale. The obtained values of C q are summarized in Table 2 . The values of C q could show the localized electron density at the atomic sites with sensitivity to effects of any perturbations to the atomic sites especially for the interacting systems [43] [44] [45] . It is important to mention that several parameters are required for analyzing electronic structure features of matters, in which computer-based works could help to provide reliable results for approaching such goal [46] [47] [48] . The major goal of this work was to perform an electronic structure analysis of the bimolecular formation of favipiravir (Fav) and a representative BNC cage for drug delivery purposes through DFT calculations. To approach this goal, singular and bimolecular models were optimized to reach stabilized structures as presented in Figs. 1-3 . Although there were only two stating singular models for participating in bimolecular formations, six bimolecular models of BNC-Fav were obtained based on various relaxing configurations of molecular counterparts towards each other. As a consequence, the optimized models were prepared for evaluating molecular and atomic scale descriptors (Tables 1 and 2) , and for visualized HOMO and LUMO distribution patterns and ESP surfaces (Fig. 4) . To this point, all required quantitative and qualitative descriptors were evaluated for the investigated models to show the impacts of such bimolecular formation on the original features of Fav to propose a possible carrier for more efficient drug delivery purposes. It is important to mention that carbon nanostructures have been always expected to work properly for pharmaceutical applications in living systems [49] [50] [51] . However, their hydrophobicity raised an idea of developing impure carbon nanostructures by substituting other atoms instead of some or all available carbon atoms to reduce hydrophobicity feature [52] [53] [54] . In this case, combinations of boron and nitrogen (BN) atoms have produced excellent single-standing boron nitride nanostructures of even heteroatomic BN-substituted carbon nanostructures [55] [56] [57] . As a consequence, the investigated cage particle of this work was a combination of boron, nitrogen, and carbon atoms to construct the BNC cage (Fig. 2) for playing a possible carrier role of Fav (Fig. 1) . Indeed, heterosurface of the BNC cage could make possible the interacting molecules to contribute to different types of interactions in comparison with the identical surface of a pure carbon cage structure. Such benefit led to the formation of six bimolecular BNC-Fav models (Fig. 3) to analyze electronic structure features to show impacts of the existence of the BNC cage on the original features of Fav. The obtained bimolecular models of BNC-Fav complexes indicated the variation of complex formation regarding the relaxed configurations of molecular counterparts towards each other. In this case, those atomic sites of Fav counterpart with a potency of contributing to interactions participated in interactions with the heterosurface of the BNC cage. To make visual such note about such interacting bimolecular Fig. 4 could help to reach brighter points. As shown by different colors for ESP surfaces of each of singular Fav and BNC, various possibilities of interactions could be expected for the formation of such bimolecular complex systems. Blue color stands for positive, green color stands for neutral, and yellow and red colors stand for negative electrostatic parts of a molecule. As a result, the models showed colorful ESP for making them possible for participating in various interacting configurations. To this point, six models were obtained for BNC-Fav complex systems, in which their strengths are more or less different from each other. As listed in Table 1 , BNC-Fav2 and BNC-Fav4 were the strongest and the weakest models for the investigated complexes, respectively. Both values of E tot and E int approved such observation for the highest and lowest strength of BNC-Fav bimolecular models. To follow the structures visually, the corresponding panels of Fig. 3 showed the formation of a covalent bond for O1 atom of Fav and B atom of the BNC cage in addition to the formation of other interactions between two counterparts of the BNC-Fav2 complex. However, there were not any appropriate atomic positions for making strong interactions between the molecular counterparts in the BNC-Fav4 model resulting in the weakest complex systems for the BNC-Fav models. There might be a question about the importance of investigating formations of such weak complexes, in which it should be mentioned that several parameters could favor the formation of interacting systems even the formation of weak complex systems. Comparing to experiments, the efficiency of reaction products could be managed by conducting the molecules to be localized in the most appropriate configurations for interactions. This is indeed an advantage of performing such computer-based works to show the possibility of all complex formations for providing more insightful results for the investigated systems. For the investigated complexes, the order of strength could be ranked regarding the obtained values of both E tot and E int parameters: BNC-Fav2 > BNC-Fav3 > BNC-Fav1 > BNC-Fav5 > BNC-Fav6 > BNC-Fav4. The obtained values of BSSE were not meaningful to make significant changes to the obtained energies of interacting complex systems. Besides the participation of atomic sites of Fav in interactions with atomic sites of the BNC cage, whole structures of two molecules were also interacting together to make the BNC-Fav6 complex with somehow parallel positioning of Fav at the partial surface of the BNC cage. One important point about obtaining such different complexes is the existence of the heterosurface of the BNC cage in comparison with the homosurface of the pure carbon cage. The presence of impurities could increase the number of configurations suitable for the interactions [58] [59] [60] . Further analysis of the models based on HOMO and LUMO features could indicate variations of such molecular orbital energy levels for the complex models in comparison with singular forms of both of Fav and BNC cage. It is known that HOMO level could imply for the tendency of the molecule for electron donating whereas LUMO level could imply for such tendency for electron accepting. By such definitions, the investigated models showed different tendencies of participating in such electron-transferring processes. In addition to such levels, energy distance between HOMO and LUMO, assigned by E gap , is also another important feature for showing the reactivity modes of molecular systems. Additionally, each of the hardness and softness (η and σ) parameters could lead to the favorability of molecules for participating in electron-transferring processes. Such variations of electric features were also observed by the obtained values of D m . Examining the visual representations of HOMO and LUMO in Fig. 4 could show a very important achievement for the investigated model systems. As expected for nanostructures to work as proper carriers for drug delivery processes, it was seen that both HOMO and LUMO patterns were moved to the cage part in the strongest BNC-Fav2 model. This is indeed showing a dominant role of the employed BNC cage for managing the participation of loaded Fav in future interactions. Localization of HOMO and LUMO at the molecules could show the tendency of molecular parts for contributing to further interactions, in which such tendency was eliminated for Fav part in the BNC-Fav2 complex. As a result, it could be suggested that the formation of BNC-Fav2 complex might somehow reduce the unwanted side effects of medication of Fav for the patients, but further investigations are still required. As a consequence, BNC-Fav2 could be proposed for the best model of complex formation of Fav and BNC cage for further analyzing of applications developments. As mentioned about the importance of atomic positions of Fav for contributing to interactions with the BNC cage, the values of quadrupole coupling constants (C q ) were evaluated for H, N, and O atoms of Fav counterpart to show the impacts of complex formation at the smallest atomic scale features [61] [62] [63] . Different atomic types actually imply for different electronic environments, which are the most important factors for atoms to determine their participation in future interactions. The magnitude of such electronic density localization at the atomic site could be determined by values of C q , in which such magnitudes could be analyzed before and after the participation of atoms in interactions to show the perturbating impacts exactly. As shown in Fig. 1 , the atoms were numbered for being detectable in the atomic analysis regarding the obtained values of C q as listed in Table 2 . A quick look at the obtained results could mean that the electronic sites of all atoms of Fav detected a more or less significant impact of complex formation in BNC-Fav models. As detected by energy features, BNC-Fav2 was the strongest complex, in which analyzing values of C q for each of H1, N1, and O1 of singular and complex models show significant effects of interactions on the features of these atoms. Moreover, slight changes in electronic features were also detected for the rest of the atoms of Fav in the BNC-Fav2 model. Interestingly, the changes of atomic features of Fav from singular to complex model were almost negligible for the BNC-Fav4 model, as indicated earlier as the weakest complex formation for the investigated models. For other models, more or less significant impacts were observed for the atomic features by comparing the singular and complex models. Indeed, this is important how to decide for future functions of drug compounds in the drug delivery processes to employ the appropriate carrier for them. As a consequence, the atomic features of Fav were almost completely managed by the BNC cage assigning a carrier with significant impacts on the loaded drug. The major goal of this work was to perform an electronic structure analysis of the bimolecular formation of Fav and BNC cage for providing more insightful information regarding drug delivery purposes. To this aim, the obtained results indicated that such an idea of BNC-Fav bimolecular formation was achievable based on the obtained values of energies for whole structures and interacting modes. To this point, the obtained complex formations were ordered regarding their strength, in which BNC-Fav2 was found to be the strongest one, whereas BNC-Fav4 was found to be the weakest one. Such strengths were obtained by the importance of atomic types and positions of Fav counterpart in addition to the appropriate heteroatomic surface of the BNC cage. As a consequence, complexes with different strengths were obtained for the investigated models. Further analysis of molecular orbital features also indicated different tendencies for participating in interactions for the complex models even different from those singular counterparts. Atomic scale C q parameters showed the exact roles of atomic sites for the investigated models, in which their values were more or less significant for atoms of the whole structure. Interestingly, the most and the least significant effects were observed for atoms of BNC-Fav2 and BNC-Fav4, respectively, as were assigned for the highest and the lowest strength of the investigated BNC-Fav complexes. As a consequence, the BNC cage was seen useful for employing in the drug delivery process of loaded Fav for further investigations. The COVID-19 epidemic The impact of COVID-19 on oil supply in the short term Influence of the epidemic COVID-19: an outlook on health, business and scientific studies The impact of COVID-19 pandemic on corporate social responsibility and marketing philosophy COVID-19: a serious warning for emergency health innovation Work and well-being during COVID-19: impact, inequalities, resilience, and the future of work Coronavirus disease 2019 (COVID-19): Immunological approaches and emerging pharmacologic treatments. Int Immunopharmacol Screening efficacy of available HIV protease inhibitors on COVID-19 protease The pharmacist's active role in combating COVID-19 medication misinformation Post-acute COVID-19 syndrome A CRISPR-based and post-amplification coupled SARS-CoV-2 detection with a portable evanescent wave biosensor The impact of the COVID-19 pandemic on subjective mental well-being: the interplay of perceived threat, future anxiety and resilience Case studies of SARS-CoV-2 treated with favipiravir among patients in critical or severe condition Structure of the SARS-CoV-2 RNA-dependent RNA polymerase in the presence of favipiravir-RTP Angioedema and acute urticaria in a patient with COVID 19 pneumonia: Favipiravir side effect or COVID-19 cutaneous manifestation Efficacy of favipiravir in COVID-19 treatment: a multi-center randomized study Favipiravir: structural analysis and activity against COVID-19 Favipiravir for the treatment of patients with COVID-19: a systematic review and meta-analysis Clinical outcomes and plasma concentrations of baloxavir marboxil and favipiravir in COVID-19 patients: an exploratory randomized, controlled trial Safety and efficacy of favipiravir versus hydroxychloroquine in management of COVID-19: a randomised controlled trial The supramolecularly complexes of calix [4] arene derivatives toward favipiravir antiviral drug (used to treatment of COVID-19): a DFT study on the geometry optimization, electronic structure and infrared spectroscopy of adsorption and sensing 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 Favipiravir antiviral efficacy against SARS-CoV-2 in a hamster model DFT studies of graphene-functionalised derivatives of capecitabine Application of the nano-drug delivery system in treatment of cardiovascular diseases DNA codon recognition by a cubane wire: in silico approach Density functional theory analyses of non-covalent complex formation of 6-thioguanine and coronene Cyclophosphamide@CNT: in silico exploration of nano drug delivery system Advances in anti-breast cancer drugs and the application of nano-drug delivery systems in breast cancer therapy 2021) sgRNACNN: identifying sgRNA on-target activity in four crops using ensembles of convolutional neural networks BP4RNAseq: a babysitter package for retrospective and newly generated RNA-seq data analyses using both alignment-based and alignment-free quantification method Computational analyses of cytidine and aza-cytidine molecular structures Modifying a graphene layer by a thymine or a uracil nucleobase: DFT studies DFT studies of single lithium adsorption on coronene Computer-based tools for structural characterizations and activity specifications of natural products: a quick review DFT studies of 5-fluorouracil tautomers on a silicon graphene nanosheet 2021) 6-Thioguanine bimolecular formation for dual chelation of iron: DFT study Cl@ Si20H20]−: parent siladodecahedrane with endohedral chloride ion Stabilizing role of halide ions in endohedral [20] silafulleranes: insights from DFT calculations toward silicon nanocages Lithium adsorption at the C20 fullerene-like cage: DFT approach Solid-state NMR and NQR spectroscopy of lead-halide perovskite materials DFT explorations of quadrupole coupling constants for planar 5-fluorouracil pairs Computational studies of the purine-functionalized graphene sheets Surface interaction of H2O and H2S onto Ca12O12 nanocluster: quantum-chemical analyses Electronic structure study of the bimetallic Cu1-xZnx alloy thin films Electronic structure study of Sn-substituted InP semiconductor Complementary oligonucleotide conjugated multicolor carbon dots for intracellular recognition of biological events Carbon nanomaterials integrated molecularly imprinted polymers for biological sample analysis: a critical review Modified (n, 0) BN nanotubes (n= 3-10) by acetic acids: DFT studies An investigation of structural and magnetic properties of Ce-Nd doped strontium hexaferrite nanoparticles as a microwave absorbent Biomimetic synthesis of CuO nanoparticle using Capparis decidua and their antibacterial activity Nafion-based ionic-polymer-metal composites: displacement rate analysis by changing electrode properties Interactions of fluorouracil by CNT and BNNT: DFT analyses DFT calculations of a cubic B4N4 cubane-like particle for CO gas adsorption Investigations of adsorption behavior and anti-cancer activity of curcumin on pure and platinum-functionalized B12N12 nanocages Beryllium oxide (BeO) nanotube provides excellent surface towards adenine adsorption: a dispersioncorrected DFT study in gas and water phases A review on eco-friendly synthesis of BiVO4 nanoparticle and its eclectic applications Synthesis of M-Nd doped Fe3O4 nanoparticles (M= Co, Ce, Cr, Ni) with tunable magnetic properties A DFT approach on tioguanine: exploring tio-tiol tautomers, frontier molecular orbitals, IR and UV spectra, and quadrupole coupling constants Non-covalent functionalisation of C30 fullerene by pyrrole-n-carboxylic acid (n= 2, 3): density functional theory studies Kukolich SG (2021) Calculations and analysis of 55Mn nuclear quadrupole coupling for asymmetric top acyl methyl manganese pentacarbonyl The authors declare no competing interests.