The Maxwell–Schrödinger equations describing the ultrashort laser pulse propagation in multilevel molecular media are shown to be analytically solvable and to reduce to nonlinear multisine-Gordon equations for certain models of the field–molecule system. Special cases of level schemes are shown to lead to soliton propagation. A numerical study is presented in order to examine the effect of nonresonant (virtual) transitions on pulse propagation and amplification in a three-level molecular system as a function of intensity and detuning.