We present a parametrization of the chiral even generalized parton distributions, $H$, $E$, $\stackrel{\texttildelow{}}{H}$, $\stackrel{\texttildelow{}}{E}$, for the quark, antiquark, and gluon, in the perturbative QCD (pQCD)-parton framework. Parametric analytic forms are given as a function of two equivalent sets of variables $x$, $\ensuremath{\xi}$, $t$ (symmetric frame) and $X$, $\ensuremath{\zeta}$, $t$ (asymmetric frame), at an initial scale, ${Q}_{o}^{2}$. In the $X>\ensuremath{\zeta}$ region, a convenient and flexible form is obtained as the product of a Regge term $\ensuremath{\propto}{X}^{\ensuremath{-}\ensuremath{\alpha}+{\ensuremath{\alpha}}^{\ensuremath{'}}t}$, describing the low $X$ behavior, times a spectator model-based functional form depending on various mass parameters; the behavior at $X<\ensuremath{\zeta}$ is determined using the generalized parton distributions symmetry and polynomiality properties. The parameters are constrained using data on the flavor separated nucleon electromagnetic elastic form factors, the axial and pseudoscalar nucleon form factors, and the parton distribution functions from both the deep inelastic unpolarized and polarized nucleon structure functions. For the gluon distributions we use, in particular, constraints provided by recent lattice QCD moments calculations. The parametrization's kinematical range of validity is $0.0001\ensuremath{\le}X\ensuremath{\le}\phantom{\rule{0ex}{0ex}}0.85$, $0.01\ensuremath{\le}\ensuremath{\zeta}\ensuremath{\le}0.85$, $0\ensuremath{\le}\ensuremath{-}t\ensuremath{\le}1\text{ }\text{ }{\mathrm{GeV}}^{2}$, $2\ensuremath{\le}{Q}^{2}\ensuremath{\le}100\text{ }\text{ }{\mathrm{GeV}}^{2}$. With the simultaneous description of the quark, antiquark, and gluon sectors, this parametrization represents a first tool enabling a global QCD analysis of deeply virtual exclusive experiments.