The mass shift parameters, the electron density at the origin, the hyperfine interaction constants and the gJ Land´e factors are computed along the beryllium, boron, carbon and nitrogen isoelectronic sequences. The calculations on B- and C-like ions are based on the wave functions described in [1-4]. The many-electron wave functions corresponding to the Be and N isoelectronic sequences are obtained using the new version of the grasp2K multiconfiguration Dirac-Fock package, following similar optimization strategies. A new program, hereafter referred as ris and designed as a module of grasp2K, calculates the mass shift parameters and the electron density at the origin within the relativistic framework. For estimating the mass shift, ris considers the expectation value of the following operator HMS = 1 2M N Xi , j pi · pj − αZ ri i + ( i · ri) ri r2 i · pj , that is more complete than the one calculated in sms92. The one-body part (i = j) of the first term of equation is responsible for the observed breakdown of the Dirac kinetic energy operator often used to evaluate the isotope normal mass shift. The second term of equation takes the nuclear recoil corrections into account at the (αZ) order. The programs ris and rhfs2 allow the storage of the angular coefficients to reduce the cpu time when calculations are performed along a given isoelectronic sequence.