Abstract:Small modifications in the chemical structure of ligands are known to dramatically change their ability to inhibit the activity of a protein. Unraveling the mechanisms that govern these dramatic changes requires scrutinizing the dynamics of protein-ligand binding and unbinding at the atomic level. As an exemplary case, we have studied Glycogen Synthase Kinase-3β (GSK-3β), a multifunctional kinase that has been implicated in a host of pathological processes. As such, there is a keen interest in identifying ligands that inhibit GSK-3β activity. One family of compounds that are highly selective and potent inhibitors of GSK-3β is exemplified by a molecule termed COB-187. COB-187 consists of a five-member heterocyclic ring with a thione at C2, a pyridine substituted methyl at N3, and a hydroxyl and phenyl at C4. We have studied the inhibition of GSK-3β by COB-187-related ligands that differ in a single heavy atom from each other (either in the location of nitrogen in their pyridine ring, or with the pyridine ring replaced by a phenyl ring), or in the length of the alkyl group joining the pyridine and the N3. The inhibition experiments show a large range of half-maximal inhibitory concentration (IC50) values from 10 nM to 10 μM, implying that these ligands exhibit vastly different propensities to inhibit GSK-3β. To explain these differences, we perform Markov State Modeling (MSM) using fully atomistic simulations. Our MSM results are in excellent agreement with the experiments in that they accurately capture differences in the binding propensities of the ligands. The simulations show that the binding propensities are related to the ligands' ability to attain a compact conformation where their two aromatic rings are spatially close. We rationalize this result by sampling numerous binding and unbinding events via funnel metadynamics simulations, which show that indeed while approaching the bound state, the ligands prefer to be in their compact conformation. We find that the presence of nitrogen in the aromatic ring increases the probability of attaining the compact conformation. Protein-ligand binding is understood to be dictated by the energetics of interactions and entropic factors, like the release of bound water from the binding pockets. This work shows that changes in the conformational distribution of ligands due to atom-level modifications in the structure play an important role in protein-ligand binding.

Dissection of the Difference Between the Group I Metal Ions in Inhibiting GSK3β: a Computational Study.

The Effect of Li+ on GSK-3 Inhibition: Molecular Dynamics Simulation

Molecular Dynamics Simulations on the Role of Structural Mg2+ Ions in Phosphoryl Transfer Catalyzed by Gsk-3 Beta

Contrasting Effects of Inhibitors Li+ and Be2+ on Catalytic Cycle of Glycogen Synthase Kinase-3β

Mechanism of Kinase Inactivation and Nonbinding of FRATide to GSK3β Due to K85M Mutation: Molecular Dynamics Simulation and Normal Mode Analysis.

Molecular modeling and molecular dynamics simulation studies of the GSK3β/ATP/substrate complex: understanding the unique P+4 primed phosphorylation specificity for GSK3β substrates.

Capturing the Effects of Single Atom Substitutions on the Inhibition Efficiency of Glycogen Synthase Kinase-3β Inhibitors via Markov State Modeling and Experiments

Unveiling the Novel Characteristics of IGPD Polymer and Inhibitors Binding Affinities Using 12-6-4 LJ-type Nonbonded Mn2+ Model

Understanding the Selectivity of Inhibitors Toward PI4KIIIα and PI4KIIIβ Based Molecular Modeling

Design, Synthesis, and Evaluation of 3‐aryl‐4‐pyrrolyl‐maleimides As Glycogen Synthase Kinase‐3β Inhibitors

Microsecond Molecular Dynamics Simulations and Dynamic Network Analysis Provide Understanding of the Allosteric Inactivation of GSK3β Induced by the L343R Mutation

Effect of Double Mutations K214/A–E215/Q of FRATide on GSK3β: Insights from Molecular Dynamics Simulation and Normal Mode Analysis

Role of Mg2+ Ions in Protein Kinase Phosphorylation: Insights from Molecular Dynamics Simulations of ATP-kinase Complexes

Structural basis for the complete loss of GSK3beta catalytic activity due to R96 mutation investigated by molecular dynamics study.

Discovery and Computational Studies of Potent Covalent Kinase Inhibitors with alpha-Substituent Electrophiles Targeting Cysteine

Discovery and Computational Studies of Potent Covalent Kinase Inhibitors with -Substituent Electrophiles Targeting Cysteine

Exploring the Selectivity of Tetrahydropyrido[1,2-a]isoindolone Derivatives to GSK3 beta and CDK5 by Computational Methods

Role of Bridging Water Molecules in GSK3β‐inhibitor Complexes: Insights from QM/MM, MD, and Molecular Docking Studies

Assessing the Potential Effects of Active Site Mg2+ Ions in the Glms Ribozyme–Cofactor Complex

Studies of glycogen synthase kinase-3 inhibitors

Design, Synthesis and Evaluation of 7-Azaindazolyl-indolyl-maleimides As Glycogen Synthase Kinase-3Β (Gsk-3Β) Inhibitors