JBCS



22:32, qua dez 4

Acesso Aberto/TP




Artigo

Graphen-phenyl-nh2 as nanocarrier: a density functional theory study

Nosrat Madadi Mahani*; Fatemeh Mostaghni; Homa Shafiekhani

Department of Chemistry, Payame Noor University (PNU), 19395-4697 Tehran, Iran

Recebido em 02/10/2021
Aceito em 13/12/2021
Publicado na web em 16/02/2022

Endereço para correspondência

*e-mail: nmmadady@gmail.com

RESUMO

Recently, graphene and modified graphene as one of the most suitable and the most important carbon nanomaterials have been introduced for drug delivery. In this paper, we have studied the binding characteristics of the EDC-NHS cross-linking process of graphene-phenyl-NH2 and 5-aminolevulinic acid (ALA) drug in both gas and solvent phases by density functional theory calculations. For describing binding properties and reaction nature between graphene-ghenyl-NH2 and ALA drug, quantum molecular descriptors, topological analysis, natural bond orbital analysis, analysis of the bond order, the density of states, and analysis bond length was investigated in solvent and gas phases. Due to the results, the complex of the graphene-phenyl-NH2 @ALA turns to absorb more electrons in water solvent than gas phase. Furthermore, the binding of graphene-phenyl-NH2 and ALA is mainly based on covalent interactions, and bond order of graphene-phenyl-NH2 @ALA complex is one in solvent and gas phases. The praphene-phenyl- NH2 @ALA complex has displayed a meaningful improvement of electronic and structural properties. Therefore, it represented that praphene-phenyl-NH2 being combined with the ALA drug is appropriate for use in drug delivery.

Palavras-chave: graphene-Phenyl-NH2; 5-aminolevulinic acid drug; density functional theory; natural bond orbital analysis; topological analysis.

INTRODUCTION

Using the nanocarriers of anticancer drugs like micelles, liposomes, solid lipid nanoparticles, nanocages, carbon nanotubes, nanocrystals, and graphene oxide indicated an appropriate therapy of cancer by targeting specific cells of the tumor and these nanocarriers have no side effects.1-5 Drug delivery systems basis on the nanoscale to causes releasing on time and in place target of drug molecule.6 Also, drug nanocarriers can reduce drug toxicity.7 They provide modifications such as biocompatibility, half-life, solubility, and its release characteristics in properties of drug.8

Nanostructures as graphene sheets due to having properties such as light-weight, electron transport, structural strength, biocompatibility, flexibility, and having an extensive surface for double face adsorption in two dimensions can be used as sensor, cellular imaging, drug delivery and gene delivery agents.9,10 Graphene sheets are valuable for drug delivery due to solubility, and stability in the biological environment.9

Recently, researchers have displayed which graphene-phenyl-NH2 have cytoprotective properties.11 Graphene-phenyl-NH2 is the securest graphene derivative because it is not endowed with hemolytic and thrombotic properties in comparison to other graphene derivatives12 and is therefore can be used in biomedicine applications. Also, The edges and defect sites are well known to be more reactive than the basal plane of graphene. Therefore, the phenyl-NH2 group onto the zigzag egde is much more reactive than the central part of graphene.13

5-Aminolevulinic acid (ALA) is a prodrug that is used in photodynamic therapy in treating non-melanoma skin cancers due to its bioavailability.14,15 ALA has two active sites which include the carboxylic acid group and the amine group which the O-side of the carboxylic acid group can be covalently binding to graphene-phenyl-NH2. In the present study, the bonding of G-NH2 with the ALA drug was theoretically investigated. The bond between graphene-phenyl-NH2 with ALA drugs was appointed through a cross-linking compound that is correlated to the creating of amides bond betwixt drug and graphene-phenyl-NH2. This study aims to access the information about electronic characteristics, the natural electronic configurations by natural bond orbital analysis, chemical bond nature due to the topological analysis of electron density, the bond order analysis due to the laplacian electron density of the reaction of graphene-Phenyl-NH2 and ALA drug applying density functional theory computational method.

The interaction of ALA with SWCNTs (5, 5) and (9, 0), C24 and C60 has been studied with hybrid density functional theory calculations.16 Ab initio calculations for investigation the electronic and structural properties of ALA functionalized BNNTs have been performed by Samanata et al..17

Adsorption of cytarabine,18 ampyra,19 tegafur,20 and cladribine21 drugs anticancer on the graphene, and functionalized graphene as nanocarriers systems have been investigated using DFT calculations. Investigation binding of Doxorubicin drug with the graphene and functionalized graphene have been done using density functional theory computational by Vovusha et al..22 A theoretical study on the reactivity and the structural properties of the Melphalan anticancer drug with PAMAM dendrimer was done in both phases of gas and solvent.23

The addition of aryl radicals to perfect and defective graphene has been studied by dispersion-corrected density functional theory.24 Also, dispersion-corrected density functional theory has been utilized to study the interaction between alkyl radicals and graphene.25 Cycloadditions reactions onto grapheme like the [2 + 2] cycloadditions,26 Diels-Alder,27 [6+4] cycloadditions28 and, heteroatom promoted cycloadditions29 have been studied by means of density functional theory. Our study's results can help a better viewpoint of the reaction between graphene-phenyl-NH2 and ALA drug and proper information about the bonding of the reaction between graphene-phenyl-NH2 with ALA drug.

 

METHODS

At first, we applied GAUSSIAN 09 package30 to optimize all structures and quantum chemical computations in the solution and gas phases at Becke's three-parameter Lee-Yang-Parr hybrid functional31 with 6-311G basis set. density functional theory computations present us excellent viewpoints into structural geometry changes, interaction energy, and load distribution in the encapsulation process. Frequency computations are done to remove the imaginary frequencies and gaining minimum energy. Quantum theory of atoms in molecules and multifunctional wavefunction analyzer32 are applied to investigate the nature of interaction and binding among the atoms of graphene-phenyl-NH2 and ALA by the topological analysis of the electron densities. Moreover, Fuzzy bond order,33 Mayer bond order,34 and natural bond orbital (NBO) for graphene-phenyl-NH2 @ ALA were calculated.

Besides, the function of electron localization as a measurement of the likelihood of finding an electron in the neighborhood space of a reference electron located at a given point and with the same spin35was achieved by Multiwfn code. For simulating the biological system, configurations are surveyed in the solvent phase. Surveying solvent impacts was done using the integral equation formalism polarized continuum model of self-consistent reaction field.36 The descriptors of quantum chemical as an energy gap, and electronegativity and global hardness can be used to assess the chemical reactivity which is computed. The bond between graphene-phenyl-NH2 and ALA is made by a cross-linking applying EDC-NHS.

The final reaction between the ALA drug and graphene-phenyl-NH2 can be written as follows:

The bonding energies resulting from the reaction between the graphene-phenyl-NH2 and ALA can be acquired from equation (2).

where X is the enthalpy (H), the Gibbs free energy (G) or total energy (E), and Xcomplex is corresponds to the values of the graphene-phenyl-NH2 @ ALA.

 

RESULTS AND DISCUSSION

We have firstly optimized the structures of graphene-phenyl-NH2 and ALA in gas and solvent phases. The geometrically optimized structures of graphene-phenyl-NH2 and ALA in the solvent phase were shown in Figure 1. Moreover, the optimized configurations of graphene-phenyl-NH2 @ ALA in the solvent phase were indicated in Figure 2.

 


Figure 1. The optimized configurations of the ALA (a), and G-Phenyl-NH2 (b) at the solvent phase

 

 


Figure 2. The optimized configurations of the ALA@ G-Phenyl-NH2 at the solvent phase

 

The chemical reactivity descriptors of ALA, graphene-phenyl-NH2 and ALA@ graphene-phenyl-NH2 have shown in Table 1. The global hardness and the energy gap of ALA in the solvent phase is less than in the gas phase. Therefore, the stability of ALA is less in the water solvent and its reactivity is higher. The global hardness and energy gap of graphene-phenyl-NH2 in the solvent phase is greater than in the gas phase. After binding, the energy gap and global hardness of the ALA@ graphene-phenyl-NH2 complex decreased compared to pristine ALA in the solvent phase and electronegativity increased compared to ALA in the water solvent.

 

 

By surveying the electronegativity, it can be understood that graphene-phenyl-NH2 turn to absorb more electron than solvent phase. Therefore, the reaction of ALA and graphene-phenyl-NH2 causes an increase in the electronegativity and reduced their energy gap and global hardness in the solvent phase. The reaction energies of ALA and graphene-phenyl-NH2 in both phases are indicated in Table 2. In the ALA@ graphene-phenyl-NH2 system, Gibbs free energies are negative which shows the adsorption is spontaneous. Moreover, the negative enthalpies indicated that this reaction is exothermic. Furthermore, the polarity of molecules is shown by the dipole moment. The polarity of ALA@ graphene-phenyl-NH2 is higher in water solvent than in the gas phase.

 

 

Grimme has taken a different approach for correction of free energy based on the rigid rotor-harmonic oscillator (RRHO) approximation37 by arguing that low-frequency modes resemble free rotations and using the corresponding entropy term for low frequency modes. The free energy of association of molecules A and B to form a bound complex C, at temperature T in solvent X is given by equation 3:

which ΔE is different in energies of complex and reactants in the gas phase, is the sum of corrections from energy to free energy in the rigid-rotor-harmonic-oscillator approximation (RRHO) also including zero-point-vibrational energy for each species in the gas phase, and is their corresponding free solvation free energy. Based on equation 3, ΔEgas = -31.7024, = -4.3767 and, =1.0905 kJ mol-1, therefore ΔGa = -34.9886 kJ mol-1 which illustrated value of free energy in the phase condensed.

The saddle point on distributing electron density is named bond critical point which has the electron density gradient of zero. The bond critical point (BCP) is a first-order saddle point in the electron density distribution. The bond critical point is described by the major diameter elements of the Hessian matrix of electron density. The laplacian of the electron density is the trace of the Hessian matrix and the sum of the diagonal elements of the second derivative matrix of the density at a point is the Laplacian of electron density at that point.38 In regions where the Laplacian is negative, the potential energy is dominant and the negative charge is concentrated and, in regions where the Laplacian is positive, the kinetic energy density dominates and a depletion of negative charge occurs.

At bond critical point (BCP), quantities as the Laplacian, ∇2ρ(r), the potential energy density, V(r), the positive definite kinetic energy density, G(r), and the total electronic energy density, H(r) = V(r) + G(r), in terms of which the bond degree, H(r)/ρ(r), can be deliberated. In addition, equation 4, displays a helpful relationship between them:

The nature of chemical bonds could be distinguished by the values of potential and kinetic energy. The characteristics of chemical bonds can be shown by the amounts of potential and kinetic energy. The chemical bonds is covalent in case Laplacian and V(r) be negative, ρ(r) > 0.02, G(r)/ρ(r) and G(r)/V(r) < 1.39 Ionic interactions displays positive Laplacian and V(r), ρ(r) < 0.1, V(r)/ρ(r) and G(r)/V(r) >1. Electron density descriptors at the bond critical points for the significant chemical bonds of ALA@ graphene-phenyl-NH2 are showed in Table 3. Laplacian and E(r) of bonds are negative in gas and solvent phases. Moreover, the electron density values them are higher than 0.02. The analysis of quantities as -G(r)/V(r), E(r), and ρ(r), showed that these bonds are covalent kind.

 

 

Analyzing the function of electron localization was done to corroborate the characteristics of bonds. The electron localization function(ELF) is an easy measurement of electron localization in the atomic and molecular systems. The electron localization function analysis can give information about the binding, the location of bonds, the electron pairs, the atomic structure of the system, and the strength of the bond by integrating the electron density of the ELF area. To prove the covalent bonding, the electron localization function can be utilized.40 The electron localization function amount is between zero and one. By enhancing the electron localization function, the electrons become more localized and conceptualized, the characteristics of the atomic bond between two atoms gradually transform into a covalent kind. electron localization function amounts are in the range from 0 to 1. By increasing the electron localization function, the characteristics of the bond become to covalent kind due to the electrons become more concentrated and localized.41 The topographical analysis of electron localization function of ALA@ graphene-phenyl-NH2 is done which is indicated in Figure 3 in both phases. The upper limit nelectron localization function = 1 corresponding to perfect localization and the value electron localization function =½ corresponding to electron-gas-like pair probability. The electron localization function for the bod between reaction ALA@ graphene-phenyl-NH2 are almost 0.42. The electron localization function analysis characterized N-C bonds are partly covalent that approved electron density analysis. Analyzing the natural bond orbitals prepare the electronic density distribution on atoms and bonds which is a beneficent instrument to understand the delocalization of electron density.

 


Figure 3. ELF plots of optimized ALA@ G-Phenyl-NH2 (a) at gas and (b) at solvent phases

 

Table 4 indicates natural bond orbital (NBO) analysis after and before the reaction of ALA and graphene-phenyl-NH2 for main atoms in both phases. Table 4 data show that the electron densities of orbitals have changed. The electron density of 2p orbitals of N102 was reduced in the gas phase while the electron density of 2s was increased in the solvent phase after the reaction of ALA with graphene-phenyl-NH2. For the O106 atom, the electron density of 2s orbital decreased after adsorption, while the electron density of the 2p orbital increased in the gas phase. The electron density of 2s and 2p orbitals decreased after the reaction in water solvent. For the C99 atom, the electron density of the 2p orbital of C85 in both phases increased after reaction but the electron density of 2s orbital decreased in gas phase and increased in the solvent phase. For C104 atom, the electron densities of 2s and 2p orbitals in the solvent phase have reduced. After the reaction, the electron densities of 1s orbitals of the H103 atom have reduced by 0.02 and 0.01 in gas and solution phases, respectively.

 

 

The advancements of Wiberg bond order can be as Mayer and Fuzzy bond order that have confirmed extremely appropriate in bonding analysis applying Mulliken population analysis and semi-empirical calculation approaches. To show ALA and graphene-phenyl-NH2 bonds, the bond lengths, Fuzzy bond order, and Mayer bond order were calculated as listed in Table 5. The bond length of key atoms in the solution phase increased slightly relative to the gas phase, and therefore the bond order decreased. The bond length of N102-C104 in the solution phase is shorter than the gas phase and a shorter bond length indicated a stronger bond. The amounts of MBO show that bond order reduced after the reaction in both phases. The amounts of FBO show that the bond order between the drug and the graphene-phenyl-NH2 reduced after the reaction in both phases.

 

 

The density of the states of a system is defined as the number of various states which are to be occupied by the system at every energy. The density of the states computations permit one to identify the general distribution of states as an energy function. The density of the states spectra help better understanding the reaction between ALA and graphene-phenyl-NH2. DOS computation showed that there is not much between ALA and graphene-phenyl-NH2 in the gas phase and the water solvent. Therefore, the difference between HOMO and LUMO energies in the two phases is low. After reaction between ALA and graphene-phenyl-NH2 (see Figure 4), the difference between HOMO energies in two phases is 0.23 eV and also the value of 0.18eV for LUMO energies. Therefore, this reaction causes that in gas phase the HOMO level shifts to lower energy compared to the solvent phase, while the LUMO level is almost unchanged.

 


Figure 4. DOS plots of optimized ALA @ G-Phenyl-NH2 (a) at gas and (b) at solvent phases

 

CONCLUSIONS

The density functional theory computations were done to investigate the reaction of graphene-phenyl-NH2 and ALA anticancer drug. In this article, the bond of graphene-phenyl-NH2 and ALA drug was made by an EDC/NHS cross-linking. Our studies reveal that the ALA @ graphene-phenyl-NH2 complex can adsorb more electrons in water solvent. As well, the characteristics of bonding between ALA and graphene-phenyl-NH2 is the covalent. The binding energy of ALA @ graphene-phenyl-NH2 displayed that this energy reduces from -31.7024 eV in the gas phase to -35.0739 eV in the solvent phase. Likewise, the computations of Natural bond orbital (NBO) analysis, bond lengths, Fuzzy bond order, and Mayer bond order analysis of compounds were done to approve the covalent bond. The density of the states spectra of ALA @ graphene-phenyl-NH2 have the same trends in their energy gaps in gas and solvent phases. These binding characteristics can be useful in the appropriate selection of carbon nanomaterial for targeted drug delivery.

 

AKNOWLEDGEMENTS

We gratefully thank Payame Noor University for financial support. The authors would like to express their appreciations to department of chemistry this university for providing research facilities.

 

REFERENCES

1. Raza, K.; Thotakura, N.; Kumar, P.; Joshi, M.; Bhushan, S.; Bhatia, A.; Kumar, V.; Malik, R.; Sharma, G.; Guru, S. K.; Int. J. Pharm. 2015, 495, 551.

2. Naiafi, M.; Chin. J. Struct. Chem. 2019, 38, 195.

3. Spencer, D. S.; Puranik, A. S.; Peppas, N. A.; Curr. Opin. Chem. Eng. 2015 , 7, 84.

4. Kamel, M.; Raissi, H.; Morsali, A.; Shahabi, M.; Appl. Surf. Sci. 2018, 434, 492.

5. Khorram, R.; Morsali, A.; Raissi, H.; Hakimi, M.; Beyramabadi, S. A.; Chin. J. Struct. Chem. 2017, 36, 1639.

6. Dass, C. R.; Su, T.; Drug Del. 2001, 8,191.

7. Dubin, C. H.; Mech. Eng. Nanotechnol. 2004, 126, 10.

8. Kuma, M. N.; J. Pharm. Pharm. Sci. 2000, 3,234.

9. Yang, K.; Feng, L.; Shi, X.; Liu, Z.; Chem. Soc. Rev. 2013,42, 530.

10. Feng, L.; Liu, Z.; Nanomedicine 2011, 6, 317.

11. Singh, S. K.; Singh, M. K.; Kulkarni, P. P.; Sonkar, V. K.; Gracio, J. J. A.; Dash, D.; ACS Nano 2012, 6, 2731.

12. Rabchinskii, M. K.; Ryzhkov, S. A.; Kirilenko, D. A.; Ulin, N. V.; Baidakova, M. V.; Shnitov, V. V.; Pavlov, S. I.; Chumakov, R. G.; Stolyarova, D. Y.; Besedina, N. A.; Shvidchenko, A. V.; Potorochin, D. V.; Roth, F.; Smirnov, D. A.; Gudkov, M. V.; Brzhezinskaya, M.; Lebedev, O. I.; Melnikov, V. P.; Brunkov, P. N.; Sci. Rep. 2020, 10, 6902.

13. Jiang, D.; J. Chem. Phys. 2007, 126, 134701.

14. Baumgartner, R.; Pottier, R.; Krammer, B.; Stepp, H.; Photodynamic therapy with ALA, RSC Publishing: London, 2006.

15. Ishizuka, M.; Abe, F.; Sano, Y.; Takahashi, K.; Inoue, K.; Nakajima, M.; Kohda, T.; Komatsu, N.; Ogura, S.;Tanaka, T.; Int. Immunopharmacol. 2011, 11, 358.

16. Kia, M.; Golzar, M.; Mahjoub, K.; Soltani, A.; Superlattices Microst. 2013, 62, 251.

17. Samanta, P. N.; Das, K. K.; RSC Adv. 2016, 95, 92547.

18. Zaboli, M.; Raissi, H.; Rahmani Moghaddam, N.; Farzad, F.; J. Mol. Liq. 2020, 301, 112458.

19. Dastani, N.; Arab, A.; Raissi, H.; Adsorption 2020, 26, 879.

20. Shahabi, M.; Raissi, H.; J. Biomol. Struct. Dyn. 2017, 36, 2517.

21. Dastani, N.; Arab, A.; Raissi, H.; Struct. Chem. 2020, 31,1691.

22. Vovusha, H.; Banerjee, D.; Kumar Yadav, M.; Perrozzi, F.; Ottaviano, L.; Sanyal, S.; Sanyal, B.; J. Phys. Chem. C 2018, 122, 21031.

23. Ehsani, E., Shojaie, F.; J. Mol. Liq. 2021, 323, 114625.

24. Denis, P. A.; Chem. Phys. Chem. 2013, 14, 3271.

25. Denis, P. A.; Iribarne, F.; Chem. Eur. J. 2012, 18, 7568.

26. Denis, P. A.; Iribarne, F. ; J. Mater. Chem. 2012, 22, 5470.

27. Denis, P. A.; Chem. Eur. 2013, 19, 15719.

28. Denis, P. A.; Chem. Select 2017, 2, 9620.

29. Denis, P. A.; Chem. Select 2016, 1, 5497.

30. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H.; Li, X.; Caricato, M.; Marenich, A. V.; Bloino, J.; Janesko, B. G.; Gomperts, R.; Mennucci, B.; Hratchian, H. P.; Ortiz, J. V.; Izmaylov, A. F.; Sonnenberg, J. L.; Williams-Young, D.; Ding, F.; Lipparini, F.; Egidi, F.; Goings, J.; Peng, B.; Petrone, A.; Henderson, T.; Ranasinghe, D.; Zakrzewski, V. G.; Gao, J.; Rega, N.; Zheng, G.; Liang, W.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Throssell, K.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M. J.; Heyd, J. J.; Brothers, E. N.; Kudin, K. N.; Staroverov, V. N.; Keith, T. A.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A. P.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Millam, J. M.; Klene, M.; Adamo, C.; Cammi, R.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B.; Fox, D. J.; Gaussian 09, Gaussian, Inc Wallingford, CT, 2017.

31. Lopez, N.; Illas, F.; J. Phys. Chem. B 1998, 102,1430.

32. Lu, T.; Chen, F.; J. Comput. Chem. 2012, 33, 580.

33. Mayer, I.; Salvador, P.; Chem. Phys. Lett. 2004, 383, 368.

34. Mayer, I.; Chem. Phys. Lett. 1983, 97, 270.

35. Lu, T.; Chen, F. W.; Acta Phys.-Chim. Sin. 2011, 27, 2786.

36. Lipparini, F.; Scalmani, G.; Mennucci, B.; Frisch, M. J.; J. Chem. Theory Comput. 2011, 7, 610.

37. Grimme, S.; Chem. - Eur. J. 2012, 18, 9955.

38. Bankiewicz, B.; Matczak, P.; Palusiak, M.; J. Phys. Chem. A 2011, 116, 452.

39. Jacobsen, H.; Jacobsen, H.; J. Comput. Chem. 2009, 30, 1093.

40. Becke, A. D.; Edgecombe, K. E.; J. Chem. Phys. 1990, 92, 5397.

41. Fuentealba, P.; Chamorro, E.; Santos, J. C.; Theor. Comput. Chem. 2007, 19, 57.

On-line version ISSN 1678-7064 Printed version ISSN 0100-4042
Qu�mica Nova
Publica��es da Sociedade Brasileira de Qu�mica
Caixa Postal: 26037 05513-970 S�o Paulo - SP
Tel/Fax: +55.11.3032.2299/+55.11.3814.3602
Free access

GN1