Universal effects of solvent species on the stabilized structure of a protein

We investigate the effects of solvent specificities on the stability of the native structure (NS) of a protein on the basis of our free-energy function (FEF). We use CPB-bromodomain (CBP-BD) and apoplastocyanin (apoPC) as representatives of the protein universe and water, methanol, ethanol, and cyclohexane as solvents. The NSs of CBP-BD and apoPC consist of 66% α-helices and of 35% β-sheets and 4% α-helices, respectively. In order to assess the structural stability of a given protein immersed in each solvent, we contrast the FEF of its NS against that of a number of artificially created, misfolded decoys possessing the same amino-acid sequence but significantly different topology and α-helix and β-sheet contents. In the FEF, we compute the solvation entropy using the morphometric approach combined with the integral equation theories, and the change in electrostatic (ES) energy upon the folding is obtained by an explicit atomistic but simplified calculation. The ES energy change is represented by the break of protein-solvent hydrogen bonds (HBs), formation of protein intramolecular HBs, and recovery of solvent-solvent HBs. Protein-solvent and solvent-solvent HBs are absent in cyclohexane. We are thus able to separately evaluate the contributions to the structural stability from the entropic and energetic components. We find that for both CBP-BD and apoPC, the energetic component dominates in methanol, ethanol, and cyclohexane, with the most stable structures in these solvents sharing the same characteristics described as an association of α-helices. In particular, those in the two alcohols are identical. In water, the entropic component is as strong as or even stronger than the energetic one, with a large gain of translational, configurational entropy of water becoming crucially important so that the relative contents of α-helix and β-sheet and the content of total secondary structures are carefully selected to achieve sufficiently close packing of side chains. If the energetic component is excluded for a protein in water, the priority is given to closest side-chain packing, giving rise to the formation of a structure with very low α-helix and β-sheet contents. Our analysis, which requires minimal computational effort, can be applied to any protein immersed in any solvent and provides robust predictions that are quite consistent with the experimental observations for proteins in different solvent environments, thus paving the way toward a more detailed understanding of the folding process.