Interactions Between Amino Acid Side Chains in Cylindrical Hydrophobic Nanopores with Applications to Peptide Stability

Confinement effects on protein stability are relevant in a number of biological applications ranging from encapsulation in the cylindrical cavity of a chaperonin, translocation through pores, and structure formation in the exit tunnel of the ribosome. Consequently, free energies of interaction between amino acid side chains in restricted spaces can provide insights into factors that control protein stability in nanopores. Using all-atom molecular dynamics simulations, we show that 3 pair interactions between side chains--hydrophobic (Ala-Phe), polar (Ser-Asn) and charged (Lys-Glu)--are substantially altered in hydrophobic, water-filled nanopores, relative to bulk water. When the pore holds water at bulk density, the hydrophobic pair is strongly destabilized and is driven to large separations corresponding to the width and the length of the cylindrical pore. As the water density is reduced, the preference of Ala and Phe to be at the boundary decreases, and the contact pair is preferred. A model that accounts for the volume accessible to Phe and Ala in the solvent-depleted region near the pore boundary explains the simulation results. In the pore, the hydrogen-bonded interactions between Ser and Asn have an enhanced dependence on their relative orientations, as compared with bulk water. When the side chains of Lys and Glu are restrained to be side by side, parallel to each other, then salt bridge formation is promoted in the nanopore. Based on these results, we argue and demonstrate that for a generic amphiphilic sequence, cylindrical confinement is likely to enhance thermodynamic stability relative to the bulk.