Modeling of laminar-turbulent transition remains one of the key challenges in the numerical simulation of boundary layers, especially at coarse grid resolutions. Transition to turbulence strongly affects practical concerns such as wall-shear stress distribution and surface heat transfer. For wall-modeled large-eddy simulations (WMLES) specifically, the thin laminar region can require 10-100 times more grid points than the turbulent region to properly capture the amplification of disturbances preceding breakdown to turbulence (Slotnick et al., 2014). We propose to use the nonlinear parabolized stability equations (NLPSE) in conjunction with WMLES to simulate transitional flows affordably and accurately.