We investigate the generation of features at large scales in the primordial power spectrum (PPS) when inflation is driven by two scalar fields. In canonical single field models of inflation, these features are often generated due to deviations from the slow-roll regime. While deviations from slow-roll can be naturally achieved in two field models due to a sharp turn in the trajectory in the field space, features at the largest scales of the types suggested by CMB temperature anisotropies are more difficult to achieve in models involving two canonical scalar fields due to the presence of isocurvature fluctuations. We show instead that a coupling between the kinetic terms of the scalar fields can easily produce such features. We discuss models whose theoretical predictions are consistent with current observations and highlight the implications of our results. © 2020 IOP Publishing Ltd and Sissa Medialab.