Computational Perspectives into Plasmepsins Structure—Function Relationship: Implications to Inhibitors Design

The development of efficient and selective antimalariais remains a challenge for the pharmaceutical industry. The aspartic proteases plasmepsins, whose inhibition leads to parasite death, are classified as targets for the design of potent drugs. Combinatorial synthesis is currently being used to generate inhibitor libraries for these enzymes, and together with computational methodologies have been demonstrated capable for the selection of lead compounds. The high structural flexibility of plasmepsins, revealed by their X-ray structures and molecular dynamics simulations, made even more complicated the prediction of putative binding modes, and therefore, the use of common computational tools, like docking and free-energy calculations. In this review, we revised the computational strategies utilized so far, for the structure-function relationship studies concerning the plasmepsin family, with special focus on the recent advances in the improvement of the linear interaction estimation (LIE) method, which is one of the most successful methodologies in the evaluation of plasmepsin-inhibitor binding affinity

Modeling the Binding of Inhibitors/Drugs to the Human Serotonin Transporter

Human serotonin transporter (hSERT), a membrane protein from the neurotransmitter sodium symporter family, is implicated in depression disorder and has been the primary target of antidepressant discovery research for several decades. Since the currently available antidepressants may cause adverse effects and have several limitations, novel drugs are highly desired. However, the efforts to develop better therapeutics are hampered by the lack of a crystal structure of hSERT. Knowledge of the binding site of the drug and its orientation in the protein is crucial in structure-based drug discovery. We employed a novel computational protocol comprised of active site detection, docking, scoring, molecular dynamics simulations, and absolute binding free energy (ABFE) calculations to elucidate the binding site and the binding mode of a dual hSERT/5HT-1A blocker SSA-426 and our in-house hSERT inhibitor DJLDU-3-79 in hSERT. Through this approach, we propose that both of these inhibitors bind in the S1 pocket of hSERT and in a similar orientation. This disproves the earlier hypothesis that both these inhibitors bind in the S2 site; however, we are in agreement with the earlier hypothesis that both of the ligands orient similarly. Further, we resolved the ambiguity in binding energies and binding trends of the tricyclic antidepressant drugs clomipramine, imipramine, and desipramine with leucine transporter (LeuT) (a bacterial homologue of hSERT) through relative binding free energy (RBFE) calculations. Based on our RBFE results, we proposed that clomipramine should have the highest affinity for LeuT, followed by imipramine and desipramine. Finally, to achieve accuracy in binding energy estimations and to perform all CHARMM simulations, we developed CHARMM general force field parameters (CGenFF) for fifteen monoamine transporter ligands