This paper proposes an approach for optimal design of multilevel systems under uncertainty. The approach utilizes the stochastic extension of the analytical target cascading formulation. The reliability of satisfying the probabilistic constraints is computed by means of the most probable point method using the hybrid mean value algorithm. A linearization technique is employed for estimating the propagation of uncertainties throughout the problem hierarchy. The proposed methodology is applied to a piston-ring/cylinder-liner engine subassembly design problem. Specifically, we assess the impact of variations in manufacturing-related properties such as surface roughness on engine attributes such as brake-specific fuel consumption. Results are compared to the ones obtained using Monte Carlo simulation.