We present a technique for the direct optimization of conservative backbone curves in nonlinear mechanical systems. The periodic orbits on the conservative backbone are computed analytically using the reduced dynamics of the corresponding Lyapunov subcenter manifold (LSM). In this manner, we avoid expensive full-system simulations and numerical continuation to approximate the nonlinear response. Our method aims at tailoring the shape of the backbone curve using a gradient-based optimization with respect to the system’s parameters. To this end, we formulate the optimization problem by imposing constraints on the frequency-amplitude relation. Sensitivities are computed analytically by differentiating the backbone expression and the corresponding LSM. At each iteration, only the reduced-order model construction and sensitivity computation are performed, making our approach robust and efficient.