Conformational Fluctuations in β2-Microglubulin Using Markov State Modeling and Molecular Dynamics

Computational exploration of allosteric inhibitors targeting CDK4/CDK6 proteins: a promising approach for multi-target drug development

Cyclin-dependent kinases (CDKs) play a pivotal role in orchestrating the intricate regulation of the cell cycle, a fundamental process governing cell growth and division. In particular, CDK4 and CDK6 are critical for the transition from the G1 phase to the S phase, where Deoxyribonucleic acid (DNA) replication occurs, and their dysregulation is linked to various diseases, notably cancer. While ATP-binding site inhibitors for CDKs are well-documented, this study focuses on uncovering allosteric inhibitors, providing a fresh perspective on CDK inhibition. Computational techniques were employed in this investigation, utilizing Molecular Operating Environment (MOE) for virtual screening of a drug-like compound library. Moreover, the stability of the most promising binding inhibitors was assessed through Molecular Dynamics (MD) simulations and MMPBSA/MMGBSA analyses. The outcome reveals that three inhibitors (C1, C2, and C3) exhibited the strongest binding affinity for CDK4/CDK6, as corroborated by docking and simulation analyses. The computed binding energies ranged from -6.1 to -7.6 kcal/mol, underscoring the potency of these allosteric inhibitors. Notably, this study identifies key residues (PHE31, HIS95, HIS100, VAL101, ASP102, ASP104, and THR107) that play pivotal roles in mediating inhibitor binding within the allosteric sites. Among the findings, the C1-CDK4 complex and C2-CDK6 complex emerge as particularly promising inhibitors, exhibiting high binding energies, favorable interaction patterns, and sustained presence within the active site. This study contributes significantly to the pursuit of multi-target drugs against CDK4/CDK6 proteins, with potential implications for the development of innovative therapies across various disorders, including cancer and other cell cycle-related conditions.

A Self-Consistent Molecular Mechanism of β2-Microglobulin Aggregation

Despite the consensus on the origin of dialysis-related amyloidosis (DRA) being β2-microglobulin (β2m) aggregation, the debate on the underlying mechanism persists because of the continuous emergence of β2m variant- and pH-dependent contradictory results. By characterizing the native monomeric (initiation) and aggregated fibrillar (termination) states of β2m via a combination of two enhanced sampling approaches, we here propose a mechanism that explains the heterogeneous behavior of wild-type (WT) and pathogenic (V27M and D76N) β2m variants in physiological and disease-pertinent acidic pH environments. It appears that the higher retainment of monomeric native folds at neutral pH (native-like) distinguishes pathogenic β2m mutants from the WT (moderate loss). However, at acidic pH, all three variants behave similarly in producing a substantial amount of partially unfolded states (conformational switch, propensity), though with different extents (WT < V27M < D76N). Whereas at the fibrillar end, all β2m variants display a pH-dependent protofilament separation pathway and a higher protofilament binding affinity (stability) at acidic pH, where the relative order of binding affinity (WT < V27M < D76N) remains consistent with pH modulation. Combining these observations, we conclude that β2m variants possibly shift from native-like aggregation to conformational switch-initiated fibrillation as the pH is altered from neutral to acidic. The combined propensity-stability approach based on the initiation and termination points of β2m aggregation not only assists us in deciphering the mechanism but also emphasizes the protagonistic roles of both terminal points in the overall aggregation process.