In the last decades, treatment planning for multicatheter interstitial breast brachytherapy has evolved considerably from fluoroscopy-based 2D to anatomy-based 3D planning. To plan the right positions of the catheters, ultrasound or computed tomography (CT) imaging can be used, but the treatment plan is always based on postimplant CT images. With CT imaging, the 3D target volume can be defined more precisely and delineation of the organs at risk volumes is also possible. Consequently, parameters calculated from dose-volume histogram can be used for quantitative plan evaluation. The catheter reconstruction is also easier and faster on CT images compared to X-ray films. In high dose rate brachytherapy, using a stepping source, a number of forward dose optimization methods (manual, geometrical, on dose points, graphical) are available to shape the dose distribution to the target volume, and these influence dose homogeneities to different extent. Currently, inverse optimization algorithms offer new possibilities to improve dose distributions further considering the requirements for dose coverage, dose homogeneity, and dose to organs at risk simultaneously and automatically. In this article, the evolvement of treatment planning for interstitial breast implants is reviewed, different forward optimization methods are discussed, and dose-volume parameters used for quantitative plan evaluation are described. Finally, some questions of the inverse optimization method are investigated and initial experiences of the authors are presented.