Protease is an enzyme that catalysed the hydrolysis of protein into peptide. Application of protease in industry has been linked with cost effective substrates and complex of enzyme-substrate stability. Molecular docking approach has identified casein as a preference substrates. However, lack of data on casein mode of binding to protease and enzyme stability represents a limitation for its production and structural optimization. In this study, we have used a molecular dynamic (MD) to examine the stability of complex enzyme-substrate of protease from Bacillus lehensis G1. The 3D structure of protease (BleG1_1979) was docked with substrate casein using AutoDock Vina. Structural analysis of the substrate-binding cleft revealed a binding site of casein was predominantly at the hydrophobic region of BleG1_1979. The MD of complex BleG1_1979-casein was tested with two temperatures; 298 K and 310 K using GROMACS v5.1.4. MD simulation showed a stable behaviour of BleG1_1979 over the 20 ns simulation period. The molecular docking and MD simulation suggested that the production of protease from B. lehensis G1 by utilization of casein and the stability of complex protease-casein could be a potential application to generate a cost effective enzyme to be develop for industrial use.