  We propose to combine a mean field approach with all atom molecular dynamics
into a multistage algorithm that can model protein folding and dynamics over
very long time periods yet with atomic level precision. As an example we
investigate an isolated monomeric Myc oncoprotein that has been implicated in
carcinomas including those in colon, breast and lungs. Under physiological
conditions a monomeric Myc is presumed to be an example of intrinsically
disordered proteins, that pose a serious challenge to existing modelling
techniques. We argue that a room temperature monomeric Myc is in a dynamical
state, it oscillates between different conformations that we identify. For this
we adopt the C-alpha backbone of Myc in a crystallographic heteromer as an
initial Ansatz for the monomeric structure. We construct a multisoliton of the
pertinent Landau free energy, to describe the C-alpha profile with ultra high
precision. We use Glauber dynamics to resolve how the multisoliton responds to
repeated increases and decreases in ambient temperature. We confirm that the
initial structure is unstable in isolation. We reveal a highly degenerate
ground state landscape, an attractive set towards which Glauber dynamics
converges in the limit of vanishing ambient temperature. We analyse the thermal
stability of this Glauber attractor using room temperature molecular dynamics.
We identify and scrutinise a particularly stable subset in which the two
helical segments of the original multisoliton align in parallel, next to each
other. During the MD time evolution of a representative structure from this
subset, we observe intermittent quasiparticle oscillations along the C-terminal
alpha-helix, some of which resemble a translating Davydov's Amide-I soliton. We
propose that the presence of oscillatory motion is in line with the expected
intrinsically disordered character of Myc.
