Characterizing the Hot Spots Involved in RON-MSPβ Complex Formation Using In Silico Alanine Scanning Mutagenesis and Molecular Dynamics Simulation

Identification of Novel Mutations in Arabidopsis thaliana DOF 4.2 Coding Gene

Purpose: DOF (DNA-binding with One Finger) proteins are plant-specific transcription factors which mediate numerous biological processes. The purpose of the current study is to report new naturally occurring mutations in the gene encoding for one of the members of DOF proteins named DOF 4.2.

Methods: The expression of zinc finger domain of DOF 4.2 (DOF 4.2-ZF) was investigated by first synthesis of cDNA library using different parts of Arabidopsis thaliana plant. Then the coding sequence for zinc finger domain of DOF 4.2 protein was prepared using nested PCR experiment and cloned into pGEX-6P-1 expression vector. Finally, the prepared construct was used for protein expression. Furthermore, molecular dynamics (MD) simulation was carried out to predict DNA binding affinity of DOF 4.2-ZF using AMBER package.

Results: For the first time a new variant of DOF 4.2-ZF protein with three mutations was detected. One of the mutations is silent while the other two mutations lead to amino acid replacement (S18G) as well as introduction of a stop codon ultimately resulting in a truncated protein production. In order to investigate whether the truncated form is able to recognize DNA binding motif, MD simulations were carried out and the results showed that the chance of binding of DOF 4.2-ZF protein to cognate DNA in its truncated form is very low.

Conclusion: The findings demonstrated that the observed mutations adversely affect the DNA binding ability of the truncated form of DOF4.2 if it is expressed in the mutant variant of A. thaliana used in this study.

Exploring the Binding Interaction of Raf Kinase Inhibitory Protein With the N-Terminal of C-Raf Through Molecular Docking and Molecular Dynamics Simulation

Protein-protein interactions are indispensable physiological processes regulating several biological functions. Despite the availability of structural information on protein-protein complexes, deciphering their complex topology remains an outstanding challenge. Raf kinase inhibitory protein (RKIP) has gained substantial attention as a favorable molecular target for numerous pathologies including cancer and Alzheimer’s disease. RKIP interferes with the RAF/MEK/ERK signaling cascade by endogenously binding with C-Raf (Raf-1 kinase) and preventing its activation. In the current investigation, the binding of RKIP with C-Raf was explored by knowledge-based protein-protein docking web-servers including HADDOCK and ZDOCK and a consensus binding mode of C-Raf/RKIP structural complex was obtained. Molecular dynamics (MD) simulations were further performed in an explicit solvent to sample the conformations for when RKIP binds to C-Raf. Some of the conserved interface residues were mutated to alanine, phenylalanine and leucine and the impact of mutations was estimated by additional MD simulations and MM/PBSA analysis for the wild-type (WT) and constructed mutant complexes. Substantial decrease in binding free energy was observed for the mutant complexes as compared to the binding free energy of WT C-Raf/RKIP structural complex. Furthermore, a considerable increase in average backbone root mean square deviation and fluctuation was perceived for the mutant complexes. Moreover, per-residue energy contribution analysis of the equilibrated simulation trajectory by HawkDock and ANCHOR web-servers was conducted to characterize the key residues for the complex formation. One residue each from C-Raf (Arg398) and RKIP (Lys80) were identified as the druggable “hot spots” constituting the core of the binding interface and corroborated by additional long-time scale (300 ns) MD simulation of Arg398Ala mutant complex. A notable conformational change in Arg398Ala mutant occurred near the mutation site as compared to the equilibrated C-Raf/RKIP native state conformation and an essential hydrogen bonding interaction was lost. The thirteen binding sites assimilated from the overall analysis were mapped onto the complex as surface and divided into active and allosteric binding sites, depending on their location at the interface. The acquired information on the predicted 3D structural complex and the detected sites aid as promising targets in designing novel inhibitors to block the C-Raf/RKIP interaction.

Identification of the retinoschisin-binding site on the retinal Na/K-ATPase

X-linked juvenile retinoschisis (XLRS) is a hereditary retinal dystrophy, caused by mutations in the RS1 gene which encodes the secreted protein retinoschisin. In recent years, several molecules have been proposed to interact with retinoschisin, including the retinal Na/K-ATPase, L-voltage gated Ca²⁺ channels, and specific sugars. We recently showed that the retinal Na/K-ATPase consisting of subunits ATP1A3 and ATP1B2 is essential for anchoring retinoschisin to plasma membranes and identified the glycosylated ATP1B2 subunit as the direct interaction partner for retinoschisin. We now aimed to precisely map the retinoschisin binding domain(s) in ATP1B2. In general, retinoschisin binding was not affected after selective elimination of individual glycosylation sites via site-directed mutagenesis as well as after full enzymatic deglycosylation of ATP1B2. Applying the interface prediction tool PresCont, two putative protein-protein interaction patches (“patch I” and “patch II”) consisting each of four hydrophobic amino acid stretches on the ATP1B2 surface were identified. These were consecutively altered by site-directed mutagenesis. Functional assays with the ATP1B2 patch mutants identified patch II and, specifically, the associated amino acid at position 240 (harboring a threonine in ATP1B2) as crucial for retinoschisin binding to ATP1B2. These and previous results led us to suggest an induced-fit binding mechanism for the interaction between retinoschisin and the Na/K-ATPase, which is dependent on threonine 240 in ATP1B2 allowing the accommodation of hyperflexible retinoschisin spikes by the associated protein-protein interaction patch on ATP1B2.

Insight Into the Binding Mechanism of p53/pDIQ-MDMX/MDM2 With the Interaction Entropy Method

The study of the p53-MDMX/MDM2 binding sites is a research hotspot for tumor drug design. The inhibition of p53-targeted MDMX/MDM2 has become an effective approach in anti-tumor drug development. In this paper, a theoretically rigorous and computationally accurate method, namely, the interaction entropy (IE) method, combined with the polarized protein-specific charge (PPC) force field, is used to explore the difference in the binding mechanism between p53-MDMX and p53-MDM2. The interaction of a 12mer peptide inhibitor (pDIQ), which is similar to p53 in structure, with MDMX/MDM2 is also studied. The results demonstrate that p53/pDIQ with MDM2 generates a stronger interaction than with MDMX. Compared to p53, pDIQ has larger binding free energies with MDMX and MDM2. According to the calculated binding free energies, the differences in the binding free energy among the four complexes that are obtained from the combination of PPC and IE are more consistent with the experimental values than with the results from the combination of the non-polarizable AMBER force field and IE. In addition, according to the decomposition of the binding free energy, the van der Waals (vdW) interactions are the main driving force for the binding of the four complexes. They are also the main source of the weaker binding affinity of p53/pDIQ-MDMX relative to p53/pDIQ-MDM2. Compared with p53-MDMX/MDM2, according to the analysis of the residue decomposition, the predicated total residue contributions are higher in pDIQ-MDMX/MDM2 than in p53-MDMX/MDM2, which explains why pDIQ has higher binding affinity than p53 with MDMX/MDM2. The current study provides theoretical guidance for understanding the binding mechanisms and designing a potent dual inhibitor that is targeted to MDMX/MDM2.

dbMPIKT: a database of kinetic and thermodynamic mutant protein interactions

Background

Protein-protein interactions (PPIs) play important roles in biological functions. Studies of the effects of mutants on protein interactions can provide further understanding of PPIs. Currently, many databases collect experimental mutants to assess protein interactions, but most of these databases are old and have not been updated for several years.

Results

To address this issue, we manually curated a kinetic and thermodynamic database of mutant protein interactions (dbMPIKT) that is freely accessible at our website. This database contains 5291 mutants in protein interactions collected from previous databases and the literature published within the last three years. Furthermore, some data analysis, such as mutation number, mutation type, protein pair source and network map construction, can be performed online.

Conclusion

Our work can promote the study on PPIs, and novel information can be mined from the new database. Our database is available in http://DeepLearner.ahu.edu.cn/web/dbMPIKT/ for use by all, including both academics and non-academics.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2493-7) contains supplementary material, which is available to authorized users.