Structural basis of valmerins as dual inhibitors of GSK3β/CDK5

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.

J, Ahmad SF, Attia SM, Emran TB, Patil RB, Ahmed SSSJ. Structural Characteristics of PON1 with Leu55Met and Gln192Arg Variants Influencing Oxidative-Stress-Related Diseases: An Integrated Molecular Modeling and Dynamics Study

Background and Objectives: PON1 is a multi-functional antioxidant protein that hydrolyzes a variety of endogenous and exogenous substrates in the human system. Growing evidence suggests that the Leu55Met and Gln192Arg substitutions alter PON1 activity and are linked with a variety of oxidative-stress-related diseases. Materials and Methods: We implemented structural modeling and molecular dynamics (MD) simulation along with essential dynamics of PON1 and molecular docking with their endogenous (n = 4) and exogenous (n = 6) substrates to gain insights into conformational changes and binding affinity in order to characterize the specific functional ramifications of PON1 variants. Results: The Leu55Met variation had a higher root mean square deviation (0.249 nm) than the wild type (0.216 nm) and Gln192Arg (0.202 nm), implying increased protein flexibility. Furthermore, the essential dynamics analysis confirms the structural change in PON1 with Leu55Met vs. Gln192Arg and wild type. Additionally, PON1 with Leu55Met causes local conformational alterations at the substrate binding site, leading to changes in binding affinity with their substrates. Conclusions: Our findings highlight the structural consequences of the variants, which would increase understanding of the role of PON1 in the pathogenesis of oxidative-stress-related diseases, as well as the management of endogenous and exogenous chemicals in the treatment of diseases.

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.

Multispectroscopic and Computational Investigations on the Binding Mechanism of Dicaffeoylquinic Acids with Ovalbumin

Recently, studies on the interactions between ovalbumin (OVA) and polyphenols have received a great deal of interest. This study explored the conformational changes and the interaction mechanism of the binding between OVA and chlorogenic acid (CGA) isomers such as 3,4-dicaffeoylquinic acids (3,4-diCQA), 4,5-dicaffeoylquinic acids (4,5-diCQA), and 3,5-dicaffeoylquinic acids (3,5-diCQA) using multispectroscopic and in silico analyses. The emission spectra show that the diCQAs caused strong quenching of OVA fluorescence under different temperatures through a static quenching mechanism with hydrogen bond (H-bond) and van der Waals (vdW) interactions. The values of binding constants (OVA-3,4-diCQA = 6.123 × 10⁵, OVA-3,5-diCQA = 2.485 × 10⁵, OVA-4,5-diCQA = 4.698 × 10⁵ dm³ mol^-1 at 298 K) suggested that diCQAs had a strong binding affinity toward OVA, among which OVA-3,4-diCQA exhibits higher binding constant. The results of UV-vis absorption and synchronous fluorescence indicated that the binding of all three diCQAs to OVA induced conformational and micro-environmental changes in the protein. The findings of molecular modeling further validate the significant role of vdW force and H-bond interactions in ensuring the stable binding of OVA-diCQA complexes. Temperature-dependent molecular dynamics simulation studies allow estimation of the individual components that contribute to the total bound free energy value, which allows evaluation of the nature of the interactions involved. This research can provide information for future investigations on food proteins' physicochemical stability and CGA bioavailability in vitro or in vivo.

Efficient synthesis of isoindolones by intramolecular cyclisation of pyridinylbenzoic acids

A straightforward one-pot method for the synthesis of unreported pyrido-[2,1-a]isoindolones in excellent yield is described. Two novel isoindolones were synthesized and fully characterized. The alkyl substituents on the pyridine play an important role in the outcome of the reaction. The mechanism, investigated through DFT calculations, features an unprecendented intramolecular cyclization reaction involving a carboxylic acid activated by tosyl chloride and an electron-poor pyridinic nitrogen. This protocol completes the known strategies to obtain functionalized isoindolones.

When Good Kinases Go Rogue: GSK3, p38 MAPK and CDKs as Therapeutic Targets for Alzheimer's and Huntington's Disease

Alzheimer's disease (AD) is a mostly sporadic brain disorder characterized by cognitive decline resulting from selective neurodegeneration in the hippocampus and cerebral cortex whereas Huntington's disease (HD) is a monogenic inherited disorder characterized by motor abnormalities and psychiatric disturbances resulting from selective neurodegeneration in the striatum. Although there have been numerous clinical trials for these diseases, they have been unsuccessful. Research conducted over the past three decades by a large number of laboratories has demonstrated that abnormal actions of common kinases play a key role in the pathogenesis of both AD and HD as well as several other neurodegenerative diseases. Prominent among these kinases are glycogen synthase kinase (GSK3), p38 mitogen-activated protein kinase (MAPK) and some of the cyclin-dependent kinases (CDKs). After a brief summary of the molecular and cell biology of AD and HD this review covers what is known about the role of these three groups of kinases in the brain and in the pathogenesis of the two neurodegenerative disorders. The potential of targeting GSK3, p38 MAPK and CDKS as effective therapeutics is also discussed as is a brief discussion on the utilization of recently developed drugs that simultaneously target two or all three of these groups of kinases. Multi-kinase inhibitors either by themselves or in combination with strategies currently being used such as immunotherapy or secretase inhibitors for AD and knockdown for HD could represent a more effective therapeutic approach for these fatal neurodegenerative diseases.

Screening potential FDA-approved inhibitors of the SARS-CoV-2 major protease 3CLpro through high-throughput virtual screening and molecular dynamics simulation

It has been confirmed that the new coronavirus SARS-CoV-2 caused the global pandemic of coronavirus disease 2019 (COVID-19). Studies have found that 3-chymotrypsin-like protease (3CL^pro) is an essential enzyme for virus replication, and could be used as a potential target to inhibit SARS-CoV-2. In this work, 3CL^pro was used as the target to complete the high-throughput virtual screening of the FDA-approved drugs, and Indinavir and other 10 drugs with high docking scores for 3CL^pro were obtained. Studies on the binding pattern of 3CL^pro and Indinavir found that Indinavir could form the stable hydrogen bond (H-bond) interactions with the catalytic dyad residues His41-Cys145. Binding free energy study found that Indinavir had high binding affinity with 3CL^pro. Subsequently, molecular dynamics simulations were performed on the 3CL^pro and 3CL^pro-Indinavir systems, respectively. The post-dynamic analyses showed that the conformational state of the 3CLpro-Indinavir system transformed significantly and the system tended to be more stable. Moreover, analyses of the residue interaction network (RIN) and H-bond occupancy revealed that the residue-residue interaction at the catalytic site of 3CL^pro was significantly enhanced after binding with Indinavir, which in turn inactivated the protein. In short, through this research, we hope to provide more valuable clues against COVID-19.

Design, synthesis, biological evaluation and molecular dynamics studies of 4-thiazolinone derivatives as protein tyrosine phosphatase 1B (PTP1B) inhibitors

Protein tyrosine phosphatase 1B (PTP1B) is a key negative regulator of insulin signaling pathway, and more and more studies have shown that it is a potential target for the treatment of type 2 diabetes mellitus (T2DM). In this study, 17 new 4-thiazolinone derivatives were designed and synthesized as novel PTP1B inhibitors, and ADMET prediction confirmed that these compounds were to be drug-like. In vitro enzyme activity experiments were performed on these compounds, and it was found that a plurality of compounds had good inhibitory activity and high selectivity against PTP1B protein. Among them, compound 7p exhibited the best inhibitory activity with an IC₅₀ of 0.92 μM. The binding mode of compound 7p and PTP1B protein was explored, revealing the reason for its high efficiency. In addition, molecular dynamics simulations for the PTP1B^WT and PTP1B^comp#7p systems revealed the effects of compound 7p on PTP1B protein at the molecular level. In summary, the study reported for the first time that 4-thiazolinone derivatives as a novel PTP1B inhibitor had good inhibitory activity and selectivity for the treatment of T2DM, providing more options for the development of PTP1B inhibitors. AbbreviationsBBBblood-brain barrierCDC25Bcell division cycle 25 homolog BCYP2D6Cytochrome P450 2D6 bindingDCCMdynamic cross-correlation mapDSDiscovery StudioH bondhydrogen bondHIAhuman intestinal absorptionLARleukocyte antigen-related phosphataseMDmolecular dynamicsMEG-2maternal-effect germ-cell defective 2MM-PBSAmolecular mechanics Poisson Boltzmann surface area)PCAprincipal component analysisPDBProtein Data BankpNPPp-nitrophenyl phosphatePPBplasma protein bindingPTP1Bprotein tyrosine phosphotase 1BRMSDroot mean square deviationRMSFroot mean square fluctuationSHP-1src homologous phosphatase-1SHP-2src homologous phosphatase-2SPCsingle-point chargeTCPTPT cell protein tyrosine phosphataseT2DMType 2 diabetes mellitusVDWvan der WaalsCommunicated by Ramaswamy H. Sarma.

How to design potent and selective DYRK1B inhibitors? Molecular modeling study

DYRK1B protein kinase is an emerging anticancer target due to its overexpression in a variety of cancers and its role in cancer chemoresistance through maintaining cancer cells in the G0 (quiescent) state. Consequently, there is a growing interest in the development of potent and selective DYRK1B inhibitors for anticancer therapy. One of the major off-targets is another protein kinase, GSK3β, which phosphorylates an important regulator of cell cycle progression on the same residue as DYRK1B and is involved in multiple signaling pathways. In the current work, we performed a detailed comparative structural analysis of DYRK1B and GSK3β ATP-binding sites and identified key regions responsible for selectivity. As the crystal structure of DYRK1B has never been reported, we built and optimized a homology model by comparative modeling and metadynamics simulations. Calculation of interaction energies between docked ligands in the ATP-binding sites of both kinases allowed us to pinpoint key residues responsible for potency and selectivity. Specifically, the role of the gatekeeper residues in DYRK1B and GSK3β is discussed in detail, and two other residues are identified as key to selectivity of DYRK1B inhibition versus GSK3β. The analysis presented in this work was used to support the design of potent and selective azaindole-quinoline-based DYRK1B inhibitors and can facilitate development of more selective inhibitors for DYRK kinases.

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.

Designing of dual inhibitors for GSK-3β and CDK5: Virtual screening and in vitro biological activities study

Alzheimer's disease is a multifactorial neurodegenerative disorder with many drug targets contributing to its etiology. Despite the devastating effects of this disease, therapeutic methods for treating Alzheimer's disease remain limited. The multifactorial nature of Alzheimer's disease strongly supports a multi-target rationale as a drug design strategy. Glycogen synthase kinase-3 beta and cyclin-dependent kinase 5 have been identified as being involved in the pathological hyperphosphorylation of tau proteins, which leads to the formation of neurofibrillary tangles and causes Alzheimer's disease. In this study, using a molecular docking method to screen a virtual library, we discovered molecules that can simultaneously inhibit Glycogen synthase kinase-3 beta and cyclin-dependent kinase 5 as lead compounds for the treatment of Alzheimer's disease. The docking results revealed the key residues in the substrate binding sites of both Glycogen synthase kinase-3 beta and cyclin-dependent kinase 5. A receiver operating characteristic curve indicated that the docking model consistently and selectively scored the majority of active compounds above decoys. The pre-treatment of cells with screened compounds protected them against Aβ_25-35- induced cell death by up to 80%. Collectively, these findings suggest that some compounds have potential to be promising multifunctional agents for Alzheimer's disease treatment.

Underlying mechanisms of cyclic peptide inhibitors interrupting the interaction of CK2α/CK2β: comparative molecular dynamics simulation studies

Protein-protein interactions (PPIs) are fundamental to all biological processes. Recently, the CK2β-derived cyclic peptide Pc has been demonstrated to efficiently antagonize the CK2α/CK2β interaction and strongly affect the phosphorylation of CK2β-dependent CK2 substrate specificity. The binding affinity of Pc to CK2α is destroyed to different extents by two single-point mutations of Tyr188 to Ala (Y188A) and Phe190 to Ala (F190A), which exert negative effects on the inhibitory activity (IC50) of Pc against the CK2α/CK2β interaction from 3.0 μM to 54.0 μM and ≫100 μM, respectively. However, the structural influences of Y188A and F190A mutations on the CK2α-Pc complex remain unclear. In this study, comparative molecular dynamics (MD) simulations, principal component analysis (PCA), domain cross-correlation map (DCCM) analysis and energy calculations were performed on wild type (WT), Y188A mutant, and F190A mutant systems. The results revealed that ordered communications between hydrophobic and polar interactions were essential for CK2α-Pc binding in the WT system. In addition to the loss of the hydrogen bond between Gln36 of CK2α and Gly189 of Pc in the two mutants, the improper recognition mechanisms occurred through different pathways. These pathways included the weakened hydrophobic interactions in the Y188A mutant as well as decreased polar and hydrophobic interactions in the F190A mutant. The energy analysis results qualitatively elucidated the instability of the two mutants and energetic contributions of the key residues. This study not only revealed the structural mechanisms for the decreased binding affinity of Y188A and F190A mutant CK2α-Pc complexes, but also provided valuable clues for the rational design of CK2α/CK2β subunit interaction inhibitors with high affinity and specificity.

The Role of Cdk5 in Alzheimer's Disease

Alzheimer's disease (AD) is known as the most fatal chronic neurodegenerative disease in adults along with progressive loss of memory and other cognitive function disorders. Cyclin-dependent kinase 5 (Cdk5), a unique member of the cyclin-dependent kinases (Cdks), is reported to intimately associate with the process of the pathogenesis of AD. Cdk5 is of vital importance in the development of CNS and neuron movements such as neuronal migration and differentiation, synaptic functions, and memory consolidation. However, when neurons suffer from pathological stimuli, Cdk5 activity becomes hyperactive and causes aberrant hyperphosphorylation of various substrates of Cdk5 like amyloid precursor protein (APP), tau and neurofilament, resulting in neurodegenerative diseases like AD. Deregulation of Cdk5 contributes to an array of pathological events in AD, ranging from formation of senile plaques and neurofibrillary tangles, synaptic damage, mitochondrial dysfunction to cell cycle reactivation as well as neuronal cell apoptosis. More importantly, an inhibition of Cdk5 activity with inhibitors such as RNA inference (RNAi) could protect from memory decline and neuronal cell loss through suppressing β-amyloid (Aβ)-induced neurotoxicity and tauopathies. This review will briefly describe the above-mentioned possible roles of Cdk5 in the physiological and pathological mechanisms of AD, further discussing recent advances and challenges in Cdk5 as a therapeutic target.