Allosteric regulation plays an important role in a myriad of bio-macromolecular processes. Specifically, in a protein, the process of allostery refers to the transmission of a local perturbation, such as ligand binding, to a distant site. Decades after the discovery of this phenomenon, models built on static images of proteins are being reconsidered with the knowledge that protein dynamics plays an important role in its function. Molecular dynamics simulations are a valuable tool for studying complex biomolecular systems, providing an atomistic description of their structure and dynamics. Unfortunately, their predictive power has been limited by the complexity of the biomolecule free-energy surface and by the length of the allosteric timescale (in the order of milliseconds). In this work, we are able to probe the origins of the allosteric changes that transcription factor mixed lineage leukemia (MLL) causes to the interactions of KIX domain of CREB-binding protein (CBP) with phosphorylated kinase inducible domain (pKID), by combing all-atom molecular dynamics with enhanced sampling methods recently developed in our group. We discuss our results in relation to previous NMR studies. We also develop a general simulations protocol to study allosteric phenomena and many other biological processes that occur in the micro/milliseconds timescale.