Disorder-to-order active site capping regulates the rate-limiting step of the inositol pathway.

Myo-inositol-1-phosphate synthase (MIPS) catalyzes the NAD+-dependent isomerization of glucose-6-phosphate (G6P) into inositol-1-phosphate (IMP), controlling the rate-limiting step of the inositol pathway. Previous structural studies focused on the detailed molecular mechanism, neglecting large-scale conformational changes that drive the function of this 240 kDa homotetrameric complex. In this study, we identified the active, endogenous MIPS in cell extracts from the thermophilic fungus Thermochaetoides thermophila. By resolving the native structure at 2.48 Å (FSC = 0.143), we revealed a fully populated active site. Utilizing 3D variability analysis, we uncovered conformational states of MIPS, enabling us to directly visualize an order-to-disorder transition at its catalytic center. An acyclic intermediate of G6P occupied the active site in two out of the three conformational states, indicating a catalytic mechanism where electrostatic stabilization of high-energy intermediates plays a crucial role. Examination of all isomerases with known structures revealed similar fluctuations in secondary structure within their active sites. Based on these findings, we established a conformational selection model that governs substrate binding and eventually inositol availability. In particular, the ground state of MIPS demonstrates structural configurations regardless of substrate binding, a pattern observed across various isomerases. These findings contribute to the understanding of MIPS structure-based function, serving as a template for future studies targeting regulation and potential therapeutic applications.

Variation of Cyclodextrin (CD) Complexation with Biogenic Amine Tyramine: Pseudopolymorphs of ß-CD Inclusion vs. α-CD Exclusion, Deep Atomistic Insights.

Tyramine (TRM) is a biogenic catecholamine neurotransmitter, which can trigger migraines and hypertension. TRM accumulated in foods is reduced and detected using additive cyclodextrins (CDs) while their association characteristics remain unclear. Here, single-crystal X-ray diffraction and density functional theory (DFT) calculation have been performed, demonstrating the elusive pseudopolymorphs in ß-CD inclusion complexes with TRM base/HCl, ß-CD·0.5TRM·7.6H2O (1) and ß-CD·TRM HCl·4H2O (2) and the rare α-CD·0.5(TRM HCl)·10H2O (3) exclusion complex. Both 1 and 2 share the common inclusion mode with similar TRM structures in the round and elliptical ß-CD cavities, belong to the monoclinic space group P21, and have similar herringbone packing structures. Furthermore, 3 differs from 2, as the smaller twofold symmetry-related, round α-CD prefers an exclusion complex with the twofold disordered TRM-H+ sites. In the orthorhombic P21212 lattice, α-CDs are packed in a channel-type structure, where the column-like cavity is occupied by disordered water sites. DFT results indicate that ß-CD remains elliptical to suitably accommodate TRM, yielding an energetically favorable inclusion complex, which is significantly contributed by the ß-CD deformation, and the inclusion complex of α-CD with the TRM aminoethyl side chain is also energetically favorable compared to the exclusion mode. This study suggests the CD implications for food safety and drug/bioactive formulation and delivery.

The Nitrofuran-Warhead-Equipped Spirocyclic Azetidines Show Excellent Activity against Mycobacterium tuberculosis.

A series of 21 new 7'H-spiro[azetidine-3,5'-furo [3,4-d]pyrimidine]s substituted at the pyrimidine ring second position were synthesized. The compounds showed high antibacterial in vitro activity against M. tuberculosis. Two compounds had lower minimum inhibitory concentrations against Mtb (H37Rv strain) compared with isoniazid. The novel spirocyclic scaffold shows excellent properties for anti-tuberculosis drug development.

Binding Affinity Determination in Drug Design: Insights from Lock and Key, Induced Fit, Conformational Selection, and Inhibitor Trapping Models.

Binding affinity is a fundamental parameter in drug design, describing the strength of the interaction between a molecule and its target protein. Accurately predicting binding affinity is crucial for the rapid development of novel therapeutics, the prioritization of promising candidates, and the optimization of their properties through rational design strategies. Binding affinity is determined by the mechanism of recognition between proteins and ligands. Various models, including the lock and key, induced fit, and conformational selection, have been proposed to explain this recognition process. However, current computational strategies to predict binding affinity, which are based on these models, have yet to produce satisfactory results. This article explores the connection between binding affinity and these protein-ligand interaction models, highlighting that they offer an incomplete picture of the mechanism governing binding affinity. Specifically, current models primarily center on the binding of the ligand and do not address its dissociation. In this context, the concept of ligand trapping is introduced, which models the mechanisms of dissociation. When combined with the current models, this concept can provide a unified theoretical framework that may allow for the accurate determination of the ligands' binding affinity.

Amlodipine Ocular Delivery Restores Ferning Patterns and Reduces Intensity of Glycosylated Peak of Carrageenan-Induced Tear Fluid: An In-Silico Flexible Docking with IL-ß1.

BACKGROUND: The tear ferning test can be an easy clinical procedure for the evaluation and characterization of the ocular tear film. OBJECTIVE: The objective of this study was to examine the restoration of tear ferning patterns and reduction of glycosylation peak after amlodipine application in carrageenan-induced conjunctivitis. METHODS: At the rabbit's upper palpebral region, carrageenan was injected for cytokine-mediated conjunctivitis. Ferning pattern and glycosylation of the tear fluid were characterized using various instrumental analyses. The effect of amlodipine was also examined after ocular instillation and flexible docking studies. RESULTS: Optical microscopy showed a disrupted ferning of the tear collected from the inflamed eye. FTIR of the induced tear fluid exhibited peaks within 1000-1200 cm-1, which might be due to the protein glycosylation absent in the normal tear spectrogram. The glycosylation peak reduced significantly in the tear sample collected from the amlodipine-treated group. Corresponding energy dispersive analysis showed the presence of sulphur, indicating protein leakage from the lacrimal gland in the induced group. The disappearance of sulphur from the treated group indicated its remedial effect. The flexible docking studies revealed a stronger binding mode of amlodipine with Interleukin-1ß (IL-1ß). The reduction in the intensity of the glycosylated peak and the restoration offering are probably due to suppression of IL-1ß. CONCLUSION: This study may be helpful in obtaining primary information for drug discovery to be effective against IL-1ß and proving tear fluid as a novel diagnostic biomarker.

Alternate recognition by dengue protease: Proteolytic and binding assays provide functional evidence beyond an induced-fit.

Proteases are key enzymes in viral replication, and interfering with these targets is the basis for therapeutic interventions. We previously introduced a hypothesis about conformational selection in the protease of dengue virus and related flaviviruses, based on conformational plasticity noted in X-ray structures. The present work presents the first functional evidence for alternate recognition by the dengue protease, in a mechanism based primarily on conformational selection rather than induced-fit. Recognition of distinct substrates and inhibitors in proteolytic and binding assays varies to a different extent, depending on factors reported to influence the protease structure. The pH, salinity, buffer type, and temperature cause a change in binding, proteolysis, or inhibition behavior. Using representative inhibitors with distinct structural scaffolds, we identify two contrasting binding profiles to dengue protease. Noticeable effects are observed in the binding assay upon inclusion of a non-ionic detergent in comparison to the proteolytic assay. The findings highlight the impact of the selection of testing conditions on the observed ligand affinity or inhibitory potency. From a broader scope, the dengue protease presents an example, where the induced-fit paradigm appears insufficient to explain binding events with the biological target. Furthermore, this protein reveals the complexity of comparing or combining biochemical assay data obtained under different conditions. This can be particularly critical for artificial intelligence (AI) approaches in drug discovery that rely on large datasets of compounds activity, compiled from different sources using non-identical testing procedures. In such cases, mismatched results will compromise the model quality and its predictive power.

Michaelis-like complex of mouse ketohexokinase isoform C.

Over the past forty years there has been a drastic increase in fructose-related diseases, including obesity, heart disease and diabetes. Ketohexokinase (KHK), the first enzyme in the liver fructolysis pathway, catalyzes the ATP-dependent phosphorylation of fructose to fructose 1-phosphate. Understanding the role of KHK in disease-related processes is crucial for the management and prevention of this growing epidemic. Molecular insight into the structure-function relationship in ligand binding and catalysis by KHK is needed for the design of therapeutic inhibitory ligands. Ketohexokinase has two isoforms: ketohexokinase A (KHK-A) is produced ubiquitously at low levels, whereas ketohexokinase C (KHK-C) is found at much higher levels, specifically in the liver, kidneys and intestines. Structures of the unliganded and liganded human isoforms KHK-A and KHK-C are known, as well as structures of unliganded and inhibitor-bound mouse KHK-C (mKHK-C), which shares 90% sequence identity with human KHK-C. Here, a high-resolution X-ray crystal structure of mKHK-C refined to 1.79 Šresolution is presented. The structure was determined in a complex with both the substrate fructose and the product of catalysis, ADP, providing a view of the Michaelis-like complex of the mouse ortholog. Comparison to unliganded structures suggests that KHK undergoes a conformational change upon binding of substrates that places the enzyme in a catalytically competent form in which the ß-sheet domain from one subunit rotates by 16.2°, acting as a lid for the opposing active site. Similar kinetic parameters were calculated for the mouse and human enzymes and indicate that mice may be a suitable animal model for the study of fructose-related diseases. Knowledge of the similarity between the mouse and human enzymes is important for understanding preclinical efforts towards targeting this enzyme, and this ground-state, Michaelis-like complex suggests that a conformational change plays a role in the catalytic function of KHK-C.

Crystal structure of NRAS Q61K with a ligand-induced pocket near switch II.

The RAS isoforms (KRAS, HRAS and NRAS) have distinct cancer type-specific profiles. NRAS mutations are the second most prevalent RAS mutations in skin and hematological malignancies. Although RAS proteins were considered undruggable for decades, isoform and mutation-specific investigations have produced successful RAS inhibitors that are either specific to certain mutants, isoforms (pan-KRAS) or target all RAS proteins (pan-RAS). While extensive structural and biochemical investigations have focused mainly on K- and H-RAS mutations, NRAS mutations have received less attention, and the most prevalent NRAS mutations in human cancers, Q61K and Q61R, are rare in K- and H-RAS. This manuscript presents a crystal structure of the NRAS Q61K mutant in the GTP-bound form. Our structure reveals a previously unseen pocket near switch II induced by the binding of a ligand to the active form of the protein. This observation reveals a binding site that can potentially be exploited for development of inhibitors against mutant NRAS. Furthermore, the well-resolved catalytic site of this GTPase bound to native GTP provides insight into the stalled GTP hydrolysis observed for NRAS-Q61K.

Time-series analysis of rhenium(I) organometallic covalent binding to a model protein for drug development.

Metal-based complexes with their unique chemical properties, including multiple oxidation states, radio-nuclear capabilities and various coordination geometries yield value as potential pharmaceuticals. Understanding the interactions between metals and biological systems will prove key for site-specific coordination of new metal-based lead compounds. This study merges the concepts of target coordination with fragment-based drug methodologies, supported by varying the anomalous scattering of rhenium along with infrared spectroscopy, and has identified rhenium metal sites bound covalently with two amino acid types within the model protein. A time-based series of lysozyme-rhenium-imidazole (HEWL-Re-Imi) crystals was analysed systematically over a span of 38 weeks. The main rhenium covalent coordination is observed at His15, Asp101 and Asp119. Weak (i.e. noncovalent) interactions are observed at other aspartic, asparagine, proline, tyrosine and tryptophan side chains. Detailed bond distance comparisons, including precision estimates, are reported, utilizing the diffraction precision index supplemented with small-molecule data from the Cambridge Structural Database. Key findings include changes in the protein structure induced at the rhenium metal binding site, not observed in similar metal-free structures. The binding sites are typically found along the solvent-channel-accessible protein surface. The three primary covalent metal binding sites are consistent throughout the time series, whereas binding to neighbouring amino acid residues changes through the time series. Co-crystallization was used, consistently yielding crystals four days after setup. After crystal formation, soaking of the compound into the crystal over 38 weeks is continued and explains these structural adjustments. It is the covalent bond stability at the three sites, their proximity to the solvent channel and the movement of residues to accommodate the metal that are important, and may prove useful for future radiopharmaceutical development including target modification.

Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch.

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.

Structural Differences at Quadruplex-Duplex Interfaces Enable Ligand-Induced Topological Transitions.

Quadruplex-duplex (QD) junctions, which represent unique structural motifs of both biological and technological significance, have been shown to constitute high-affinity binding sites for various ligands. A QD hybrid construct based on a human telomeric sequence, which harbors a duplex stem-loop in place of a short lateral loop, is structurally characterized by NMR. It folds into two major species with a (3+1) hybrid and a chair-type (2+2) antiparallel quadruplex domain coexisting in a K+ buffer solution. The antiparallel species is stabilized by an unusual capping structure involving a thymine and protonated adenine base AH+ of the lateral loop facing the hairpin duplex to form a T·AH+·G·C quartet with the interfacial G·C base pair at neutral pH. Addition and binding of Phen-DC3 to the QD hybrid mixture by its partial intercalation at corresponding QD junctions leads to a topological transition with exclusive formation of the (3+1) hybrid fold. In agreement with the available experimental data, such an unprecedented discrimination of QD junctions by a ligand can be rationalized following an induced fit mechanism.

Optimizing Iodine Enrichment through Induced-Fit Transformations in a Flexible Ag(I)-Organic Framework: From Accelerated Adsorption Kinetics to Record-High Storage Density.

Efficient capture and storage of radioactive I2 is a prerequisite for developing nuclear power but remains a challenge. Here, two flexible Ag-MOFs (FJI-H39 and 40) with similar active sites but different pore sizes and flexibility are prepared; both of them can capture I2 with excellent removal efficiencies and high adsorption capacities. Due to the more flexible pores, FJI-H39 not only possesses the record-high I2 storage density among all the reported MOFs but also displays a very fast adsorption kinetic (124 times faster than FJI-H40), while their desorption kinetics are comparable. Mechanistic studies show that FJI-H39 can undergo induced-fit transformations continuously (first contraction then expansion), making the adsorbed iodine species enrich near the Ag(I) nodes quickly and orderly, from discrete I- anion to the dense packing of various iodine species, achieving the very fast adsorption kinetic and the record-high storage density simultaneously. However, no significant structural transformations caused by the adsorbed iodine are observed in FJI-H40. In addition, FJI-H39 has excellent stability/recyclability/obtainability, making it a practical adsorbent for radioactive I2. This work provides a useful method for synthesizing practical radioactive I2 adsorbents.

Multi-Channel smFRET study reveals a Compact conformation of EF-G on the Ribosome.

While elongation factor G (EF-G) is crucial for ribosome translocation, the role of its GTP hydrolysis remains ambiguous. EF-G's indispensability is further exemplified by the phosphorylation of human eukaryotic elongation factor 2 (eEF2) at Thr56, which inhibits protein synthesis globally, but its exact mechanism is not clear. In this study, we developed a multi-channel single-molecule FRET (smFRET) microscopy methodology to examine the conformational changes of E. coli EF-G induced by mutations that closely aligned with eEF2's Thr56 residue. We utilized Alexa 488/594 double-labeled EF-G to catalyze the translocation of fMet-Phe-tRNAPhe-Cy3 inside Cy5-L27 labeled ribosomes, allowing us to probe both processes within the same complex. Our findings indicate that in the presence of either GTP or GDPCP, wild-type EF-G undergoes a conformational extension upon binding to the ribosome to promote normal translocation. On the other hand, T48E and T48V mutations did not affect GTP/GDP binding or GTP hydrolysis, but impeded Poly(Phe) synthesis and caused EF-G to adopt a unique compact conformation, which wasn't observed when the mutants interact solely with the sarcin/ricin loop. This study provides new insights into EF-G's adaptability and sheds light on the modification mechanism of human eEF2.

Structural basis for the allosteric behaviour and substrate specificity of Lactococcus lactis Prolidase.

Prolidase (EC 3.4.13.9) is an enzyme that specifically hydrolyzes Xaa-Pro dipeptides into free amino acids. We previously studied kinetic behaviours and solved the crystal structure of wild-type (WT) Lactococcus lactis prolidase (Llprol), showing that this homodimeric enzyme has unique characteristics: allosteric behaviour and substrate inhibition. In this study, we focused on solving the crystal structures of three Llprol mutants (D36S, H38S, and R293S) which behave differently in v-S plots. The D36S and R293S Llprol mutants do not show allosteric behaviour, and the Llprol mutant H38S has allosteric behaviour comparable to the WT enzyme (Hill constant 1.52 and 1.58, respectively). The crystal structures of Llprol variants suggest that the active site of Llprol formed with amino acid residues from both monomers, i.e., located in an interfacial area of dimer. The comparison between the structure models of Llprol indicated that the two monomers in the dimers of Llprol variants have different relative positions among Llprol variants. They showed different interatomic distances between the amino acid residues bridging the two monomers and varied sizes of the solvent-accessible interface areas in each Llprol variant. These observations indicated that Llprol could adapt to different conformational states with distinctive substrate affinities. It is strongly speculated that the domain movements required for productive substrate binding are restrained in allosteric Llprol (WT and H38S). At low substrate concentrations, only one out of the two active sites at the dimer interface could accept substrate; as a result, the asymmetrical activated dimer leads to allosteric behaviour.

"Covalent-Disassembly"-Based Approaches For Sensing Applications.

The detection of analytes with small molecular probes is crucial for the analysis and understanding of chemical, medicinal, environmental and biological situations as well as processes. Classic detection approaches rely on the concept of molecular recognition and bond formation reactions. Bond breakage reactions have been less explored in similar contexts. This concept article introduces metal-salen and metal-imine complexes as "covalent-disassembly"-based (DB)-probes for detecting polyoxophosphates, thiols, amino acids, HCN and changes in pH. It discusses the roles, importance and combinations of structurally functionalized molecular building blocks in the construction of DB-probes. Applications of optimized DB-probes for analyte detection in live cells and foodstuff are also discussed. Furthermore, the mechanism of the disassembly of a Fe(III)-salen probe upon pyrophosphate binding is presented. Extraordinary selectivity for this analyte was achieved by a multistep disassembly sequence including an unprecedented structural change of the metal complex (i. e. "induced-fit" principle). Design principles of probes for sensing applications following the "covalent-disassembly" approach are summarized, which will help improving current systems, but will also facilitate the development of new DB-probes for challenging analytic targets.

Discovery of new thymol-3,4-disubstituted thiazole hybrids as dual COX-2/5-LOX inhibitors with in vivo proof.

New thymol-3,4-disubstitutedthiazole hybrids were synthesised as dual COX-2/5-LOX inhibitors. Compounds 6b, 6d, 6e, and 6f displayed in vitro inhibitory activity against COX-2 (IC50= 0.037, 0.042, 0.046, and 0.039 µM) nearly equal to celecoxib (IC50= 0.045 µM). 6b, 6d, and 6f showed SI (379, 341, and 374, respectively) higher than that of celecoxib (327). 6a-l elicited in vitro 5-LOX inhibitory activity higher than quercetin. 6a-f, 6i-l, 7a, and 7c possessed in vivo inhibition of formalin induced paw edoema higher than celecoxib. 6a, 6b, 6f, 6h-l, and 7b showed gastrointestinal safety profile as celecoxib and diclofenac sodium in the population of fasted rats. Induced fit docking and molecular dynamics simulation predicted good fitting of 6b and 6f without changing the packing and globularity of the apo protein. In conclusion, 6b and 6f achieved the target goal as multitarget inhibitors of inflammation.

Dynamics and interactions of intrinsically disordered proteins.

Intrinsically disordered proteins (IDPs) are widespread in eukaryotes and participate in a variety of important cellular processes. Numerous studies using state-of-the-art experimental and theoretical methods have advanced our understanding of IDPs and revealed that disordered regions engage in a large repertoire of intra- and intermolecular interactions through their conformational dynamics, thereby regulating many intracellular functions in concert with folded domains. The mechanisms by which IDPs interact with their partners are diverse, depending on their conformational propensities, and include induced fit, conformational selection, and their mixtures. In addition, IDPs are implicated in many diseases, and progress has been made in designing inhibitors of IDP-mediated interactions. Here we review these recent advances with a focus on the dynamics and interactions of IDPs.

Unveiling the therapeutic potential of a mutated paraoxonase 2 in diabetic retinopathy: Defying glycation, mitigating oxidative stress, ER stress and inflammation.

Paraoxonase 2 (PON2) is an intracellular anti-oxidant protein ubiquitously expressed in all cells and reduces reactive oxygen species, endoplasmic reticulum (ER) stress, further improves mitochondrial function and thereby shows anti-apoptotic function. In diabetes and its complications this PON gets glycated and becomes in effective. The PON activity is reported to be reduced in diabetic retinopathy and we have earlier showed Carboxy methyl lysine (AGE) decreased PON2 expression and activity in Human retinal endothelial cells (HREC) . In this study, we have designed and developed a mutated PON2 by in silico and in vitro approach which can resist glycation. Where in glycation-prone residues in PON2 was predicted using in silico analyses and a mutated PON2 was developed using in vitro site directed mutagenesis (SDM) assay mPON2 (mutant PON2-PON2-K70A) and its efficacy was compared with wPON2 (wild type PON2). CML glycated wPON2 and reduced its activity when compared with mPON2 in HREC confirmed by immunoprecipitation and in vitro experiments. Additionally, mPON2 interaction efficiency with its substrates was higher than wPON2 by insilico assay and demonstrated enhanced inhibition against CML-induced oxidative stress, ER stress, pro-inflammation, and mitochondrial fission than wPON2 by invitro assay. Further mPON2 showed increased inhibition of phosphorylation of NFĸB induced by CML. Our investigation establishes that the over expression of mPON2 in HREC can defy glycation and therefore mitigate ER stress and inflammation against CML than endogenous wPON2. These findings imply that mPON2 can be a beneficial therapeutic target against diabetic retinopathy.

A Kinetic Transition Network Model Reveals the Diversity of Protein Dimer Formation Mechanisms.

Protein homodimers have been classified as three-state or two-state dimers depending on whether a folded monomer forms before association, but the details of the folding-binding mechanisms are poorly understood. Kinetic transition networks of conformational states have provided insight into the folding mechanisms of monomeric proteins, but extending such a network to two protein chains is challenging as all the relative positions and orientations of the chains need to be included, greatly increasing the number of degrees of freedom. Here, we present a simplification of the problem by grouping all states of the two chains into two layers: a dissociated and an associated layer. We combined our two-layer approach with the Wako-Saito-Muñoz-Eaton method and used Transition Path Theory to investigate the dimer formation kinetics of eight homodimers. The analysis reveals a remarkable diversity of dimer formation mechanisms. Induced folding, conformational selection, and rigid docking are often simultaneously at work, and their contribution depends on the protein concentration. Pre-folded structural elements are always present at the moment of association, and asymmetric binding mechanisms are common. Our two-layer network approach can be combined with various methods that generate discrete states, yielding new insights into the kinetics and pathways of flexible binding processes.

Benzimidazoles Containing Piperazine Skeleton at C-2 Position as Promising Tubulin Modulators with Anthelmintic and Antineoplastic Activity.

Benzimidazole anthelmintic drugs hold promise for repurposing as cancer treatments due to their interference with tubulin polymerization and depolymerization, manifesting anticancer properties. We explored the potential of benzimidazole compounds with a piperazine fragment at C-2 as tubulin-targeting agents. In particular, we assessed their anthelmintic activity against isolated Trichinella spiralis muscle larvae and their effects on glioblastoma (U-87 MG) and breast cancer (MDA-MB-231) cell lines. Compound 7c demonstrated exceptional anthelmintic efficacy, achieving a 92.7% reduction in parasite activity at 100 µg/mL after 48 hours. In vitro cytotoxicity analysis of MDA-MB 231 and U87 MG cell lines showed that derivatives 7b, 7d, and 7c displayed lower IC50 values compared to albendazole (ABZ), the control. These piperazine benzimidazoles effectively reduced cell migration in both cell lines, with compound 7c exhibiting the most significant reduction, making it a promising candidate for further study. The binding mode of the most promising compound 7c, was determined using the induced fit docking-molecular dynamics (IFD-MD) approach. Regular docking and IFD were also employed for comparison. The IFD-MD analysis revealed that 7c binds to tubulin in a unique binding cavity near that of ABZ, but the benzimidazole ring was fitted much deeper into the binding pocket. Finally, the absolute free energy of perturbation technique was applied to evaluate the 7c binding affinity, further confirming the observed binding mode.