Interplay of secondary structures and side-chain contacts in the denatured state of BBA1

The denatured state of a miniprotein BBA1 is studied under the native condition with the AMBER/Poisson-Boltzmann energy model and with the self-guided enhanced sampling technique. Forty independent trajectories are collected to sample the highly diversified denatured structures. Our simulation data show that the denatured BBA1 contains high percentage of native helix and native turn, but low percentage of native hairpin. Conditional population analysis indicates that the native helix formation and the native hairpin formation are not cooperative in the denatured state. Side-chain analysis shows that the native hydrophobic contacts are more preferred than the non-native hydrophobic contacts in the denatured BBA1. In contrast, the salt-bridge contacts are more or less nonspecific even if their populations are higher than those of hydrophobic contacts. Analysis of the trajectories shows that the native helix mostly initiates near the N terminus and propagates to the C terminus, and mostly forms from 3(10)-helix/turn to alpha helix. The same analysis shows that the native turn is important but not necessary in its formation in the denatured BBA1. In addition, the formations of the two strands in the native hairpin are rather asymmetric, demonstrating the likely influence of the protein environment. Energetic analysis shows that the native helix formation is largely driven by electrostatic interactions in denatured BBA1. Further, the native helix formation is associated with the breakup of non-native salt-bridge contacts and the accumulation of native salt-bridge contacts. However, the native hydrophobic contacts only show a small increase upon the native helix formation while the non-native hydrophobic contacts stay essentially the same, different from the evolution of hydrophobic contacts observed in an isolated helix folding.