We show how generalized parton distributions (GPDs) can be determined in the case where hadrons are described in terms of their partonic degrees of freedom through solutions of dynamical equations. We demonstrate our approach on the example of two-quark bound states described by the Bethe-Salpeter equation, and three-quark bound states described by three-and four-dimensional Faddeev-like equations. Within the model of strong interactions defined by the dynamical equations, all possible mechanisms contributing to the GPDs are taken into account, and all GPD sum rules are satisfied automatically. The formulation is general and can be applied to determine generalized quark distributions, generalized gluon distributions, transition GPDs, nucleon distributions in nuclei, etc. Our approach is based on the gauging of equations method.