From matrix scaling to Brascamp-Lieb and back