Molecular docking studies, structural and spectroscopic properties of monomeric and dimeric species of benzofuran-carboxylic acids derivatives: DFT calculations and biological activities

Computational-based investigation of antioxidative potential polyphenolic compounds of Salvia officinalis L.: combined DFT and molecular docking approaches

CONTEXT

The antioxidant properties of the three polyphenolic compounds (carnosol, cirsiliol, and luteolin) of Salvia officinalis L. were investigated employing the density functional theory (DFT) calculations at the B3LYP of basis set at 6-311 +  + G (d, p) in order to evaluate their antioxidant activity. The enthalpies of reactions associated with the SET-PT, SPLET, and HAT mechanisms were analyzed in gas and in different solvents using the CPCM (conductor-like polarizable continuum) model. For all possible hydrogen donor sites, the corresponding parameters (BDE, AIP, PDE, PA, ETE, HOMOs, and LUMOs) and reactivity indices (IPE, EA, Χ, η, S, and ω) were also evaluated. The calculated results showed that derivatives 12-OH, 11-OH, 4'-OH, and 3'-OH had the lowest antioxidant activity. The results showed as well that carnosol, cirsiliol, and luteolin have higher reactivity compared to ascorbic acid and could be considered better antioxidants. According to research, the catechol group is crucial in influencing the studied compounds antioxidant activity. The theoretically predicted order of antioxidant efficiencies in this work agrees well with the QSAR (quantitative structure-activity relationship) data. The findings show that in the vacuum as well as benzene media. HAT would be the most effective mechanism; in contrast, the thermodynamic equilibrium approach in polar media is the SPLET mechanism. Likewise, the outcomes of the docking modeling confirm that the selected molecules have high inhibitory activity to glutathione-S-transferases (GSTs) receptors. Moreover, they have very important pharmacokinetic, chemical, and biological profiles. Finally, all the results show that the three natural molecules have good pharmacokinetic profiles, particularly the bioavailability and permeability toward biological membranes.

METHODS

The software packages used in this investigation are Gaussian 16, Discovery studio Visualizer, and AutoDock vina. The three compounds (carnosol, cirsiliol, and luteolin) of Salvia officinalis L. were optimized with DFT/B3LYP of basis set at 6-311 +  + G (d, p). The optimized structures were established via vibrational analysis (i.e., no imaginary frequencies in the frequency set). All enthalpies were zero-point (ZPE) corrected. Vibrational frequency calculations were performed at 298.15 K and 1 atmosphere pressure to determine the thermodynamic characteristics of the investigated reactions. The descriptors were associated with the antioxidant mechanisms for investigated molecules in vacuum and in various solvents. The molecular docking was used by AutoDock vina to estimate and evaluate the title compounds compatibility as potential antioxidant drugs utilizing appropriate receptor proteins. The solvation effect in the medium of benzene (ɛ = 2.27) and water (ɛ = 78.39) was taken into account. Furthermore, a methanol solvent (ɛ = 32.61) was also taken into consideration to compare with the empirical data.

Identification of potential anti-hyperglycemic compounds in Cordyceps militaris ethyl acetate extract: in vitro and in silico studies

Cordyceps militaris has been long known for valuable health benefits by folk experience and was recently reported with diabetes-tackling evidences, thus deserving extending efforts on screening for component-activity relationship. In this study, experiments were carried out to find the evidence, justification, and input for computations on the potential against diabetes-related protein structures: PDB-4W93, PDB-3W37, and PDB-4A3A. Liquid chromatography identified 14 bioactive compounds in the ethyl acetate extract (1-14) and quantified the contents of cordycepin (0.11%) and adenosine (0.01%). Bioassays revealed the overall potential of the extract against α-amylase (IC50 = 6.443 ± 0.364 mg.mL-1) and α-glucosidase (IC50 = 2.580 ± 0.194 mg.mL-1). A combination of different computational platforms was used to select the most promising candidates for applications as anti-diabetic bio-inhibitors, i.e. 1 (ground state: -888.49715 a.u.; dipole moment 3.779 Debye; DS ¯ -12.3 kcal.mol-1; polarizability 34.7 Å3; logP - 1.30), 10 (ground state: -688.52406 a.u.; dipole moment 5.487 Debye; DS ¯ -12.6 kcal.mol-1; polarizability 24.9 Å3; logP - 3.39), and 12 (ground state: -1460.07276 a.u.; dipole moment 3.976 Debye; DS ¯ -12.5 kcal.mol-1; polarizability 52.4 Å3; logP - 4.39). The results encourage further experimental tests on cordycepin (1), mannitol (10), and adenosylribose (12) to validate their in-practice diabetes-related activities, thus conducive to hypoglycemic applications.Communicated by Ramaswamy H. Sarma.

Ternary Mixture of Pentanamide in Solvent Analogy with Halogenated Phenol: Experimental, Theoretical, and In Silico Biological Studies

This research describes the preparation of mixtures of new halogen-substituted phenol derivatives and their effects due to linkages with a fatty amide (pentanamide). The molecules were optimized using DFT, and the vibrational and electronic analysis was done subsequently. The energies of frontier molecular orbitals (FMOs) were used to estimate the global chemical reactivity parameters as we suggest that hydrogen-bonded networks may have contributed to the stability and reactivity of the compound. In addition to the experimental investigation, dielectric parameters were calculated. Fukui functions were analyzed to study the chemical reactivity. To get insight into interactions of σ → π* orbitals, natural bond orbital calculations were done. Additionally, surface analysis of the MEP and Hirshfeld charges were performed at the equivalent DFT levels. The research also indicated that both (interaction region indicator) IRI and (electron delocalize range) EDR would proficiently identify chemical-bonding and weak interaction regions, providing a significant advantage in exploring diverse chemical systems and reactions. This indicated that compounds could diffuse through noncovalent interactions, including intramolecular hydrogen bonding. Dielectric relaxation studies taken at five distinct molar ratios identified significant dielectric properties such as ε', ε″, ε0, and ε∞. The PA with FP, CP, BP, and IP molecules has potential antiviral and antioxidant benefits for carbonic anhydrase, with favorable drug-like features and diverse biological benefits. Pharmacological effects were forecasted using the PASS server, and these molecules exhibited favorable pharmacokinetic properties.

Identification of potential inhibitors against Alzheimer-related proteins in Cordyceps militaris ethanol extract: experimental evidence and computational analyses

Laboratory experiments were carried out to identify the chemical composition of Cordyceps militaris and reveal the first evidence of their Alzheimer-related potential. Liquid chromatography-mass spectrometry analysis identified 21 bioactive compounds in the ethanol extract (1-21). High-performance liquid chromatography quantified the content of cordycepin (0.32%). Bioassays revealed the overall anti-Alzheimer potential of the extract against acetylcholinesterase (IC50 = 115.9 ± 11.16 µg mL-1). Multi-platform computations were utilized to predict the biological inhibitory effects of its phytochemical components against Alzheimer-related protein structures: acetylcholinesterase (PDB-4EY7) and β-amyloid protein (PDB-2LMN). In particular, 7 is considered as a most effective inhibitor predicted by its chemical stability in dipole-based environments (ground state - 467.26302 a.u.; dipole moment 11.598 Debye), inhibitory effectiveness (DS ¯ - 13.6 kcal mol-1), polarized compatibility (polarizability 25.8 Å3; logP - 1.01), and brain penetrability (logBB - 0.244; logPS - 3.047). Besides, 3 is promising as a brain-penetrating agent (logBB - 0.257; logPS - 2.400). The results preliminarily suggest further experimental attempts to verify the pro-cognitive effects of l(-)-carnitine (7).

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03714-9.

In-silico assay of a dosing vehicle based on chitosan-TiO2 and modified benzofuran-isatin molecules against Pseudomonas aeruginosa

A high priority of the World Health Organization (WHO) is the study of drugs against Pseudomonas aeruginosa, which has developed antibiotic resistance. In this order, recent research is analyzing biomaterials and metal oxide nanoparticles, such as chitosan (QT) and TiO2 (NT), which can transport molecules with biological activity against bacteria, to propose them as drug carrier candidates. In the present work, 10 modified benzofuran-isatin molecules were studied through computational simulation using density functional theory (DFT) and molecular docking assays against Hfq and LpxC (proteins of P. aeruginosa). The results show that the ligand efficiency of commercial drugs C-CP and C-AZI against Hfq is low compared with the best-designed molecule MOL-A. However, we highlight that the influence of NT promotes a better interaction of some molecules, where MOL-E generates a better interaction by 0.219 kcal/mol when NT is introduced in Hfq, forming the system Hfq-NT (Target-NT). Similar behavior is observed in the LpxC target, in which MOL-J is better at 0.072 kcal/mol. Finally, two pharmacophoric models for Hfq and LpxC implicate hydrophobic and aromatic-hydrophobic fragments.

Vibrational Spectroscopies, Global Reactivity, Molecular Docking, Thermodynamic Properties and Linear and Nonlinear Optical Parameters of Monohydrate Arsonate Salt of 4-Aminopyridine

In this work, a structural and electronic properties of a novel organic arsenate template by 4-aminopyridine, with the general formula (C₅H₇N₂)(C₅H₈N₂)[AsO₄]·H₂O ((4-APH)(4-APH₂)[AsO₄]·H₂O) have been presented. The density functional theory (DFT) along with B3LYP hybrid functional is employed. The optimized structure was found to be in well consistent with the X-ray diffraction geometry. The examination of the vibrational spectrum was correlated by DFT calculation using the unit cell parameters obtained from the experiment data. Besides, the thermodynamic functions (heat capacity, entropy, enthalpy) from spectroscopic data by statistical methods were obtained for the range of temperature 100–1000 K. In addition, the molecular orbital calculations such as Natural Bond Orbitals (NBOs), AIM approach, HOMO–LUMO energy gap, NLO characteristic and Hirshfeld surface analysis were also performed with the same level of DFT. Electronic stability of the compound arising from hyper conjugative interactions and charge delocalization were also investigated based on the natural bond orbital (NBO) analysis. Molecular docking studies were also conducted as part of this study. The theoretical results showed an excellent agreement with the experimental values.

Structure-activity relationships and the underlying mechanism of α-amylase inhibition by hyperoside and quercetin: Multi-spectroscopy and molecular docking analyses

Inhibiting the activity of α-amylase has been considered an effective strategy to manage hyperglycemia. Hyperoside and quercetin are the main natural flavonoids in various plants, and the inhibition mechanism on α-amylase remains unclear. In this study, the structure-activity relationships between hyperoside/quercetin and α-amylase were evaluated by enzyme kinetic analysis, multi-spectroscopic techniques, and molecular docking analysis. Results showed that hyperoside and quercetin exhibited significant α-amylase inhibitory activities with IC₅₀ values of 0.491 and 0.325 mg/mL, respectively. The α-amylase activity decreased in the presence of hyperoside and quercetin in a competitive and noncompetitive manner, respectively. UV-vis spectra suggested that the aromatic amino acid residues (Trp and Tyr) microenvironment of α-amylase changed in the presence of these two flavonoids. FTIR and CD spectra showed the vibrations of the amide bands and the secondary structure content changes. The fluorescence quenching mechanism of α-amylase by hyperoside and quercetin belonged to the static quenching type. Finally, molecular docking intuitively showed that hyperoside/quercetin formed hydrogen bonds with the key active site residues (Asp197, Glu233, and Asp300) in α-amylase. MD simulation indicated hyperoside/quercetin-α-amylase docked complexes had good stability. Taken together, this research provides new sights to developing potent drugs or functional foods with hyperoside and quercetin, offering new avenues for hyperglycemia treatment.

The antioxidative potential of procyanidin B1: DFT (density functional theory) and docking approaches

Procyanidin B1 is one of the natural dimeric flavonoids. It has established a great role in antioxidative activity. In the current study, we wish to provide crucial information on its antioxidative action by the DFT computational and docking approaches. From point of thermodynamic view, at the M062X/6-311G(d,p) level, the HAT (hydrogen atom transfer) and SPL-ET (sequential proton loss-electron transfer) are principal antioxidative routes of this compound in gas and methanol, respectively. OH groups of two phenyl rings of this molecule are likely to be the best antiradical sites. In the kinetics of the interactions with HOO^• radicals, OH groups of phenyl rings have also generated the best ΔG^# (Gibbs free energy of activation) and rate constant K. The antioxidative action of procyanidin B1 is further confirmed by its chelation to metal ions, in which complex formation with Cu²⁺ having lower binding energy is more stable than complex formation with Zn²⁺. Docking study revealed that the antioxidative activity of procyanidin B1 involved human tyrosinase enzyme inhibition through interaction with essential residues, focusing on the OH groups of two phenyl rings.

Facile synthesis, crystal structure, biological evaluation, and molecular modeling studies of N-((4-acetyl phenyl) carbamothioyl) pivalamide as the multitarget-directed ligand

The crystal structure of N-((4-acetylphenyl)carbamothioyl)pivalamide (3) was synthesized by inert refluxing pivaloyl isothiocyanate (2) and 4-aminoacetophenone in dry acetone. The spectroscopic characterization (¹H-NMR, ¹³CNMR, FT-IR) and single crystal assays determined the structure of synthesized compound (3). Systematic experimental and theoretical studies were conducted to determine the molecular characteristics of the synthesized crystal. The biological examination of (3) was conducted against a variety of enzymes i.e., acetyl cholinesterase (AChE), butyl cholinesterase (BChE), alpha amylase, and urease enzyme were evaluated. The crystal exhibited approximately 85% enzyme inhibition activity against BChE and AChE, but only 73.8 % and 57.9% inhibition activity against urease and alpha amylase was observed respectively. The theoretical calculations were conducted using density functional theory studies (DFTs) with the 6–31G (d, p) basis set and B3LYP functional correlation. The Frontier molecular orbital analysis revealed that the HOMO/LUMO energy gap was smaller, which corresponds to the molecule’s reactivity. In terms of reactivity, the chemical softness value was found to be in good agreement with experimental values. In Crystal structure analysis, the intramolecular N—H•••O hydrogen bond generates a S 6) ring motif and N—H•••O interactions exist in crystal structure between the centroids of neighboring parallel aromatic (C4-C9) rings with a centroid to centroid distance of 3.9766 (7)Å. These intermolecular interactions were useful in structural stabilization. The Hirshfeld surfaces and their related two-dimensional fingerprint plots were used for thorough investigation of intermolecular interactions. According to Hirshfeld surface analysis of the crystal structure the most substantial contributions to the crystal packing are from H ••• O and H ••• N/N ••• H interactions. Molecular docking studies were conducted to evaluate the binding orientation of synthesized crystal with multiple targets. The compound exhibited stronger interactions with AChE and BChE with binding energies of -7.5 and -7.6 kcal/mol, respectively. On the basis of in-vitro and in-silico findings, it is deduced that N-((4-acetylphenyl)carbamothioyl)pivalamide 3) possesses reactive and potent multiple target inhibitory properties.

Antioxidant activity and α-glucosidase inhibitability of Distichochlamys citrea M.F. Newman rhizome fractionated extracts: in vitro and in silico screenings

Distichochlamys citrea M.F. Newman (commonly known as “Black Ginger”) is an endemic plant to Vietnam and has been extensively exploited by folk medication for treatments of infection-related diseases and diabetes. In this work, its rhizomes were subjected to fractionated extraction, phytochemical examination, evaluation of antioxidant effect by DDPH free radical neutralization, and inhibitory activity toward α-glucosidase. The compositional components were subjected to in silico screening, including density functional theory calculation, molecular docking simulation, physicochemical analysis, and pharmacokinetic regression. In the trials, EtOAc fraction is found as the bioactive part of most effectiveness, regarding both antioxidant effect (IC₅₀ = 90.27 µg mL⁻¹) and α-glucosidase inhibitory activity (IC₅₀ = 115.75 μg mL⁻¹). Chemical determination reveals there are 13 components of its composition. DFT-based calculations find no abnormal constraints in their structures. Docking-based simulation provides order of inhibitory effectiveness: 3-P53341 > 12-P53341 > 7-P53341 > 4-P53341 > 11-P53341 > 10-P53341. QSARIS-based investigations implicate their biocompatibility. ADMET-based regressions indicate that all candidates are generally safe for medicinal applications. The findings would contribute to the basis for further studies on the chemical compositions of Distichochlamys citrea and their biological activities.

Modelling the DFT structural and reactivity study of feverfew and evaluation of its potential antiviral activity against COVID-19 using molecular docking and MD simulations

Abstract

The unavailability of a proper drug against SARS-CoV-2 infections and the emergence of various variants created a global crisis. In the present work, we have studied the antiviral behavior of feverfew plant in treating COVID-19. We have reported a systematic in silico study with the antiviral effects of various phytoconstituents Borneol (C₁₀H₁₈O), Camphene (C₁₀H₁₆), Camphor (C₁₀H₁₆O), Alpha-thujene (C₁₀H₁₆), Eugenol (C₁₀H₁₄O), Carvacrol (C₁₀H₁₄O) and Parthenolide (C₁₅H₂₀O₃) of feverfew on the viral protein of SARS-CoV-2. Parthenolide shows the best binding affinity with both main protease (M^pro) and papain-like protease (PL^pro). The molecular electrostatic potential and Mulliken atomic charges of the Parthenolide molecule shows the high chemical reactivity of the molecule. The docking of Parthenolide with PL^pro give score of −8.0 kcal/mol that validates the good binding of Parthenolide molecule with PL^pro. This complex was further considered for molecular dynamics simulations. The binding energy of the complex seems to range in between −3.85 to −11.07 kcal/mol that is high enough to validate the stability of the complex. Free energy decomposition analysis have been also performed to understand the contribution of residues that reside into the binding site. Good binding affinity and reactivity response suggested that Parthenolide can be used as a promising drug against the COVID-19.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s11696-022-02067-6.

DFT and MD simulation investigation of favipiravir as an emerging antiviral option against viral protease (3CLpro) of SARS-CoV-2

As per date, around 20 million COVID-19 cases reported from across the globe due to a tiny 125 nm sized virus: SARS-CoV-2 which has created a pandemic and left an unforgettable impact on our world. Besides vaccine, medical community is in a race to identify an effective drug, which can fight against this disease effectively. Favipiravir (F) has recently attracted too much attention as an effective repurposed drug against COVID-19. In the present study, the pertinency of F has been tested as an antiviral option against viral protease (3CL^pro) of SARS-CoV-2 with the help of density functional theory (DFT) and MD Simulation. Different electronic properties of F such as atomic charges, molecular electrostatic properties (MEP), chemical reactivity and absorption analysis have been studied by DFT. In order to understand the interaction and stability of inhibitor F against viral protease, molecular docking and MD simulation have been performed. Various output like interaction energies, number of intermolecular hydrogen bonding, binding energy etc. have established the elucidate role of F for the management of CoV-2 virus for which there is no approved therapies till now. Our findings highlighted the need to further evaluate F as a potential antiviral against SARS-CoV-2.

Wojcik M. Impact of non-covalent interactions on FT-IR spectrum and properties of 4-methylbenzylammonium nitrate. A DFT and molecular docking study

In this research, the impact of non-covalent interactions on the FT-IR spectrum and structural, electronic, topological and vibrational properties of hybrid 4-methylbenzylammonium nitrate (4MBN) have been studied combining B3LYP/CC-PVTZ calculations with molecular docking. 4MBN was synthesized and characterized by using the FT-IR spectrum while the optimized structures in gas phase and in ethanol and aqueous solutions have evidenced monodentate coordination between the nitrate and methylbenzylammonium groups, in agreement with that experimental determined for this species by X-ray diffraction. Here, non-covalent interactions were deeply analyzed in terms of topological parameters (AIM), electron localization function (ELF), localized orbital locator (LOL), Hirshfeld surface and reduced density gradient (RDG) method. Weak interactions such as H-bonds, VDW and steric effect in 4MBN were visualized and quantified by the independent gradient density (IGM) based on the promolecular density. The hyper-conjugative and the delocalization of charge in 4MBN have been elucidated by natural bonding orbital (NBO) while its chemical reactivity was studied and discussed by using molecular electrostatic potential surface (MESP), frontier molecular orbital (FMOs), density of state (DOS) and partial density of state (PDOS). The complete vibrational assignments of 69 vibration modes expected for 4MBN are reported together with the scaled force constants while the electronic transitions were evaluated by TD-DFT calculations in ethanol solution. Thermal analysis (DTA and DSC) was also determined. Molecular docking calculations have suggested that 4MBN presents biological activity and could act as a good inhibitor against schizophrenia disease.

Steric and energetic characterizations of mouse and human musk receptors activated by nitro musk smelling compounds at molecular level: Statistical physics treatment and molecular docking analysis

Understanding olfaction process at a microscopic or molecular level needs more elucidation of the multiple stages involved in the olfaction mechanism. A worth full elucidation and a better understanding of this molecular mechanism, a necessary preamble should be achieved. The content of this work is a preamble for that. A study of the mouse and human olfactory receptors activation in response to two nitro musks stimuli, which are the musk xylol and the musk ketone, are considered here, first, for their wide expanded use in perfumery, but also to show some particular aspects of this process in the case of these two stimuli, which could help to deduce more details and more general aspects in the global olfactory mechanism. A statistical physics modeling using the monolayer model with two independent types of receptor binding sites of the response of the mouse olfactory receptor MOR215-1 and the human olfactory receptor OR5AN1, which are identified as specifically responding to musk compounds, is used to characterize the interaction between the two nitro musk molecules, the mouse and the human olfactory receptors and to determine the olfactory band of these two odorants through the determination of the molar adsorption energies and the adsorption energy distributions. The physico-chemical model parameters can be used for the steric characterization via the calculation of the receptor site size distributions. The docking computation between these two nitro musks and the human olfactory receptor OR5AN1 is performed demonstrating a large similarity in receptor-ligand detection process. Thus, docking finding results prove that the calculated binding affinities were belonging to the spectrum of adsorption energies.

Quantum chemical studies on molecular structure, AIM, ELF, RDG and antiviral activities of hybrid hydroxychloroquine in the treatment of COVID-19: Molecular docking and DFT calculations

Structure-activity relationships for hydroxychloroquine compound and its derivatives resulted in a potent antiviral activity. Where hydroxychloroquine derivatives showed an apparent efficacy against coronavirus related pneumonia. For this reason, the current study is focused on the structural properties of hydroxychloroquine and hydroxychloroquine sulfate. Optimized structures of these molecules have been reported by using DFT method at B3LYP/6-31G* level of theory. The geometric were determined and compared with the experimental crystal structure. The intra and intermolecular interactions which exist within these compounds are analyzed by different methods namely the topological analysis AIM, ELF and the reduced gradient of the density. These approaches make it possible in particular to study the properties of hydrogen bonds. The highest occupied molecular orbital and the lowest unoccupied molecular orbital energy levels are constructed and the corresponding frontier energy gaps are determined to realize the charge transfer within the molecule. The densities of state diagrams were determined to calculate contributions to the molecular orbitals. The molecular electrostatic potential surfaces are determined to give a visual representation of charge distribution of these ligands and to provide information linked to electrophilic and nucleophilic sites localization. Finally, these derivatives were evaluated for the inhibition of COVID-19 activity by using the molecular docking method.

Theoretical and experimental study of guar gum sulfation

The synthesis of guar gum sulfates by a complex of sulfur trioxide with 1,4-dioxane was studied. The influence of temperature, process duration, and the volume of chlorosulfonic acid on the degree of substitution of guar gum sulfates was studied. The sulfation process has been optimized using the Box-Behnken design. It was shown that the optimal conditions for sulfation of guar gum with a complex of sulfur trioxide-1.4-dioxane: temperature 60 °C, duration 2.9 h, and a volume of chlorosulfonic acid of 3.1 ml. Sulfate groups embedding into the structure of guar gum was confirmed by elemental analysis and FTIR. The initial and sulfated guar gum were also characterized by methods: X-ray diffraction, scanning electron microscopy, and gel permeation chromatography. Using X-ray diffraction, it was shown that amorphization of guar gum occurs during sulfation. Using scanning electron microscopy, it was shown that the morphology of guar gum changes in the process of sulfation. Using gel permeation chromatography, it was shown in the process of guar gum sulfation by a complex of sulfur trioxide with 1,4-dioxane, the molecular weight decreases from 600 to 176 kDa. The geometric parameters of all complexes were carried out by using the DFT/B3PW91 method with a 6-31 + G (d,p) basis set. These structures are optimized to predict the important properties of a theme. MEP with contour map has been performed to obtain the electronic properties. Frontier molecular orbital HOMO-LUMO orbital diagram has been obtained for different energy levels and their band gap energies have been computed.

DFT and molecular docking study of chloroquine derivatives as antiviral to coronavirus COVID-19

The recently emerged COVID-19 virus caused hundreds of thousands of deaths and instigated a widespread fear, threatening the world’s most advanced health security. In 2020, chloroquine derivatives are among the drugs tested against the coronavirus pandemic and showed an apparent efficacy. In the present work, the chloroquine and the chloroquine phosphate molecules have been proposed as potential antiviral for the treatment of COVID-19 diseases combining DFT and molecular docking calculations. Molecular geometries, electronic properties and molecular electrostatic potential were investigated using density functional theory (DFT) at the B3LYP/6-31G* method. As results, we found a good agreement between the theoretical and the experimental geometrical parameters (bond lengths and bond angles). The frontier orbitals analysis has been calculated at the same level of theory to determine the charge transfer within the molecule. In order to perform a better description of the FMOs, the density of states was determined. The molecular electrostatic potential maps were calculated to provide information on the chemical reactivity of molecule and also to describe the intermolecular interactions. All these studies help us a lot in determining the reactivity of the mentioned compounds. Finally, docking calculations were carried out to determine the pharmaceutical activities of the chloroquine derivatives against coronavirus diseases. The choice of these ligands was based on their antiviral activities.