Understanding the structural and energetic requisites of ligand binding toward its molecular target is of paramount relevance in drug design. In recent years, atomistic free energy calculations have proven to be a valid tool to complement experiments in characterizing the thermodynamic and kinetic properties of protein/ligand interaction. Here, we investigate, through a recently developed metadynamics-based protocol, the unbinding mechanism of an inhibitor of the pharmacologically relevant target p38 MAP kinase. We provide a thorough description of the ligand unbinding pathway identifying the most stable binding mode and other thermodynamically relevant poses. From our simulations, we estimated the unbinding rate as k_off = 0.020 ± 0.011 s^–1. This is in good agreement with the experimental value (k_off = 0.14 s^–1). Next, we developed a Markov state model that allowed identifying the rate-limiting step of the ligand unbinding process. Our calculations further show that the solvation of the ligand and that of the active site play crucial roles in the unbinding process. This study paves the way to investigations on the unbinding dynamics of more complex p38 inhibitors and other pharmacologically relevant inhibitors in general, demonstrating that metadynamics can be a powerful tool in designing new drugs with engineered binding/unbinding kinetics.

Significance: Converging evidence points to the build-up of phosphorylated α-synuclein (α-syn) at residue serine 129 (pS129) in Lewy body disease, suggesting its central role in the regulation of α-syn aggregation and neuronal degeneration. However, a comprehensive understanding of the role of α-syn phosphorylation at pS129 in α-synuclenopathies pathogenesis is still lacking. Herein, we study the phosphorylation incidence and its effect on α-syn aggregation propensity and cellular toxicity. Collectively, our data suggest that pS129 occurred subsequent to initial α-syn aggregation, lessened aggregation propensity, and attenuated cytotoxicity through diverse assays. Our findings highlight major implications for a better understanding of the role of a molecular modification on protein aggregation.

We present a novel approach to explore the conformational space of globular proteins near their native state. It combines the advantages of coarse-grained models with those of all-atoms simulations, required to treat molecular recognition processes. The comparison between calculated structural properties with those obtained with all-atoms molecular dynamics simulations establishes the accuracy of the model. Our method has the potential to be extended to molecular recognition processes in systems whose characteristic size and time scale prevent an analysis based on all-atoms molecular dynamics.

The prolyl isomerase PIN1, a critical modifier of multiple signalling pathways, is overexpressed in the majority of cancers and its activity strongly contributes to tumour initiation and progression. Inactivation of PIN1 function conversely curbs tumour growth and cancer stem cell expansion, restores chemosensitivity and blocks metastatic spread, thus providing the rationale for a therapeutic strategy based on PIN1 inhibition. Notwithstanding, potent PIN1 inhibitors are still missing from the arsenal of anti-cancer drugs. By a mechanism-based screening, we have identified a novel covalent PIN1 inhibitor, KPT-6566, able to selectively inhibit PIN1 and target it for degradation. We demonstrate that KPT-6566 covalently binds to the catalytic site of PIN1. This interaction results in the release of a quinone-mimicking drug that generates reactive oxygen species and DNA damage, inducing cell death specifically in cancer cells. Accordingly, KPT-6566 treatment impairs PIN1-dependent cancer phenotypes in vitro and growth of lung metastasis in vivo.

Fast lateral proton migration along membranes is of vital importance for cellular energy homeostasis and various proton-coupled transport processes. It can only occur if attractive forces keep the proton at the interface. How to reconcile this high affinity to the membrane surface with high proton mobility is unclear. Here, we tested whether a minimalistic model interface between an apolar hydrophobic phase ( n -decane) and an aqueous phase mimics the biological pathway for lateral proton migration. The observed diffusion span, on the order of tens of micrometers, and the high proton mobility were both similar to the values previously reported for lipid bilayers. Extensive ab initio simulations on the same water/ n -decane interface reproduced the experimentally derived free energy barrier for the excess proton. The free energy profile G _{H ⁺} adopts the shape of a well at the interface, having a width of two water molecules and a depth of 6 ± 2 RT . The hydroniums in direct contact with n -decane have a reduced mobility. However, the hydroniums in the second layer of water molecules are mobile. Their in silico diffusion coefficient matches that derived from our in vitro experiments, (5.7 ± 0.7) × 10 ^-5 cm ² s ^-1 . Conceivably, these are the protons that allow for fast diffusion along biological membranes.