Simulations of reversible protein aggregate and crystal structure.

We simulated the structure of reversible protein aggregates as a function of protein surface characteristics, protein-protein interaction energies, and the entropic penalty accompanying the immobilization of protein in a solid phase. These simulations represent an extension of our previous work on kinetically irreversible protein aggregate structure and are based on an explicit accounting of the specific protein-protein interactions that occur within reversible aggregates and crystals. We considered protein monomers with a mixture of hydrophobic and hydrophilic surface regions suspended in a polar solvent; the energetic driving force for aggregation is provided by the burial of solvent-exposed hydrophobic surface area. We analyzed the physical properties of the generated aggregates, including density, protein-protein contact distributions, solvent accessible surface area, porosity, and order, and compared our results with the protein crystallization literature as well as with the kinetically irreversible case. The physical properties of reversible aggregates were consonant with those observed for the irreversible aggregates, although in general, reversible aggregates were more stable energetically and were more crystal-like in their order content than their irreversible counterparts. The reversible aggregates were less dense than the irreversible aggregates, indicating that the increased energetic stability is derived primarily from the optimality rather than the density of the packing in the solid phase. The extent of hydrophobic protein-protein contacts and solvent-exposed surface area within the aggregate phase depended on the aggregation pathway: reversible aggregates tended to have a greater proportion of hydrophobic-hydrophobic contacts and a smaller fraction of hydrophobic solvent-exposed surface area. Furthermore, the arrangement of hydrophobic patches on the protein surface played a major role in the distribution of protein contacts and solvent content. This was readily reflected in the order of the aggregates: the greater the contiguity of the hydrophobic patches on the monomer surface, the less ordered the aggregates became, despite the opportunities for rearrangement offered by a reversible pathway. These simulations have enhanced our understanding of the impact of protein structural motifs on aggregate properties and on the demarcation between aggregation and crystallization.