Thermodynamics and structure of a salmon cold active goose-type lysozyme

Molecular characteristics of lysozyme G in Euryglossa orientalis; cDNA cloning, phylogenic analysis, physicochemical properties and tissue gene expression

Lysozymes are the key molecules of innate immune system against bacterial infections. In the present study, we identified the molecular characteristics, physicochemical properties, antibacterial activity, evolutionary relationship and tissue expression pattern of g-type lysozyme in Euryglossa orientalis (EuOr LysG). The full-length EuOr LysG cDNA is composed of 588 nucleotides and an open reading frame encoding a protein with 195 amino acids with more than 65 % identity to g-type lysozyme of Solea senegalensis (73 %) and Scophthalmus rhombus (64 %). Amino acid sequence alignment showed that EuOr LysG protein possessed a conserved catalytic motif (Glu71-Asp84-Asp101) and one predicted disulfide bond between Cys32 and Cys122. Phylogenetic analysis based on the g-type lysozyme sequences indicated that E. orientalis and other fish of Pleuronectiformes were diverged together in the evolutionary history. The K _m and V _max values of the recombinant EuOr LysG were 0.266 mg/ml of Micrococcus lysodeikticus as substrate and 667 U/mg of protein, respectively. The optimum temperature and pH of recombinant EuOr LysG were 45 and 6 °C, respectively. Real-time PCR analysis showed that EuOr LysG transcript was most abundant in head kidney and gill and less in muscle. We also showed that the EuOr LysG had potent lytic activity against major fish bacterial pathogens with the highest activity against Bacillus cereus and Aeromonas hydrophila. Bacterial challenge with Vibrio parahaemolyticus could upregulate LysG in immune-related tissues. Our results help to understand the molecular and physicochemical characteristics of g-type lysozyme in E. orientalis which might play an important role in host defense against the bacterial infection.

Inclusion of the orientational entropic effect and low-resolution experimental information for protein-protein docking in Critical Assessment of PRedicted Interactions (CAPRI)

Inclusion of entropy is important and challenging for protein-protein binding prediction. Here, we present a statistical mechanics-based approach to empirically consider the effect of orientational entropy. Specifically, we globally sample the possible binding orientations based on a simple shape-complementarity scoring function using an FFT-type docking method. Then, for each generated orientation, we calculate the probability through the partition function of the ensemble of accessible states, which are assumed to be represented by the set of nearby binding modes. For each mode, the interaction energy is calculated using our ITScorePP scoring function that was developed in our laboratory based on principles of statistical mechanics. Using the above protocol, we present the results of our participation in Rounds 22-27 of the Critical Assessment of PRedicted Interactions (CAPRI) experiment for 10 targets (T46-T58). Additional experimental information, such as low-resolution small-angle X-ray scattering data, was used when available. In the prediction (or docking) experiments of the 10 target complexes, we achieved correct binding modes for six targets: one with high accuracy (T47), two with medium accuracy (T48 and T57), and three with acceptable accuracy (T49, T50, and T58). In the scoring experiments of seven target complexes, we obtained correct binding modes for six targets: one with high accuracy (T47), two with medium accuracy (T49 and T50), and three with acceptable accuracy (T46, T51, and T53).

Extending RosettaDock with water, sugar, and pH for prediction of complex structures and affinities for CAPRI rounds 20-27

Rounds 20-27 of the Critical Assessment of PRotein Interactions (CAPRI) provided a testing platform for computational methods designed to address a wide range of challenges. The diverse targets drove the creation of and new combinations of computational tools. In this study, RosettaDock and other novel Rosetta protocols were used to successfully predict four of the 10 blind targets. For example, for DNase domain of Colicin E2-Im2 immunity protein, RosettaDock and RosettaLigand were used to predict the positions of water molecules at the interface, recovering 46% of the native water-mediated contacts. For α-repeat Rep4-Rep2 and g-type lysozyme-PliG inhibitor complexes, homology models were built and standard and pH-sensitive docking algorithms were used to generate structures with interface RMSD values of 3.3 Å and 2.0 Å, respectively. A novel flexible sugar-protein docking protocol was also developed and used for structure prediction of the BT4661-heparin-like saccharide complex, recovering 71% of the native contacts. Challenges remain in the generation of accurate homology models for protein mutants and sampling during global docking. On proteins designed to bind influenza hemagglutinin, only about half of the mutations were identified that affect binding (T55: 54%; T56: 48%). The prediction of the structure of the xylanase complex involving homology modeling and multidomain docking pushed the limits of global conformational sampling and did not result in any successful prediction. The diversity of problems at hand requires computational algorithms to be versatile; the recent additions to the Rosetta suite expand the capabilities to encompass more biologically realistic docking problems.

Expanding the frontiers of protein-protein modeling: from docking and scoring to binding affinity predictions and other challenges

In addition to protein-protein docking, this CAPRI edition included new challenges, like protein-water and protein-sugar interactions, or the prediction of binding affinities and ΔΔG changes upon mutation. Regarding the standard protein-protein docking cases, our approach, mostly based on the pyDock scheme, submitted correct models as predictors and as scorers for 67% and 57% of the evaluated targets, respectively. In this edition, available information on known interface residues hardly made any difference for our predictions. In one of the targets, the inclusion of available experimental small-angle X-ray scattering (SAXS) data using our pyDockSAXS approach slightly improved the predictions. In addition to the standard protein-protein docking assessment, new challenges were proposed. One of the new problems was predicting the position of the interface water molecules, for which we submitted models with 20% and 43% of the water-mediated native contacts predicted as predictors and scorers, respectively. Another new problem was the prediction of protein-carbohydrate binding, where our submitted model was very close to being acceptable. A set of targets were related to the prediction of binding affinities, in which our pyDock scheme was able to discriminate between natural and designed complexes with area under the curve = 83%. It was also proposed to estimate the effect of point mutations on binding affinity. Our approach, based on machine learning methods, showed high rates of correctly classified mutations for all cases. The overall results were highly rewarding, and show that the field is ready to move forward and face new interesting challenges in interactomics.

Enzyme characterisation and gene expression profiling of Atlantic salmon chicken- and goose-type lysozymes

Lysozymes represent important innate immune components against bacteria. In this study, Atlantic salmon (Salmo salar) goose (g-) and chicken (c-) types of lysozyme were subjected to protein characterisations and tissue expression analyses. Specific bacterial protein inhibitors of g- and c-type lysozymes were employed to discriminate between respective enzyme activities. Blood, gills and liver contained activities exclusive for the g-type lysozyme. Only haematopoietic organs (head kidney and spleen) contained enzyme activities of both g- and c-lysozyme enzymes and c-type activity was not found outside these organs. Gene transcript levels proportional to enzyme activity levels were detected for the g-type lysozyme but not for the c-type. In vitro studies revealed significant induction of c-type gene expression and enzyme activity in macrophages after incubation with lipopolysaccharide (LPS) while expression of the g-type lysozyme gene was unaffected. The activity of purified native c-type enzyme was profoundly reduced by divalent cations and displayed low tolerance to monovalent cations, while the native g-type lysozyme was stimulated by monovalent cations and tolerated low concentrations of divalent cations. Activities of both enzymes increased with temperature elevations up to 60°C. The native g-type lysozyme responses to temperature in particular are in apparent conflict to the ones for the recombinant salmon g-lysozyme. Our results imply separate expression regulations and different functions of c- and g-type lysozymes in salmon. LPS-induced expression of c-type lysozyme and broad constitutive tissue distribution of g-type lysozyme in salmon is different from findings in other studied fish species.

Structural basis of bacterial defense against g-type lysozyme-based innate immunity

Gram-negative bacteria can produce specific proteinaceous inhibitors to defend themselves against the lytic action of host lysozymes. So far, four different lysozyme inhibitor families have been identified. Here, we report the crystal structure of the Escherichia coli periplasmic lysozyme inhibitor of g-type lysozyme (PliG-Ec) in complex with Atlantic salmon g-type lysozyme (SalG) at a resolution of 0.95 Å, which is exceptionally high for a complex of two proteins. The structure reveals for the first time the mechanism of g-type lysozyme inhibition by the PliG family. The latter contains two specific conserved regions that are essential for its inhibitory activity. The inhibitory complex formation is based on a double 'key-lock' mechanism. The first key-lock element is formed by the insertion of two conserved PliG regions into the active site of the lysozyme. The second element is defined by a distinct pocket of PliG accommodating a lysozyme loop. Computational analysis indicates that this pocket represents a suitable site for small molecule binding, which opens an avenue for the development of novel antibacterial agents that suppress the inhibitory activity of PliG.

Genomic characterization, phylogeny and gene regulation of g-type lysozyme in sole (Solea senegalensis)

The g-type lysozyme is a key protein of the innate immune system to fight bacterial infections. In this study we cloned and characterized the gene encoding for g-type lysozyme in Senegalese sole (Solea senegalensis). The deduced amino acid sequence comprised 195 residues containing the three conserved catalytic residues and two cysteines. A BAC analysis revealed that the gene is structured in 5 exons and 4 introns. Also, two polyadenylation signals that generate two cDNAs differing in 3'-UTR length were detected. Promoter analysis showed the presence of the main cis-acting elements involved in the transcriptional regulation of the gene. At genomic level, the g-type lysozyme was associated with mucolipin 1 and the peptidoglycan recognition protein 2 conforming a cluster of antidefensive genes with a well-conserved synteny across Percomorpha. FISH analysis using the BAC clone revealed a single hybridization signal located in an acrocentric chromosome pair. The phylogenetic analysis confirmed that the g-type lysozyme represents a complex group in fish that has been shaped by gene duplications and diversification with several positions under Darwinian selection. Expression analysis in juvenile tissues indicated that transcript levels were higher in gills, spleen and heart. During development, gene expression activated just at the beginning of metamorphosis, increasing progressively until climax. Hormonal treatments demonstrated that this gene was regulated positively by thyroid hormones during development and negatively by dexamethasone. In contrast, no response was observed after all-trans retinoic acid or 4-diethylaminobenzaldehyde treatments. Finally, treatments using lipopolysaccharide, lipoteichoic acid, peptidoglycan, zymosan and poly(I:C) activated gene expression in a time- and tissue-specific manner. Taken together, data indicate that g-type lysozyme is a high evolutionary conserved gene that diversified to adapt to changing environment and pathogen conditions. Gene expression can be activated by diverse pathogen stimuli and modulated by physiological factors with important consequences for the aquaculture of this species.