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

Proteases influence colony aggregation behavior in Vibrio cholerae

Aggregation behavior provides bacteria protection from harsh environments and threats to survival. Two uncharacterized proteases, LapX and Lap, are important for Vibrio cholerae liquid-based aggregation. Here, we determined that LapX is a serine protease with a preference for cleavage after glutamate and glutamine residues in the P1 position, which processes a physiologically based peptide substrate with a catalytic efficiency of 180 ± 80 M-1s-1. The activity with a LapX substrate identified by a multiplex substrate profiling by mass spectrometry screen was 590 ± 20 M-1s-1. Lap shares high sequence identity with an aminopeptidase (termed VpAP) from Vibrio proteolyticus and contains an inhibitory bacterial prepeptidase C-terminal domain that, when eliminated, increases catalytic efficiency on leucine p-nitroanilide nearly four-fold from 5.4 ± 4.1 × 104 M-1s-1 to 20.3 ± 4.3 × 104 M-1s-1. We demonstrate that LapX processes Lap to its mature form and thus amplifies Lap activity. The increase is approximately eighteen-fold for full-length Lap (95.7 ± 5.6 × 104 M-1s-1) and six-fold for Lap lacking the prepeptidase C-terminal domain (11.3 ± 1.9 × 105 M-1s-1). In addition, substrate profiling reveals preferences for these two proteases that could inform in vivo function. Furthermore, purified LapX and Lap restore the timing of the V. cholerae aggregation program to a mutant lacking the lapX and lap genes. Both proteases must be present to restore WT timing, and thus they appear to act sequentially: LapX acts on Lap, and Lap acts on the substrate involved in aggregation.

Biocontrol Potential of Bacillus subtilis and Bacillus tequilensis against Four Fusarium Species

The use of biological control agents as opposed to synthetic agrochemicals to control plant pathogens has gained momentum, considering their numerous advantages. The aim of this study is to investigate the biocontrol potential of plant bacterial isolates against Fusarium oxysporum, Fusarium proliferatum, Fusarium culmorum, and Fusarium verticillioides. Isolation, identification, characterization, and in vitro biocontrol antagonistic assays of these isolates against Fusarium species were carried out following standard protocols. The bacterial endophytes were isolated from Glycine max. L leaves (B1), Brassica napus. L seeds (B2), Vigna unguiculata seeds (B3), and Glycine max. L seeds (B4). The bacterial isolates were identified using 16S rRNA PCR sequencing. A phylogenetic analysis shows that the bacterial isolates are closely related to Bacillus subtilis (B1) and Bacillus tequilensis (B2-B4), with an identity score above 98%. All the bacterial isolates produced a significant amount (p < 0.05) of indole acetic acid (IAA), siderophores, and protease activity. In vitro antagonistic assays of these isolates show a significant (p < 0.05) growth inhibition of the fungal mycelia in the following order: F. proliferatum > F. culmorum > F. verticillioides > F. oxysporum, compared to the control. The results suggest that these bacterial isolates are good biocontrol candidates against the selected Fusarium species.

Crystal structure of vespid phospholipase A(1) reveals insights into the mechanism for cause of membrane dysfunction

Vespid phospholipase A1 (vPLA1) from the black-bellied hornet (Vespa basalis) catalyzes the hydrolysis of emulsified phospholipids and shows potent hemolytic activity that is responsible for its lethal effect. To investigate the mechanism of vPLA1 towards its function such as hemolysis and emulsification, we isolated vPLA1 from V. basalis venom and determined its crystal structure at 2.5 Å resolution. vPLA1 belongs to the α/β hydrolase fold family. It contains a tightly packed β-sheet surrounded by ten α-helices and a Gly-X-Ser-X-Gly motif, characteristic of a serine hydrolyase active site. A bound phospholipid was modeled into the active site adjacent to the catalytic Ser-His-Asp triad indicating that Gln95 is located at hydrogen-bonding distance from the substrate's phosphate group. Moreover, a hydrophobic surface comprised by the side chains of Phe53, Phe62, Met91, Tyr99, Leu197, Ala167 and Pro169 may serve as the acyl chain-binding site. vPLA1 shows global similarity to the N-terminal domain of human pancreatic lipase (HPL), but with some local differences. The lid domain and the β9 loop responsible for substrate selectivity in vPLA1 are shorter than in HPL. Thus, solvent-exposed hydrophilic residues can easily accommodate the polar head groups of phospholipids, thereby accounting for the high activity level of vPLA1. Our result provides a potential explanation for the ability of vPLA1 to hydrolyze phospholipids of cell membrane.

Amsacta moorei entomopoxvirus encodes a functional esterase (amv133) with protease activity

OBJECTIVES

Lipolytic genes have been investigated in several viral genomes, and some of them show enzyme activity which can be used for various functions including the production of DNA replication metabolites, rescue from endosomes, and membrane fusion. Amsacta moorei entomopoxvirus (AMEV) replicates in nearly the entire insect body, especially in the adipose tissue. One of the open reading frames (ORFs) in the AMEV genome, amv133, encodes a putative lipase enzyme. In this study, we therefore investigate the enzyme activity of amv133.

METHODS

amv133 was aligned with known lipase genes and their homologs in entomopoxviruses. Expressed proteins were partially purified and assayed for lipase, esterase and protease.

RESULTS

We found that amv133 contains all the domains required for a functional lipase enzyme and that it shows a significant similarity with homologs in other entomopoxviruses. Since there is a similarity of the catalytic triad between lipases and serine proteases, we also investigated the protease activity of amv133. Lipase, esterase and protease assays showed that amv133 encodes a functional esterase enzyme with protease activity.

CONCLUSION

The current data show that amv133 is a conserved gene in all entomopoxvirus genomes sequenced so far and might contribute greatly to degrading the lipids or proteins and hence improve the virus infection.

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.

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.

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.