Abstract
Excipients are included within protein biotherapeutic solution formulations to improve colloidal and conformational stability but are generally not designed for the specific purpose of preventing aggregation and improving cryoprotection in solution. In this work, we have explored the relationship between the structure and antiaggregation activity of excipients by utilizing coarse-grained molecular dynamics modeling of protein-excipient interaction. We have studied human serum albumin as a model protein, and we report the interaction of 41 excipients (polysorbates, fatty alcohol ethoxylates, fatty acid ethoxylates, phospholipids, glucosides, amino acids, and others) in terms of the reduction of solvent accessible surface area of aggregation-prone regions, proposed as a mechanism of aggregation prevention. Polyoxyethylene sorbitan had the greatest degree of interaction with aggregation-prone regions, decreasing the solvent accessible surface area of APRs by 20.7 nm(2) (40.1%). Physicochemical descriptors generated by Mordred are employed to probe the structure-property relationship using partial least-squares regression. A leave-one-out cross-validated model had a root-mean-square error of prediction of 4.1 nm(2) and a mean relative error of prediction of 0.077. Generally, longer molecules with a large number of alcohol-terminated PEG units tended to interact more, with qualitatively different protein interactions, wrapping around the protein. Shorter or less ethoxylated compounds tend to form hemimicellar clusters at the protein surface. We propose that an improved design would feature many short chains of 5 to 10 PEG units in many distinct branches and at least some hydrophobic content in the form of medium-length or greater aliphatic chains (i.e., six or more carbon atoms). The combination of molecular dynamics simulation and quantitative modeling is an important first step in an all-purpose protein-independent model for the computer-aided design of stabilizing excipients.