Inferring Pathways of Oxidative Folding from Prefolding Free Energy Landscapes of Disulfide-Rich Toxins

Short, cysteine-rich peptides can exist in stable or metastable structural ensembles due to the number of possible patterns of formation of their disulfide bonds. One interesting subset of this peptide group is the coonotoxins, which are produced by aquatic snails in the family Conidae. The µ conotoxins, which are antagonists and blockers of the voltage-gated sodium channel, exist in a folding spectrum: on one end of the spectrum are more hirudin-like folders, which form disulfide bonds and then reshuffle them, leading to an ensemble of kinetically trapped isomers–and on the other end are more BPTI-like folders–which form the native disulfide bonds one by one in a particular order, leading to a preponderance of conformations existing in a single stable state. In this article, we employ the composite diffusion map approach to study the unified free energy surface of pre-folding µ-conotoxin equilibrium. We identify the two most important nonlinear collective modes of the unified folding landscape and demonstrate that in the absence of their disulfides, the conotoxins can be thought of as largely disordered polymers. A small increase in the number of hydrophobic residues in the protein shifts the free energy landscape towards hydrophobically collapsed coil conformations responsible for cysteine proximity in hirudin-like folders, compared to semi-extended coil conformations with more distal cysteines in BPTI-like folders. Overall, this work sheds important light on the folding processes and free energy landscapes of cysteinerich peptides and demonstrates the extent to which sequence and length contribute to these landscapes.