A computational analysis of the binding affinities of FKBP12 inhibitors using the MM-PB/SA method

The FK506-binding proteins have been targets of pharmaceutical interests over years. We have studied the binding of a set of 12 nonimmunosuppressive small-molecule inhibitors to FKBP12 through molecular dynamics simulations. Each complex was subjected to 1-ns MD simulation conducted in an explicit solvent environment under constant temperature and pressure. The binding free energy of each complex was then computed by the MM-PB/SA method in the AMBER program. Our MM-PB/SA computation produced a good correlation between the experimentally determined and the computed binding free energies with a correlation coefficient (R(2)) of 0.93 and a standard deviation as low as 0.30 kcal/mol. The vibrational entropy term given by the normal mode analysis was found to be helpful for achieving this correlation. Moreover, an adjustment to one weight factor in the PB/SA model was essential to correct the absolute values of the final binding free energies to a reasonable range. A head-to-head comparison of our MM-PB/SA model with a previously reported Linear Response Approximation (LRA) model suggested that the MM-PB/SA method is more robust in binding affinity prediction for this class of compounds.     

The FKBP12 subunit modifies the long-range allosterism of the ryanodine receptor

Ryanodine receptors (RyRs) are large conductance intracellular channels controlling intracellular calcium homeostasis in myocytes, neurons, and other cell types. Loss of RyR’s constitutive cytoplasmic partner FKBP results in channel sensitization, dominant subconductance states, and increased cytoplasmic Ca²⁺. FKBP12 binds to RyR1’s cytoplasmic assembly 130?Å away from the ion gate at four equivalent sites in the RyR1 tetramer. To understand how FKBP12 binding alters RyR1’s channel properties, we studied the 3D structure of RyR1 alone in the closed conformation in the context of the open and closed conformations of FKBP12-bound RyR1. We analyzed the metrics of conformational changes of existing structures, the structure of the ion gate, and carried out multivariate statistical analysis of thousands of individual cryoEM RyR1 particles. We find that under closed state conditions, in the presence of FKBP12, the cytoplasmic domain of RyR1 adopts an upward conformation, whereas absence of FKBP12 results in a relaxed conformation, while the ion gate remains closed. The relaxed conformation is intermediate between the RyR1-FKBP12 complex closed (upward) and open (downward) conformations. The closed-relaxed conformation of RyR1 appears to be consistent with a lower energy barrier separating the closed and open states of RyR1-FKBP12, and suggests that FKBP12 plays an important role by restricting conformations within RyR1’s conformational landscape.     

Docking and molecular dynamics studies of new potential inhibitors of the human epidermal receptor 2

Compounds similar to lapatinib and gefitinib have been investigated as potential inhibitors of the intracellular receptor tyrosine kinase (RTK) domain of the human epidermal receptor 2 (HER2), which is a promising molecular target to the drug design of new chemotherapies for breast, lung, ovarian and colorectal cancers. In this study, we have searched potential HER2 inhibitors used for treatment of other illnesses such as hepatitis, bacterial infections and sexual impotence screened in the DrugBank. The compounds selected were subjected to virtual screening docking in order to evaluate the main interactions between them and the RTK domain of HER2. The selected compounds were investigated by flexible docking, molecular dynamics studies and ΔG _bind calculations. The results suggest that antrafenine, saprisartan, reserpine, irinotecan and udenafil are potential candidates to inhibit the RTK domain of HER2.     

Insertion of a chaperone domain converts FKBP12 into a powerful catalyst of protein folding

The catalytic activity of human FKBP12 as a prolyl isomerase is high towards short peptides, but very low in proline-limited protein folding reactions. In contrast, the SlyD proteins, which are members of the FKBP family, are highly active as folding enzymes. They contain an extra "insert-in-flap" or IF domain near the prolyl isomerase active site. The excision of this domain did not affect the prolyl isomerase activity of SlyD from Escherichia coli towards short peptide substrates but abolished its catalytic activity in proline-limited protein folding reactions. The reciprocal insertion of the IF domain of SlyD into human FKBP12 increased its folding activity 200-fold and generated a folding catalyst that is more active than SlyD itself. The IF domain binds to refolding protein chains and thus functions as a chaperone module. A prolyl isomerase catalytic site and a separate chaperone site with an adapted affinity for refolding protein chains are the key elements for a productive coupling between the catalysis of prolyl isomerization and conformational folding in the enzymatic mechanisms of SlyD and other prolyl isomerases, such as trigger factor and FkpA.     

Structural insights into the inhibitor binding and new inhibitor design to Polo-like kinase-1 Polo-box domain using computational studies

Polo box domain (PBD) from Polo-Like Kinase-1 (PLK-1) a cell cycle regulator is one of the important non-kinase targets implicated in various cancers. The crystal structure of PLK-1 PBD bound to phosphopeptide inhibitor is available and acylthiourea derivatives have been reported as potent PBD inhibitors. In this work, structure and ligand-based pharmacophore methods have been used to identify new PBD inhibitors. The binding of acylthiourea analogs and new inhibitors to PBD were assessed using molecular docking and molecular dynamics simulations to understand their binding interactions, investigate the complex stability and reveal the molecular basis for inhibition. This study provides the binding free energies and residue-wise contributions to decipher the essential interactions in the protein-inhibitor complementarity for complex formation and the design of new PBD inhibitors with better binding.

Communicated by Ramaswamy H. Sarma     

Identification of potential inhibitors of Fasciola gigantica thioredoxin1: computational screening,molecular dynamics simulation,and binding free energy studies

Fasciola gigantica is the causative organism of fascioliasis and is responsible for major economic losses in livestock production globally. F. gigantica thioredoxin1 (FgTrx1) is an important redox-active enzyme involved in maintaining the redox homeostasis in the cell. To identify a potential anti-fasciolid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,740) against the FgTrx1 structure. The ligands were docked against FgTrx1 and 309 ligands were found to have better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 30 compounds were found to fit well for re-docking studies. After refinement by molecular docking and drug-likeness analysis, three potential inhibitors (ZINC15970091, ZINC9312362, and ZINC9312661) were identified. These three ligands were further subjected to molecular dynamics simulation (MDS) to compare the dynamics and stability of the protein structure after binding of the ligands. The binding free energy analyses were calculated to determine the intermolecular interactions. The results suggested that the two compounds had a binding free energy of –82.237, and –109.52 kJ.mol^?1 for compounds with IDs ZINC9312362 and ZINC9312661, respectively. These predicted compounds displayed considerable pharmacological and structural properties to be drug candidates. We concluded that these two compounds could be potential drug candidates to fight against F. gigantica parasites.     

Molecular modeling,dynamics, and an insight into the structural inhibition of cofactor independent phosphoglycerate mutase isoform 1 from Wuchereria bancrofti using cheminformatics and mutational studies

Phosphoglycerate mutase catalyzes the interconversion between 2-phosphoglycerate and 3-phosphoglycerate in the glycolytic and gluconeogenic pathways. They exist in two unrelated forms, that is either cofactor (2,3-diphosphoglycerate) dependent or cofactor-independent. These two enzymes have no similarity in amino acid sequence, tertiary structure, and in catalytic mechanism. Wuchereria bancrofti (WB) contains the cofactor-independent form, whereas other organisms can possess the dependent form or both. Since, independent phosphoglycerate mutase (iPGM) is an essential gene for the survival of nematodes, and it has no sequence or structural similarity to the cofactor-dependent phosphoglycerate mutase found in mammals, it represents an attractive drug target for the filarial nematodes. In this current study, a putative cofactor-iPGM gene was identified in the protein sequence of the WB. In the absence of crystal structure, a three-dimensional structure was determined using the homology modeling approximation, and the most stable protein conformation was identified through the molecular dynamics simulation studies, using GROMACS 4.5. Further, the functional or characteristic residues were identified through the sequence analysis, potential inhibitors were short-listed and validated, and potential inhibitors were ranked using the cheminformatics and molecular dynamics simulations studies, Prime MM-GBSA approach, respectively.