Telomestatin, a natural product isolated from Streptomyces anulatus, stabilizes telomeric DNA G-quadruplexes. Treatment with this ligand induces apoptosis of various cancer cells with a relatively low effect on somatic cells because of its high selectivity toward G-quadruplex over duplex DNA. A high-resolution structure of a G-quadruplex in complex with telomestatin does not yet exist because of its low solubility, and the binding nature of this ligand remains elusive. In this study, we utilized molecular binding simulations and MMGBSA binding energy analysis to decipher the nature of the binding of telomestatin to a telomeric G-quadruplex. We identified three major binding poses: bottom intercalation, top stacking, and groove binding. The top stacking mode resembles the pose observed in an NMR complex of the same G-quadruplex with the telomestatin analogue L2H. The bottom intercalation and groove binding poses were not observed in the previous studies of L2H. The bottom intercalation mode exhibited the most favorable binding energy among the three modes, while also partially intercalating into the telomeric quadruplex. The dynamic and energetic properties of these three binding modes are thoroughly examined. "Flip insertion" and "slide insertion" were observed in the bottom intercalation mode. Our findings also provide insight into the design of more selective DNA quadruplex ligands as anticancer agents in the future.

RHPS4, a potent binder to human telomeric DNA G-quadruplex, shows high efficacy in tumor cell growth inhibition. However, it's preferential binding to DNA G-quadruplex over DNA duplex (about 10 fold) remains to be improved toward its clinical application. A high resolution structure of the single-stranded telomeric DNA G-quadruplexes, or B-DNA duplex, in complex with RHPS4 is not available yet, and the binding nature of this ligand to these DNA forms remains to be elusive. In this study, we carried out 40 μs molecular dynamics binding simulations with a free ligand to decipher the binding pathway of RHPS4 to a DNA duplex and three G-quadruplex folders (parallel, antiparallel and hybrid) of the human telomeric DNA sequence. The most stable binding mode identified for the duplex, parallel, antiparallel and hybrid G-quadruplexes is an intercalation, bottom stacking, top intercalation and bottom intercalation mode, respectively. The intercalation mode with similar binding strength to both the duplex and the G-quadruplexes, explains the lack of binding selectivity of RHPS4 to the G-quadruplex form. Therefore, a ligand modification that destabilizes the duplex intercalation mode but stabilizes the G-quadruplex intercalation mode will improve the binding selectivity toward G-quadruplex. The intercalation mode of RHPS4 to both the duplex and the antiparallel and the hybrid G-quadruplex follows a base flipping-insertion mechanism rather than an open-insertion mechanism. The groove binding, the side binding and the intercalation with flipping out of base were observed to be intermediate states before the full intercalation state with paired bases.

DNA G-quadruplex structures are emerging cancer-specific targets for chemotherapeutics. Ligands that bind to and stabilize DNA G-quadruplexes have the potential to be anti-cancer drugs. Lack of binding selectivity to DNA G-quadruplex over DNA duplex remains a major challenge when attempting to develop G-quadruplex ligands into successful anti-cancer drugs. Thorough understanding of the binding nature of existing non-selective ligands that bind to both DNA quadruplex and DNA duplex will help to address this challenge. Daunomycin and doxorubicin, two commonly used anticancer drugs, are examples of non-selective DNA ligands. In this study, we extended our early all-atom binding simulation studies between doxorubicin and a DNA duplex (d(CGATCG)) to probe the binding between daunomycin and a parallel DNA quadruplex (d(TGGGGT)) and DNA duplex. In addition to the end stacking mode, which mimics the mode in the crystal structure, a pure groove binding mode was observed in our free binding simulations. The dynamic and energetic properties of these two binding modes are thoroughly examined, and a detailed comparison is made between DNA quadruplex binding modes and DNA duplex binding modes. Implications on the design of more selective DNA quadruplex ligands are also discussed. Graphical abstract Top stacking and groov binding modes from the MD simulations.

Elucidating the complex microbial interactions in biological environments requires the identification and characterization of not only the bacterial component but also the eukaryotic viruses, bacteriophage, and fungi. In a proof of concept experiment, next generation sequencing approaches, accompanied by the development of novel computational and bioinformatics tools, were utilized to examine the evolution of the microbial ecology of the avian trachea during the growth of a healthy commercial broiler flock. The flock was sampled weekly, beginning at placement and concluding at 49 days, the day before processing. Metagenomic sequencing of DNA and RNA was utilized to examine the bacteria, virus, bacteriophage, and fungal components during flock growth. The utility of using a metagenomic approach to study the avian respiratory virome was confirmed by detecting the dysbiosis in the avian respiratory virome of broiler chickens diagnosed with infection with infectious laryngotracheitis virus. This study provides the first comprehensive analysis of the ecology of the avian respiratory microbiome and demonstrates the feasibility for the use of this approach in future investigations of avian respiratory diseases.