Alkaline Phosphatase Superfamily Research Articles

Unveiling the molecular mechanisms of the type IX secretion system's response regulator: Structural and functional insights.

The type IX secretion system (T9SS) is a nanomachinery utilized by bacterial pathogens to facilitate infection. The system is regulated by a signaling cascade serving as its activation switch. A pivotal member in this cascade, the response regulator protein PorX, represents a promising drug target to prevent the secretion of virulence factors. Here, we provide a comprehensive characterization of PorX both in vitro and in vivo. First, our structural studies revealed PorX harbors a unique enzymatic effector domain, which, surprisingly, shares structural similarities with the alkaline phosphatase superfamily, involved in nucleotide and lipid signaling pathways. Importantly, such pathways have not been associated with the T9SS until now. Enzymatic characterization of PorX's effector domain revealed a zinc-dependent phosphodiesterase activity, with active site dimensions suitable to accommodate a large substrate. Unlike typical response regulators that dimerize via their receiver domain upon phosphorylation, we found that zinc can also induce conformational changes and promote PorX's dimerization via an unexpected interface. These findings suggest that PorX can serve as a cellular zinc sensor, broadening our understanding of its regulatory mechanisms. Despite the strict conservation of PorX in T9SS-utilizing bacteria, we demonstrate that PorX is essential for virulence factors secretion in Porphyromonas gingivalis and affects metabolic enzymes secretion in the nonpathogenic Flavobacterium johnsoniae, but not for the secretion of gliding adhesins. Overall, this study advances our structural and functional understanding of PorX, highlighting its potential as a druggable target for intervention strategies aimed at disrupting the T9SS and mitigating virulence in pathogenic species.

Mechanistic and evolutionary insights into alkaline phosphatase superfamily through structure-function studies on Sphingomonas alkaline phosphatase

Alkaline phosphatases (APs), represented by E. coli AP (ECAP), employ an arginine residue to stabilize the phosphoryl group in the active site; whereas, AP from Sphingomonas (SPAP) shows a unique combination of substrate-binding residues; Thr89, Asn110, Lys171, and Arg173. Although such combination has been observed only in SPAP, these residues are present separately in different members of the AP superfamily. Here, we establish the presence of two distinct classes of APs; ECAP-type and SPAP-type. Bioinformatic analyses show that SPAP-type of APs are widely distributed in the bacterial kingdom. The role of active site residues in the catalytic mechanism has been delineated through a set of crystal structures reported here. These structures, representing different stages of the reaction pathway provide wealth of information for the catalytic mechanism. Despite critical differences in the substrate binding residues, SPAP follows a mechanism similar to that of ECAP-type of APs. Structure-based phylogenetic analysis suggests that SPAP and ECAP may have diverged very early during the evolution from a common ancestor. Moreover, it is proposed that the SPAP-type of APs are fundamental members of the AP superfamily and are more closely related to other members of the superfamily as compared to the ECAP-type of APs.

Response regulator PorX coordinates oligonucleotide signalling and gene expression to control the secretion of virulence factors.

The PglZ family of proteins belongs to the alkaline phosphatase superfamily, which consists of metallohydrolases with limited sequence identity but similar metal-coordination architectures in otherwise divergent active sites. Proteins with a well-defined PglZ domain are ubiquitous among prokaryotes as essential components of BREX phage defence systems and two-component systems (TCSs). Whereas other members of the alkaline phosphatase superfamily are well characterized, the activity, structure and biological function of PglZ family proteins remain unclear. We therefore investigated the structure and function of PorX, an orphan response regulator of the Porphyromonas gingivalis TCS containing a putative PglZ effector domain. The crystal structure of PorX revealed a canonical receiver domain, a helical bundle, and an unprecedented PglZ domain, similar to the general organization of the phylogenetically related BREX-PglZ proteins. The PglZ domain of PorX features an active site cleft suitable for large substrates. An extensive search for substrates revealed that PorX is a phosphodiesterase that acts on cyclic and linear oligonucleotides, including signalling molecules such as cyclic oligoadenylates. These results, combined with mutagenesis, biophysical and enzymatic analysis, suggest that PorX coordinates oligonucleotide signalling pathways and indirectly regulates gene expression to control the secretion of virulence factors.

Structural insights into choline-O-sulfatase reveal the molecular determinants for ligand binding.

Choline-O-sulfatase (COSe; EC 3.1.6.6) is a member of the alkaline phosphatase (AP) superfamily, and its natural function is to hydrolyze choline-O-sulfate into choline and sulfate. Despite its natural function, the major interest in this enzyme resides in the landmark catalytic/substrate promiscuity of sulfatases, which has led to attention in the biotechnological field due to their potential in protein engineering. In this work, an in-depth structural analysis of wild-type Sinorhizobium (Ensifer) meliloti COSe (SmeCOSe) and its C54S active-site mutant is reported. The binding mode of this AP superfamily member to both products of the reaction (sulfate and choline) and to a substrate-like compound are shown for the first time. The structures further confirm the importance of the C-terminal extension of the enzyme in becoming part of the active site and participating in enzyme activity through dynamic intra-subunit and inter-subunit hydrogen bonds (Asn146A-Asp500B-Asn498B). These residues act as the `gatekeeper' responsible for the open/closed conformations of the enzyme, in addition to assisting in ligand binding through the rearrangement of Leu499 (with a movement of approximately 5 Å). Trp129 and His145 clamp the quaternary ammonium moiety of choline and also connect the catalytic cleft to the C-terminus of an adjacent protomer. The structural information reported here contrasts with the proposed role of conformational dynamics in promoting the enzymatic catalytic proficiency of an enzyme.

Alkaline Phosphatases: in Silico Study on the Catalytic Effect of Conserved Active Site Residues Using Human Placental Alkaline Phosphatase (PLAP) As a Model Protein.

The metalloenzymes from the alkaline phosphatase (AP) superfamily catalyze the hydrolysis and transphosphorylation of phosphate monoesters. The role of several amino acids highly conserved in the active site of this family of enzymes was examined, using human placental AP (PLAP) as a model protein. By employing an active-site model based on the X-ray crystal structure of PLAP, mutations of several key residues were modeled by quantum mechanical methods in order to determine their impact on the catalytic activity. Kinetic and thermodynamic estimations were achieved for each reaction step of the catalytic mechanism by characterization of the intermediates and transition states on the reaction pathway, and the effects of mutations on the activation barriers were analyzed. A good accordance was observed between the present computational results and experimental measurements reported in the literature.

The Escherichia coli cellulose synthase subunit G (BcsG) is a Zn2+-dependent phosphoethanolamine transferase

Bacterial biofilms are cellular communities that produce an adherent matrix. Exopolysaccharides are key structural components of this matrix and are required for the assembly and architecture of biofilms produced by a wide variety of microorganisms. The human bacterial pathogens Escherichia coli and Salmonella enterica produce a biofilm matrix composed primarily of the exopolysaccharide phosphoethanolamine (pEtN) cellulose. Once thought to be composed of only underivatized cellulose, the pEtN modification present in these matrices has been implicated in the overall architecture and integrity of the biofilm. However, an understanding of the mechanism underlying pEtN derivatization of the cellulose exopolysaccharide remains elusive. The bacterial cellulose synthase subunit G (BcsG) is a predicted inner membrane-localized metalloenzyme that has been proposed to catalyze the transfer of the pEtN group from membrane phospholipids to cellulose. Here we present evidence that the C-terminal domain of BcsG from E. coli (EcBcsGΔN) functions as a phosphoethanolamine transferase in vitro with substrate preference for cellulosic materials. Structural characterization of EcBcsGΔN revealed that it belongs to the alkaline phosphatase superfamily, contains a Zn2+ ion at its active center, and is structurally similar to characterized enzymes that confer colistin resistance in Gram-negative bacteria. Informed by our structural studies, we present a functional complementation experiment in E. coli AR3110, indicating that the activity of the BcsG C-terminal domain is essential for integrity of the pellicular biofilm. Furthermore, our results established a similar but distinct active-site architecture and catalytic mechanism shared between BcsG and the colistin resistance enzymes.

A High-Throughput Assay Platform for Next-Generation Mechanistic Enzymology and Applications

Twentieth-century biochemistry dissected cellular metabolism, identified enzymes responsible for its regulation, and revealed reaction pathways. Nevertheless, enhanced understanding is needed to dissect evolutionary pathways, rationally engineer novel catalysts, and predict disease effects of human alleles. Deep mechanistic studies of active site residues have revealed much about how they contribute to catalysis, but the enzyme scaffold is required for these residues to participate in catalysis. Deep mutational scanning and other current high-throughput approaches have revealed functional effects throughout enzyme scaffolds, but these studies intrinsically lack mechanistic depth. We have developed HT-MEK (High-Throughput Microfluidic Enzyme Kinetics), a technology for expression, purification, and assay of >1000 rationally-selected enzyme variants in a single experiment, thereby merging the power of traditional and recent high-throughput technologies. HT-MEK yields quantitative measurements of multiple kinetic and thermodynamic parameters for each variant. As a first demonstration of HT-MEK, we performed enzyme-wide mutational scans of a well-characterized Alkaline Phosphatase superfamily member, PafA. We made thousands of kinetic and thermodynamic measurements for over 1000 PafA variants with a battery of substrates and inhibitors to probe specific catalytic features. Consistent with the simplest expectation, the affecting mutations cluster around the active site regions previously shown to be directly involved in each mechanism. Nevertheless, the extent and direction of these clusters was not predictable, nor was the observation that several clusters extend to the protein surface, demonstrating that HT-MEK can be used to identify surfaces that may allow allosteric control of particular functions. Parallel approaches provide analogous quantitative stability measurements, some of which we will also present. The efficient and systematic dissection of any enzyme or system of enzymes allowed by HT-MEK enables next-generation applications to engineer allostery, novel catalysts, and efficient biosynthetic pathways for use in industry and medicine.

Resistance to the "last resort" antibiotic colistin: a single-zinc mechanism for phosphointermediate formation in MCR enzymes.

MCR (mobile colistin resistance) enzymes catalyse phosphoethanolamine (PEA) addition to bacterial lipid A, threatening the "last-resort" antibiotic colistin. Molecular dynamics and density functional theory simulations indicate that monozinc MCR supports PEA transfer to the Thr285 acceptor, positioning MCR as a mono- rather than multinuclear member of the alkaline phosphatase superfamily.

Bringing Enzymology into the Genomic Era: Developing and Deploying New Tools to Quantitatively Map Functional Connections Throughout an Enzyme

Despite enormous advances in our understanding of enzymes and their mechanisms over the past decades, how they achieve their enormous rate enhancements and exquisite specificities remains elusive, as evidenced by our limited ability to engineer new enzymes. A fundamental reason for this lack of understanding is the large number of interconnected residues and long-range interactions that cannot be revealed via traditional, low-throughput quantitative assays. While high-throughput screens and selections can identify enzymes with desired changes, the improvements obtained are typically modest and the qualitative or limited quantitative data obtainable by these approaches is insufficient to reveal functional interconnections between remote and active site residues. To meet this unmet need, we have developed a new technology, HT-MEK (High-throughput Microfluidic Enzyme Kinetics), that allows measurement of myriad quantitative biochemical parameters for 1000 rationally chosen enzyme variants in a single experiment. Using this platform, we carried out deep functional analysis of multiple substitutions at each position within a single enzyme, the Alkaline Phosphatase superfamily member PafA. For each substitution, we measured multiple biochemical parameters, including Michaelis-Menten constants for cognate and non-cognate substrates and inhibition constants for multiple inhibitors and transition state analogs. The results provide a detailed functional signature for each position within PafA. We observe distinct physically contiguous ‘functional units’ throughout the enzyme that tune specific aspects of function, yielding blueprints for future rational design efforts. HT-MEK is applicable to enzyme function and stability, provides data that rivals or exceeds the precision typically obtained from standard assays, and reduces the time, effort, and cost to quantitatively profile the effects of individual substitutions on enzyme function by orders of magnitude, providing a technology to access previously inaccessible questions across multiple systems and disciplines.

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 1011- to 1017-fold for their primary and 109- to 1012-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.

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.

Differential catalytic promiscuity of the alkaline phosphatase superfamily bimetallo core reveals mechanistic features underlying enzyme evolution

Members of enzyme superfamilies specialize in different reactions but often exhibit catalytic promiscuity for one another's reactions, consistent with catalytic promiscuity as an important driver in the evolution of new enzymes. Wanting to understand how catalytic promiscuity and other factors may influence evolution across a superfamily, we turned to the well-studied alkaline phosphatase (AP) superfamily, comparing three of its members, two evolutionarily distinct phosphatases and a phosphodiesterase. We mutated distinguishing active-site residues to generate enzymes that had a common Zn2+ bimetallo core but little sequence similarity and different auxiliary domains. We then tested the catalytic capabilities of these pruned enzymes with a series of substrates. A substantial rate enhancement of ∼1011-fold for both phosphate mono- and diester hydrolysis by each enzyme indicated that the Zn2+ bimetallo core is an effective mono/di-esterase generalist and that the bimetallo cores were not evolutionarily tuned to prefer their cognate reactions. In contrast, our pruned enzymes were ineffective sulfatases, and this limited promiscuity may have provided a driving force for founding the distinct one-metal-ion branch that contains all known AP superfamily sulfatases. Finally, our pruned enzymes exhibited 107-108-fold phosphotriesterase rate enhancements, despite absence of such enzymes within the AP superfamily. We speculate that the superfamily active-site architecture involved in nucleophile positioning prevents accommodation of the additional triester substituent. Overall, we suggest that catalytic promiscuity, and the ease or difficulty of remodeling and building onto existing protein scaffolds, have greatly influenced the course of enzyme evolution. Uncovering principles and properties of enzyme function, promiscuity, and repurposing provides lessons for engineering new enzymes.

Structural snapshots of the catalytic cycle of the phosphodiesterase Autotaxin.

Autotaxin (ATX) is a secreted phosphodiesterase that produces the signalling lipid lysophosphatidic acid (LPA). The bimetallic active site of ATX is structurally related to the alkaline phosphatase superfamily. Here, we present a new crystal structure of ATX in complex with orthovanadate (ATX-VO5), which binds the Oγ nucleophile of Thr209 and adopts a trigonal bipyramidal conformation, following the nucleophile attack onto the substrate. We have now a portfolio of ATX structures we discuss as intermediates of the catalytic mechanism: the new ATX-VO5 structure; a unique structure where the nucleophile Thr209 is phosphorylated (ATX-pThr). Comparing these to a complex with the LPA product (ATX-LPA) and with a complex with a phosphate ion (ATX-PO4), that represent the Michaelis complex of the reaction, we observe movements of Thr209, changes in the relative displacement of the zinc ions, and a water molecule that likely fulfils the second nucleophilic attack. We propose that ATX follows the associative two-step in-line displacement mechanism.

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.

Understanding Functional Evolution in the Alkaline Phosphatase Superfamily

Over the past 40 years, it has been demonstrated that many enzymes are capable of promiscuous catalytic activities, facilitating the turnover of more than one chemically distinct substrate. This ha ...

Site-Directed Mutagenesis Maps Interactions That Enhance Cognate and Limit Promiscuous Catalysis by an Alkaline Phosphatase Superfamily Phosphodiesterase

Catalytic promiscuity, an evolutionary concept, also provides a powerful tool for gaining mechanistic insights into enzymatic reactions. Members of the alkaline phosphatase (AP) superfamily are highly amenable to such investigation, with several members having been shown to exhibit promiscuous activity for the cognate reactions of other superfamily members. Previous work has shown that nucleotide pyrophosphatase/phosphodiesterase (NPP) exhibits a >10⁶-fold preference for the hydrolysis of phosphate diesters over phosphate monoesters, and that the reaction specificity is reduced 10³-fold when the size of the substituent on the transferred phosphoryl group of phosphate diester substrates is reduced to a methyl group. Here we show additional specificity contributions from the binding pocket for this substituent (herein termed the R' substituent) that account for an additional ~250-fold differential specificity with the minimal methyl substituent. Removal of four hydrophobic side chains suggested on the basis of structural inspection to interact favorably with R' substituents decreases phosphate diester reactivity 10⁴-fold with an optimal diester substrate (R' = 5'-deoxythymidine) and 50-fold with a minimal diester substrate (R' = CH₃). These mutations also enhance the enzyme's promiscuous phosphate monoesterase activity by nearly an order of magnitude, an effect that is traced by mutation to the reduction of unfavorable interactions with the two residues closest to the nonbridging phosphoryl oxygen atoms. The quadruple R' pocket mutant exhibits the same activity toward phosphate diester and phosphate monoester substrates that have identical leaving groups, with substantial rate enhancements of ~10¹¹-fold. This observation suggests that the Zn²⁺ bimetallo core of AP superfamily enzymes, which is equipotent in phosphate monoester and diester catalysis, has the potential to become specialized for the hydrolysis of each class of phosphate esters via addition of side chains that interact with the substrate atoms and substituents that project away from the Zn²⁺ bimetallo core.

Stabilization of Different Types of Transition States in a Single Enzyme Active Site: QM/MM Analysis of Enzymes in the Alkaline Phosphatase Superfamily

The first step for the hydrolysis of a phosphate monoester (pNPP(2-)) in enzymes of the alkaline phosphatase (AP) superfamily, R166S AP and wild-type NPP, is studied using QM/MM simulations based on an approximate density functional theory (SCC-DFTBPR) and a recently introduced QM/MM interaction Hamiltonian. The calculations suggest that similar loose transition states are involved in both enzymes, despite the fact that phosphate monoesters are the cognate substrates for AP but promiscuous substrates for NPP. The computed loose transition states are clearly different from the more synchronous ones previously calculated for diester reactions in the same AP enzymes. Therefore, our results explicitly support the proposal that AP enzymes are able to recognize and stabilize different types of transition states in a single active site. Analysis of the structural features of computed transition states indicates that the plastic nature of the bimetallic site plays a minor role in accommodating multiple types of transition states and that the high degree of solvent accessibility of the AP active site also contributes to its ability to stabilize diverse transition-state structures without the need of causing large structural distortions of the bimetallic motif. The binding mode of the leaving group in the transition state highlights that vanadate may not always be an ideal transition state analog for loose phosphoryl transfer transition states.

Multiple Catalytic Promiscuity: Specificity and Efficiency in the Alkaline Phosphatase Superfamily

‘Promiscuous’ enzymes possess additional activities in addition to their native ones, challenging the textbook adage “one enzyme – one activity”. In particular the question arises, whether efficiency and specialisation are inextricably linked. The observation of strong promiscuous activities in the alkaline phosphatase (AP) superfamily ‐ where one active site can catalyse up to six chemically distinct hydrolytic reactions with promiscuous second order rate accelerations between 109 and 1017‐ suggests that even broadly promiscuous catalysis can be rather efficient. Mechanistic and structural studies are used to delineate how efficiency and specificity can be combined and managed. We demonstrate that crosswise promiscuity relationships in the AP superfamily indicate that an enzyme is ‘pregnant’ with another activity, i.e. has the potential to be mutated or evolved into a new catalyst. These catalysts are multifunctional generalists that have won additional activities, at low trade‐off cost to the other existing activities. In particular, we use picoliter water‐in–oil emulsion droplets produced in microfluidic devices as high‐throughput screening reactors and present a workflow that allows screening of >;106 clones and allows successful selections, where lower throughput approaches have failed.

Harnessing Promiscuity Patterns to Map Evolution in the Alkaline Phosphatase Superfamily

In the past 30 yeas scientists have demonstrated that many enzymes are capable of promiscuous catalytic activity, in which one enzyme catalyses the turnover of chemically distinct substrates. It has been suggested that such promiscuity can play an important role in the evolution of enzyme function. This phenomenon is particularly pronounced in the alkaline phosphatase superfamily, in which the native reaction for one member of the superfamily is often a promiscuous side reaction for another. Moreover, despite the similarity of the substrates involved, the reactions catalysed (cleavage of P-O and S-O bonds) proceed with very different solvation and protonation requirements. We present here a detailed study of the promiscuous catalytic activity of three evolutionarily related but structurally different members of this superfamily: an arylsulfatase from Pseudonomas aeruginosa (PAS) and two phosphonate monoester hydrolases (PMH) from Rhizobium leguminosarum and Burkholderia caryophylli. The latter enzymes are particularly interesting due to the fact that, in addition to their broad promiscuity, they are also capable of catalysing a xenobiotic substrate for which no natural catalyst is known. We demonstrate that subtle changes in active site residues can have a significant impact on mechanism. However, despite such mechanistic variation, the promiscuity simply arises out of the ability of the non-native substrates to exploit the pre-existing electrostatic pre-organization of the active site towards the native substrate. We believe this provides an example of chemically-driven protein evolution within an enzyme superfamily.

Modeling catalytic promiscuity in the alkaline phosphatase superfamily.

In recent years, it has become increasingly clear that promiscuity plays a key role in the evolution of new enzyme function. This finding has helped to elucidate fundamental aspects of molecular evolution. While there has been extensive experimental work on enzyme promiscuity, computational modeling of the chemical details of such promiscuity has traditionally fallen behind the advances in experimental studies, not least due to the nearly prohibitive computational cost involved in examining multiple substrates with multiple potential mechanisms and binding modes in atomic detail with a reasonable degree of accuracy. However, recent advances in both computational methodologies and power have allowed us to reach a stage in the field where we can start to overcome this problem, and molecular simulations can now provide accurate and efficient descriptions of complex biological systems with substantially less computational cost. This has led to significant advances in our understanding of enzyme function and evolution in a broader sense. Here, we will discuss currently available computational approaches that can allow us to probe the underlying molecular basis for enzyme specificity and selectivity, discussing the inherent strengths and weaknesses of each approach. As a case study, we will discuss recent computational work on different members of the alkaline phosphatase superfamily (AP) using a range of different approaches, showing the complementary insights they have provided. We have selected this particular superfamily, as it poses a number of significant challenges for theory, ranging from the complexity of the actual reaction mechanisms involved to the reliable modeling of the catalytic metal centers, as well as the very large system sizes. We will demonstrate that, through current advances in methodologies, computational tools can provide significant insight into the molecular basis for catalytic promiscuity, and, therefore, in turn, the mechanisms of protein functional evolution.