Computational studies of proteins have significantly improved our understanding of protein folding. These studies are normally carried out by using chains in isolation. However, in many systems of practical interest, proteins fold in the presence of other molecules. To obtain insight into folding in such situations, we compare the thermodynamics of folding for a Miyazawa-Jernigan model 64-mer in isolation to results obtained in the presence of additional chains. The melting temperature falls as the chain concentration increases. In multichain systems, free-energy landscapes for folding show an increased preference for misfolded states. Misfolding is accompanied by an increase in interprotein interactions; however, near the folding temperature, the transition from folded chains to misfolded and associated chains is entropically driven. A majority of the most probable interprotein contacts are also native contacts, suggesting that native topology plays a role in early stages of aggregation.computer simulation ͉ protein aggregation ͉ lattice model L attice-model proteins have played a key role in developing our understanding of protein folding. These model proteins contain enough detail to capture the essential physics of the folding process, yet are amenable to rigorous calculation of free-energy landscapes used to describe the folding pathway. Such calculations have provided a conceptual solution to the Levinthal Paradox, which ponders the ability of a protein to navigate a vast amount of conformational space to reach its native state on time scales of seconds or less (1-3). The funnellike nature of the free-energy landscapes, first calculated from lattice models, shows that proteins only need to sample a small fraction of conformations to reach the native state. The energetic bias toward the native state exists because native interactions are, on average, more stable than nonnative ones. Similar features are observed in folding landscapes generated from experiments, validating results from model calculations (1).Most computational studies of protein folding have examined a single chain in isolation (4-7). However, in many systems of practical interest, including in vivo folding, proteins fold in crowded environments. In such situations, interactions with other biological molecules compete with the intraprotein interactions that bias a protein's conformation toward its native state. Some biological molecules, such as molecular chaperones, promote folding (8). However, intermolecular interactions can also induce misfolding and aggregation (9, 10), resulting in loss of protein function (11). Further, protein aggregates can be toxic. Protein association has been linked to Ͼ20 human diseases, including Alzheimer's, Parkinson's, and Huntington's diseases (12, 13).We report simulations for systems containing one-, two-, or four-lattice model 64-mers. This chain length is greater than that in most model studies of multichain systems and correspondingly provides a more realistic surface area͞volume ratio than that for sm...