Functional interrelationships in the alkaline phosphatase superfamily: phosphodiesterase activity of Escherichia coli alkaline phosphatase

The microbial phosphorus cycle in aquatic ecosystems

Phosphorus is an essential element for life, and phosphorus cycling is crucial to planetary habitability. In aquatic environments, microorganisms are a major component of phosphorus cycling and rapidly transform the diverse chemical forms of phosphorus through various uptake, assimilation and release pathways. Recent discoveries have revealed a more dynamic and complex aquatic microbial phosphorus cycle than previously understood. Some microorganisms have been shown to use and produce new phosphorus compounds, including those in reduced forms. New findings have also raised numerous unanswered questions that warrant further investigation. There is an expanding influence of human activity on aquatic ecosystems. Advancements in understanding the phosphorus biogeochemistry of evolving aquatic environments offer a unique opportunity to comprehend, anticipate and mitigate the effect of human activities. In this Review, I discuss the wealth of new aquatic phosphorus cycle research, spanning disciplines from omics and physiology to global biogeochemical modelling, with a focus on the current comprehension of how aquatic microorganisms sense, transport, assimilate, store, produce and release phosphorus. Of note, I delve into cellular phosphorus allocation, an underexplored topic with wide-ranging implications for energy and element flux in aquatic ecosystems.

Selection of a promiscuous minimalist cAMP phosphodiesterase from a library of de novo designed proteins

The ability of unevolved amino acid sequences to become biological catalysts was key to the emergence of life on Earth. However, billions of years of evolution separate complex modern enzymes from their simpler early ancestors. To probe how unevolved sequences can develop new functions, we use ultrahigh-throughput droplet microfluidics to screen for phosphoesterase activity amidst a library of more than one million sequences based on a de novo designed 4-helix bundle. Characterization of hits revealed that acquisition of function involved a large jump in sequence space enriching for truncations that removed >40% of the protein chain. Biophysical characterization of a catalytically active truncated protein revealed that it dimerizes into an α-helical structure, with the gain of function accompanied by increased structural dynamics. The identified phosphodiesterase is a manganese-dependent metalloenzyme that hydrolyses a range of phosphodiesters. It is most active towards cyclic AMP, with a rate acceleration of ~109 and a catalytic proficiency of >1014 M-1, comparable to larger enzymes shaped by billions of years of evolution.

Genetic Perturbation Alters Functional Substates in Alkaline Phosphatase

Enzymes inherently exhibit molecule-to-molecule heterogeneity in their conformational and functional states, which is considered to be a key to the evolution of new functions. Single-molecule enzyme assays enable us to directly observe such multiple functional states or functional substates. Here, we quantitatively analyzed functional substates in the wild-type and 69 single-point mutants of Escherichia coli alkaline phosphatase by employing a high-throughput single-molecule assay with a femtoliter reactor array device. Interestingly, many mutant enzymes exhibited significantly heterogeneous functional substates with various types, while the wild-type enzyme showed a highly homogeneous substate. We identified a correlation between the degree of functional substates and the level of improvement in promiscuous activities. Our work provides much comprehensive evidence that the functional substates can be easily altered by mutations, and the evolution toward a new catalytic activity may involve the modulation of the functional substates.

Charge-reversal nanoemulsions: A systematic investigation of phosphorylated PEG-based surfactants

Surfactants bearing monophosphate esters with PEG of increasing chain length and different lipophilic tail structures were investigated to improve the effectiveness of enzyme triggered charge-converting nanoemulsions. The surfactants PEG-8-stearate, PEG-22-tocopheryl succinate (TPGS), PEG-3-oleate, PEG-9-oleate and PEG-9-lauryl ether were phosphorylated and incorporated in a self-emulsifying drug delivery system (SEDDS) exhibiting a defined PEG corona. To provide a positive zeta potential increasing amounts of the cationic surfactant benzalkonium chloride (BA) were incorporated. The effect of these PEG monophosphate esters (P-PEG-surfactants) was evaluated based on enzyme induced phosphate release and change in zeta potential. Significant enzyme induced charge conversion was observed for all P-PEG-surfactants, showing shifts from Δ3 mV to Δ31 mV. Surfactants comprising the shortest and longest PEG chain showed similar amplitudes (P-PEG-3-oleate: Δ11.9 mV; P-PEG-22-TPGS Δ10.2 mV), whereas P-PEG-8-stearate, P-PEG-9-oleate and P-PEG-9-lauryl ether bearing similarly long PEG chains but different lipophilic tail structures resulted in pronounced differences in amplitudes of Δ10.3 mV, Δ14.5 mV and Δ18.1 mV, respectively. Furthermore, an indirect correlation between the lipophilicity of P-PEG-surfactants and the obtained charge-reversing effect was observed. With the exception of P-PEG-lauryl ether, this charge-reversal effect decreased with increasing BA concentrations. In conclusion, the enzyme induced amplitude of charge conversion of P-PEG-surfactants depends to a high extent on their lipophilic tail structure. Based on this knowledge potent charge-reversal nanoemulsions can be designed.

High throughput and quantitative enzymology in the genomic era

Accurate predictions from models based on physical principles are the ultimate metric of our biophysical understanding. Although there has been stunning progress toward structure prediction, quantitative prediction of enzyme function has remained challenging. Realizing this goal will require large numbers of quantitative measurements of rate and binding constants and the use of these ground-truth data sets to guide the development and testing of these quantitative models. Ground truth data more closely linked to the underlying physical forces are also desired. Here, we describe technological advances that enable both types of ground truth measurements. These advances allow classic models to be tested, provide novel mechanistic insights, and place us on the path toward a predictive understanding of enzyme structure and function.

Inhibition of the activated sludge-associated enzyme phosphatase by transition metal oxyanions

Organic esters of phosphoric acid and other organophosphorous compounds are enzymatically hydrolyzed during wastewater treatment by microbial phosphoesterases, especially by phosphomonoesterase (phosphatase). For physiological reasons, the enzyme is inhibited by its main inorganic reaction product, ortho-phosphate. It is known that oxyanions of transition metals, resembling the molecular topology of ortho-phosphate, e.g. vanadate and tungstate, are more potent inhibitors for microbial alkaline phosphatase than phosphate. To proof this effect for activated sludge, a multitude of samples from a communal wastewater treatment plant was exposed at pH values from 7.00 to 8.50 to tungstate, vanadate, and molybdate. Inhibition effects were determined by a sensitive fluorimetric microplate assay and characteristic parameters (IC₅₀ and IC₂₀ concentrations) were deduced from modelled dose-response functions. Mean inhibitor concentrations (in brackets: ranges) causing 50% inactivation (IC₅₀) at pH 7.50 were 2.5 (1.3-4.1) μM tungstate, 2.9 (1.6-5.5) μM vanadate, and 41.4 (33.6-56.7) μM molybdate. Vanadate and tungstate concentrations between 0.6 and 0.7 μM provoked a 20% (IC₂₀) inhibition. The inhibition efficiency of tungstate and molybdate decreased with increasing pH, whereas vanadate reacted pH independently. These results underline the necessity to consider enzyme inhibition assessing the limitations and potentials of biological wastewater treatment processes.

Mechanism and Inhibition of Human Methionine Adenosyltransferase 2A

S-Adenosyl-l-methionine (AdoMet) is synthesized by the MAT2A isozyme of methionine adenosyltransferase in most human tissues and in cancers. Its contribution to epigenetic control has made it a target for anticancer intervention. A recent kinetic isotope effect analysis of MAT2A demonstrated a loose nucleophilic transition state. Here we show that MAT2A has a sequential mechanism with a rate-limiting step of formation of AdoMet, followed by rapid hydrolysis of the β-γ bond of triphosphate, and rapid release of phosphate and pyrophosphate. MAT2A catalyzes the slow hydrolysis of both ATP and triphosphate in the absence of other reactants. Positional isotope exchange occurs with ¹⁸O as the 5'-oxygen of ATP. Loss of the triphosphate is sufficiently reversible to permit rotation and recombination of the α-phosphoryl group of ATP. Adenosine (α-β or β-γ)-imido triphosphates are slow substrates, and the respective imido triphosphates are inhibitors. The hydrolytically stable (α-β, β-γ)-diimido triphosphate (PNPNP) is a nanomolar inhibitor. The MAT2A protein structure is highly stabilized against denaturation by binding of PNPNP. A crystal structure of MAT2A with 5'-methylthioadenosine and PNPNP shows the ligands arranged appropriately in the ATP binding site. Two magnesium ions chelate the α- and γ-phosphoryl groups of PNPNP. The β-phosphoryl oxygen is in contact with an essential potassium ion. Imidophosphate derivatives provide contact models for the design of catalytic site ligands for MAT2A.

Origin of the Phosphoprotein Phosphatase (PPP) sequence family in Bacteria: Critical ancestral sequence changes, radiation patterns and substrate binding features

Background

Phosphoprotein phosphatases (PPP) belong to the PPP Sequence family, which in turn belongs to the broader metallophosphoesterase (MPE) superfamily. The relationship between the PPP Sequence family and other members of the MPE superfamily remains unresolved, in particular what transitions took place in an ancestral MPE to ultimately produce the phosphoprotein specific phosphatases (PPPs).

Methods

We use structural and sequence alignment data, phylogenetic tree analysis, sequence signature (Weblogo) analysis, in silico protein-peptide modeling data, and in silico mutagenesis to trace a likely route of evolution from MPEs to the PPP Sequence family. Hidden Markov Model (HMM) based iterative database search strategies were utilized to identify PPP Sequence Family members from numerous bacterial groups.

Results

Using Mre11 as proxy for an ancestral nuclease-like MPE we trace a possible evolutionary route that alters a single active site substrate binding His-residue to yield a new substrate binding accessory, the "2-Arg-Clamp". The 2-Arg-Clamp is not found in MPEs, but is present in all PPP Sequence family members, where the phosphomonesterase reaction predominates. Variation in position of the clamp arginines and a supplemental sequence loop likely provide substrate specificity for each PPP Sequence family group.

Conclusions

Loss of a key substrate binding His-in MPEs opened the path to bind novel substrates and evolution of the 2-Arg-Clamp, a sequence change seen in both bacterial and eukaryotic phosphoprotein phosphatases.General significance: We establish a likely evolutionary route from nuclease-like MPE to PPP Sequence family enzymes, that includes the phosphoprotein phosphatases.

Discovery of the First Human Arylsulfatase A Reversible Inhibitor Impairing Mouse Oocyte Fertilization

Arylsulfatase A (ARSA) plays a crucial role in the reproduction of mammals due to its involvement in the specific gamete interaction preceding sperm and egg fusion leading to fertilization. Recently, it has been shown that zona pellucida (ZP) sperm binding and in vivo fertilization in mice are markedly hampered by using a specific anti-ARSA antibody. Herein, the design and discovery of the first ARSA small molecule inhibitor based on a coumarin-containing polycycle are presented. Through a structure-based approach applied on our in-house library, compound 1r was identified as an ARSA reversible inhibitor (ARSAi); then its activity was validated through both surface plasmon resonance and biochemical inhibition experiments, the first providing a K_D value of 21 μM and the latter an IC₅₀ value of 13.2 μM. Further investigations highlighted that compound 1r induced 20% sperm death at 25 μM and also impaired sperm motility; nevertheless both the effects were mediated by ROS production, since they were rescued by the cotreatment of 1r and N-acetyl cysteine (NAC). Interestingly, while 1r was not able to hamper the ZP/sperm binding, it markedly decreased the in vitro oocyte fertilization by mouse sperm up to 60%. Notably, this effect was not hampered by 1r/NAC coadministration, hence allowing the ruling out of an ROS-dependent mechanism. In conclusion, herein is reported the first ever hit of ARSAi as a chemical tool that will enable better exploration of ARSA's biological role in fertilization as well as provide a starting point for developing 1r structure optimization aimed at increasing enzyme inhibition potency but also providing a deeper understanding of the involvement of ARSA in the fertilization pathway mechanism.

The physical basis and practical consequences of biological promiscuity

Proteins interact with metabolites, nucleic acids, and other proteins to orchestrate the myriad catalytic, structural and regulatory functions that support life from the simplest microbes to the most complex multicellular organisms. These molecular interactions are often exquisitely specific, but never perfectly so. Adventitious "promiscuous" interactions are ubiquitous due to the thousands of macromolecules and small molecules crowded together in cells. Such interactions may perturb protein function at the molecular level, but as long as they do not compromise organismal fitness, they will not be removed by natural selection. Although promiscuous interactions are physiologically irrelevant, they are important because they can provide a vast reservoir of potential functions that can provide the starting point for evolution of new functions, both in nature and in the laboratory.

A Novel Alkaline Phosphatase/Phosphodiesterase, CamPhoD, from Marine Bacterium Cobetia amphilecti KMM 296

A novel extracellular alkaline phosphatase/phosphodiesterase from the structural protein family PhoD that encoded by the genome sequence of the marine bacterium Cobetia amphilecti KMM 296 (CamPhoD) has been expressed in Escherichia coli cells. The calculated molecular weight, the number of amino acids, and the isoelectric point (pI) of the mature protein’s subunit are equal to 54832.98 Da, 492, and 5.08, respectively. The salt-tolerant, bimetal-dependent enzyme CamPhoD has a molecular weight of approximately 110 kDa in its native state. CamPhoD is activated by Co²⁺, Mg²⁺, Ca²⁺, or Fe³⁺ at a concentration of 2 mM and exhibits maximum activity in the presence of both Co²⁺ and Fe³⁺ ions in the incubation medium at pH 9.2. The exogenous ions, such as Zn²⁺, Cu²⁺, and Mn²⁺, as well as chelating agents EDTA and EGTA, do not have an appreciable effect on the CamPhoD activity. The temperature optimum for the CamPhoD activity is 45 °C. The enzyme catalyzes the cleavage of phosphate mono- and diester bonds in nucleotides, releasing inorganic phosphorus from p-nitrophenyl phosphate (pNPP) and guanosine 5′-triphosphate (GTP), as determined by the Chen method, with rate approximately 150- and 250-fold higher than those of bis-pNPP and 5′-pNP-TMP, respectively. The Michaelis–Menten constant (K_m), V_max, and efficiency (k_cat/K_m) of CamPhoD were 4.2 mM, 0.203 mM/min, and 7988.6 S⁻¹/mM; and 6.71 mM, 0.023 mM/min, and 1133.0 S⁻¹/mM for pNPP and bis-pNPP as the chromogenic substrates, respectively. Among the 3D structures currently available, in this study we found only the low identical structure of the Bacillus subtilis enzyme as a homologous template for modeling CamPhoD, with a new architecture of the phosphatase active site containing Fe³⁺ and two Ca²⁺ ions. It is evident that the marine bacterial phosphatase/phosphidiesterase CamPhoD is a new structural member of the PhoD family.

Serum alkaline phosphatase level is positively associated with metabolic syndrome: A nationwide population-based study

BACKGROUND

Serum alkaline phosphatase (ALP), a useful marker of hepatobiliary or bone disorder, has recently emerged as a biomarker of chronic low-grade inflammation and cardiometabolic disease. This study aimed to examine the association of serum ALP level with metabolic syndrome (MetS) in apparently healthy adults.

METHODS

A cross-sectional study was conducted to examine the relationship between serum ALP level and MetS in 7,101 men and 8,873 women aged 19 to 75 years using data from the 2008 to 2011 Korean National Health and Nutrition Examination Survey. The odds ratios (ORs) and 95% confidence intervals (CIs) for MetS were calculated using multiple logistic regression analyses across serum ALP quartiles (Q1: ≤190 U/L; Q2: 191-224 U/L; Q3: 225-263 U/L; and Q4: ≥264 U/L for men and Q1: ≤163 U/L; Q2: 164-201 U/L; Q3: 202-251 U/L; and Q4: ≥252 U/L for women).

RESULTS

The mean values of most cardiometabolic variables, HOMA-IR, and leukocyte count gradually increased with serum ALP quartile. The prevalence of MetS significantly increased in accordance with serum ALP quartile. In comparison with those of individuals in the lowest quartile, the OR (95% CI) for MetS in the highest quartile was 1.32 (1.05-1.64) in men and 1.99 (1.42-3.81) in women after adjusting for age, cigarette smoking, alcohol intake, regular exercise, household income, education level, occupation, AST, ALT, and GGT levels.

CONCLUSION

Serum ALP level was positively and independently associated with MetS in men and women.

Mechanism of catalysis and inhibition of Mycobacterium tuberculosis SapM, implications for the development of novel antivirulence drugs

Mycobacterium tuberculosis (Mtb) SapM is a secreted virulence factor critical for intracellular survival of the pathogen. The role of SapM in phagosome maturation arrest in host macrophages suggests its potential as a drug target to assist in the clearance of tuberculosis infection. However, the mechanism of action of SapM at the molecular level remains unknown. In this study, we provide new insights into the mechanism of catalysis, substrate specificity and inhibition of SapM, and we identify the critical residues for catalysis and substrate binding. Our findings demonstrate that SapM is an atypical monoester alkaline phosphatase, with a serine-based mechanism of catalysis probably metal-dependent. Particularly relevant to SapM function and pathogenesis, is its activity towards PI(4,5)P₂ and PI3P, two phosphoinositides that function at the early stages of microbial phagocytosis and phagosome formation. This suggests that SapM may have a pleiotropic role with a wider importance on Mtb infection than initially thought. Finally, we have identified two inhibitors of SapM, L-ascorbic acid and 2-phospho-L-ascorbic, which define two different mechanisms by which the catalytic activity of this phosphatase could be regulated. Critically, we demonstrate that 2-phospho-L-ascorbic reduces mycobacterial survival in macrophage infections, hence confirming the potential of SapM as a therapeutic drug target.

Catalytic Hydrolysis of Phosphate Monoester by Supramolecular Complexes Formed by the Self-Assembly of a Hydrophobic Bis(Zn2+-cyclen) Complex, Copper, and Barbital Units That Are Functionalized with Amino Acids in a Two-Phase Solvent System

We previously reported on the preparation of supramolecular complexes by the 2:2:2 assembly of a dinuclear Zn²⁺-cyclen (cyclen = 1,4,7,10-tetraazacyclododecane) complex having a 2,2′-bipyridyl linker equipped with 0~2 long alkyl chains (Zn₂L¹~Zn₂L³), 5,5-diethylbarbituric acid (Bar) derivatives, and a copper(II) ion (Cu²⁺) in aqueous solution and two-phase solvent systems and their phosphatase activities for the hydrolysis of mono(4-nitrophenyl) phosphate (MNP). These supermolecules contain Cu₂(μ-OH)₂ core that mimics the active site of alkaline phosphatase (AP), and one of the ethyl groups of the barbital moiety is located in close proximity to the Cu₂(μ-OH)₂ core. The generally accepted knowledge that the amino acids around the metal center in the active site of AP play important roles in its hydrolytic activity inspired us to modify the side chain of Bar with various functional groups in an attempt to mimic the active site of AP in the artificial system, especially in two-phase solvent system. In this paper, we report on the design and synthesis of new supramolecular complexes that are prepared by the combined use of bis(Zn²⁺-cyclen) complexes (Zn₂L¹, Zn₂L², and Zn₂L³), Cu²⁺, and Bar derivatives containing amino acid residues. We present successful formation of these artificial AP mimics with respect to the kinetics of the MNP hydrolysis obeying Michaelis–Menten scheme in aqueous solution and a two-phase solvent system and to the mode of the product inhibition by inorganic phosphate.

Catalytic Hydrolysis of Phosphate Monoester by Supramolecular Phosphatases Formed from a Monoalkylated Dizinc(II) Complex, Cyclic Diimide Units, and Copper(II) in Two-Phase Solvent System

Design and synthesis of enzyme mimic with programmed molecular interaction among several building blocks including metal complexes and metal chelators is of intellectual and practical significance. The preparation of artificial enzymes that mimic the natural enzymes such as hydrolases, phosphatases, etc. remains a great challenge in the field of supramolecular chemistry. Herein we report on the design and synthesis of asymmetric (nonsymmetric) supermolecules by the 2:2:2 self-assembly of an amphiphilic zinc(II)-cyclen complex containing a 2,2'-bipyridyl linker and one long alkyl chain (Zn₂L³), barbital analogues, and Cu²⁺ as model compounds of an enzyme alkaline phosphatase that catalyzes the hydrolysis of phosphate monoesters such as mono(4-nitrophenyl)phosphate at neutral pH in two-phase solvent system (H₂O/CHCl₃) in pH 7.4 and 37 °C. Hydrolytic activity of these complexes was found to be catalytic, and their catalytic turnover numbers are 3-4. The mechanistic studies based on the UV/vis and emission spectra of the H₂O and CHCl₃ phases of the reaction mixtures suggest that the hydrophilicity/hydrophobicity balance of the supramolecular catalysts is an important factor for catalytic activity.

Balancing Specificity and Promiscuity in Enzyme Evolution: Multidimensional Activity Transitions in the Alkaline Phosphatase Superfamily

Highly proficient, promiscuous enzymes can be springboards for functional evolution, able to avoid loss of function during adaptation by their capacity to promote multiple reactions. We employ a systematic comparative study of structure, sequence, and substrate specificity to track the evolution of specificity and reactivity between promiscuous members of clades of the alkaline phosphatase (AP) superfamily. Construction of a phylogenetic tree of protein sequences maps out the likely transition zone between arylsulfatases (ASs) and phosphonate monoester hydrolases (PMHs). Kinetic analysis shows that all enzymes characterized have four chemically distinct phospho- and sulfoesterase activities, with rate accelerations ranging from 10¹¹- to 10¹⁷-fold for their primary and 10⁹- to 10¹²-fold for their promiscuous reactions, suggesting that catalytic promiscuity is widespread in the AP-superfamily. This functional characterization and crystallography reveal a novel class of ASs that is so similar in sequence to known PMHs that it had not been recognized as having diverged in function. Based on analysis of snapshots of catalytic promiscuity "in transition", we develop possible models that would allow functional evolution and determine scenarios for trade-off between multiple activities. For the new ASs, we observe largely invariant substrate specificity that would facilitate the transition from ASs to PMHs via trade-off-free molecular exaptation, that is, evolution without initial loss of primary activity and specificity toward the original substrate. This ability to bypass low activity generalists provides a molecular solution to avoid adaptive conflict.

Analysis of Phosphoryl-Transfer Enzymes with QM/MM Free Energy Simulations

We discuss the application of quantum mechanics/molecular mechanics (QM/MM) free energy simulations to the analysis of phosphoryl transfers catalyzed by two enzymes: alkaline phosphatase and myosin. We focus on the nature of the transition state and the issue of mechanochemical coupling, respectively, in the two enzymes. The results illustrate unique insights that emerged from the QM/MM simulations, especially concerning the interpretation of experimental data regarding the nature of enzymatic transition states and coupling between global structural transition and catalysis in the active site. We also highlight a number of technical issues worthy of attention when applying QM/MM free energy simulations, and comment on a number of technical and mechanistic issues that require further studies.

Structural and Mechanistic Analysis of the Choline Sulfatase from Sinorhizobium melliloti: A Class I Sulfatase Specific for an Alkyl Sulfate Ester

Hydrolysis of organic sulfate esters proceeds by two distinct mechanisms, water attacking at either sulfur (S-O bond cleavage) or carbon (C-O bond cleavage). In primary and secondary alkyl sulfates, attack at carbon is favored, whereas in aromatic sulfates and sulfated sugars, attack at sulfur is preferred. This mechanistic distinction is mirrored in the classification of enzymes that catalyze sulfate ester hydrolysis: arylsulfatases (ASs) catalyze S-O cleavage in sulfate sugars and arylsulfates, and alkyl sulfatases break the C-O bond of alkyl sulfates. Sinorhizobium meliloti choline sulfatase (SmCS) efficiently catalyzes the hydrolysis of alkyl sulfate choline-O-sulfate (kcat/KM=4.8×103s-1M-1) as well as arylsulfate 4-nitrophenyl sulfate (kcat/KM=12s-1M-1). Its 2.8-Å resolution X-ray structure shows a buried, largely hydrophobic active site in which a conserved glutamate (Glu386) plays a role in recognition of the quaternary ammonium group of the choline substrate. SmCS structurally resembles members of the alkaline phosphatase superfamily, being most closely related to dimeric ASs and tetrameric phosphonate monoester hydrolases. Although >70% of the amino acids between protomers align structurally (RMSDs 1.79-1.99Å), the oligomeric structures show distinctly different packing and protomer-protomer interfaces. The latter also play an important role in active site formation. Mutagenesis of the conserved active site residues typical for ASs, H218O-labeling studies and the observation of catalytically promiscuous behavior toward phosphoesters confirm the close relation to alkaline phosphatase superfamily members and suggest that SmCS is an AS that catalyzes S-O cleavage in alkyl sulfate esters with extreme catalytic proficiency.

Challenges and advances in the computational modeling of biological phosphate hydrolysis

Phosphate ester hydrolysis is fundamental to many life processes, and has been the topic of substantial experimental and computational research effort. However, even the simplest of phosphate esters can be hydrolyzed through multiple possible pathways that can be difficult to distinguish between, either experimentally, or computationally. Therefore, the mechanisms of both the enzymatic and non-enzymatic reactions have been historically controversial. In the present contribution, we highlight a number of technical issues involved in reliably modeling these computationally challenging reactions, as well as proposing potential solutions. We also showcase examples of our own work in this area, discussing both the non-enzymatic reaction in aqueous solution, as well as insights obtained from the computational modeling of organophosphate hydrolysis and catalytic promiscuity amongst enzymes that catalyze phosphoryl transfer.

An atypical phosphodiesterase capable of degrading haloalkyl phosphate diesters from Sphingobium sp. strain TCM1

Sphingobium sp. strain TCM1 can degrade tris(2-chloroethyl) phosphate (TCEP) to inorganic phosphate and 2-chloroethanol. A phosphotriesterase (PTE), phosphodiesterase (PDE) and phosphomonoesterase (PME) are believed to be involved in the degradation of TCEP. The PTE and PME that respectively catalyze the first and third steps of TCEP degradation in TCM1 have been identified. However, no information has been reported on a PDE catalyzing the second step. In this study, we identified, purified, and characterized a PDE capable of hydrolyzing haloalkyl phosphate diesters. The final preparation of the enzyme had a specific activity of 29 µmol min⁻¹ mg⁻¹ with bis(p-nitrophenyl) phosphate (BpNPP) as the substrate. It also possessed low PME activity with p-nitrophenyl phosphate (pNPP) as substrate. The catalytic efficiency (k_cat/K_m) with BpNPP was significantly higher than that with pNPP, indicating that the enzyme prefers the organophosphorus diester to the monoester. The enzyme degraded bis(2,3-dibromopropyl) phosphate, bis(1,3-dichloro-2-propyl) phosphate and bis(2-chloroethyl) phosphate, suggesting that it is involved in the metabolism of haloalkyl organophosphorus triesters. The primary structure of the PDE from TCM1 is distinct from those of typical PDE family members and the enzyme belongs to the polymerase and histidinol phosphatase superfamily.

Lysophosphatidylcholine elicits intracellular calcium signaling in a GPR55-dependent manner

The GPR55 signaling is fertile ground for drug discovery, however despite considerable research progress during the past 10 years, many open questions remain. The GPR55 pharmacology remains controversial, as many ligands have been reported with inconsistent results. Here, we show that various molecular species of lysophosphatidylcholine (LPC) elicit intracellular Ca²⁺ mobilization in GPR55-expressing PC-3 human prostate carcinoma cells. The response was even stronger than [Ca²⁺]_i flux evoked by endogenous (OEA) and synthetic (Abn-CBD) agonists. Treatment with GPR55 antagonists CID16020046 and ML193 as well as the lipid raft disrupter methyl-β-cyclodextrin strongly blunted LPC-induced calcium signal. Additionally, molecular modeling analysis revealed that LPC 16:0 and LPC 18:1 interact stronger with the receptor than to OEA. Identified electrostatic interactions between GPR55 residues and the ligands overlap with the binding site identified previously for lysophosphatidylinositol. Therefore, we prove that LPC is another GPR55-sensitive ligand. This finding is relevant in understanding lysophospolipids-mediated signaling and opens new avenues to develop therapeutic approach based on GPR55 targeting.

Genome-wide TOP2A DNA cleavage is biased toward translocated and highly transcribed loci

Type II topoisomerases orchestrate proper DNA topology, and they are the targets of anti-cancer drugs that cause treatment-related leukemias with balanced translocations. Here, we develop a high-throughput sequencing technology to define TOP2 cleavage sites at single-base precision, and use the technology to characterize TOP2A cleavage genome-wide in the human K562 leukemia cell line. We find that TOP2A cleavage has functionally conserved local sequence preferences, occurs in cleavage cluster regions (CCRs), and is enriched in introns and lincRNA loci. TOP2A CCRs are biased toward the distal regions of gene bodies, and TOP2 poisons cause a proximal shift in their distribution. We find high TOP2A cleavage levels in genes involved in translocations in TOP2 poison–related leukemia. In addition, we find that a large proportion of genes involved in oncogenic translocations overall contain TOP2A CCRs. The TOP2A cleavage of coding and lincRNA genes is independently associated with both length and transcript abundance. Comparisons to ENCODE data reveal distinct TOP2A CCR clusters that overlap with marks of transcription, open chromatin, and enhancers. Our findings implicate TOP2A cleavage as a broad DNA damage mechanism in oncogenic translocations as well as a functional role of TOP2A cleavage in regulating transcription elongation and gene activation.

Mechanistic and Evolutionary Insights from Comparative Enzymology of Phosphomonoesterases and Phosphodiesterases across the Alkaline Phosphatase Superfamily

Naively one might have expected an early division between phosphate monoesterases and diesterases of the alkaline phosphatase (AP) superfamily. On the contrary, prior results and our structural and biochemical analyses of phosphate monoesterase PafA, from Chryseobacterium meningosepticum, indicate similarities to a superfamily phosphate diesterase [Xanthomonas citri nucleotide pyrophosphatase/phosphodiesterase (NPP)] and distinct differences from the three metal ion AP superfamily monoesterase, from Escherichia coli AP (EcAP). We carried out a series of experiments to map out and learn from the differences and similarities between these enzymes. First, we asked why there would be independent instances of monoesterases in the AP superfamily? PafA has a much weaker product inhibition and slightly higher activity relative to EcAP, suggesting that different metabolic evolutionary pressures favored distinct active-site architectures. Next, we addressed the preferential phosphate monoester and diester catalysis of PafA and NPP, respectively. We asked whether the >80% sequence differences throughout these scaffolds provide functional specialization for each enzyme's cognate reaction. In contrast to expectations from this model, PafA and NPP mutants with the common subset of active-site groups embedded in each native scaffold had the same monoesterase:diesterase specificities; thus, the >10⁷-fold difference in native specificities appears to arise from distinct interactions at a single phosphoryl substituent. We also uncovered striking mechanistic similarities between the PafA and EcAP monoesterases, including evidence for ground-state destabilization and functional active-site networks that involve different active-site groups but may play analogous catalytic roles. Discovering common network functions may reveal active-site architectural connections that are critical for function, and identifying regions of functional modularity may facilitate the design of new enzymes from existing promiscuous templates. More generally, comparative enzymology and analysis of catalytic promiscuity can provide mechanistic and evolutionary insights.

QM/MM Analysis of Transition States and Transition State Analogues in Metalloenzymes

Enzymology is approaching an era where many problems can benefit from computational studies. While ample challenges remain in quantitatively predicting behavior for many enzyme systems, the insights that often come from computations are an important asset for the enzymology community. Here we provide a primer for enzymologists on the types of calculations that are most useful for mechanistic problems in enzymology. In particular, we emphasize the integration of models that range from small active-site motifs to fully solvated enzyme systems for cross-validation and dissection of specific contributions from the enzyme environment. We then use a case study of the enzyme alkaline phosphatase to illustrate specific application of the methods. The case study involves examination of the binding modes of putative transition state analogues (tungstate and vanadate) to the enzyme. The computations predict covalent binding of these ions to the enzymatic nucleophile and that they adopt the trigonal bipyramidal geometry of the expected transition state. By comparing these structures with transition states found through free energy simulations, we assess the degree to which the transition state analogues mimic the true transition states. Technical issues worth treating with care as well as several remaining challenges to quantitative analysis of metalloenzymes are also highlighted during the discussion.

Leaving Group Ability Observably Affects Transition State Structure in a Single Enzyme Active Site

A reaction's transition state (TS) structure plays a critical role in determining reactivity and has important implications for the design of catalysts, drugs, and other applications. Here, we explore TS structure in the enzyme alkaline phosphatase using hybrid Quantum Mechanics/Molecular Mechanics simulations. We find that minor perturbations to the substrate have major effects on TS structure and the way the enzyme stabilizes the TS. Substrates with good leaving groups (LGs) have little cleavage of the phosphorus-LG bond at the TS, while substrates with poor LGs have substantial cleavage of that bond. The results predict nonlinear free energy relationships for a single rate-determining step, and substantial differences in kinetic isotope effects for different substrates; both trends were observed in previous experimental studies, although the original interpretations differed from the present model. Moreover, due to different degrees of phosphorus-LG bond cleavage at the TS for different substrates, the LG is stabilized by different interactions at the TS: while a poor LG is directly stabilized by an active site zinc ion, a good LG is mainly stabilized by active site water molecules. Our results demonstrate the considerable plasticity of TS structure and stabilization in enzymes. Furthermore, perturbations to reactivity that probe TS structure experimentally (i.e., substituent effects) may substantially perturb the TS they aim to probe, and thus classical experimental approaches such as free energy relations should be interpreted with care.

Promiscuity in the Enzymatic Catalysis of Phosphate and Sulfate Transfer

The enzymes that facilitate phosphate and sulfate hydrolysis are among the most proficient natural catalysts known to date. Interestingly, a large number of these enzymes are promiscuous catalysts that exhibit both phosphatase and sulfatase activities in the same active site and, on top of that, have also been demonstrated to efficiently catalyze the hydrolysis of other additional substrates with varying degrees of efficiency. Understanding the factors that underlie such multifunctionality is crucial both for understanding functional evolution in enzyme superfamilies and for the development of artificial enzymes. In this Current Topic, we have primarily focused on the structural and mechanistic basis for catalytic promiscuity among enzymes that facilitate both phosphoryl and sulfuryl transfer in the same active site, while comparing this to how catalytic promiscuity manifests in other promiscuous phosphatases. We have also drawn on the large number of experimental and computational studies of selected model systems in the literature to explore the different features driving the catalytic promiscuity of such enzymes. Finally, on the basis of this comparative analysis, we probe the plausible origins and determinants of catalytic promiscuity in enzymes that catalyze phosphoryl and sulfuryl transfer.

Catalytic and substrate promiscuity: distinct multiple chemistries catalysed by the phosphatase domain of receptor protein tyrosine phosphatase

The presence of latent activities in enzymes is posited to underlie the natural evolution of new catalytic functions. However, the prevalence and extent of such substrate and catalytic ambiguity in evolved enzymes is difficult to address experimentally given the order-of-magnitude difference in the activities for native and, sometimes, promiscuous substrate/s. Further, such latent functions are of special interest when the activities concerned do not fall into the domain of substrate promiscuity. In the present study, we show a special case of such latent enzyme activity by demonstrating the presence of two mechanistically distinct reactions catalysed by the catalytic domain of receptor protein tyrosine phosphatase isoform δ (PTPRδ). The primary catalytic activity involves the hydrolysis of a phosphomonoester bond (C─O─P) with high catalytic efficiency, whereas the secondary activity is the hydrolysis of a glycosidic bond (C─O─C) with poorer catalytic efficiency. This enzyme also displays substrate promiscuity by hydrolysing diester bonds while being highly discriminative for its monoester substrates. To confirm these activities, we also demonstrated their presence on the catalytic domain of protein tyrosine phosphatase Ω (PTPRΩ), a homologue of PTPRδ. Studies on the rate, metal-ion dependence, pH dependence and inhibition of the respective activities showed that they are markedly different. This is the first study that demonstrates a novel sugar hydrolase and diesterase activity for the phosphatase domain (PD) of PTPRδ and PTPRΩ. This work has significant implications for both understanding the evolution of enzymatic activity and the possible physiological role of this new chemistry. Our findings suggest that the genome might harbour a wealth of such alternative latent enzyme activities in the same protein domain that renders our knowledge of metabolic networks incomplete.

Tungstate as a Transition State Analog for Catalysis by Alkaline Phosphatase

The catalytic mechanisms underlying Escherichia coli alkaline phosphatase's (AP) remarkable rate enhancement have been probed extensively. Past work indicated that whereas the serine nucleophile (Ser102) electrostatically repels the product phosphate, another oxyanion, tungstate, binds more strongly in the presence of Ser102. These results predict a covalent bond between the serine nucleophile and tungstate, a model that we test herein. The crystal structure of tungstate-bound alkaline phosphatase provides evidence for a covalent adduct model and further shows that the ligand adopts trigonal bipyramidal geometry, which is infrequently observed for tungstate in small molecules and other active sites but mirrors the geometry of the presumed phosphoryl transfer transition state. The AP active site is known to stabilize another oxyanion, vanadate, in trigonal bipyramidal geometry, but the extent to which binding of either ligand reproduces the energetics of the transition state cannot be deduced from structural inspection alone. To test for transition state analog behavior, we determined the relationship between catalytic activity and affinity for tungstate and vanadate for a series of 20 AP variants. Affinity and activity were highly correlated for tungstate (r(2) = 0.89) but not vanadate (r(2) = 0.23), indicating that the tungstate•AP complex may better mimic this enzyme's transition state properties. The results herein suggest that tungstate will be a valuable tool for further dissecting AP catalysis and may prove helpful in mechanistic studies of other phosphoryl transfer enzymes.

High-throughput analysis and protein engineering using microcapillary arrays

We describe a multipurpose technology platform, termed μSCALE (microcapillary single-cell analysis and laser extraction), that enables massively parallel, quantitative biochemical and biophysical measurements on millions of protein variants expressed from yeast or bacteria. μSCALE spatially segregates single cells within a microcapillary array, enabling repeated imaging, cell growth and protein expression. We performed high-throughput analysis of cells and their protein products using a range of fluorescent assays, including binding-affinity measurements and dynamic enzymatic assays. A precise laser-based extraction method allows rapid recovery of live clones and their genetic material from microcapillaries for further study. With μSCALE, we discovered a new antibody against a clinical cancer target, evolved a fluorescent protein biosensor and engineered an enzyme to reduce its sensitivity to its inhibitor. These protein analysis and engineering applications each have unique assay requirements and different host organisms, highlighting the flexibility and technical capabilities of the μSCALE platform.

Promiscuity and electrostatic flexibility in the alkaline phosphatase superfamily

Catalytic promiscuity, that is, the ability of single enzymes to facilitate the turnover of multiple, chemically distinct substrates, is a widespread phenomenon that plays an important role in the evolution of enzyme function. Additionally, such pre-existing multifunctionality can be harnessed in artificial enzyme design. The members of the alkaline phosphatase superfamily have served extensively as both experimental and computational model systems for enhancing our understanding of catalytic promiscuity. In this Opinion, we present key recent computational studies into the catalytic activity of these highly promiscuous enzymes, highlighting the valuable insight they have provided into both the molecular basis for catalytic promiscuity in general, and its implications for the evolution of phosphatase activity.

Cooperative Electrostatic Interactions Drive Functional Evolution in the Alkaline Phosphatase Superfamily

It is becoming widely accepted that catalytic promiscuity, i.e., the ability of a single enzyme to catalyze the turnover of multiple, chemically distinct substrates, plays a key role in the evolution of new enzyme functions. In this context, the members of the alkaline phosphatase superfamily have been extensively studied as model systems in order to understand the phenomenon of enzyme multifunctionality. In the present work, we model the selectivity of two multiply promiscuous members of this superfamily, namely the phosphonate monoester hydrolases from Burkholderia caryophylli and Rhizobium leguminosarum. We have performed extensive simulations of the enzymatic reaction of both wild-type enzymes and several experimentally characterized mutants. Our computational models are in agreement with key experimental observables, such as the observed activities of the wild-type enzymes, qualitative interpretations of experimental pH-rate profiles, and activity trends among several active site mutants. In all cases the substrates of interest bind to the enzyme in similar conformations, with largely unperturbed transition states from their corresponding analogues in aqueous solution. Examination of transition-state geometries and the contribution of individual residues to the calculated activation barriers suggest that the broad promiscuity of these enzymes arises from cooperative electrostatic interactions in the active site, allowing each enzyme to adapt to the electrostatic needs of different substrates. By comparing the structural and electrostatic features of several alkaline phosphatases, we suggest that this phenomenon is a generalized feature driving selectivity and promiscuity within this superfamily and can be in turn used for artificial enzyme design.

Extensive site-directed mutagenesis reveals interconnected functional units in the alkaline phosphatase active site

Enzymes enable life by accelerating reaction rates to biological timescales. Conventional studies have focused on identifying the residues that have a direct involvement in an enzymatic reaction, but these so-called 'catalytic residues' are embedded in extensive interaction networks. Although fundamental to our understanding of enzyme function, evolution, and engineering, the properties of these networks have yet to be quantitatively and systematically explored. We dissected an interaction network of five residues in the active site of Escherichia coli alkaline phosphatase. Analysis of the complex catalytic interdependence of specific residues identified three energetically independent but structurally interconnected functional units with distinct modes of cooperativity. From an evolutionary perspective, this network is orders of magnitude more probable to arise than a fully cooperative network. From a functional perspective, new catalytic insights emerge. Further, such comprehensive energetic characterization will be necessary to benchmark the algorithms required to rationally engineer highly efficient enzymes.

Identification and characterization of a novel phosphodiesterase from the metagenome of an Indian coalbed

Phosphoesterases are involved in the degradation of organophosphorus compounds. Although phosphomonoesterases and phosphotriesterases have been studied in detail, studies on phosphodiesterases are rather limited. In our search to find novel phosphodiesterases using metagenomic approach, we cloned a gene encoding a putative phosphodiesterase (PdeM) from the metagenome of the formation water collected from an Indian coal bed. Bioinformatic analysis showed that PdeM sequence possessed the characteristic signature motifs of the class III phosphodiesterases and phylogenetic study of PdeM enabled us to identify three distinct subclasses (A, B, and C) within class III phosphodiesterases, PdeM clustering in new subclass IIIB. Bioinformatic, biochemical and biophysical characterization of PdeM further revealed some of the characteristic features of the phosphodiesterases belonging to newly described subclass IIIB. PdeM is a monomer of 29.3 kDa, which exhibits optimum activity at 25°C and pH 8.5, but low affinity for bis(pNPP) as well as pNPPP. The recombinant PdeM possessed phosphodiesterase, phosphonate-ester hydrolase and nuclease activity. It lacked phosphomonoesterase, phosphotriesterase, and RNAse activities. Overexpression of PdeM in E.coli neither affected catabolite respression nor did the recombinant protein hydrolyzed cAMP in vitro, indicating its inability to hydrolyze cAMP. Although Mn²⁺ was required for the activity of PdeM, but addition of metals (Mn²⁺ or Fe³⁺) did not induce oligomerization. Further increase in concentration of Mn²⁺ upto 3 mM, increased α-helical content as well as the phosphodiesterase activity. Structural comparison of PdeM with its homologs showed that it lacked critical residues required for dimerization, cAMP hydrolysis, and for the high affinity binding of bis(pNPP). PdeM, thus, is a novel representative of new subclass of class III phosphodiesterases.

Probing the Origins of Catalytic Discrimination between Phosphate and Sulfate Monoester Hydrolysis: Comparative Analysis of Alkaline Phosphatase and Protein Tyrosine Phosphatases

Catalytic promiscuity, the ability of enzymes to catalyze multiple reactions, provides an opportunity to gain a deeper understanding of the origins of catalysis and substrate specificity. Alkaline phosphatase (AP) catalyzes both phosphate and sulfate monoester hydrolysis reactions with a ∼10¹⁰-fold preference for phosphate monoester hydrolysis, despite the similarity between these reactions. The preponderance of formal positive charge in the AP active site, particularly from three divalent metal ions, was proposed to be responsible for this preference by providing stronger electrostatic interactions with the more negatively charged phosphoryl group versus the sulfuryl group. To test whether positively charged metal ions are required to achieve a high preference for the phosphate monoester hydrolysis reaction, the catalytic preference of three protein tyrosine phosphatases (PTPs), which do not contain metal ions, were measured. Their preferences ranged from 5 × 10⁶ to 7 × 10⁷, lower than that for AP but still substantial, indicating that metal ions and a high preponderance of formal positive charge within the active site are not required to achieve a strong catalytic preference for phosphate monoester over sulfate monoester hydrolysis. The observed ionic strength dependences of k_cat/K_M values for phosphate and sulfate monoester hydrolysis are steeper for the more highly charged phosphate ester with both AP and the PTP Stp1, following the dependence expected based on the charge difference of these two substrates. However, the dependences for AP were not greater than those of Stp1 and were rather shallow for both enzymes. These results suggest that overall electrostatics from formal positive charge within the active site is not the major driving force in distinguishing between these reactions and that substantial discrimination can be attained without metal ions. Thus, local properties of the active site, presumably including multiple positioned dipolar hydrogen bond donors within the active site, dominate in defining this reaction specificity.

Alkaline phosphatases in microbialites and bacterioplankton from Alchichica soda lake, Mexico

Dissolved organic phosphorus utilization by different members of natural communities has been closely linked to microbial alkaline phosphatases whose affiliation and diversity is largely unknown. Here we assessed genetic diversity of bacterial alkaline phosphatases phoX and phoD, using highly diverse microbial consortia (microbialites and bacterioplankton) as study models. These microbial consortia are found in an oligo-mesotrophic soda lake with a particular geochemistry, exhibiting a low calcium concentration and a high Mg : Ca ratio relative to seawater. In spite of the relative low calcium concentration in the studied system, our results highlight the diversity of calcium-based metallophosphatases phoX and phoD-like in heterotrophic bacteria of microbialites and bacterioplankton, where phoX was the most abundant alkaline phosphatase found. phoX and phoD-like phylotypes were more numerous in microbialites than in bacterioplankton. A larger potential community for DOP utilization in microbialites was consistent with the TN : TP ratio, suggesting P limitation within these assemblages. A cross-system comparison indicated that diversity of phoX in Lake Alchichica was similar to that of other aquatic systems with a naturally contrasting ionic composition and trophic state, although no phylotypes were shared among systems.

Crystallographic study of the phosphoethanolamine transferase EptC required for polymyxin resistance and motility in Campylobacter jejuni

The foodborne enteric pathogen Campylobacter jejuni decorates a variety of its cell-surface structures with phosphoethanolamine (pEtN). Modifying lipid A with pEtN promotes cationic antimicrobial peptide resistance, whereas post-translationally modifying the flagellar rod protein FlgG with pEtN promotes flagellar assembly and motility, which are processes that are important for intestinal colonization. EptC, the pEtN transferase required for all known pEtN cell-surface modifications in C. jejuni, is a predicted inner-membrane metalloenzyme with a five-helix N-terminal transmembrane domain followed by a soluble sulfatase-like catalytic domain in the periplasm. The atomic structure of the catalytic domain of EptC (cEptC) was crystallized and solved to a resolution of 2.40 Å. cEptC adopts the α/β/α fold of the sulfatase protein family and harbors a zinc-binding site. A phosphorylated Thr266 residue was observed that was hypothesized to mimic a covalent pEtN-enzyme intermediate. The requirement for Thr266 as well as the nearby residues Asn308, Ser309, His358 and His440 was ascertained via in vivo activity assays on mutant strains. The results establish a basis for the design of pEtN transferase inhibitors.

Differential effects of a mutation on the normal and promiscuous activities of orthologs: implications for natural and directed evolution

Neutral drift occurring over millions or billions of years results in substantial sequence divergence among enzymes that catalyze the same reaction. Although natural selection maintains the primary activity of orthologous enzymes, there is, by definition, no selective pressure to maintain physiologically irrelevant promiscuous activities. Thus, the levels and the evolvabilities of promiscuous activities may vary among orthologous enzymes. Consistent with this expectation, we have found that the levels of a promiscuous activity in nine gamma-glutamyl phosphate reductase (ProA) orthologs vary by about 50-fold. Remarkably, a single amino acid change from Glu to Ala near the active site appeared to be critical for improvement of the promiscuous activity in every ortholog. The effects of this change varied dramatically. The improvement in the promiscuous activity varied from 50- to 770-fold, and, importantly, was not correlated with the initial level of the promiscuous activity. The decrease in the original activity varied from 190- to 2,100-fold. These results suggest that evolution of a novel enzyme may be possible in some microbes, but not in others. Further, these results underscore the importance of using multiple orthologs as starting points for directed evolution of novel enzyme activities.