A Monte Carlo study of the dynamics of G-protein activation

To link quantitatively the cell surface binding of ligand to receptor with the production of cellular responses, it may be necessary to explore early events in signal transduction such as G-protein activation. Two different model frameworks relating receptor/ligand binding to G-protein activation are examined. In the first framework, a simple ordinary differential equation model is used to describe receptor/ligand binding and G-protein activation. In the second framework, the events leading to G-protein activation are simulated using a dynamic Monte Carlo model. In both models, reactions between ligand-bound receptors and G-proteins are assumed to be diffusion-limited. The Monte Carlo model predicts two regimes of G-protein activation, depending upon whether the lifetime of a receptor/ligand complex is long or short compared with the time needed for diffusional encounters of complexes and G-proteins. When the lifetime of a complex is relatively short compared with the diffusion time, the movement of ligand among free receptors by binding and unbinding ("switching") significantly enhances G-protein activation. Receptor antagonists dramatically reduce G-protein activation and, thus, signal transduction in this case, and significant clustering of active G-proteins near receptor/ligand complexes results. The simple ordinary differential equation model poorly predicts G-protein activation for this situation. In the alternative case, when diffusion is relatively fast, ligand movement among receptors is less important and the simple ordinary differential equation model and Monte Carlo model results are similar. In this case, there is little clustering of active G-proteins near receptor/ligand complexes. Results also indicate that as the GTPase activity of the alpha-subunit decreases, the steady-state level of alpha-GTP increases, although temporal sensitivity is compromised.