SciELO - Scientific Electronic Library Online

 
vol.64 número3Phycocyanin Thermo-photostability: An Accelerated Life-test AnalysisConceptual DFT Reactivity Descriptors Computational Study of Graphene and Derivatives Flakes: Doped Graphene, Graphane, Fluorographene, Graphene Oxide, Graphyne, and Graphdiyne índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • No hay artículos similaresSimilares en SciELO

Compartir


Journal of the Mexican Chemical Society

versión impresa ISSN 1870-249X

J. Mex. Chem. Soc vol.64 no.3 Ciudad de México jul./sep. 2020  Epub 04-Jun-2021

https://doi.org/10.29356/jmcs.v64i3.1213 

Articles

Cloroquine and Hydroxychloroquine: The Yin-yang of these Drugs from a Theoretical Study

Ana Martínez1  * 

1Departamento de Materiales de Baja Dimensionalidad, Instituto de Investigaciones en Materiales, Universidad Nacional Autónoma de México, Circuito Exterior SN, Ciudad Universitaria, CP 04510, CDMX, México.


Abstract

On November 12, 2019 the first case of infection of a novel coronavirus (COVID-19) was identified at Wuhan, the capital of Hubei province in China. Regardless of the origin of this virus, which is uncertain, it has produced a pandemic that has been the cause of human deaths worldwide. Two drugs are being used as antiviral against this virus; cloroquine and hydroxychloroquine, although controversy exists concerning their effectiveness. The main objective of this report is to present the electronic properties of cloroquine and hydroxychloroquine that might help explain the action mechanisms against virus. The idea that emerges from this study is that acid-base equilibrium is not the only criteria of importance to explain the action mechanism, but that the oxide-reduction balance may also help explain the toxicity or effectiveness of these drugs. Being molecules able to oxidize other molecules is similar to yin-yang; a dualism that describes contrary forces, as oxidation may produce dysfunction and affect the conditions needed for viral infection, replication and propagation of the virus, but also contribute to increasing oxidative stress. These results offer a further step along the path of understanding these action mechanisms.

Keywords: COVID-19; SARS-CoV2; coronavirus; oxidation; electroaccepting power

Resumen

El 12 de noviembre de 2019 se identificó el primer caso de infección de un nuevo coronavirus (COVID-19) en Wuhan, la capital de la provincia de Hubei en China. Independientemente del origen de este virus, que es incierto, ha producido una pandemia que ha sido la causa de muertes humanas en todo el mundo. Se están usando dos medicamentos como antivirales contra este virus; cloroquina e hidroxicloroquina, aunque existe controversia sobre su efectividad. El objetivo principal de este trabajo es presentar las propiedades electrónicas de la cloroquina y la hidroxicloroquina que podrían ayudar a explicar los mecanismos de acción contra los virus. La idea que surge de este estudio es que el equilibrio ácido-base no es el único criterio de importancia para explicar el mecanismo de acción, sino que el equilibrio óxido-reducción también puede ayudar a explicar la toxicidad o efectividad de estos medicamentos. Ser moléculas capaces de oxidar otras moléculas es similar al yin-yang; un dualismo que describe fuerzas contrarias, ya que la oxidación puede producir disfunción y afectar las condiciones necesarias para la infección viral, la replicación y la propagación del virus, pero también contribuye a aumentar el estrés oxidativo. Estos resultados ofrecen un paso más en el camino hacia la comprensión del mecanismo de acción.

Palabras clave: COVID-19; SARS-CoV2; coronavirus; oxidatción; poder electroaceptor

Introduction

On November 12, 2019, the first case of infection of a novel coronavirus (COVID-19) was identified at Wuhan, the capital of Hubei province in China. [1] Regardless of the origin of this virus, which is uncertain, it has produced a pandemic that has been the cause of human deaths worldwide. Governments have attempted to control the pandemic by requesting that the population maintain a healthy distance and stay at home. This produces collateral damage, as it affects the economy of the population. Social inequality in Mexico is particularly shameful. Fifty percent of the population lives in poverty and survives on what they can earn each day. For these people, staying at home has many repercussions, making it difficult to stop the pandemic solely with these actions (social distancing and remaining at home). Therefore, it is very important to develop the vaccine and/or find drugs to help control this infection.

The Coronaviridae family is constructed from enveloped viruses with a positive-strand RNA genome. [2,3] These coronaviruses mostly contain four structural proteins: spike protein, envelope protein, membrane protein, and nucleocapsid protein. [4] Spike proteins are relevant during virus infection, as they promote host attachment and virus-cell membrane fusion. [5] COVID-19 is a coronavirus that belongs to this family. Several compounds have been investigated for the purpose of developing antiviral drugs or a vaccine against human COVID-19 infection, but an effective and definitive antiviral strategy to halt this emergency has still not been devised.

Two drugs that have been used for this purpose are cloroquine and hydroxychloroquine. First synthetized in 1934, cloroquine is a medicine for the treatment and prevention of malaria and autoimmune conditions (rheumatoid arthritis and lupus erythematosus for example). Hydroxychloroquine was produced later in 1955, and was reported to be as effective as cloroquine against malaria, but less toxic. It was later stated that based on its physicochemical properties, hydroxychloroquine was much less active than cloroquine. Both were studied for the treatment of SARS-CoV, which is the virus responsible for the SARS coronavirus epidemic that we suffered in 2002. [6-9]

Up until now, it has not been clear whether these drugs are effective for the treatment of the COVID-19 infection. [9-11] French studies reported that hydroxychloroquine combined with azithromycin might reduce the viral load in patients, [12,13] but other studies reported that there was insufficient evidence to verify the effectiveness of these drugs against COVID-19. [14] However, in previous studies authors suggested that the lack of effectiveness of these drugs may relate to other patient conditions, such as renal or hepatic damage. Although there are many clinical studies, the action mechanism of chloroquine and hydroxychloroquine against COVID-19 is not yet fully understood. One of the hypotheses is that both drugs can be incorporated into endosomes and lysosomes, producing an increment in the pH of the intracellular compartments. This produces dysfunction and affects the conditions required for viral infection, replication and propagation of the virus. [15]

In order to understand the action mechanism, it is important to have accurate information about any molecules that might be used in the treatment of COVID-19. In spite of all the studies reported until now, there are no theoretical studies concerning the electronic properties of these drugs. For this reason, the main objective of this report is to present electronic properties of chloroquine and hydroxychloroquine that may help elucidate the action mechanisms against this virus. Even though now it is known that these two molecules do not have positive effect on the disease produced by the virus and apparently they are not longer useful for the treatment of COVID-19, it is important to have an electronic characterization of these drugs, since they were used with success with several patients (Plaquenil ® was used in Mexico for patients with COVID-19) and they were also useful for the treatment of other illness.

Computational details

Gaussian09 was used for all electronic calculations. [16] Geometry optimizations without symmetry constraints were implemented at M062x/6-311+G(2d,p) level of theory [17-21] while applying the continuum solvation model density (SMD) with water in order to mimic polar environments. [22] Harmonic analyses were calculated to verify local minima (zero imaginary frequencies). We considered protonated states of all drugs following the available experimental evidence. PubChem configurations were used as initial geometries.

Conceptual density functional theory (CDFT) is a chemical reactivity theory founded on Density Functional Theory based concepts. [23-27] Within this theory there are response functions that help us to understand the chemicial reactivity. The response functions that we used in this investigation are the electro-donating (ω-) and electro-accepting (ω+) powers, previously reported by Gázquez et al. [28, 29] These authors defined the propensity to donate charge or ω- as follows:

ω- = (3I + A)2 / 16 (I-A) (1)

whereas the propensity to accept charge or ω+ is defined as

ω+ = (I + 3A)2 / 16 (I-A) (2)

I and A are vertical ionization energy and vertical electron affinity, respectively. Lower values of ω- imply greater capacity for donating charge. Higher values of ω+ imply greater capacity for accepting charge. In contrast to I and A, ω- and ω+ refer to charge transfers, not necessarily from one electron. This definition is based on a simple charge transfer model expressed in terms of chemical potential and hardness. These chemical descriptors have been used successfully in many different chemical systems. [30-34]

Results and discussion

Cloroquine and hydroxychloroquine represent diprotic weak bases that are protonated at physiological conditions, depending on pKa values. It was reported that these drugs should be administered in tablets in its basic form. For this reason, it is important to calculate neutral and protonated molecules. Fig. 1 reports the optimized structures of cloroquine (neutral and protonated). Bond distances between the nitrogen atom of pyridine and a carbon atom of one methyl group are included in each case. This indirectly indicates the doubling of the molecule. Evidently, the di-protonated molecule is less folded than the neutral and the bond distance of the mono-protonated structure is the shortest. For hydroxychloroquine, the optimized structures (see Figure 2) indicate that molecules are also folded. The distance between the nitrogen atom of the pyridine and the oxygen atom is an index of the molecular angle. The mono-protonate structure has longer distance than the neutral, with the di-protonated molecule representing the shortest. These modifications to structure may relate to the action mechanism. In terms of geometry, there is no significant difference between these two drugs. Under physiological conditions, both protonated species are present and both may contribute to incrementing the pH values in the intracellular compartments. These optimized structures should be taken into account, as the geometry of the molecule is important when there is interaction with the proteins of the virus. These optimized structure are folded, and this is in agreement with previous experimental results of the X-Ray structure [35]. Figures 1 and 2 also report the Mulliken atomic charges of some atoms. In all cases, the carbon atom that is bonded to the chlorine is positive (almost +1), whereas the chlorine atom is slightly negative. The oxygen atom of hydroxychloroquine has a negative atomic charge, as anticipated.

Fig. 1 Schematic representation and optimized structures of cloroquine, neutral and protonated. N-CH3 bond distances in Å (in bold) and Mulliken atomic charges (in italics) are included. 

Fig. 2 Schematic representation and optimized structures of hydroxychloroquine, neutral and protonated. N-O bond distances in Å (in bold) and Mulliken atomic charges (in italics) are reported. 

As the basic forms are used as drugs against diseases, it is important to analyze the di-protonated species in more detail. Fig. 3 reports the molecular frontier orbitals for the di-protonated molecules.

Fig. 3 Frontier molecular orbitals of di-protonated molecules. Eigenvalues in a.u. 

Cloroquine and hydroxychloroquine have practically the same molecular frontier orbitals. All molecular orbitals are located at the benzene and pyridine rings and the rest of the atoms do not participate. These orbitals are π-bonding orbitals. There is also a p non-bonding orbital of the chlorine atom and there are π-antibonding orbitals between chlorine and carbon atoms. In the case of hydroxychloroquine, the HOMO-3 orbital forms a π-bond above the region where the oxygen atom is placed. This orbital is different from the corresponding orbital of cloroquine. However, HOMO and LUMO are similar for both molecules. The electrons located in these orbitals can be said to be the active electrons, which may participate in charge transfer processes. Benzene and pyridine rings and chlorine atoms will accept electrons and the electrons will be ionized from the π-bonding orbitals of the rings.

Vertical ionization energy and vertical electron affinity were obtained for all the molecules and are reported in Table 1. As expected, protonated species have higher ionization energies and higher electron affinities than non-protonated molecules. The ionization energy reveals that it is energetically costly to donate electrons. The electron affinity indicates that it is energetically favorable to accept electrons, so that these molecules are oxidants. For comparison, we included the values for two oxidant free radical molecules in Table 1. Di-protonated molecules have higher electron affinity than NO, as they are better electron acceptors. This is important as these molecules may oxidize other molecules, such as proteins and DNA-nitrogen bases. Apparently, the oxidative capacity of these molecules is high, a factor which may affect the toxicity or effectiveness of these molecules.

Table 1 Vertical ionization energy and vertical electron affinity (IE and EA in eV) are reported. Electrodonating and electroaccepting powers (ω- and ω+ respectively) are also included.  

Molecule IE EA ω- ω+
Cloroquine 5.93 1.83 5.88 1.99
[Cloroquine-H]+1 5.99 1.89 6.01 2.07
[Cloroquine-2H]+2 6-77 2.58 7.81 3.15
Hydroxychloroquine 5.91 1.81 5.83 1.97
[Hydroxychloroquine-H]+1 5.98 1.87 5.96 2.04
[Hydroxychloroquine-2H]+2 6.75 2.55 7.74 3.09
NO 6.32 1.88 6.11 2.01
OCH3 8.30 4.75 15.46 8.94

Table 1 also reports the electrodonating and electroaccepting powers of the molecules under study. These quantities refer to the partial charge transfer processes. Both quantities are increasing according to the protonation. Non-protonated molecules are better electron donors and worse electron acceptors than protonated species.

Conclusion

Expressed simply, the hypothesis about the action mechanism of these molecules indicates that they increment the basicity of the environment, thus curtailing the reproduction of the virus. This is clear from the pKa values, as under physiological conditions these molecules are mono and di-protonated. However, these molecules may also interact by transferring electrons or oxidizing other molecules of the organism, and may also oxidize proteins in the virus. The idea that emerges from this study is that not only the acid-base equilibrium is important. Besides this, the reduction-oxidation equilibrium may explain the toxicity or the effectiveness of these drugs. What is certain is that these drugs may contribute to increased oxidative stress in the body. Being a molecule able to oxidize other molecules is similar to yin-yang; a dualism that describes contrary forces, as oxidation may contribute to attacking the virus, but also contributes to increasing oxidative stress. These results should help improve understanding of the action mechanisms.

Acknowledgment

This work was carried out using Miztli HP Cluster 3000 supercomputer (LANCAD-UNAM-DGTIC-141) provided by Dirección General de Cómputo y Tecnologías de Información y Comunicación (DGTIC), Universidad Nacional Autónoma de México (UNAM). Author would like to acknowledge Oralia L Jiménez, María Teresa Vázquez and Caín González for their technical support.

References

1. WHO. “International Health Regulations, WHA58.3” (World Health Organization, Geneva, ed. 2, 2005). [ Links ]

2. Weiss S. R.; Navas-Martin S.; Microbiol. Mol. Biol. Rev. 2005, 69, 635-664. [ Links ]

3. Masters P. S. Adv. Virus Res. 2006, 66, 193-292. [ Links ]

4. Bosch B.J.; van der Zee R.; de Haan C.A.; Rottier P.J. J. Virol. 2003, 77, 8801-8811. [ Links ]

5. Wu C.; Liu Y.; Yang Y.; Zhang P.; Zhong W.; Wang Y.; Wang O.; Xu Y.; Li M.; Li X.; Zheng M.; Chen L.; Li H. Acta Pharm. Sinica B. (2020) https://doi.org/10.1016/j.apsb.2020.02.008 [ Links ]

6. Tonnesmann E.; Kandolf R.; Lewalter T. Immunopharmacol Immunotoxicol 2013, 35, 434-442. [ Links ]

7. Ben-Zvi I.; Kivity S.; Langevitz P.; Shoenfeld Y. Clin. Rev. Allergy Immunol. 2012, 42, 145-153. [ Links ]

8. Vincent M.J.; Bergeron E.; Benjannet S.; et al. Virol. J. 2005, 2, 69, 1-10. [ Links ]

9. Warhurst D. C.; Steele J. C. P.; Adagu I.S.; Craig J.C.; Cullander C. J. Antimicrobial Chemotherapy 2003, 52, 188-193. [ Links ]

10. Taccone F.S.; Gorham J.; Vincent J.L. Lancet Respir. Med. 2020 Published Online April 15, 2020 https://doi.org/10.1016/ S2213-2600(20)30172-7. [ Links ]

11. Benjamin N.; Rome M.D.; Jerry Avorn M.D. New England J. of Med. This article was published on April 14, 2020, DOI: 10.1056/NEJMp2009457. [ Links ]

12. Gautret P.; Lagier J.C.; Parola P. et al. Int J Antimicrob Agents 2020; published online March 20. DOI:10.1016/j.ijantimicag.2020.105949. [ Links ]

13. Gautret P. ; Lagier J.C. ; Parola P. et al. 2020, https://www.mediterranee-infection.com/wp-content/uploads/2020/03/COVID-IHU-2-1.pdf (accessed April 2, 2020) [ Links ]

14. Alhazzani W.; Moller M.H.; Arabi Y.M. et al. Intensive Care Med 2020; published online March 28. DOI:10.1007/s00134-020-06022-5. [ Links ]

15. Al-Bari M.A. J. Antimicrob. Chemother 2015, 70, 1608-1621. [ Links ]

16. Frisch M. J.; Trucks G. W.; Schlegel H. B.; Scuseria G. E.; Robb M. A.; Cheeseman J. R.; Montgomery J. J. A.; Vreven T.; Kudin K. N.; Burant J. C.; Millam J. M.; Iyengar S. S.; Tomasi J.; Barone V.; Mennucci B.; Cossi M.; Scalmani G.; Rega N.; Petersson G. A.; Nakatsuji H.; Hada M.; Ehara M.; Toyota K.; Fukuda R.; Hasegawa J.; Ishida M.; Nakajima T.; Honda Y.; Kitao O.; Nakai H.; Klene M.; Li X. ; Knox J. E.; Hratchian H. P.; Cross J. B.; Bakken V.; Adamo C.; Jaramillo J.; Gomperts R.; Stratmann R.E.; Yazyev O.; Austin A. J.; Cammi R.; Pomell, C.; Ochterski J. W.; Ayala P. Y.; Morokuma K.; Voth G. A.; Salvador P.; Dannenberg J. J.; Zakrzewski V. G.; Dapprich S.; Daniels A. D.; Strain M. C.; Farkas O.; Malick D. K.; Rabuck A. D.; Raghavachari K.; Foresman J. B.; Ortiz J. V.; Cui Q.; Baboul A. G.; Clifford S.; Cioslowski J.; Stefanov B. B; Liu G.; Liashenko A.; Piskorz P.; Komaromi I.; Martin R. L.; Fox D. J.; Keith T.; Al-Laham M. A.; Peng C. Y.; Nanayakkara A.; Challacombe M.; Gill P. M. W.; Johnson B.; Chen W.; Wong M. W.; Gonzalez C.J.; Pople A. Gaussian 09, Revision A.08, Inc. Wallingford, CT, 2009. [ Links ]

17. Zhao Y.; Truhlar D. G. Theor. Chem. Acc. 2008, 120, 215-241. [ Links ]

18. Petersson G. A. ; Bennett A.; Tensfeldt T. G.; Al-Laham M. A. ; Shirley W. A.; Mantzaris J. A. J. Chem. Phys. 1988, 89, 2193-2218. [ Links ]

19. Petersson G. A. ; Al-Laham M. A. J. Chem. Phys. 1991,94, 6081-6090. [ Links ]

20. McLean A. D.; Chandler G. S. J. Chem. Phys. 1980, 72, 5639-5648. [ Links ]

21. Krishnan R.; Binkley J. S.; Seeger R; Pople J. J. Chem. Phys. 1980, 72, 650-654. [ Links ]

22. A) Mennucci B. ; Tomasi J. J. Chem. Phys. 1997, 106, 5151; (b) Tomasi J. ; Mennucci B. ; Cances E. J. Mol. Struct. (THEOCHEM) 1999, 464, 211. [ Links ]

23. Geerlings P.; Chattaraj P.K.; De Proft F.; Gázquez J.L.; Liu S.; Morell Ch.; Toro-Labbé A.; Vela A.; Ayers P. Theor. Chem. Acc. 2020, 139, 5-18. [ Links ]

24. Geerlings P. ; De Proft F. ; Langenaeke W. Chem. Rev. 2003, 103, 1793-1874. [ Links ]

25. Chermette H. J. Comput. Chem. 1999, 20, 129-154. [ Links ]

26. Liu SB. Acta Phys. Chem. Sin. 2009, 25, 590-600. [ Links ]

27. Wang B.; Rong Ch.; Chattaraj P.K. ; Liu SB. Theor. Chem. Acc. 2019, 138, 124. [ Links ]

28. Gázquez J. L.; Cedillo A.; Vela A. J. Phys. Chem. A. 2007, 111, 1966-1970. [ Links ]

29. Gázquez J. L. J. Mex. Chem. Soc. 2008, 52, 3-10. [ Links ]

30. Martínez A.; Rodríguez-Gironés M. A.; Barbosa A.; Costas M. J. Phys. Chem. A. 2008, 112, 9037-9042. [ Links ]

31. Martínez, A. J. Phys. Chem. B. 2009, 113, 4915-4921. [ Links ]

32. Cerón-Carrasco J. P.; Bastida A.; Requena A.; Zuñiga J.; Miguel B. J. Phys. Chem. B. 2010, 114, 4366-4372. [ Links ]

33. Alfaro R. A. D.; Gómez-Sandoval Z.; Mammino L. J. Mol. Model. 2014, 20, 2337. [ Links ]

34. García-Hernández E.; Garza J. J. Mex. Chem. Soc. 2017, 61(3), 222-228 [ Links ]

35. Semeniuk A.; Kalinowska-Tluscik J.; Nitek W.; Oleksyn B. J. J. Chem. Crystallogr. 2008, 38, 333-338. [ Links ]

Received: May 06, 2020; Accepted: June 22, 2020

*Corresponding author: Ana Martínez, email: martina@unam.mx

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License