Search

The calculation of relative free energies (ΔΔ G ) for charge-changing mutations at protein-protein interfaces through alchemical methods remains challenging due to variations in the system's net charge during charging steps, the possibility of mutated and contacting ionizable residues occurring in various protonation states, and undersampling issues. In this study, we present a set of strategies, collectively termed TIRST/TIRST-H + , to address some of these challenges. Our approaches combine thermodynamic integration (TI) with the prediction of p K a shifts to calculate ΔΔ G values. Moreover, special sets of restraints are employed to keep the alchemically transformed molecules separated. The accuracy of the devised approaches was assessed on a large and diverse data set comprising 164 point mutations of charged residues (Asp, Glu, Lys, and Arg) to Ala at the protein-protein interfaces of complexes with known three-dimensional structures. Mean absolute and root-mean-square errors ranging from 1.38 to 1.66 and 1.89 to 2.44 kcal/mol, respectively, and Pearson correlation coefficients of ∼0.6 were obtained when testing the approaches on the selected data set using the GPU-TI module of Amber18 suite and the ff14SB force field. Furthermore, the inclusion of variable protonation states for the mutated acid residues improved the accuracy of the predicted ΔΔ G values. Therefore, our results validate the use of TIRST/TIRST-H + in prospective studies aimed at evaluating the impact of charge-changing mutations to Ala on the stability of protein-protein complexes.

Spiders of Loxosceles genus are widely distributed and their venoms contain phospholipases D (PLDs), which degrade phospholipids and trigger inflammatory responses, dermonecrosis, hematological changes, and renal injuries. Biochemical, functional, and structural properties of three recombinant PLDs from L. intermedia , L. laeta, and L. gaucho , the principal species clinically relevant in South America, were analyzed. Sera against L. gaucho and L. laeta PLDs strongly cross-reacted with other PLDs, but sera against L. intermedia PLD mostly reacted with homologous molecules, suggesting underlying structural and functional differences. PLDs presented a similar secondary structure profile but distinct melting temperatures. Different methods demonstrated that all PLDs cleave sphingomyelin and lysophosphatidylcholine, but L. gaucho and L. laeta PLDs excelled. L. gaucho PLD showed greater "in vitro" hemolytic activity. L. gaucho and L. laeta PLDs were more lethal in assays with mice and crickets. Molecular dynamics simulations correlated their biochemical activities with differences in sequences and conformations of specific surface loops, which play roles in protein stability and in modulating interactions with the membrane. Despite the high similarity, PLDs from L. gaucho and L. laeta venoms are more active than L. intermedia PLD, requiring special attention from physicians when these two species prevail in endemic regions.

Brown spider envenomation results in dermonecrosis, characterized by an intense inflammatory reaction. The principal toxins of brown spider venoms are phospholipase-D isoforms, which interact with different cellular membrane components, degrade phospholipids, and generate bioactive mediators leading to harmful effects. The Loxosceles intermedia phospholipase D, LiRecDT1, possesses a loop that modulates the accessibility to the active site and plays a crucial role in substrate. In vitro and in silico analyses were performed to determine aspects of this enzyme's substrate preference. Sphingomyelin d18:1/6:0 was the preferred substrate of LiRecDT1 compared to other Sphingomyelins. Lysophosphatidylcholine 16:0/0:0 was preferred among other lysophosphatidylcholines, but much less than Sphingomyelin d18:1/6:0. In contrast, phosphatidylcholine d18:1/16:0 was not cleaved. Thus, the number of carbon atoms in the substrate plays a vital role in determining the optimal activity of this phospholipase-D. The presence of an amide group at C2 plays a key role in recognition and activity. In silico analyses indicated that a subsite containing the aromatic residues Y228 and W230 appears essential for choline recognition by cation-π interactions. These findings may help to explain why different cells, with different phospholipid fatty acid compositions exhibit distinct susceptibilities to brown spider venoms.

Staphylococcal exfoliative toxins (ETs) are glutamyl endopeptidases that specifically cleave the Glu381-Gly382 bond in the ectodomains of desmoglein 1 (Dsg1) via complex action mechanisms. To date, four ETs have been identified in different Staphylococcus aureus strains and ETE is the most recently characterized. The unusual properties of ETs have been attributed to a unique structural feature, i.e., the 180° flip of the carbonyl oxygen (O) of the nonconserved residue 192/186 (ETA/ETE numbering), not conducive to the oxyanion hole formation. We report the crystal structure of ETE determined at 1.61 Å resolution, in which P186(O) adopts two conformations displaying a 180° rotation. This finding, together with free energy calculations, supports the existence of a dynamic transition between the conformations under the tested conditions. Moreover, enzymatic assays showed no significant differences in the esterolytic efficiency of ETE and ETE/P186G, a mutant predicted to possess a functional oxyanion hole, thus downplaying the influence of the flip on the activity. Finally, we observed the formation of ETE homodimers in solution and the predicted homodimeric structure revealed the participation of a characteristic nonconserved loop in the interface and the partial occlusion of the protein active site, suggesting that monomerization is required for enzymatic activity.

The SARS-CoV-2 main protease, also known as 3-chymotrypsin-like protease (3CL pro ), is a cysteine protease responsible for the cleavage of viral polyproteins pp1a and pp1ab, at least, at eleven conserved sites, which leads to the formation of mature nonstructural proteins essential for the replication of the virus. Due to its essential role, numerous studies have been conducted so far, which have confirmed 3CL pro as an attractive drug target to combat Covid-19 and have reported a vast number of inhibitors and their co-crystal structures. Despite all the ongoing efforts, D-peptides, which possess key advantages over L-peptides as therapeutic agents, have not been explored as potential drug candidates against 3CL pro . The current work fills this gap by reporting an in silico approach for the discovery of D-peptides capable of inhibiting 3CL pro that involves structure-based virtual screening (SBVS) of an in-house library of D-tripeptides and D-tetrapeptides into the protease active site and subsequent rescoring steps, including Molecular Mechanics Generalized-Born Surface Area (MM-GBSA) free energy calculations and molecular dynamics (MD) simulations. In vitro enzymatic assays conducted for the four top-scoring D-tetrapeptides at 20 μM showed that all of them caused 55-85% inhibition of 3CL pro activity, thus highlighting the suitability of the devised approach. Overall, our results present a promising computational strategy to identify D-peptides capable of inhibiting 3CL pro , with broader application in problems involving protein inhibition., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Hernández González, Eberle, Willbold and Coronado.)

Allosteric inhibitors regulate enzyme activity from remote and usually specific pockets. As they promise an avenue for less toxic and safer drugs, the identification and characterization of allosteric inhibitors has gained great academic and biomedical interest in recent years. Research on falcipain-2 (FP-2), the major papain-like cysteine hemoglobinase of Plasmodium falciparum , might benefit from this strategy to overcome the low selectivity against human cathepsins shown by active site-directed inhibitors. Encouraged by our previous finding that methacycline inhibits FP-2 noncompetitively, here we assessed other five tetracycline derivatives against this target and characterized their inhibition mechanism. As previously shown for methacycline, tetracycline derivatives inhibited FP-2 in a noncompetitive fashion, with K i values ranging from 121 to 190 μM. A possible binding to the S' side of the FP-2 active site, similar to that described by X-ray crystallography (PDB: 6SSZ) for the noncompetitive inhibitor E- chalcone 48 (EC48), was experimentally discarded by kinetic analysis using a large peptidyl substrate spanning the whole active site. By combining lengthy molecular dynamics (MD) simulations that allowed methacycline to diffuse from solution to different FP-2 surface regions and free energy calculations, we predicted the most likely binding mode of the ligand. Of note, the proposed binding pose explains the low differences in K i values observed for the tested tetracycline derivatives and the calculated binding free energies match the experimental values. Overall, this study has implications for the design of novel allosteric inhibitors against FP-2 and sets the basis for further optimization of the tetracycline scaffold to produce more potent and selective inhibitors.

Falcipain-2 (FP-2) is a Plasmodium falciparum hemoglobinase widely targeted in the search for antimalarials. FP-2 can be allosterically modulated by various noncompetitive inhibitors that have been serendipitously identified. Moreover, the crystal structures of two inhibitors bound to an allosteric site, termed site 6, of the homolog enzyme human cathepsin K (hCatK) suggest that the equivalent region in FP-2 might play a similar role. Here, we conduct the rational identification of FP-2 inhibitors through virtual screenings (VS) of compounds into several pocket-like conformations of site 6, sampled during molecular dynamics (MD) simulations of the free enzyme. Two noncompetitive inhibitors, ZINC03225317 and ZINC72290660, were confirmed using in vitro enzymatic assays and their poses into site 6 led to calculated binding free energies matching the experimental ones. Our results provide strong evidence about the allosteric inhibition of FP-2 through binding of small molecules to site 6, thus opening the way toward the discovery of new inhibitors against this enzyme., (© 2021. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

During their life cycle, Leishmania parasites display a fine-tuned regulation of the mRNA translation through the differential expression of isoforms of eukaryotic translation initiation factor 4E (LeishIF4Es). The interaction between allosteric modulators such as 4E-interacting proteins (4E-IPs) and LeishIF4E affects the affinity of this initiation factor for the mRNA cap. Here, several computational approaches were employed to elucidate the molecular bases of the previously-reported allosteric modulation in L. major exerted by 4E-IP1 (Lm4E-IP1) on eukaryotic translation initiation factor 4E 1 (LmIF4E-1). Molecular dynamics (MD) simulations and accurate binding free energy calculations (Δ G bind ) were combined with network-based modeling of residue-residue correlations. We also describe the differences in internal motions of LmIF4E-1 apo form, cap-bound, and Lm4E-IP1-bound systems. Through community network calculations, the differences in the allosteric pathways of allosterically-inhibited and active forms of LmIF4E-1 were revealed. The Δ G bind values show significant differences between the active and inhibited systems, which are in agreement with the available experimental data. Our study thoroughly describes the dynamical perturbations of LmIF4E-1 cap-binding site triggered by Lm4E-IP1. These findings are not only essential for the understanding of a critical process of trypanosomatids' gene expression but also for gaining insight into the allostery of eukaryotic IF4Es, which could be useful for structure-based design of drugs against this protein family., Competing Interests: The authors declare no competing interests., (© 2021 The Authors.)

Falcipain-2 (FP-2) is hemoglobinase considered an attractive drug target of Plasmodium falciparum . Recently, it has been shown that peptidomimetic nitriles containing a 3-pyridyl (3Pyr) moiety at P2 display high affinity and selectivity for FP-2 with respect to human cysteine cathepsins (hCats), outperforming other P2-Pyr isomers and analogs. Further characterization demonstrated that certain P3 variants of these compounds possess micromolar inhibition of parasite growth in vitro and no cytotoxicity against human cell lines. However, the structural determinants underlying the selectivity of the 3Pyr-containing nitriles for FP-2 remain unknown. In this work, we conduct a thorough computational study combining MD simulations and free energy calculations to decipher the bases of the selectivity of the aforementioned nitriles. Our results reveal that water bridges involving the nitrogen and one carboxyl oxygen of I85 and D234 of FP-2, respectively, and the nitrogen of the neutral 3Pyr moiety, which are either less prevalent or nonexistent in the other complexes, explain the experimental activity profiles. The presence of crystallographic waters close to the bridging water positions in the studied proteases strongly supports the occurrence of such interactions. Overall, our findings suggest that selective FP-2 inhibitors can be designed by promoting water bridge formation at the bottom of the S2 subsite and/or by introducing complementary groups that displace the bridging water.

The arginine repressor (ArgR) regulates the expression of genes involved in arginine biosynthesis. Upon attaining a threshold concentration of arginine in the cytoplasm, the trimeric C-terminal domain of ArgR binds three arginines in a shallow surface cleft and subsequently hexamerizes forming a dimer of trimers containing six Arg co-repressor molecules which are buried at the subunit interfaces. The N-terminal domains of this complex bind to the DNA promoter thereby interrupting the transcription of the genes related to Arg biosynthesis. The crystal structures of the wild type and mutant Pro115Gln ArgR from Corynebacterium pseudotuberculosis determined at 1.7 Å demonstrate that a single amino acid substitution switches co-repressor specificity from Tyr to Arg. Molecular dynamics simulations indicate that the first step, i.e., the binding of the co-repressor, occurs in the trimeric state and that Pro115Gln ArgR preferentially binds Arg. It was also shown that, in Pro115 ArgR hexamers, the concomitant binding of sodium ions shifts selectivity to Tyr. Structural data combined with phylogenetic analyses of ArgR from C. pseudotuberculosis suggest that substitutions in the binding pocket at position 115 may alter its specificity for amino acids and that the length of the protein interdomain linker can provide further functional flexibility. These results support the existence of alternative ArgR regulatory mechanisms in this pathogenic bacterium., Competing Interests: Declaration of Competing Interest None., (Copyright © 2020 Elsevier B.V. All rights reserved.)

Falcipain-2 (FP-2) is a Plasmodium falciparum cysteine protease that has been extensively targeted to identify novel antimalarials. Remarkably, previous reports have shown that FP-2 can be allosterically modulated and, for a particular noncompetitive chalcone inhibitor, the existing lines of experimental evidence can guide the prediction of its unknown binding mode to the enzyme in a reliable fashion. In this work, we propose a structure of FP-2 in complex with the aforementioned compound that fulfills all of the experimental data, by employing a combination of molecular modeling tools, such as pocket volume measurements, docking, molecular dynamics (MD) simulations, and free energy calculations. Our results show that the studied inhibitor binds a transient pocket occluded in all of the available FP-2 crystal structures and lying in a region previously characterized as a potential allosteric site in related cysteine proteases. In addition, we detected in silico the occurrence of significant community reorganization in FP-2, increased signal transmission between the allosteric pocket and the active site, and change in loop motions and residue p K a values upon the compound binding, thus providing insight into the uncharacterized allosteric mechanism. Overall, this study yields valuable predictions for the design of novel allosteric inhibitors against FP-2 and other cysteine proteases.

Trypanosoma cruzi is the causative agent of Chagas disease, a neglected infection affecting millions of people in tropical regions. There are several chemotherapeutic agents for the treatment of this disease, but most of them are highly toxic and generate resistance. Currently, the development of allosteric inhibitors constitutes a promising research field, since it can improve the accessibility to more selective and less toxic medicines. To date, the allosteric drugs prediction is a state-of-the-art topic in rational structure-based computational design. In this work, a simulation strategy was developed for computational discovery of allosteric inhibitors, and it was applied to cruzain, a promising target and the major cysteine protease of T. cruzi. Molecular dynamics simulations, binding free energy calculations and network-based modelling of residue interactions were combined to characterize and compare molecular distinctive features of the apo form and the cruzain-allosteric inhibitor complexes. By using geometry-based criteria on trajectory snapshots, we predicted two main allosteric sites suitable for drug targeting. The results suggest dissimilar mechanisms exerted by the same allosteric site when binding different potential allosteric inhibitors. Finally, we identified the residues involved in suboptimal paths linking the identified site and the orthosteric site. The present study constitutes the first approximation to the design of cruzain allosteric inhibitors and may serve for future pharmacological intervention. Here, no major effects on active site structure were observed due to compound binding (modification of distance and angles between catalytic residues), which indicates that allosteric regulation in cruzain might be mediated via alterations of its dynamical properties similarly to allosteric inhibition of human cathepsin K (HCatK). The current findings are particularly relevant for the design of allosteric modulators of papain-like cysteine proteases., Competing Interests: The authors have declared that no competing interests exist.

Background: Falcipain 2 (FP-2) is the hemoglobin-degrading cysteine protease of Plasmodium falciparum most extensively targeted to develop novel antimalarials. However, no commercial antimalarial drugs based on FP-2 inhibition are available yet due to the low selectivity of most FP-2 inhibitors against the human cysteine proteases., Methods: A structure-based virtual screening (SVBS) using Maybridge HitFinder™ compound database was conducted to identify potential FP-2 inhibitors. In vitro enzymatic and cell-growth inhibition assays were performed for the top-scoring compounds. Docking, molecular dynamics (MD) simulations and free energy calculations were employed to study the interaction of the best hits with FP-2 and other related enzymes., Results and Conclusions: Two hits based on 4-(9H-fluoren-9-yl) piperazin-1-yl) methanone scaffold, HTS07940 and HTS08262, were identified as inhibitors of FP-2 (half-maximal inhibitory concentration (IC 50 ) = 64 μM and 14.7 μM, respectively) without a detectable inhibition against the human off-target cathepsin K (hCatK). HTS07940 and HTS08262 inhibited the growth of the multidrug-resistant P. falciparum strain FCR3 in culture (half-maximal inhibitory concentrations (IC50) = 2.91 μM and 34 μM, respectively) and exhibited only moderate cytotoxicity against HeLa cells (Half-maximal cytotoxic concentration (CC50) = 133 μM and 350 μM, respectively). Free energy calculations reproduced the experimental affinities of the hits for FP-2 and explained the selectivity with respect to hCatK., General Significance: To the best of our knowledge, HTS07940 stands among the most selective FP-2 inhibitors identified by SBVS reported so far, displaying moderate antiplasmodial activity and low cytotoxicity against human cells. Hence, this compound constitutes a promising lead for the design of more potent and selective FP-2 inhibitors., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Falcipain-2 (FP-2) is a major hemoglobinase of Plasmodium falciparum, considered an important drug target for the development of antimalarials. A previous study reported a novel series of 20 reversible peptide-based inhibitors of FP-2. However, the lack of tridimensional structures of the complexes hinders further optimization strategies to enhance the inhibitory activity of the compounds. Here we report the prediction of the binding modes of the aforementioned inhibitors to FP-2. A computational approach combining previous knowledge on the determinants of binding to the enzyme, docking, and postdocking refinement steps, is employed. The latter steps comprise molecular dynamics simulations and free energy calculations. Remarkably, this approach leads to the identification of near-native ligand conformations when applied to a validation set of protein-ligand structures. Overall, we proposed substrate-like binding modes of the studied compounds fulfilling the structural requirements for FP-2 binding and yielding free energy values that correlated well with the experimental data. Proteins 2017; 85:1666-1683. © 2017 Wiley Periodicals, Inc., (© 2017 Wiley Periodicals, Inc.)

Crystallographic data of the dimeric and octameric forms of fragaceatoxin C (FraC) suggested the key role of a small hydrophobic protein-protein interaction surface for actinoporins oligomerization and pore formation in membranes. However, site-directed mutagenesis studies supporting this hypothesis for others actinoporins are still lacking. Here, we demonstrate that disrupting the key hydrophobic interaction between V60 and F163 (FraC numbering scheme) in the oligomerization interface of FraC, equinatoxin II (EqtII), and sticholysin II (StII) impairs the pore formation activity of these proteins. Our results allow for the extension of the importance of FraC protein-protein interactions in the stabilization of the oligomeric intermediates of StII and EqtII pointing out that all of these proteins follow a similar pathway of membrane disruption. These findings support the hybrid pore proposal as the universal model of actinoporins pore formation. Moreover, we reinforce the relevance of dimer formation, which appears to be a functional intermediate in the assembly pathway of some different pore-forming proteins., (© 2016 The Protein Society.)

The Kunitz-type protease inhibitor ShPI-1 inhibits human neutrophil elastase (HNE, Ki = 2.35·10-8 M) but does not interact with the porcine pancreatic elastase (PPE); whereas its P1 site variant, ShPI-1/K13L, inhibits both HNE and PPE (Ki = 1.3·10-9 M, and Ki = 1.2·10-8 M, respectively). By employing a combination of molecular modeling tools, e.g., structural alignment, molecular dynamics simulations and Molecular Mechanics Generalized-Born/Poisson-Boltzmann Surface Area free energy calculations, we showed that D226 of HNE plays a critical role in the interaction of this enzyme with ShPI-1 through the formation of a strong salt bridge and hydrogen bonds with K13 at the inhibitor's P1 site, which compensate the unfavorable polar-desolvation penalty of the latter residue. Conversely, T226 of PPE is unable to establish strong interactions with K13, thereby precluding the insertion of K13 side-chain into the S1 subsite of this enzyme. An alternative conformation of K13 site-chain placed at the entrance of the S1 subsite of PPE, similar to that observed in the crystal structure of ShPI-1 in complex with chymotrypsin (PDB: 3T62), is also unfavorable due to the lack of stabilizing pair-wise interactions. In addition, our results suggest that the higher affinity of ShPI-1/K13L for both elastases mainly arises from the lower polar-desolvation penalty of L13 compared to that of K13, and not from stronger pair-wise interactions of the former residue with those of each enzyme. These results provide insights into the PPE and HNE inhibition and may contribute to the design of more potent and/or specific inhibitors toward one of these proteases.

Background: Fasciola hepatica is the causative agent of fascioliasis, a disease affecting grazing animals, causing economic losses in global agriculture and currently being an important human zoonosis. Overuse of chemotherapeutics against fascioliasis has increased the populations of drug resistant parasites. F. hepatica cathepsin L3 is a protease that plays important roles during the life cycle of fluke. Due to its particular collagenolytic activity it is considered an attractive target against the infective phase of F. hepatica., Methodology/principal Findings: Starting with a three dimensional model of FhCL3 we performed a structure-based design of novel inhibitors through a computational study that combined virtual screening, molecular dynamics simulations, and binding free energy (ΔGbind) calculations. Virtual screening was carried out by docking inhibitors obtained from the MYBRIDGE-HitFinder database inside FhCL3 and human cathepsin L substrate-binding sites. On the basis of dock-scores, five compounds were predicted as selective inhibitors of FhCL3. Molecular dynamic simulations were performed and, subsequently, an end-point method was employed to predict ΔGbind values. Two compounds with the best ΔGbind values (-10.68 kcal/mol and -7.16 kcal/mol), comparable to that of the positive control (-10.55 kcal/mol), were identified. A similar approach was followed to structurally and energetically characterize the interface of FhCL3 in complex with a peptidic substrate. Finally, through pair-wise and per-residue free energy decomposition we identified residues that are critical for the substrate/ligand binding and for the enzyme specificity., Conclusions/significance: The present study is the first computer-aided drug design approach against F. hepatica cathepsins. Here we predict the principal determinants of binding of FhCL3 in complex with a natural substrate by detailed energetic characterization of protease interaction surface. We also propose novel compounds as FhCL3 inhibitors. Overall, these results will foster the future rational design of new inhibitors against FhCL3, as well as other F. hepatica cathepsins.