Basith S, Manavalan B, Lee G. Amyotrophic lateral sclerosis disease-related mutations disrupt the dimerization of superoxide dismutase 1 - A comparative molecular dynamics simulation study.
Comput Biol Med 2022;
151:106319. [PMID:
36446187 DOI:
10.1016/j.compbiomed.2022.106319]
[Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2022] [Revised: 10/31/2022] [Accepted: 11/13/2022] [Indexed: 11/27/2022]
Abstract
More than 150 genes are involved in amyotrophic lateral sclerosis (ALS), with superoxide dismutase 1 (SOD1) being one of the most studied. Mutations in SOD1 gene, which encodes the enzyme SOD1 is the second most prevalent and studied cause of familial ALS. SOD1 is a ubiquitous, homodimeric metalloenzyme that forms a critical component of the cellular defense against reactive oxygen species. Several mutations in the SOD1 enzyme cause misfolding, dimerization instability, and increased aggregate formation in ALS. However, there is a lack of information on the dimerization of SOD1 monomers and the mechanistic underpinnings on how the pathogenic mutations disrupt the dimerization mechanism. Here, we presented microsecond-scale molecular dynamics (MD) simulations to unravel how interface-based mutations compromise SOD1 dimerization and provide mechanistic understanding into the corresponding process using WT and three interface-based mutant systems (A4V, T54R, and I113T). Structural stability analysis showed that the mutant systems displayed disparate variations in the catalytic sites which may directly alter the stability and activity of the SOD1 enzyme. Based on the dynamic network analysis and principal component analysis, it has been identified that the mutations weakened the correlated motions along the dimer interface and altered the protein conformational behavior, thus weakening the stability of dimer formation. Moreover, the simulation results identified crucial residues such as G51, D52, G114, I151, and Q153 in establishing the dimerization interaction network, which were weakened or absent in the presence of interfacial mutants. Surface potential analysis on mutant systems also displayed changes in the dimerization potential, thus showing the unfavorable dimer formation. Furthermore, network analysis identified the hotspot residues necessary for SOD1 signal transduction which were surprisingly found in the catalytic sites rather than the anticipated dimerization interface.
Collapse