Thermal stability limits of proteins in solution and adsorbed on a hydrophobic surface

A coarse-grained Monte Carlo simulation is used to study thermal denaturation of small proteins in an infinitely dilute solution and adsorbed on a flat hydrophobic surface. Intermolecular interactions are modeled using the Miyazawa-Jernigan (MJ) knowledge-based potential for implicit solvent with the BULDG hydrophobicity scale. We analyze the thermal behavior of lysozyme for its prevalence of alpha-helices, fibronectin for its prevalence of beta-sheets, and a short single helical peptide. Protein dimensions and contact maps are studied in detail before and during isothermal adsorption and heating. The MJ potential is shown to correctly predict the native conformation in solution under standard conditions, and the anticipated thermal stabilization of adsorbed proteins is observed when compared with heating in solution. The helix of the peptide is found to be much less stable thermally than the helices of lysozyme, reinforcing the importance of long-range forces in defining the protein structure. Contact map analysis of the adsorbed proteins shows correlation between the hydrophobicity of the secondary structure and their thermal stability on the surface.