  Over the last decades, there has been a tremendous increase in research on
extrasolar planets. Many exosolar systems, which consist of a Star and two
inclined Planets, seem to be locked in 4/3, 3/2, 2/1, 5/2, 3/1 and 4/1 mean
motion resonance (MMR). We herewith present the model used to simulate three
dimensional planetary systems and provide planar families of periodic orbits
(PO), which belong to all possible configurations that each MMR has, along with
their linear horizontal and vertical stability. We focus on depicting stable
spatial families (most of them up to mutual inclination of $60^\circ$)
generated by PO of planar circular families, because the trapping in MMR could
be a consequence of planetary migration process. We attempt to connect the
linear stability of PO with long-term stability of a planetary system close to
them. This can stimulate the search of real planetary systems in the vicinity
of stable spatial PO-counterbalanced by the planets' orbital elements, masses
and MMR; all of which could constitute a suitable environment convenient to
host them.
