A stable molecular structure is important in the development of a protein candidate into a therapeutic product. A therapeutic protein often contains many different variants; some of them may have an impact on the conformational stability of the protein. Conventionally, to evaluate the impact of a variant on stability, the variant must be enriched to a reasonable purity, and then its stability characterized by chromatographic or biophysical techniques. However, it is often impractical to purify and characterize each variant in a therapeutic protein. A workflow, based on limited proteolysis followed by MS detection, was established to simultaneously assess the impact of a large number of variants on conformational stability without enrichment. Because a less stable domain is more susceptible to proteolytic degradation, conformational stability of the domain can be reported from the release rate of a proteolytic peptide. A kinetic model is established to quantitatively determine the extent of domain stabilization/destabilization of different variants. The methodology is demonstrated by examining variants known to affect the stability of immunoglobulin domains, such as different N-glycoforms, methionine oxidations, and sequence variants. With this methodology, near 100 variants may be evaluated within 2 days in a single experiment. Insights into the sequence−stability relationship will be obtained by monitoring the large number of low-level sequence variants, facilitating engineering of more stable molecules.