Molecular interactions between mesenchymal-derived Keratinocyte growth factor (KGF) and Kit ligand (KITLG) are essential for follicular development. These factors are expressed by theca and granulosa cells. We determined full length coding sequence of buffalo KGF and KITLG proteins having 194 and 274 amino acids, respectively. The recombinant KGF and KITLG proteins were solubilized in 10 mM Tris, pH 7.5 and 50 mM Tris, pH 7.4 and purified using Ni-NTA column and GST affinity chromatography, respectively. The purity and molecular weight of His-KGF (~23 kDa) and GST-KITLG (~57 kDa) proteins were confirmed by SDS-PAGE and western blotting. The co-immunoprecipitation assay accompanied with computational analysis demonstrated the interaction between KGF and KITLG proteins. We deduced 3D structures of the candidate proteins and assessed their binding based on protein docking. In the process, KGF specific residues, Lys123, Glu135, Lys140, Lys155 and Trp156 and KITLG specific ones, Ser226, Phe233, Gly234, Ala235, Phe236, Trp238 and Lys239 involved in the formation of KGF-KITLG complex were detected. The hydrophobic interactions surrounding KGF-KITLG complex affirmed their binding affinity and stability to the interacting interface. Additionally, in-silico site directed mutagenesis enabled the assessment of changes that occurred in the binding energies of mutated KGF-KITLG protein complex. Our results demonstrate that in the presence of KITLG, KGF mimics its native binding mode suggesting all the KGF residues are specific to their binding complex. This study provides an insight on the critical amino acid residues participating in buffalo ovarian folliculogenesis.