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

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.

Identifying Potential Ligand Binding Sites on Glycogen Synthase Kinase 3 Using Atomistic Cosolvent Simulations.

Glycogen synthase kinase 3 ß (GSK3ß) is a serine/threonine kinase that phosphorylates several protein substrates in crucial cell signaling pathways. Owing to its therapeutic importance, there is a need to develop GSK3ß inhibitors with high specificity and potency. One approach is to find small molecules that can allosterically bind to the GSK3ß protein surface. We have employed fully atomistic mixed-solvent molecular dynamics (MixMD) simulations to identify three plausible allosteric sites on GSK3ß that can facilitate the search for allosteric inhibitors. Our MixMD simulations narrow down the allosteric sites to precise regions on the GSK3ß surface, thereby improving upon the previous predictions of the locations of these sites.

Stability of Water Confined between Supported Self-Assembled Monolayers.

We present a thermodynamic argument showing that the evaporation and condensation free-energy barriers of water confined between two hydrophobic self-assembled monolayers (SAMs) vary more gradually with the SAM hydrophobicity as compared to the case of water confined between two bare hydrophobic surfaces (no SAMs). We validate our theory by calculating the free-energy profiles of water confined between two SAMs and between two bare surfaces of different hydrophobicities. An implication of our findings is the existence of three regimes of stability of confined water as a function of the hydrophobicity of the SAMs. In comparison to bare planar surfaces with no SAMs, the highly hydrophobic SAMs act to stabilize the liquid state, whereas weakly hydrophobic SAMs stabilize the vapor state of confined water. For intermediate hydrophobicities, the SAMs reduce both the evaporation and the condensation free-energy barriers. These results imply that the effects of SAM hydrophobicity on the behavior of confined water are nontrivial and richer than previously thought.