A soft ellipsoid contact pair potential model for biaxial molecules is proposed considering the configuration dependent energy anisotropy explicitly alongwith their geometrical aspects. Analytic expressions of forces and torques for this coarse-grained potential are provided to make it convenient for molecular dynamics simulation study. Additionally, some molecular dynamics simulations are performed to generate biaxial smectic and nematic liquid crystal phases.