Calcium calmodulin dependent kinase II (CaMKII) is sequestered in dendritic spines within seconds upon synaptic stimulation. The program Smoldyn was used to develop scenarios of single molecule CaMKII diffusion and binding in virtual dendritic spines. We first validated simulation of diffusion as a function of spine morphology. Additional cellular structures were then incorporated to simulate binding of CaMKII to the post-synaptic density (PSD); binding to cytoskeleton; or their self-aggregation. The distributions of GFP tagged native and mutant constructs in dissociated hippocampal neurons were measured to guide quantitative analysis. Intra-spine viscosity was estimated from fluorescence recovery after photo-bleach (FRAP) of red fluorescent protein. Intra-spine mobility of the GFP-CaMKII? constructs was measured, with hundred-millisecond or better time resolution, from FRAP of distal spine tips in conjunction with fluorescence loss (FLIP) from proximal regions. Different FRAP \ FLIP profiles were predicted from our Scenarios and provided a means to differentiate binding to the PSDs from self-aggregation. The predictions were validated by experiments. Simulated fits of the Scenarios provided estimates of binding and rate constants. We utilized these values to assess the role of self-aggregation during the initial response of native CaMKII holoenzymes to stimulation. The computations revealed that self-aggregation could provide a concentration-dependent switch to amplify CaMKII sequestration and regulate its activity depending on its occupancy of the actin cytoskeleton.