Pnictogen bond-driven control of the molecular interaction between organophosphorus and acetylcholinesterase enzyme.

This study addresses a comprehensive assessment of the interaction between chemical warfare agents (CWA) and acetylcholinesterase (AChE) systems, focus on the intriguing pnictogen-bond interaction (PnB). Utilizing the crystallographic data from the Protein Data Bank pertaining to the AChE-CWA complex involving Sarin (GB), Cyclosarin (GF), 2-[fluoro(methyl)phosphoryl]oxy-1,1-dimethylcyclopentane (GP) and venomous agent X (VX) agents, the CWA is systematically displaced by increments of 0.1 Å along the PO bond axis, extending its distance by 4 Å from the original position. The AIM analysis was carried out and consistently revealed the presence of a significant interaction along the PO bond. Investigating the intrinsic nature of the PnB, the NBO and the EDA analysis unearthed the contribution of orbital factors to the overall energy of the system. Strikingly, this observation challenges the conventional σ-hole explanation commonly associated with such interactions. This finding adds a layer of complexity to understanding of PnB, encouraging further exploration into the underlying mechanisms governing these intriguing chemical phenomena.

Enhanced Sampling in Molecular Dynamics Simulations: How Many MD Snapshots can be Needed to Reproduce the Biological Behavior?

Since its early days in the 19th century, medicinal chemistry has concentrated its efforts on the treatment of diseases, using tools from areas such as chemistry, pharmacology, and molecular biology. The understanding of biological mechanisms and signaling pathways is crucial information for the development of potential agents for the treatment of diseases mainly because they are such complex processes. Given the limitations that the experimental approach presents, computational chemistry is a valuable alternative for the study of these systems and their behavior. Thus, classical molecular dynamics, based on Newton's laws, is considered a technique of great accuracy, when appropriated force fields are used, and provides satisfactory contributions to the scientific community. However, as many configurations are generated in a large MD simulation, methods such as Statistical Inefficiency and Optimal Wavelet Signal Compression Algorithm are great tools that can reduce the number of subsequent QM calculations. Accordingly, this review aims to briefly discuss the importance and relevance of medicinal chemistry allied to computational chemistry as well as to present a case study where, through a molecular dynamics simulation of AMPK protein (50 ns) and explicit solvent (TIP3P model), a minimum number of snapshots necessary to describe the oscillation profile of the protein behavior was proposed. For this purpose, the RMSD calculation, together with the sophisticated OWSCA method was used to propose the minimum number of snapshots.

Combining computational tools and experimental studies towards endocrine disruptors mitigation: A review of biocatalytic and adsorptive processes.

The endocrine disruptors (EDCs) are an important group of emerging contaminants, and their mitigation has been a huge challenge due to their chemistry complexity and variety of these compounds. The traditional treatments are inefficient to completely remove EDCs, and adsorptive processes are the major alternative investigated on their removal. Also, the use of EDCs degrading enzymes has been encouraged due to ecofriendly approach of biocatalytic processes. This paper highlights the occurrence, classification, and toxicity of EDCs with special focus in the use of enzyme-based and adsorptive technologies in the elimination of EDCs from ambiental matrices. Numerous prior reviews have focused on the discussions toward these technologies. However, the literature lacks theoretical discussions about important aspects of these methods such as the mechanisms of EDCs adsorption on the adsorbent surface or the interactions between degrading enzymes - EDCs. In this sense, theoretical calculations combined to experimental studies may help in the development of more efficient technologies to EDCs mitigation. In this review, we point out how computational tools such as molecular docking and molecular dynamics have to contribute to the design of new adsorbents and efficient catalytic processes towards endocrine disruptors mitigation.

A theoretical protocol for the rational design of the bioinspired multifunctional hybrid material MIP@cercosporin.

CONTEXT: Rational design of polymeric materials prepared with the molecular imprinting technology is gaining even more space, as it can provide the optimal conditions to direct the laboratory molecularly imprinting polymer (MIP) preparation, maximizing their efficiency while reducing costs and preparation time, when compared to try-and-error approaches. We perform a rational design of an MIP with specific cavities for cercosporin accommodation by means of computational tools. The main steps of an MIP preparation were simulated and it was found that the most appropriated functional monomer to be used in the MIP preparation for cercosporin is the acrylamide, while the most suitable crosslinking agent is found to be p-divinylbenzene. Also, the most suitable solvents to remove cercosporin from the cavity are those with low dielectric constant, such as chloroform. This kind of solvent can then be used in washing step, in the case of use the MIP for sensing destinations. On the other hand, solvents like water, which has high dielectric constants, can efficiently improve the interactions between cercosporin and the functional monomer acrylamide, being indicated when the objective is to attract or maintain the cercosporin inside the MIP cavity. Thus, a MIP@cercosporin hybrid material can be used in aqueous solutions more reliably, or even the cercosporin detection in this media can be favoured. In the selectivity analysis of the material prepared in this specific condition, the results point that this MIP can also detect elsinochrome A with high efficiency, and could be more selective for hypericin, altertoxin, hypocrelin A, and phleichrome mycotoxins. METHOD: The main steps of a MIP synthesis were theoretically simulated trough density functional theory (DFT) calculations aiming to direct and optimize the synthesis and applications of the material before the bench tests. Initially, in order to choose the most suitable functional to be employed for cercosporin calculations, eight of the DFT functionals that had been previously used for cercosporin calculations in literature were tested, which were the LCWPBE, B3LYP, CAM-B3LYP, M062-X, mPW1PW91, PBE0, TPSSh, and ωb97Xd. The theoretical 1H NMR chemical shifts for cercosporin molecule were calculated and compared with experimental results to analyze the performance of the functionals. Of all these, the best results were obtained with the TPSSh functional, employing the 6-31G(d,p) basis set, and this level of theory was then used for all the following steps. All the simulations were performed by means of geometry optimizations and frequency calculations. Additionally, AIM calculations were employed for further analysis of the interactions between the chosen functional monomer and cercosporin template in step 1, which was functional monomer selection. In washing step, the calculations were done using implicit solvation model, and finally, in selectivity tests, the putative "solid" MIP was simulated by freezing the positions of the monomers after the template remotion, and then other structurally similar toxins were placed in its cavity for the geometry optimizations and frequency calculations.

Vanadium complex as a potential modulator of the autophagic mechanism through proteins PI3K and ULK1: development, validation and biological implications of a specific force field for [VO(bpy)2Cl].

The modulation of autophagy has been presented as a very useful strategy in anticancer treatments. In this sense, the vanadium complex (VC) bis(2,2'-bipyridine)chlorooxovanadium(IV), [VO(bpy)2Cl], is known for its ability to induce autophagy in triple-negative breast cancer cells (TNBC). An excellent resource to investigate the role of VC in the induction of autophagy is to make use of Molecular Dynamics (MD) simulations. However, until now, the scarcity of force field parameters for the VC prevented a reliable analysis. The autophagy signaling pathway starts with the PI3K protein and ends with ULK1. Therefore, in the first stage of this work, we developed a new AMBER force field for the VC (VCFF) from a quantum structure, obtained by DFT calculations. In the second stage, the VCFF was validated through structural analyses. From this, it was possible to investigate, through docking and MD (200 ns), the performance of the PI3K-VC and ULK1-VC systems (third stage). The analyses of this last stage involved RMSD, hydrogen bonds, RMSF and two pathways for the modulation of autophagy. In general, this work fills in the absence of force field parameters (FF) for VC by proposing an efficient and new FF, in addition to investigating, at the molecular level, how VC is able to induce autophagy in TNBC cells. This study encourages new parameterizations of metallic complexes and contributes to the understanding of the duality of autophagic processes.Communicated by Ramaswamy H. Sarma.

Silver(I) and Copper(II) 1,10-Phenanthroline-5,6-dione Complexes as Promising Antivirulence Strategy against Leishmania: Focus on Gp63 (Leishmanolysin).

Leishmaniasis, caused by protozoa of the genus Leishmania, encompasses a group of neglected diseases with diverse clinical and epidemiological manifestations that can be fatal if not adequately and promptly managed/treated. The current chemotherapy options for this disease are expensive, require invasive administration and often lead to severe side effects. In this regard, our research group has previously reported the potent anti-Leishmania activity of two coordination compounds (complexes) derived from 1,10-phenanthroline-5,6-dione (phendione): [Cu(phendione)3].(ClO4)2.4H2O and [Ag(phendione)2].ClO4. The present study aimed to evaluate the effects of these complexes on leishmanolysin (gp63), a virulence factor produced by all Leishmania species that plays multiple functions and is recognized as a potential target for antiparasitic drugs. The results showed that both Ag-phendione (-74.82 kcal/mol) and Cu-phendione (-68.16 kcal/mol) were capable of interacting with the amino acids comprising the active site of the gp63 protein, exhibiting more favorable interaction energies compared to phendione alone (-39.75 kcal/mol) or 1,10-phenanthroline (-45.83 kcal/mol; a classical gp63 inhibitor) as judged by molecular docking assay. The analysis of kinetic parameters using the fluorogenic substrate Z-Phe-Arg-AMC indicated Vmax and apparent Km values of 0.064 µM/s and 14.18 µM, respectively, for the released gp63. The effects of both complexes on gp63 proteolytic activity were consistent with the in silico assay, where Ag-phendione exhibited the highest gp63 inhibition capacity against gp63, with an IC50 value of 2.16 µM and the lowest inhibitory constant value (Ki = 5.13 µM), followed by Cu-phendione (IC50 = 163 µM and Ki = 27.05 µM). Notably, pretreatment of live L. amazonensis promastigotes with the complexes resulted in a significant reduction in the expression of gp63 protein, including the isoforms located on the parasite cell surface. Both complexes markedly decreased the in vitro association indexes between L. amazonensis promastigotes and THP-1 human macrophages; however, this effect was reversed by the addition of soluble gp63 molecules to the interaction medium. Collectively, our findings highlight the potential use of these potent complexes in antivirulence therapy against Leishmania, offering new insights for the development of effective treatments for leishmaniasis.

Parameterization and validation of a new AMBER force field for an oxovanadium (IV) complex with therapeutic potential implications in Alzheimer's disease.

The scarcity of efficient force fields to describe metal complexes may be a problem for new advances in medicinal chemistry. Thus, the development of force fields for these compounds can be valuable for the scientific community, especially when it comes to molecules that show interesting outputs regarding potential treating of diseases. Vanadium complexes, for instance, have shown promising results towards therapeutics of Alzheimer's Disease, most notably the bis(maltolato)oxovanadium (IV). Therefore, the mainly goal of this work is to develop and validate a new set of parameters for this vanadium complex from a minimum energy structure, obtained by DFT calculations, where great results of the new force field are found when confronted with experimental and quantum reference values. Moreover, the new force field showed to be quite effective to describe the molecule of under study whilst GAFF could not describe it effectively. In addition, a case study points out hydrogen bonds in the vanadium complex-PTP1B system.

Evaluation of autophagy inhibition to combat cancer: (vanadium complex)-protein interactions, parameterization, and validation of a new force field.

Autophagy has drawn attention from the scientific community, mainly because of its significant advantages over chemotherapeutic processes. One of these advantages is its direct action on cancer cells, avoiding possible side effects, unlike chemotherapy, which reaches tumor cells and affects healthy cells in the body, leading to a great loss in the quality of life of patients. In this way, it is known that vanadium complex (VC) [VO(oda)(phen)] has proven inhibition effect on autophagy process in pancreatic cancer cells. Keeping that in mind, molecular dynamics (MD) simulations can be considered excellent strategies to investigate the interaction of metal complexes and their biological targets. However, simulations of this type are strongly dependent on the appropriate choice of force field (FF). Therefore, this work proposes the development of AMBER FF parameters for VC, having a minimum energy structure as a starting point, obtained through DFT calculations with B3LYP/def2-TZVP level of theory plus ECP for the vanadium atom. An MD simulation in vacuum was performed to validate the developed FF. From the structural analyses, satisfying values of VC bond lengths and angles were obtained, where a good agreement with the experimental data and the quantum reference was found. The RMSD analysis showed an average of only 0.3%. Finally, we performed docking and MD (120 ns) simulations with explicit solvent between VC and PI3K. Overall, our findings encourage new parameterizations of metal complexes with significant biological applications, as well as allow to contribute to the elucidation of the complex process of autophagy.

Molecular Dynamics-Assisted Interaction of Vanadium Complex-AMPK: From Force Field Development to Biological Application for Alzheimer's Treatment.

A large part of the world's population is affected by Alzheimer's disease (AD) and diabetes mellitus type 2, which cause both social and economic impacts. These two conditions are associated with one protein, AMPK. Studies have shown that vanadium complexes, such as bis(N',N'-dimethylbiguanidato)-oxovanadium(IV), VO(metf)2·H2O, are potential agents against AD. A crucial step in drug design studies is obtaining information about the structure and interaction of these complexes with the biological targets involved in the process through molecular dynamics (MD) simulations. However, MD simulations depend on the choice of a good force field that could present reliable results. Moreover, general force fields are not efficient for describing the properties of metal complexes, and a VO(metf)2·H2O-specific force field does not yet exist; thus, the proper development of a parameter set is necessary. Furthermore, this investigation is essential and relevant given the importance for both the scientific community and the population that is affected by this neurodegenerative disease. Therefore, the present work aims to develop and validate the AMBER force field parameters for VO(metf)2·H2O since the literature lacks such information on metal complexes and investigate through classical molecular dynamics the interactions made by the complex with the protein. The proposed force field proved to be effective for describing the vanadium complex (VC), supported by different analyses and validations. Moreover, it had a great performance when compared to the general AMBER force field. Beyond that, MD findings provided an in-depth perspective of vanadium complex-protein interactions that should be taken into consideration in future studies.

Exploring electronic, structural and dynamics parameters of phenylbenzothiazole complexes with Mn2+, Cu2+ and Zn2+ for designing new magnetic resonance imaging (MRI) probes: congruence between computation and spectroscopic data.

Cancer is one of the leading causes of human death worldwide, being one of the most serious problems faced by mankind. For the diagnosis, Magnetic Resonance Imaging (MRI), through effective contrast agents (Cas), has greatly helped in the diagnosis at the initial stages. However, it is necessary to include new compounds more effective and selective for cancer diagnosis. The complexes with Mn2+, Cu2+ and Zn2+ have received great attention due to their applications as CAs for MRI. Those materials can shorten the T2 and T2* transverse relaxation times. Thus, the representative structures for hyperfine coupling constants (HFCCs) were selected from docking results by frequency of occupancy calculations. From the Multivariate Analysis to obtain the PCA graphs in the choice of a representative conformations. it is possible to notice that the variable energy does not present a high correlation with the other variables, and structural factors, such as the spatial positions of the metal atoms, seem to be important in the reactivity of the complexes. Structural factors, such as the spatial positions of the metal atoms, seem to be important in the reactivity of the complexes. Theoretical findings suggest that the compounds are capable of increasing the Aiso values of the water molecules, but the complex [Zn(H2O)(NNO)] shows a greater influence, being more sensitive to the Electron paramagnetic resonance parameters than the complexes [CuCl(H2O)NNO] and [MnCl2(H2O)(NNO)] with the explicit solvent and the enzyme. MRI contrast agents have generated various problems due to their high toxicity. In this perspective, this compound may be a promising alternative for transporting the CAs into diseased tissue.Communicated by Ramaswamy H. Sarma.

Synergistic effect of intra- and intermolecular hydrogen bond in 2-haloethanols probed by infrared.

Fluorine is often considered the only halogen to effectively engage in hydrogen bonds, while the other halogens, particularly iodine, are not electronegative enough to participate as hydrogen bond acceptors in electrostatic interactions. 2-Fluoroethanol and 2-iodoethanol have been studied herein to test this assumption, since a highly stable gauche conformation can experience the intramolecular hydrogen bond. However, the infrared O  H stretching frequency indicates that the hydroxyl group in 2-fluoroethanol is not engaged in intramolecular hydrogen bond, while the corresponding vibration mode for 2-iodoethanol suggests that not only the O  H is engaged in such interaction, but also that intramolecular hydrogen bond may drive the conformational equilibrium in this molecule. Theoretical calculations support the covalent nature of this interaction, and provide evidence that intermolecular hydrogen bond with a water molecule, and probably with the polar solvents tested experimentally, occurs with the hydroxyl rather than with the iodine substituent, as conventionally, in order to keep the intramolecular hydrogen bond effective.

Intermolecular Covalent Interactions: Nature and Directionality.

Quantum chemical methods were employed to analyze the nature and the origin of the directionality of pnictogen (PnB), chalcogen (ChB), and halogen bonds (XB) in archetypal Fm Z⋅⋅⋅F- complexes (Z=Pn, Ch, X), using relativistic density functional theory (DFT) at ZORA-M06/QZ4P. Quantitative Kohn-Sham MO and energy decomposition analyses (EDA) show that all these intermolecular interactions have in common that covalence, that is, HOMO-LUMO interactions, provide a crucial contribution to the bond energy, besides electrostatic attraction. Strikingly, all these bonds are directional (i.e., F-Z⋅⋅⋅F- is approximately linear) despite, and not because of, the electrostatic interactions which, in fact, favor bending. This constitutes a breakdown of the σ-hole model. It was shown how the σ-hole model fails by neglecting both, the essential physics behind the electrostatic interaction and that behind the directionality of electron-rich intermolecular interactions. Our findings are general and extend to the neutral, weaker ClI⋅⋅⋅NH3 , HClTe⋅⋅⋅NH3 , and H2 ClSb⋅⋅⋅NH3 complexes.

Assessing the Therapeutic and Toxicological Profile of Novel Acetylcholinesterase Reactivators: Value of In Silico And In Vitro Data.

Organophosphorus compounds (OP) make up an important class of inhibitors, mostly employed as pesticides, even as chemical weapons. These toxic substances act through the inhibition of the acetylcholinesterase (AChE) enzyme, which results in elevated synaptic acetylcholine (ACh) levels, leading to serious adverse effects under the cholinergic syndrome. Many reactivators have been developed to combat the toxic effects of these AChE inhibitors. In this line, the oximes highlight because of their good reactivating power of cholinesterase enzymes. To date, no universal antidotes can reactivate AChE inhibited by any OP agent. This review summarizes the intoxication process by neurotoxic OP agents, along with the development of reactivators capable of reversing their effects, approaching aspects like the therapeutic and toxicological profile of these antidotes. Computational methods and conscious in vitro studies, capable of significantly predicting the toxicological profile of these drug candidates, might support the process of development of these reactivators before entering in vivo studies in animals, and then clinical trials. These approaches can assist in the design of safer and more effective molecules, reducing related cost and time for the process.

Exploring 129Xe NMR parameters for structural investigation of biomolecules: relativistic, solvent, and thermal effects.

In recent years, the study of new probes has aroused great interest in the scientific community around the world. Therefore, in the present work, we present a potential candidate for a new spectroscopic probe, the Xe(CO)3(NNO) conjugated to 2-(4'-aminophenyl) benzothiazole complex, XeABT. For this proposal, chemical shift calculations at the DFT level were performed; thus, a factorial design was carried out in order to choose the best computational method. The best combination was the base function ZORA-def2-TZVP, with the functional PBE0 and considering the relativistic effects with the ZORA implementation. Our findings reveal that the 129Xe chemical shifts are affected by thermal and solvent effects, and considering an enzymatic environment, a significant decrease in δ(129Xe) values is observed, suggesting with the XeABT complex it may be a promising spectroscopic probe.

Improved Protocol for the Selection of Structures from Molecular Dynamics of Organic Systems in Solution: The Value of Investigating Different Wavelet Families.

Wavelets are mathematical tools used to decompose and represent another function described in the time domain, allowing the study of each component of the original function with a scale-compatible resolution. Thus, these transforms have been used to select conformations from molecular dynamics (MD) trajectories in systems of fundamental and technological interest. Recently, our research group has used wavelets to develop and validate a method, meant to select structures from MD trajectories, which we named OWSCA (optimal wavelet signal compression algorithm). Here, we moved forward on this project by demonstrating the efficacy of this method on the study of three different systems (non-flexible organic, flexible organic, and protein). For each system, 93 wavelets were investigated to verify which is the best one for a given organic system. The results show that the best wavelets were different for each system and, also, very close to the experimental values, with the wavelets db1, rbio 3.1, and bior1.1 being selected for the non-flexible, flexible organic, and protein systems, respectively. This reinforces our OWSCA as a very efficient and promising method for the selection of structures from MD trajectories of different classes of compounds. Our findings also point out that additional studies considering wavelet families are needed for defining the best wavelet for representing each system under study.

Insights into the value of statistical models, solvent, and relativistic effects for investigating Re complexes of 2-(4'-aminophenyl)benzothiazole: a potential spectroscopic probe.

Cancer affects a major part of the worldwide population, and, to minimize deaths, the diagnosis in the early stages of the disease is fundamental. Thus, to improve diagnosis and treatment new potential spectroscopic probes are crucial. Benzothiazole derivates present antitumor properties and are highly selective and interact strongly with the enzyme phosphoinositide 3-kinase (PI3K), which was associated with cell proliferation and breast cancer cells. In this paper, the rhenium shielding tensors (187Re(σ)) and hydrogen and carbon chemical shifts (1H(δ) and 13C(δ)) of the Re(CO)3(NNO) complex conjugated with 2-(4'-aminophenyl)benzothiazole (ReABT) were evaluated. A statistical HCA model was used to analyze the best DFT protocol to compute σ and δ values and to evaluate the relativistic effects, both in the basis set and Hamiltonian as well as the functionals M06L or PBE0. The best protocol was applied to obtain 187Re(σ) of the ReABT complex in different environments (gas phase, solution, and in the active site of the PI3K enzyme). The results point out that 187Re(σ) values of the ReABT complex change significantly when the complex is docked in the PI3K enzyme.

Value of Contrast-enhanced Magnetic Resonance Imaging (MRI) in the Diagnosis of Breast Cancer.

This review article aims to address the main features of breast cancer. Thus, the general aspects of this disease have been shown since the first evidence of breast cancer in the world until the numbers today. In this way, there are some ways to prevent breast cancer, such as the woman's lifestyle (healthy eating habits and physical activities) that helps to reduce the incidence of this anomaly. The first noticeable symptom of this anomaly is typically a lump that feels different from the rest of the breast tissue. More than 80% of breast cancer are discovered when the woman feels a lump being present and about 90% of the cases, the cancer is noticed by the woman herself. Currently, the most used method for the detection of cancer and other injuries is the Magnetic Resonance Imaging (MRI) technique. This technique has been shown to be very effective, however, for a better visualization of the images, Contrast Agents (CAs) are used, which are paramagnetic compounds capable of increasing the relaxation of the hydrogen atoms of the water molecules present in the body tissues. The most used CAs are Gd3+ complexes, although they are very efficient, they are toxic to the organism. Thus, new contrast agents have been studied to replace Gd3+ complexes; we can mention iron oxides as a promising substitute.

Deciphering inhibitory activity of flavonoids against tau protein kinases: a coupled molecular docking and quantum chemical study.

Today, Alzheimer's disease (AD) is one of the most important neurodegenerative disorders that affected millions of people worldwide. Hundreds of academic investigations highlighted the potential roles of natural metabolites in the cornerstone of AD prevention. Nevertheless, alkaloids are only metabolites that successfully showed promising clinical therapeutic effects on the prevention of AD. In this regard, other plant metabolites such as flavonoids are also considered as promising substances in the improvement of AD complications. The lack of data on molecular mode of action of flavonoids inside brain tissues, and their potential to transport across the blood-brain barrier, a physical hindrance between bloodstream and brain tissues, limited the large-scale application of these compounds for AD therapy programs. Herein, a coupled docking and quantum study was applied to determine the binding mode of flavonoids and three protein kinases involved in the pathogenesis of AD. The results suggested that all docked metabolites showed considerable binding affinity to interact with target receptors, but some compounds possessed higher binding energy values. Because docking simulation cannot entirely reveal the potential roles of ligand substructures in the interaction with target residues, quantum chemical analyses (QCAs) were performed to cover this drawback. Accordingly, QCAs determined that distribution of molecular orbitals have a pivotal function in the determination of the type of reaction between ligands and receptors; therefore, using such quantum chemical descriptors may correct the results of virtual docking outcomes to highlight promising backbones for further developments.Communicated by Ramaswamy H. Sarma.

In Silico Studies of Potential Selective Inhibitors of Thymidylate Kinase from Variola virus.

Continuing the work developed by our research group, in the present manuscript, we performed a theoretical study of 10 new structures derived from the antivirals cidofovir and ribavirin, as inhibitor prototypes for the enzyme thymidylate kinase from Variola virus (VarTMPK). The proposed structures were subjected to docking calculations, molecular dynamics simulations, and free energy calculations, using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method, inside the active sites of VarTMPK and human TMPK (HssTMPK). The docking and molecular dynamic studies pointed to structures 2, 3, 4, 6, and 9 as more selective towards VarTMPK. In addition, the free energy data calculated through the MM-PBSA method, corroborated these results. This suggests that these compounds are potential selective inhibitors of VarTMPK and, thus, can be considered as template molecules to be synthesized and experimentally evaluated against smallpox.