Active site detection by spatial conformity and electrostatic analysis--unravelling a proteolytic function in shrimp alkaline phosphatase

Transcriptomic analysis of the antioxidant responses and immunomodulatory effects of dietary melatonin in red swamp crayfish (Procambarus clarkii)

This study aimed to investigate the effects of dietary melatonin (MT) levels on the antioxidant capacity, immunomodulatory, and transcriptional regulation of red swamp crayfish. Six experimental diets with different levels of MT (0, 22.5, 41.2, 82.7, 165.1, and 329.2 mg/kg diet) were fed to juvenile crayfish for 60 d. The transcriptome data of the control group and the group supplemented with dietary MT at 165.1 mg/kg were obtained using RNA-seq. In total, 3653 differentially expressed genes (2082 up-regulated and 1571 down-regulated) were identified. Pathways and genes related to antioxidant immune and growth performance were verified by qRT-PCR. The total hemocyte count, phagocytosis rate, and respiratory burst were significantly increased in the MT (165.1 mg/kg) group compared to the control group. Analysis of antioxidant immune-related enzymes in the hepatopancreas demonstrated that dietary MT (165.1 mg/kg) significantly increased activities of catalase, superoxide dismutase, glutathione reductase, and glutathione peroxidase and significantly decreased aspartate aminotransferase and alanine aminotransferase activity. At the transcriptional level, dietary MT up-regulated expression levels of genes associated with antioxidant immune and development, which included toll-like receptors, Crustin, C-type lectin, and so on. To conclude, MT could be used as a supplement in crayfish feed to increase immunity and antioxidant capacity and according to the broken line regression, the ideal MT concentration was the 159.02 mg/kg. Overall, this study demonstrates the role of melatonin in the antioxidant responses and immunomodulatory of Procambarus clarkii, laying the foundation for the development of melatonin as a feed additive in the aquaculture of this species.

De Novo Arginine Synthesis Is Required for Full Virulence of Xanthomonas arboricola pv. juglandis During Walnut Bacterial Blight Disease

Walnut blight (WB) disease caused by Xanthomonas arboricola pv. juglandis (Xaj) threatens orchards worldwide. Nitrogen metabolism in this bacterial pathogen is dependent on arginine, a nitrogen-enriched amino acid that can either be synthesized or provided by the plant host. The arginine biosynthetic pathway uses argininosuccinate synthase (argG), associated with increased bacterial virulence. We examined the effects of bacterial arginine and nitrogen metabolism on the plant response during WB by proteomic analysis of the mutant strain Xaj argG^-. Phenotypically, the mutant strain produced 42% fewer symptoms and survived in the plant tissue with 2.5-fold reduced growth compared with wild type, while showing itself to be auxotrophic for arginine in vitro. Proteomic analysis of infected tissue enabled the profiling of 676 Xaj proteins and 3,296 walnut proteins using isobaric labeling in a data-dependent acquisition approach. Comparative analysis of differentially expressed proteins revealed distinct plant responses. Xaj wild type (WT) triggered processes of catabolism and oxidative stress in the host under observed disease symptoms, while most of the host biosynthetic processes triggered by Xaj WT were inhibited during Xaj argG^- infection. Overall, the Xaj proteins revealed a drastic shift in carbon and energy management induced by disruption of nitrogen metabolism while the top differentially expressed proteins included a Fis transcriptional regulator and a peptidyl-prolyl isomerase. Our results show the critical role of de novo arginine biosynthesis to sustain virulence and minimal growth during WB. This study is timely and critical as copper-based control methods are losing their effectiveness, and new sustainable methods are urgently needed in orchard environments.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

The plant-based chimeric antimicrobial protein SlP14a-PPC20 protects tomato against bacterial wilt disease caused by Ralstonia solanacearum

Cecropin-B (CecB) is a peptide with well-established antimicrobial properties against different phytopathogenic bacteria. Despite modest action against Ralstonia solanacearum, its animal source limits the acceptance in transgenic applications. To overcome this, we selected eight alpha-helical (AH) cationic peptides derived from plant protein sequences and investigated their antimicrobial properties against R. solanacearum. Remarkably, PPC20 (a linear AH-peptide present in phosphoenolpyruvate carboxylase) has a three-fold lower lethal dose on R. solanacearum than CecB and lower toxicity to human intestinal epithelial cells. Linking PPC20 to SlP14a (part of a pathogenesis-related protein) established an apoplast-targeted protein providing a means of secreting and stabilizing the antimicrobial peptide in the plant compartment colonized by the pathogen. SlP14a is also a potential antimicrobial, homologous to a human elastase which likely targets outer membrane proteins in Gram-negative bacteria. Recombinant SlP14a-PPC20 showed antibacterial activity against R. solanacearum in vitro, making it a promising candidate for plant protection. This was confirmed with genetically-modified tomato plants engineered to express SlP14a-PPC20, in which bacterial populations in stems were reduced compared to inoculated wild-type control plants. Disease symptoms were also markedly less severe in SlP14a-PPC20-expressing plants, demonstrating a viable strategy to improve resistance against bacterial wilt in tomato.

Computational insights of K1444N substitution in GAP-related domain of NF1 gene associated with neurofibromatosis type 1 disease: a molecular modeling and dynamics approach

The NF1 gene encodes for neurofibromin protein, which is ubiquitously expressed, but most highly in the central nervous system. Non-synonymous SNPs (nsSNPs) in the NF1 gene were found to be associated with Neurofibromatosis Type 1 disease, which is characterized by the growth of tumors along nerves in the skin, brain, and other parts of the body. In this study, we used several in silico predictions tools to analyze 16 nsSNPs in the RAS-GAP domain of neurofibromin, the K1444N (K1423N) mutation was predicted as the most pathogenic. The comparative molecular dynamic simulation (MDS; 50 ns) between the wild type and the K1444N (K1423N) mutant suggested a significant change in the electrostatic potential. In addition, the RMSD, RMSF, Rg, hydrogen bonds, and PCA analysis confirmed the loss of flexibility and increase in compactness of the mutant protein. Further, SASA analysis revealed exchange between hydrophobic and hydrophilic residues from the core of the RAS-GAP domain to the surface of the mutant domain, consistent with the secondary structure analysis that showed significant alteration in the mutant protein conformation. Our data concludes that the K1444N (K1423N) mutant lead to increasing the rigidity and compactness of the protein. This study provides evidence of the benefits of the computational tools in predicting the pathogenicity of genetic mutations and suggests the application of MDS and different in silico prediction tools for variant assessment and classification in genetic clinics.

Structural characterization of a novel KH-domain containing plant chloroplast endonuclease

Chlamydomonas reinhardtii is a single celled alga that undergoes apoptosis in response to UV-C irradiation. UVI31+, a novel UV-inducible DNA endonuclease in C. reinhardtii, which normally localizes near cell wall and pyrenoid regions, gets redistributed into punctate foci within the whole chloroplast, away from the pyrenoid, upon UV-stress. Solution NMR structure of the first putative UV inducible endonuclease UVI31+ revealed an α₁–β₁–β₂–α₂–α₃–β₃ fold similar to BolA and type II KH-domain ubiquitous protein families. Three α−helices of UVI31+ constitute one side of the protein surface, which are packed to the other side, made of three-stranded β–sheet, with intervening hydrophobic residues. A twenty-three residues long polypeptide stretch (D54-H76) connecting β₁ and β₂ strands is found to be highly flexible. Interestingly, UVI31+ recognizes the DNA primarily through its β–sheet. We propose that the catalytic triad residues involving Ser114, His95 and Thr116 facilitate DNA endonuclease activity of UVI31+. Further, decreased endonuclease activity of the S114A mutant is consistent with the direct participation of Ser114 in the catalysis. This study provides the first structural description of a plant chloroplast endonuclease that is regulated by UV-stress response.

Computational analysis of perturbations in the post-fusion Dengue virus envelope protein highlights known epitopes and conserved residues in the Zika virus

The dramatic transformation of the Zika virus (ZIKV) from a relatively unknown virus to a pathogen generating global-wide panic has exposed the dearth of detailed knowledge about this virus. Decades of research in the related Dengue virus (DENV), finally culminating in a vaccine registered for use in endemic regions (CYD-TDV) in three countries, provides key insights in developing strategies for tackling ZIKV, which has caused global panic to microcephaly and Guillain-Barre Syndrome. Dengue virus (DENV), a member of the family Flaviviridae, the causal agent of the self-limiting Dengue fever and the potentially fatal hemorrhagic fever/dengue shock syndrome, has been a scourge in tropical countries for many centuries. The recently solved structure of mature ZIKV (PDB ID:5IRE) has provided key insights into the structure of the envelope (E) and membrane (M) proteins, the primary target of neutralizing antibodies. The previously established MEPP methodology compares two conformations of the same protein and identifies residues with significant spatial and electrostatic perturbations. In the current work, MEPP analyzed the pre-and post-fusion DENV type 2 envelope (E) protein, and identified several known epitopes (His317, Tyr299, Glu26, Arg188, etc.) (MEPPitope). These residues are overwhelmingly conserved in ZIKV and all DENV serotypes, and also enumerates residue pairs that undergo significant polarity reversal. Characterization of α-helices in E-proteins show that α1 is not conserved in the sequence space of ZIKV and DENV. Furthermore, perturbation of α1 in the post-fusion DENV structure includes a known epitope Asp215, a residue absent in the pre-fusion α1. A cationic β-sheet in the GAG-binding domain that is stereochemically equivalent in ZIKV and all DENV serotypes is also highlighted due to a residue pair (Arg286-Arg288) that has a significant electrostatic polarity reversal upon fusion. Finally, two highly conserved residues (Thr32 and Thr40), with little emphasis in existing literature, are found to have significant electrostatic perturbation. Thus, a combination of different computational methods enable the rapid and rational detection of critical residues as epitopes in the search for an elusive therapy or vaccine that neutralizes multiple members of the Flaviviridae family. These secondary structures are conserved in the related Dengue virus (DENV), and possibly rationalize isolation techniques particle adsorption on magnetic beads coated with anionic polymers and anionic antiviral agents (viprolaxikine) for DENV. These amphipathic α-helices could enable design of molecules for inhibiting α-helix mediated protein-protein interactions. Finally, comparison of these secondary structures in proteins from related families illuminate subtle changes in the proteins that might render them ineffective to previously successful drugs and vaccines, which are difficult to identify by a simple sequence or structural alignment. Finally, conflicting results about residues that are involved in neutralizing a DENV-E protein by the potent antibody 5J7 (PDB ID:3J6U) are reported.

The Secreted Protease PrtA Controls Cell Growth, Biofilm Formation and Pathogenicity in Xylella fastidiosa

Pierce's disease (PD) is a deadly disease of grapevines caused by the Gram-negative bacterium Xylella fastidiosa. Though disease symptoms were formerly attributed to bacteria blocking the plant xylem, this hypothesis is at best overly simplistic. Recently, we used a proteomic approach to characterize the secretome of X. fastidiosa, both in vitro and in planta, and identified LesA as one of the pathogenicity factors of X. fastidiosa in grapevines that leads to leaf scorching and chlorosis. Herein, we characterize another such factor encoded by PD0956, designated as an antivirulence secreted protease "PrtA" that displays a central role in controlling in vitro cell proliferation, length, motility, biofilm formation, and in planta virulence. The mutant in X. fastidiosa exhibited reduced cell length, hypermotility (and subsequent lack of biofilm formation) and hypervirulence in grapevines. These findings are supported by transcriptomic and proteomic analyses with corresponding plant infection data. Of particular interest, is the hypervirulent response in grapevines observed when X. fastidiosa is disrupted for production of PrtA, and that PD-model tobacco plants transformed to express PrtA exhibited decreased symptoms after infection by X. fastidiosa.

Computational analysis of perturbations in the post-fusion Dengue virus envelope protein highlights known epitopes and conserved residues in the Zika virus

The dramatic transformation of the Zika virus (ZIKV) from a relatively unknown virus to a pathogen generating global-wide panic has exposed the dearth of detailed knowledge about this virus. Decades of research in the related Dengue virus (DENV), finally culminating in a vaccine registered for use in endemic regions (CYD-TDV) in three countries, provides key insights in developing strategies for tackling ZIKV, which has caused global panic to microcephaly and Guillain-Barre Syndrome. Dengue virus (DENV), a member of the family Flaviviridae, the causal agent of the self-limiting Dengue fever and the potentially fatal hemorrhagic fever/dengue shock syndrome, has been a scourge in tropical countries for many centuries. The recently solved structure of mature ZIKV (PDB ID:5IRE) has provided key insights into the structure of the envelope (E) and membrane (M) proteins, the primary target of neutralizing antibodies. The previously established MEPP methodology compares two conformations of the same protein and identifies residues with significant spatial and electrostatic perturbations. In the current work, MEPP analyzed the pre-and post-fusion DENV type 2 envelope (E) protein, and identified several known epitopes (His317, Tyr299, Glu26, Arg188, etc.) (MEPPitope). These residues are overwhelmingly conserved in ZIKV and all DENV serotypes, and also enumerates residue pairs that undergo significant polarity reversal. Characterization of α-helices in E-proteins show that α1 is not conserved in the sequence space of ZIKV and DENV. Furthermore, perturbation of α1 in the post-fusion DENV structure includes a known epitope Asp215, a residue absent in the pre-fusion α1. A cationic β-sheet in the GAG-binding domain that is stereochemically equivalent in ZIKV and all DENV serotypes is also highlighted due to a residue pair (Arg286-Arg288) that has a significant electrostatic polarity reversal upon fusion. Finally, two highly conserved residues (Thr32 and Thr40), with little emphasis in existing literature, are found to have significant electrostatic perturbation. Thus, a combination of different computational methods enable the rapid and rational detection of critical residues as epitopes in the search for an elusive therapy or vaccine that neutralizes multiple members of the Flaviviridae family. These secondary structures are conserved in the related Dengue virus (DENV), and possibly rationalize isolation techniques particle adsorption on magnetic beads coated with anionic polymers and anionic antiviral agents (viprolaxikine) for DENV. These amphipathic α-helices could enable design of molecules for inhibiting α-helix mediated protein-protein interactions. Finally, comparison of these secondary structures in proteins from related families illuminate subtle changes in the proteins that might render them ineffective to previously successful drugs and vaccines, which are difficult to identify by a simple sequence or structural alignment. Finally, conflicting results about residues that are involved in neutralizing a DENV-E protein by the potent antibody 5J7 (PDB ID:3J6U) are reported.

Sequence/structural analysis of xylem proteome emphasizes pathogenesis-related proteins, chitinases and β-1, 3-glucanases as key players in grapevine defense against Xylella fastidiosa

Background. Xylella fastidiosa, the causative agent of various plant diseases including Pierce’s disease in the US, and Citrus Variegated Chlorosis in Brazil, remains a continual source of concern and economic losses, especially since almost all commercial varieties are sensitive to this Gammaproteobacteria. Differential expression of proteins in infected tissue is an established methodology to identify key elements involved in plant defense pathways.

Methods. In the current work, we developed a methodology named CHURNER that emphasizes relevant protein functions from proteomic data, based on identification of proteins with similar structures that do not necessarily have sequence homology. Such clustering emphasizes protein functions which have multiple copies that are up/down-regulated, and highlights similar proteins which are differentially regulated. As a working example we present proteomic data enumerating differentially expressed proteins in xylem sap from grapevines that were infected with X. fastidiosa.

Results. Analysis of this data by CHURNER highlighted pathogenesis related PR-1 proteins, reinforcing this as the foremost protein function in xylem sap involved in the grapevine defense response to X. fastidiosa. β-1, 3-glucanase, which has both anti-microbial and anti-fungal activities, is also up-regulated. Simultaneously, chitinases are found to be both up and down-regulated by CHURNER, and thus the net gain of this protein function loses its significance in the defense response.

Discussion. We demonstrate how structural data can be incorporated in the pipeline of proteomic data analysis prior to making inferences on the importance of individual proteins to plant defense mechanisms. We expect CHURNER to be applicable to any proteomic data set.

Sequence/structural analysis of xylem proteome emphasizes pathogenesis-related proteins, chitinases and β-1, 3-glucanases as key players in grapevine defense against Xylella fastidiosa

Background. Xylella fastidiosa, the causative agent of various plant diseases including Pierce's disease in the US, and Citrus Variegated Chlorosis in Brazil, remains a continual source of concern and economic losses, especially since almost all commercial varieties are sensitive to this Gammaproteobacteria. Differential expression of proteins in infected tissue is an established methodology to identify key elements involved in plant defense pathways. Methods. In the current work, we developed a methodology named CHURNER that emphasizes relevant protein functions from proteomic data, based on identification of proteins with similar structures that do not necessarily have sequence homology. Such clustering emphasizes protein functions which have multiple copies that are up/down-regulated, and highlights similar proteins which are differentially regulated. As a working example we present proteomic data enumerating differentially expressed proteins in xylem sap from grapevines that were infected with X. fastidiosa. Results. Analysis of this data by CHURNER highlighted pathogenesis related PR-1 proteins, reinforcing this as the foremost protein function in xylem sap involved in the grapevine defense response to X. fastidiosa. β-1, 3-glucanase, which has both anti-microbial and anti-fungal activities, is also up-regulated. Simultaneously, chitinases are found to be both up and down-regulated by CHURNER, and thus the net gain of this protein function loses its significance in the defense response. Discussion. We demonstrate how structural data can be incorporated in the pipeline of proteomic data analysis prior to making inferences on the importance of individual proteins to plant defense mechanisms. We expect CHURNER to be applicable to any proteomic data set.

DOCLASP - Docking ligands to target proteins using spatial and electrostatic congruence extracted from a known holoenzyme and applying simple geometrical transformations

The ability to accurately and effectively predict the interaction between proteins and small drug-like compounds has long intrigued researchers for pedagogic, humanitarian and economic reasons. Protein docking methods (AutoDock, GOLD, DOCK, FlexX and Glide to name a few) rank a large number of possible conformations of protein-ligand complexes using fast algorithms. Previously, it has been shown that structural congruence leading to the same enzymatic function necessitates the congruence of electrostatic properties (CLASP). The current work presents a methodology for docking a ligand into a target protein, provided that there is at least one known holoenzyme with ligand bound - DOCLASP (Docking using CLASP). The contact points of the ligand in the holoenzyme defines a motif, which is used to query the target enzyme using CLASP. If there are significant matches, the holoenzyme and the target protein are superimposed based on congruent atoms. The same linear and rotational transformations are also applied to the ligand, thus creating a unified coordinate framework having the holoenzyme, the ligand and the target enzyme. In the current work, the dipeptidyl peptidase-IV inhibitor vildagliptin was docked to the PI-PLC structure complexed with myo-inositol using DOCLASP. Also, corroboration of the docking of phenylthiourea to the modelled structure of polyphenol oxidase (JrPPO1) from walnut is provided based on the subsequently solved structure of JrPPO1 (PDBid:5CE9). Analysis of the binding of the antitrypanosomial drug suramin to nine non-homologous proteins in the PDB database shows a diverse set of binding motifs, and multiple binding sites in the phospholipase A2-likeproteins from the Bothrops genus of pitvipers. The conformational changes in the suramin molecule on binding highlights the challenges in docking flexible ligands into an already 'plastic' binding site. Thus, DOCLASP presents a method for 'soft docking' ligands to proteins with low computational requirements.

DOCLASP - Docking ligands to target proteins using spatial and electrostatic congruence extracted from a known holoenzyme and applying simple geometrical transformations

The ability to accurately and effectively predict the interaction between proteins and small drug-like compounds has long intrigued researchers for pedagogic, humanitarian and economic reasons. Protein docking methods (AutoDock, GOLD, DOCK, FlexX and Glide to name a few) rank a large number of possible conformations of protein-ligand complexes using fast algorithms. Previously, it has been shown that structural congruence leading to the same enzymatic function necessitates the congruence of electrostatic properties (CLASP). The current work presents a methodology for docking a ligand into a target protein, provided that there is at least one known holoenzyme with ligand bound - DOCLASP (Docking using CLASP). The contact points of the ligand in the holoenzyme defines a motif, which is used to query the target enzyme using CLASP. If there are significant matches, the holoenzyme and the target protein are superimposed based on congruent atoms. The same linear and rotational transformations are also applied to the ligand, thus creating a unified coordinate framework having the holoenzyme, the ligand and the target enzyme. In the current work, the dipeptidyl peptidase-IV inhibitor vildagliptin was docked to the PI-PLC structure complexed with myo-inositol using DOCLASP. Also, corroboration of the docking of phenylthiourea to the modelled structure of polyphenol oxidase (JrPPO1) from walnut is provided based on the subsequently solved structure of JrPPO1 (PDBid:5CE9). Analysis of the binding of the antitrypanosomial drug suramin to nine non-homologous proteins in the PDB database shows a diverse set of binding motifs, and multiple binding sites in the phospholipase A2-likeproteins from the Bothrops genus of pitvipers. The conformational changes in the suramin molecule on binding highlights the challenges in docking flexible ligands into an already ’plastic’ binding site. Thus, DOCLASP presents a method for ’soft docking’ ligands to proteins with low computational requirements.

Directed evolution induces tributyrin hydrolysis in a virulence factor of Xylella fastidiosa using a duplicated gene as a template

Duplication of genes is one of the preferred ways for natural selection to add advantageous functionality to the genome without having to reinvent the wheel with respect to catalytic efficiency and protein stability. The duplicated secretory virulence factors of Xylella fastidiosa (LesA, LesB and LesC), implicated in Pierce's disease of grape and citrus variegated chlorosis of citrus species, epitomizes the positive selection pressures exerted on advantageous genes in such pathogens. A deeper insight into the evolution of these lipases/esterases is essential to develop resistance mechanisms in transgenic plants. Directed evolution, an attempt to accelerate the evolutionary steps in the laboratory, is inherently simple when targeted for loss of function. A bigger challenge is to specify mutations that endow a new function, such as a lost functionality in a duplicated gene. Previously, we have proposed a method for enumerating candidates for mutations intended to transfer the functionality of one protein into another related protein based on the spatial and electrostatic properties of the active site residues (DECAAF). In the current work, we present in vivo validation of DECAAF by inducing tributyrin hydrolysis in LesB based on the active site similarity to LesA. The structures of these proteins have been modeled using RaptorX based on the closely related LipA protein from Xanthomonas oryzae. These mutations replicate the spatial and electrostatic conformation of LesA in the modeled structure of the mutant LesB as well, providing in silico validation before proceeding to the laborious in vivo work. Such focused mutations allows one to dissect the relevance of the duplicated genes in finer detail as compared to gene knockouts, since they do not interfere with other moonlighting functions, protein expression levels or protein-protein interaction.

DOCLASP - Docking ligands to target proteins using spatial and electrostatic congruence extracted from a known holoenzyme and applying simple geometrical transformations

The ability to accurately and effectively predict the interaction between proteins and small drug-like compounds has long intrigued researchers for pedagogic, humanitarian and economic reasons. Protein docking methods (AutoDock, GOLD, DOCK, FlexX and Glide to name a few) rank a large number of possible conformations of protein-ligand complexes using fast algorithms. Previously, it has been shown that structural congruence leading to the same enzymatic function necessitates the congruence of electrostatic properties (CLASP). The current work presents a methodology for docking a ligand into a target protein, provided that there is at least one known holoenzyme with ligand bound - DOCLASP (Docking using CLASP). The contact points of the ligand in the holoenzyme defines a motif, which is used to query the target enzyme using CLASP. If there are significant matches, the holoenzyme and the target protein are superimposed based on congruent atoms. The same linear and rotational transformations are also applied to the ligand, thus creating a unified coordinate framework having the holoenzyme, the ligand and the target enzyme. In the current work, the dipeptidyl peptidase-IV inhibitor vildagliptin was docked to the PI-PLC structure complexed with myo-inositol using DOCLASP. Also, corroboration of the docking of phenylthiourea to the modelled structure of polyphenol oxidase (JrPPO1) from walnut is provided based on the subsequently solved structure of JrPPO1 (PDBid:5CE9). Analysis of the binding of the antitrypanosomial drug suramin to nine non-homologous proteins in the PDB database shows a diverse set of binding motifs, and multiple binding sites in the phospholipase A2-likeproteins from the Bothrops genus of pitvipers. The conformational changes in the suramin molecule on binding highlights the challenges in docking flexible ligands into an already 'plastic' binding site. Thus, DOCLASP presents a method for 'soft docking' ligands to proteins with low computational requirements.

The dipeptidyl peptidase IV inhibitors vildagliptin and K-579 inhibit a phospholipase C: a case of promiscuous scaffolds in proteins

The long term side effects of any newly introduced drug is a subject of intense research, and often raging controversies. One such example is the dipeptidyl peptidase-IV (DPP4) inhibitor used for treating type 2 diabetes, which is inconclusively implicated in increased susceptibility to acute pancreatitis. Previously, based on a computational analysis of the spatial and electrostatic properties of active site residues, we have demonstrated that phosphoinositide-specific phospholipase C (PI-PLC) from Bacillus cereus is a prolyl peptidase using in vivo experiments. In the current work, we first report the inhibition of the native activity of PI-PLC by two DPP4 inhibitors - vildagliptin (LAF-237) and K-579. While vildagliptin inhibited PI-PLC at micromolar concentrations, K-579 was a potent inhibitor even at nanomolar concentrations. Subsequently, we queried a comprehensive, non-redundant set of 5000 human proteins (50% similarity cutoff) with known structures using serine protease (SPASE) motifs derived from trypsin and DPP4. A pancreatic lipase and a gastric lipase are among the proteins that are identified as proteins having promiscuous SPASE scaffolds that could interact with DPP4 inhibitors. The presence of such scaffolds in human lipases is expected since they share the same catalytic mechanism with PI-PLC. However our methodology also detects other proteins, often with a completely different enzymatic mechanism, that have significantly congruent domains with the SPASE motifs. The reported elevated levels of serum lipase, although contested, could be rationalized by inhibition of lipases reported here. In an effort to further our understanding of the spatial and electrostatic basis of DPP4 inhibitors, we have also done a comprehensive analysis of all 76 known DPP4 structures liganded to inhibitors till date. Also, the methodology presented here can be easily adopted for other drugs, and provide the first line of filtering in the identification of pathways that might be inadvertently affected due to promiscuous scaffolds in proteins.

The dipeptidyl peptidase IV inhibitors vildagliptin and K-579 inhibit a phospholipase C: a case of promiscuous scaffolds in proteins

The long term side effects of any newly introduced drug is a subject of intense research, and often raging controversies. One such example is the dipeptidyl peptidase-IV (DPP4) inhibitor used for treating type 2 diabetes, which is inconclusively implicated in increased susceptibility to acute pancreatitis. Previously, based on a computational analysis of the spatial and electrostatic properties of active site residues, we have demonstrated that phosphoinositide-specific phospholipase C (PI-PLC) from Bacillus cereus is a prolyl peptidase using in vivo experiments. In the current work, we first report the inhibition of the native activity of PI-PLC by two DPP4 inhibitors - vildagliptin (LAF-237) and K-579. While vildagliptin inhibited PI-PLC at micromolar concentrations, K-579 was a potent inhibitor even at nanomolar concentrations. Subsequently, we queried a comprehensive, non-redundant set of 5000 human proteins (50% similarity cutoff) with known structures using serine protease (SPASE) motifs derived from trypsin and DPP4. A pancreatic lipase and a gastric lipase are among the proteins that are identified as proteins having promiscuous SPASE scaffolds that could interact with DPP4 inhibitors. The presence of such scaffolds in human lipases is expected since they share the same catalytic mechanism with PI-PLC. However our methodology also detects other proteins, often with a completely different enzymatic mechanism, that have significantly congruent domains with the SPASE motifs. The reported elevated levels of serum lipase, although contested, could be rationalized by inhibition of lipases reported here. In an effort to further our understanding of the spatial and electrostatic basis of DPP4 inhibitors, we have also done a comprehensive analysis of all 76 known DPP4 structures liganded to inhibitors till date. Also, the methodology presented here can be easily adopted for other drugs, and provide the first line of filtering in the identification of pathways that might be inadvertently affected due to promiscuous scaffolds in proteins.

Binding site comparison for function prediction and pharmaceutical discovery

While structural genomics resulted in thousands of new protein crystal structures, we still do not know the functions of most of these proteins. One reason for this shortcoming is their unique sequences or folds, which leaves them assigned as proteins of 'unknown function'. Recent advances in and applications of cutting edge binding site comparison algorithms for binding site detection and function prediction have begun to shed light on this problem. Here, we review these algorithms and their use in function prediction and pharmaceutical discovery. Finding common binding sites in weakly related proteins may lead to the discovery of new protein functions and to novel ways of drug discovery.

The electrostatic profile of consecutive Cβ atoms applied to protein structure quality assessment

The structure of a protein provides insight into its physiological interactions with other components of the cellular soup. Methods that predict putative structures from sequences typically yield multiple, closely-ranked possibilities. A critical component in the process is the model quality assessing program (MQAP), which selects the best candidate from this pool of structures. Here, we present a novel MQAP based on the physical properties of sidechain atoms. We propose a method for assessing the quality of protein structures based on the electrostatic potential difference (EPD) of Cβ atoms in consecutive residues. We demonstrate that the EPDs of Cβ atoms on consecutive residues provide unique signatures of the amino acid types. The EPD of Cβ atoms are learnt from a set of 1000 non-homologous protein structures with a resolution cuto of 1.6 Å obtained from the PISCES database. Based on the Boltzmann hypothesis that lower energy conformations are proportionately sampled more, and on Annsen's thermodynamic hypothesis that the native structure of a protein is the minimum free energy state, we hypothesize that the deviation of observed EPD values from the mean values obtained in the learning phase is minimized in the native structure. We achieved an average specificity of 0.91, 0.94 and 0.93 on hg_structal, 4state_reduced and ig_structal decoy sets, respectively, taken from the Decoys `R' Us database. The source code and manual is made available at https://github.com/sanchak/mqap and permanently available on 10.5281/zenodo.7134.

Protein structure quality assessment based on the distance profiles of consecutive backbone Cα atoms

Predicting the three dimensional native state structure of a protein from its primary sequence is an unsolved grand challenge in molecular biology. Two main computational approaches have evolved to obtain the structure from the protein sequence - ab initio/de novo methods and template-based modeling - both of which typically generate multiple possible native state structures. Model quality assessment programs (MQAP) validate these predicted structures in order to identify the correct native state structure. Here, we propose a MQAP for assessing the quality of protein structures based on the distances of consecutive Cα atoms. We hypothesize that the root-mean-square deviation of the distance of consecutive Cα (RDCC) atoms from the ideal value of 3.8 Å, derived from a statistical analysis of high quality protein structures (top100H database), is minimized in native structures. Based on tests with the top100H set, we propose a RDCC cutoff value of 0.012 Å, above which a structure can be filtered out as a non-native structure. We applied the RDCC discriminator on decoy sets from the Decoys 'R' Us database to show that the native structures in all decoy sets tested have RDCC below the 0.012 Å cutoff. While most decoy sets were either indistinguishable using this discriminator or had very few violations, all the decoy structures in the fisa decoy set were discriminated by applying the RDCC criterion. This highlights the physical non-viability of the fisa decoy set, and possible issues in benchmarking other methods using this set. The source code and manual is made available at https://github.com/sanchak/mqap and permanently available on 10.5281/zenodo.7134.

Protein structure quality assessment based on the distance profiles of consecutive backbone Cα atoms

Predicting the three dimensional native state structure of a protein from its primary sequence is an unsolved grand challenge in molecular biology. Two main computational approaches have evolved to obtain the structure from the protein sequence - ab initio/de novo methods and template-based modeling - both of which typically generate multiple possible native state structures. Model quality assessment programs (MQAP) validate these predicted structures in order to identify the correct native state structure. Here, we propose a MQAP for assessing the quality of protein structures based on the distances of consecutive Cα atoms. We hypothesize that the root-mean-square deviation of the distance of consecutive Cα (RDCC) atoms from the ideal value of 3.8 Å, derived from a statistical analysis of high quality protein structures (top100H database), is minimized in native structures. Based on tests with the top100H set, we propose a RDCC cutoff value of 0.012 Å, above which a structure can be filtered out as a non-native structure. We applied the RDCC discriminator on decoy sets from the Decoys 'R' Us database to show that the native structures in all decoy sets tested have RDCC below the 0.012 Å cutoff. While most decoy sets were either indistinguishable using this discriminator or had very few violations, all the decoy structures in the fisa decoy set were discriminated by applying the RDCC criterion. This highlights the physical non-viability of the fisa decoy set, and possible issues in benchmarking other methods using this set. The source code and manual is made available at https://github.com/sanchak/mqap and permanently available on 10.5281/zenodo.7134.

Bioengineering of plant (tri)terpenoids: from metabolic engineering of plants to synthetic biology in vivo and in vitro

Terpenoids constitute a large and diverse class of natural products that serve many functions in nature. Most of the tens of thousands of the discovered terpenoids are synthesized by plants, where they function as primary metabolites involved in growth and development, or as secondary metabolites that optimize the interaction between the plant and its environment. Several plant terpenoids are economically important molecules that serve many applications as pharmaceuticals, pesticides, etc. Major challenges for the commercialization of plant-derived terpenoids include their low production levels in planta and the continuous demand of industry for novel molecules with new or superior biological activities. Here, we highlight several synthetic biology methods to enhance and diversify the production of plant terpenoids, with a foresight towards triterpenoid engineering, the least engineered class of bioactive terpenoids. Increased or cheaper production of valuable triterpenoids may be obtained by 'classic' metabolic engineering of plants or by heterologous production of the compounds in other plants or microbes. Novel triterpenoid structures can be generated through combinatorial biosynthesis or directed enzyme evolution approaches. In its ultimate form, synthetic biology may lead to the production of large amounts of plant triterpenoids in in vitro systems or custom-designed artificial biological systems.

A computational module assembled from different protease family motifs identifies PI PLC from Bacillus cereus as a putative prolyl peptidase with a serine protease scaffold

Proteolytic enzymes have evolved several mechanisms to cleave peptide bonds. These distinct types have been systematically categorized in the MEROPS database. While a BLAST search on these proteases identifies homologous proteins, sequence alignment methods often fail to identify relationships arising from convergent evolution, exon shuffling, and modular reuse of catalytic units. We have previously established a computational method to detect functions in proteins based on the spatial and electrostatic properties of the catalytic residues (CLASP). CLASP identified a promiscuous serine protease scaffold in alkaline phosphatases (AP) and a scaffold recognizing a β-lactam (imipenem) in a cold-active Vibrio AP. Subsequently, we defined a methodology to quantify promiscuous activities in a wide range of proteins. Here, we assemble a module which encapsulates the multifarious motifs used by protease families listed in the MEROPS database. Since APs and proteases are an integral component of outer membrane vesicles (OMV), we sought to query other OMV proteins, like phospholipase C (PLC), using this search module. Our analysis indicated that phosphoinositide-specific PLC from Bacillus cereus is a serine protease. This was validated by protease assays, mass spectrometry and by inhibition of the native phospholipase activity of PI-PLC by the well-known serine protease inhibitor AEBSF (IC50 = 0.018 mM). Edman degradation analysis linked the specificity of the protease activity to a proline in the amino terminal, suggesting that the PI-PLC is a prolyl peptidase. Thus, we propose a computational method of extending protein families based on the spatial and electrostatic congruence of active site residues.

A quantitative measure of electrostatic perturbation in holo and apo enzymes induced by structural changes

Biological pathways are subject to subtle manipulations that achieve a wide range of functional variation in differing physiological niches. In many instances, changes in the structure of an enzyme on ligand binding germinate electrostatic perturbations that form the basis of its changed catalytic or transcriptional efficiency. Computational methods that seek to gain insights into the electrostatic changes in enzymes require expertise to setup and computing prowess. In the current work, we present a fast, easy and reliable methodology to compute electrostatic perturbations induced by ligand binding (MEPP). The theoretical foundation of MEPP is the conserved electrostatic potential difference (EPD) in cognate pairs of active site residues in proteins with the same functionality. Previously, this invariance has been used to unravel promiscuous serine protease and metallo-β-lactamase scaffolds in alkaline phosphatases. Given that a similarity in EPD is significant, we expect differences in the EPD to be significant too. MEPP identifies residues or domains that undergo significant electrostatic perturbations, and also enumerates residue pairs that undergo significant polarity change. The gain in a certain polarity of a residue with respect to neighboring residues, or the reversal of polarity between two residues might indicate a change in the preferred ligand. The methodology of MEPP has been demonstrated on several enzymes that employ varying mechanisms to perform their roles. For example, we have attributed the change in polarity in residue pairs to be responsible for the loss of metal ion binding in fructose 1,6-bisphosphatases, and corroborated the pre-organized state of the active site of the enzyme with respect to functionally relevant changes in electric fields in ketosteroid isomerases.

Promiscuity-based enzyme selection for rational directed evolution experiments

Error-prone PCR, DNA shuffling, and saturation mutagenesis are techniques used by protein engineers to mimic the natural "evolutionary walk" that conjures new enzymes. Rational design is often critical in efforts to accelerate this "random walk" into a "resolute sprint." Previous work by our group established a computational method for detecting active sites (CLASP) based on spatial and electrostatic properties of catalytic residues, and a method to quantify promiscuous activities in a wide range of proteins (PROMISE). Here, we describe a rational design flow (DECAAF) based on the PROMISE methodology to choose a protein which, when subjected to minimal mutations, is most likely to mirror the scaffold of a desired enzymatic function. Modeling the diversity in catalytic sites and providing precise user control to guide the search is a key goal of our implementation. The flow details have been worked out in a real-life example to select a plant protein to substitute for human neutrophil elastase in a chimeric antimicrobial enzyme designed to bolster the innate immune defense system in plants.

The electrostatic profile of consecutive Cβ atoms applied to protein structure quality assessment

The structure of a protein provides insight into its physiological interactions with other components of the cellular soup. Methods that predict putative structures from sequences typically yield multiple, closely-ranked possibilities. A critical component in the process is the model quality assessing program (MQAP), which selects the best candidate from this pool of structures. Here, we present a novel MQAP based on the physical properties of sidechain atoms. We propose a method for assessing the quality of protein structures based on the electrostatic potential difference (EPD) of Cβ atoms in consecutive residues. We demonstrate that the EPDs of Cβ atoms on consecutive residues provide unique signatures of the amino acid types. The EPD of Cβ atoms are learnt from a set of 1000 non-homologous protein structures with a resolution cuto of 1.6 Å obtained from the PISCES database. Based on the Boltzmann hypothesis that lower energy conformations are proportionately sampled more, and on Annsen's thermodynamic hypothesis that the native structure of a protein is the minimum free energy state, we hypothesize that the deviation of observed EPD values from the mean values obtained in the learning phase is minimized in the native structure. We achieved an average specificity of 0.91, 0.94 and 0.93 on hg_structal, 4state_reduced and ig_structal decoy sets, respectively, taken from the Decoys `R' Us database. The source code and manual is made available at https://github.com/sanchak/mqap and permanently available on 10.5281/zenodo.7134.

Structural phylogeny by profile extraction and multiple superimposition using electrostatic congruence as a discriminator

Phylogenetic analysis of proteins using multiple sequence alignment (MSA) assumes an underlying evolutionary relationship in these proteins which occasionally remains undetected due to considerable sequence divergence. Structural alignment programs have been developed to unravel such fuzzy relationships. However, none of these structure based methods have used electrostatic properties to discriminate between spatially equivalent residues. We present a methodology for MSA of a set of related proteins with known structures using electrostatic properties as an additional discriminator (STEEP). STEEP first extracts a profile, then generates a multiple structural superimposition providing a consolidated spatial framework for comparing residues and finally emits the MSA. Residues that are aligned differently by including or excluding electrostatic properties can be targeted by directed evolution experiments to transform the enzymatic properties of one protein into another. We have compared STEEP results to those obtained from a MSA program (ClustalW) and a structural alignment method (MUSTANG) for chymotrypsin serine proteases. Subsequently, we used PhyML to generate phylogenetic trees for the serine and metallo-β-lactamase superfamilies from the STEEP generated MSA, and corroborated the accepted relationships in these superfamilies. We have observed that STEEP acts as a functional classifier when electrostatic congruence is used as a discriminator, and thus identifies potential targets for directed evolution experiments. In summary, STEEP is unique among phylogenetic methods for its ability to use electrostatic congruence to specify mutations that might be the source of the functional divergence in a protein family. Based on our results, we also hypothesize that the active site and its close vicinity contains enough information to infer the correct phylogeny for related proteins.

A measure of the broad substrate specificity of enzymes based on 'duplicate' catalytic residues

The ability of an enzyme to select and act upon a specific class of compounds with unerring precision and efficiency is an essential feature of life. Simultaneously, these enzymes often catalyze the reaction of a range of similar substrates of the same class, and also have promiscuous activities on unrelated substrates. Previously, we have established a methodology to quantify promiscuous activities in a wide range of proteins. In the current work, we quantitatively characterize the active site for the ability to catalyze distinct, yet related, substrates (BRASS). A protein with known structure and active site residues provides the framework for computing ‘duplicate’ residues, each of which results in slightly modified replicas of the active site scaffold. Such spatial congruence is supplemented by Finite difference Poisson Boltzmann analysis which filters out electrostatically unfavorable configurations. The congruent configurations are used to compute an index (BrassIndex), which reflects the broad substrate profile of the active site. We identify an acetylhydrolase and a methyltransferase as having the lowest and highest BrassIndex, respectively, from a set of non-homologous proteins extracted from the Catalytic Site Atlas. The acetylhydrolase, a regulatory enzyme, is known to be highly specific for platelet-activating factor. In the methyltransferase (PDB: 1QAM), various combinations of glycine (Gly38/40/42), asparagine (Asn101/11) and glutamic acid (Glu59/36) residues having similar spatial and electrostatic profiles with the specified scaffold (Gly38, Asn101 and Glu59) exemplifies the broad substrate profile such an active site may provide. ‘Duplicate’ residues identified by relaxing the spatial and/or electrostatic constraints can be the target of directed evolution methodologies, like saturation mutagenesis, for modulating the substrate specificity of proteins.

An automated flow for directed evolution based on detection of promiscuous scaffolds using spatial and electrostatic properties of catalytic residues

The aspiration to mimic and accelerate natural evolution has fueled interest in directed evolution experiments, which endow or enhance functionality in enzymes. Barring a few de novo approaches, most methods take a template protein having the desired activity, known active site residues and structure, and proceed to select a target protein which has a pre-existing scaffold congruent to the template motif. Previously, we have established a computational method (CLASP) based on spatial and electrostatic properties to detect active sites, and a method to quantify promiscuity in proteins. We exploit the prospect of promiscuous active sites to serve as the starting point for directed evolution and present a method to select a target protein which possesses a significant partial match with the template scaffold (DECAAF). A library of partial motifs, constructed from the active site residues of the template protein, is used to rank a set of target proteins based on maximal significant matches with the partial motifs, and cull out the best candidate from the reduced set as the target protein. Considering the scenario where this ‘incubator’ protein lacks activity, we identify mutations in the target protein that will mirror the template motif by superimposing the target and template protein based on the partial match. Using this superimposition technique, we analyzed the less than expected gain of activity achieved by an attempt to induce β-lactamase activity in a penicillin binding protein (PBP) (PBP-A from T. elongatus), and attributed this to steric hindrance from neighboring residues. We also propose mutations in PBP-5 from E. coli, which does not have similar steric constraints. The flow details have been worked out in an example which aims to select a substitute protein for human neutrophil elastase, preferably related to grapevines, in a chimeric anti-microbial enzyme which bolsters the innate immune defense system of grapevines.

Enumerating pathways of proton abstraction based on a spatial and electrostatic analysis of residues in the catalytic site

The pathways of proton abstraction (PA), a key aspect of most catalytic reactions, is often controversial and highly debated. Ultrahigh-resolution diffraction studies, molecular dynamics, quantum mechanics and molecular mechanic simulations are often adopted to gain insights in the PA mechanisms in enzymes. These methods require expertise and effort to setup and can be computationally intensive. We present a push button methodology--Proton abstraction Simulation (PRISM)--to enumerate the possible pathways of PA in a protein with known 3D structure based on the spatial and electrostatic properties of residues in the proximity of a given nucleophilic residue. Proton movements are evaluated in the vicinity of this nucleophilic residue based on distances, potential differences, spatial channels and characteristics of the individual residues (polarity, acidic, basic, etc). Modulating these parameters eliminates their empirical nature and also might reveal pathways that originate from conformational changes. We have validated our method using serine proteases and concurred with the dichotomy in PA in Class A β-lactamases, both of which are hydrolases. The PA mechanism in a transferase has also been corroborated. The source code is made available at www.sanchak.com/prism.

A measure of the promiscuity of proteins and characteristics of residues in the vicinity of the catalytic site that regulate promiscuity

Promiscuity, the basis for the evolution of new functions through 'tinkering' of residues in the vicinity of the catalytic site, is yet to be quantitatively defined. We present a computational method Promiscuity Indices Estimator (PROMISE)--based on signatures derived from the spatial and electrostatic properties of the catalytic residues, to estimate the promiscuity (PromIndex) of proteins with known active site residues and 3D structure. PromIndex reflects the number of different active site signatures that have congruent matches in close proximity of its native catalytic site, the quality of the matches and difference in the enzymatic activity. Promiscuity in proteins is observed to follow a lognormal distribution (μ = 0.28, σ = 1.1 reduced chi-square = 3.0E-5). The PROMISE predicted promiscuous functions in any protein can serve as the starting point for directed evolution experiments. PROMISE ranks carboxypeptidase A and ribonuclease A amongst the more promiscuous proteins. We have also investigated the properties of the residues in the vicinity of the catalytic site that regulates its promiscuity. Linear regression establishes a weak correlation (R(2)∼0.1) between certain properties of the residues (charge, polar, etc) in the neighborhood of the catalytic residues and PromIndex. A stronger relationship states that most proteins with high promiscuity have high percentages of charged and polar residues within a radius of 3 Å of the catalytic site, which is validated using one-tailed hypothesis tests (P-values∼0.05). Since it is known that these characteristics are key factors in catalysis, their relationship with the promiscuity index cross validates the methodology of PROMISE.