Decoding the TAome and computational insights into parDE toxin-antitoxin systems in Pseudomonas aeruginosa

Toxin-antitoxin (TA) modules are widely found in the genomes of pathogenic bacteria. They regulate vital cellular functions like transcription, translation, and DNA replication, and are therefore essential to the survival of bacteria under stress. With a focus on the type II parDE modules, this study thoroughly examines TAome in Pseudomonas aeruginosa, a bacterium well-known for its adaptability and antibiotic resistance. We explored the TAome in three P. aeruginosa strains: ATCC 27,853, PAO1, and PA14, and found 15 type II TAs in ATCC 27,853, 12 in PAO1, and 13 in PA14, with significant variation in the associated mobile genetic elements. Five different parDE homologs were found by further TAome analysis in ATCC 27,853, and their relationships were confirmed by sequence alignments and precise genomic positions. After comparing these ParDE modules' sequences to those of other pathogenic bacteria, it was discovered that they were conserved throughout many taxa, especially Proteobacteria. Nucleic acids were predicted as potential ligands for ParD antitoxins, whereas ParE toxins interacted with a wide range of small molecules, indicating a diverse functional repertoire. The interaction interfaces between ParDE TAs were clarified by protein-protein interaction networks and docking studies, which also highlighted important residues involved in binding. This thorough examination improves our understanding of the diversity, evolutionary dynamics, and functional significance of TA systems in P. aeruginosa, providing insights into their roles in bacterial physiology and pathogenicity.