Stellarator optimization for nested magnetic surfaces at finite β and toroidal current