In order to distinguish between regular and chaotic planetary orbits we apply a new technique called MEGNO in a wide neighbourhood of orbital parameters determined using standard two-body Keplerian fits for HD 12661, HD 38529, HD 37124 and HD 160691 planetary systems. We show that the currently announced orbital parameters place these systems in very different situations from the point of view of dynamical stability. While HD 38529 and HD 37124 are located within large stability zones in the phase space around their determined orbits, the preliminary orbits in HD 160691 are highly unstable. The orbital parameters of the HD 12661 planets are located in a border region between stable and unstable dynamical regimes, so while its currently determined orbital parameters produce stable regular orbits, a minor change within the margin of error of just one parameter may result in a chaotic dynamical system.

12 pages, 3 figures, accepted ApJ, revised version following the referee's report