key: cord-0835489-jpn3nnct
authors: Salah, M.; Belghiti, M.E.; Aitouna, A.O.; Zeroual, A.; Jorio, S.; El Alaoui Abdellaoui, H.; El Hadki, H.; Marakchi, K.; Komiha, N.
title: MEDT Study of the 1,3-DC Reaction of Diazomethane with Psilostachyin and investigation about the interactions of some pyrazoline derivatives with Protease (M(pro)) of nCoV-2
date: 2020-09-24
journal: J Mol Graph Model
DOI: 10.1016/j.jmgm.2020.107763
sha: b32223ebf4303817f5ddfe277ea6b1c2b6ba2a8a
doc_id: 835489
cord_uid: jpn3nnct

The molecular electronic density theory (MEDT) was invested to elucidate the chemo-, regio- and stereo-selectivity of the 1,3-dipolar cycloaddition between Diazomethane (DZM) and Psilostachyin (PSH). The DFT method at B3LYP/6-31+G(d,p) level of theory was used. Reactivity indices, transition structures theory, IGM and ELF analysis were employed to reveal the mechanism of the reaction. The addition of DZM to PSH takes place through a one-step mechanism and an asynchronous transition states. Eight possible addition channels of reaction were investigated (addition of C(sp2) to Diazomethane at C4, C5, C6 or C7). The addition of C (sp2) at C5 leading to P1 product is the preferred channel . The addition of ether does not affect the chemo-, regio- and stereo-selectivity of the reaction. Analysis of transfer of charges along the IRC path associated with the P1 product shows a polar character for the studied reaction. We have also used the noncovalent interaction (NCI) which is very helpful to reveal the most favored addition channel of the reaction, by analyzing the weak interactions in different TSs. Finally, we investigate about the potential of inhibition of some pyrazoline compounds against COVID-19-M(pro) by performing a molecular docking calculations.

1,3-Dipolar cycloaddition reactions play a key role in organic chemistry [1] . They are routes for the synthesis of a considerable diversity of functionalized five-membered heterocyclic systems, which present a perfect compounds of "non heterocyclic" derivatives [2] . The effect of the 1,3-dipolar cycloaddition reactions in domain of heterocyclic synthesis is comparable to that of Diels-Alder reactions on synthesis of carbocyclic products [3] [4] . In fact, the abundance of various classes of dipoles and dipolarophiles offer a lot of perspectives. Indeed, the dipolar cycloaddition reactions serve as main key in the total synthesis of natural products or analogues [5] . 1,3-Dipolar cycloaddition is a main synthesis strategy to the regio-and stereoselective of five membered heterocycles and their ring opened acyclic derivatives, more than that, it serves as a chemical test [6] . Indeed, the 1,3-dipolar cycloaddition of Diazomethane (DZM) leading to pyrazolines, was used to test the presence of an exocyclic methylene α, β-unsaturated-γ-lactone functionnality [7] . Large experimental and theoretical studies focused at understanding the reaction mechanism and the origin of the observed chemo-, regio-and stereo-selectivity observed in many reactions. The earlier work of Huisgen and coworkers revealed that the mechanism of 1,3-dipolar cycloadditions which lead to 5-membered rings are managed by the stereoselectivity of dipolarophiles, the nature of solvent and substituents, the activation parameters, and by orientation phenomena [8] [9] [10] . Whereas Houk et al. [11] [12] explain that the frontier orbital energies and coefficients allow a specific qualitative treatment of reactivity of individual 1,3 dipoles. Recently, molecular electronic density theory (MEDT) accompagnied with the IGM and ELF analysis have been invested with a great success to reveal the mechanism of many reactions [13] [14] [15] [16] [17] . In the present work, we have used the MEDT, IGM and ELF tools to study the 1,3-dipolar cycloaddition reaction between Diazomethane (DZM) and Psilostachyin (PSH). This reaction has been realized by the Adriana Ortiz-León group [7] in ether solution, leading to mexicanin I acetate pyrazoline product (Scheme 1). However, to date, the stereochemistry of this reaction has not been explained nor a mechanism associated with the previous reaction had been elaborated. The DZM, known as a 1,3-dipole of the propargyl/allenyl anion type presents an electron donating capacity of the Catom with double Carbone-Carbon bonds in dipolarophile, involving its highest occupied molecular orbital (HOMO).

with the PSH reagent.

Theoretically, there are eight possibilities for the C-atom of DZM to add up at the double Carbone-Carbon bonds of PSH reagent (addition of C(sp2) of DZM at C4, C5, C6 or C7). So, this interaction leads a priori to eight products: P1, P2, P3, P4, P5, P6, P7 and P8 (Scheme 2). This reaction evolves through the pathways specified, where eight transition states (TS1, TS2, TS3, TS4, TS5, TS6, TS7 and TS8) are involved. So, our goal in this work is to investigate the different possibilities of addition between the two reagents, to detect the different transition structures, to predict the more favored channel of addition, in order to reveal the mechanism of reaction with a good electronic density description.

The different possibilities of addition of Diazomethane to carbon-carbon double bonds of PSH reagent leading to eight possible products The obtained products P1, P2, P3, P4, P5, P6, P7 and P8 belong to a large family of pyrazoline derivatives (pseudoguaianolide pyrazolines). These compounds are an important biological agents and many studies in literature have been focused on this class, due to their relevant biological properties, that include antiviral, antibacterial, antitumor, anti-tubercular, antifungal, antiparasitic, anti-inflammatory, anti-diabetic, anaesthetic, analgesic and insecticidal agents [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] . Considering the large anti-viral properties of pyrazoline products [28] [29] [30] [31] , they deserve to be tested against COVID-19, seen till this date, there is no vaccine and no antiviral for the treatment of this virus [32] [33] . The appearance of new structural data of key proteins associated with SARSCoV-2 and the host, such as the main protease (M pro ) [34] , the RNA-dependent RNA polymerase (RdRp) [35] , human angiotensin-converting enzyme 2 (hACE2) [36] and spike glycoprotein (S) [37] , make the operation of search of the drug against COVID-19 from the knowledge of the three-dimensional structure of the protein targets and to identify small molecules that exhibit the desired chemical characteristics wanted to bind to the proteins targeted by molecular docking. This should be promising [38] [39] [40] . In this context, The last part of this work aims to investigated about the potential of inhibition of some derivatives of the family of pyrazoline products against COVID-19, by performing a molecular docking calculations. In this paper, we have focused on J o u r n a l P r e -p r o o f main protease protein (M pro ) who intervenes at viral replication and transcription, as key to inhibit this virus [41] [42] [43] [44] [45] [46] . Indeed, the M pro possesses an important role in cutting down polyproteins into functional pieces. Thus, we have targeted this protein in the docking calculations in order to check the ability of our pyrazoline products to inhibit the activity of M pro which could block the replication of the SARS-COV virus.

The Density functional theory (DFT) has become a useful tool with a notable successes in most branches of chemistry [47] [48] [49] . Indeed, many theoretical works on chemical reactivity have adopted this theory for their calculations to providing theoretical insights for chemo-, regio-and stereo-selectivity. In this work, we have performed our DFT calculations by using the GAUSSIAN 09 [50] set of programs along with the graphical interface GaussView-5.0. Geometry optimization of different stationary points were performed at DFT/B3LYP/6-31+G(d,p) level and visualized by CYLview program [51] . Calculations of frequencies have been performed at 298.15°K and 1.0 atm, in order to verify that the transition structures (TSs) possess one imaginary frequency and that the minima have no imaginary frequency [52] . The energy profiles have been developed by the IRC calculation. The electronic structures of stationary points and the Wiberg bond indices [53] were analyzed through the natural bond orbital (NBO) method [54] . The solvent effects of ether were considered using a self-consistent reaction field (SCRF) based on the polarizable continuum model (PCM) of Thomasi [55] . The local and global reactivity have been analyzed through the Parr function using population analysis (NPA) [56] . The reactivity of molecular systems has been unveiled through several concepts [13] , [57] , [58] the electronic chemical potential μ, the chemical hardness η, the global electrophilicity ω and global nucleophilicity N. The previous global indices in the framework of the conceptual Density Functional Theory, are calculated by the next expressions, μ = (ε HOMO +ε LUMO )/2 , η = ε LUMO -ε HOMO and ω = μ 2 /2 η. The global nucleophilicity N, is calculated referring to tetracyanoethylene, N = E HOMO (nucleophile) -E HOMO (TCE). The local indices ω k , N k at atomic site K are obtained using Parr functions, respectively P + (r) and P -(r), as follows ω k = ωP + (r) and N k = NP -(r). P + (r) = ρ s ra (r) and P -(r) = ρ s rc (r). ρ s rc (r) and ρ s ra (r) are respectively the atomic spin density (ASD) of the radical cation and anion. According to Hohenberg-Kohn [59] , the total molecular information is in the electron density (ED) and the models based on the ED are more physically significant for chemical bonding analysis. Indeed, the independent gradient model approach (IGM) [60] , is more efficient to evaluate the evolution of the formation of new bonds. The IGM is more qualified for description of the covalent bonding pattern and revealing the weak interactions where other tools were obviously failing. It offers direct 3-D depiction in the form of closed domains associated with the spatial localization of the interactions. This is done through the mono determinantal wave functions of the ground state (GS) and that of the transition structures calculated at B3LYP/6-31+G(d,p) level, using the Multiwfn software and visualized by the VMD program [61] . The physical descriptor associated with IGM [60] is symbolized by δg, δg=|∇ρ IGM |-|∇ρ|. ρ represents the electron density (it's a local function), and is an experimentally accessible scalarfield, |∇ρ IGM | represents the virtual upper limit of the electron density gradient and |∇ρ| represents a non-interacting system electron density gradient. δg is a measure of the ED gradient collapse caused by electron sharing between contra gradient ED sources present in the chemical system. The ELF calculation [56] , [62] , [63] has been adopted, in order to realize a good description of the electronic structures of our molecular systems and reveal the mechanism of the favored addition along the IRC path. The ELF analysis had a great success in literature. It's defined as follow,

and D 0 (r) = (3 ) ( ) . The D(r) reveals the excess kinetic energy density caused by Pauli repulsion. The D 0 (r) can be interpreted as Thomas-Fermi kinetic energy density, which is the exact kinetic energy density of non-interacting, uniform electron gas. Hence D 0 (r) is introduced as a reference. A large ELF value means that electrons are greatly localized, indicating that there is a covalent bond, a lone pair or inner shells of the atom involved. Finally, we have used the noncovalent interaction (NCI) method [64] , which is also known as reduced density gradient (RDG) method. It will be very helpful to reveal the most favored addition channel, by analyzing the weak interactions in different TSs. ∇ρ(r) (non-interacting system electron density gradient). The region where the electron density is small, reveals the presence of weak interaction region [65] .

Concerning the docking part, the biological activity of our compounds has been evaluated by using Maestro software (Schrödinger, LLC, New York, NY, 2018) [66], Against Main protease (M pro ) of SARS-CoV2 (DOI: 10.2210/pdb6LU7/pdb) [66] which is a potential target of COVID-19 disease. Our compounds were initially docked to the binding site Main protease receptor using Glide SP (Schrödinger Suite 2018) [66] with standard settings. The procedure was composed of the preparation of the ligand and receptor, grid generation, and docking. The binding pocket of the proteins was predicted as to be the nearest space with distance less than 3 Å from ligand crystallized in the Main protease (N3). Concerning the preparation of the ligand in the Schrodinger software, we have used the LigPrep [66] to optimize the geometry of ligand and offer its three dimensional structure with correct chiralities. The LigPrep can also produce a number of structures from each input structure with various ionization states, tautomers, stereochemistry, and ring conformations.

The presentation of the results concerning the reactivity study will be structured as follows. We start byan analysis of the conceptual DFT (CDFT) reactivity indices [67] and a topological analysis of ELF. Secondly, we present a mechanistic study for the eight competitive reaction paths associated with the 1,3-dipolar cycloaddition between DZM and PSH (Scheme 2). In the next we will be focused on analysis of the asynchronicity of the formation of the new sigma bonds using Wiberg bond order and IGM analysis [13] . Then, we show the ELF topological analysis along the IRC path associated with the favored product. This allows an accurate description for the mechanism of our reaction. Finally, we mention an analysis of non conventional interactions (NCI) in the different TSs associated with the addition of DZM at the C4=C5 double bond area, which will be very helpful to reveal the most favored path for addition by analyzing the weak interactions.

To investigate the local atomic properties of reagents we will start by CDFT reactivity indices. Thus, the global reactivity indices defined via the CDFT are a useful theoretical tool to explain the reactivity in cycloaddition reactions [67] [68] . In table 1 are displayed the global properties of the different reagents, the values of the electronic chemical potential, μ, the chemical hardness, η, electrophilicity, ω, and nucleophilicity, N, for the different reagents. Table 1 . HOMO, LUMO energies (eV), electronic chemical potential (µ,eV), chemical hardness (η,eV), global electrophilicity (ω,eV) and global nucleophilicity (N,eV) at B3LYP/6-31+G(d,p) E HOMO (TCE) = -9.41 eV at B3LYP/6-31+G(d,p)

As presented in Table 1 , the electronic chemical potential of DZM (-4.015eV) is higher than that of dipolarophile PSH (-4.635eV), consequently, the transfer of charge will be more favored from dipole to dipolarophile (NED: Normal Electronic Demand). The values of η for the two reagents confirm the previous results of μ, indeed the PSH (η=5.01 eV) is more hard to give electrons than DZM (η=4.65 eV). The electrophilicity ω scale, permits to classifying the organic compounds as strong electrophiles if ω > 1. 5 . Therefore, the values of (ω;Ν) related to the dipolarophile: PSH (2.14;2.27) and dipole DZM (1.73;3.07) indicate that the dipolarophile will act as a strong electrophile and for DZM will act as powerful nucleophile , in agreement with other previous parameters. When a nucleophile approaches an electrophile along a polar procedure, a specific electron density is moved from the nucleophile to the electrophilic species. As a result, the nucleophilic reagent loses a certain electronic density, while the electrophilic reagent gains this electronic density. These changes require a reorganization of the electronic density in two frames. In this case, the electrophilic and nucleophilic Parr functions have proved to be the most precise and insightful tools for studying the redistribution of electron density along polar processes, and are widely used to analyze localities in polar processes reactivity. The electrophilic Parr function of PSH and the nucleophilic Parr function of DZM are represented in Figure 1 . ; the most favorable regioisomeric channel will be linked with the two-center interaction between the most nucleophilic center in the DZM framework, the C1 carbon, and the most electrophilic center of electrophilic PSH C6 carbon, in disagreement with the experimental results and energetic analysis of the different competitive reaction paths. The failure of the Parr indices to predict the more favored zone of addition is explained by the fact that these indices are obtained from separated reagents without any consideration of energetic profile (TSs and products), and the interaction between the two reagents (steric effects and weak interactions). For these reasons we investigate in the next, the mechanistic factors beside the steric effects and the weak interactions associated with our reaction (NCI analysis).

For more details and deep investigation about the local atomic properties of reagents we have realized a topological analysis of ELF and NPA of reagents. Thus, the quantum-chemical evaluation of ELF is an attractive procedure, that provides a simple link between the chemical structures and the distribution of electron density. To describe the simplest electronic structure of reagents, and determine their reactivity in 1,3-dipolar cycloaddition reactions, a topological evaluation of the ELF was first carried out. The ELF basin attractor positions, ELF localization domains as well as the populations of the valence basin, together with the Lewis type structures proposed and the natural atomic charges are illustrated in Figure 2 . (Figure 2) . After realizing an investigation about the local atomic properties of reagents by CDFT and ELF analysis, we have investigated about the kinetic and thermodynamic factors which influence our reaction. Thus we have performed a mechanistic study of our reaction.

Analysis of the mechanistic results in gas-phase indicates that the reaction between DZM and PSH occur along asynchronous concerted processes, therefore the different products (P1, P2, P3, P4, P5, P6, P7 and P8) and their corresponding TSs (TS1, TS2, TS3, TS4, TS5, TS6, TS7 and TS8) were obtained and verified. The total and relative energies are presented in table S1. The geometries of the different TSs are depicted in figure  S1 . The structure of the transition states are characterized by a vibrational frequency analysis. Indeed each transition state possess a single imaginary frequency.

J o u r n a l P r e -p r o o f

According to the relative energies values in gas-phase summarized in Table S1 , Scheme 3a and Scheme 3b, the activation barriers associated with the different possible products P1; P2; P3; P4; P5; P6; P7 and P8 are equals to 13,29; 13,77; 19,09 ; 24,18; 18,64; 18,61; 18,90 and 17 .92 (Kcal/mol) respectively. Consequently, the P1 product has the lowest value of the activation barrier. Kenetically, the P1 product is slightly favored than P2 product, but much more favored than the rest of products, in agreement with the experimental outcomes [7] . Thermodynamically, the most stable product is P2. The difference of stability between P2 (-27,12 Kcal/mol) and P1(-23.29 Kcal/mol) is equal to 3.93 kcal/mol. Hence, the studied reaction is controlled by the kinetic factors. To test the solvent effect in reaction, we have performed a total optimization of all stationary points in ether at the B3LYP/6-31+G(d,p) level by using the self-consistent reaction field (SCRF) model based on the polarizable continuum model (PCM) of Tomasi [55] . All the results of energies are collected in table S1. The addition of ether causes an increase of activation barrier energy for the different products by 0.48 to 2.30 Kcal/mol. The difference between the activation barrier energies of TS1(13.77 Kcal/mol) and TS2(15.39 Kcal/mol) in ether, increase to 1.62 Kcal/mol, so the addition of ether is helpful for the formation of P1 product. The stability of the products in ether, decreases slightly by 0.124 to 1.67 Kcal/mol (Scheme 3). In summary, we can conclude that the addition of ether has a weak influence on the kinetic and thermodynamic factors, without any influence on chemo-, regio-and stereo-selectivity of reaction.

Scheme 3a. The stereo-and regio-isomeric reaction paths associated with the addition of Diazomethane to carbon-carbon double bonds (C4=C5) of PSH. Reagent relative energies at B3LYP/6-31+G(d,p) level, in gaz phase and in ether (in parentheses), are given in Kcal/mol relative to the separated reagents. The thermodynamic factors associated with the reaction of DZM with PSH in ether had been studied at B3LYP/6-31+G(d,p) level by calculating the thermal enthalpy (H) and the thermal free energy (G). The different results of G and H have been mapped at Table S2 . The enthalpy profiles of the reaction paths associated with the 1,3-dipolar cycloaddition of Diazomethane to the PSH reagent, in ether at B3LYP/6-31+G(d,p) (Figure 3) , show that the P1 product need a lowest energy than other products. Indeed it need only 14.96 Kcal/mol as energy for activation and release a low energy of -18.98 Kcal/mol. Generally, all possibilities of reaction leading to different products P1; P2; P3; P4; P5; P6; P7 and P8 release energy in favor of the exterior, which shows that the 1,3-dipolar cycloaddition of Diazomethane to PSH reagentis exothermic. The energy released for different products varies from 16.99 Kcal/mol to 22.02 Kcal/mol.

The optimized geometries of the TSs corresponding to the eights channels of reaction (Scheme 2) are reported in figure S1. The length of the new sigma bonds in evolution C1-C5 and N3-C4 at TS1 (lowest value of activation barrier) are 2.1223 Å and 2.5766 Å respectively. The degree of synchronicity (∆d, in Å) can be measured as the difference between the lengths of the two new sigma bonds in evolution. The value of ∆d between the C1-C5 and N3-C4 bonds in evolution of TS1 is equal to 0.4543 Å, indicating that the process of formation of the new sigma bonds in evolution occurs asynchronously. For TS2, the values of length of the new bonds in evolution are d(C1-C5) = 2.24399 Å and d(N3-C4) = 2.25172 Å, hence the formation of the two bonds is more or less synchronous. For more details about the degree of synchronicity ∆d for the rest of TSs, see Table  S3 . To confirm the previous results we use, in the next, the Wiberg bond indices and IGM analysis.

The results of TVs allow us to collect useful information about the chemical process linked to the different transition structures detected in the studied reactions. The imaginary frequencies, the main TVs components and their corresponding geometric parameters for TSs are presented in Table 2 . J o u r n a l P r e -p r o o f the Diels Alder cycloadditions (500 cm −1 ), indicating that these processes are associated with the motions of the heavy atoms and are also related to the earlier transition states.

Beside the wiberg bond indices which have been computed by using the NBO population analysis, the novel IGM descriptor δg=|∇ρ IGM |-|∇ρ|, can be also invested as a real and an effective physical index, based on the electron density to obtain a deep analysis of the evolution of the cycloadducts process. The IGM analysis can be obtained by evaluating the electron density using the corresponding mono-determinantal wave functions of the ground state, corresponding to the TSs. The δg descriptor measures the interaction between atoms in the system, so we have focused on the atoms that will participate in the formation After detecting the different TSs and analyzing their geometries and frequencies, we have then elaborate the mechanism of our reaction. Hence the necessity to an ELF topological analysis along the IRC path associated with the favored product (P1).

The ELF topological analysis of sorted structures along the IRC profile are listed in table S5. The most significant points of the IRC path were throughout analyzed and characterized, which lead to obtaining an important information concerning the mechanism of the reaction. The absence of an intermediate complex along IRC (Scheme S1), indicates that the formation of the new sigma bonds N3-C4 and C1-C5 takes place through a one-step molecular mechanism. Topological analysis of the ELF highlight that the formation of the N3-C4 and C1-C5 bonds is done in two stages and one-step. The ELF analysis reveals an interesting valence basins from S0 to S18 (see supplementary).

The topological analysis of the ELF of the first structure of the IRC (S0), displays similar bonding pattern to those of the seperated reagents. Indeed, at S0 there are two important disynaptic basins V(C,C)=8 and V(C,C)=9 with a sum of electronic population equal to 3.45e, indicating that C4-C5 is a double bond, in front of this, we have two other disynaptic basins V(N,N)=3 and V(N-N)=4 integrating a total of population of 3.59e, which indicates a double bond between the two N atoms of DZM fragment. Furthermore, we have a disynaptic basin V(N,C)=5 between N atom and C atom of DZM integrating 3.09e which implied a pseudo double bond between the two previous atoms. Concerning the most important monosynaptic basins we have V(N)=1; V(N)=2; V(C)=6 and V(C)=7 which integrate 2.04e; 1.83e; 0.50e; 0.47e respectively. The values of V(N)=1; V(N)=2 indicate the existence of a lone pair and negative charge (-1.83e) on N atom. The change in basins is more remarkable at the structure S3, where we observe the disappearance of V(C)=6 and V(C)=7 basins and appearance of a novel basin V(C)=e with 0.72e, which increases to 1.49e at S6. In the same time the two disynaptic basins V(C-C)=8 and V(C-C)=9 disappear leading to a novel disynaptic basin V(C-C)=13 at S4 including 3.36e. At S11 we notice the formation of the novel C-C sigma bond (1.75e) which increases its electronic population to 1.85e at S18. Consequently the C-C bond is more advanced than the N-C sigma bond. The next novel sigma bond (N-C) begin to be formed at S15 by the appearance of monosynaptic basin V(C)=17 integrating 0.75e , which increases to 0.87e at S16, leading then to appearance of monosynaptic basin V(C)=14 and a disynaptic basin V(C-C)=18 at S17 integrating 1.58e as electronic population which increases to 1.71e at S18 leading to the formation of the second new sigma bond N3-C4.

Analysis of charge transfer along the IRC path associated with the P1 product, indicates that it takes place from DMZ fragment to PSH fragment in conformity with the results of CDFT reactivity indices, so we confirm the NED character of the reaction. The values of global electronic density (table S5) transferred from DZM to PSH varied between 0.044e and 0.205e, which indicates the polar process of our reaction.

As we have already reported the local indices have failed to predict the more favored zone of addition. This could been explained by the fact that these indices are obtained from separated reagents without any consideration of steric effects and weak interactions. For these reasons we have investigated about the steric effects and the weak interactions associated with our reaction by realizing a non covalent interactions analysis. Thus, the NCI index developed by Yang's group using quantum-mechanical electron density [69] is a visualization index, which allow an identification of non-covalent interactions (favorable and unfavorable interactions). It provides a rich representation of the weak interactions where other methods fail to describe them. The visualization for different intramolecular and intermolecular interactions (attractive and repulsive forces) is performed by coloration ( Figure 5) . The more blue implies the stronger attractive interaction, the interaction region marked by green can be identified as VdW interactions region and the regions marked by red, implied a zone of strong repulsion effect [58] . In our case, we have focused on the impact of the weak interactions on the regio-and stereo-selectivity at C4=C5 double bond area, so we have performed the NCI calculations with NCI plot [70] for the different TSs,TS1; TS2;TS3 and TS4 ( figure 6) . The visualization has been realized by the VMD program [61] . Herein, the 2D scatter graph and 3D RDG isosurface for the different TSs are represented in figure 6 and figure 7 . The RDG and sign (λ 2 )ρ are a pair of very important functions for revealing weak interaction regions, they are collectively employed in NCI method. The scatter plots of RDG versus the second largest eigenvalue of the Hessian matrix of electron density functions are obtained by using Multiwfn 3.6 program [71] , based on the obtained wave functions calculated at the B3LYP/6-31+G(d,p) level and plotted using gnuplot program [72]. The figure 6 presented the different scatter graphs associated with the eight TSs. These scatter graphs, inform us about the type of dominant interactions by calculation of RDG in function with the sign(λ 2 )ρ. According to the obtained Scatter graphs, it exists tree type of interactions: strong attractions (blue area), Van Der Waals interactions (VdW) (green area) and strong repulsion (red area). The scatter graph of TS1 presents a large blue area (-0.05 a.u to -0.02 a.u), a narrow green area (-0.02 a.u to 0.01 a.u) and a red area (0.01 a.u to 0.05 a.u). Generaly, for the rest of scatter graphs associated with other TSs we remark that the green area expands slightly at the expense of blue area, furthermore the red area gets slightly bigger. The analysis of color-mapped RDG isosurfaces ( figure 7) show a strong attraction (blue color) of N and C atoms of DZM with the two C atoms (Double bond C=C) of PSH in different TSs. However, the addition approach associated withTS1 is more favored because the DZM fragment can get near the double bond easily without any other effects that can hinder it. On the other hand, for the rest of the TSs we can observe more VdW interactions (green color) which weaken the forces of addition on the double bond. So, the approach of DZM to the double bond of PSH in these cases is more difficult for the rest of the TSs. Consequently, the P1 pathway is more favored than the other products in conformity with the experimental outcomes [7] . 

The outbreak of the novel coronavirus disease COVID-19, caused by the SARS-CoV-2 virus, has spread rapidly around the world, the epidemiological characteristics of this virus allow it a high contagiousness risks. This invincible virus, causes enormous damage every day across the globe, and destabilize the human society and economy. As today, there is no vaccine and no antiviral for the treatment of this virus. However, some compounds can been promising for treatment, as pyrazoline molecules. So, we have tested some compounds of pyrazoline family (Table 3) with the Main protease (M pro ) of SARS-CoV2 which is an essential protein for the proteolytic maturation of the virus. Indeed, it has been examined as a potential target protein to prevent the spread of infection by inhibiting the cleavage of the viral polyprotein [73]. 

Molecular docking study is a helpful cheminformatics tool [41] [42] [43] that gives compounds interacting with the protein (M pro ), related to a specific biological activity, the manner in which they interact, and also the way a ligand (pyrazoline compounds) and M pro protein can bind. The crystal structure of COVID-19 main protease in complex with an inhibitor N3 has been realized by Jin Zheming et al. [46] the global symmetry of the complex is Cyclic-C2 and its global stoichiometry is Homo 2-mer-A2 (Figure 8) . with standard settings. We choose the compounds who have a binding affinity higher than -5kcal/mol for further analysis. We evaluate their interaction with the binding site of the target.

The 3D interacting site of the ligands with the binding site of the protein are in the figure 9 , the Hydrogen bond interactions for all the compounds who have a binding activities less than -5kcal/mol with the protein are shown in figure 10 . The molecular docking scores of the compounds and the amino acids found in the active site pockets of M pro protein are showed in the table 4. Main protease has the hydrogen bond attractors and hydrogen bond donors with the amino acids in the binding site at a distance less than 2.73 Å. The 2D structure of the most active compounds interacting with the amino acids in the binding site of the Main protease show for complex Spiro G-M pro a hydrogen bond donor between oxygen atom of Spiro-G and CYS145, GLU166, hydrogen bond acceptor with the hydoxyd group (OH) of Spiro-G and ASN142. For Spiro-F, there are two hydrogen bonds donors with oxygen atom and CYS145 , GLN189, hydrogen bond donor between nitrogen atom and GLY143. For the Spiro C there are two hydrogen bonds between oxygen atom and SER144, THR26. Concerning the P1 compound there are three hydrogen bonds donors between oxygen atom and GLN189, GLY143, GLU 166. Finally, for P3, there are two hydrogen bonds donors between oxygen atom and GLU166, CYS145. We noted also, that the most common amino acids in the binding site for the different complexes are CYS145, GLU166 and GLY143. The comparison of the binding energy of our products (Spiro-C, Spiro G and P3) is greater than that of Lopinavir and the mixture of Lopinavir+Oseltamivir which have been used to cure MERS and SARS [73]. For all these compounds must be predicted their toxicity, druglikness and ADME (Data of ADME/ToX) using the preADMET server [74] . According to results associated with toxicity, druglikness and ADME we note the following. For the different pyrazolline products the Ames_test is non-mutagen. The risk associated with the hERG_inhibition is low and their CMC_like_Rule is qualified. The solubility of P1, P3, spiro-C, spiro-F and spiro-G in pure water is between 71-144 mg/l. The value of binding with the plasma protein is from 69.94 to 76.80, consequently, these products possess a good absorption. For more details about the prediction of ADME/Tox see suplementary (Data of ADME/ToX). In summary, the pyrazoline compounds may act as inhibitors of the COVID-M pro , and they deserve more investigations. 

The B3LYP/6-31+G(d,p) level of computation was used to study the 1,3-DC reaction of diazomethane to PSH. The kinetic, thermodynamic and geometrical parameters associated with the formation of the different possible products were analysed in the gas phase and ether. The global and local reactivity indices at the ground state of the reagents were evaluated to understand the reactivity in the polar processes. Whereas the electrophilic and nucleophilic character of the reagents is are responsible for the high polar character of this 1,3-DC reaction, the energetic results for 1,3-DC reaction of diazomethane to PSH, indicate that the reaction path, leading to P1 product, is kinetically preferred. The addition of ether as a solvent has a weak influence on the kinetic and thermodynamic of the studied reaction without any influence on the chemo-, regio-and stereo-selectivity. ELF calculations reveal that the formation of the N3-C4 and C1-C5 bonds is done in two stages and one-step and the C1-C5 bond is more advanced than N3-C4 bond. The weak interactions associated with the VdW forces unfavored the addition of DZM in the cumbersome zone by weakening the forces of addition, leading to the formation of P1 product. Finally, the investigation about the inhibition of Mpro-COVID-19 by pyrazoline products have given encouraging results. More accurate researches on this family of products are worth undertaking for the struggle against COVID-19 virus.

Supplementary data. Table S1 . Energies of different stationary points and relative energies of TS and products in gas phase and solvent at B3LYP/6-31+G(d,p) Table S2 . Thermochemistry parameters H and G in (u.a) of different stationary points in ether Table S3 . Geometric parameters of the new sigma bonds (lengths in Angstrom) of TSs in gas phase and ether. Table S4 . Wiberg bond orders in TSs and products, percentage of formation and IGM values δg in TSs associated with the new sigma bonds . Table S5 . Populations calculated from the ELF of some selected points, along the IRC curve in gaz phase at B3LYP/6-31+G(d,p) associated with the channel leading to the formation of P1 product of reaction between DZM and PSH. GEDT (global electronic density transferred) from DZM to PSH. 

An efficient stereoselective synthesis of 3-spirocyclopentene-and 3-spiropyrazole-2-oxindoles via 1,3-dipolar cycloaddition reaction

Modern Heterocyclic Chemistry, 1re éd

Isoxazoline-carbocyclic aminols for nucleoside synthesis through aza-Diels-Alder reactions

Carbocyclic Ring Construction via an Intramolecular Diels-Alder Reaction of an in Situ-Generated, Heteroatom-Stabilized Allyl Cation: Total Synthesis of (±)-Lycopodine

Current Trends towards the Synthesis of Bioactive Heterocycles and Natural Products Using 1,3-Dipolar Cycloadditions (1,3-DC) with Azomethine Ylides

Theoretical Prediction and Experimental Tests of Conformational Switches in Transition States of Diels−Alder and 1,3-Dipolar Cycloadditions to Enol Ethers

The stereochemistry of the 1,3-dipolar cycloadditions of diazomethane to pseudoguaianolides

-Dipolar Cycloadditions. Past and Future

Concerted nature of 1,3-dipolar cycloadditions and the question of diradical intermediates

Kinetics and Mechanism of 1,3-Dipolar Cycloadditions

Origin of reactivity, regioselectivity, and periselectivity in 1,3-dipolar cycloadditions

Frontier molecular orbitals of 1,3 dipoles and dipolarophile

Theoretical study of the 1,3-DC reaction between fluorinated alkynes and azides: Reactivity indices, transition structures, IGM and ELF analysis

Reactivity Indices related to DFT Theory, the Electron Localization Function (ELF) and Non-Covalent Interactions (NCI) Calculations in the Formation of the non-Halogenated Pyruvic Esters in Solution

A molecular electron density theory study of the enhanced reactivity of aza aromatic compounds participating in Diels-Alder reactions

A molecular electron density theory investigation of the molecular mechanism, regioselectivity, stereoselectivity and chemoselectivity of cycloaddition reaction between acetonitrile Noxide and 2,5-dimethyl-2H

Deciphering the role of acid additives in chiral diaminecatalyzed asymmetric aldol reactions of cyclohexanones with aldehydes

Biological Activities of Pyrazoline Derivatives-A Recent Development

Antifungal and antimycobacterial activity of 1-(3, 5-diaryl-4,5-dihydro-1H-pyrazol-4-yl)-1H-imidazole derivatives

Synthesis and antimicrobial studies on novel chloro-fluorine containing hydroxy pyrazolines

Bis-pyrazolines: Synthesis, characterization and antiamoebic activity as inhibitors of growth of Entamoeba histolytica

Synthesis and antiinflammatory activity of some new 1,3,5-trisubstituted pyrazolines bearing benzene sulfonamide

Design, synthesis, and biological evaluation of 1-(4-sulfamylphenyl)-3-trifluoromethyl-5-indolyl pyrazolines as cyclooxygenase-2 (COX-2) and lipoxygenase (LOX) inhibitors

Synthesis of novel thiazolone-based compounds containing pyrazoline moiety and evaluation of their anticancer activity

Design, synthesis, and biological testing of pyrazoline derivatives of combretastatin-A4

Synthesis and evaluation of pyrazolignans. A new class of cytotoxic agents

Synthesis and antidepressant activity of some 1,3,5-triphenyl-2-pyrazolines and 3-(2-hydroxy naphthalen-1-yl)-1,5-diphenyl-2-pyrazolines

Triaryl Pyrazoline Compound Inhibits Flavivirus RNA Replication

Ecofriendly one-pot synthesis and antiviral evaluation of novel pyrazolyl pyrazolines of medicinal interest

Anti-HIV-1 and cytotoxicity studies of piperidyl-thienyl chalcones and their 2-pyrazoline derivatives

Synthesis of 5-(4-hydroxy-3-methylphenyl)-5-(substituted phenyl)-4, 5-dihydro-1H-1-pyrazolyl-4-pyridylmethanone derivatives with anti-viral activity

COVID-19): current status and future perspectives

Ongoing Clinical Trials for the Management of the COVID-19 Pandemic

The crystal structure of COVID-19 main protease in complex with an inhibitor 11b

Crystal structure of the RNA-dependent RNA polymerase of a human picorbirnavirus

Structure of SARS-CoV-2 chimeric receptor-binding domain complexed with its receptor human ACE2

Crystal structure of the SARS-CoV Spike protein fusion core

A search for medications to treat COVID-19 via in silico molecular docking models of the SARS-CoV-2 spike glycoprotein and 3CL protease

In silico studies on therapeutic agents for COVID-19: Drug repurposing approach

Molecular docking, simulation and MM-PBSA studies of nigella sativa compounds: a computational quest to identify potential natural antiviral for COVID-19 treatment

Computational screening of antagonists against the SARS-CoV-2 (COVID-19) coronavirus by molecular docking

Identification of potential inhibitors of SARS-CoV-2 main protease from Aloe vera compounds: A molecular docking study

Identification of phytochemical inhibitors against main protease of COVID-19 using molecular modeling approaches

Identification of potential drug candidates to combat COVID-19: a structural study using the main protease (mpro) of SARS-CoV-2

Guo Emerging coronaviruses: Genome structure, replication, and pathogenesis

Stucture of Mpro from SARS-Cov-2 and discovery of its inhibitors

DFT computations and molecular dynamics investigations on conformers of some pyrazinamide derivatives as corrosion inhibitors for aluminum

A DFT study on the adsorption of benzodiazepines to vermiculite surfaces

A theoretical study on the regioselectivity of 1,3-dipolar cycloadditions using DFT-based reactivity indexes

Gaussian 09, revision A

CYLview User Manual

Comparison between Optimized Geometries and Vibrational Frequencies Calculated by the DFT Methods

The Mayer bond order as a tool in inorganic chemistry

Natural bond orbital methods

Polarizable Continuum Model (PCM) Calculations of Solvent Effects on Optical Rotations of Chiral Molecules

Understanding the local reactivity in polar organic reactions through electrophilic and nucleophilic Parr functions

CDFT-Based Reactivity Descriptors as a Useful MEDT Chemoinformatics Tool for the Study of the Virotoxin Family of Fungal Peptides

Computational study of the 1,3-dipolar cycloaddition between methyl 2-trifluorobutynoate and substituted azides in terms of reactivity indices and activation parameters

Density-functional theory beyond the Hohenberg-Kohn theorem

The Independent Gradient Model: A New Approach for Probing Strong and Weak Interactions in Molecules from Wave Function Calculations

Origin of the synchronicity in bond formation in polar Diels-Alder reactions: an ELF analysis of the reaction between cyclopentadiene and tetracyanoethylene

Understanding the regioselectivity in hetero Diels-Alder reactions. An ELF analysis of the reaction between nitrosoethylene and 1-vinylpyrrolidine

Accurately extracting the signature of intermolecular interactions present in the NCI plot of the reduced density gradient versus electron density

Utilizing Hirshfeld surface calculations, non-covalent inter-action (NCI) plots and the calculation of inter-action energies in the analysis of mol-ecular packing

Schrödinger Release 2018-4: Prime, Schrödinger, LLC

Schrödinger Release 2018-4: Epik, Schrödinger, LLC

Applications of the Conceptual Density Functional Theory Indices to Organic Chemistry Reactivity

DFT-Based Regioselectivity Criteria for Cycloaddition Reactions

Analysis of Hydrogen-Bond Interaction Potentials from the Electron Density: Integration of Noncovalent Interaction Regions

NCIPLOT: a program for plotting non-covalent interaction regions

Multiwfn: a multifunctional wavefunction analyzer

Computational studies of drug repurposing and synergism of lopinavir, oseltamivir and ritonavir binding with SARS-CoV-2 Protease against COVID-19

We would like to thank MENFPESRS-DESRS and UCD for their financial support of this work.J o u r n a l P r e -p r o o f