In 1938, H. L. Krall found a necessary and sufficient condition for an orthogonal polynomial set {P(n)(x)}0infinity to satisfy a linear differential equation of the form SIGMA(i=0)N l(i)(x)y(i)(x) = lambda(n)y(x). Here the authors give a new simple proof of Krall's theorem as well as some other characterizations of such orthogonal polynomial sets based on the symmetrizability of the differential operator. In particular it is shown that such orthogonal polynomial sets are characterized by a certain Sobolev-type orthogonality, which generalizes Hahn's charaterization of classical orthogonal polynomials.