Hybrid Bernstein Normalizing Flows for Flexible Multivariate Density Regression with Interpretable Marginals