Geometry optimization
The structure of the title compound Mhadhbi et al. [15], crystallizes in the centrosymmetric triclinic system with P space group; the unit cell parameters are a = 6.570 (5) Å, b = 8.528 (2) Å, c = 12.529 (3) Å, α = 95.83 (3)°, β = 97.149 (4)°, γ = 95.116 (7)° and V = 689.2 (6) Å3. The unit formula (Fig. 1a) is made of two 4-aminopyridinium cations mono and diprotonated ([C5H7N2] +, [C5H8N2]2+), an arsenate anion [AsO4]3- and one water molecule. The different components are connected by a network of N/O–H···O hydrogen bonds. In order to find the optimized geometry, theoretical calculations were carried out by the DFT/B3LYP with LanL2DZ and 6-311++G(d,p) basis sets. The experimental and calculated molecular geometry parameters (bond lengths, bond angles of the title compound) were given in Table S1, showing the observed and calculated bond parameters are compatible. The Fig. 1b represents the corresponding optimized geometry by the DFT calculation. The difference between the optimized and experimental geometries was analyzed in term of Root Mean Square Deviation (RMSD) overlay. The RMSD fits of the bond lengths and bond angles of experimental and calculated geometries are 0.026 Å and 2.946° respectively, which indicate the agreement between the two geometries. As shown in Table S1, a set of calculated geometrical parameters are slightly deviated from the experimental values; this deviation might be due to that fact that the calculations are carried out for a single molecule in the gaseous phase and the experimental results are obtained in the solid state.
The details of the molecular bond parameters are briefly summarized below. The experimentally obtained values of the C-N and C-C bond lengths are in the range of 1.311–1.390 and 1.363–1.423 Å, respectively. It can be seen that the N-C and C-C bond lengths predicted by the B3LYP basis are found to be between 1.341–1.472 and 1.373–1.420 Å, respectively. The experimental values of the C-C-C bond angles in the phenyl rings and N-C-C angles vary in the range of 117.038–120.023° and 117.169–122.805°, respectively. However, the corresponding theoretical values are in the range of 117.441– 121.309° and 119.303–122.945°, respectively. As for the anionic group, each arsenic atom at the center of the tetrahedron is connected to four oxygen atoms. The experimental values of the large As-O bond lengths are 1.631–1.651 Å, whereas the theoretical values of the corresponding bond lengths are 1.719–1.828 Å. On one hand, the experimental O-As-O angles vary between 108.239 and 111.234°. These values are in agreement with the theoretical corresponding angles which vary between 100.085 and 115.848°, respectively. The molecular geometry in gas phase may differ from the solid phase owing to the extended hydrogen bonding and stacking interactions (attractive non-covalent interactions between two aromatic rings). The difference between the theoretical and experimental geometry might be to the fact that calculations were performed using isolated molecule in the gaseous phase to obtain theoretical results and in solid state for experimental results. Thus, it is found that most of the optimized bond lengths and the bond angles are in reasonable agreement with the corresponding experimental values. As seen, the calculated geometric parameters represent good approximation and can be used as a foundation to calculate the other parameters for the compound. There are slight differences since the experimental results belong to the solid phase, whereas the calculated values are related to the gas phase.
Molecular electrostatic potential (MESP) maps
The molecular electrostatic potential (MEP) surface was drawn by the DFT method. The electrostatic potential is used to exhibit the molecular outline, dimension, and electrostatic potential in terms of color grading [31]. The potential is presented in the map of MEP and visualized with different colors as shown in Fig. 2. At DFT level the MEP surface is plotted for the title compound. The color code of these maps ranges from -2.07 10-1 (red) and 2.07 10-1 a. u. (blue). The mapping represents the donor atoms in blue regions with positive potential; the most positive electrostatic potential (nucleophilic site) is localized hydrogen atoms. The acceptor atoms with negative potential in red regions represent the most electronegative potential (electrophilic sites) and it is mainly localized over anionic group. The green or yellow color indicates the neutral region. This graph shows that the electrophilic sites are located around the hydrogen atoms, whereas the nucleophilic sites are located around the oxygen atoms. MEP surface represents a proof of the creation of hydrogen bond [32]; it is defined as a donor acceptor interaction. The negative region plays an important role in forming a hydrogen bond interaction with the protein. This is well correlated with molecular docking study as well as molecular dynamics study. The electrostatic potential of the donor becomes less negative as the acceptor becomes more negative. The contour map of electrostatic potential confirms that the different negative and positive potential sites of the molecule are in accordance with the total electron density surface map.
Natural Bond Orbital (NBO) Analysis
The strength of the interaction between electron donors and the electron acceptors, or the donating tendency which is from electron donors [33] to electron acceptors and hence the degree of conjugation of the system is measured by the value of energy of the hyperconjugative interactions, E(2). In order to evaluate the donor–acceptor interactions, the second-order Fock matrix was carried out in NBO analysis [34]. For each donor (i) and acceptor (j), the stabilization energy E(2) associated with the delocalization i→j is estimated as:
E(2) = qi F(i,j) /( εi - εj)
Where qi is the donor orbital occupancy, εi and εj are diagonal elements and F(i,j) is the off diagonal NBO Fock matrix element between the natural bonding orbitals.
Delocalization of electron density between occupied Lewis type (bond or lone pair) NBO orbitals and formally unoccupied (antibonding and Rydberg) non- Lewis NBO orbitals correspond to a stabilizing donor–acceptor interaction. The NBO analysis has been performed on the compound using NBO 3.1 program as implemented in the Gaussian 09W package in order to elucidate the intramolecular interaction, re-hybridization and delocalization of electron density within the molecule are presented in Table S2. For the title compound, the stabilization energy E(2) values are chemically important and are used to estimate the intramolecular delocalization. Interactions between the bonding orbital s(C15-C16), s(C15-C24) and s(C22 - C24) and the anti- bonding orbital s*(C15 - C24), s*(C15- C16) and s*(C15-N 26) possess medium stabilization energy roughly equal to 5.5 kcal/mol. In addition, the interactions between the p (C4 - C6), p (C9 - C11), p (C16 - C18) and p (C22 - C24) bonding orbitals and antibonding orbitals LP*(1) C8, LP*(1) C8, LP*(1) C15 and LP*(1) C15 lead to an average stabilization energy around 52.5 kcal/ mol. Similarly, the strongest stabilization energies are occurred in p (C4 - C6)→LP (1) N1 and p (C9 - C11)→LP (1) N1 with 118.91 and 119.09 kcal/ mol, respectively. Furthermore, the resonance interaction energies for the donor-acceptor excitations of the LP (N) bonding orbitals to p* antibonding orbitals are LP (1) N1→p*(C4 - C6) 48.1 kcal/ mol, LP (1) N1→p*(C9 - C11) 47.51 kcal/ mol, LP (1) N 20→p*(C16 - C18) 42.93 kcal/ mol and LP (1) N 20→p*(C22 - C24) 41.07 kcal/ mol. On the other hand, the NBO analysis of ((4-APH)(4-APH2)[AsO4]·H2O) clearly explains the evidences of the formations of strong interaction between LP(N) and σ* (As - O) anti-bonding orbitals and are given in Table S2. The strongest stabilization energies are occurred in s*(As30 - O31)→s*(As30 - O33) and s*(As30 - O34)→ s*(As30 - O33) with 165.72 and 167.76 kcal/ mol, respectively. NBO analysis confirms the existence of X-H…Y hydrogen bonds between the different interactions donor-acceptor, as, for example, the hyperconiugative interactions energy between O31-H2 involving the lone pair (LP(1), and LP(3)) on the oxygen atom (O31). The stabilization energy of [LP(1)O31-σ*(N1-H2)] and [LP(3)O31-σ*(N1-H2)] is 8.56 and 12 kcal/ mol, respectively. This result may be justified by the extent of overlap between the lone pair orbitals (LP (1) and LP (3)) of O31 and the σ*(N1-H2). According to NBO results, the strong donor-acceptor interactions occurred in the title compound implies good nonlinear optical properties. Finally, we can conclude that the NBO analysis confirms that the hydrogen bonding in the (4-APH)(4-APH2)[AsO4]·H2O compound leads to its crystal structure.
AIM approach: topological parameters at bond critical points (BCP)
In the AIM theory the nature of the bonding interaction can be determined through an analysis of the properties [35] of the charge density, ρ, and its Laplacian ∇2ρ, local potential energy density (V(r)), local gradient kinetic energy density (G(r)), total energy density (H(r) = (V(r) + G(r)) and the bond energy E at the BCPs, and through the properties of the atoms, which are obtained by integrating the charge density over the atom orbitals [35]. In addition, the ellipticity (ε) at the BCP is a sensitive index to monitor the π -character of bond. The ε is related to λ1 and λ2 which correspond to the eigen values of Hessian and is defined by the relationship: ε = (λ1/ λ2)-1. Hence, the AIM graph displaying all the BCPs in the asymmetric unit is presented in Fig. 3. The topological parameters calculated at BCPs are detailed in Table 1. According to the BCP analysis of Table 1 and Based on the Rozas et al. [36] criterion all of the seven hydrogen bonds are considered weak since the Laplacian and the energy density values are positive. The electron density ρ is low and its corresponding Laplacian ∇2ρ is positive at all BCPs and are in the range of 0.0040 - 0.0367 a.u. and 0.0135 - 0.1217 a.u., respectively. According to Koch and Popelier criteria [37], the hydrogen bonding interaction should have ρ and ∇2ρ values within the range 0.0020 - 0.0400 and 0.0240 - 0.1390 a.u, respectively. Referring to Table 2, the most strong interaction is found for O31···H2 contact (ρ = 0.0367 and ∇2ρ = 0.1217 a.u.) with an energy that is equal to 41.78 kJ. mol-1. In order to investigate the effect of π-electron delocalization in bonds associated with N and O atoms of N-H···O, the analysis of the bond ellipticity is performed. The ellipticities (ε) of bond of aromatic ring at BCP are in range of 0.028-1.629. The ε values confirm the presence of resonance assisted intermolecular hydrogen bonds.
Noncovalent interaction-Reduced Density Gradient (RDG) analysis
To prove the existence of weak interactions [38,39], NCI plots are generated with the plots of the RDG S versus (sign λ2)ρ, where (sign λ2)ρ is the electron density multiplied by the sign of the second Hessian eigenvalue (λ2). The value of (sign λ2)ρ is useful to predict the nature of interaction; for a repulsive interaction (sign λ2)ρ > 0 (steric effect in ring and cage) and for attractive interaction (sign λ2)ρ < 0 (hydrogen bonds). The values close to zero indicating Van der Waals interactions.
To study non-covalent interactions, Johnson et al. [40] developed an approach based on the electronic density and its derivatives named RDG. The RDG were defined using the following equation.
The strength and nature of interactions can be interpreted from the product sign λ2 times ρ, as blue, green, and red color codes are used to describe stabilizing H-bonding, Van der Waals, as well as destabilizing steric interaction, respectively. The results were collected in Fig. 4a and illustrate the attractive, Van der Waals and repulsive interactions in the present structure. The NCI-RDG graph of the title compound is plotted in Fig. 4b. The observed red areas (spots) indicate steric repulsions and are localized mostly in the middle of benzene rings. The observed green and blue areas between the organic and the inorganic parts indicate the existence of Van der Waals and steric effects.
Frontier Molecular Orbital Analysis
In order to evaluate the energetic behavior [41], we carried out the calculations in gas phase and in solvent. The calculated energy values of HOMO and LUMO are around -6.60 eV and 1.40 eV, respectively, in the gas phase and in water. The energies and the pictorial illustration of HOMO, LUMO, HOMO-1 and LUMO+1 frontier molecular orbitals are shown in Fig. 5. The brown color refers to the positive phase, which indicates the nucleophilic site. The negative phase is indicated by green color, and is related to the electrophilic site. Accordingly, the energy gap of HOMO-LUMO explains the eventual charge transfer interaction within the molecule. Furthermore, in going from the gas phase to the solvent phase, the increasing value of the energy gap and molecule becomes more stable and the low values of frontier orbital gap in (4-APH)(4-APH2)[AsO4]·H2O make it more reactive. Additionally, a molecule with a small frontier orbital gap is readily polarizable and normally exhibits high chemical reactivity and low kinetic stability [42-44]. Furthermore, the global chemical reactivity descriptors (GCRD) are another vital technique to understand the chemical properties of a molecule such as chemical hardness (η), chemical potential (μ), chemical softness (S), electronegativity (χ) and electrophilic index (ω). The formula (Eq. 1 to Eq. 5) to obtain GCRD is derived from HOMO-LUMO energies, taking the HOMO energy as ionization potential (I) and LUMO as electron affinity (A) and the data are tabulated in Table 2. The concept of these parameters is related to each other [34,45-47] where:
Chemical hardness (η) = [(ELUMO - EHOMO) / 2] Eq. 1
Chemical softness (S) = [1/2η] Eq. 2
Chemical potential (μ) = [(EHOMO + ELUMO) / 2] Eq. 3
Electronegativity (χ) = [(I + A) /2] Eq. 4
Electrophilicity index (ω) = [μ2 /2η] Eq. 5
Using the above equations, the chemical potential, hardness and electrophilicity index have been calculated and their values are shown in Table 2. The ionization potential value indicates that energy of 6.60 eV in gas phase is required to remove an electron from the HOMO. The lower value of electron affinity shows higher molecular reactivity with the nucleophiles. Higher hardness and lower softness values confirm the higher molecular hardness associated with the molecule. The electrophilicity index helps in describing the biological activity of the compound.
Non-linear optical effects (NLO)
The application of (4-APH)(4-APH2)[AsO4]·H2O in the field of non linear optics was investigated through the hyperpolarizability calculation by analyzing the vibrational modes of the optimized molecular geometry. The calculated values of dipole moment, polarizability and first order hyperpolarizability and their components in gas phase and water solvent are listed in Table 3. However, the values of α and β from Gaussian output are given in atomic units (a.u.) and converted into electronic units (esu) (α; 1 a.u. = 0.1482 ×10−24 esu, β; 1 a.u. = 8.6393×10−33 esu). The total dipole moment (μt) and first order hyperpolarizability (βt) are found to be 39.1985 Debye and 2205.2506×10−33 esu, respectively. According to the NLO results, the crystal has a total dipole moment of 28.54 times greater than that of urea. The first order hyperpolarizability is more than 5.91 times greater than that of urea (μ and β of urea are 1.3732 Debye and 0.3728×10−30 esu) and 3.22 times that of the reference crystal KDP (βKDP = 6.85×10-31 esu) [48]. In addition, it is associated with intra-molecular charge transfer that is attributed to electron cloud movement through π - conjugated framework from electron donor to electron acceptor groups. Therefore, the components of hyperpolarizability are useful to understand charge delocalization in the molecule. The maximum charge delocalization occurs along βxxx in the molecule under investigation. According to this result, this compound may exhibit excellent nonlinear optical properties.
Thermodynamic analysis
The computed parameters were summarized in Table 4. On the basis of vibrational analysis at B3LYP/ 6-31G (d, p), the standard statistical thermodynamic functions: heat capacity Cp,mo entropy Smo and enthalpy change Hmo for the title compound were calculated within the temperature range from 100 to 1000 K and as observed in Table S3. From Fig. S1, it can be observed that these thermodynamic functions are increasing with temperature due to the fact that the molecular vibrational intensities increase with temperature [49]. The correlation equations between the thermodynamic functions and temperatures were fitted by quadratic formulas, and R is the corresponding fitting factors (R2) for the thermodynamic properties. The corresponding R2 fitting factor of each thermodynamic quantities Cp,mo, Smo, and Hmo are 0.99956, 0.9998 and 0.99955, respectively. The corresponding fitting equations between thermodynamic properties and temperature are as follows:
Cp,mo = 57.73802+ 1.14384T -5.10116 ´ 10-4 T2
Smo = 295.48957+ 1.44264T -4.09224 ´ 10-4 T2
Hmo = -17.75073+0.19291T+2.8839 ´ 10-4 T2
All the thermodynamic data afford helpful information for the study of thermodynamic energies and estimate directions of chemical reactions according to the second law of thermodynamics in thermo chemical field [50]. Please note that all thermodynamic calculations were done in gas phase and they could not be used in solution.
Hirshfeld Surface Analyses
The molecular Hirshfeld surface [51,52]: dnorm, shape index, curvedness, de and di for the title compound are illustrated in Figs. S2(a-e) and mapped over: dnorm ranges (-0.632 to 1.120 Å), shape index ranges (- 1.000 to 1.000 Å), curvedness (-4.000 to 0.400 Å), de ranges (0.718 to 2.411 Å) and di ranges (0.719 to 2.430 Å), respectively. The examination of the full distribution of intermolecular interactions is shown in 2D map fingerprint plots (Fig. S3). The asymmetric property of the compound is due to the hydrogen bonds and the intermolecular interactions within the crystal. The 2D fingerprints (Fig. S3a) show that the intermolecular, O···H and H···H interactions are the dominant contacts. Fig. S3b represents the O···H/ H···O contacts between the hydrogen atoms located inside the Hirshfeld surface and the oxygen atoms characterized by two symmetrical points located at the top left and right. The largest contribution to the total Hirshfeld area is about (42.4%). Fig. S3c illustrates the interaction H···H. The contribution to the surface hirshfeld about (37.6%) in the surface grouped points which occupy half of surface. Fig. S3d shows the interaction between groups C···H/ H···C. These interactions have a contribution to the surface hirshfeld of (12.4%). The contribution of other intermolecular interactions is found to be in decreasing order: [C···C (3.7%), N···H (2.5%), N···C (1%), C···O (0.2%), N···O (0.1%) and H···As (0.1%)]. Fig. S4 summarized the distribution of individual intermolecular interactions on the basis of HS analysis for the title compound.
Vibrational studies
The experimental and theoretical FT-IR and FT-Raman spectra were shown in Figs. 6 and 7, respectively. Theoretically computed frequencies with their relative intensities, probable assignments and potential energy distribution (PED) are summarized in Table S4. To the assignments, only potential energy distributions (PED) contributions ≥ 10% and the experimental IR spectrum were considered. The theoretical intensities show slight deviations from the experimental values since the theoretical wavenumbers are obtained from the isolated molecule in the gaseous phase and the experimental wavenumbers are obtained from the isolated molecule in solid state.
Arsenate groups vibrations
According to the literature, the isolated arsenate anion [AsO4]3- with Td point group symmetry led to four Raman active normal modes: n1(A1), n2(E), n3(F2) and n4(F2) with average wavenumbers 837, 349, 887 and 463 cm-1, respectively [53]. n1 and n3 involve the symmetric and the asymmetric stretching mode of the As–O bond, whereas n2 and n4 involve mainly O–As–O symmetric and asymmetric bending modes. We have distinguished the four fundamental modes (ν1, ν2, ν3 and ν4) of the [AsO4]3- anion in the wave number region from 300 to 1000 cm-1. In the region, the [AsO4]3- stretching and bending vibrations expected to appear, as well as the modes associated with the internal modes of the organic cations. The band at 850 cm-1 in Raman spectrum is assigned to the asymmetric stretching 𝜈3 mode of the arsenate ion. In IR spectrum, this mode also appears as a single band at 865 cm-1 (with PED contribution of 93%). The band related to the symmetric stretching mode ν1 vibration is located at 817 cm-1 in the IR spectrum while in the Raman spectrum it is absent. The bands concerned with the asymmetric and symmetric bending vibration ν4 and ν2 of the arsenate ion are identified in the 550-700 cm-1 and 300-450 cm-1 frequency regions respectively. The very strong band observed in the Raman spectrum at 667 cm-1 and the strong band at 677 cm-1 in the IR spectrum can be easily assigned to 𝜈4 mode (with PED contribution of 88%). The same bending vibration is predicted at 665 and 677 cm-1 respectively, by DFT method. The symmetric mode 𝜈2 appears in the Raman spectrum as a shoulder band at 410 cm-1. This mode is theoretically calculated at 374 cm-1 by DFT method. As seen from Table S4, all modes associated to the arsenates groups are well reproduced by theoretical method.
C-H vibrations
Heterocyclic aromatic compounds and its derivatives are structurally very close to benzene. The C- H stretching frequency of such compounds falls very nearly in the region 3000 - 3100 cm-1 which is the characteristic region for the ready identification of C–H stretching vibrations. Hence, the both FT-IR and FT-Raman bands were observed at 3095 cm-1 and 3059 cm-1 assigned to CH asymmetric stretching. It is noticeable that asymmetric stretching vibrations occurred at higher wave number. The theoretically calculated harmonic wavenumber at 3230 cm-1 and 3094 cm-1 in B3LYP/6-31G(d,p) basis set respectively with PED contribution of 86%.
N-H vibrations
Generally, NH stretching modes cause to vibrational bands in absorption region above 3000-3500 cm-1. This band in the FT-IR spectrum is observed at 3432 cm-1, but the mode is not observed Raman spectra of the molecule. The theoretically calculated harmonic wavenumber at 3452 cm-1 and 3307 cm-1 in B3LYP/6-31G(d,p) basis set respectively with PED contribution of 97%. On the other hand, the NH2 in plane bending (scissoring) mode are observed as an individual band at 1592-1602 (IR-Raman). The theoretically calculated harmonic wavenumber at 1627 cm-1 and 1599 cm-1 in B3LYP/6-31G(d,p) basis set respectively with PED contribution of 89%.
C=C vibrations
The C=C skeletal stretching vibrations in aromatic rings generally appear at 1430-1625 cm-1
[54-55]. While the C=C stretching modes in the benzene rings of the molecule are observed at 1436 cm-1 in FT-IR and 1435 cm-1 in FT-Raman spectra. The wavenumber values calculated for these C=C vibrations are calculated at 1467 and 1442 cm-1 with 45% contributions of PED, respectively. Also it is observed at 1648 cm-1 in FT-IR spectrum. The theoretically calculated harmonic wavenumber at 1697 and 1668 cm-1 with B3LYP/6- 31G(d,p) basis sets in both FT-IR and FT-Raman spectra respectively.
C-N vibrations
The identification of C-N stretching modes in the side chains is a rather difficult task since there are problems in identifying these frequencies from other vibrations. The vibrational bands observed at 1336 and 1332 cm-1 in FT-IR and FT-Raman spectra are assigned to C-N vibrations. The theoretically calculated harmonic wavenumber at 1333 and 1310 cm-1 respectively with PED contribution of 56%. Additionally, the C-NH and C-NH2 stretching band are observed at 1040 and 1057 cm-1 in IR and Raman spectra. The wavenumber values calculated for these CN vibrations are calculated at 1074 and 1057 cm-1 with 67% contributions of PED, respectively.
O-H vibrations
At high wavenumbers, the IR and Raman spectra show an absorption band centered at 3042 and 3036 cm-1, respectively, which is assignable to the stretching vibration mode of water molecules (O–H). The theoretically computed values using DFT calculations are 3083 cm-1 and 2954 cm-1 (with PED contribution of 91%) are in good agreement with the experimental values.
Molecular docking simulation, druggability and pharmacokinetics
Molecular docking is a powerful computational tool for predicting the binding affinity of a ligand with proteins which is very useful in modern structure based drug designing [26,28]. The binding interactions of the ligand with the target protein or enzyme are generally through intermolecular bonds [27,29,30]. This method is crucial in drug design and it’s is frequently the preferred way for estimating the bioactive conformation of the ligand in the protein binding region [28,56,57]. Using iGEMDOCK program, we can determine the interaction energy of a ligand with a protein (also called free energy binding). Generally, this interaction energy is not found from a single structure. In fact, it corresponds to the average energy of several stable structures. The results reveal four docked poses, so four energy values. The best position that has the lowest energy corresponds to the most stable protein/ligand complex [58,59]. Docking calculations of interaction energies were subsequently clustered in Table 5. Fig. 8 represents the molecular docking of Alpha-synuclein protein (3q25) with the two ligands that are Aresnic and Maltose. Molecular docking led to the following results: the total energies equal to -93.6683 and -107.384 kcal/mol for Aresnic and Maltose interacting with Alpha-synuclein protein, respectively. Likewise, for Alpha-synuclein protein, the Maltose ligand is the strongest binding since it possesses the strongest energy (in absolute value) -107.384 kcal/mol, the van der Waals interaction -55.243 kcal/mol, the strongest averconpair 26.3043 kcal/mol and H-bond interaction -52.1406 kcal/mol. The ligand Aresnic is found with energy equal to -93.6683 kcal/mol. Following the previous observations, a new protein was studied called Rho-kinase (2f2u) in interaction with the following Ligand: Aresnic and 5-(1,4-diazepan-1-sulfonyl) isoquinoline (M77). Localization of the best docked poses of these complexes estimated via docking calculations is given in Fig. 9. Table 5 analysis shows that our compound 5-(1,4-diazepan-1-sulfonyl) isoquinoline takes the highest energy value which is equal to -95.4741 kcal/mol, also it presents the strongest van der Waals interaction -75.5128 kcal/mol) and averconpair 28.2 kcal/mol. According to the docking results, the interactions observed between Alpha-synuclein (and also between Rho-kinase) with the Aresnic ligand show the good efficacy of inhibition of the studied compound, indicating that they are promising inhibitors in the treatment of Alzheimer's diseases and in general way of psychic diseases. Later on, the antiviral potential of the compound was assessed. Four different proteins were targeted: 6LU7, 6M0J, 6M2N and 7BTF. The binding affinities of the complex, its interaction bonds, the closest interacting residues and the distance to closest interacting residues are shown in table 6. The results showed negative binding energies ranging from -3.1 to -4.1 kcal/mol and up to five conventional hydrogen bonds, particularly with 6LU7, which is the SARS-CoV-2 main protease [58,60]. The latter is highly involved in SARS-CoV-2 inducing COVID-19 pathogenesis and even its complication with acute respiratory distress syndrome, which might lead to death [60]. Regarding the measure of the distance to the closest interaction residue, it has been predicted that it ranges between 2.267 A to 2.714 A for 6M0J and 6LU7 respectively, which suggest that the compound is strongly embedded in the different protein pocket regions. These data prove that the compound filled well into the pocket region of the targeted proteins and established conventional hydrogen bonds, electrostatic and hydrophobic bonds, which satisfactory support the potential biological effect of the compound, particularly its antiviral effect as compared with previous published reports [26,59,60]. In fact, it has been reported that several biological and synthesized compounds possessed potential ameliorative effect against SARS-CoV-2 inducing COVID-19 [26,59,60,61]. To rationalize these results and the possible use of its novel organic arsenate templated by 4-aminopyridine compound, its lipophilicity, druggability, medicinal chemistry and pharmacokinetics were assessed. The findings are shown in table 7. The compound showed acceptable consensus Log Po/w and bioavailability score (0.97). While its absorption via both the gastro-intestinal tract and the skin were low (Log Kp = -10.51 cm/s), it was not categorized as blood-brain-barrier (BBB) permeant but it was a substrate of the p-glycoprotein. Cytochrome P450 (CYP) isoforms, such as CYP1A2, CYP2C19, CYP2C9, CYP2D6 and CYP3A4, are commonly used in biotransformation of drugs and xenobiotics [27,30,58]. In this study, the screened CYP enzyme isoforms data indicate that the compound is inhibiting none of these isoenzymes. Thus, suggesting promising pharmacokinetic properties.