Microsecond MD Simulations to Explore the Structural and Energetic Differences between the Human RXRα-PPARγ vs. RXRα-PPARγ-DNA

The heterodimeric complex between retinoic X receptor alpha (RXRα) and peroxisome proliferator-activated receptor gamma (PPARγ) is one of the most important and predominant regulatory systems, controlling lipid metabolism by binding to specific DNA promoter regions. X-ray and molecular dynamics (MD) simulations have revealed the average conformation adopted by the RXRα-PPARγ heterodimer bound to DNA, providing information about how multiple domains communicate to regulate receptor properties. However, knowledge of the energetic basis of the protein-ligand and protein-protein interactions is still lacking. Here we explore the structural and energetic mechanism of RXRα-PPARγ heterodimer bound or unbound to DNA and forming complex with co-crystallized ligands (rosiglitazone and 9-cis-retinoic acid) through microsecond MD simulations, molecular mechanics generalized Born surface area binding free energy calculations, principal component analysis, the free energy landscape, and correlated motion analysis. Our results suggest that DNA binding alters correlated motions and conformational mobility within RXRα–PPARγ system that impact the dimerization and the binding affinity on both receptors. Intradomain correlated motions denotes a stronger correlation map for RXRα-PPARγ-DNA than RXRα-PPARγ, involving residues at the ligand binding site. In addition, our results also corroborated the greater role of PPARγ in regulation of the free and bound DNA state.