Objective To investigate the flexibility and mobility of the Bacillus thuringiensis toxin Cry1 Aa. Methods The graph theory-based program Constraint Network Analysis and normal mode-based program NMsim were used to analyze the global and local flexibility indices as well as the fluctuation of individual residues in detail. Results The decrease in Cry1 Aa network rigidity with the increase of temperature was evident. Two phase transition points in which the Cry1 Aa structure lost rigidity during the thermal simulation were identified. Two rigid clusters were found in domains I and II. Weak spot...