Opioid receptors, a kind of G protein-coupled receptors (GPCRs), mainly mediate an analgesic response via allosterically transducing the signal of endogenous ligand binding in the extracellular domain to couple to effector proteins in the intracellular domain. The δ opioid receptor (DOP) is associated with emotional control besides pain control, which makes it an attractive therapeutic target. However, its allosteric mechanism and key residues responsible for the structural stability and signal communication are not completely clear. Here we utilize the Gaussian network model (GNM) and amino acid network (AAN) combined with perturbation methods to explore the issues. The constructed fcfGNMMD, where the force constants are optimized with the inverse covariance estimation based on the correlated fluctuations from the available DOP molecular dynamics (MD) ensemble, shows a better performance than traditional GNM in reproducing residue fluctuations and cross-correlations and in capturing functionally low-frequency modes. Additionally, fcfGNMMD can consider implicitly the environmental effects to some extent. The lowest mode can well divide DOP segments and identify the two sodium ion (important allosteric regulator) binding coordination shells, and from the fastest modes, the key residues important for structure stabilization are identified. Using fcfGNMMD combined with a dynamic perturbation-response method, we explore the key residues related to the sodium ion binding. Interestingly, we identify not only the key residues in sodium ion binding shells but also the ones far away from the perturbation sites, which are involved in binding with DOP ligands, suggesting the possible long-range allosteric modulation of sodium binding for the ligand binding to DOP. Furthermore, utilizing the weighted AAN combined with attack perturbations, we identify the key residues for allosteric communication. This work helps strengthen the understanding of the allosteric communication mechanism in δ opioid receptor and can provide valuable information for drug design.

Protein–deoxyribonucleic acid (DNA) interactions are important in a variety of biological processes. Accurately predicting protein-DNA binding affinity has been one of the most attractive and challenging issues in computational biology. However, the existing approaches still have much room for improvement. In this work, we propose an ensemble model for Protein-DNA Binding Affinity prediction (emPDBA), which combines six base models with one meta-model. The complexes are classified into four types based on the DNA structure (double-stranded or other forms) and the percentage of interface residues. For each type, emPDBA is trained with the sequence-based, structure-based and energy features from binding partners and complex structures. Through feature selection by the sequential forward selection method, it is found that there do exist considerable differences in the key factors contributing to intermolecular binding affinity. The complex classification is beneficial for the important feature extraction for binding affinity prediction. The performance comparison of our method with other peer ones on the independent testing dataset shows that emPDBA outperforms the state-of-the-art methods with the Pearson correlation coefficient of 0.53 and the mean absolute error of 1.11 kcal/mol. The comprehensive results demonstrate that our method has a good performance for protein-DNA binding affinity prediction. Availability and implementation: The source code is available at https://github.com/ChunhuaLiLab/emPDBA/.

No abstract

The YTH domain of YTHDF3 belongs to a class of protein "readers" recognizing the N6-methyladenosine (m 6 A) modification in mRNA. Although static crystal structure reveals m 6 A recognition by a conserved aromatic cage, the dynamic process in recognition and importance of aromatic cage residues are not completely clear. Here, molecular dynamics (MD) simulations are performed to explore the issues and negative selectivity of YTHDF3 toward unmethylated substrate. Our results reveal that there exist conformation selectivity and induced-fit in YTHDF3 binding with m 6 Amodified RNA, where recognition loop and loop6 play important roles in the specific recognition. m 6 A modification enhances the stability of YTHDF3 in complex with RNA. The methyl group of m 6 A, like a warhead, enters into the aromatic cage of YTHDF3, where Trp492 anchors the methyl group and constraints m 6 A, making m 6 A further stabilized by π-π stacking interactions from Trp438 and Trp497. In addition, the methylation enhances the hydrophobicity of adenosine, facilitating water molecules excluded out of the aromatic cage, which is another reason for the specific recognition and stronger intermolecular interaction. Finally, the comparative analyses of hydrogen bonds and binding free energy between the methylated and unmethylated complexes reveal the physical basis for the preferred recognition of m 6 A-modified RNA by YTHDF3. This study sheds light on the mechanism by which YTHDF3 specifically recognizes m 6 A-modified RNA and can provide important information for structure-based drug design.

No abstract