Computational screening of potential inhibitors targeting MurF of Brugia malayi Wolbachia through multi-scale molecular docking, molecular dynamics and MM-GBSA analysis

Lymphatic filariasis is a parasitic disease caused by the worms Wuchereria bancrofti, Brugia malayi and Brugia timori. Three anti-filarial drugs namely Diethylcarbamazine, Ivermectin and Albendazole and their combinations are used as the control strategies for filariasis. The disease has received much attention in drug discovery due to the unavailability of vaccines and the toxic pharmaceutical properties of the existing drugs. In Wolbachia endosymbiont Brugia malayi, the UDP-N-acetylmuramoyl-tripeptide-d-alanyl-d-alanine ligase (MurF) plays a key role in peptidoglycan biosynthesis pathway and therefore can be considered as effective drug target against filariasis disease. Therefore, in the present study, MurF was selected as the therapeutic target to identify specific inhibitors against filariasis. Homology modeling was performed to predict the three-dimensional structure of MurF due to the absence of the experimental structure. Further molecular dynamics simulation and structure-based high throughput virtual screening with three different chemical databases (Zinc, Maybridge and Specs) were carried out to identify potent inhibitors and also to check their conformations inside the binding site of MurF, respectively. Top three compounds with high docking score and high relative binding affinity against MurF were selected. Further, validation studies, including predicted ADME (Absorption, Distribution, Metabolism, Excretion) assessment, binding free energy using MM-GBSA (Molecular Mechanics Generalized Born Surface Area) and DFT (Density Functional Theory) calculations were performed for the top three compounds. From the results, it was observed that all the three compounds were predicted to show high reactivity, acceptable range of pharmacokinetic properties and high binding affinity with the drug target MurF. Overall, the results could provide more understanding on the inhibition of MurF enzyme and the screened compounds could lead to the development of new specific anti-filarial drugs.

Water Mapping and Scoring approaches to predict the role of Hydration sites in Binding Affinity of PAK1 inhibitors

AIM

This study aims to develop and establish a computational model that can identify potent molecules for p21-activating kinase 1 (PAK1).

BACKGROUND

PAK1 is a well-established drug target that has been explored for various therapeutic interventions. Control of this protein requires an indispensable inhibitor to curb the structural changes and subsequent activation of signalling effectors responsible for the progression of diseases, such as cancer, inflammatory, viral, and neurological disorders.

OBJECTIVE

To establish a computational model that could identify active molecules which will further provide a platform for developing potential PAK1 inhibitors.

METHOD

A congeneric series of 27 compounds was considered for this study with Ki (nm) covering a minimum of 3 log range. The compounds were developed based on a previously reported Group-I PAK inhibitor, namely G-5555. The 27 compounds were subjected to the SP and XP mode of docking, to understand the binding mode, its conformation and interaction patterns. To understand the relevance of biological activity from computational approaches, the compounds were scored against generated water maps to obtain WM/MM ΔG binding energy. Moreover, molecular dynamics analysis was performed for the highly active compound, to understand the conformational variability and complex's stability. We then evaluate the predictable binding pose obtained from the docking studies.

RESULT

From the SP and XP modes of docking, the common interaction pattern with the amino acid residues Arg299 (cation-π), Glu345 (Aromatic hydrogen bond), hinge region Leu347, salt bridges Asp393 and Asp407 was observed, among the congeneric compounds. The interaction pattern was compared with the co-crystal inhibitor FRAX597 of the PAK1 crystal structure (PDB id: 4EQC). The correlation with different docking parameters in the SP and XP modes was insignificant and thereby revealed that the SP and XP's scoring functions could not predict the active compounds. This was due to the limitations in the docking methodology that neglected the receptor flexibility and desolvation parameters. Hence, to recognise the desolvation and explicit solvent effects, as well as to study the Structure-Activity Relationships (SARs) extensively, WaterMap (WM) calculations were performed on the congeneric compounds. Based on displaceable unfavourable hydration sites (HS) and their associated thermodynamic properties, the WM calculations facilitated to understand the significance of correlation in the folds of activity of highly (19 and 17), moderate (16 and 21) and less active (26 and 25) compounds. Furthermore, the scoring function from WaterMap, namely WM/MM, led to a significant R2 value of 0.72, due to a coupled conjunction with MM treatment and displaced unfavourable waters at the binding site. To check the "optimal binding conformation", molecular dynamics simulation was carried out with the highly active compound 19 to explain the binding mode, stability, interactions, solvent accessible area, etc., which could support the predicted conformation with bioactive conformation.

CONCLUSION

This study determined the best scoring function, established SARs and predicted active molecules through a computational model. This will contribute towards development of the most potent PAK1 inhibitors.

Identification of Pak1 inhibitors using water thermodynamic analysis

p21-activated kinases (Paks) play an integral component in various cellular diverse processes. The full activation of Pak is dependent upon several serine residues present in the N-terminal region, a threonine present at the activation loop, and finally the phosphorylation of these residues ensure the complete activation of Pak1. The present study deals with the identification of novel potent candidates of Pak1 using computational methods as anti-cancer compounds. A diverse energy based pharmacophore (e-pharmacophore) was developed using four co-crystal inhibitors of Pak1 having pharmacophore features of 5 (DRDRR), 6 (DRHADR), and 7 (RRARDRP and DRRDADH) hypotheses. These models were used for rigorous screening against e-molecule database. The obtained hits were filtered using ADME/T and molecular docking to identify the high affinity binders. These hits were subjected to hierarchical clustering using dendritic fingerprint inorder to identify structurally diverse molecules. The diverse hits were scored against generated water maps to obtain WM/MM ΔG binding energy. Furthermore, molecular dynamics simulation and density functional theory calculations were performed on the final hits to understand the stability of the complexes. Five structurally diverse novel Pak1 inhibitors (4835785, 32198676, 32407813, 76038049, and 32945545) were obtained from virtual screening, water thermodynamics and WM/MM ΔG binding energy. All hits revealed similar mode of binding pattern with the hinge region residues replacing the unstable water molecules in the binding site. The obtained novel hits could be used as a platform to design potent drugs that could be experimentally tested against cancer patients having increased Pak1 expression.