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.

The formation of a heterogeneous set of advanced glycation end products (AGEs) is the final outcome of a non-enzymatic process that occurs in vivo on long-life biomolecules. This process, known as glycation, starts with the reaction between reducing sugars, or their autoxidation products, with the amino groups of proteins, DNA, or lipids, thus gaining relevance under hyperglycemic conditions. Once AGEs are formed, they might affect the biological function of the biomacromolecule and, therefore, induce the development of pathophysiological events. In fact, the accumulation of AGEs has been pointed as a triggering factor of obesity, diabetes-related diseases, coronary artery disease, neurological disorders, or chronic renal failure, among others. Given the deleterious consequences of glycation, evolution has designed endogenous mechanisms to undo glycation or to prevent it. In addition, many exogenous molecules have also emerged as powerful glycation inhibitors. This review aims to provide an overview on what glycation is. It starts by explaining the similarities and differences between glycation and glycosylation. Then, it describes in detail the molecular mechanism underlying glycation reactions, and the bio-molecular targets with higher propensity to be glycated. Next, it discusses the precise effects of glycation on protein structure, function, and aggregation, and how computational chemistry has provided insights on these aspects. Finally, it reports the most prevalent diseases induced by glycation, and the endogenous mechanisms and the current therapeutic interventions against it.

This article reviews different formulations of the thermodynamic cycles used for the prediction of pK_a values, their advantages, and disadvantages with special emphasis on the limitations resulting from the necessity of gas‐phase calculations, which allow introducing some difficult cases that motivated alternative strategies. Before introducing the protocols that do not consider gas‐phase calculations, the two current opinions available in the literature on the debate about the correct formalism for the calculation of free energies in solution are briefly introduced. Then, the isodesmic proton exchange reaction in solution is reviewed by analyzing its performance on difficult cases for thermodynamic cycles such as carbon acids and amino acids. The pK_a values predicted by the isodesmic reaction for common acid species are also reviewed to compare their accuracy results in relation with those of thermodynamic cycles. Linear regressions between experimental pK_a values and the calculated free energies obtained with the isodesmic reaction provide expressions for the dependence of the error in the calculated pK_as on the pK_a difference between the studied acid and the reference species. Finally, it is shown that linear regressions correct the calculated free energies of the isodesmic reaction, when high constant precision is required in a large pK_a range. The deviations from the expected behavior are equivalent to those reported previously for different pK_a calculation protocols and are determined by the inaccuracies of continuum solvent models on the interactions with ionic species.