Cooperative Assembly of Co-Smad4 MH1 with R-Smad1/3 MH1 on DNA: A Molecular Dynamics Simulation Study

Background: Smads, the homologs of Sma and MAD proteins, play a key role in gene expression regulation in the transforming growth factor-β (TGF-β) signaling pathway. Recent experimental studies have revealed that Smad4/R-Smad heterodimers bound on DNA are energetically more favorable than homodimeric R-Smad/R-Smad complexes bound on DNA, which indicates that Smad4 might act as binding vehicle to cooperatively assemble with activated R-Smads on DNA in the nucleus. However, the details of interaction mechanism for cooperative recruitment of Smad4 protein to R-Smad proteins on DNA, and allosteric communication between the Smad4-DNA and R-Smad-DNA interfaces via DNA mediating are not yet clear so far. Methodology: In the present work, we have constructed a series of Smadn+DNA+Smadn (n = 1, 3, 4) models and carried out molecular dynamics simulations, free energy calculations and DNA dynamics analysis for them to study the interaction properties of Smadn (n = 1, 3, 4) with DNA molecule. Results: The results revealed that the binding of Smad4 protein to DNA molecule facilitates energetically the formation of the heteromeric Smad4+DNA+Smad1/3 complex by increasing the affinity of Smad1/3 with DNA molecule. Further investigations through the residue/base motion correlation and DNA dynamics analyses predicted that the binding of Smad4 protein to DNA molecule in the heteromeric Smad4+DNA+Smad1/3 model induces an allosteric communication from the Smad4-DNA interface to Smad1/Smad3-DNA interface via DNA base-pair helical motions, surface conformation changes and new hydrogen bond formations. The present work theoretically explains the mechanism of cooperative recruitment of Smad4 protein to Smad1/3 protein via DNA-mediated indirect readout mode in the nucleus.