One of the advantages of the chemically meaningful and reliable Z-matrix is the possibility to interpolate internal coordinates (IC) in an automated fashion. Using the same construction table for both minimum structures, their Z-matrices can be interpolated just as easily as their cartesian coordinates and if the Zmatrices were defined in a chemically meaningful way, the interpolated structures will be usually much more chemically reasonable. We demonstrate the superior automated interpolation by generating starting guesses for the Nudged Elastic Band (NEB) method. NEB approximates minimum energy paths (MEPs) by connecting minimum structures as a string of discrete intermediate structures ("images"). At the time of writing, most NEB implementations generate the initial pathway by interpolation in cartesian coordinates. These interpolations will often generate structures with compressed bonds or overlapping atoms when bond angles change or molecule fragments rotate against each other. There has been development of methods that obviate this problem, such as the Image-Dependent Pair Potential (IDPP), which performs a pre-optimization of the initial cartesian interpolation, reducing atom overlaps and bond deformations. In this work we demonstrate that starting guesses from interpolation in internal coordinates, that use black-box-like automatically generated Z-matrices, performs similar or better compared to cartesian interpolation + IDPP.