Structural and Dynamic Properties of the YTH Domain in Complex with N6-methyladenosine RNA Studied by Accelerated Molecular Dynamics Simulations

Mingwei Li,Guanglin Chen,Zhiyong Zhang

DOI: https://doi.org/10.15302/j-qb-022-0297

2023-01-01

Quantitative Biology

Abstract:Background: N6-methyl adenosine (m6A) modifications of mRNA and long non-coding RNA (lncRNAs) are known to play a significant role in regulation of gene expression and organismal development. Besides writer and eraser proteins of this dynamic modification, the YT521-B homology (YTH) domain can recognize the modification involved in numerous cellular processes. The function of proteins containing YTH domain and its binding mode with N6-Methyladenosine RNA has attracted considerable attention. However, the structural and dynamic characteristics of the YTH domain in complex with m6A RNA is still unknown.Method: This work presents results of accelerated molecular dynamics (aMD) simulations at the timescale of microseconds. Principal component analysis (PCA), molecular mechanics generalized Born surface area (MM/GBSA) calculations, contact analysis and contact-based principal component analysis (conPCA) provide new insights into structure and dynamics of the YTH-RNA complex.Results: The aMD simulations indicate that the recognition loop has a larger movement away from the binding pocket in the YTH-A3 RNA than that in the YTH-m6A3 RNA. In aMD trajectories of the apo YTH, there is a significant close-open transition of the recognition loop, that is to say, the apo YTH can take both the closed and open structure. We have found that the YTH domain binds more favorably to the methylated RNA than the non-methylated RNA. The per-residue free energy decomposition and conPCA suggest that hydrophobic residues including W380, L383-V385, W431-P434, M437, and M441-L442, may play important roles in favorable binding of the m6A RNA to the YTH domain, which is also supported by aMD simulations of a double mutated system (L383A/M437A).Conclusion: The results are in good agreement with higher structural stability of the YTH-m6A RNA than that of the YTH-A3 RNA. The addition of a methylation group on A3 can enhance its binding to the hydrophobic pocket in the YTH domain. Our simulations support a ‘conformational selection’ mechanism between the YTH-RNA binding. This work may aid in our understanding of the structural and dynamic characteristics of the YTH protein in complex with the methylated RNA.