Background: Peptidoglycan (PGN) recognition proteins (PGRPs) are important pattern recognition receptors of the host innate immune system that are involved in the immune defense against bacterial pathogens. PGRPs have been characterized in several fish species. The PGN-binding ability is important for the function of PGRPs. However, the PGRP-PGN interaction mechanism in fish remains unclear. In the present study, the 3-D model of a long PGRP of half-smooth tongue sole (Cynoglossus semilaevis) (csPGRP-L), a marine teleost with great economic value, was constructed through the comparative modeling method, and the key amino acids involved in the interaction with Lys-type PGNs and Dap-type PGNs were analyzed by molecular dynamics and molecular docking methods.

Results: csPGRP-L possessed a typical PGRP structure, consisting of five β-sheets and four α-helices. Molecular docking showed that the van der Waals forces had a slightly larger contribution than Coulombic interaction in the csPGRP-L-PGN complex. Moreover, the binding energies of csPGRP-L-PGNs computed by MM-PBSA method revealed that csPGRP-L might selectively bind both types of MTP-PGNs and MPP-PGNs. In addition, the binding energy of each residue of csPGRP-L was also calculated, revealing that the residues involved in the interaction with Lys-type PGNs were different from that with Dap-type PGNs.

Conclusions: The 3-D structure of csPGRP-L possessed typical PGRP structure and might selectively bind both types of MTP- and MPP-PGNs, which provided useful insights to understanding the functions of fish PGRPs.