The nematocyst walls of Hydra are formed by proteins containing small cysteine-rich domains (CRDs) of ~25 amino acids. The first CRD of nematocyst outer all antigen (NW1) and the C-terminal CRD of minicollagen-1 (Mcol1C) contain six cysteines at identical sequence positions, however adopt different disulfide bonded structures. NW1 shows the disulfide connectivities C2-C14/C6-C19/C10-C18 and Mcol1C C2-C18/C6-C14/C10-C19. To analyze if both show structural preferences in the open, non-disulfide bonded form, which explain the formation of either disulfide connectivity pattern, molecular dynamics (MD) simulations at different temperatures were performed. NW1 maintained in the 100-ns MD simulations at 283 K a rather compact fold that is stabilized by specific hydrogen bonds. The Mcol1C structure fluctuated overall more, however stayed most of the time also rather compact. The analysis of the backbone Φ/ψ angles indicated different turn propensities for NW1 and Mcol1C, which mostly can be explained based on published data about the influence of different amino acid side chains on the local backbone conformation. Whereas a folded precursor mechanism may be considered for NW1, Mcol1C may fold according to the quasi-stochastic folding model involving disulfide bond reshuffling and conformational changes, locking the native disulfide conformations. The study further demonstrates the power of MD simulations to detect local structural preferences in rather dynamic systems such as the open, non-disulfide bonded forms of NW1 and Mcol1C, which complement published information from NMR backbone residual dipolar couplings. Because the backbone structural preferences encoded by the amino acid sequence embedding the cysteines influence which disulfide connectivities are formed, the data are generally interesting for a better understanding of oxidative folding and the design of disulfide stabilized therapeutics.
Read full abstract