Structural Characterization of Functionally Important Chloride Binding Sites in the Marine Vibrio Alkaline Phosphatase.

Enzyme stability and function can be affected by various environmental factors, such as temperature, pH, and ionic strength. Enzymes that are located outside the relatively unchanging environment of the cytosol, such as those residing in the periplasmic space of bacteria or extracellularly secreted, are challenged by more fluctuations in the aqueous medium. Bacterial alkaline phosphatases (APs) are generally affected by ionic strength of the medium, but this varies substantially between species. An AP from the marine bacterium Vibrio splendidus (VAP) shows complex pH-dependent activation and stabilization in the 0-1.0 M range of halogen salts and has been hypothesized to specifically bind chloride anions. Here, using X-ray crystallography and anomalous scattering, we have located two chloride binding sites in the structure of VAP, one in the active site and another one at a peripheral site. Further characterization of the binding sites using site-directed mutagenesis and small-angle X-ray scattering showed that upon binding of chloride to the peripheral site, structural dynamics decreased locally, resulting in thermal stabilization of the VAP active conformation. Binding of the chloride ion in the active site did not displace the bound inorganic phosphate product, but it may promote product release by facilitating rotational stabilization of the substrate-binding Arg129. Overall, these results reveal the complex nature and dynamics of chloride binding to enzymes through long-range modulation of electronic potential in the vicinity of the active site, resulting in increased catalytic efficiency and stability.

pH-Dependent Binding of Chloride to a Marine Alkaline Phosphatase Affects the Catalysis, Active Site Stability, and Dimer Equilibrium.

The effect of ionic strength on enzyme activity and stability varies considerably between enzymes. Ionic strength is known to affect the catalytic activity of some alkaline phosphatases (APs), such as Escherichia coli AP, but how ions affect APs is debated. Here, we studied the effect of various ions on a cold-adapted AP from Vibrio splendidus (VAP). Previously, we have found that the active form of VAP is extremely unstable at low ionic strengths. Here we show that NaCl increased the activity and stability of VAP and that the effect was pH-dependent in the range of pH 7-10. The activity profile as a function of pH formed two maxima, indicating a possible conformational change. Bringing the pH from the neutral to the alkaline range was accompanied by a large increase in both the Ki for inorganic phosphate (product inhibition) and the KM for p-nitrophenyl phosphate. The activity transitions observed as the pH was varied correlated with structural changes as monitored by tryptophan fluorescence. Thermal and urea-induced inactivation was shown to be accompanied by neither dissociation of the active site metal ions nor dimer dissociation. This would suggest that the inactivation involved subtle changes in active site conformation. Furthermore, the VAP dimer equilibrium was studied for the first time and shown to highly favor dimerization, which was dependent on pH and NaCl concentration. Taken together, the data support a model in which anions bind to some specific acceptor in the active site of VAP, resulting in great stabilization and catalytic rate enhancement, presumably through a different mechanism.

Cold-active alkaline phosphatase is irreversibly transformed into an inactive dimer by low urea concentrations.

Alkaline phosphatase is a homodimeric metallo-hydrolase where both Zn(2+) and Mg(2+) are important for catalysis and stability. Cold-adapted alkaline phosphatase variants have high activity at low temperatures and lower thermal stability compared with variants from mesophilic hosts. The instability, and thus inactivation, could be due to loose association of the dimers and/or loosely bound Mg(2)(+) in the active site, but this has not been studied in detail for the cold-adapted variants. Here, we focus on using the intrinsic fluorescence of Trp in alkaline phosphatase from the marine bacterium Vibrio splendidus (VAP) to probe for dimerization. Trp→Phe substitutions showed that two out of the five native Trp residues contributed mostly to the fluorescence emission. One residue, 15Šaway from the active site (W460) and highly solvent excluded, was phosphorescent and had a distant role in substrate binding. An additional Trp residue was introduced to the dimer interface to act as a possible probe for dimerization. Urea denaturation curves indicated that an inactive dimer intermediate, structurally equivalent to the native state, was formed before dimer dissociation took place. This is the first example of the transition of a native dimer to an inactive dimer intermediate for alkaline phosphatase without using mutagenesis, ligands, or competitive inhibition.

Dynamics fingerprint and inherent asymmetric flexibility of a cold-adapted homodimeric enzyme. A case study of the Vibrio alkaline phosphatase.

BACKGROUND: Protein dynamics influence protein function and stability and modulate conformational changes. Such motions depend on the underlying networks of intramolecular interactions and communicating residues within the protein structure. Here, we provide the first characterization of the dynamic fingerprint of the dimeric alkaline phosphatase (AP) from the cold-adapted Vibrio strain G15-21 (VAP), which is among the APs with the highest known kcat at low temperatures. METHODS: Multiple all-atom explicit solvent molecular dynamics simulations were employed in conjunction with different metrics to analyze the dynamical patterns and the paths of intra- and intermolecular communication. RESULTS: Interactions and coupled motions at the interface between the two VAP subunits have been characterized, along with the networks of intramolecular interactions. It turns out a low number of intermolecular interactions and coupled motions, which result differently distributed in the two monomers. The paths of long-range communication mediated from the catalytic residues to distal sites were also characterized, pointing out a different information flow in the two subunits. CONCLUSIONS: A pattern of asymmetric flexibility is evident in the two identical subunits of the VAP dimer that is intimately linked to a different distribution of intra- and intermolecular interactions. The asymmetry was also evident in pairs of cross-correlated residues during the dynamics. GENERAL SIGNIFICANCE: The results here discussed provide a structural rationale to the half-of-site mechanism previously proposed for VAP and other APs, as well as a general framework to characterize asymmetric dynamics in homomeric enzymes.

The high catalytic rate of the cold-active Vibrio alkaline phosphatase requires a hydrogen bonding network involving a large interface loop.

The role of surface loops in mediating communication through residue networks is still a relatively poorly understood part in the study of cold adaptation of enzymes, especially in terms of their quaternary interactions. Alkaline phosphatase (AP) from the psychrophilic marine bacterium Vibrio splendidus (VAP) is characterized by an analogous large surface loop in each monomer, referred to as the large loop, that hovers over the active site of the other monomer. It presumably has a role in the high catalytic efficiency of VAP which accompanies its extremely low thermal stability. Here, we designed several different variants of VAP with the aim of removing intersubunit interactions at the dimer interface. Breaking the intersubunit contacts from one residue in particular (Arg336) reduced the temperature stability of the catalytically potent conformation and caused a 40% drop in catalytic rate. The high catalytic rates of enzymes from cold-adapted organisms are often associated with increased dynamic flexibility. Comparison of the relative B-factors of the R336L crystal structure to that of the wild-type confirmed surface flexibility was increased in a loop on the opposite monomer, but not in the large loop. The increase in flexibility resulted in a reduced catalytic rate. The large loop increases the area of the interface between the subunits through its contacts and may facilitate an alternating structural cycle demanded by a half-of-sites reaction mechanism through stronger ties, as the dimer oscillates between high affinity (active) or low phosphoryl group affinity (inactive).

The 1.4 A crystal structure of the large and cold-active Vibrio sp. alkaline phosphatase.

Alkaline phosphatase (AP) from the cold-adapted Vibrio strain G15-21 is among the AP variants with the highest known k(cat) value. Here the structure of the enzyme at 1.4 A resolution is reported and compared to APs from E. coli, human placenta, shrimp and the Antarctic bacterium strain TAB5. The Vibrio AP is a dimer although its monomers are without the long N-terminal helix that embraces the other subunit in many other APs. The long insertion loop, previously noted as a special feature of the Vibrio AP, serves a similar function. The surface does not have the high negative charge density as observed in shrimp AP, but a positively charged patch is observed around the active site that may be favourable for substrate binding. The dimer interface has a similar number of non-covalent interactions as other APs and the "crown"-domain is the largest observed in known APs. Part of it slopes over the catalytic site suggesting that the substrates may be small molecules. The catalytic serines are refined with multiple conformations in both monomers. One of the ligands to the catalytic zinc ion in binding site M1 is directly connected to the crown-domain and is closest to the dimer interface. Subtle movements in metal ligands may help in the release of the product and/or facilitate prior dephosphorylation of the covalent intermediate. Intersubunit interactions may be a major factor for promoting active site geometries that lead to the high catalytic activity of Vibrio AP at low temperatures.

X-ray crystal structure of Vibrio alkaline phosphatase with the non-competitive inhibitor cyclohexylamine.

BACKGROUND: Para-nitrophenyl phosphate, the common substrate for alkaline phosphatase (AP), is available as a cyclohexylamine salt. Here, we report that cyclohexylamine is a non-competitive inhibitor of APs. METHODS: Cyclohexylamine inhibited four different APs. Co-crystallization with the cold-active Vibrio AP (VAP) was performed and the structure solved. RESULTS: Inhibition of VAP fitted a non-competitive kinetic model (Km unchanged, Vmax reduced) with IC50 45.3 mM at the pH optimum 9.8, not sensitive to 0.5 M NaCl, and IC50 27.9 mM at pH 8.0, where the addition of 0.5 M NaCl altered the inhibition to the level observed at pH 9.8. APs from E. coli and calf intestines were less sensitive to cyclohexylamine, whereas an Antarctic bacterial AP was similar to VAP in this respect. X-ray crystallography at 2.3 Å showed two binding sites, one in the active site channel and another at the surface close to dimer interface. Antarctic bacterial AP and VAP have Trp274 in common in their active-sites, that takes part in binding cyclohexylamine. VAP variants W274A, W274K, and W274H gave IC50 values of 179 mM, 188 mM and 187 mM, respectively, at pH 9.8. CONCLUSIONS: The binding of cyclohexylamine in locations at the dimeric interface and/or in the active site of APs may delay product release or reduce the rate of catalytic step(s) involving conformational changes and intersubunit communications. GENERAL SIGNIFICANCE: Cyclohexylamine is a common chemical in industries and used as a counterion in substrates for alkaline phosphatase, a clinically important and common enzyme in the biosphere.

Chloride promotes refolding of active Vibrio alkaline phosphatase through an inactive dimeric intermediate with an altered interface.

Most enzymes are homodimers or higher order multimers. Cold-active alkaline phosphatase from Vibrio splendidus (VAP) transitions into a dimer with very low activity under mild denaturation conditions. The desire to understand why this dimer fails to efficiently catalyse phosphomonoester hydrolysis led us to investigate interfacial communication between subunits. Here, we studied in detail the unfolding mechanism at two pH values and in the presence or absence of sodium chloride. At pH 8.0, the denaturation model had to include an inactive dimer intermediate and follow the pathway: N2 → I2 → 2U. At pH 10.5, the model was of a two-state nature. Enzyme activity was not recovered under several examined refolding conditions. However, in the presence of 0.5 m NaCl, the enzyme was nearly fully reactivated after urea treatment. Thermal inactivation experiments were biphasic where the inactivation could be detected using CD spectroscopy at 190-200 nm. By incorporating a bimane fluorescence probe at the dimer interface, we could monitor inactivation/denaturation at two distinct sites at the dimer interface. A change in bimane fluorescence at both sites was observed during inactivation, but prior to the global unfolding event. Furthermore, the rate of change in bimane fluorescence correlated with inactivation rates at 40 °C. These results indicate and support the hypothesis that the subunits of VAP are only functional in the dimeric state due to the cooperative nature of the reaction mechanism when proper crosstalk between subunits is facilitated.

Engineered disulfide bonds increase active-site local stability and reduce catalytic activity of a cold-adapted alkaline phosphatase.

Alkaline phosphatase is an extracellular enzyme that is membrane-bound in eukaryotes but resides in the periplasmic space of bacteria. It normally carries four cysteine residues that form two disulfide bonds, for instance in the APs of Escherichia coli and vertebrates. An AP variant from a Vibrio sp. has only one cysteine residue. This cysteine is second next to the nucleophilic serine in the active site. We have individually modified seven residues to cysteine that are on two loops predicted to be within a 5 A radius. Four of them formed a disulfide bond to the endogenous cysteine. Thermal stability was monitored by circular dichroism and activity measurements. Global stability was similar to the wild-type enzyme. However, a significant increase in heat-stability was observed for the disulfide-containing variants using activity as a measure, together with a large reduction in catalytic rates (k(cat)) and a general decrease in Km values. The results suggest that a high degree of mobility near the active site and in the helix carrying the endogenous cysteine is essential for full catalytic efficiency in the cold-adapted AP.

Effects of replacing active site residues in a cold-active alkaline phosphatase with those found in its mesophilic counterpart from Escherichia coli.

Alkaline phosphatase (AP) from a North Atlantic marine Vibrio bacterium was previously characterized as being kinetically cold-adapted. It is still unknown whether its characteristics originate locally in the active site or are linked to more general structural factors. There are three metal-binding sites in the active site of APs, and all three metal ions participate in catalysis. The amino acid residues that bind the two zinc ions most commonly present are conserved in all known APs. In contrast, two of the residues that bind the third metal ion (numbered 153 and 328 in Escherichia coli AP) are different in various APs. This may explain their different catalytic efficiencies, as the Mg2+ most often present there is important for both structural stability and the reaction mechanism. We have mutated these key residues to the corresponding residues in E. coli AP to obtain the double mutant Asp116/Lys274, and both single mutants. All these mutants displayed reduced substrate affinity and lower overall reaction rates. The Lys274 and Asp116/Lys274 mutants also displayed an increase in global heat stability, which may be due to the formation of a stabilizing salt bridge. Overall, the results show that a single amino acid substitution in the active site is sufficient to alter the structural stability of the cold-active Vibrio AP both locally and globally, and this influences kinetic properties.

Reversible inactivation of alkaline phosphatase from Atlantic cod (Gadus morhua) in urea.

Alkaline phosphatase (AP) from Atlantic cod (Gadus morhua) is a zinc and magnesium containing homodimer that requires the oligomeric state for activity. Its kinetic properties are indicative of cold-adaptation. Here, the effect of urea on the structural stability was studied in order to correlate the activity with metal content, the microenvironment around tryptophan residues, and events at the subunit interface. At the lowest concentrations of urea, the first detected alteration in properties was an increase in the activity of the enzyme. This was followed by inactivation, and the release of half of the zinc content when the amount of urea reached levels of 2 M. Intrinsic tryptophan fluorescence and circular dichroism ellipticity changed in the range 2.5 to 8 M urea, signaling dimer dissociation, followed by one major monomer unfolding transition at 6-8 M urea as indicated by ANS fluorescence and KI fluorescence quenching. Gibbs free energy was estimated by the linear extrapolation method using a three-state model as 8.6 kcal/mol for dimer stability and 11.6 kcal/mol for monomer unfolding giving a total of 31.8 kcal/mol. Dimer association had a very small ionic contribution. Dimers were stable in relatively high concentration of urea, whereas the immediate vicinity around the active site was vulnerable to low concentrations of urea. Thus, inactivation did not coincide with dimer dissociation, suggesting that the active site is the most dynamic part of the molecule and closest related to cold-adaptation of its enzymatic activity.

Microscopic rate-constants for substrate binding and acylation in cold-adaptation of trypsin I from Atlantic cod.

Temperature imposes limits on where life can thrive and this is evident in the evolution of the basic structural properties of proteins. Cold-adaptation of enzymes is one example, where the catalytic rate constant (k(cat)) is increased compared with hot-acclimated homologous under identical assay conditions. Trypsin I from Atlantic cod (Gadus morhua) has catalytic efficiency (k(cat)/K(m)) for amide hydrolysis that is 17-fold larger than observed for bovine trypsin. Here, the individual rate-constants for association of substrate (k(1)), dissociation of substrate (k(-1)), and acylation of the enzyme (k(2)) have been determined using benzoyl-Arg-p-nitroanilide or benzyloxycarbonyl-Gly-Pro-Arg-p-nitroanilide as substrates. Rather unexpectedly, by far the largest difference (37-fold increase) was observed in k(1), the rate constant for binding of substrate. The cold-adaptation of the dissociation and catalytic steps were not as prominent (increased by 3.7-fold). The length of substrate did have an effect by increasing the reaction rate by 70-fold, and again, the step most affected was the initial binding-step.

Amino acid sequence of the cold-active alkaline phosphatase from Atlantic cod (Gadus morhua).

Atlantic cod is a marine fish that lives at low temperatures of 0-10 degrees C and contains a cold-adapted alkaline phosphatase (AP). Preparations of AP from either the lower part of the intestines or the pyloric caeca area were subjected to proteolytic digestion, mass spectrometry and amino acid sequencing by Edman degradation. The primary structure exhibits greatest similarity to human tissue non-specific AP (80%), and approximately 30% similarity to AP from Escherichia coli. The key residues required for catalysis are conserved in the cod AP, except for the third metal binding site, where cod AP has the same variable residues as mammalian APs (His153 and His328 by E. coli AP numbering). General comparison of the amino acid composition with mammalian APs showed that cod AP contains fewer Cys, Leu, Met and Ser, but proportionally more Asn, Asp, Ile, Lys, Trp and Tyr residues. Three N-linked glycosylation sites were found. The glycan structure was determined as complex biantennary in type with fucose and sialic acid attached, although a trace of complex tri-antennary structure was also observed. A three-dimensional model was obtained by homology modelling using the human placental AP scaffold. Cod AP has fewer charged and hydrophobic residues, but more polar residues at the intersubunit surface. The N-terminal helix arm that embraces the second subunit in dimeric APs may be more flexible due to a replaced Pro at its base. One disulfide bridge was found instead of the two present in most other APs. This may invoke greater movement in the structure that together with weaker subunit contacts leads to improved catalytic efficiency.

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.

Characterizing alpha helical properties of Ebola viral proteins as potential targets for inhibition of alpha-helix mediated protein-protein interactions.

Ebola, considered till recently as a rare and endemic disease, has dramatically transformed into a potentially global humanitarian crisis. The genome of Ebola, a member of the Filoviridae family, encodes seven proteins. Based on the recently implemented software (PAGAL) for analyzing the hydrophobicity and amphipathicity properties of alpha helices (AH) in proteins, we characterize the helices in the Ebola proteome. We demonstrate that AHs with characteristically unique features are involved in critical interactions with the host proteins. For example, the Ebola virus membrane fusion subunit, GP2, from the envelope glycoprotein ectodomain has an AH with a large hydrophobic moment. The ability of this AH to bind to other host proteins is disrupted by a neutralizing antibody derived from a human survivor of the 1995 Kikwit outbreak, emphasizing the critical nature of this secondary structure in the virulence of the Ebola virus. Our method ensures a comprehensive list of such `hotspots'. These helices probably are or can be the target of molecules designed to inhibit AH mediated protein-protein interactions. Further, by comparing the AHs in proteins of the related Marburg viruses, we are able to elicit subtle changes in the proteins that might render them ineffective to previously successful drugs. Such differences are difficult to identify by a simple sequence or structural alignment. Thus, analyzing AHs in the small Ebola proteome can aid rational design aimed at countering the `largest Ebola epidemic, affecting multiple countries in West Africa' ( http://www.cdc.gov/vhf/ebola/outbreaks/2014-west-africa/index.html).

The PDB database is a rich source of alpha-helical anti-microbial peptides to combat disease causing pathogens.

The therapeutic potential of α-helical anti-microbial peptides (AH-AMP) to combat pathogens is fast gaining prominence. Based on recently published open access software for characterizing α-helical peptides (PAGAL), we elucidate a search methodology (SCALPEL) that leverages the massive structural data pre-existing in the PDB database to obtain AH-AMPs belonging to the host proteome. We provide in vitro validation of SCALPEL on plant pathogens ( Xylella fastidiosa, Xanthomonas arboricola and Liberibacter crescens) by identifying AH-AMPs that mirror the function and properties of cecropin B, a well-studied AH-AMP. The identified peptides include a linear AH-AMP present within the existing structure of phosphoenolpyruvate carboxylase (PPC20), and an AH-AMP mimicing the properties of the two α-helices of cecropin B from chitinase (CHITI25). The minimum inhibitory concentration of these peptides are comparable to that of cecropin B, while anionic peptides used as control failed to show any inhibitory effect on these pathogens. Substitute therapies in place of conventional chemotherapies using membrane permeabilizing peptides like these might also prove effective to target cancer cells. The use of native structures from the same organism could possibly ensure that administration of such peptides will be better tolerated and not elicit an adverse immune response. We suggest a similar approach to target Ebola epitopes, enumerated using PAGAL recently, by selecting suitable peptides from the human proteome, especially in wake of recent reports of cationic amphiphiles inhibiting virus entry and infection.

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.

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.

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.