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.

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 pK 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.

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.Abbreviations: AA, acrylamide; CAS, computational alanine scanning; CD, circular dichroism; CF, carboxyfluorescein; EM, energy minimization; EqtII, equinatoxin II; FraC, fragaceatoxin C; DG bind , binding free energy; DG polar , polar desolvation energy; DG nonpolar , nonpolar desolvation energy; DG solvat , solvation free energy; HA, hemolytic activity; Ksv, SternVolmer constant; m, slope of regression fit/line; MD, molecular dynamic; MLV, multilamellar vesicles; PFP, pore-forming protein; PFT, pore-forming toxins; POPC, 1-palmitoyl-2-oleylphosphatidylcholine; SM, sphingomyelin; StII, sticholysin II; SUV, small unilamellar vesicles; DV ele , electrostatic energy; DV vdw , van der Waals Energy; p, surface pressure; p 0 , initial surface pressure; Dp, increment in surface pressure. 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.