Theoretical study for regulatory property of scaffold protein on MAPK cascade: A qualitative modeling

Understand the Functions of Scaffold Proteins in Cell Signaling by a Mesoscopic Simulation Method

Scaffold proteins are central players in regulating the spatial-temporal organization of many important signaling pathways in cells. They offer physical platforms to downstream signaling proteins so that their transient interactions in a crowded and heterogeneous environment of cytosol can be greatly facilitated. However, most scaffold proteins tend to simultaneously bind more than one signaling molecule, which leads to the spatial assembly of multimeric protein complexes. The kinetics of these protein oligomerizations are difficult to quantify by traditional experimental approaches. To understand the functions of scaffold proteins in cell signaling, we developed a, to our knowledge, new hybrid simulation algorithm in which both spatial organization and binding kinetics of proteins were implemented. We applied this new technique to a simple network system that contains three molecules. One molecule in the network is a scaffold protein, whereas the other two are its binding targets in the downstream signaling pathway. Each of the three molecules in the system contains two binding motifs that can interact with each other and are connected by a flexible linker. By applying the new simulation method to the model, we show that the scaffold proteins will promote not only thermodynamics but also kinetics of cell signaling given the premise that the interaction between the two signaling molecules is transient. Moreover, by changing the flexibility of the linker between two binding motifs, our results suggest that the conformational fluctuations in a scaffold protein play a positive role in recruiting downstream signaling molecules. In summary, this study showcases the capability of computational simulation in understanding the general principles of scaffold protein functions.

Spatial distribution and dose-response relationship for different operation modes in a reaction-diffusion model of the MAPK cascade

The mitogen-activated protein kinase (MAPK) cascade plays a critical role in the control of cell growth. Deregulation of this pathway contributes to the development of many cancers. To better understand its signal transduction, we constructed a reaction-diffusion model for the MAPK pathway. We modeled the three layers of phosphorylation-dephosphorylation reactions and diffusion processes from the cell membrane to the nucleus. Based on different types of feedback in the MAPK cascade, four operation modes are introduced. For each of the four modes, spatial distributions and dose-response curves of active kinases (i.e. ppMAPK) are explored by numerical simulation. The effects of propagation length, diffusion coefficient and feedback strength on the pathway dynamics are investigated. We found that intrinsic bistability in the MAPK cascade can generate a traveling wave of ppMAPK with constant amplitude when the propagation length is short. ppMAPK in this mode of intrinsic bistability decays more slowly than it does in all other modes as the propagation length increases. Moreover, we examined the global and local responses to Ras-GTP of these four modes, and demonstrated how the shapes of these dose-response curves change as the propagation length increases. Also, we found that larger diffusion constant gives a higher response level on the zero-order regime and makes the ppMAPK profiles flatter under strong Ras-GTP stimulus. Furthermore, we observed that spatial responses of ppMAPK are more sensitive to negative feedback than to positive feedback in the broader signal range. Finally, we showed how oscillatory signals pass through the kinase cascade, and found that high frequency signals are damped faster than low frequency ones.