Developing efficient and robust computational models is essential to improve our understanding of protein solution behavior. This becomes particularly important to tackle the high-concentration regime. In this context, the main challenge is to put forward coarse-grained descriptions able to reduce the level of detail, while retaining key features and relevant information. In this work, we develop an efficient strategy that can be used to investigate and gain insight into monoclonal antibody solutions under different conditions. We use a multi-scale numerical approach, which connects information obtained at all-atom and amino-acid levels to bead models. The latter has the advantage of reproducing the properties of interest while being computationally much faster. Indeed, these models allow us to perform many-protein simulations with a large number of molecules. We can, thus, explore conditions not easily accessible with more detailed descriptions, perform effective comparisons with experimental data up to very high protein concentrations, and efficiently investigate protein–protein interactions and their role in phase behavior and protein self-assembly. Here, a particular emphasis is given to the effects of charges at different ionic strengths.