advanced usage of PDE toolbox - prescribing a list of points for the boundaries and assigning BCs

Hi folks, we have been back and forth with Matlab support for several weeks/months on this issue and I still don't have a solution or an indication if it is at all possible.

We are solving moving boundary problem using Matlab's pde toolbox. The algorithm is quite simple. Prescribe the geometry, mesh, solve, move the boundaries (update the co-ordinates of the boundary nodes), mesh, solve, update boundary, mesh, solve, ....

Unfortunately although simple to describe there are issues. We have different boundary conditions on different parts of the domain. Applying boundary condition using edge labels is not simple because for some reason Matlab changes the edge labels as the simulation proceeds. We end up running into a mess.

Here is what I would like to do. I would like to be able to give the mesher and solver a list of points for the boundaries, and I would like to tell solver which boundary condition to use on each edge. At the next step the co-ordinates of the boundary points would move but we would again just update which boundary condition to use on each edge and repeat. Is this possible?