Computational prediction of small molecules with predicted binding to FGFR3 and testing biological effects in bone cells

Activating anabolic receptor-mediated signaling is essential for stimulating new bone formation and for promoting bone healing in humans. Fibroblast growth factor receptor (FGFR) 3 is reported to be an important positive regulator of osteogenesis. Presently, recombinant proteins are used to stimulate FGFR3 function but have limitations for therapy due to expense and stability. Therefore, there is a need for identification of novel small molecules binding to FGFR3 that promote biological function. In silico molecular docking and high-throughput virtual screening on zinc database identified seven compounds predicted to bind to an active site within the βCʹ-βE loop, specific to FGFR3. All seven compounds fall within an acceptable range of ADME/T properties. Four compounds showed a 30–65% oral absorption rate. Density functional theory analysis revealed a high HOMO-LUMO gap, reflecting high molecular stability for compounds 14977614 and 13509082. Five compounds exhibited mutagenicity, while the other three compounds presented irritability. Computational mutagenesis predicted that mutating G322 affected compound binding to FGFR3. Molecular dynamics simulation revealed compound 14977614 is stable in binding to FGFR3. Furthermore, compound 14977614, with an oral absorption rate of 60% and high molecular stability, produced significant increases in both proliferation and differentiation of bone marrow stromal cells in vitro. Anti-FGFR3 treatment completely blocked the stimulatory effect of 14977614 on BMSC proliferation. Ex vivo treatment of mouse calvaria in organ culture for seven days with 14977614 increased mineralization and expression levels of bone formation markers. In conclusion, computational analyses identified seven compounds that bind to the FGFR3, and in vitro studies showed that compound 14977614 exerts significant biological effects on osteogenic cells.