Determination of structural properties in the adsorption of drugs on chitosan-hydrogels for type 2 diabetes by means of the PM6 method
Determinación de propiedades estructurales en la adsorción de fármacos en hidrogeles de quitosano para diabetes tipo 2 mediante el método PM6
Determinação das propriedades estruturais en adsorção de drogas em hidrogelos de quitosan para diabetes do tipo 2 pelo método PM6
DOI:
https://doi.org/10.15446/rev.colomb.quim.v49n2.84723Palabras clave:
Adsorption, diabetes, hydrogels, glibenclamide, metformin (en)adsorción, diabetes, hidrogeles, glibenclamida, metformina (es)
adsorption, diabetes, hydrogels, glibenclamide, metformin (pt)
Descargas
Computational chemistry performs the modeling and calculation of physicochemical properties that allow understanding of the different molecular interactions at the nanometric scale in medical applications such as the design of controlled release systems. The PM6 model was used to analyze metformin and glibenclamide. First, the energy properties as the Gibbs free energy and enthalpies were obtained. The results showed the affinity of both drugs with water (glibenclamide: -7.96 and metformin: -11.49) due to the formation of hydrogen bonds, which were verified by the electronegativities corresponding to the dipole moment and to the partition coefficient (Log P).
Subsequently, the main properties for the design of a release system using the metformin/glibenclamide complex in the chitosan hydrogel were determined. In this process it was appreciated that the Gibbs free energy (-2157.60 kcal/mol) determined the thermodynamic stability of the adsorption. In addition, the Log P (-25.82) indicated an instantaneous solubility through the formation of hydrogen bonds and were verified by the electronic distribution and the change in dipole moment.
La química computacional realiza el modelado y el cálculo de propiedades fisicoquímicas que permiten comprender las diferentes interacciones moleculares a escala nanométrica en aplicaciones médicas como por ejemplo, el diseño de sistemas de liberación controlada. El modelo PM6 se utilizó para analizar metformina y glibenclamida. Primero, se obtuvieron las propiedades energéticas como la energía libre de Gibbs y las entalpías. Los resultados mostraron la afinidad de ambos fármacos con el agua (glibenclamide: -7.96 y metformina: -11.49) debido a la formación de enlaces de hidrógeno, que fueron verificados por las electronegatividades correspondientes al momento dipolar y al coeficiente de partición (Log P).
Posteriormente, se determinaron las principales propiedades para el diseño de un sistema de liberación que usa el complejo metformina/glibenclamida en el hidrogel de quitosano. En este proceso se apreció que la energía libre de Gibbs (-2157.60 kcal/mol) determinó la estabilidad termodinámica de la adsorción. Además, el Log P (-25.82) indicó una solubilidad instantánea a través de la formación de enlaces de hidrógeno y se verificó mediante la distribución electrónica y el cambio en el momento dipolar.
A química computacional realiza a modelagem e o cálculo das propriedades físico-químicas que permitem compreender as diferentes interações moleculares em escala nanométrica em aplicações médicas, como o projeto de sistemas de liberação controlada. O modelo PM6 é usado para analisar metformina e glibenclamida. Primeiro, obtenha as propriedades energéticas como a energia livre de Gibbs e as entalpias. Os resultados mostram a afinidade de ambos os componentes com água (glibenclamida: -7.96 e forma: -11.49) debitados na forma de ligações de hidrogênio, que verificam por eletronegatividades correspondentes no momento dipolar e coeficiente de participação (Log P).
Posteriormente, selecione as principais propriedades para o projeto de um sistema de liberação que usa o método completo/glibenclamida no hidrogel de quitosano. Neste processo, aprecie a energia livre de Gibbs (-2157.60 kcal/mol) que determina a estabilidade termodinâmica da adsorção. Além disso, o Log P (-25.82) indica uma solução instantânea para a passagem de forma de hidrogênio e é verificada usando a distribuição eletrônica e o câmbio no momento dipolar.

Química Orgánica y Bioquímica
Determination of structural properties in the adsorption of drugs on chitosan-hydrogels for type 2 diabetes by means of the PM6 method
Determinación de propiedades estructurales en la adsorción de fármacos en hidrogeles de quitosano para diabetes tipo 2 mediante el método PM6
Determinação das propriedades estruturais en adsorção de drogas em hidrogelos de quitosan para diabetes do tipo 2 pelo método PM6
Determination of structural properties in the adsorption of drugs on chitosan-hydrogels for type 2 diabetes by means of the PM6 method
Revista Colombiana de Química, vol. 49, núm. 2, 2020
Universidad Nacional de Colombia
Recepción: 24 Enero 2020
Aprobación: 18 Abril 2020
Abstract:
Computational chemistry performs the modeling and calculation of physicochemical properties that allow understanding of the different molecular interactions at the nanometric scale in medical applications such as the design of controlled release systems. The PM6 model was used to analyze metformin and glibenclamide. First, the energy properties as the Gibbs free energy and enthalpies were obtained. The results showed the affinity of both drugs with water (glibenclamide: -7.96 and metformin: -11.49) due to the formation of hydrogen bonds, which were verified by the electronegativities corresponding to the dipole moment and to the partition coefficient (Log P). Subsequently, the main properties for the design of a release system using the metformin/glibenclamide complex in the chitosan hydrogel were determined. In this process it was appreciated that the Gibbs free energy (-2157.60 kcal/mol) determined the thermodynamic stability of the adsorption. In addition, the Log P (-25.82) indicated an instantaneous solubility through the formation of hydrogen bonds and were verified by the electronic distribution and the change in dipole moment.
Keywords: Adsorption, diabetes, hydrogels, glibenclamide, metformin.
Resumen:
La química computacional realiza el modelado y el cálculo de propiedades fisicoquímicas que permiten comprender las diferentes interacciones moleculares a escala nanométrica en aplicaciones médicas como por ejemplo, el diseño de sistemas de liberación controlada. El modelo PM6 se utilizó para analizar metformina y glibenclamida. Primero, se obtuvieron las propiedades energéticas como la energía libre de Gibbs y las entalpías. Los resultados mostraron la afinidad de ambos fármacos con el agua (glibenclamide: -7.96 y metformina: -11.49) debido a la formación de enlaces de hidrógeno, que fueron verificados por las electronegatividades correspondientes al momento dipolar y al coeficiente de partición (Log P). Posteriormente, se determinaron las principales propiedades para el diseño de un sistema de liberación que usa el complejo metformina/glibenclamida en el hidrogel de quitosano. En este proceso se apreció que la energía libre de Gibbs (-2157.60 kcal/mol) determinó la estabilidad termodinámica de la adsorción. Además, el Log P (-25.82) indicó una solubilidad instantánea a través de la formación de enlaces de hidrógeno y se verificó mediante la distribución electrónica y el cambio en el momento dipolar.
Palabras clave: adsorción, diabetes, hidrogeles, glibenclamida, metformina.
Resumo:
A química computacional realiza a modelagem e o cálculo das propriedades físico-químicas que permitem compreender as diferentes interações moleculares em escala nanométrica em aplicações médicas, como o projeto de sistemas de liberação controlada. O modelo PM6 é usado para analisar metformina e glibenclamida. Primeiro, obtenha as propriedades energéticas como a energia livre de Gibbs e as entalpias. Os resultados mostram a afinidade de ambos os componentes com água (glibenclamida: -7.96 e forma: -11.49) debitados na forma de ligações de hidrogênio, que verificam por eletronegatividades correspondentes no momento dipolar e coeficiente de participação (Log P). Posteriormente, selecione as principais propriedades para o projeto de um sistema de liberação que usa o método completo/glibenclamida no hidrogel de quitosano. Neste processo, aprecie a energia livre de Gibbs (-2157.60 kcal/mol) que determina a estabilidade termodinâmica da adsorção. Além disso, o Log P (-25.82) indica uma solução instantânea para a passagem de forma de hidrogênio e é verificada usando a distribuição eletrônica e o câmbio no momento dipolar.
Introduction
Nanotechnology has allowed for the development of technologies for the early identification of diseases, a field where nanomedicine represents the main area of application. The applications of nanomedicine consist in the formation of images, the detection, and diagnosis of diseases as well as controlled delivery systems which increase the effectiveness of the release, the tolerability, the specificity of the drugs, and the reduction of side effects, which has increased the availability of new medical treatments [1-4]. Nanomedicine allows the study of different types of molecules such as proteins, drugs, vaccines or nucleic acids to carry out the diagnosis and treatment for cancer, diabetes and tissue engineering [5-6].
The World Health Organization (WHO) predicts a population of 500 million adults affected with type 2 diabetes by the year 2030 [7]. The design of new drug administration systems allows the penetration in the cellular tissues to carry out the glucose measurement and later perform the controlled release of antidiabetic drugs, such as glibenclamide and metformin, in complete doses. These delivery systems are synthesized using polymers, micelles, liposomes, and biodegradable polymer nanoparticles [8]. Release systems represent the main application of hydrogels, which are obtained by the dispersion in water of a polymer solution and are characterized by the degree of swelling, viscosity, pH, biodegradation, and degree of crosslinking [9].
Hydrogels are used as intelligent transporters due to their swelling capacity, release kinetics, and biocompatibility, which have generated more research at the molecular level in order to modify their structure and optimize their efficiency in pharmaceutical applications [10]. The main polymers include chitosan, poly(lactic acid-co-glycolic acid), and dextran [11].
Chitosan is a polysaccharide that allows for the release of a drug in a hydrogel due to the hydrophilic nature of functional groups such as OH, COOH, and NH. through the formation of hydrogels [8]. It also allows the adsorption of nanoparticles due to the rupture of the internal structure of enterocyte to reach the bloodstream, obtaining a controlled release system sensitive to external stimuli [11]. Chitosan hydrogels have been used in anticancer therapies as well as in controlled release systems [12-14].
Currently, computer models are used to analyze, design, and evaluate the effectiveness of new drugs with physical-chemical properties and specific biological activities, reducing the problems of effectiveness in medical treatments [15]. Thus, the pharmaceutical industry has increased the design of release systems using classical force fields or quantum mechanical models for the modeling of biological systems or the design of drugs where the calculation of Gibbs free energy is obtained for several thousand atoms simultaneously [16].
Quantitative Structure-Activity Relationship (QSAR) properties allow the calculation of non–covalent interactions, that is, they allow predicting the structure and biological behavior of macromolecules and their interactions with drugs, receptors or drug–receptors using descriptors such as thermodynamic, lipophilicity, physicochemical, electronic, and biological properties, solubility, and diffusion coefficient [3, 17-20]. The partition coefficient (Log P) is one of the most relevant QSAR properties because it provides information about the effectiveness of a drug to cross biological barriers and to reach the specific site for drug release. Log P defines the intermolecular behavior of the drugs in the human body, being widely used in the pharmaceutical sector [21].
The calculation time in DFT and AB-initio results in a higher computational cost due to the size and type of systems that are analyzed. However, the semi-empirical models of quantum mechanics simplify these times due to the incorporation of parametrizations as well as the exclusion of functions of different atoms. The PM6 model of quantum mechanics was developed for 9 000 compounds from experimental bases in conjunction with the AB–initio method, which predicts the heat of formation and the geometry of optimization more accurately because information is obtained from the interactions of atomic, hydrogen bonds as well as the use of d orbitals that allow the study of transition metals [16, 22-24].
In this work, the PM6 model was used to determine the optimization geometry, QSAR properties and the molecular electrostatic potential (MESP) in the adsorption of a metformin/glibenclamide complex on chitosan hydrogels.
Materials and Methods
All the simulations were performed in an Intel(R) Core (TM) i7 CPU @ 2.53 GHz with 4 GB RAM, having windows XP environment. The simulations involved metformin, gliblenclamide, metformin/water, and glibenclamide/water, as well as the complex (metformin/glibenclamide) where the calculation of the energetic properties and the molecular orbitals were made by selecting the PM6 semi–empirical method. Additionally, for determining the solvation energy (Aq) in the different simulations, the Cramer–Truhlar SM3 method, a SM5.4 solvation of the Spartan software was used [24].
Subsequently, using the Compute menu the QSAR properties such as surface area, volume, mass, and polar surface area were obtained. The Log P was calculated using the Spartan software by selecting the Ghose–Crippen method. This method is independent of the wave energy function, that is, the calculation is not influenced whether it is obtained using a quantum mechanics, molecular mechanics or a semi-empirical method. HOMO and LUMO molecular orbitals were obtained by solving the wave function by means of the PM6 semi–empirical method.
Tables 1 and 2 show the thermodynamic and QSAR properties of water, metformin, and glibenclamide, where a negative Gibbs free energy was obtained. Besides, the solvation energy indicates the dissolution of both drugs in water through hydrogen bonds. Electrostatic Potential maps (MESP) were calculated in the Spartan program by selecting the Hartree-Fock method with a 3-21G set. Figure 1 shows the MESP of water, metformin, and glibenclamide. It can be observed that the regions of high electronic density (in red) are presented in N-H, S=O, C=O and O-H bonds, while the regions of low electronic density are presented in the CH bonds.
| Properties | Units | Water | Metformin | Glibenclamide |
| Gibbs free energy (ΔG) | kcal/mol | - 47426 | - 268288 | - 1414298 |
| Solvation energy | kcal/mol | - 12.55 | - 13.44 | - 28.01 |
| Enthalpy of formation (H°f) | kJ/mol | - 54.27 | 46.46 | - 191.84 |
| Enthalpy of aqueous formation (H°f) | kJ/mol | - 66.82 | 33.01 | - 219.79 |
| Standard entropy (S°) | kJ/mol | 0.450 | 0.094 | 0.167 |
| Dipolar moment | Debye | 1.86 | 2.99 | 11.23 |
| E HOMO | eV | - 11.91 | - 9.2 | - 9.44 |
| E LUMO | eV | 4.070 | 0.35 | - 0.91 |
Thermodynamic properties of different molecules.
| Properties | Units | Water | Metformin | Glibenclamide |
| Mass | amu | 18 | 129.16 | 494.012 |
| Log P | ------- | - 0.38 | - 0.50 | 3.80 |
| Surface area | Ų | 36.36 | 172.2 | 512.04 |
| Polar surface area | A² | 36.35 | 78.27 | 99 |
| Polar surface | % | 99.9 | 44 | 19.33 |
| Volume | ų | 19.22 | 137.53 | 471.83 |
QSAR properties of different molecules.

Figure 1.
MESP of (a) water, (b) metformin, and (c) glibenclamide using the PM6 model. The red color corresponds to regions of high electronic density and the blue color to regions of low electronic density.
Results and Discussions
Drugs/water
Tables 3 and 4 show the energetic and QSAR properties of the metformin/water and glibenclamide/water solutions at 298.15 K, where the Gibbs free energy (ΔG) indicates the spontaneous formation of both individual aqueous solutions attributed to a strong interaction of the OH bonds of the water with the receptor groups of both drugs and to the drugs polar surface originated by the presence of atoms with different electronegativity in each molecule.
The existence of an electrostatic attraction between water and both drugs caused by the rupture of the individual water ions to obtain a thermodynamic equilibrium was also observed, and then corroborated by the negative partition coefficient shown in Table 4 and the dipole moment indicating a conformational change in the structure of both solutions caused by the electronegativity of OH, CO, SO in both drugs [25-26].
| Properties | Units | Metformin/water | Glibenclamide/wáter |
| Gibbs free energy (ΔG) | kcal/mol | - 1429 | - 1503 |
| Enthalpy of formation (H°f) | kJ/mol | - 1818 | - 2076 |
| Enthalpy of aqueous formation (H°f) | kJ/mol | - 1324 | - 1378 |
| Standard entropy (S°) | kJ/mol | 0.354 | 0.42 |
| Dipolar moment | Debye | 6.48 | 11.46 |
| E HOMO | eV | - 9.49 | - 9.37 |
| E LUMO | eV | 0.26 | - 1.18 |
Energetic properties of antidiabetic drugs in water.
| Properties | Units | Metformin/water | Glibenclamide/water |
| Mass | amu | 669.617 | 1,034.462 |
| Log P | ------- | - 11.49 | - 7.96 |
| Surface area | Ų | 811.53 | 1,077.61 |
| Polar surface area | A² | 727.94 | 706.07 |
| Polar surface | % | 89.70 | 65.55 |
| Volume | ų | 638.55 | 961.47 |
| Ovality | -------- | 2.26 | 2.29 |
QSAR properties of antidiabetic drugs in water.
The MESP of the dissolution of both drugs in water are depicted in Figure 2. The electronic distribution shows that the regions of high electronic density were located in the OH, C=O and S=O bonds, while the regions of low electronic density were observed in the CH bonds. Figure 2(a) metformin shows an attraction between hydrogen in water and oxygen in the OH, while in Figure 2(b) it is observed that hydrogen bonds originate between hydrogen in water and the carbonyl group. This behavior corroborates the formation of a new bond that generates the dissolution of both drugs with water.

Figure 2.
MESP of the (a) Metformin/water and (b) Glibenclamide/water using the PM6 model.
Hydrogel/Complex (metformin/glibenclamide)
Table 5 shows the adsorption of the complex (metformin/glibenclamide) on the chitosan hydrogel because the Gibbs free energy and enthalpy indicate an exothermic process, that is, a spontaneity in the adsorption process produced by the increase of interactions in the hydrogel/complex. Also, the dipole moment reveals that the decrease of the dipole-dipole forces produces a greater solubility due to the increase of the intermolecular forces through the formation of hydrogen bonds, generating geometry with the greater structural stability of the hydrogel/complex [27-28].
The formation of hydrogen bonds is attributed to the increase in the polar surface (Table 6). Also, the existence of the affinity of the complex (metformin/glibenclamide) with the water present in the hydrogel is observed due to the negative value of the partition coefficient that indicates that the drugs are related to the polar behavior of the water. Additionally, the surface area contributes to penetrate the cell membrane and may be in contact with the bloodstream to establish a drug delivery system [29-30]. The cationic character of chitosan allows obtaining stimuli sensitive nanohydrogels to the glucose [14] due to the protonation of the amine groups that produce the adsorption of drugs in the hydrogel/complex.
Additionally, a new controlled release system is obtained since the nanometric particles of the complex have greater stability, bioavailability and longer life time, being thus more effective causing the rupture of the permeability barriers, achieving the controlled release [31].
| Properties | Units | Hydrogel/Complex |
| Gibbs free energy (ΔG) | kcal/mol | - 2157.60 |
| Enthalpy of formation (H°f) | kJ/mol | - 3784.42 |
| Enthalpy of aqueous formation-Aq (H°f) | kJ/mol | - 1920.15 |
| Standard entropy (S°) | kJ/mol | 0.7945 |
| Dipolar moment | Debye | 5.97 |
| E HOMO | eV | - 8.81 |
| E LUMO | eV | - 0.83 |
| Ovality | ---------- | 2.86 |
Energetic properties of hydrogel/complex (metformin/glibenclamide).
| Properties | Units | Hydrogel/Complex |
| Mass | amu | 2749.16 |
| Log P | ------- | - 25.82 |
| Surface area | Ų | 2545.02 |
| Polar surface area | A² | 1348.59 |
| Polar surface | % | 52.96 |
| Volume | ų | 2494.83 |
| Polarizability | ------- | 39.49 |
QSAR properties of hydrogel/complex.
Figure 3 shows the map of the electronic distribution of the adsorption of the drugs (glibenclamide and metformin) in the chitosan hydrogel, where the deformations present in the rings of glibenclamide are attributed to the functional groups OH and NH of the chitosan. Additionally, the regions of high electronic density zones correspond to the polar groups on the surface of the hydrogel, attributing the adsorption of both drugs to the formation of hydrogen bonds. In addition, in the sulphonylurea group of glibenclamide, a positive electronic density is observed, while in metformin, ionic interactions producing electrostatic forces that corroborate the transfer of charge with the hydrogel are observed [32-34].

Figure 3.
MESP of the hydrogel/(metformin/glibenclamide) complex using PM6 model.
Conclusions
The computational analysis showed that there is a process of adsorption of metformine/glibenclamide complex using chitosan hydrogels due to the thermodynamic values as a Gibbs free energy of -2157.60 kcal/mol. The negative value of Gibbs free energy determined a spontaneous reaction which corresponded to the formation of hydrogen bonds. This was corroborated by the negative value of the coefficient partition (Log P) of -25.82, which indicated the hydrophilic capacity increases in the adsorption and the dipole moment change as a product of a greater electronegativity attributed to hydrogel.
Moreover, the surface area indicated that both drugs are dispersed molecularly, that is, there would be no problems of crystallinity in the hydrogel attributable to the interactions and the size of the pores in the hydrogel that produce thermodynamic stability. Finally, through the enthalpy of formation and the enthalpy of aqueous formation it was possible to determine that glibenclamide has a greater affinity for hydrogel due to the volume occupied by the molecule in the hydrogel. The OH and SO functional groups of glibenclamide increased the water solubility of the hydrogel as well as the presence of secondary bonds such as the van der Waals forces.
References
[1] R. Gupta, <Diabetes Treatment by Nanotechnology,> J. Biotechnol. Biomater. vol. 7, n° 3, pp. 1-2, 2017. DOI: https://doi.org/10.4172/2155-952X.1000268
[2] B. B. Kadhim, H. O. Muhsen, <Structural and Electronic Properties of SWGaPNT Drug Carrier,> Nanosci. Nanotechnol. vol. 7, n° 1, pp. 9-13, 2017. DOI: https://doi.org/ 10.5923/j.nn.20170701.03
[3] N. Haddish-Berhane, J. L. Rickus, K. Haghigh, <The role of multiscale computational approaches for rational design of conventional and nanoparticle oral drug delivery systems,> Int. J. Nanomed. vol. 2, n°3, pp. 315–331, 2007.
[4] P. S. Appasaheb, S. D. Manohar, S. R. Bhanudas, <A review on intranasal drug delivery system,> J. Adv. Pharm. Edu. & Res. vol. 3, n° 4, pp. 333-346, 2013. DOI: https://doi.org/10.5958/2321-5844.2018.00004.3
[5] G. Sharma, A. R. Sharma, J. S. Nam, G. P. C. Doss, S. S. Lee, C. Chakraborty, <Nanoparticle based insulin delivery system: the next generation efficient therapy for Type 1 diabetes,> J. Nanobiotechnol. vol. 13, pp. 1-13, 2015. DOI: https://doi.org/10.1186/s12951-015-0136-y
[6] R. M. Di-Santo, V. Subramanian, Z. Gu, <Recent advances in nanotechnology for diabetes treatment,>” Wiley Interdiscip. Rev. Nanomed. Nanobiotechnol.vol. 7, pp. 548–564, 2015. DOI: https://doi.org/10.1002/wnan.1329
[7] O. Veiseh, B. C. Tang, K. A. Whitehead, D. G. Anderson, R. Langer, <Managing diabetes with nanomedicine: challenges and opportunities,> Nat. Rev. Drug Discov. vol. 14, pp. 45–57, 2015. DOI: https://doi.org/10.1038/nrd4477
[8] C. S. Satish, K. P. Satish, H. G. Shivakumar, <Hydrogels as controlled drug delivery systems: Synthesis, crosslinking water and drug transport mechanism,> Indian J. Pharma. Sci. vol. 68, pp. 133-140, 2006. DOI: https://doi.org/10.4103/0250-474X.25706
[9] I. S. Raja, N. N. Fathima, <Porosity and dielectric properties as tools to predict drug release trends from hydrogels,> Springerplus. vol. 3, pp. 393-402, 2014. DOI: https://doi.org/ 10.1186/2193-1801-3-393
[10] P. A. Arredondo, M. E. Londoño, <Hydrogels: potentials biomaterials for controlled drug delivery,> Rev. Ing. Biomédica. vol. 3, n° pp. 83-94, 2009. DOI: https://doi.org/10.24050/19099762.n5.2009.70
[11] V. R. Devadasu, T. M. Alshammari, M. Aljofan, <Current advances in the utilization of nanotechnology for the diagnosis and treatment of diabetes,> Int. J. Diabetes Dev. Ctries. vol. 38, n° 1, pp. 11-19, 2018. DOI: https://doi.org/10.1007/s13410-017-0558-1
[12] C. Dalwadi, P. Gayatri, <Application of nanohydrogels in drug delivery systems: recent patents review,> Rec. Pat. Nanotech. vol. 9, n° 1, pp. 17-25, 2015. DOI: https://doi.org/ 10.2174/1872210509666150101151521
[13] S. K. Basu, A. Rajendran, <Studies in the development of nateglinide loaded calcium alginate and chitosan coated calcium alginate beads,> Chem. Pharm. Bull. vol. 56, n° 8, pp. 1077-1084, 2008. DOI: https://doi.org/10.1248/cpb.56.1077
[14] L. J. Del-Valle, A. Díaz, J. Puiggalí, <Hydrogels for biomedical applications: cellulose, chitosan, and protein/peptide derivatives,> Gels. Vol. 3, n° 3, 2017. DOI: https:/doi.org/10.3390/gels3030027
[15] R. Galeazzi, <Molecular dynamics as a tool in rational drug design: Current status and some major applications,> Curr. Comput-Aid Drug. vol. 5, pp. 225-240, 2009. DOI: https:/doi.org/10.2174/157340909789577847
[16] A. V. Sulimov, D. C. Kutov, E. V. Katkova, V. B. Sulimov, <Combined docking with classical force field and quantum chemical semiempirical method PM7,> Adv. Bioinformat. vol. 2017, n° 5, pp. 1-6, 2017. DOI: https:/doi.org/10.1155/2017/7167691
[17] J. L. Medina-Franco, F. López-Vallejo, R. Castillo, <Diseño de fármacos asistido por computadora,> Educ. Química. vol. 17, pp. 452-457, 2006. DOI: https://doi.org/10.22201/fq.18708404e.2006.4.66027
[18] J. Liu, H. Zheng, X. Dai, S. Sun, H. G. Machens, A. F. Schilling, <Biomaterials for promoting wound healing in diabetes,> J. Tissue Sci. Eng. vol. 8, pp. 193, 2017. DOI: https:/doi.org/10.4172/2157-7552.1000193
[19] L. G. Ferreira, R. N. Dos-Santos, G. Oliva, A. D. Andricopulo, <Molecular docking and structure-based drug design strategies,> Molecules.vol. 20, pp. 13384-13421, 2015. DOI: https:/doi.org/10.3390/molecules200713384
[20] M. V. Vellayappan, J. R. Venugopal, S. Ramakrishna, S. Ray, A. F. Ismail, M. Mandal, A. Manikandan, S. Seal, S. K. Jaganathan, <Electrospinning applications from diagnosis to treatment of diabetes,> RSC Adv. vol. 6, pp. 83638–83655, 2016. DOI: https:/doi.org/ 10.1039/C6RA15252J
[21] A. Kovalenko, <Molecular theory of solvation: Methodology summary and illustrations,> Condens. Matter. Phys. vol. 18, n° 3, pp. 1-24, 2015. DOI: https://doi.org/10.5488/CMP.18.32601
[22] J. J. P- Stewart, <Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements,> J. Mol. Model. vol. 13, n° 12, pp. 1173-1213, 2007. DOI: https://doi.org/10.1007/s00894-007-0233-4
[23] J. J. P. Stewart, <Application of the PM6 method to modeling the solid state,> J. Mol. Model. vol. 14,n° 6, pp. 499–535, 2008. DOI: https://doi.org/10.1007/s00894-008-0299-7
[24] L. O. Olasunkanmi, J. Ige, G. O. Ogunlusi, <Theoretical study of the molecular geometries, electronic and thermodynamic properties of chlorinated dipyrido-(3,2-a:2 ′ ,3 ′ -c)-phenazine,> J. Chem. vol. 2013, pp. 1-7, 2013. DOI: https://doi.org/10.1155/2013/640649
[25] G. Schaftenaar, J. Vlieg, <Quantum mechanical polar surface area,> J. Comput. Aid. Mol. Des. vol. 26, n° 3, pp. 311–318, 2012. DOI: https://doi.org/10.1007/s10822-012-9557-y
[26] C. W. Pouton, <Formulation of poorly water-soluble drugs for oral administration: Physicochemical and physiological issues and the lipid formulation classification system,> Eur. J. Pharm. Sci. vol. 29, n° 3-4, pp. 278-287, 2006. DOI: https://doi.org/ 10.1016/j.ejps.2006.04.016
[27] E. Mah, R. Ghosh, <Thermo-responsive hydrogels for stimuli-responsive membranes,> Processes. vol. 1, pp. 238-262, 2013. DOI: https://doi.org/10.3390/pr1030238
[28] Q. Chai, Y. Jiao, X. Yu, <Hydrogels for biomedical applications: Their characteristics and the mechanisms behind them,> Gels. 2017. DOI: https://doi.org/10.3390/gels3010006
[29] C. F. Santa, B. L. López-Osorio, <Materiales poliméricos en nanomedicina: transporte y liberación controlada de fármacos,> Rev. Acad. Colomb. Cienc. Exact. Fis. Nat. vol. 37, pp. 115-124, 2013.
[30] S. Hurst, C. M. Loi, J. Brodfuehrer, A. El-Kattan, <Impact of physiological, physicochemical and biopharmaceutical factors in absorption and metabolism mechanisms on the drug oral bioavailability of rats and humans,> Expert. Opin. Drug Metab. Toxicol.vol. 3, n° 4, pp. 469-489, 2007. DOI: https://doi.org/10.1517/17425225.3.4.469
[31] S. Jana, S. Maiti, <Chitosan-based nanoparticulate systems for oral drug delivery,> in Nanostructures for oral medicine: Nanostructures in therapeutic medicine series, E. Andronescu, A. M. Grumezescu, Edits., Elsevier, 2017, pp. 607-638.
[32] I. Amăriuţei, M. Florea, M. V. Ghica, L. Popa, <Microemulsion-based strategy for the capillary electrophoretic analysis of some antidiabetic drugs,> Romanian J. Biophys.vol. 25, n° 3, pp. 173–184, 2015.
[33] D. Horter, J. B. Dressman, <Influence of physicochemical properties on dissolution of drugs in the gastrointestinal tract,> Adv. Drug. Deliv. Rev. vol. 46, pp. 75-87, 2001. DOI: https://doi.org/10.1016/s0169-409x(00)00130-7
[34] R. Galeazzi, <Molecular dynamics as a tool in rational drug design: Current status and some major applications,> Curr. Comput-Aid Drug, vol. 5, pp. 225-240, 2009. DOI: https://doi.org/10.2174/157340909789577847
Notas de autor
polymer1979@yahoo.com
Recibido: 24 de enero de 2020; Aceptado: 18 de abril de 2020
Abstract
Computational chemistry performs the modeling and calculation of physicochemical properties that allow understanding of the different molecular interactions at the nanometric scale in medical applications such as the design of controlled release systems. The PM6 model was used to analyze metformin and glibenclamide. First, the energy properties as the Gibbs free energy and enthalpies were obtained. The results showed the affinity of both drugs with water (glibenclamide: -7.96 and metformin: -11.49) due to the formation of hydrogen bonds, which were verified by the electronegativities corresponding to the dipole moment and to the partition coefficient (Log P).
Subsequently, the main properties for the design of a release system using the metformin/glibenclamide complex in the chitosan hydrogel were determined. In this process it was appreciated that the Gibbs free energy (-2157.60 kcal/mol) determined the thermodynamic stability of the adsorption. In addition, the Log P (-25.82) indicated an instantaneous solubility through the formation of hydrogen bonds and were verified by the electronic distribution and the change in dipole moment.
Keywords
Adsorption, diabetes, hydrogels, glibenclamide, metformin.Resumen
La química computacional realiza el modelado y el cálculo de propiedades fisicoquímicas que permiten comprender las diferentes interacciones moleculares a escala nanométrica en aplicaciones médicas como por ejemplo, el diseño de sistemas de liberación controlada. El modelo PM6 se utilizó para analizar metformina y glibenclamida. Primero, se obtuvieron las propiedades energéticas como la energía libre de Gibbs y las entalpías. Los resultados mostraron la afinidad de ambos fármacos con el agua (glibenclamide: -7.96 y metformina: -11.49) debido a la formación de enlaces de hidrógeno, que fueron verificados por las electronegatividades correspondientes al momento dipolar y al coeficiente de partición (Log P).
Posteriormente, se determinaron las principales propiedades para el diseño de un sistema de liberación que usa el complejo metformina/glibenclamida en el hidrogel de quitosano. En este proceso se apreció que la energía libre de Gibbs (-2157.60 kcal/mol) determinó la estabilidad termodinámica de la adsorción. Además, el Log P (-25.82) indicó una solubilidad instantánea a través de la formación de enlaces de hidrógeno y se verificó mediante la distribución electrónica y el cambio en el momento dipolar.
Palabras clave
adsorción, diabetes, hidrogeles, glibenclamida, metformina.Resumo
A química computacional realiza a modelagem e o cálculo das propriedades físico-químicas que permitem compreender as diferentes interações moleculares em escala nanométrica em aplicações médicas, como o projeto de sistemas de liberação controlada. O modelo PM6 é usado para analisar metformina e glibenclamida. Primeiro, obtenha as propriedades energéticas como a energia livre de Gibbs e as entalpias. Os resultados mostram a afinidade de ambos os componentes com água (glibenclamida: -7.96 e forma: -11.49) debitados na forma de ligações de hidrogênio, que verificam por eletronegatividades correspondentes no momento dipolar e coeficiente de participação (Log P).
Posteriormente, selecione as principais propriedades para o projeto de um sistema de liberação que usa o método completo/glibenclamida no hidrogel de quitosano. Neste processo, aprecie a energia livre de Gibbs (-2157.60 kcal/mol) que determina a estabilidade termodinâmica da adsorção. Além disso, o Log P (-25.82) indica uma solução instantânea para a passagem de forma de hidrogênio e é verificada usando a distribuição eletrônica e o câmbio no momento dipolar.
Introduction
Nanotechnology has allowed for the development of technologies for the early identification of diseases, a field where nanomedicine represents the main area of application. The applications of nanomedicine consist in the formation of images, the detection, and diagnosis of diseases as well as controlled delivery systems which increase the effectiveness of the release, the tolerability, the specificity of the drugs, and the reduction of side effects, which has increased the availability of new medical treatments [1-4]. Nanomedicine allows the study of different types of molecules such as proteins, drugs, vaccines or nucleic acids to carry out the diagnosis and treatment for cancer, diabetes and tissue engineering [5-6].
The World Health Organization (WHO) predicts a population of 500 million adults affected with type 2 diabetes by the year 2030 [7]. The design of new drug administration systems allows the penetration in the cellular tissues to carry out the glucose measurement and later perform the controlled release of antidiabetic drugs, such as glibenclamide and metformin, in complete doses. These delivery systems are synthesized using polymers, micelles, liposomes, and biodegradable polymer nanoparticles [8]. Release systems represent the main application of hydrogels, which are obtained by the dispersion in water of a polymer solution and are characterized by the degree of swelling, viscosity, pH, biodegradation, and degree of crosslinking [9].
Hydrogels are used as intelligent transporters due to their swelling capacity, release kinetics, and biocompatibility, which have generated more research at the molecular level in order to modify their structure and optimize their efficiency in pharmaceutical applications [10]. The main polymers include chitosan, poly(lactic acid-co-glycolic acid), and dextran [11].
Chitosan is a polysaccharide that allows for the release of a drug in a hydrogel due to the hydrophilic nature of functional groups such as OH, COOH, and NH. through the formation of hydrogels [8]. It also allows the adsorption of nanoparticles due to the rupture of the internal structure of enterocyte to reach the bloodstream, obtaining a controlled release system sensitive to external stimuli [11]. Chitosan hydrogels have been used in anticancer therapies as well as in controlled release systems [12-14].
Currently, computer models are used to analyze, design, and evaluate the effectiveness of new drugs with physical-chemical properties and specific biological activities, reducing the problems of effectiveness in medical treatments [15]. Thus, the pharmaceutical industry has increased the design of release systems using classical force fields or quantum mechanical models for the modeling of biological systems or the design of drugs where the calculation of Gibbs free energy is obtained for several thousand atoms simultaneously [16].
Quantitative Structure-Activity Relationship (QSAR) properties allow the calculation of non–covalent interactions, that is, they allow predicting the structure and biological behavior of macromolecules and their interactions with drugs, receptors or drug–receptors using descriptors such as thermodynamic, lipophilicity, physicochemical, electronic, and biological properties, solubility, and diffusion coefficient [3, 17-20]. The partition coefficient (Log P) is one of the most relevant QSAR properties because it provides information about the effectiveness of a drug to cross biological barriers and to reach the specific site for drug release. Log P defines the intermolecular behavior of the drugs in the human body, being widely used in the pharmaceutical sector [21].
The calculation time in DFT and AB-initio results in a higher computational cost due to the size and type of systems that are analyzed. However, the semi-empirical models of quantum mechanics simplify these times due to the incorporation of parametrizations as well as the exclusion of functions of different atoms. The PM6 model of quantum mechanics was developed for 9 000 compounds from experimental bases in conjunction with the AB–initio method, which predicts the heat of formation and the geometry of optimization more accurately because information is obtained from the interactions of atomic, hydrogen bonds as well as the use of d orbitals that allow the study of transition metals [16, 22-24].
In this work, the PM6 model was used to determine the optimization geometry, QSAR properties and the molecular electrostatic potential (MESP) in the adsorption of a metformin/glibenclamide complex on chitosan hydrogels.
Materials and Methods
All the simulations were performed in an Intel(R) Core (TM) i7 CPU @ 2.53 GHz with 4 GB RAM, having windows XP environment. The simulations involved metformin, gliblenclamide, metformin/water, and glibenclamide/water, as well as the complex (metformin/glibenclamide) where the calculation of the energetic properties and the molecular orbitals were made by selecting the PM6 semi–empirical method. Additionally, for determining the solvation energy (Aq) in the different simulations, the Cramer–Truhlar SM3 method, a SM5.4 solvation of the Spartan software was used [24].
Subsequently, using the Compute menu the QSAR properties such as surface area, volume, mass, and polar surface area were obtained. The Log P was calculated using the Spartan software by selecting the Ghose–Crippen method. This method is independent of the wave energy function, that is, the calculation is not influenced whether it is obtained using a quantum mechanics, molecular mechanics or a semi-empirical method. HOMO and LUMO molecular orbitals were obtained by solving the wave function by means of the PM6 semi–empirical method.
Tables 1 and 2 show the thermodynamic and QSAR properties of water, metformin, and glibenclamide, where a negative Gibbs free energy was obtained. Besides, the solvation energy indicates the dissolution of both drugs in water through hydrogen bonds. Electrostatic Potential maps (MESP) were calculated in the Spartan program by selecting the Hartree-Fock method with a 3-21G set. Figure 1 shows the MESP of water, metformin, and glibenclamide. It can be observed that the regions of high electronic density (in red) are presented in N-H, S=O, C=O and O-H bonds, while the regions of low electronic density are presented in the CH bonds.
Table 1.:
Properties
Units
Water
Metformin
Glibenclamide
Gibbs free energy (ΔG)
kcal/mol
- 47426
- 268288
- 1414298
Solvation energy
kcal/mol
- 12.55
- 13.44
- 28.01
Enthalpy of formation (H°
f
)
kJ/mol
- 54.27
46.46
- 191.84
Enthalpy of aqueous formation (H°
f
)
kJ/mol
- 66.82
33.01
- 219.79
Standard entropy (S°)
kJ/mol
0.450
0.094
0.167
Dipolar moment
Debye
1.86
2.99
11.23
E HOMO
eV
- 11.91
- 9.2
- 9.44
E LUMO
eV
4.070
0.35
- 0.91
Table 2.:
Properties
Units
Water
Metformin
Glibenclamide
Mass
amu
18
129.16
494.012
Log P
-------
- 0.38
- 0.50
3.80
Surface area
Ų
36.36
172.2
512.04
Polar surface area
A²
36.35
78.27
99
Polar surface
%
99.9
44
19.33
Volume
ų
19.22
137.53
471.83
Figure 1.:
Results and Discussions
Drugs/water
Tables 3 and 4 show the energetic and QSAR properties of the metformin/water and glibenclamide/water solutions at 298.15 K, where the Gibbs free energy (ΔG) indicates the spontaneous formation of both individual aqueous solutions attributed to a strong interaction of the OH bonds of the water with the receptor groups of both drugs and to the drugs polar surface originated by the presence of atoms with different electronegativity in each molecule.
The existence of an electrostatic attraction between water and both drugs caused by the rupture of the individual water ions to obtain a thermodynamic equilibrium was also observed, and then corroborated by the negative partition coefficient shown in Table 4 and the dipole moment indicating a conformational change in the structure of both solutions caused by the electronegativity of OH, CO, SO in both drugs [25-26].
Table 3.:
Properties
Units
Metformin
/water
Glibenclamide
/wáter
Gibbs free energy (ΔG)
kcal/mol
- 1429
- 1503
Enthalpy of formation (H°
f
)
kJ/mol
- 1818
- 2076
Enthalpy of aqueous formation (H°
f
)
kJ/mol
- 1324
- 1378
Standard entropy (S°)
kJ/mol
0.354
0.42
Dipolar moment
Debye
6.48
11.46
E HOMO
eV
- 9.49
- 9.37
E LUMO
eV
0.26
- 1.18
Table 4.:
Properties
Units
Metformin
/water
Glibenclamide
/water
Mass
amu
669.617
1,034.462
Log P
-------
- 11.49
- 7.96
Surface area
Ų
811.53
1,077.61
Polar surface area
A²
727.94
706.07
Polar surface
%
89.70
65.55
Volume
ų
638.55
961.47
Ovality
--------
2.26
2.29
The MESP of the dissolution of both drugs in water are depicted in Figure 2. The electronic distribution shows that the regions of high electronic density were located in the OH, C=O and S=O bonds, while the regions of low electronic density were observed in the CH bonds. Figure 2(a) metformin shows an attraction between hydrogen in water and oxygen in the OH, while in Figure 2(b) it is observed that hydrogen bonds originate between hydrogen in water and the carbonyl group. This behavior corroborates the formation of a new bond that generates the dissolution of both drugs with water.
Figure 2.:
Hydrogel/Complex (metformin/glibenclamide)
Table 5 shows the adsorption of the complex (metformin/glibenclamide) on the chitosan hydrogel because the Gibbs free energy and enthalpy indicate an exothermic process, that is, a spontaneity in the adsorption process produced by the increase of interactions in the hydrogel/complex. Also, the dipole moment reveals that the decrease of the dipole-dipole forces produces a greater solubility due to the increase of the intermolecular forces through the formation of hydrogen bonds, generating geometry with the greater structural stability of the hydrogel/complex [27-28].
The formation of hydrogen bonds is attributed to the increase in the polar surface (Table 6). Also, the existence of the affinity of the complex (metformin/glibenclamide) with the water present in the hydrogel is observed due to the negative value of the partition coefficient that indicates that the drugs are related to the polar behavior of the water. Additionally, the surface area contributes to penetrate the cell membrane and may be in contact with the bloodstream to establish a drug delivery system [29-30]. The cationic character of chitosan allows obtaining stimuli sensitive nanohydrogels to the glucose [14] due to the protonation of the amine groups that produce the adsorption of drugs in the hydrogel/complex.
Additionally, a new controlled release system is obtained since the nanometric particles of the complex have greater stability, bioavailability and longer life time, being thus more effective causing the rupture of the permeability barriers, achieving the controlled release [31].
Table 5.:
Properties
Units
Hydrogel
/Complex
Gibbs free energy (ΔG)
kcal/mol
- 2157.60
Enthalpy of formation (H°
f
)
kJ/mol
- 3784.42
Enthalpy of aqueous formation-Aq (H°
f
)
kJ/mol
- 1920.15
Standard entropy (S°)
kJ/mol
0.7945
Dipolar moment
Debye
5.97
E HOMO
eV
- 8.81
E LUMO
eV
- 0.83
Ovality
----------
2.86
Table 6.:
Properties
Units
Hydrogel
/Complex
Mass
amu
2749.16
Log P
-------
- 25.82
Surface area
Ų
2545.02
Polar surface area
A²
1348.59
Polar surface
%
52.96
Volume
ų
2494.83
Polarizability
-------
39.49
Figure 3 shows the map of the electronic distribution of the adsorption of the drugs (glibenclamide and metformin) in the chitosan hydrogel, where the deformations present in the rings of glibenclamide are attributed to the functional groups OH and NH of the chitosan. Additionally, the regions of high electronic density zones correspond to the polar groups on the surface of the hydrogel, attributing the adsorption of both drugs to the formation of hydrogen bonds. In addition, in the sulphonylurea group of glibenclamide, a positive electronic density is observed, while in metformin, ionic interactions producing electrostatic forces that corroborate the transfer of charge with the hydrogel are observed [32-34].
Figure 3.:
Conclusions
The computational analysis showed that there is a process of adsorption of metformine/glibenclamide complex using chitosan hydrogels due to the thermodynamic values as a Gibbs free energy of -2157.60 kcal/mol. The negative value of Gibbs free energy determined a spontaneous reaction which corresponded to the formation of hydrogen bonds. This was corroborated by the negative value of the coefficient partition (Log P) of -25.82, which indicated the hydrophilic capacity increases in the adsorption and the dipole moment change as a product of a greater electronegativity attributed to hydrogel.
Moreover, the surface area indicated that both drugs are dispersed molecularly, that is, there would be no problems of crystallinity in the hydrogel attributable to the interactions and the size of the pores in the hydrogel that produce thermodynamic stability. Finally, through the enthalpy of formation and the enthalpy of aqueous formation it was possible to determine that glibenclamide has a greater affinity for hydrogel due to the volume occupied by the molecule in the hydrogel. The OH and SO functional groups of glibenclamide increased the water solubility of the hydrogel as well as the presence of secondary bonds such as the van der Waals forces.
References
Referencias
R. Gupta, J. Biotechnol. Biomater. vol. 7, n° 3, pp. 1-2, 2017. DOI: https://doi.org/10.4172/2155-952X.1000268
B. B. Kadhim, H. O. Muhsen, Nanosci. Nanotechnol. vol. 7, n° 1, pp. 9-13, 2017. DOI: https://doi.org/ 10.5923/j.nn.20170701.03
N. Haddish-Berhane, J. L. Rickus, K. Haghigh, Int. J. Nanomed. vol. 2, n°3, pp. 315–331, 2007.
P. S. Appasaheb, S. D. Manohar, S. R. Bhanudas, J. Adv. Pharm. Edu. & Res. vol. 3, n° 4, pp. 333-346, 2013. DOI: https://doi.org/10.5958/2321-5844.2018.00004.3
G. Sharma, A. R. Sharma, J. S. Nam, G. P. C. Doss, S. S. Lee, C. Chakraborty, J. Nanobiotechnol. vol. 13, pp. 1-13, 2015. DOI: https://doi.org/10.1186/s12951-015-0136-y
R. M. Di-Santo, V. Subramanian, Z. Gu, ” Wiley Interdiscip. Rev. Nanomed. Nanobiotechnol. vol. 7, pp. 548–564, 2015. DOI: https://doi.org/10.1002/wnan.1329
O. Veiseh, B. C. Tang, K. A. Whitehead, D. G. Anderson, R. Langer, Nat. Rev. Drug Discov. vol. 14, pp. 45–57, 2015. DOI: https://doi.org/10.1038/nrd4477
C. S. Satish, K. P. Satish, H. G. Shivakumar, Indian J. Pharma. Sci. vol. 68, pp. 133-140, 2006. DOI: https://doi.org/10.4103/0250-474X.25706
I. S. Raja, N. N. Fathima, Springerplus. vol. 3, pp. 393-402, 2014. DOI: https://doi.org/ 10.1186/2193-1801-3-393
P. A. Arredondo, M. E. Londoño, Rev. Ing. Biomédica. vol. 3, n° pp. 83-94, 2009. DOI: https://doi.org/10.24050/19099762.n5.2009.70
V. R. Devadasu, T. M. Alshammari, M. Aljofan, Int. J. Diabetes Dev. Ctries. vol. 38, n° 1, pp. 11-19, 2018. DOI: https://doi.org/10.1007/s13410-017-0558-1
C. Dalwadi, P. Gayatri, Rec. Pat. Nanotech. vol. 9, n° 1, pp. 17-25, 2015. DOI: https://doi.org/ 10.2174/1872210509666150101151521
S. K. Basu, A. Rajendran, Chem. Pharm. Bull. vol. 56, n° 8, pp. 1077-1084, 2008. DOI: https://doi.org/10.1248/cpb.56.1077
L. J. Del-Valle, A. Díaz, J. Puiggalí, Gels. Vol. 3, n° 3, 2017. DOI: https:/doi.org/10.3390/gels3030027
R. Galeazzi, Curr. Comput-Aid Drug. vol. 5, pp. 225-240, 2009. DOI: https:/doi.org/10.2174/157340909789577847
A. V. Sulimov, D. C. Kutov, E. V. Katkova, V. B. Sulimov, Adv. Bioinformat. vol. 2017, n° 5, pp. 1-6, 2017. DOI: https:/doi.org/10.1155/2017/7167691
J. L. Medina-Franco, F. López-Vallejo, R. Castillo, Educ. Química. vol. 17, pp. 452-457, 2006. DOI: http://dx.doi.org/10.22201/fq.18708404e.2006.4.66027
J. Liu, H. Zheng, X. Dai, S. Sun, H. G. Machens, A. F. Schilling, J. Tissue Sci. Eng. vol. 8, pp. 193, 2017. DOI: https:/doi.org/10.4172/2157-7552.1000193
L. G. Ferreira, R. N. Dos-Santos, G. Oliva, A. D. Andricopulo, Molecules. vol. 20, pp. 13384-13421, 2015. DOI: https:/doi.org/10.3390/molecules200713384
M. V. Vellayappan, J. R. Venugopal, S. Ramakrishna, S. Ray, A. F. Ismail, M. Mandal, A. Manikandan, S. Seal, S. K. Jaganathan, RSC Adv. vol. 6, pp. 83638–83655, 2016. DOI: https:/doi.org/ 10.1039/C6RA15252J
A. Kovalenko, Condens. Matter. Phys. vol. 18, n° 3, pp. 1-24, 2015. DOI: https://doi.org/10.5488/CMP.18.32601
J. J. P- Stewart, J. Mol. Model. vol. 13, n° 12, pp. 1173-1213, 2007. DOI: https://doi.org/10.1007/s00894-007-0233-4
J. J. P. Stewart, J. Mol. Model. vol. 14,n° 6, pp. 499–535, 2008. DOI: https://doi.org/10.1007/s00894-008-0299-7
L. O. Olasunkanmi, J. Ige, G. O. Ogunlusi, J. Chem. vol. 2013, pp. 1-7, 2013. DOI: https://doi.org/10.1155/2013/640649
G. Schaftenaar, J. Vlieg, J. Comput. Aid. Mol. Des. vol. 26, n° 3, pp. 311–318, 2012. DOI: https://doi.org/10.1007/s10822-012-9557-y
C. W. Pouton, Eur. J. Pharm. Sci. vol. 29, n° 3-4, pp. 278-287, 2006. DOI: https://doi.org/ 10.1016/j.ejps.2006.04.016
E. Mah, R. Ghosh, Processes. vol. 1, pp. 238-262, 2013. DOI: https://doi.org/10.3390/pr1030238
Q. Chai, Y. Jiao, X. Yu, Gels. 2017. DOI: https://doi.org/10.3390/gels3010006
C. F. Santa, B. L. López-Osorio, Rev. Acad. Colomb. Cienc. Exact. Fis. Nat. vol. 37, pp. 115-124, 2013.
S. Hurst, C. M. Loi, J. Brodfuehrer, A. El-Kattan, Expert. Opin. Drug Metab. Toxicol. vol. 3, n° 4, pp. 469-489, 2007. DOI: https://doi.org/10.1517/17425225.3.4.469
S. Jana, S. Maiti, in Nanostructures for oral medicine: Nanostructures in therapeutic medicine series, E. Andronescu, A. M. Grumezescu, Edits., Elsevier, 2017, pp. 607-638.
I. Amăriuţei, M. Florea, M. V. Ghica, L. Popa, Romanian J. Biophys. vol. 25, n° 3, pp. 173–184, 2015.
D. Horter, J. B. Dressman, Adv. Drug. Deliv. Rev. vol. 46, pp. 75-87, 2001. DOI: https://doi.org/10.1016/s0169-409x(00)00130-7
R. Galeazzi, Curr. Comput-Aid Drug, vol. 5, pp. 225-240, 2009. DOI: https://doi.org/10.2174/157340909789577847
Cómo citar
IEEE
ACM
ACS
APA
ABNT
Chicago
Harvard
MLA
Turabian
Vancouver
Descargar cita
CrossRef Cited-by
1. Marisol Ibarra-Rodríguez, Víctor Víctor Jiménez Pérez, Blanca M. Muñoz-Flores, Mario Sánchez. (2025). Metformin Adsorption on Binuclear Boron Schiff Complexes. Journal of the Mexican Chemical Society, , p.652. https://doi.org/10.29356/jmcs.v69i3.2207.
2. Olawale Folorunso Akinyele, Emmanuel Gabriel Fakola, Omolara Olubunmi Adeboye, Sunday Chimela Chinuomah. (2023). Investigation of Methyl-5-(pentan-3-yloxy)-7-oxabicyclo[4.1.0]hept-3-ene-3-carboxyhydrazide Derivatives as Potential Inhibitors of COVID-19 Main Protease: DFT and Molecular Docking Study. Bioinformatics and Biology Insights, 17 https://doi.org/10.1177/11779322231182050.
Dimensions
PlumX
Visitas a la página del resumen del artículo
Descargas
Licencia
Derechos de autor 2020 Revista Colombiana de Química

Esta obra está bajo una licencia internacional Creative Commons Atribución 4.0.
Los autores/as conservarán sus derechos de autor y garantizarán a la revista el derecho de primera publicación de su obra, el cuál estará simultáneamente sujeto a la Licencia de reconocimiento de Creative Commons (CC. Atribución 4.0) que permite a terceros compartir la obra siempre que se indique su autor y su primera publicación en esta revista.
Los autores/as podrán adoptar otros acuerdos de licencia no exclusiva de distribución de la versión de la obra publicada (p. ej.: depositarla en un archivo telemático institucional o publicarla en un volumen monográfico) siempre que se indique la publicación inicial en esta revista.
Se permite y recomienda a los autores/as difundir su obra a través de Internet (p. ej.: en archivos telemáticos institucionales o en su página web) antes y durante el proceso de envío, lo cual puede producir intercambios interesantes y aumentar las citas de la obra publicada. (Véase El efecto del acceso abierto).



