Deciphering the Molecular Interaction Process of Gallium Maltolate on SARS-CoV-2 Main and Papain-Like Proteases: A Theoretical Study

: This study explored the inhibitory potential of gallium maltolate against severe acute respiratory syndrome coronavirus 2 and main and papain-like proteases. Computational methods, including density functional theory and molecular docking, were used to assess gallium maltolate reactivity and binding interactions. Density functional theory calculations revealed gallium maltolate’s high electron-capturing capacity, particularly around the gallium metal atom, which may contribute to their activity. Molecular docking demonstrated that gallium maltolate can form strong hydrogen bonds with key amino acid residues like glutamate-166 and cysteine-145, tightly binding to main and papain-like proteases. The binding energy and interactions of gallium maltolate were comparable to known SARS-CoV-2 inhibitors like N-[(5-methyl-1,2-oxazol-3-yl)carbonyl]-L-alanyl-L-valyl-N- {(2S,3E)-5-(benzyloxy)-5-oxo-1-[(3S)-2-oxopyrrolidin-3-yl]pent-3-en-2-yl}-L-leucinamide, indicating its potential as an antiviral agent. However, further experimental validation is required to confirm its effectiveness in inhibiting SARS-CoV-2 replication and treating COVID-19.


Introduction
Over the years, the pharmaceutical industry has been dominated mainly by organic chemistry; however, with the emergence of metallated drugs, a remarkable growth in inorganic medicinal chemistry has developed.In its beginnings, a well-known anticancer drug, cis-platinum (cis-[PtCl 2 (NH 3 ) 2 ]), a compound that forms adducts with DNA inhibiting its replication and transcription processes [1,2], leading to cell death by apoptosis, stood out [3].However, this drug presented side effects such as nausea and neurotoxicity, which motivated different researchers to develop derivatives of this drug in order to improve its effectiveness and reduce its side effects [4,5].Thus, over the years, different In this comprehensive research endeavor, we present an extensive in silico investigation focused on unraveling the inhibitory potential of gallium maltolate (GaM) against the SARS-CoV-2 main proteases, M pro and PL pro .This multifaceted investigation was undertaken to shed light on the therapeutic prospects of GaM in combating the COVID-19 pandemic.This computational study was initiated by employing advanced reactivity descriptors, which are meticulously defined within the purview of the density functional theory (DFT) [31].These descriptors were instrumental in scrutinizing and delineating the reactive attributes of the organometallic GaM molecule.It is worth noting that the use of DFT-based reactivity descriptors has been a proven and widely accepted methodology in the scientific community.It has been successfully utilized by numerous research groups, underscoring its efficacy in elucidating the reactivity patterns of various molecular systems.Following the rigorous assessment of GaM's reactivity, we employed molecular docking techniques to assess its potential as an inhibitor of SARS-CoV-2.The results of this computational analysis yielded highly promising outcomes regarding GaM's ability to inhibit the virus.Notably, the binding energies exhibited by GaM closely approximated those of N3, a well-established reference molecule known for its potent inhibitory activity against SARS-CoV-2.This suggests that GaM holds substantial promise as an antiviral agent.The basis for GaM's remarkable inhibitory activity lies in its ability to forge strong interactions with the key amino acid residues Glu166 and Cys145, which are pivotal for the activity of the M pro enzyme.These interactions are primarily facilitated by the formation of hydrogen bonds, underscoring the strategic significance of GaM in targeting the vital components of the virus's replication machinery.

Koopmans' Theorem References
Global hardness (η) Electron acceptor (ω Net electrophilicity (∆ω ± ) Additionally, three-dimensional molecular electrostatic potential (MEP) surfaces were generated to depict the electronic charge distributions of GaM, TTT, and N3 molecules.These MEP surfaces provide a qualitative picture of chemical reactivity and site-specific interactions, as they provide a visual representation of electron density variations.Electrophilic (electron-deficient) and nucleophilic (electron-rich) attack sites are readily identified, enhancing our understanding of the potential noncovalent interactions, particularly within the active sites of protein-ligand complexes.The color-coded MEPs function as a predictive framework for electrostatic interaction potential: regions with a positive poten-tial (shown in shades of blue) indicate likely targets for nucleophilic attacks, whereas areas with a negative potential (transitioning from green to red) are suggestive of electrophilic character, predisposing them to nucleophilic encounters.This nuanced visualization facilitates the discernment of key interactive sites that govern ligands' binding affinity and specificity to their protein counterparts.

Molecular Docking
Molecular docking simulations were carried out employing AutoDock 4.0 [47,48], a widely recognized computational tool in medicinal chemistry for predicting the preferred orientation of small molecules to their protein targets.The coordination compound of interest, gallium maltolate (GaM), was procured from the Cambridge Crystallography Data Centre (CCDC) [49].Benchmarking its interactions, we included two established co-crystallized ligands, specifically 5-amino-2-methyl-N-[(1R)-1-naphthalen-1-ylethyl] benzamide (TTT) and the complex N-[(5-Methylisoxazol-3-yl)Carbonyl]Alanyl-L-Valyl-N1-((1R,2Z)-4-(Benzyloxy)-4-Oxo-1-{[(3R)-2-Oxopyrrolidin-3-yl]methyl}But-2-Enyl)-L-Leucinamide (N3), both of which were retrieved from the Protein Data Bank (PDB) [50].To tailor the computational environment for the unique characteristics of GaM, particularly the gallium ion, modifications were incorporated into the parameter file of AutoDock 4.0 using Autodock Tools (ADT) [51].This step ensured accurate grid box delineation around the enzymatic binding site, as specified in the Grid Parameter File (GPF).Such precision in the computational setup is crucial for reflecting the physicochemical environment of the metal-ligand interaction within the active site.The ligands' geometries, including those of GaM, TTT, and N3, were optimized at the PBE0/Def2-TZVP level using the Gaussian16 program, with Gasteiger charges applied to approximate the molecular electrostatic potential, informing the likely intermolecular interactions.For the receptor proteins, SARS-CoV-2 main protease M pro (PDB code 6LU7) [52], and SARS-CoV-2 papain-like protease PL pro (PDB code 7CMD) [53], datasets were meticulously curated.A series of bioinformatics protocols were implemented to validate and refine the structural models using the Protein Preparation Wizard [54] from the Schrödinger suite, adopting the OPLS force field [55] to ensure computational modeling fidelity.The binding sites were inferred from the centroid of the co-crystallized ligands [56-58], which is a critical determinant for the generation of the grid maps by AutoGrid 4.0.These maps were constructed with a resolution of 0.375 Å across a 60 3 -point grid, establishing the spatial context for the docking calculations.The Lamarckian Genetic Algorithm (LGA) underpinned the search strategy, simulating an evolutionary process to optimize the ligand conformations.Initial conditions included a diverse population of 150 random individuals, with the algorithm permitted to iterate over 2.5 × 10 7 energy evaluations and 27,000 generations, subjected to a mutation rate of 0.02 and a crossover rate of 0.80, to converge on the energetically favorable binding poses.The end-state of this computational endeavor was the elucidation of the lowest-energy bound conformations for GaM in complex with M pro and PL pro , offering structural insights that are instrumental for the design and development of potential therapeutic agents against SARS-CoV-2.

Non-Covalent Interactions
To gain deeper insights into the intermolecular interactions, we probed the noncovalent binding between the gallium maltolate (GaM) ligand and the SARS-CoV-2 main protease (M pro ) as well as the papain-like protease (PL pro ).The analysis was centered on the energetically most favorable GaM-protein complexes, employing the Non-Covalent Interaction (NCI) index [59] computed with the NCIPLOT software [60].This methodological framework is based on the exploration of the electron density (ρ(r)) landscape and its gradient properties, principally through the reduced density gradient (RDG, represented as s).The RDG was computed through the following equation: The visualization of NCI isosurfaces, corresponding to areas of low RDG values and low electron density, allows for the spatial localization of non-covalent interactions.Additionally, the qualitative nature of these interactions is discerned by examining the sign of the second eigenvalue of the electron density Hessian matrix (λ 2 ) in conjunction with the electron density ρ(λ 2 ).The resulting NCI isosurfaces are color-coded: blue indicates regions of strong and attractive interactions, typically hydrogen bonds; green denotes weaker interactions, such as van der Waals forces; and red highlights areas with strong, repulsive steric interactions.This chromatic coding provides a direct and interpretable mapping of interaction strengths and types, enhancing understanding of the molecular interplay at the binding interface.

Reactivity Indexes and Molecular Electrostatic Potential Analysis
The therapeutic efficacy of a pharmacological agent is contingent upon its molecular affinity for binding to a biological target.This affinity can manifest as reversible or irreversible interactions based on the drug's inherent chemical properties [61].Reversible agents typically engage in non-covalent bonding with the target, encompassing a spectrum of weak intermolecular forces such as hydrogen bonds, ionic interactions, van der Waals interactions, and hydrophobic contacts.These interactions, although transient, are crucial for the reversible nature of the drug-target complex.In contrast, irreversible inhibitors form stable covalent bonds with specific amino acid residues within the target site, a process commonly facilitated by an electrophilic fragment within the molecule that reacts with a nucleophilic cysteine residue present in the enzyme's active site.
The propensity of a molecule to form such bonds, indicative of its reactivity or binding potential, can be theoretically appraised through frontier molecular orbital (FMO) analysis.This entails an evaluation of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) [62].The FMO theory suggests that a molecule's chemical reactivity is established by the orbital interactions wherein a nucleophilic HOMO of the pharmacophore can interact with an electrophilic LUMO of the biological target.The HOMO-LUMO gap serves as an index of reactivity, with a smaller gap indicative of heightened reactivity due to more accessible orbital interactions.Further, FMOs provide insights into the potential electron flow from the drug to the target, aiding in the elucidation of reactivity through various chemical descriptors [62].This analytical perspective is pivotal in guiding the structural optimization of both reversible and irreversible inhibitors.It enables the strategic positioning of noncovalent interaction points and the incorporation of electrophilic moieties designed to form covalent bonds with the biological target [63].Therefore, integrating noncovalent interaction analyses and covalent bonding potential via a molecular orbital approach could substantially refine the process of rational drug design within the domain of medicinal chemistry.
The results in Table 2 provide insights into the electronic structure and reactivity of the GaM, TTT, and N3 compounds in vacuum and water environments.Overall, the data indicate that solvation in water leads to the stabilization of all three compounds, with the HOMO being stabilized more than the LUMO.This increases the HOMO-LUMO gap, reflecting greater stability in the solvent.When these compounds are placed in a water environment, the HOMO energy becomes more negative and the LUMO energy becomes more positive compared to their vacuum counterparts.For N3, the HOMO energy ranges from −7.0 to −7.3 eV, indicating that it readily donates electrons, making it an effective electron donor.On the other hand, GaM exhibits the least negative HOMO energy, ranging from −5.7 to −6.2 eV, suggesting that it is less prone to donating electrons.Regarding LUMO energy, N3 exhibits the lowest LUMO energy, ranging from −1.5 to −1.6 eV, indicating its high electron-accepting ability.GaM, with its higher LUMO energy ranging from −1.2 to −1.5 eV, is less effective as an electron acceptor.The fact that the GAP increases in water for all compounds suggests that water solvation enhances the stability of these compounds.N3, with the largest GAP ranging from 5.4 to 5.8 eV, should be the most stable among the three, i.e., less prone to undergoing chemical reactions.In contrast, GaM and TTT with smaller GAPs (4.5 to 4.7 eV) should be comparatively more reactive.Electronegativity (χ) represents a molecule's ability to attract electrons, while global hardness (η) indicates its resistance to changes in electronic configuration.The slight increase in electronegativity and global hardness when the compounds are in a water environment suggests that they exhibit a stronger electron-attracting tendency and enhanced stability in water.N3, with the highest electronegativity (4.3 to 4.4 eV) and the highest global hardness range (2.7 to 2.9 eV), is the most electron-attracting and stable compound.In contrast, GaM, with the lowest electronegativity ranging from 3.5 to 3.9 eV and a global hardness of 2.3 eV, is less electron-attracting and less stable than the other compounds.The variations in electrophilicity (ω), electrodonating power (ω − ), and electroaccepting power (ω + ) in the presence of water are noteworthy.For N3, they slightly decrease in water, while GaM and TTT experience an increase in these descriptor values.N3 maintains the highest values for these parameters even after the slight decrease, signifying its strong reactivity and its ability to both donate and accept electrons.GaM exhibits the lowest ω − and ω + in both solvents, indicating its reduced propensity to donate or accept electrons.The increase in net electrophilicity (∆ω ± ) in the presence of water implies that water enhances the overall reactivity of these compounds.Once again, N3 stands out with the highest ∆ω ± , indicating that it is the most reactive compound among the three.GaM, with the lowest ∆ω ± , exhibits lower reactivity compared to the other compounds.The results show N3 is the most reactive compound with the highest electronegativity, hardness, and electrophilic power, while GaM is the least reactive with the lowest values.Solvation in water stabilizes all three compounds.
According to the observed FMO calculations, GaM exhibits a reactivity profile quite comparable to that of the TTT molecule, a well-known reference compound.However, GaM distinguishes itself with a notably heightened capacity for electron capture, as indicated by its net electrophilicity.This unique ability to capture electrons seems to have its epicenter on the metal atom within GaM, but it does not confine itself to just the metal atom; instead, it extends its influence to neighboring atoms by drawing in the surrounding electronic cloud.In order to visualize and precisely quantify the magnitude of this electron-attracting effect, we employed molecular electrostatic potential (MEP) calculations for the three compounds under investigation.These MEP calculations provide a detailed map of the electron distribution around the molecules, shedding light on the areas where electrons are either drawn toward or repelled from the molecular structure.The insights gained from these MEP calculations will help us better understand the electron-capturing behavior of GaM and its implications for potential applications.
The electrostatic molecular potential serves as a valuable and insightful analytical instrument for evaluating drug reactivity when exposed to electrophilic and nucleophilic attacks.This assessment is crucial for understanding how pharmaceutical compounds interact with their target receptors and biological systems [64].As depicted in Figure 1, the electrostatic molecular potential highlights specific charge distributions and regions of interest.Notably, GaM exhibits a distinctive charge distribution pattern.It is characterized by a pronounced positive charge (depicted in blue) concentrated around the methyl groups, indicating an electrophilic nature in these areas.Simultaneously, there is a negative charge predominantly localized on the metal-bound carboxylic groups, signifying susceptibility to nucleophilic attacks.This dual charge distribution suggests that GaM may engage in both electrophilic and nucleophilic interactions, making it a versatile candidate for pharmaceutical reactivity.In contrast, another compound, TTT, presents a distinct electrostatic profile.In Figure 1, TTT shows positive charge regions that surround its aromatic ring, depicted in green.These positive regions indicate a hydrophobic nature, which might suggest a preference for interactions in non-polar environments.Understanding these electrostatic characteristics is vital for predicting how TTT may interact with hydrophobic binding sites within biological systems or receptor proteins.
The electrostatic molecular potential serves as a valuable and insightful analytical instrument for evaluating drug reactivity when exposed to electrophilic and nucleophilic attacks.This assessment is crucial for understanding how pharmaceutical compounds interact with their target receptors and biological systems [64].As depicted in Figure 1, the electrostatic molecular potential highlights specific charge distributions and regions of interest.Notably, GaM exhibits a distinctive charge distribution pattern.It is characterized by a pronounced positive charge (depicted in blue) concentrated around the methyl groups, indicating an electrophilic nature in these areas.Simultaneously, there is a negative charge predominantly localized on the metal-bound carboxylic groups, signifying susceptibility to nucleophilic attacks.This dual charge distribution suggests that GaM may engage in both electrophilic and nucleophilic interactions, making it a versatile candidate for pharmaceutical reactivity.In contrast, another compound, TTT, presents a distinct electrostatic profile.In Figure 1, TTT shows positive charge regions that surround its aromatic ring, depicted in green.These positive regions indicate a hydrophobic nature, which might suggest a preference for interactions in non-polar environments.Understanding these electrostatic characteristics is vital for predicting how TTT may interact with hydrophobic binding sites within biological systems or receptor proteins.
Figure 1.A schematic representation of molecular electrostatic potential (MEP) maps (isovalue is 0.03 a.u.) for GaM, TTT, and N3 ais shown in both a vacuum and a water environment.In these representations, areas shaded in red, orange, and yellow denote electrophilic reactivity, while those in green and blue correspond to nucleophilic reactivity.The following elements have been highlighted in this figure: carbon (gray), hydrogen (light gray), oxygen (red), gallium (purple), and nitrogen (blue).

Molecular Docking Analysis
The structural features of two crucial enzymes associated with SARS-CoV-2 are the main protease (M pro ) and the papain-like protease (PL pro ).These enzymes are pivotal in mediating viral replication and transcription [52,65], making them prime targets for drug development and therapeutic interventions.M pro , also known as 3CLpro, is a dimeric enzyme with each protomer consisting of three distinct domains.The first two domains, Domain I and Domain II, form an antiparallel β-barrel structure, while Domain III is a Figure 1.A schematic representation of molecular electrostatic potential (MEP) maps (isovalue is 0.03 a.u.) for GaM, TTT, and N3 ais shown in both a vacuum and a water environment.In these representations, areas shaded in red, orange, and yellow denote electrophilic reactivity, while those in green and blue correspond to nucleophilic reactivity.The following elements have been highlighted in this figure: carbon (gray), hydrogen (light gray), oxygen (red), gallium (sapphire), and nitrogen (blue).

Molecular Docking Analysis
The structural features of two crucial enzymes associated with SARS-CoV-2 are the main protease (M pro ) and the papain-like protease (PL pro ).These enzymes are pivotal in mediating viral replication and transcription [52,65], making them prime targets for drug development and therapeutic interventions.M pro , also known as 3CLpro, is a dimeric enzyme with each protomer consisting of three distinct domains.The first two domains, Domain I and Domain II, form an antiparallel β-barrel structure, while Domain III is a globular cluster of five α-helices connected to Domain II by a long loop.The substrate binding site, where the catalytic action occurs, is between Domain I and Domain II and involves a catalytic dyad of Cys145 and His41 [52,65].The substrate binding pocket of M pro is highly conserved among coronaviruses, making it an attractive target for the development of broad-spectrum antiviral drugs.Inhibitors such as N3 and N1 covalently bind to Cys145, effectively occupying this pocket.The binding of inhibitors can induce significant conformational changes in M pro , effectively sealing the entrance to the substrate binding pocket.This structural flexibility is crucial for accommodating different substrates and inhibitors [52,65,66].PL pro is a monomeric enzyme characterized by four subdomains: the N-terminal ubiquitin-like domain (Ubl), the α-helical thumb domain, the β-stranded finger domain, and the palm domain.The catalytic triad responsible for PL pro 's enzymatic activity includes Cys111, His272, Tyr264, Tyr268, and Asp286 [53,56,67,68].PL pro is highly similar to its counterpart in the SARS-CoV virus, sharing conserved structural features and an impressive 82% sequence identity.This similarity suggests that research on SARS-CoV PL pro can inform studies on SARS-CoV-2 PL pro .Crystal structures have shown that the binding of inhibitors causes the closure of the BL2 loop, sealing the cleft entrance where substrate cleavage occurs.The S3-S4 subsites within the substrate binding cleft are critical for accommodating key inhibitor groups, such as the zinc-binding moiety of TTT.Inhibitors targeting cysteine residues in the zinc-binding domain are particularly promising for drug development [53,56,67].
To investigate the possible mechanism by which selected GaM acts, molecular docking analysis of GaM was carried out at the active sites of SARS-CoV-2 M pro and PL pro .The molecular docking scores are displayed in Table 3.The results of molecular docking studies are presented in terms of interaction energy (∆E binding ), which reflects the strength of binding between the ligands and the target proteins.The ligands GaM, TTT, and N3 were evaluated for their binding interactions with two key proteins of the SARS-CoV-2 virus, M pro and PL pro .GaM demonstrated moderate binding affinities with both M pro and PL pro .With ∆E binding values of −6.98 kcal/mol for M pro and −5.90 kcal/mol for PL pro , GaM appears to form reasonably favorable interactions with both proteins.TTT exhibited a significant interaction with SARS-CoV-2 PL pro , with an impressive ∆E binding value of −10.34 kcal/mol.This indicates a strong and favorable binding affinity.On the other hand, TTT has a higher ∆E binding value specified for its interaction with M pro , which indicates that the calculations performed are correct since it is not its original protein.N3, showed an ∆E binding of −7.29 kcal/mol with M pro , similar to the GaM interaction with M pro .However, N3 has a higher ∆E binding value specified for its interaction with PL pro , which indicates that the calculations performed are correct since it is not its original protein.Table 4 presents a detailed overview of the interactions between GaM organometallic molecules and the amino acid residues of SARS-CoV-2 M pro and PL pro proteases.This interaction data sheds light on the molecular mechanisms through which GaM interacts with these proteases, providing valuable insights into potential inhibitory effects.In the case of M pro , several amino acid residues are involved in interacting with GaM.Notably, Cys145, Glu166, and Gln189 are identified as forming hydrogen bonds with GaM.These hydrogen bonds imply a strong and specific binding between GaM and M pro , indicating a potentially robust inhibitory effect.Additionally, other amino acids like Phe140, His163, Ser144, and Ser305 participate in the interaction primarily through van der Waals forces and hydrophobic interactions.This diverse array of interactions suggests that GaM can engage with M pro through multiple mechanisms, further reinforcing the likelihood of effective inhibition.In contrast, GaM interacts with PL pro through a different set of amino acid residues.Only two amino acids, Arg166 and Tyr264, form hydrogen bonds with GaM in the case of PL pro .This limited number of hydrogen bonding interactions in comparison to M pro implies a comparatively weaker binding affinity for GaM with PL pro .Other amino acids in PL pro contribute to the interaction mainly through van der Waals forces and hydrophobic interactions, which is similar to the interaction pattern observed with M pro .GaM forms a greater number of hydrogen bonds with M pro as compared to PL pro .This difference in hydrogen bonding interactions suggests that GaM likely exhibits a stronger binding affinity for M pro , making it a more potent inhibitor for this protease.The distinct sets of interacting residues also indicate that GaM binds to different regions or sites on the two proteases.These insights underscore the potential of GaM as an inhibitor of these proteases and suggest that it might exert a stronger inhibitory effect on M pro due to a greater number of hydrogen bonding interactions.This knowledge can be instrumental in the development of therapeutic strategies targeting SARS-CoV-2 proteases for antiviral drug development.All these interactions are complemented by the analyses performed using the noncovalent interaction index (Figure 2).This index provides us with a deeper and more accurate view of the relationships that are established between the components (in terms of their non-covalent interactions).These relationships can include hydrogen bonds, dispersion forces, ionic interactions, and other types of non-covalent attractions.The results of these analyses allow for a more complete understanding of how the molecules or particles involved are intertwined and affect each other, which in turn sheds light on the properties and behavior of the system under study.Attractive interactions appear at λ 2 < 0, whereas in cases where λ 2 is positive (as in rings or cages), usually several atoms interact but are not bonded, which corresponds to steric crowding according to classical chemistry.The blue regions indicate pronounced attractive interactions, whereas the green regions indicate relatively weaker interactions.

Conclusions
In conclusion, the computational studies conducted on gallium maltolate (GaM) demonstrate its promising potential as an antiviral agent against SARS-CoV-2.Density functional theory (DFT) calculations have revealed GaM's distinct reactivity, showcasing its increased electron-accepting ability in comparison to reference inhibitors like N3.Moreover, molecular docking studies have indicated GaM's effective binding to the active sites of SARS-CoV-2 proteases M pro and PL pro , with binding energies comparable to established potent inhibitors such as N3.A detailed analysis has uncovered that GaM forms important bonding interactions, particularly hydrogen bonds, with crucial amino acid residues like Glu166 and Cys145 in M pro .These interactions likely play a pivotal role in facilitating GaM's ability to inhibit viral replication by targeting essential components of the viral machinery.Consequently, based on the wealth of computational evidence, GaM is supported as a promising candidate for further investigation as an antiviral therapeutic against COVID-19 or as a foundational basis for delving into the broader realm of metallodrugs, particularly those involving gallium.

Biophysica 2024, 4 , 10 Figure 2 .
Figure 2. Schematic representations of the NCIPLOT isosurface gradient (0.5 a.u.) of GaM on the structure of SARS-CoV-2 (A,B) M pro and (C,D) PL pro .The sign of λ2 enables the identification of the interaction type.Attractive interactions appear at λ2 < 0, whereas in cases where λ2 is positive (as in rings or cages), usually several atoms interact but are not bonded, which corresponds to steric crowding according to classical chemistry.The blue regions indicate pronounced attractive interactions, whereas the green regions indicate relatively weaker interactions.

Figure
Figure Schematic representations of the NCIPLOT isosurface gradient (0.5 a.u.) of GaM on the structure of SARS-CoV-2 (A,B) M pro and (C,D) PL pro .The sign of λ 2 enables the identification of the

Table 1 .
Summary of equations used to calculate various global reactivity indexes in the TAFF[36]

Table 2 .
Global reactivity descriptors for GaM, N3, and TTT molecules optimized both in vacuum and water solvents at the PBE0/Def2-TZVP level.All units are measured in eV.

Table 3 .
Molecular docking study between selected ligands, GaM, TTT, and N3, and SARS-CoV-2 M pro and PL pro .Intermolecular docking values (*), presented with their interaction energy (∆E binding ).All units are measured in kcal/mol.Values are listed as a three-colored scheme from green (high energy) to red (low energy). *

Table 4 .
Amino acid residues of SARS-CoV-2 M pro and PL pro and hydrogen bonding (*) with the GaM molecule within a distance of 3.5 Å.