The Impact of Small Molecule Binding on the Energy Landscape of the Intrinsically Disordered Protein C-Myc

Intrinsically disordered proteins are attractive therapeutic targets owing to their prevalence in several diseases. Yet their lack of well-defined structure renders ligand discovery a challenging task. An intriguing example is provided by the oncoprotein c-Myc, a transcription factor that is over expressed in a broad range of cancers. Transcriptional activity of c-Myc is dependent on heterodimerization with partner protein Max. This protein-protein interaction is disrupted by the small molecule 10058-F4 (1), that binds to monomeric and disordered c-Myc. To rationalize the mechanism of inhibition, structural ensembles for the segment of the c-Myc domain that binds to 1 were computed in the absence and presence of the ligand using classical force fields and explicit solvent metadynamics molecular simulations. The accuracy of the computed structural ensembles was assessed by comparison of predicted and measured NMR chemical shifts. The small molecule 1 was found to perturb the composition of the apo equilibrium ensemble and to bind weakly to multiple distinct c-Myc conformations. Comparison of the apo and holo equilibrium ensembles reveals that the c-Myc conformations binding 1 are already partially formed in the apo ensemble, suggesting that 1 binds to c-Myc through an extended conformational selection mechanism. The present results have important implications for rational ligand design efforts targeting intrinsically disordered proteins

Analyzing and Biasing Simulations with PLUMED

This chapter discusses how the PLUMED plugin for molecular dynamics can be used to analyze and bias molecular dynamics trajectories. The chapter begins by introducing the notion of a collective variable and by then explaining how the free energy can be computed as a function of one or more collective variables. A number of practical issues mostly around periodic boundary conditions that arise when these types of calculations are performed using PLUMED are then discussed. Later parts of the chapter discuss how PLUMED can be used to perform enhanced sampling simulations that introduce simulation biases or multiple replicas of the system and Monte Carlo exchanges between these replicas. This section is then followed by a discussion on how free-energy surfaces and associated error bars can be extracted from such simulations by using weighted histogram and block averaging techniques