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.

A Combination of Solid-State NMR and MD Simulations Reveals the Binding Mode of a Rhomboid Protease Inhibitor

Structure and Dynamics of the Rhomboid Protease GlpG in Liposomes Studied by Solid-State NMR.

Multi-Scaled Explorations of Binding-Induced Folding of Intrinsically Disordered Protein Inhibitor IA3 to Its Target Enzyme

Molecular Dynamics Simulations Investigate the Pathway of Substrate Entry Active Site of Rhomboid Protease

An in vitro platform for the enzymatic characterization of the rhomboid protease RHBDL4

pMD-Membrane: A Method for Ligand Binding Site Identification in Membrane-Bound Proteins

Lipid-polymer nanoparticles to probe the native-like environment of intramembrane rhomboid protease GlpG and its activity

From docking to dynamics: Unveiling the potential non-peptide and non-covalent inhibitors of M pro from natural products

Structural Analysis of a Rhomboid Family Intramembrane Protease Reveals a Gating Mechanism for Substrate Entry

A Combined STD-NMR/Molecular Modeling Protocol for Predicting the Binding Modes of the Glycosidase Inhibitors Kifunensine and Salacinol to Golgi Α-Mannosidase II

Discovery and Structure-Based Optimization of Fragments Binding the Mixed Lineage Kinase Domain-like Protein Executioner Domain

Probing the interaction mechanism of small molecule inhibitors with matriptase based on molecular dynamics simulation and free energy calculations.

Probing Ligand-Binding Modes and Binding Mechanisms of Benzoxazole-Based Amide Inhibitors with Soluble Epoxide Hydrolase by Molecular Docking and Molecular Dynamics Simulation.

Structural visualization of inhibitor binding in prolyl oligopeptidase

Probing Protein-Ligand Methyl-π Interaction Geometries through Chemical Shift Measurements of Selectively Labeled Methyl Groups

Discovery and Validation of the Binding Poses of Allosteric Fragment Hits to Protein Tyrosine Phosphatase 1b: From Molecular Dynamics Simulations to X-ray Crystallography

A multiscale approach to predict the binding mode of metallo beta‐lactamase inhibitors

Probing the binding mode and interactions of proteinase K and glutathione: molecular simulation and experiments

Importance of Three-Body Problems and Protein-Protein Interactions in Proteolysis-Targeting Chimera Modeling: Insights from Molecular Dynamics Simulations

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

Ensemble Docking as a Tool for the Rational Design of Peptidomimetic Staphylococcus aureus Sortase A Inhibitors