The prediction of neutron stars properties is strictly connected to the employed nuclear interactions. The appearance of hyperons in the inner core of the star is strongly dependent on the details of the underlying hypernuclear force. We summarize our recent quantum Monte Carlo results on the development of realistic two- and three-body hyperon-nucleon interactions based on the available experimental data for light- and medium-heavy hypernuclei.