Influence of DNA sequences on thermodynamic and structural stability of ZTA transcription factor - DNA complex: An all-atom molecular dynamics study

Abstract

The Epstein-Barr virus (EBV) is one of the cancer-causing gamma type viruses. Although more than 99% people are infected by this virus at some point, it remains in the body in a latent state, typically causing only minor symptoms. Our current understanding is that a known transcription factor (TF), the ZTA protein, binds with dsDNA (double stranded deoxyribonucleic acid) and plays crucial role in mediating the viral latent-to-lytic cycle through binding of specific ZTA responsive elements (ZREs). However, there is no clear understanding of the effect of DNA sequences on the structural stability and quantitative estimation of the binding affinity between the ZTA TF and DNA, along with their mechanistic details. In this study, we employ integrated classical all-atom molecular dynamics (MD) and enhanced sampling simulations to study the ZTA-dsDNA structural properties, thermodynamics, and mechanistic details for the ZTA protein and for two different dsDNA systems: core motif and core motif with flanking end sequences. For each system, we studied three different ZTA responsive elements (ZREs) sequences: ZRE 1, ZRE 2 and ZRE 3. We performed structural analyses, including RMSD and RMSF calculations, to assess conformational stability, along with detailed interaction profiles and hydrogen bond analysis. We conducted residue-level and nucleic acid-level analyses to assess the important protein residues and DNA bases forming interactions between the ZTA and dsDNA systems. We also explored the effect of adding flanking end sequences to the core motif on DNA groove lengths and interstrand hydrogen bonds. Our results indicate that the flanking sequences surrounding the core motif significantly influence the structural stability and binding affinity of the ZTA-dsDNA complex. Among ZRE 1, ZRE 2, and ZRE 3, particularly when paired with their naturally occurring flanking ends, ZRE 3 exhibits higher stability and binding affinity. These findings provide insights into the molecular mechanisms underlying EBV pathogenesis and may indicate potential targets for therapeutic intervention. A detailed of the binding mechanisms will allow for the design of better-targeted therapies against EBV-associated cancers. This study will serve as a holistic benchmark for future studies on these viral protein interactions.