Conservation analysis and decomposition of residue correlation networks in the phospholipase A2 superfamily (PLA2s): Insights into the structure-function relationships of snake venom toxins

Phospholipases A2 (PLA(2)s) comprise a superfamily of glycerophospholipids hydrolyzing enzymes present in many organisms in nature, whose catalytic activity was majorly unveiled by analysis of snake venoms. The latter have pharmaceutical and biotechnological interests and can be divided into different functional sub-classes. Our goal was to identify important residues and their relation to the functional and class-specific characteristics in the PLA(2)s family with special emphasis on snake venom PLA(2)s (svPLA(2)s). We identified such residues by conservation analysis and decomposition of residue coevolution networks (DRCN), annotated the results based on the available literature on PLA(2)s, structural analysis and molecular dynamics simulations, and related the results to the phylogenetic distribution of these proteins. A filtered alignment of PLA(2)s revealed 14 highly conserved positions and 3 sets of coevolved residues, which were annotated according to their structural or functional role. These residues are mostly involved in ligand binding and catalysis, calcium-binding, the formation of disulfide bridges and a hydrophobic cluster close to the binding site. An independent validation of the inference of structure-function relationships from our co-evolution analysis on the svPLA2s family was obtained by the analysis of the pattern of selection acting on the Viperidae and Elapidae lineages. Additionally, a molecular dynamics simulation on the Lys49 PLA(2) from Agkistrodon contortrix laticinctus was carried out to further investigate the correlation of the Lys49-Glu69 pair. Our results suggest this configuration can result in a novel conformation where the binding cavity collapses due to the approximation of two loops caused by a strong salt bridge between Glu69 and Arg34. Finally, phylogenetic analysis indicated a correlation between the presence of residues in the coevolved sets found in this analysis and the Glade localization. The results provide a guide for important positions in the family of PLA(2)s, and potential new objects of investigation. (C) 2018 Elsevier Ltd. All rights reserved.