Abstract
In this work, we investigate the role of solvent in the binding reaction of p53 transactivation domain (TAD) peptide to its receptor MDM2. Previously, our group generated 831 μs of explicit-solvent aggregate molecular simulation trajectory data for the MDM2-p53 peptide binding reaction using large-scale distributed computing, and subsequently built a Markov State Model (MSM) of the binding reaction (Zhou et al. 2017). Here, we perform tICA anlaysis and construct an MSM using only solvent-based structural features, and find a remarkably similar landscape and implied timescales for the slowest motions. The solvent shells contributing most to the first tICA eigenvector are those centered on Lys24 and Thr18 of p53 TAD peptide in the range of 3–6 Å. Important solvent shells were visualized to reveal wetting and dewetting transitions along peptide-protein binding trajectories. Our results provide a solvent-centric view of the hydrophobic effect in action for a realistic peptide-protein binding scenario.