Generation of the configurational ensemble of an intrinsically disordered protein from unbiased molecular dynamics simulation

Abstract

Significance A major challenge in biology is characterizing the structural flexibility of intrinsically disordered proteins (IDPs). Ensemble-averaged experimental data do not provide the underlying protein structures. Here, we performed independently small-angle neutron and X-ray scattering experiments and unbiased molecular dynamics simulations to probe the solution structure of an IDP. We report that enhancing the sampling of the simulations can generate an ensemble of IDP structures in quantitative agreement with scattering and NMR, without the need for biasing the simulation or reweighting the results. The demonstration of established simulation technology that produces accurate physical models of flexible biosystems may pave the way to relating conformational flexibility to biological function. Intrinsically disordered proteins (IDPs) are abundant in eukaryotic proteomes, play a major role in cell signaling, and are associated with human diseases. To understand IDP function it is critical to determine their configurational ensemble, i.e., the collection of 3-dimensional structures they adopt, and this remains an immense challenge in structural biology. Attempts to determine this ensemble computationally have been hitherto hampered by the necessity of reweighting molecular dynamics (MD) results or biasing simulation in order to match ensemble-averaged experimental observables, operations that reduce the precision of the generated model because different structural ensembles may yield the same experimental observable. Here, by employing enhanced sampling MD we reproduce the experimental small-angle neutron and X-ray scattering profiles and the NMR chemical shifts of the disordered N terminal (SH4UD) of c-Src kinase without reweighting or constraining the simulations. The unbiased simulation results reveal a weakly funneled and rugged free energy landscape of SH4UD, which gives rise to a heterogeneous ensemble of structures that cannot be described by simple polymer theory. SH4UD adopts transient helices, which are found away from known phosphorylation sites and could play a key role in the stabilization of structural regions necessary for phosphorylation. Our findings indicate that adequately sampled molecular simulations can be performed to provide accurate physical models of flexible biosystems, thus rationalizing their biological function.