Interior-boundary conditions (IBCs) have been suggested as a possibility to circumvent the problem of ultraviolet divergences in quantum field theories. In the IBC approach, particle creation and annihilation is described with the help of linear conditions that relate the wave functions of two sectors of Fock space: ψ (n) (p) at an interior point p and ψ (n+m) (q) at a boundary point q, typically a collision configuration. Here, we extend IBCs to the relativistic case. To do this, we make use of Dirac's concept of multi-time wave functions, i.e., wave functions ψ(x 1 , ..., x N ) depending on N space-time coordinates x i for N particles. This provides the manifestly covariant particle-position representation that is required in the IBC approach. In order to obtain rigorous results, we construct a model for Dirac particles in 1+1 dimensions that can create or annihilate each other when they meet. Our main results are an existence and uniqueness theorem for that model, and the identification of a class of IBCs ensuring local probability conservation on all Cauchy surfaces. Furthermore, we explain how these IBCs relate to the usual formulation with creation and annihilation operators. The Lorentz invariance is discussed and it is found that apart from a constant matrix (which is required to transform in a certain way) the model is manifestly Lorentz invariant. This makes it clear that the IBC approach can be made compatible with relativity.