Proteins in vacuo: A molecular dynamics study of the unfolding behavior of highly charged disulfide-bond-intact lysozyme subjected to a temperature pulse

Molecular dynamics simulations were used to interpret a variety of experimental data on highly charged disulfide-bond-intact lysozyme in vacuo. The simulation approach involved submitting a model of the protein [Reimann, Velázquez, and Tapia, J. Phys. Chem. B $102,$ 9344 (1998)] in a given charge state to a 3-ns-long heat pulse (usually at 500 K) followed by cooling or relaxation for 1 ns back to room temperature (293 K). This treatment yielded a charge threshold around $Q_0=8+$ for obtaining significant unfolding, as indicated by an enhancement in collision cross section and conformer length. The collision cross sections and lengths theoretically obtained, along with the threshold charge state for initiating unfolding, were compatible with experimental results on lysozyme in vacuo. The unfolded, highly elongated conformations obtained for $Q>~9+$ displayed a significant level of non-native $\beta$-sheet content which appeared to be additionally stabilized by charge self-solvation.