In this study, we investigate the atomistic details of Keap1-Nrf2 inhibitors by in-depth modeling techniques, including molecular dynamics (MD) simulations, and the path-based free energy method of umbrella sampling (US). The protein–protein interaction (PPI) of Keap1-Nrf2 is implicated in several neurodegenerative diseases like cancer, diabetes, and cardiomyopathy. A better understanding of the five sub-pocket binding sites for Nrf2 (ETGE and DLG motifs) inside the Kelch domain would expedite the inhibitor design process. We selected four protein–ligand complexes with distinct co-crystal ligands and binding occupancies inside the Nrf2 binding site. We performed 100 ns of MD simulation for each complex and analyzed the trajectories. From the results, it is evident that one ligand (1VV) has flipped inside the binding pocket, whereas the remaining three were stable. We found that Coulombic (Arg483, Arg415, Ser363, Ser508, and Ser602) and Lennard–Jones (Tyr525, Tyr334, and Tyr572) interactions played a significant role in complex stability. The obtained binding free energy values from US simulations were consistent with the potencies of simulated ligands. US simulation highlight the importance of basic and aromatic residues in the binding pocket. A detailed description of the dissociation process brings valuable insight into the interaction of the four selected protein–ligand complexes, which could help in the future to design more potent PPI inhibitors.