Molecular dynamics study on inhibition mechanism of CDK-2 and GSK-3β by CHEMBL272026 molecule

Identifying potential Alzheimer's disease therapeutics through GSK-3β inhibition: A molecular docking and dynamics approach

Emerging as a promising drug target for Alzheimer's disease (AD) therapy, glycogen synthase kinase 3β (GSK-3β) has garnered attention. This study sought to rigorously scrutinize a compendium of natural compounds retrieved from the ZINC database through pharmacodynamic experiments, employing a 1 H-indazole-3-carboxamide (INDZ) scaffold, to identify compounds capable of inhibiting the GSK-3β protein. Utilizing a multi-step approach, the study involved pharmacophore analysis, followed by molecular docking to select five promising ligands for further investigation. Subsequently, ESMACS simulations were employed to assess the stability of the ligand-protein interactions. Evaluation of the binding modes and free energy of the ligands revealed that five compounds (2a-6a) exhibited crucial interactions with the active site residues. Furthermore, various methodologies, including hydrogen bond and clustering analyses, were utilized to ascertain their inhibitory potential and elucidate the factors contributing to ligand binding in the protein's active site. The findings from MMPBSA/GBSA analysis indicated that these five selected small molecules closely approached the IC50 value of the reference ligand (OH8), yielding energy values of -34.85, -32.58, -31.71, and -30.39 kcal/mol, respectively. Additionally, an assessment of the interactions using hydrogen bond and dynamic analyses delineated the effective binding of the ligands with the binding pockets in the protein. Through computational analysis, we obtained valuable insights into the molecular mechanisms of GSK-3β, aiding in the development of more potent inhibitors.

Elucidating the binding mechanisms of GABA and Muscimol as an avenue to discover novel GABA-mimetic small molecules

Gamma-aminobutyric acid (GABA) signaling is the principal inhibitory pathway in the central nervous system. It is critical in neuronal cell proliferation and fate determination. Any aberration in GABA inhibition results in psychiatric and neurological diseases. Thus, modulating GABAergic neurotransmission has become the basis of drug therapy for psychiatric and several neurological diseases. Though GABA and muscimol are classical inhibitors of GABA receptors, the search for novel inhibitors continues unabated. In this study, the binding mechanism of GABA and muscimol was elucidated and applied in the search for small molecule GABAergic inhibitors using comprehensive computational techniques. It was revealed that a high-affinity binding of GABA and muscimol was mediated by a water molecule involving α1Thr129 and then stabilized by strong interactions including salt bridges with β2Glu155 and α1Arg66 amidst hydrogen bonds, π-π stacking, and π -cation interactions with other residues. The binding of GABA and muscimol was also characterized by stability and deeper penetration into the hydrophobic core of the protein which resulted in conformational changes of the binding pocket and domain, by inducing correlated motions of the residues. Thermodynamics analysis showed GABA and muscimol exhibited total binding free energies of -19.85 ± 8.83 Kcal/mol and -26.55 ± 3.42 Kcal/mol, respectively. A pharmacophore model search, based on the energy contributions of implicating binding residues, resulted in the identification of ZINC68604167, ZINC19735138, ZINC04202466, ZINC00901626, and ZINC01532854 as potential GABA-mimetic compounds from metabolites and natural products libraries. This study has elucidated the binding mechanisms of GABA and muscimol and successfully applied in the identification of GABA-mimetic compounds.Communicated by Ramaswamy H. Sarma.

Pinpointing top inhibitors for GSK3β from pool of indirubin derivatives using rigorous computational workflow and their validation using molecular dynamics (MD) simulations

Glycogen synthase kinase-3β (GSK3β) is a pivotal protein kinase implicated in a spectrum of debilitating diseases, encompassing cancer, diabetes, and neurodegenerative disorders. While the therapeutic potential of GSK3β inhibition is widely recognized, there remains an unmet need for a rigorous, systematic analysis probing the theoretical inhibition dynamics of a comprehensive library of indirubin derivatives against GSK3β using advanced computational methodologies. Addressing this gap, this study embarked on an ambitious endeavor, leveraging indirubin-a renowned scaffold-as a template to curate a vast library of 1000 indirubin derivatives from PubChem. These were enriched with varied substitutions and modifications, identified via a structure similarity search with a Tanimoto similarity threshold of 85%. Harnessing a robust virtual screening workflow, we meticulously identified the top 10 contenders based on XP docking scores. Delving deeper, we gauged the binding free energy differentials (ΔGBind) of these hits, spotlighting the top three compounds that showcased unparalleled binding prowess. A comparative pharmacophore feature mapping with the reference inhibitor OH8, co-crystallized with GSK3β (PDB ID: 6Y9R), was undertaken. The binding dynamics of these elite compounds were further corroborated with 100 ns molecular dynamics simulations, underlining their stable and potent interactions with GSK3β. Remarkably, our findings unveil that these indirubin derivatives not only match but, in certain scenarios, surpass the binding affinity and specificity of OH8. By bridging this research chasm, our study amplifies the therapeutic promise of indirubin derivatives, positioning them as frontrunners in the quest for groundbreaking GSK3β inhibitors, potentially revolutionizing treatments for a myriad of ailments.

Multi-dimensional structural footprint identification for the design of potential scaffolds targeting METTL3 in cancer treatment from natural compounds

CONTEXT

[Formula: see text]-adenosine-methyltransferase (METTL3) is the catalytic domain of the 'writer' proteins which is involved in the post modifications of [Formula: see text]-methyladinosine ([Formula: see text]). Though its activities are essential in many biological processes, it has been implicated in several types of cancer. Thus, drug developers and researchers are relentlessly in search of small molecule inhibitors that can ameliorate the oncogenic activities of METTL3. Currently, STM2457 is a potent, highly selective inhibitor of METTL3 but is yet to be approved.

METHODS

In this study, we employed structure-based virtual screening through consensus docking by using AutoDock Vina in PyRx interface and Glide virtual screening workflow of Schrodinger Glide. Thermodynamics via MM-PBSA calculations was further used to rank the compounds based on their total free binding energies. All atom molecular dynamics simulations were performed using AMBER 18 package. FF14SB force fields and Antechamber were used to parameterize the protein and compounds respectively. Post analysis of generated trajectories was analyzed with CPPTRAJ and PTRAJ modules incorporated in the AMBER package while Discovery studio and UCSF Chimera were used for visualization, and origin data tool used to plot all graphs.

RESULTS

Three compounds with total free binding energies higher than STM2457 were selected for extended molecular dynamics simulations. The compounds, SANCDB0370, SANCDB0867, and SANCDB1033, exhibited stability and deeper penetration into the hydrophobic core of the protein. They engaged in relatively stronger intermolecular interactions involving hydrogen bonds with resultant increase in stability, reduced flexibility, and decrease in the surface area of the protein available for solvent interactions suggesting an induced folding of the catalytic domain. Furthermore, in silico pharmacokinetics and physicochemical analysis of the compounds revealed good properties suggesting these compounds could serve as promising MEETL3 entry inhibitors upon modifications and optimizations as presented by natural compounds. Further biochemical testing and experimentations would aid in the discovery of effective inhibitors against the berserk activities of METTL3.

A Computational Study of the Immobilization of New 5-Nitroisatine Derivatives with the Use of C60-Based Functionalized Nanocarriers

Isatin-based compounds are a large group of drugs used as competitive inhibitors of ATP. The 5-nitroisatin derivatives studied in this work are inhibitors of the CDK2 enzyme, which can be used in the development of new anti-cancer therapies. One of the basic activities that often allows for an increase in biological activity while reducing the undesirable effects associated with the toxicity of medicinal substances is immobilization based on carriers. In this work, fifty nanocarriers derived from C60 fullerene, containing a bound phenyl ring on their surfaces, were used in the process of the immobilization of isatin derivatives. Based on flexible docking methods, the binding capacities of the drugs under consideration were determined using a wide range of nanocarriers containing symmetric and asymmetric modifications of the phenyl ring, providing various types of interactions. Based on the data collected for each of the tested drugs, including the binding affinity and the structure and stability of complexes, the best candidates were selected in terms of the type of substituent that modified the nanoparticle and its location. Among the systems with the highest affinity are the dominant complexes created by functionalized fullerenes containing substituents with a symmetrical location, such as R2-R6 and R3-R5. Based on the collected data, nanocarriers with a high potential for immobilization and use in the development of targeted therapies were selected for each of the tested drugs.

Designing and Synthesis of New Isatin Derivatives as Potential CDK2 Inhibitors

Tumors are still one of the main causes of death; therefore, the search for new therapeutic agents that will enable the implementation of effective treatment is a significant challenge for modern pharmacy. One of the important factors contributing to the development of neoplastic diseases is the overexpression of enzymes responsible for the regulation of cell division processes such as cyclin-dependent kinases. Numerous studies and examples of already-developed drugs confirm that isatin is a convenient basis for the development of new groups of inhibitors for this class of enzyme. Therefore, in this work, a new group of potential inhibitors of the CDK2 enzyme, utilizing isatin derivatives and substituted benzoylhydrazines, has been designed based on the application of computational chemistry methods, such as docking and molecular dynamics, and their inhibiting ability was assessed. In the cases of the selected compounds, a synthesis method was developed, and the selected physicochemical properties of the newly synthesized derivatives were estimated. As part of the completed project, new compounds are developed which are potential inhibitors of the CDK2 enzyme.

Multifunctional Ligands with Glycogen Synthase Kinase 3 Inhibitory Activity as a New Direction in Drug Research for Alzheimer's Disease

Alzheimer's disease (AD) belongs to the most common forms of dementia that causes a progressive loss of brain cells and leads to memory impairment and decline of other thinking skills. There is yet no effective treatment for AD; hence, the search for new drugs that could improve memory and other cognitive functions is one of the hot research topics worldwide. Scientific efforts are also directed toward combating behavioral and psychological symptoms of dementia, which are an integral part of the disease. Several studies have indicated that glycogen synthase kinase 3 beta (GSK3β) plays a crucial role in the pathogenesis of AD. Moreover, GSK3β inhibition provided beneficial effects on memory improvement in multiple animal models of AD. The present review aimed to update the most recent reports on the discovery of novel multifunctional ligands with GSK3β inhibitory activity as potential drugs for the symptomatic and disease-modifying therapy of AD. Compounds with GSK3β inhibitory activity seem to be an effective pharmacological approach for treating the causes and symptoms of AD as they reduced neuroinflammation and pathological hallmarks in animal models of AD and provided relief from cognitive and neuropsychiatric symptoms. These compounds have the potential to be used as drugs for the treatment of AD, but their precise pharmacological, pharmacokinetic, toxicological and clinical profiles need to be defined.

The Oxindole Derivatives, New Promising GSK-3β Inhibitors as One of the Potential Treatments for Alzheimer's Disease-A Molecular Dynamics Approach

Simple Summary

Enzymatic overexpression is a determinant of the development of many diseases. Increased activity of the GSK-3β enzyme is a factor that manifests itself in the development of numerous disease entities such as Alzheimer’s disease, schizophrenia, diabetes and cancers. An important medical procedure in such cases is the inhibition of enzyme activity. Based on the comprehensive use of computational chemistry methods, a group of new compounds derived from 2-oxindole was designed. The conducted research allowed the assessment of the conformational properties of the ligand molecules in the GSK-3β active site, the dynamic stability of the obtained complexes and their exact energetic characteristics. Taking into account the obtained data, a narrow group of derivatives showing an affinity for the active site of the GSK-3β enzyme was selected. The comparison of binding properties of selected 2-oxindole derivatives with an inhibitor with confirmed pharmacological activity indicates the high application potential of the newly developed compounds.

Abstract

The glycogen synthase kinase 3β (GSK-3β) is a protein kinase involved in regulating numerous physiological processes such as embryonic development, transcription, insulin action, cell division cycle and multiple neuronal functions. The overexpression of this enzyme is related to many diseases such as schizophrenia, Alzheimer’s disease, diabetes and cancer. One of the basic methods of treatment in these cases is the usage of ATP-competitive inhibitors. A significant group of such compounds are indirubin and its analogs, e.g., oxindole derivatives. The compounds considered in this work are 112 newly designed oxindole derivatives. In the first stage, such molecular properties of considered compounds as toxicity and LogP were estimated. The preliminary analysis of the binding capabilities of considered compounds towards the GSK-3β active site was conducted with the use of the docking procedure. Based on obtained molecular properties and docking simulations, a selected group of complexes that were analyzed in the molecular dynamics stage was nominated. The proposed procedure allowed for the identification of compounds such as Oxind_4_9 and Oxind_13_10, which create stable complexes with GSK-3β enzyme and are characterized by the highest values of binding affinity. The key interactions responsible for stabilization of considered ligand–protein complexes were identified, and their dynamic stability was also determined. Comparative analysis including analyzed compounds and reference molecule 3a, which is also an oxindole derivative with a confirmed inhibitory potential towards GSK3B protein, clearly indicates that the proposed compounds exhibit an analogous binding mechanism, and the obtained binding enthalpy values indicate a slightly higher binding potential than the reference molecule.

Screening of potential drug for Alzheimer's disease: a computational study with GSK-3 β inhibition through virtual screening, docking, and molecular dynamics simulation

The global impact of Alzheimer's disease (AD) necessitates intensive research to find appropriate and effective drugs. Many studies in AD suggested beta-amyloid plaques and neurofibrillary tangles-associated tau protein as the key targets for drug development. On the other hand, it is proved that triggering of Glycogen Synthase Kinase-3β (GSK-3β) also cause AD, therefore, GSK-3β is a potential drug target to combat AD. We, in this study, investigated the ability of small molecules to inhibit GSK-3β through virtual screening, Absorption, Distribution, Metabolism, and Excretion (ADME), induced-fit docking (IFD), molecular dynamics simulation, and binding free energy calculation. Besides, molecular docking was performed to reveal the binding and interaction of the ligand at the active site of GSK-3β. We found two compounds such as 6961 and 6966, which exhibited steady-state interaction with GSK-3β for 30 ns in molecular dynamics simulation. The compounds (6961 and 6966) also achieved a docking score of -9.05 kcal/mol and -8.11 kcal/mol, respectively, which is relatively higher than the GSK-3β II inhibitor (-6.73 kcal/mol). The molecular dynamics simulation revealed that the compounds have a stable state during overall simulation time, and lesser root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) values compared with co-crystal. In conclusion, we suggest the two compounds (6966 and 6961) as potential leads that could be utilized as effective inhibitors of GSK-3β to combat AD.

The Immobilization of ChEMBL474807 Molecules Using Different Classes of Nanostructures

Indirubin derivatives and analogues are a large group of compounds which are widely and successfully used in treatment of many cancer diseases. In particular, the ChEMBL474807 molecule, which has confirmed inhibiting abilities against CDK2 and GSK3B enzymes, can be included in this group. The immobilization of inhibitors with the use of nanocarriers is an often used strategy in creation of targeted therapies. Evaluations were made of the possibility of immobilizing ligand molecules on different types of nanocarrier, such as carbon nanotubes (CNT), functionalized fullerene C60 derivatives (FF_X), and functionalized cube rhombellanes, via the use of docking methods. All results were compared with a reference system, namely C60 fullerene. The realized calculations allowed indication of a group of compounds that exhibited significant binding affinity relative to the ligand molecule. Obtained data shows that structural modifications, such as those related to the addition of functional groups or changes of structure symmetry, realized in particular types of considered nanostructures, can contribute to increases of their binding capabilities. The analysis of all obtained nano complexes clearly shows that the dominant role in stabilization of such systems is played by stacking and hydrophobic interactions. The realized research allowed identification of potential nanostructures that, together with the ChEMBL474807 molecule, enable the creation of targeted therapy.

Investigation of the Inhibition Potential of New Oxindole Derivatives and Assessment of Their Usefulness for Targeted Therapy

Oxindole derivatives are a large group of compounds that can play the role of Adenosine triphosphate (ATP) competitive inhibitors. The possibility of modification of such compounds by addition of active groups to both cyclic systems of oxindole allows the obtaining of derivatives showing significant affinity toward cyclin-dependent kinase (CDK) proteins. Overexpression of that enzyme is observed in the case of most cancers. The discovery of new efficient inhibitors, which could be used in the development of targeted therapies, is one of the current goals setting trends in recent research. In this research, an oxindole molecular core was used, which was modified by the addition of different substituents to both side chains. The realized procedure allowed the creation of a set of oxindole derivatives characterized by binding affinity values and molecular descriptors evaluated during docking procedures and QSAR calculations. The most promising structures characterized by best sets of parameters were used during the molecular dynamics stage. The analysis of structural and energetic properties of systems obtained during this stage of computation gives an indication of inhibitors creating the most stable complexes, characterized by the highest affinity. During this stage, two structures were selected, where affinity towards potential nanocarriers was evaluated. Realized calculations confirmed a significant role of stacking interactions in the stabilization of ligand complexes with fullerene molecules. Obtained data indicates that complexes of oxindole derivatives and considered nanocarriers exhibit significant potential in the creation of immobilized drugs, and can be used in the development of targeted therapies.

The Immobilization of Oxindole Derivatives with Use of Cube Rhombellane Homeomorphs

A key aspect of modern drug research is the development of delivery methods that ensure the possibility of implementing targeted therapy for a specific biological target. The use of nanocarriers enables to achieve this objective, also allowing to reduce the toxicity of used substances and often extending their bioavailability. Through the application of docking methods, the possibility of using cube rhombellanes as potential carriers for two oxindole derivatives was analyzed. In the studies, compounds identified as inhibitors of the CDK2 enzyme and a set of nanostructures proposed by the Topo Cluj Group were used. The popular fullerene molecule C60 was used as the reference system. The estimated binding affinities and structures of obtained complexes show that use of functionalized cube rhombellanes containing hydrogen bond donors and acceptors in their external molecular shell significantly increases ligand affinity toward considered nanocariers, compared to classic fullerenes. The presented values also allow to state that an important factor determining the mutual affinity of the tested ligands and nanostructures is the symmetry of the analyzed nanocarriers and its influence on the distribution of binding groups (aromatic systems, donors and acceptors of hydrogen bonds) on the surface of nanoparticles.

Recent Developments and Applications of the MMPBSA Method

The Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) approach has been widely applied as an efficient and reliable free energy simulation method to model molecular recognition, such as for protein-ligand binding interactions. In this review, we focus on recent developments and applications of the MMPBSA method. The methodology review covers solvation terms, the entropy term, extensions to membrane proteins and high-speed screening, and new automation toolkits. Recent applications in various important biomedical and chemical fields are also reviewed. We conclude with a few future directions aimed at making MMPBSA a more robust and efficient method.

Inhibition mechanism of CDK-2 and GSK-3β by a sulfamoylphenyl derivative of indoline-a molecular dynamics study

A good understanding of the inhibition mechanism of enzymes exhibiting high levels of similarity is the first step to the discovery of new drugs with selective potential. Examples of such proteins include glycogen synthase kinase-3 (GSK-3β) and cyclin-dependent kinase 2 (CDK-2). This article reports the mechanism of such enzyme inhibition as analyzed by an indoline sulfamylophenyl derivative (CHEMBL410072). Previous work has shown that such compounds exhibit selective properties towards their biological targets. This study used a combined procedure involving docking and molecular dynamics simulations, which allowed identification of interactions responsible for stabilization of complexes, and analysis of the dynamic stability of the systems obtained. The initial data obtained during the molecular docking stage show that the ligand molecule exhibits a similar affinity towards both active sites, which was confirmed by quantification of identified interactions and energy values. However, the data do not cover dynamic aspects of the considered systems. Molecular dynamics simulations realized for both complexes indicate significant dissimilarities in dynamics properties of both side chains of the considered ligands, especially in the case of the part containing the sulfamide group. Such increased mobility of the analyzed systems disrupts the stability of binding in the stabilized complex with GSK-3β protein, which finally affects also the binding affinity of the ligand molecule towards this enzyme.

Potential inhibitory effect of indolizine derivatives on the two enzymes: nicotinamide phosphoribosyltransferase and beta lactamase, a molecular dynamics study

Nicotinamide phosphoribosyl-transferases (NAMPT) are enzymes that play a role in targeting cancer metabolism, while beta lactamases are involved in bacterial resistance to beta-lactam antibiotics. Many protein inhibitors exhibit such property which is often correlated with their cellular potency. In order to understand such a phenomenon, the present article conducts an analysis of the dynamic behavior of complexes formed by the inhibitors, that is indolizine derivatives, with the studied enzymes. Both docking and molecular dynamics led to identification of their interactions and showed the mechanism of inhibition of the two studied enzymes. The differences in the behavior of ligand at the active sites of beta lactamases and nicotinamide phosphoribosyl-transferases are indicated by structural and enthalpy values.