Molecular dynamics simulation and steered molecular dynamics simulation on irisin dimers

Irisin is found closely associated with promoting the browning of beige fat cells in white adipose tissue. The crystal structure reveals that irisin forms a continuous inter-subunit beta-sheet dimer. Here, molecular dynamics (MD) simulation and steered molecular dynamics (SMD) simulation were performed to investigate the dissociation process and the intricate interactions between the two irisin monomers. In the process of MD, the interactions between the monomers were roughly analyzed through the average numbers of both hydrophobic contacts and H-bonds. Then, SMD was performed to investigate the accurate interaction energy between the monomers. By the analysis of dissociation energy, the van derWaals (vdW) force was identified as the major energy to maintain the dimer structure, which also verified the results of MD simulation. Meanwhile, 11 essential residues were discovered by the magnitude of rupture force during dissociation. Among them, residues Arg75, Glu79, Ile77, Ala88, and Trp90 were reported in a previous study using the method of mutagenesis and size exclusion chromatography, and several new important residues (Arg72, Leu74, Phe76, Gln78, Val80, and Asp91) were also identified. Interestingly, the new important residues that we discovered and the important residues that were reported are located in the opposite side of the beta-sheet of the dimer.