JBCS

INTRODUCCION

Las ácidos húmicos (HA) aumentan la fertilidad, la retención de minerales, la acidez del agua adsorbida y el transporte de contaminantes a través del suelo, también afectan la química de la atmósfera a través del ciclo del carbono.^1-3 La gran mayoría de los estudios realizados con sustancias húmicas y fúlvicas se han enfocado en suelos⁴ debido a la propiedad que tienen estas sustancias para recuperar suelos erosionados e infértiles; otro aspecto positivo son sus propiedades ambientales, bioquímicas, y terapéuticas. Aunque existen discrepancias entre las estructuras propuestas para los HA, estas poseen semejanzas en cuanto al número y tipo de grupos funcionales presentes.⁵

Los estudios reportados sobre sustancias húmicas son relativamente escasos, debido a su gran diversidad y ambigüedad estructural. Asimismo, las interpretaciones de las propiedades químicas teóricas de dichas sustancias han sido muy sesgadas en cuanto a predicción de su utilidad e interacción con los principales cationes metálicos de interés vital. Entre los estudios más sobresalientes encontramos: En el año 2007 a Troul y colaboradores, los cuales estudiaron las interacciones de acidos fulvicos con benceno y Al (III).⁶ En el mismo año, De Castro y colaboradores reportaron un estudio teórico de la complejación de Zn (II) y Cd (II) con ácidos fulvicos.⁷ En el 2008, André y colaboradores estudiaron la regioselectividad de Al (III) frente a sustancias con grupos catecol.⁸ En el año 2009, Burk y colaboradores estudiaron las interacciones de cesio con sustancias orgánicas de naturaleza aniónica provenientes del suelo.⁹ En el año 2014, Aquino y colaboradores reportaron las regiones de transferencia protónica de ácidos húmicos.¹⁰ Recientemente, en el año 2017, Kimura y colaboradores estudiaron el efecto de agregación de suelos y su relación con ácido húmicos como indicadores de restauración de un bosque ribereño.¹¹ En el mismo año, Xu y sus colaboradores estudiaron las estructuras de los acidos húmicos presente en dos suelos típicos de China.¹²

En la presente investigación se estudiaron las interacciones de dos modelos de ácidos húmicos con los cationes Al³⁺, Ca²⁺, Mg²⁺, Zn²⁺, K⁺ y NH4⁺. El estudio teórico se realizó bajo el marco conceptual de la Teoría del Funcional de la Densidad (DFT) con el fin de conocer la naturaleza de las interacciones de los ácidos húmicos con distintos cationes y su afinidad. Los modelos de HA: fueron el TNB (Temple-Northeastern-Birmingham)¹³ y el de Kolla^1,3 (ver Figura 1); se utilizaron descriptores de reactividad global como el potencial químico, la electronegatividad y la dureza; y descriptores locales como la suavidad química y las Funciones de Fukui^14-18 para describir su comportamiento.

Figura 1 Estructura de los ácidos húmicos según el modelo de TNB y Kolla

METODOLOGIA

Para optimizar la geometría de los HA y sus complejos metálicos se usó la siguiente metodología computacional:

Primero se optimizaron las estructuras HA a un nivel de cálculo semiempírico PM6. Posteriormente, estas geometrías fueron corridas a un nivel DFT utilizando el funcional B3LYP y la función de base 6-31+G (d),¹⁷ para simular la solvatación se usó el modelo continuo polarizable PCM. Las estructuras optimizadas de HA fueron usadas para construir los posibles complejos, los cuales se construyeron teniendo en cuenta los sitios de afinidad electrónica de acuerdo a los descriptores moleculares locales. Por último, se optimizaron las estructuras de los complejos formados entre los sitios activos de HA (sHA) y los cationes mencionados anteriormente, utilizando niveles de cálculos PM6 y DFT/LanL2dz iniciando en una distancia de enlace M-sHA de 2 Å.^11,19 La primera optimización de la geometría de los complejos se realizó a nivel PM6.¹³ Luego, esta geometría se tomó como punto de partida para el cálculo de la energía a un nivel DFT/LanL2dz.

Finalmente, los resultados de los diferentes optimizaciones se utilizaron para calcular los descriptores de reactividad química tales como. Suavidad, dureza y funciones de Fukui para los sectores de afinidad electrónica. Todas las geometrías iniciales fueron dibujadas usando el programa Gaussview y los cálculos descritos anteriormente fueron realizados usando el programa Gaussian 09.^19,20

RESULTADOS Y DISCUSION

Optimización de las estructuras

Las estructuras iniciales de TNB y Kolla fueron optimizadas exitosamente y no se observaron frecuencias negativas, lo cual asegura que son las estructuras más estables o con la energía más baja. También se observaron varios enlaces de hidrogeno intramolecular, lo cual le da mayor estabilidad y menor energía al sistema molecular HA. Las Figuras 2 y 3 muestran las estructuras, correspondiente a los estereoisómeros RSSRSSR para el modelo TNB y RRRSRSR para el modelo de Kolla. Los isómeros de menor energía fueron los SSRSRSR, para Kolla y RRRSRSR para TNB.

Figura 2. Enlaces de hidrógeno intramolecular y las tres interacciones tipo Van der Waals O---H-C en el stereoisomero RSSRSSR del modelo TNB (Temple-Northeastern-Birmingham)

Figura 3. Enlaces de hidrógeno intramolecular y las tres interacciones tipo Van der Waals O---H-C en el estereoisomero RRRSRSR del modelo de Kolla

En las estructuras optimizadas de ambos isómeros se observaron dos diferencias marcadas, las cuales fueron el número de enlaces de hidrogeno y la posición de la región acido carboxílico, lo que le da el carácter acido a las moléculas HA. Mientras que el modelo TNB presentó un enlace de hidrógeno intramolecular y dos interacciones tipo Van der Waals entre los átomos O---H-C, la estructura de Kolla mostro dos enlaces de hidrógeno intramolecular y tres interacciones tipo Van der Waals entre los átomos O---H-C., por lo que se puede considerar una estructura más estable y menos susceptible a transportar cationes.

Análisis de los sitios de interacción con cationes

Una vez optimizadas los diferentes isómeros, se usaron los de menor energía para estudiar la interacción con los cationes de mayor importancia en los procesos que ocurren en el suelo, entre los cuales están Al³⁺, Ca²⁺, Mg²⁺, Zn²⁺, K⁺ y NH₄⁺.

Siguiendo el código de colores que muestra el mapa de potencial electrostático (Figuras 4 y 5), las regiones de mayor afinidad nucleófilica o los centros de mayor probabilidad de interacción con cationes para el modelo de Kolla son los átomos de oxigeno 2, 3, 12, 13, 17 y 18 y los átomos de nitrógeno 4 y 16; mientras que para el modelo TNB son los átomos de oxigeno 2, 3, 8, 9, 11, 12, 13, 17 y 18 y los nitrógenos 4 y 16. La ilustración cualitativa de los MEP y los indicadores de las funciones de Fukui para ataques electrófilicos (f+) permitió estimar el punto de partida para simular los diferentes complejos catión-HA. Teniendo en cuenta que valores altos de estos descriptores moleculares indican una alta reactividad. El análisis de los descriptores también mostró que el modelo TNB posee mayor número de sitios activos para interaccionar con cationes. Por su parte los descriptores moleculares suavidad local y funciones de Fukui para un ataque electrófilico mostraron concordancia a lo previsto en los MEP's (ver Tabla 1)

Figura 4. MEP para los átomos más representativos y susceptibles a la interacción con los cationes metálicos del modelo de Kolla

Figura 5. MEP para los átomos más representativos y susceptibles a la interacción con los cationes metálicos del modelo TNB

De acuerdo a los átomos de mayor afinidad electrónica y la estructura global de las modelos moleculares, se pudieron establecer las zonas de mayor probabilidad a interaccionar con cationes. Las huecos nucleofilicos quedaron distribuidas en cuatro zonas desde la A hasta la D. (ver Figura 6). De acuerdo a las zonas se simularon los distintos acercamientos catiónicos.

Figura 6. Regiones de afinidad nucleófilica de acuerdo a los indicadores de reactividad local en las estructuras del modelo (TNB) y (Kolla) respectivamente. Las letras A, B, C y D indican los sectores que presentaron mayor probabilidad de interacción frente a un catión

Específicamente, para el modelo TNB, en la región A se encuentran los átomos O1, O2, O3 y N4; los átomos con mayor valor en las Funciones de Fukui, N4, O2, y O3 (Tabla 1 y 2) tienen la mayor posibilidad de recibir un ataque electrófilico los cationes. El átomo O1 tiene menos posibilidades de interactuar con los cationes debido a su bajo valor en las Funciones de Fukui. En la región B, la cual presentó mayor reactividad se encuentran los átomos O5, O6 y O7; el O6,. En la región C los átomos, O8, O9 y O11, los cuales mostraron un valor muy parecido en las funciones de Fukui lo que hace a esta región muy reactiva. En la región D, los átomos N16, O12, O13, O17 y O18 mostraron, en promedio, los valores más altos en las Funciones de Fukui; esto se debe a que estos átomos se encuentran libres de interacciones o repulsiones con otros átomos o grupos funcionales de la estructura por lo que no comparten su carga electrónica, además los grupos hidroxilo y amino son excelentes nucleófilos. Sin embargo, para predecir que átomo es más propenso a sufrir un ataque electrófilico deben usarse en conjunto las Funciones de Fukui y el potencial electrostático molecular.

Para el modelo de Kolla, en la región A los átomos con mayor valor en las Funciones de Fukui fueron: N4, O2 y O3; en la región B el átomo con mayor valor en las Funciones de Fukui fue el O6, en la región C no hay átomos que presenten un valor significativo en este parámetro, y en la región D los átomos con mayor valor fueron N16, O17, O12, O13 y O18. Altos valores de suavidad indican que los átomos son susceptibles a la deformación de la nube electrónica o son fácilmente polarizable. Los átomos con mayores valores de suavidad local y función de Fukui fueron los mencionados arriba, ya que estas funciones se relacionan directamente (Figura 6).

Por otro lado, se observó que los grupos fenólicos, hidroxilos, carboxilatos, carbonilos y amino son más reactivos en medio acuoso que en el vacío.

Generalmente, los valores de dureza fueron bajos y los de suavidad altos,^21-23 lo cual las hace estructuras susceptibles a la deformación de la nube electrónica tendiendo así a la interacción con un agente electrófilico o a la disociación, por lo que estas estructuras podrían interactuar fácilmente con los metales y otros aceptores de electrones; ambos modelos presentaron valores muy bajos de electrofilicidad y mucho más el de Kolla, por lo que este puede tener más afinidad por los agentes electrófilicos que el modelo TNB.

Por su parte, los cationes tendieron a ser blandos (Tabla 3) excepto el aluminio. Según el principio dureza-suavidad acido-base (HSAB), las sustancias blandas como los cationes reportados prefieren unirse a los grupos carboxilatos, los cuales también son blandos.²⁴ Igualmente sustancias duras como el aluminio se unen con átomos duros como los oxígenos de los grupos hidroxilos. Los valores bajos de electronegatividad para los ácidos húmicos y altos para los cationes nos indican que hay transferencias de carga desde los ácidos húmicos hacia los cationes. Igualmente, se observaron valores de electrofilicidad bajos de los ácidos húmicos y altos para los cationes, lo cual mostró la dirección del movimiento de las cargas.

Análisis energético de los ácidos húmicos por los cationes

Para determinar la afinidad de los modelos HA de Kolla y TNB por cationes de distinta naturaleza se evaluó la estabilidad de los complejos formados entre los cationes Al³⁺, Ca²⁺, Mg²⁺, Zn²⁺, K⁺ y NH₄⁺ y cada cavidad nucleófilica de los isómeros RRRSRSR y RSSRSSR (Figuras 2 y 3). Para este análisis se calculó la energía siguiendo la reacción 1, las cuales están presentadas en la Tabla 4:

Las energías de interacción (δE) entre HA y los respectivos cationes a un nivel de cálculo DFT B3LYP/LANL2DZ al vacío mostraron valores negativos indicando una fuerte afinidad entre ellos, sugiriendo además que dichas reacciones son espontáneas a condiciones normales. En algunas regiones específicas los cationes mostraron mayor afinidad por los modelos AH; por ejemplo, con el modelo de Kolla, los Al³⁺, Mg²⁺ y Zn²⁺ mostraron mayor afinidad en las regiones A y D de la estructura, ya que las regiones B y C presentaron mayor conectividad por enlaces de hidrógeno. Consecuentemente, los descriptores moleculares indicaron que la mayoría de los átomos involucrados en la interacción con estos cationes provienen de grupos hidroxilos, unos pocos carbonilos y aminos, y muy pocos carboxilos.

El zinc mostró un comportamiento especial al ser el más fuertemente atraído por el modelo de Kolla en todas sus regiones, como lo indican las energías de interacción con los distintos átomos al vacío

Para el modelo TNB al vacío, los de mayor energía de interacción en promedio fueron Al³⁺, Mg²⁺ y Ca²⁺. Las regiones del estereoisomero que mostraron mayor afinidad por el Al³⁺ fueron las regiones A y D pero para Mg²⁺ y Ca²⁺ fueron las regiones B y C.

En términos generales, ambos modelos presentaron fuerte afinidad por los cationes metálicos. El catión aluminio mostró mayor energía de interacción en ambos modelos, debido que este se compleja de manera tridentada con los estereoisomeros. Los cationes con dos cargas también presentaron valores altos en las energías de interacción. Los cationes monovalentes, debido a sus interacciones monodentadas, presentaron los valores más bajos (ver Tablas 4 y 5):

Con el fin de corroborar los resultados arrojados en vacío se calcularon las energías de los complejos y de interacción de la reacción en este medio acuoso. Siguiendo el modelo de reacción siguiente:

En términos generales, debido a que el solvente estabiliza las cargas, las energías de interacción en ambos modelos fueron más bajas que al vacío, aunque estas también indicaron procesos espontáneos. Las energías de interacción de los cationes con los grupos hidroxilo, amino y carbonilo siguieron la tendencia siendo altas y la de los carboxilatos un poco más bajas, lo cual puede deberse a la capacidad de este grupo de neutralizar la carga eléctrica remanente de un átomo.¹³ Siguiendo la misma tendencia, las energías de interacción indicaron que, los complejos con el Al³⁺ fueron más estables que los de los demás cationes debido a la mayor carga de este catión (ver Tablas 6 y 7).

CONCLUSIONES

Globalmente los sitios más reactivos fueron los oxígenos de los grupos hidroxilos, fenólicos, carbonilos y los nitrógenos tanto al vacío como en medio acuoso, especialmente los grupos carbonilo. Los isómeros reportados para los modelos de TNB y de Kolla mostraron afinidad por los cationes de diferentes estados de oxidación. Se notó influencia de los enlaces de hidrogeno en la interacción con el catión en alguna cavidades nucleófilica o sitios reactivos. Los MEP's, la suavidad local, las funciones de Fukui y el principio HSAB describieron adecuadamente los sitios o cavidades reactivas y las interacciones complejantes HA-catión guardando relación en ambos medios. Los resultados encontrados muestran las propiedades quelantes de los HA's y su función en el suelo, lo cual los hacen atractivos como componentes influyentes en la actividad agrícola.

REFERENCIAS BIBLIOGRAFICAS

1. Al-Faiyz, Y. S. S.; Arab. J. Chem. 2012, doi.org/10.1016/j.arabjc. 2012.12.018.

2. Humus chemistry: genesis, composition and reactions. En: Stevenson F.J., editor. New York: Wiley, 1994.

3. Senesi, N.; Soil Sci. 1999, 164, 841.

4. Jiménez Montoya, J. A.; Revista Científica 1999, 1, 143.

5. Le Person, A.; Moncomble, A.; Cornard, J. P.; J. Phys. Chem. A 2014, 118, 2646.

6. Trout, C.C.; Kubicki, J. D.; Geochim. Cosmochim. Acta 2007, 71, 3859.

7. de Castro, T. R.; da Cunha, F. F. E.; Bicca, R.; Espínola, A.; Water Air Soil Pollut. 2007, 183, 467.

8. André, E.; Lapouge, C.; Cornard, J. P.; J. Phys. Chem. A 2008, 112, 9829.

9. Burk, P.; Tammiku-Taul, J.; Tamp, S.; Sikk, L.; Sillar, K.; Mayeux, C.; Gal, J. F.; Maria, J. P.; J. Phys. Chem. A 2009, 113, 10734.

10. Aquino, A. J. A.; Tunega, D.; Schaumann, G. E.; Haberhauer, G.; Gerzabek, M. H.; Lischka, H.; Geoderma 2014, 213, 115.

11. Kimura, A.; Baptista, M. B.; Scotti, M. R.; Ecological Engineering 2017, 98, 146.

12. Xu, J.; Zhao, B.; Chu, W.; Mao, J.; Zhang, J.; Ecological Engineering 2017, 576, 444.

13. Vivas-Reyes, R.; Arguello, E.; Ayazo. H.; Revista de la Real Academia de Ciencias de Zaragoza 2007, 62, 99.

14. Geerlings, P.; De Proft, F.; Langenaeker, W.; Density Functional Theory: A bridge between Chemistry and Physics, VUB Press: Brussel, 1998.

15. Jassen, F.; Introduction to Computational Chemistry, Wiley: Hoboken, 2006.

16. Szabo, A.; Ostlund, N. S.; Modern Quantum Chemistry. Introduction to Advanced Electronic Structure Theory. Dover Publications: Mineola, 1996.

17. Hehre, W. J.; Radom, L.; Schleyer P. V. R.; Pople, J.; Ab initio molecular orbital Theory, Wiley: New York, 1986.

18. Foresman, J. B.; Frisch, A.; Exploring Chemistry with Electronic Structure Methods, 2^nd ed., Gaussian Inc., 1996.

19. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, Ö.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J.; Gaussian 09, Revision B.01, Gaussian, Inc., Wallingford CT, 2009.

20. Kolla, S.; Paciolla, M. D.; Sein, L. T.; Moyer, J.; Walia, D.; Heaton H.; Jansen, S. A. In Humic Substances: Structures, Properties and Uses; Ghabbour, E.A., Davies G., eds.; Royal Society of Chemistry: Cambridge, 1998. pp. 215-225.

21. Vivas-Reyes, R.; Zapata, J.; Rev. Colomb. Quim. 2008, 37, 145.

22. Vivas-Reyes, R.; Espinosa-Fuentes, E.; Forigua, J.; Arias, H.; Gaitán, R.; Arguello, E.; J. Mol. Struct.: THEOCHEM 2008, 862, 92.

23. Pumachagua, R.; Pecho, R. H.; Pino, R. H.; Nagles, E. O.; Hurtado, J. J.; Rev. Soc. Quim. Peru 2010, 76, 92.

24. Chattaraj, P. K.; Lee, H.; Parr, R. G.; J. Am. Chem. Soc. 1991, 113, 1855.