Computational modeling-driven study of signaling pathway affecting cleft formation during submandibular salivary gland development