Molecular Simulations of a Dynamic Protein Complex: Role of Salt-Bridges and Polar Interactions in Configurational Transitions

Ion charge pairs and hydrogen bonds have been extensively studied for their roles in stabilizing protein complexes and in steering the process of protein association. Recently, it has become clear that some protein complexes are dynamic in that they interconvert between several alternate configurations. We have previously characterized one such system: the EphA2:SHIP2 SAM-SAM heterodimer by solution NMR. Here we carried out extensive all-atom molecular-dynamics simulations on a microsecond time-scale starting with different NMR-derived structures for the complex. Transitions are observed between several discernible configurations at average time intervals of 50–100 ns. The domains reorient relative to one another by substantial rotation and a slight shifting of the interfaces. Bifurcated and intermediary salt-bridge and hydrogen-bond interactions play a role in the transitions in a process that can be described as moving along a “monkey-bar”. We notice an increased density of salt bridges near protein interaction surfaces that appear to enable these transitions, also suggesting why the trajectories can become kinetically hindered in regions where fewer of such interactions are possible. In this context, even microsecond molecular-dynamics simulations are not sufficient to sample the energy landscape unless the structures remain close to their experimentally derived low-energy configurations.