Dear Label_Name,

Greetings! I hope you are well.

I am writing on behalf of an international publishing house, Science and Education Publishing, USA (http://www.sciepub.com).

Science and Education Publishing (SciEP) is a multidisciplinary open access publisher of Journals covering the fields of Science, Technology and Medicine. All SciEP publications are subject to high-quality peer review, editorial and production processes and are freely available online.

Welcome to the April, 2013 edition of SciEP newsletter, providing monthly summaries of articles on science & education related topics of interest to scientist and professionals. The mission of our journals is to provide an open access forum for publishing new findings on issues pertinent to scientific and academic field. We also invite you to submit original research papers, review articles, short reports, or case studies are considered for publication in our Multidisciplinary Journals.

American Journal of Public Health Research
DOI: 10.12691/ajphr-1-1-3

The Implementation of a Tele-Homecare System with Service Oriented Architecture and HL7 Message Transmission Standard
Ching-Sung Wang

American Journal of Medical and Biological Research
DOI: 10.12691/ajmbr-1-1-2

The Influence of Ethanol on Pyruvate Kinases Activity in Vivo, in Vitro, in Silico
Lelevich Sergey Vladimirovich, Khrustalev Vladislav Victorovich, Barkovsky Eugene Victorovich, Shedogubova Tatyana Aleksandrovna

American Journal of Medical Sciences and Medicine
DOI: 10.12691/ajmsm-1-1-1

Clinical and Histological Features in Wegener Granulomatosis
Li Cavoli Gioacchino, Passantino Rita, Bono Luisa, Tortorici Calogera, Ferrantelli Angelo, Giammarresi Carlo, Rotolo Ugo

American Journal of Infectious Diseases and Microbiology
DOI: 10.12691/ajidm-1-1-2

The Curse of Lymphatic Filariasis: Would the Continual Use of Diethylcarbamazine Eliminate this Scourge in Papua New Guinea?
Anne Lanham, Lillian Mwanri

We have released issues of some journals on February 2013.
American Journal of Public Health Research
American Journal of Medical Sciences and Medicine
American Journal of Medical and Biological Research
American Journal of Pharmacological Sciences
American Journal of Clinical Medicine Research
American Journal of Nursing Research
American Journal of Cancer Prevention
American Journal of Medicine Studies
American Journal of Infectious Diseases and Microbiology
American Journal of Cardiovascular Disease Research
American Journal of Biomedical Research
American Journal of Mechanical Engineering
American Journal of Microbiological Research
Journal of Preventive Medicine
American Journal of Food and Nutrition
American Journal of Hypertension Research
American Journal of Sports Science and Medicine
Journal of Cancer Research and Treatment
American Journal of Zoological Research
Journal of Biomedical Engineering and Technology
American Journal of Microbiological Research
Biomedical Science and Engineering

We are pleased to announce that submissions are now being accepted for peer review for all the journals across all our subject areas in SciEP.

Please submit your manuscript via our online submission system:
http://www.sciepub.com/author/signin

To unsubscribe from our mailing list, please click on the following link: unsubscribe

I am looking forward to hearing from you.

Kind regards,

David Rain
Email:join@sciepub.com
Science and Education Publishing (http://www.sciepub.com)

Plasmodium Resistance to Antimalarial Drugs: Functional Theory of Density (DFT) Study of the Stability and Reactivity of Heme-Artemisinin Adducts

Affoué Lucie Bédé, Amon Benjamine Assoma, El Hadji Sawaliho Bamba*, Maftei Dan, Ionel Humelnicu

Laboratoire de Chimie Organique Structuralle, Université Félix Houphouët-Boigny 22 BP582 Abidjan22 (Côte-d’Ivoire)

Abstract: This study is part of the search for solutions to the resistance of the parasite to antimalarials. The aim of this work is to analyze the relative stability and reactivity of adducts in order to understand the resistance of the parasite to artemisinin. The DFT/B3LYP/GenECPs method associated with bases 6-31G, 6-31G ** and LANL2DZ is used to perform these calculations. On one hand analysis of the thermodynamic quantities of formation and reaction revealed that β adduct is the most stable among the studied adducts. In other words, the latter is the majority product of the interaction between heme and artemisinin. On the other hand, the adduct γ is the least stable. The values of the pameters such as dipole moments, intramolecular hydrogen bond lengths and thermodynamic parameters (formation and reaction quantities) have shown that it can exist a correlation between the following three parameters those are the relative stability of adducts (α, β, δ), their dipole moment and the binding intramolecular hydrogen. The analysis of the Frontier Molecular Orbitals also showed that the β adduct is susceptible to electrophilic attack, unlike the γ adduct, which is rather favorable to nucleophilic attack and can be the most excited component. Moreover, the determination of the alkylation rate will be considered in order to specify the order of complexation of the meso positions α, β and δ of the heme.

Keywords: malaria, adducts, heme, artemisinin, DFT/B3LYP/genECPs

Citation: Affoué Lucie Bédé, Amon Benjamine Assoma, El Hadji Sawaliho Bamba, Maftei Dan, Ionel Humelnicu, Plasmodium Resistance to Antimalarial Drugs: Functional Theory of Density (DFT) Study of the Stability and Reactivity of Heme-Artemisinin Adducts. Journal of Materials Physics and Chemistry, 6 (1). 17-22. January 2018. doi:10.12691/jmpc-6-1-3

1. Introduction

Artemisinin is an antimalarial which is the basis of a lot of antimalarial drugs currently used in first-line treatment [1,2,3,4,5]. This antimalarial alkylates the heme on its meso positions α, β, γ and δ to form adducts [6,7]. This alkylation is weak in γ and strong in α, β and δ positions. In 2001, WHO (World Health Organization) considered artemisinin to be "the world's biggest hope in the struggle against malaria" [8]. Unfortunately, it was a short hope because since 2007, artemisinin-resistant parasites have appeared in Southeast Asia [9,10]. This molecule does not completely eliminate all malaria parasites, hence the need to combine it with other antimalarials [11,12]. But the first cases of clinical failures with artemisinin-based combination therapies (ACTs) have been identified in the same area (Southeast Asia) [13,14,15]. The problem of parasite resistance to antimalarials remains despite the introduction of new molecules. This study is a side of the research for solutions to the resistance of the parasite to antimalarials. It assumes that the instability of the adducts can explain the resistance of the parasite to artemisinin. The aim of this work is to analyze the relative stability of heme-artemisinin adducts (Figure 1) in order to understand the ineffectiveness of this antimalarial. The theoretical approach of adduct reactivity is based on the theory of Frontier Molecular Orbitals (FMO). The calculations are performed in gaseous phase and in solution using the DFT/B3LYP/GenECPs method associated with bases 6-31+G (d,p), 6-31G and LanL2DZ.

Figure 1. Molecular structures of α, β, γ and δ adducts.

2. Materials and Methods

2.1. Level of Calculation Theory

Theoretical studies are carried out with Gaussian 09 calculation software [16]. The different structures were optimized using density functional theory (DFT) with B3LYP / GenECPs as the level of calculation [17,18]. The atoms of carbon, nitrogen and oxygen (C, N and O) are optimized at B3LYP/6-31 + G (d, p) theory level. The atoms of hydrogen (H) and iron (Fe) are optimized respectively at the B3LYP/6-31G and at B3LYP/LanL2DZ level. LanL2DZ is a pseudo-potential base [19,20]. All these calculations are made by using the same method (DFT). An artemisinin derivative was added at one of the meso positions of the heme unit to form adducts in order to improve the representation of the coordination environment found in the parasite vacuole. As the alkylation reaction takes place in an acid environment, two solvents were used for the solution calculation, namely water and acetic acid. The Conductor-like Polarizable Continuum Model (CPCM) is the solvation model used for solution calculations [21]. The complete optimizations were carried out without any symmetry constraint. The harmonic vibration frequencies have been calculated to confirm that the optimized geometry corresponds exactly to a local minimum which does not possess a negative frequency.

2.2. Thermodynamic Parameters of Formation and Reaction

The knowledge of the variations of energy contributions to the internal energy at 0 K and at 298.15 K under normal atmospheric pressure (P = 1atm.) between products and reagents contributes to the energetic characterization of a chemical reaction [22,23]. For a given energy parameter X, its variation is determined according to the following relation (1):

[picture] | (1)

The studied energetic parameters are the variations of the electronic energy of reaction (internal energy of this reaction), the variation of the enthalpy and the variation of the free enthalpy . To get access to the internal energy, the contributions of different movements such as translation, rotation and vibration of the nuclei must therefore be taken into account in order to evaluate them. In the ideal gas approximation, the rotational and translational contributions are given by relation (2):

[picture] | (2)

The contribution ZPVE (Zero Point Vibrational Energy) which means the energy of the lowest vibrational level due to the 3N-6 normal modes of vibration (3N-5 for the linear molecules) of frequencies of N nuclei at 0 K, is defined by the relation below (3):

[picture] | (3)

To obtain the energy corresponding to 298 K, it is necessary to take into account the additional energy due to the population of the vibration levels during the temperature rising from 0 to 298 K, defined by the relation (4) below:

[picture] | (4)

At the end, the variation of internal energy at 298 K can be written like relation (5):

[picture] | (5)

The formation’s enthalpy of reaction at 298 K corresponds to the internal energy variation which has been corrected from the term Δ (PV), ie (–RT) because the variation in the number of moles is equal to (-1). This variation is given by the following expression (6):

[picture] | (6)

The estimation of the free reaction enthalpy is then obtained by equations (7) and (8).

[picture] | (7)

[picture] | (8)

2.3. Frontier Molecular Orbital Theory (FMO)

HOMO’s orbitals usually act as electron donors and LUMO’s ones as electron acceptors. They can give a reasonable qualitative prediction of the excited state properties of a molecule and the capability of a molecule to transport electron. They are also quantum chemical parameters those determine molecular reactivity [24]. The energy gap noted ΔE is the energy difference between the energy levels of HOMO (the highest occupied molecular orbital) and LUMO (lowest vacant molecular orbital). This difference (ΔE) serves as a measure of the molecule’s excitability. Thus, the lower the energy gap is, the more the molecule can interact with the environment. A large gap HOMO-LUMO implies a high stability for the molecule and low reactivity in terms of chemical reactions. Similarly, a small difference implies a high reactivity of the molecule [25,26,27].

3. Results and Discussion

3.1. Relative Stability of Adducts

Table 1 contains the enthalpies of formation (ΔfH298K) in atomic unit, the reaction enthalpies (ΔrH298K) and the free reaction enthalpies (ΔrG298K) in Kcal.mol⁻¹, the lengths of intramolecular hydrogen bonds (d) in angström and the dipolar moments (μ) in Debye of adducts α, β, γ and δ.As with tables and equations, figures should be set in one column if possible unless two-column display is essential. The resolution of graphics and image should be adequate to reveal the important detail in the figure.

Table 1. Formation Enthalpies, Reaction Enthalpies, Free Reaction Enthalpies, Dipole Moments and Intramolecular Hydrogen Bond Lengths Calculated at the B3LYP/GenECPs Theory Level

The comparison of the enthalpies of formation mentioned in Table 1 indicates that the alkylation of heme on the meso position β gives the lowest energy and the γ position leads to the highest energetic value. α and δ adducts’ energies are very closed to the same energetic value so they have a similar relative stability. The enthalpies of formation of these adducts increase in acetic acid and decrease in water which means that these adducts are more stable in water than in the acidic area. β adduct remains the most stable molecule and the least stable is γ adduct because they possess respectively the lowest and the highest values of energy. These results permit us to establish the following descending order of stability in acidic environment: β > δ > α > γ. In water, α adduct is more stable than δ adduct and it allow us to get the following order of decreasing stability: β > α > δ > γ. The instability of the adduct γ can be explained by a steric gene (trouble) due to the presence of carboxyl groups. Robert et al. [6,7] showed that the alkylation rate of heme is high in α, β and δ positions, and low in the γ’s one. However, their work does not give any comment on the relative stability of these adducts. The results of our work allow us to establish the order of relative stability of these molecules. The alkylation of heme in the β position leads to the most stable molecule. We can argue that the major product of the interaction between heme and artemisinin is β adduct.

The dipole moment (μ) indicates the stability of a molecule in water, in particular in aqueous solution. Thus, a strong dipole moment will result low solubility in organic solvents and high solubility in water [28,29]. Like enthalpies of formation, the dipole moment increases in the acid environment. α adduct has the highest dipole moment. The relative stability of α, β and δ adducts increases with the dipole moment that is to say more the value of the dipole moment of the adduct is high, less the adduct is stable. Thus, α adduct with the highest dipole moment value (6.1887 and 6.9855 Debye respectively in water and in acid), is the least stable among the studied adducts. The most stable is β adduct because it has the smallest value of dipole moment (3.9948 and 4.1574 Debye respectively in water and in acid). There is a correlation between the relative stability of these adducts and their dipole moment.

The free enthalpy values and calculated enthalpies of reaction (ΔrH298K and ΔrG298K) are all positive (Table 1). The alkylation at 298.15 K is an endothermic and non-spontaneous reaction at this temperature. In other words, the alkylation reaction takes place at high temperatures.

Figure 2. (A): Hydrogen bond establishment between hydroxyl groups of α and δ adducts; (B): hydrogen bond between carbonyl and hydroxyl groups of the β adduct

Figure 3. Structures of the C1, C2, C3 and C4 heme’s conformers

Concerning H bond lengths (distance d), the practice consist to consider a contact as a H bond if the distance d is smaller than the sum of Van der Waals radii, by taking 1.52 Å for the oxygen atom and 1.0 Å for hydrogen one; that means d ≤ 2.52 Å [30,31]. It is also known that the (link) H bond is strong when its length is short [32,33,34]. Moreover, studies have shown that the hydrogen bond is stronger when the energy of the reaction is weak [35,36,37]. The analysis of the intramolecular hydrogen bond across the H bond lengths of the studied adducts shows that the length of the H bond is shorter in the β adduct. It varies from 1.688 to 1.704 Å (Table 1). In α and δ compounds, the values of the hydrogen bond are closer and are comprised between 1.726 Å and 1.733 Å. Concerning the standard free reaction enthalpy, low values are also recorded around the β adduct in water and in acetic acid, which are respectively 21.062 and 22.406. K.cal.mol⁻¹. All these observations show that there is a strong existence of intramolecular hydrogen bond around the β adduct. This H bond is established between a carboxyl group and the hydrogen atom of the hydroxyl group. Contrary to the β adduct, the hydrogen bond in the other compounds occurs between the two hydroxyl groups (Figure 2). The groups that interact in hydrogen bond formation differ from one adduct to another. This observation assumes the existence of heme conformers based on intramolecular hydrogen bonding between hydroxyl and carbonyl groups. The different structures of the heme conformers are shown in Figure 3. The β-adduct which is the most stable molecule has the smallest distances of hydrogen bond lengths. α and δ compounds with neighboring energy values also have very close hydrogen bonding values. The less stable γ compound does not have hydrogen bond. This observation supposes that there is a correlation between the relative stability of α, β and δ adducts and the hydrogen bond.

3.2. Relative Stability of Heme Conformers

Figure 4 shows the different structures of the heme. Calculations of heme’s conformers are carried out in gaseous phase. The energies of enthalpies of formation, boundary orbitals, dipole moments and hydrogen bond lengths are summarized in Table 2.

The results show that the C3 conformer has the lowest value of formation enthalpy (Table 2). It is therefore the most stable. The least stable is the C4 conformer. This assumes that in the gas phase, the heme is in the form of the C3 conformer. The order of relative decreasing stability of these molecules is as follows: C3 > C1 > C2 > C4. The conformers C2 and C4 have the highest dipole moment values. The energy gaps indicate that the conformer C3 is the most stable. The dipole moment of C3 conformer is relatively weak (2.160 Debye) compared with those of the C1, C2 and C4 molecules whose dipole moments are 6.876, 10.081 and 10.408 Debye respectively (Table 2). The relative stability of the conformers increases with their dipole moment. The mentioned values in Table 2 also show that the shorter the hydrogen bond length, the lower the energy of the molecule. These results highlight the probable existence of a correlation between the relative stability of these molecules, their dipole moment and their intramolecular hydrogen bonding.

Figure 4. Probable transformation reaction of heme-artemisinin adducts with A the structure of the alkylated heme by an artemisinin group in position β (β adduct) and B a probable structure of alkylated heme having in one of its meso positions a less bulky group

Table 2. Formation Enthalpies, Boundary Orbitals, Dipole Moments and Hydrogen Bond Lengths Calculated at the B3LYP/genECPs Level

3.3. Reactivity by the Molecular Orbitals Theory Frontiers

The energies of the HOMO and LUMO boundary molecular orbitals were calculated at the DFT/B3LYP/genECPs level. The gap ΔE = ELUMO - EHOMO is estimated from these energies. The results are shown in Table 3.

Table 3. Energies of HOMO, LUMO Boundary Orbitals and HOMO-LUMO Energy Gaps of Adducts Calculated at B3LYP/genECPs Level

HOMO levels of adducts increase and LUMO levels decrease with acidic solvent (Table 3). The molecular orbital boundary HOMO reflects the electro-donor (nucleophilic) character of the molecule. The higher the energy of this OM, the easier the molecule will yield electrons. As for LUMO, it reflects the electro-acceptor (electrophile) character of the molecule. The lower the energy of this molecular orbital, the easier the molecule will accept electrons. The β adduct has the highest energy of HOMO level so this molecule will be susceptible to electrophilic attack. The γ adduct’s LUMO level energy is the lowest one. This molecule can be favorable to a nucleophilic attack. The smallest value of ΔE is obtained with the adduct γ. The β adduct leads to the largest HOMO-LUMO gap. The γ adduct is susceptible to be more reactive than the other adducts and β adduct is the most stable among the studied molecules. The gap values (ΔE) are close to one another (Table 3) and close to that of the adduct γ. This reflects the neigh bourhood of reactivity of these molecules. The continuation of the alkylation reaction which leads to an alkylated heme molecule with a less bulky group (Figure 4) can be said. Thus this last structure could more easily undergo a dimerization. This could explain the resistance of the parasite to artemisinin. These may be able to reduce the steric hindrance (trouble) effect of the artemisinin moiety and thus facilitate the formation of heme dimers.

4. Conclusion

Analysis of formation enthalpies showed that the β-position alkylation of heme gives the most stable adduct. This adduct seems to be the major product of the heme-artemisinin interaction. That in position γ leads to the least stable adduct and the most reactive. This could justify the low level of alkylation of heme observed experimentally in this position. Our results showed a probable correlation between the relative stability of α, β, δ adducts, their dipole moment and intramolecular hydrogen bonding. The energy gap values indicate that the reactivity of these molecules are all closed. In this sense, the possibility of a reaction that would transform these adducts into other molecules is possible. These may be able to reduce the steric hindrance effect of the artemisinin moiety and thus facilitate the formation of heme dimers. Our results also showed that the groups involved in hydrogen bond formation differ from one adduct to another. The results of this work have opened other avenues of scientific investigation such as the probable existence of a transformation of adducts into other molecules. The determination of the alkylation rate in the framework of the study of the mechanism of the reaction of the heme-artemisinin interaction could bring more precision to this first conclusion. All these lines of investigation would contribute to a probable explanation of the resistance of Plasmodium to artemisinin.

Acknowledgments

The calculations were carried out by Dr. Affoué Lucie Bédé during a post-doctoral research internship in the Department of Physical and Theoretical Chemistry of Alexandru Ioan Cuza University of Iasi-Romania, led by Professor Ionel HUMELNICU, as beneficiary of a postdoctoral fellowship under the "Eugen Ionescu" program 2015-2016 of AUF BECO. We thank them for allowing the calculations in this publication.

References

[1]   | Timms, A.R. and Bueding, E., “Studies of a proteolytic enzyme from Schistosoma mansoni”, Brit. J.Pharmacol., 14 (1), 68-73, March 1959.
[2]   | Vroman, J.A. and Alirm, G.M., “Current progress in the chemistry, medicinal chemistry and drug design of artemisinin based antimalarials”, Curr. Pharm. Design, 5, 101-138, February 1999.
[3]   | Touze, J.E., Fourcade, L., Pradines, B., Hovette, P., Paule, P. and Heno, Ph., “Les modes d’action des antipaludiques. Intérêt de l’association atovaquone-proguanil”, Symposium, Medecine Tropicale 2002, 62, 219-224.
[4]   | Pradines, B., Vial, H. and Olliaro, P., “Prophylaxie et traitement du paludisme : problèmes, récents développements et perspectives”. Synthèse, Medecine Tropicale 2003, 63, 79-98.
[5]   | Ebiamadon, A., Umoren E.U., Fraideh, B., Pedro M.M., Jorge, F.S.F., Devanand, L., Xianli, W. and Ronald L.P., “Nutritional characterisation and antioxidant capacity of different tissues of Artemisia annua L.”, Food Chemistry, 115(4), 1240-1246, August 2009.
[6]   | Robert, A., Cazelles, J. and Meunier, B., “Characterization of the alkylation product of heme by the antimalarial drug artemisinin”. Angew. Chem. Int. Ed., 40, 1954-1957, May 2001.
[7]   | Robert, A., Coppel, Y. and Meunier, B., “Alkylation of heme by the antimalarial drug artemisinin”, Chem. Commun, 0, 414-416, Janvier 2002.
[8]   | Pradinesa, B., Dormoi, J., Briolant S., Bogreau H. and Rogier, C., “La resistance aux antipaludiques”, Revue Francophone Des Laboratoires, 2010 (422), 51-62, May 2010.
[9]   | Dondorp, A.M., Fairhurst, R.M., Slutsker, L., John, R.M.A., Breman, J.G., Guerin, P.J., Wellems, T.E., Ringwald, P., Newman, R.D. and Plowe, C.V., “The threat of artémisinine-resistant malaria”, N. Engl. J. Med., 365, 1073-1075, September 2011.
[10]   | Didier, M., Frédéric, A. and Odile, M.P., “Étude de la résistance de Plasmodium falciparum aux antipaludiques au sein du réseau international des Instituts Pasteur (RIIP-Palu)”, Med. Sci. (Paris), 29(6-7), 647–655, July 2013.
[11]   | Nosten, F. and White, N.J, “Artemisinin-based combination treatment of falciparum malaria”, Am J Trop Med Hyg , 77(6), 181-192, July 2007.
[12]   | Eastman, R.T. and Fidock, D.A., “Artemisinin-based combination therapies: a vital tool in efforts to eliminate malaria”, Nat. Rev. Microbiol., 7, 864-874, November 2009.
[13]   | Noedl, H., Se, Y., Schaecher, K. and Smith, B.L., “Evidence of artemisinin-resistant malaria in western Cambodia”, N. Engl. J. Med., 359, 2619-2620, December 2008.
[14]   | Wongsrichanalai, C. and Meshnick, S.R., “Declining artesunate-mefloquine efficacy against falciparum malaria on the Cambodia-Thailand border”, Emerg. Infect. Dis., 14 (5), 716–719, May 2008.

[15]  
Dondorp, A.M., Nosten, F., Yi, P., Das, D., Phyo, A.P., Tarning, J., Lwin, K.M., Ariey, F., Hanpithakpong, W., Lee, S.J., Ringwald, P., Silamut, K., Imwong, M., Chotivanich, K., Lim, P., Herdman, T., An, S.S., Yeung, S., Singhasivanon, P., Nicholas, P.J., Lindegardh, N., Socheat, D. and Nicholas, J., “Artemisinin Resistance In Plasmodium Falciparum Malaria”, N. Engl. J. Med., 361, 455-467, July 2009.

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

[17]   | Becke, A.D.,. “A new mixing of Hartree–Fock and local density-functional theories”, J. Chem. Phys., 98(2), 1372-1377, October 1993.
[18]   | Lee, C., Yang, W. and Parr, R.G., “Development of the Colle–Salvetti correlationenergy formula into a functional of the electron density”, Phys. Rev. B, 37(2), 785-789, January 1988.
[19]   | Hay, P.J. and Willard, R.W., “Ab initio effective core potentials for molecular calculations. Potentials for K to Au including the outermost core orbitals”, J. Chem. Phys., 82(1), pp 299-310, September 1985.
[20]   | Lovallo, C.C. and Klobukowski, M., “Development of new pseudopotential methods: Improved model core potentials for the first-row transition metals”, J. Comput. Chem., 24, 1009-1015, May 2003.
[21]   | Cossi, M., Rega, N., Scalmani, G. and Barone, V., “Energies, structures, and electronic properties of molecules in solution with the C-PCM solvation model”, J. Comput. Chem., 24, 669-681, March 2003.
[22]   | J. W. Ochterski, “Thermochemistry in Gaussian”, 2000, Gaussian, Inc., 1-19, June 2, 2000. help@gaussian.com
[23]   | Becke, A. D., “Density‐Functional Thermochemistry. III. The Role of Exact Exchange”, J. Chem. Phys., 98, 5648-5652, 1993.
[24]   | Belletête, M., Morin, J.F., Leclerc, M. and Durocher, G., “A Theoretical, Spectroscopic, and Photophysical Study of 2, 7-Carbazolenevinylene-Based Conjugated Derivatives”, J. Phys. Chem. A, 109(31), 6953-6959, July 2005.
[25]   | Lewis, D.F.V., Ioannides, C. and Parke, D.V., “Interaction of a series of nitriles with the alcohol-inducible isoform of P450: Computer analysis of structure-activity relationships”, Journal Xenobiotica, 24(5), 401-408, September 2008.
[26]   | Zhou, Z. and Parr, R.G., “Activation hardness: new index for describing the orientation of electrophilic aromatic substitution”, J. Am. Chem. Soc., 112(15), 5720-5724, July 1990.
[27]   | Kikuchi, O., “Systematic QSAR Procedures with Quantum Chemical Descriptors”, Molecular Informatics, 6(4), 179–184, 1987.
[28]   | Rivelino, R., Canuto, S. and Canuto, K., “Conformational Stability of Furfural in Aqueous Solution: The Role of Hydrogen Bonding”, Braz. J. Phys., 34(1), 84-89, March 2004.
[29]   | Wong, M.W., Frisch, M.J. and Wiberg, K. B., “Solvent effects. 1. The mediation of electrostatic effects by solvents”, J. Am. Chem. Soc., 113 (13), 4776-4782, June 1991.
[30]   | Bondi, A., “van der Waals Volumes and Radii”, J. Phys. Chem., 68 (3), 441-451, March 1964.
[31]   | Rowland, R.S., “Taylor, R.,.Intermolecular nonbonded contact distances inorganic crystal structures: comparison with distances expected from van der Waals Radii”, J. Phys. Chem., 100 (18), 7384-7391, May 1996.
[32]   | Ziao, N., Laurence, C. and Le Questel, J.Y., “Amino nitrogen and carbonyl oxygen in competitive situations: which is the best hydrogen-bond acceptor site?”, Cryst. Eng. Comm., 4, 326-335, Apr 2002.
[33]   | Desiraju, G.R. and Steiner, T., The Weak Hydrogen Bond in Structural Chemistry and Biology, Oxford, University Press, 1999.
[34]   | Kassi, K.F., Koné, M.G-R, Affi, S.T. and Ziao, N., “ONIOM Method Characterization of Hydrogen Bonding Sites of Mycolactone A/B, a Buruli Ulcer Toxin”, Computational Chemistry, 5, 103-112, July2017.
[35]   | Kone, M.G.R., Affi, S.T., Ziao, N., Bamba, K. and Assanvo, E.F., “Hydrogen bonding sites in benzimidazolyl-chalcones molecules: An ab initio and DFT investigation”, Journal of Chemical and Pharmaceutical Research, 7(12), 805-812, 2015.
[36]   | Traore, Y., Bamba, K., Ziao, N., Affi, S.T. and Koné, M.G-R., “Quantum Chemical Characterization of Hydrogen Bonding Sites in Three 4-(4-Halo-Phenyl)-6-(Furan-2-yl) Pyrimidin-2-Amine Derivatives”, Computational Chemistry, 5, 91-102, July 2017.
[37]   | Perrin, C.L. and Nielson, J.B., “Strong hydrogen bonds in chemistry and biology”, Annu. Rev. Phys. Chem. 48, 511-544, October 1997.