Molecular modeling and dynamic simulation of chicken Mx protein with the S631N polymorphism

Impact of nsSNPs in human AIM2 and IFI16 gene: a comprehensive in silico analysis

AIM2 and IFI16 are the most studied members of AIM2-like receptors (ALRs) in humans and share a common N-Terminal PYD domain and C-terminal HIN domain. The HIN domain binds to dsDNA in response to the invasion of bacterial and viral DNA, and the PYD domain directs apoptosis-associated speck-like protein via protein-protein interactions. Hence, activation of AIM2 and IFI16 is crucial for protection against pathogenic assaults, and any genetic variation in these inflammasomes can dysregulate the human immune system. In this study, different computational tools were used to identify the most deleterious and disease-causing non-synonymous single nucleotide polymorphisms (nsSNPs) in AIM2 and IFI16 proteins. Molecular dynamic simulation was performed for the top damaging nsSNPs to study single amino acid substitution-induced structural alterations in AIM2 and IFI16. The observed results suggest that the variants G13V, C304R, G266R, and G266D for AIM2, and G13E and C356F are deleterious and affect structural integrity. We hope that the suggested deleterious nsSNPs and structural dynamics of AIM2 and IFI16 variants will guide future research to better understand the function of these variants with large-scale studies and may assist in fresher therapeutics focusing on these polymorphisms.Communicated by Ramaswamy H. Sarma.

Deciphering the intrinsic dynamics of unphosphorylated IRAK4 kinase bound to type I and type II inhibitors

Interleukin-1 receptor-associated kinase 4 (IRAK4) is a vital protein involved in Toll-like and interleukin-1 receptor signal transduction. Several studies have reported regarding the crystal structure, dynamic properties, and interactions with inhibitors of the phosphorylated form of IRAK4. However, no dynamic properties of inhibitor-bound unphosphorylated IRAK4 have been previously studied. Herein, we report the intrinsic dynamics of unphosphorylated IRAK4 (uIRAK4) bound to type I and type II inhibitors. The corresponding apo and inhibitor-bound forms of uIRAK4 were subjected to three independent simulations of 500 ns (total 1.5 μs) each, and their trajectories were analyzed. The results indicated that all three systems were relatively stable, except for the type II inhibitor-bound form of uIRAK4, which exhibited less compact folding and higher solvent surface area. The intra-hydrogen bonds corroborated the structural deformation of the type-II inhibitor-bound complex, which could be attributed to the long molecular structure of the type-II inhibitor. Moreover, the type II inhibitor bound to uIRAK4 showed higher binding free energy with uIRAK4 than the type I inhibitor. The free energy landscape analysis showed a reorientation of Phe330 side chain from the DFG motif at different metastable states for all the systems. The intra-residual distance between residues Lys213, Glu233, Tyr262, and Phe330 suggests a functional interplay when the inhibitors are bound to uIRAK4, thereby hinting at their crucial role in the inhibition mechanism. Ultimately, the intrinsic dynamics study observed between type I/II inhibitor-bound forms of uIRAK4 may assist in better understanding the enzyme and designing therapeutic compounds.

Computational Insights into the Structural Dynamics of MDA5 Variants Associated with Aicardi-Goutières Syndrome and Singleton-Merten Syndrome

Melanoma differentiation-associated protein 5 (MDA5) is a crucial RIG-I-like receptor RNA helicase enzyme encoded by IFIH1 in humans. Single nucleotide polymorphisms in the IFIH1 results in fatal genetic disorders such as Aicardi–Goutières syndrome and Singleton–Merten syndrome, and in increased risk of type I diabetes in humans. In this study, we chose four different amino acid substitutions of the MDA5 protein responsible for genetic disorders: MDA5^L372F, MDA5^A452T, MDA5^R779H, and MDA5^R822Q and analyzed their structural and functional relationships using molecular dynamic simulations. Our results suggest that the mutated complexes are relatively more stable than the wild-type MDA5. The radius of gyration, interaction energies, and intra-hydrogen bond analysis indicated the stability of mutated complexes over the wild type, especially MDA5^L372F and MDA5^R822Q. The dominant motions exhibited by the wild-type and mutant complexes varied significantly. Moreover, the betweenness centrality of the wild-type and mutant complexes showed shared residues for intra-signal propagation. The observed results indicate that the mutations lead to a gain of function, as reported in previous studies, due to increased interaction energies and stability between RNA and MDA5 in mutated complexes. These findings are expected to deepen our understanding of MDA5 variants and may assist in the development of relevant therapeutics against the disorders.