The enzyme phenylalanine ammonia lyase (PAL) is of special importance for the treatment of phenylketonuria patients. The aim of this study was to find a stable recombinant PAL with suitable kinetic properties among all natural PAL producing species using in silica and experimental approaches. To find such a stable PAL among 481 natural isoforms, 48,000 of 3-D models were predicted using the Modeller 9.10 program and evaluated by Ramachandran plot. Correlation analysis between Ramachandran plot and the energy of different thermodynamic components indicated that this plot could be an appropriate tool to predict protein stability. Hence, PAL6 from Lotus japonicus (LjPAL6) was selected as a stable isoform. Molecular dynamic (MD) simulation for 50 ns and docking has been conducted for LjPAL6-phenylalanine complex. The best PAL-phenylalanine frame was selected by re-docking with L-phenylalanine (L.Phe) and root-mean-square deviation (RMSD) value. MD simulation showed that the complex has a good stability, depicted by the low RMSD value, binding free energy and hydrogen bindings. Docking results showed that LjPAL6 has a high affinity toward L-Phe according to the low level of binding free energy. By overexpressing Ljpal6 in E. coli BL21, a total of 33.5 mg/l of protein was obtained, which has been increased to 83.7 mg/l via the optimization of LjPAL6 production using response surface methodology. The optimal pH and temperature were 8.5 and 50 degrees C, respectively. LjPAL6 showed a specific activity of 42 nkat/mg protein, with K-m,K-cat and K-cat/K-m values of 0.483 mM, 7 S-1 and 14.5 S-1 mM(-1) for L-phe, respectively. In conclusion, finding models with the most reasonable stereo-chemical quality and lowest numbers of steric clashes would result in easier folding. Hence, in silica analyses of bulk data from natural origin will lead one to find an optimal model for in vitro studies and drug design.