We propose an affine term structure model which accommodates nonlinearities in the drift and volatility function of the short-term interest rate. Such nonlinearities are a consequence of discrete beta-distributed regime shifts constructed on multiple thresholds. We derive iterative closed-form formula for the whole yield curve dynamics that can be estimated using a linearized Kalman filter. Fitting the model on US data, we collect empirical evidence of its potential in estimating conditional volatility and correlation across yields.