The p38α mitogen-activated protein kinase (MAPK) is a central regulator of cellular responses to stress and inflammation, making it an attractive therapeutic target for a broad range of diseases. Like other kinases, several types of inhibitors have been reported for p38α MAPK. These inhibitors have been classified according to the binding site they occupy on the target protein, as well as the specific conformations of key structural motifs involved in their interaction. Type II inhibitors targeting p38α MAPK bind to the inactive-like (DFG-out) conformation of the enzyme and engage adjacent allosteric pockets contiguous to the ATP-binding site, exhibiting a promising profile in terms of potency and on-target residence time. Increasing evidence indicates that drug-target residence time is a critical determinant of in vivo efficacy, and its consideration is helping to shape a new paradigm in rational drug design. In this work, we compute on-target residence times of type II p38α MAPK inhibitors using τ-Random Acceleration Molecular Dynamics (τ-RAMD) simulations. Conventional molecular dynamics simulations were first used to assess the structural stability and interaction profiles of nine type II inhibitors bound to p38α MAPK. Our results showed that τ-RAMD-derived residence times exhibited reasonable correspondence with experimental data, with a good correlation between calculated (τcalc) and experimental (τexp) residence time values. Overall, these results suggest that τ-RAMD can capture relative residence-time trends within mechanistically homogeneous inhibitor types. This study underscores the importance of binding-mode-specific validation when using enhanced sampling approaches to model drug-target kinetics.
This work describes an example of using τ-Random Acceleration Molecular Dynamics (τRAMD) in kinetic calculations.