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.

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.

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. 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 (ΔG_bind) 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 ΔG_bind values. Two compounds with the best ΔG_bind 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. 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. Fascioliosis is considered an emerging disease in humans, causing important losses in global agriculture through the infection of livestock animals. The outcome of resistant parasites has increased the search for new drugs which may contribute to disease control. In recent decades, Fasciola cathepsins (FhCs) have been defined as the principal virulence factors of this parasite. Despite being in the same protein family, they have different specificities and, thus, distinct roles throughout the fluke life cycle. Differences in specificity have been attributed to a few variations in the sequence of key FhCs subsites. Currently, the structure-based drug design of inhibitors against Fasciola cathepsin Ls (FhCLs) with unknown structures is possible due to the availability of the three-dimensional structure of FhCL1. Our detailed structural analysis of the major infective juvenile enzyme (FhCL3) identifies the molecular determinants for protein binding. Also, novel potential inhibitors against FhCL3 are proposed, which might reduce host invasion and penetration processes. These compounds are predicted to interact with the binding site of the enzyme, therefore they could prevent substrate processing by competitive inhibition. The structure-based drug design strategy described here will be useful for the development of new potent and selective inhibitors against other FhCs.

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.