Interactions of anionic, cationic and metal phthalocyanine with adenosine deaminase were studied by molecular dynamics and docking simulation. Structural parameters such as solvent accessible surface area (SAS), mid-point of transition temperature (T m ), radial distribution function (RDF) and hydrogen bond, helix, coil, beta percentage and other physical parameters were obtained. The denaturation of adenosine deaminase (ADA) by heat, anionic and cationic phthalocyanines was compared. A series of 20 ns simulation performed at temperatures ranging from 275 to 450 K, starting from the ADA native structure. Results of radial distribution functions (RDFs) showed that metallic derivative at low concentration behaves the same as osmolytes that increases the beta form and increases the enzyme stability. Molecular docking studies have been carried out to confirm the simulation results. Investigation of binding site and free energy confirmed that the efficiency of interaction with adenosine deaminase depends on metal core. Binding energy of non-metallic form is more negative than metallic form and it significantly decreases for phthalocyanine. Self-aggregation of anionic phthalocyanine decreases in comparison with cationic derivative, therefore enzyme denaturation in the presence of anionic form is higher than the other. Furthermore, thermal stability of the enzyme also depends on temperature in presence of phthalocyanine. Binding site of phthalocyanine on the enzyme has been identified by docking analysis.