Pouring SIRAH on NAMD.

Molecular dynamics (MD) simulations provide an invaluable platform for exploring the dynamics of complex biomolecular systems at atomic resolution. However, compatibility issues between force fields and MD software engines can limit the interoperability and transferability of simulations. This work demonstrates the successful use of the coarse-grained SIRAH force field on the widely used NAMD MD engine across a range of increasingly complex biomolecular systems. By leveraging NAMD's ability to read AMBER input files, SIRAH simulations can be run seamlessly on NAMD, including its recently released GPU-accelerated version, NAMD3. The benchmark systems demonstrate consistent results across AMBER, NAMD2, and NAMD3. Thus, these data highlight the enhanced simulation throughput achievable on GPU-accelerated desktop computers using all three engines along with SIRAH. Overall, this study expands the range of the SIRAH force field by utilizing advanced GPU computing resources and high-performance supercomputing facilities, which are particularly effective with NAMD.

Structural Characterization and Rat Aortic Vascular Reactivity of Bradykinin-Potentiating Peptides (BPPs) from the Snake Venom of Bothrops moojeni from Delta do Parnaíba Region, Brazil.

Snake venoms contain various bradykinin-potentiating peptides (BPPs). First studied for their vasorelaxant properties due to angiotensin converting enzyme (ACE) inhibition, these molecules present a range of binding partners, among them the argininosuccinate synthase (AsS) enzyme. This has renewed interest in their characterization from biological sources and the evaluation of their pharmacological activities. In the present work, the low molecular weight fraction of Bothrops moojeni venom was obtained and BPPs were characterized by mass spectrometry. Eleven BPPs or related peptides were sequenced, and one of them, BPP-Bm01, was new. Interestingly, some oxidized BPPs were detected. The three most abundant peptides were BPP-Bm01, BPP-Bax12, and BPP-13a, and their putative interactions with the AsS enzyme were investigated in silico. A binding cavity for these molecules was predicted, and docking studies allowed their ranking. Three peptides were synthesized and submitted to vasorelaxation assays using rat aortic rings. While all BPPs were active, BPP-Bm01 showed the highest potency in this assay. This work adds further diversity to BPPs from snake venoms and suggests, for the first time, a putative binding pocket for these molecules in the AsS enzyme. This can guide the design of new and more potent AsS activators.

Evaluating Known Zika Virus NS2B-NS3 Protease Inhibitor Scaffolds via In Silico Screening and Biochemical Assays.

The NS2B-NS3 protease (NS2B-NS3pro) is regarded as an interesting molecular target for drug design, discovery, and development because of its essential role in the Zika virus (ZIKV) cycle. Although no NS2B-NS3pro inhibitors have reached clinical trials, the employment of drug-like scaffolds can facilitate the screening process for new compounds. In this study, we performed a combination of ligand-based and structure-based in silico methods targeting two known non-peptide small-molecule scaffolds with micromolar inhibitory activity against ZIKV NS2B-NS3pro by a virtual screening (VS) of promising compounds. Based on these two scaffolds, we selected 13 compounds from an initial library of 509 compounds from ZINC15's similarity search. These compounds exhibited structural modifications that are distinct from previously known compounds yet keep pertinent features for binding. Despite promising outcomes from molecular docking and initial enzymatic assays against NS2B-NS3pro, confirmatory assays with a counter-screening enzyme revealed an artifactual inhibition of the assessed compounds. However, we report two compounds, 9 and 11, that exhibited antiviral properties at a concentration of 50 µM in cellular-based assays. Overall, this study provides valuable insights into the ongoing research on anti-ZIKV compounds to facilitate and improve the development of new inhibitors.

Structure-based discovery of thiosemicarbazones as SARS-CoV-2 main protease inhibitors.

Aim: Discovery of novel SARS-CoV-2 main protease (Mpro) inhibitors using a structure-based drug discovery strategy. Materials & methods: Virtual screening employing covalent and noncovalent docking was performed to discover Mpro inhibitors, which were subsequently evaluated in biochemical and cellular assays. Results: 91 virtual hits were selected for biochemical assays, and four were confirmed as reversible inhibitors of SARS CoV-2 Mpro with IC50 values of 0.4-3 µM. They were also shown to inhibit SARS-CoV-1 Mpro and human cathepsin L. Molecular dynamics simulations indicated the stability of the Mpro inhibitor complexes and the interaction of ligands at the subsites. Conclusion: This approach led to the discovery of novel thiosemicarbazones as potent SARS-CoV-2 Mpro inhibitors.

Rad5 HIRAN domain: Structural insights into its interaction with ssDNA through molecular modeling approaches.

The Rad5 protein is an SWI/SNF family ubiquitin ligase that contains an N-terminal HIRAN domain and a RING C3HC4 motif. The HIRAN domain is critical for recognition of the stalled replication fork during the replication process and acts as a sensor to initiate the damaged DNA checkpoint. It is a conserved domain widely distributed in eukaryotic organisms and is present in several DNA-binding proteins from all kingdoms. Here we showed that distant species have important differences in key residues that affect affinity for ssDNA. Based on these findings, we hypothesized that different HIRAN domains might affect fork reversal and translesion synthesis through different metabolic processes. To address this question, we predicted the tertiary structure of both yeast and human HIRAN domains using molecular modeling. Structural dynamics experiments showed that the yeast HIRAN domain exhibited higher structural denaturation than its human homolog, although both domains became stable in the presence of ssDNA. Analysis of atomic contacts revealed that a greater number of interactions between the ssDNA nucleotides and the Rad5 domain are electrostatic. Taken together, these results provide new insights into the molecular mechanism of the HIRAN domain of Rad5 and may guide us to further elucidate differences in the ancient eukaryotes HIRAN sequences and their DNA affinity.Communicated by Ramaswamy H. Sarma.

Structure-Based Identification of Naphthoquinones and Derivatives as Novel Inhibitors of Main Protease Mpro and Papain-like Protease PLpro of SARS-CoV-2.

The worldwide COVID-19 pandemic caused by the coronavirus SARS-CoV-2 urgently demands novel direct antiviral treatments. The main protease (Mpro) and papain-like protease (PLpro) are attractive drug targets among coronaviruses due to their essential role in processing the polyproteins translated from the viral RNA. In this study, we virtually screened 688 naphthoquinoidal compounds and derivatives against Mpro of SARS-CoV-2. Twenty-four derivatives were selected and evaluated in biochemical assays against Mpro using a novel fluorogenic substrate. In parallel, these compounds were also assayed with SARS-CoV-2 PLpro. Four compounds inhibited Mpro with half-maximal inhibitory concentration (IC50) values between 0.41 µM and 9.0 µM. In addition, three compounds inhibited PLpro with IC50 ranging from 1.9 µM to 3.3 µM. To verify the specificity of Mpro and PLpro inhibitors, our experiments included an assessment of common causes of false positives such as aggregation, high compound fluorescence, and inhibition by enzyme oxidation. Altogether, we confirmed novel classes of specific Mpro and PLpro inhibitors. Molecular dynamics simulations suggest stable binding modes for Mpro inhibitors with frequent interactions with residues in the S1 and S2 pockets of the active site. For two PLpro inhibitors, interactions occur in the S3 and S4 pockets. In summary, our structure-based computational and biochemical approach identified novel naphthoquinonal scaffolds that can be further explored as SARS-CoV-2 antivirals.

Bacterial 2'-Deoxyguanosine Riboswitch Classes as Potential Targets for Antibiotics: A Structure and Dynamics Study.

The spread of antibiotic-resistant bacteria represents a substantial health threat. Current antibiotics act on a few metabolic pathways, facilitating resistance. Consequently, novel regulatory inhibition mechanisms are necessary. Riboswitches represent promising targets for antibacterial drugs. Purine riboswitches are interesting, since they play essential roles in the genetic regulation of bacterial metabolism. Among these, class I (2'-dG-I) and class II (2'-dG-II) are two different 2'-deoxyguanosine (2'-dG) riboswitches involved in the control of deoxyguanosine metabolism. However, high affinity for nucleosides involves local or distal modifications around the ligand-binding pocket, depending on the class. Therefore, it is crucial to understand these riboswitches' recognition mechanisms as antibiotic targets. In this work, we used a combination of computational biophysics approaches to investigate the structure, dynamics, and energy landscape of both 2'-dG classes bound to the nucleoside ligands, 2'-deoxyguanosine, and riboguanosine. Our results suggest that the stability and increased interactions in the three-way junction of 2'-dG riboswitches were associated with a higher nucleoside ligand affinity. Also, structural changes in the 2'-dG-II aptamers enable enhanced intramolecular communication. Overall, the 2'-dG-II riboswitch might be a promising drug design target due to its ability to recognize both cognate and noncognate ligands.

pH and non-covalent ligand binding modulate Zika virus NS2B/NS3 protease binding site residues: Discoveries from MD and constant pH MD simulations.

Zika virus (ZIKV) is a global health concern and has been linked to severe neurological pathologies. Although no medication is available yet, many efforts to develop antivirals and host cell binding inhibitors led to attractive drug-like scaffolds, mainly targeting the nonstructural NS2B/NS3 protease (NS2B/NS3pro). NS2B/NS3pro active site has several titratable residues susceptible to pH changes and ligand binding; hence, understanding these residues' protonation is essential to drug design efforts targeting the active site. Here we use in silico methods to probe non-covalent binding and its effect on pKa shifts of the active site residues on a ligand-free protease and with a non-peptidic competitive inhibitor (Ki=13.5 µM). By applying constant pH molecular dynamics, we found that the catalytic residues of the unbound NS2B/NS3pro achieved the protonation needed for the serine protease mechanism over the pH value of 8.5. Nevertheless, the protease in the holo state achieved this same scenario at lower pH values. Also, non-covalent binding affected the catalytic triad (H51, D75, and S135) by stabilizing their distances and interaction network. Thus, NS2B/NS3pro residues configuration for activity might be both pH-dependent and influenced by ligand binding. However, compound presence within the binding site destabilized the NS2B, interfering with the closed and active conformation necessary for substrate binding and catalysis. Our outcomes provide valuable insights into non-covalent inhibitor behavior and its effect on protease active site residues, impacting optimization and design of novel compounds. Communicated by Ramaswamy H. Sarma.

Propedia: a database for protein-peptide identification based on a hybrid clustering algorithm.

BACKGROUND: Protein-peptide interactions play a fundamental role in a wide variety of biological processes, such as cell signaling, regulatory networks, immune responses, and enzyme inhibition. Peptides are characterized by low toxicity and small interface areas; therefore, they are good targets for therapeutic strategies, rational drug planning and protein inhibition. Approximately 10% of the ethical pharmaceutical market is protein/peptide-based. Furthermore, it is estimated that 40% of protein interactions are mediated by peptides. Despite the fast increase in the volume of biological data, particularly on sequences and structures, there remains a lack of broad and comprehensive protein-peptide databases and tools that allow the retrieval, characterization and understanding of protein-peptide recognition and consequently support peptide design. RESULTS: We introduce Propedia, a comprehensive and up-to-date database with a web interface that permits clustering, searching and visualizing of protein-peptide complexes according to varied criteria. Propedia comprises over 19,000 high-resolution structures from the Protein Data Bank including structural and sequence information from protein-peptide complexes. The main advantage of Propedia over other peptide databases is that it allows a more comprehensive analysis of similarity and redundancy. It was constructed based on a hybrid clustering algorithm that compares and groups peptides by sequences, interface structures and binding sites. Propedia is available through a graphical, user-friendly and functional interface where users can retrieve, and analyze complexes and download each search data set. We performed case studies and verified that the utility of Propedia scores to rank promissing interacting peptides. In a study involving predicting peptides to inhibit SARS-CoV-2 main protease, we showed that Propedia scores related to similarity between different peptide complexes with SARS-CoV-2 main protease are in agreement with molecular dynamics free energy calculation. CONCLUSIONS: Propedia is a database and tool to support structure-based rational design of peptides for special purposes. Protein-peptide interactions can be useful to predict, classifying and scoring complexes or for designing new molecules as well. Propedia is up-to-date as a ready-to-use webserver with a friendly and resourceful interface and is available at: https://bioinfo.dcc.ufmg.br/propedia.

Shared Binding Mode of Perrottetinene and Tetrahydrocannabinol Diastereomers inside the CB1 Receptor May Incentivize Novel Medicinal Drug Design: Findings from an in Silico Assay.

In recent years, therapeutic compounds derived from phytocannabinoids have brought renewed attention to the benefits they offer to ameliorate chronic disease symptoms. Among cannabinoids, tetrahydrocannabinol (THC) is a well-known component of the Cannabis plant, whose active principles have been studied through the years. Another psychoactive phytocannabinoid, derived from liverworts Radula, perrottetinene (PET), has created interest, especially as a pharmaceutical product and for its legal recreational use. Unfortunately, so far, the interaction mode of these compounds at the type 1 cannabinoid receptors (CB1R) binding site remains unknown, and no experimental three-dimensional structure in complex with THC or PET is available in the Protein Data Bank. Today, many computational methodologies can assist in this crusade and help unveil how these molecules bind, based on the already known pose of a structurally similar compound. In this work, we aim to elucidate the binding mode of THC and PET molecules in both cis and trans conformers, using a combination of several computational methodologies, including molecular docking, molecular dynamics, free energy calculations, and protein-energy network studies. We found that THC and PET interact similarly with the CB1R, in a different conformation depending on the considered diastereomer. We have observed that cis ligands adopted a half-chair conformation of the cycle ring containing the dimethyl group, assuming an axial or equatorial conformation producing a different induced fitting of the surrounding residues compared with trans ligands, with higher interaction energy than the trans conformer. For PET, we have seen that Trp-279 and Trp-356 have a marked influence on the binding. After binding, Trp-279 accommodates its side chain to better interact with the PET's terminal phenyl group, disturbing CB1R residues communication. The interaction with Trp-356 might impair the activation of CB1R and can influence the binding of PET as a partial agonist. Understanding the PET association with CB1R from a molecular perspective can offer a glimpse of preventing potential toxicological or recreational effects since it is an attractive lead for drug development with fewer side effects than trans-THC.

Proteus: An algorithm for proposing stabilizing mutation pairs based on interactions observed in known protein 3D structures.

BACKGROUND: Protein engineering has many applications for industry, such as the development of new drugs, vaccines, treatment therapies, food, and biofuel production. A common way to engineer a protein is to perform mutations in functionally essential residues to optimize their function. However, the discovery of beneficial mutations for proteins is a complex task, with a time-consuming and high cost for experimental validation. Hence, computational approaches have been used to propose new insights for experiments narrowing the search space and reducing the costs. RESULTS: In this study, we developed Proteus (an acronym for Protein Engineering Supporter), a new algorithm for proposing mutation pairs in a target 3D structure. These suggestions are based on contacts observed in other known structures from Protein Data Bank (PDB). Proteus' basic assumption is that if a non-interacting pair of amino acid residues in the target structure is exchanged to an interacting pair, this could enhance protein stability. This trade is only allowed if the main-chain conformation of the residues involved in the contact is conserved. Furthermore, no steric impediment is expected between the proposed mutations and the surrounding protein atoms. To evaluate Proteus, we performed two case studies with proteins of industrial interests. In the first case study, we evaluated if the mutations suggested by Proteus for four protein structures enhance the number of inter-residue contacts. Our results suggest that most mutations proposed by Proteus increase the number of interactions into the protein. In the second case study, we used Proteus to suggest mutations for a lysozyme protein. Then, we compared Proteus' outcomes to mutations with available experimental evidence reported in the ProTherm database. Four mutations, in which our results agree with the experimental data, were found. This could be initial evidence that changes in the side-chain of some residues do not cause disturbances that harm protein structure stability. CONCLUSION: We believe that Proteus could be used combined with other methods to give new insights into the rational development of engineered proteins. Proteus user-friendly web-based tool is available at < http://proteus.dcc.ufmg.br >.

Glutantßase: a database for improving the rational design of glucose-tolerant ß-glucosidases.

Β-glucosidases are key enzymes used in second-generation biofuel production. They act in the last step of the lignocellulose saccharification, converting cellobiose in glucose. However, most of the ß-glucosidases are inhibited by high glucose concentrations, which turns it a limiting step for industrial production. Thus, ß-glucosidases have been targeted by several studies aiming to understand the mechanism of glucose tolerance, pH and thermal resistance for constructing more efficient enzymes. In this paper, we present a database of ß-glucosidase structures, called Glutantßase. Our database includes 3842 GH1 ß-glucosidase sequences collected from UniProt. We modeled the sequences by comparison and predicted important features in the 3D-structure of each enzyme. Glutantßase provides information about catalytic and conserved amino acids, residues of the coevolution network, protein secondary structure, and residues located in the channel that guides to the active site. We also analyzed the impact of beneficial mutations reported in the literature, predicted in analogous positions, for similar enzymes. We suggested these mutations based on six previously described mutants that showed high catalytic activity, glucose tolerance, or thermostability (A404V, E96K, H184F, H228T, L441F, and V174C). Then, we used molecular docking to verify the impact of the suggested mutations in the affinity of protein and ligands (substrate and product). Our results suggest that only mutations based on the H228T mutant can reduce the affinity for glucose (product) and increase affinity for cellobiose (substrate), which indicates an increment in the resistance to product inhibition and agrees with computational and experimental results previously reported in the literature. More resistant ß-glucosidases are essential to saccharification in industrial applications. However, thermostable and glucose-tolerant ß-glucosidases are rare, and their glucose tolerance mechanisms appear to be related to multiple and complex factors. We gather here, a set of information, and made predictions aiming to provide a tool for supporting the rational design of more efficient ß-glucosidases. We hope that Glutantßase can help improve second-generation biofuel production. Glutantßase is available at http://bioinfo.dcc.ufmg.br/glutantbase .

Structural insights into NS5B protein of novel equine hepaciviruses and pegiviruses complexed with polymerase inhibitors.

Infections produced by hepaciviruses have been associated with liver disease in horses. Currently, at least three viruses belonging to the Flaviviridae family are capable of producing a chronic infection in equines: non-primate hepacivirus (NPHV), Theiler's disease-associated virus (TDAV), and equine pegivirus (EPgV). The RNA-dependent RNA polymerases of viruses (RdRp) (NS5 protein), from the flavivirus family, use de novo RNA synthesis to initiate synthesis. The two antiviral drugs currently used to treat hepatitis C (HCV), sofosbuvir and dasabuvir, act on the viral NS5B polymerase as nucleoside and non-nucleoside inhibitors, respectively. Both drugs have shown significant clinical inhibition of viral response. In this work, we aimed to model the NS5B polymerase of the equine hepacivirus (EHCV) subtypes 1 and 2, TDAV and EPgV, to assess whether current direct-acting antiviral drugs against HCV interact with these proteins. Crystal structures of HCV-NS5B were used as templates for modeling target sequences in both conformations (open and closed). Also, molecular docking of sofosbuvir and dasabuvir were performed to predict their possible binding modes at the modeled NS5B polymerase binding sites. We observed that the NS5B models of the EHCV and EPgV shared well-conserved 3D structures to HCV-NS5B and other RdRps, suggesting functional conservation. Interactions of EHCV subtypes 1, 2 and TDAV polymerases with sofosbuvir showed a similar molecular interaction pattern compared to HCV-NS5B, while interactions with dasabuvir were less conserved. In silico studies of molecular interactions between these modeled structures and sofosbuvir suggest that this compound could be efficient in combating equine pathogens, thus contributing to animal welfare.

Studying effects of different protonation states of His11 and His102 in ribose-5-phosphate isomerase of Trypanosoma cruzi: an example of cooperative behavior.

The Trypanosoma cruzi ribose-5-phosphate isomerase B (TcRpiB) is a crucial piece in the pentose phosphate pathway and thus is a potential drug target for treatment of Chagas' disease. TcRpiB residues, such as Cys69, Asp45, Glu149 and Pro47, have confirmed their roles in substrate recognition, catalytic reaction and binding site conformation. However, the joint performance of His11 and His102, in the D-ribose-5-phosphate (R5P) in the catalysis is not well understood. In this work, we probed the influence of different protonation states of His11 and His102 on the behavior of the ligand R5P using molecular dynamics simulations, network analysis and thermodynamic integration. Simulations revealed that a protonated His11 combined with a neutral His102 (His11+‒His102) was able to stabilize the ligand R5P in the binding site. Moreover, calculated relative free energy differences showed that when protonated His11 was coupled to a neutral His102 an exergonic process takes place. On the other hand, neutral His11 combined with a protonated His102 (His11‒His102+), sampled conformations that resembled the catalyzed product D-ribulose-5-phosphate (Ru5P). Network analysis also demonstrated some peculiarities for these systems with some negatively correlated nodes in the binding site for His11‒His102+, and exclusive suboptimal paths for His11+‒His102. Therefore, the combined approach presented in this paper proposes two suitable protonation states for the TcRpiB catalytic mechanism, where an extra proton in either histidines might favor R5P binding or influence isomerization reaction to Ru5P. Our results may guide further in silico drug discovery studies. Communicated by Ramaswamy H. Sarma.

Integrating Molecular Docking and Molecular Dynamics Simulations.

Computational methods, applied at the early stages of the drug design process, use current technology to provide valuable insights into the understanding of chemical systems in a virtual manner, complementing experimental analysis. Molecular docking is an in silico method employed to foresee binding modes of small compounds or macromolecules in contact with a receptor and to predict their molecular interactions. Moreover, the methodology opens up the possibility of ranking these compounds according to a hierarchy determined using particular scoring functions. Docking protocols assign many approximations, and most of them lack receptor flexibility. Therefore, the reliability of the resulting protein-ligand complexes is uncertain. The association with the costly but more accurate MD techniques provides significant complementary with docking. MD simulations can be used before docking since a series of "new" and broader protein conformations can be extracted from the processing of the resulting trajectory and employed as targets for docking. They also can be utilized a posteriori to optimize the structures of the final complexes from docking, calculate more detailed interaction energies, and provide information about the ligand binding mechanism. Here, we focus on protocols that offer the docking-MD combination as a logical approach to improving the drug discovery process.

Benzimidazole inhibitors of the major cysteine protease of Trypanosoma brucei.

Aim: Limitations in available therapies for trypanosomiases indicate the need for improved medicines. Cysteine proteases cruzain and rhodesain are validated targets for treatment of Chagas disease and human African trypanosomiasis. Previous studies reported a benzimidazole series as potent cruzain inhibitors. Results & methodology: Considering the high similarity between these proteases, we evaluated 40 benzimidazoles against rhodesain. We describe their structure-activity relationships (SAR), revealing trends similar to those observed for cruzain and features that lead to enzyme selectivity. This series comprises noncovalent competitive inhibitors (best Ki = 0.21 µM against rhodesain) and micromolar activity against Trypanosoma brucei brucei. A cheminformatics analysis confirms scaffold novelty, and the inhibitors described have favorable predicted physicochemical properties. Conclusion: Our results support this series as a starting point for new human African trypanosomiasis medicines.

Thiophenacetamide as a potential modulator to NF-κB: structure and dynamics study using in silico and molecular biology assays.

Nuclear factor kappa B (NF-κB) plays critical roles in the regulation of many pathophysiological processes, including inflammation and immune responses, cell growth and apoptosis. This DNA-binding protein receptor is considered an important molecular target to treat many diseases through host-directed therapy. In this line, several drugs containing thiophene cores have been extensively evaluated due to their ability to interfere on NF-κB translocation to the nucleus. In this work, assays using drug affinity responsive target stability (DARTS) revealed that the parent compound N-(Aryl)-2-thiophen-2-ylacetamide referred to as thiophenacetamide (TAA) specifically binds to the p65 subunit of the NF-κB. Since no experimental binding mode of TAA with p65 is available, we explored TAA within putative sites in silico to gain insights into its possible binding mode and behavior. The binding mode of TAA found in Site 1 formed hydrogen bonds with Lys37 and Asp125 on p65, important residues near DNA-binding region. Molecular dynamics simulations showed the stability of this mode of binding in contrast to the other also tested modes. Our results suggest that TAA binding could occur in regions close to residues responsible for DNA binding, increasing NF-κB protein rigidity and affecting the association between DNA and NF-κB. Communicated by Ramaswamy H. Sarma.

Understanding Structure-Activity Relationships for Trypanosomal Cysteine Protease Inhibitors by Simulations and Free Energy Calculations.

The protozoan cysteine proteases cruzain in Trypanosoma cruzi and rhodesain in Trypanosoma brucei are therapeutic targets for Chagas disease and Human African Trypanosomiasis (HAT), respectively. A benzimidazole series was previously characterized as potent noncovalent competitive cruzain and rhodesain inhibitors with activity against trypanosomes. Common structure-activity relationships (SAR) trends and structural modifications leading to selectivity against each enzyme were described. However, some of these trends could not be understood based on the reported binding mode of lead compound 1. Therefore, we employed microsecond molecular dynamics simulations and free energy calculations to understand qualitative SAR trends and to quantitatively recapitulate them. Simulations revealed the most stable protein-ligand interactions and provided insights concerning enzyme selectivity. Calculated relative binding free energies of compound 1 analogs exhibited deviations of 1.1 and 2.2 kcal/mol from the experimental values for cruzain and rhodesain, respectively. These data encourage prospective thermodynamic integration (TI) studies to optimize this series and facilitate the prioritization of compounds for synthesis.