Molecular dynamics simulation to investigate the impact of disulfide bond formation on conformational stability of chicken cystatin I66Q mutant

Abstract

Chicken cystatin (cC) mutant I66Q is located in the hydrophobic core of the protein and increases the propensity for amyloid formation. Here, we demonstrate that under physiological conditions, the replacement of Ile with the Gln in the I66Q mutant increases the susceptibility for the disulfide bond Cys71–Cys81 to be reduced when compared to the wild type (WT) cC. Molecular dynamics (MD) simulations under conditions favoring cC amyloid fibril formation are in agreement with the experimental results. MD simulations were also performed to investigate the impact of disrupting the Cys71–Cys81 disulfide bond on the conformational stability of cC at the atomic level, and highlighted major disruption to the cC appendant structure. Domain swapping and extensive unfolding has been proposed as one of the possible mechanisms initiating amyloid fibril formation by cystatin. Our in silico studies suggest that disulfide bond formation between residues Cys95 and Cys115 is necessary to maintain conformational stability of the I66Q mutant following breakage of the Cys71–Cys81 disulfide bridge. Subsequent breakage of disulfide bond Cys95–Cys115 resulted in large structural destabilization of the I66Q mutant, which increased the α–β interface distance and expanded the hydrophobic core. These experimental and computational studies provide molecular-level insight into the relationship between disulfide bond formation and progressive unfolding of amyloidogenic cC mutant I66Q. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:23