Proteins aggregate and precipitate from high concentration solutions in a wide variety of problems of natural and technological interest. Consequently, there is a broad interest in developing new ways to model the thermodynamic and kinetic aspects of protein stability in these crowded cellular or solution environments. We use a coarse-grained modeling approach to study the effects of different crowding agents on the conformational equilibria of proteins and the thermodynamic phase behavior of their solutions. At low to moderate protein concentrations, we find that crowding species can either stabilize or destabilize the native state, depending on the strength of their attractive interaction with the proteins. At high protein concentrations, crowders tend to stabilize the native state due to excluded volume effects, irrespective of the strength of the crowder-protein attraction. Crowding agents reduce the tendency of protein solutions to undergo a liquid-liquid phase separation driven by strong protein-protein attractions. The aforementioned equilibrium trends represent, to our knowledge, the first simulation predictions for how the properties of crowding species impact the global thermodynamic stability of proteins and their solutions.