This paper presents a method for multidisciplinary design optimization of offshore wind turbines at system level. The formulation and implementation that enable the integrated aerodynamic and structural design of the rotor and tower simultaneously are detailed. The objective function to be minimized is the levelized cost of energy. The model includes various design constraints: stresses, deflections, modal frequencies and fatigue limits along different stations of the blade and tower. The rotor design variables are: chord and twist distribution, blade length, rated rotational speed and structural thicknesses along the span. The tower design variables are: tower thickness and diameter distribution, as well as the tower height. For the other wind turbine components, a representative mass model is used to include their dynamic interactions in the system. To calculate the system costs, representative cost models of a wind turbine located in an offshore wind farm are used. To show the potential of the method and to verify its usefulness, the 5 MW NREL wind turbine is used as a case study. The result of the design optimization process shows 2.3% decrease in the levelized cost of energy for a representative Dutch site, while satisfying all the design constraints. (C) 2014 Elsevier Ltd. All rights reserved.