Optimization of Molecular Dynamics Simulations of c-MYC(1-88)—An Intrinsically Disordered System

Many of the proteins involved in key cellular regulatory events contain extensive intrinsically disordered regions that are not readily amenable to conventional structure/function dissection. The oncoprotein c-MYC plays a key role in controlling cell proliferation and apoptosis and more than 70% of the primary sequence is disordered. Computational approaches that shed light on the range of secondary and tertiary structural conformations therefore provide the only realistic chance to study such proteins. Here, we describe the results of several tests of force fields and water models employed in molecular dynamics simulations for the N-terminal 88 amino acids of c-MYC. Comparisons of the simulation data with experimental secondary structure assignments obtained by NMR establish a particular implicit solvation approach as highly congruent. The results provide insights into the structural dynamics of c-MYC(1-88), which will be useful for guiding future experimental approaches. The protocols for trajectory analysis described here will be applicable for the analysis of a variety of computational simulations of intrinsically disordered proteins.