Author Correction: Higher-order epistasis shapes the fitness landscape of a xenobiotic-degrading enzyme.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Higher-order epistasis shapes the fitness landscape of a xenobiotic-degrading enzyme.

Characterizing the adaptive landscapes that encompass the emergence of novel enzyme functions can provide molecular insights into both enzymatic and evolutionary mechanisms. Here, we combine ancestral protein reconstruction with biochemical, structural and mutational analyses to characterize the functional evolution of methyl-parathion hydrolase (MPH), an organophosphate-degrading enzyme. We identify five mutations that are necessary and sufficient for the evolution of MPH from an ancestral dihydrocoumarin hydrolase. In-depth analyses of the adaptive landscapes encompassing this evolutionary transition revealed that the mutations form a complex interaction network, defined in part by higher-order epistasis, that constrained the adaptive pathways available. By also characterizing the adaptive landscapes in terms of their functional activities towards three additional organophosphate substrates, we reveal that subtle differences in the polarity of the substrate substituents drastically alter the network of epistatic interactions. Our work suggests that the mutations function collectively to enable substrate recognition via subtle structural repositioning.

FAD-sequestering proteins protect mycobacteria against hypoxic and oxidative stress.

The ability to persist in the absence of growth triggered by low oxygen levels is a critical process for the survival of mycobacterial species in many environmental niches. MSMEG_5243 (fsq), a gene of unknown function in Mycobacterium smegmatis, is up-regulated in response to hypoxia and regulated by DosRDosS/DosT, an oxygen- and redox-sensing two-component system that is highly conserved in mycobacteria. In this communication, we demonstrate that MSMEG_5243 is a flavin-sequestering protein and henceforth refer to it as Fsq. Using an array of biochemical and structural analyses, we show that Fsq is a member of the diverse superfamily of flavin- and deazaflavin-dependent oxidoreductases (FDORs) and is widely distributed in mycobacterial species. We created a markerless deletion mutant of fsq and demonstrate that fsq is required for cell survival during hypoxia. Using fsq deletion and overexpression, we found that fsq enhances cellular resistance to hydrogen peroxide treatment. The X-ray crystal structure of Fsq, solved to 2.7 Å, revealed a homodimeric organization with FAD bound noncovalently. The Fsq structure also uncovered no potential substrate-binding cavities, as the FAD is fully enclosed, and electrochemical studies indicated that the Fsq:FAD complex is relatively inert and does not share common properties with electron-transfer proteins. Taken together, our results suggest that Fsq reduces the formation of reactive oxygen species (ROS) by sequestering free FAD during recovery from hypoxia, thereby protecting the cofactor from undergoing autoxidation to produce ROS. This finding represents a new paradigm in mycobacterial adaptation to hypoxia.

Evolution of cyclohexadienyl dehydratase from an ancestral solute-binding protein.

The emergence of enzymes through the neofunctionalization of noncatalytic proteins is ultimately responsible for the extraordinary range of biological catalysts observed in nature. Although the evolution of some enzymes from binding proteins can be inferred by homology, we have a limited understanding of the nature of the biochemical and biophysical adaptations along these evolutionary trajectories and the sequence in which they occurred. Here we reconstructed and characterized evolutionary intermediate states linking an ancestral solute-binding protein to the extant enzyme cyclohexadienyl dehydratase. We show how the intrinsic reactivity of a desolvated general acid was harnessed by a series of mutations radiating from the active site, which optimized enzyme-substrate complementarity and transition-state stabilization and minimized sampling of noncatalytic conformations. Our work reveals the molecular evolutionary processes that underlie the emergence of enzymes de novo, which are notably mirrored by recent examples of computational enzyme design and directed evolution.

Mechanistic Basis of pH-Dependent 5-Flucytosine Resistance in Aspergillus fumigatus.

The antifungal drug 5-flucytosine (5FC), a derivative of the nucleobase cytosine, is licensed for the treatment of fungal diseases; however, it is rarely used as a monotherapeutic to treat Aspergillus infection. Despite being potent against other fungal pathogens, 5FC has limited activity against Aspergillus fumigatus when standard in vitro assays are used to determine susceptibility. However, in modified in vitro assays where the pH is set to pH 5, the activity of 5FC increases significantly. Here we provide evidence that fcyB, a gene that encodes a purine-cytosine permease orthologous to known 5FC importers, is downregulated at pH 7 and is the primary factor responsible for the low efficacy of 5FC at pH 7. We also uncover two transcriptional regulators that are responsible for the repression of fcyB and, consequently, mediators of 5FC resistance, the CCAAT binding complex (CBC) and the pH regulatory protein PacC. We propose that the activity of 5FC might be enhanced by the perturbation of factors that repress fcyB expression, such as PacC or other components of the pH-sensing machinery.

Correction: Sterol Biosynthesis and Azole Tolerance Is Governed by the Opposing Actions of SrbA and the CCAAT Binding Complex.

[This corrects the article DOI: 10.1371/journal.ppat.1005775.].

Sterol Biosynthesis and Azole Tolerance Is Governed by the Opposing Actions of SrbA and the CCAAT Binding Complex.

Azole drugs selectively target fungal sterol biosynthesis and are critical to our antifungal therapeutic arsenal. However, resistance to this class of drugs, particularly in the major human mould pathogen Aspergillus fumigatus, is emerging and reaching levels that have prompted some to suggest that there is a realistic probability that they will be lost for clinical use. The dominating class of pan-azole resistant isolates is characterized by the presence of a tandem repeat of at least 34 bases (TR34) within the promoter of cyp51A, the gene encoding the azole drug target sterol C14-demethylase. Here we demonstrate that the repeat sequence in TR34 is bound by both the sterol regulatory element binding protein (SREBP) SrbA, and the CCAAT binding complex (CBC). We show that the CBC acts complementary to SrbA as a negative regulator of ergosterol biosynthesis and show that lack of CBC activity results in increased sterol levels via transcriptional derepression of multiple ergosterol biosynthetic genes including those coding for HMG-CoA-synthase, HMG-CoA-reductase and sterol C14-demethylase. In agreement with these findings, inactivation of the CBC increased tolerance to different classes of drugs targeting ergosterol biosynthesis including the azoles, allylamines (terbinafine) and statins (simvastatin). We reveal that a clinically relevant mutation in HapE (P88L) significantly impairs the binding affinity of the CBC to its target site. We identify that the mechanism underpinning TR34 driven overexpression of cyp51A results from duplication of SrbA but not CBC binding sites and show that deletion of the 34 mer results in lack of cyp51A expression and increased azole susceptibility similar to a cyp51A null mutant. Finally we show that strains lacking a functional CBC are severely attenuated for pathogenicity in a pulmonary and systemic model of aspergillosis.

The role of protein dynamics in the evolution of new enzyme function.

Enzymes must be ordered to allow the stabilization of transition states by their active sites, yet dynamic enough to adopt alternative conformations suited to other steps in their catalytic cycles. The biophysical principles that determine how specific protein dynamics evolve and how remote mutations affect catalytic activity are poorly understood. Here we examine a 'molecular fossil record' that was recently obtained during the laboratory evolution of a phosphotriesterase from Pseudomonas diminuta to an arylesterase. Analysis of the structures and dynamics of nine protein variants along this trajectory, and three rationally designed variants, reveals cycles of structural destabilization and repair, evolutionary pressure to 'freeze out' unproductive motions and sampling of distinct conformations with specific catalytic properties in bi-functional intermediates. This work establishes that changes to the conformational landscapes of proteins are an essential aspect of molecular evolution and that change in function can be achieved through enrichment of preexisting conformational sub-states.

Synthesis, X-ray Structural Characterization, and DFT Calculations of Mononuclear Nickel(II) Complexes Containing Diamine and Methacrylate Ligands.

The mononuclear Ni(II) complexes [Ni(en)2(H2O)2](MAA)2 (1) and [Ni(pn)2(MAA)2] (2), where MAA, en and pn are methacrylate, ethylendiamine and 1,3-propylendiamine, respectively, have been synthesized and characterized by elemental analysis, FT-IR and UV�Vis spectroskopy. Structures of the complexes have been determined by single-crystal X-ray diffraction analyses. In the nickel(II) complexes 1 and 2 nickel(II) ion is six-coordinate and has a distorted octahedral geometry. Ni(II) is bonded to four nitrogen atoms of the two diamines and additionally to two oxygen atoms of aqua ligand in 1, and two oxygen atoms of methacrylate ligands in 2. The theoretical geometries of the studied compounds have been calculated by means of density functional theory (DFT) at the B3LYP/6-311+G(d,p)/LanL2DZ level and considering effective core potential (ECP). The comparison of the results indicates that the employed DFT method yields good agreement with experimental data.

Structure of an Insecticide Sequestering Carboxylesterase from the Disease Vector Culex quinquefasciatus: What Makes an Enzyme a Good Insecticide Sponge?

Carboxylesterase (CBE)-mediated metabolic resistance to organophosphate and carbamate insecticides is a major problem for the control of insect disease vectors, such as the mosquito. The most common mechanism involves overexpression of CBEs that bind to the insecticide with high affinity, thereby sequestering them before they can interact with their target. However, the absence of any structure for an insecticide-sequestering CBE limits our understanding of the molecular basis for this process. We present the first structure of a CBE involved in sequestration, Cqestß21, from the mosquito disease vector Culex quinquefasciatus. Lysine methylation was used to obtain the crystal structure of Cqestß21, which adopts a canonical α/ß-hydrolase fold that has high similarity to the target of organophosphate and carbamate insecticides, acetylcholinesterase. Sequence similarity networks of the insect carboxyl/cholinesterase family demonstrate that CBEs associated with metabolic insecticide resistance across many species share a level of similarity that distinguishes them from a variety of other classes. This is further emphasized by the structural similarities and differences in the binding pocket and active site residues of Cqestß21 and other insect carboxyl/cholinesterases. Stopped-flow and steady-state inhibition studies support a major role for Cqestß21 in organophosphate resistance and a minor role in carbamate resistance. Comparison with another isoform associated with insecticide resistance, Cqestß1, showed both enzymes have similar affinity to insecticides, despite 16 amino acid differences between the two proteins. This provides a molecular understanding of pesticide sequestration by insect CBEs and could facilitate the design of CBE-specific inhibitors to circumvent this resistance mechanism in the future.

A Palladium-Catalyzed Ullmann Cross-Coupling/Reductive Cyclization Route to the Carbazole Natural Products 3-Methyl-9H-carbazole, Glycoborine, Glycozoline, Clauszoline K, Mukonine, and Karapinchamine A.

The title natural products 2-7 have been prepared by reductive cyclization of the relevant 2-arylcyclohex-2-en-1-one (e.g. 20) to the corresponding tetrahydrocarbazole and dehydrogenation (aromatization) of this to give the target carbazole (e.g. 4). Compounds such as 20 were prepared using a palladium-catalyzed Ullmann cross-coupling reaction between the appropriate 2-iodocyclohex-2-en-1-one and o-halonitrobenzene.

A Unified Approach to the Isomeric α-, ß-, γ-, and δ-Carbolines via their 6,7,8,9-Tetrahydro Counterparts.

A cross-coupling/reductive cyclization protocol has been employed in a unified approach to all four carbolines. So, for example, the 2-nitropyridine 8, which is readily prepared through an efficient palladium-catalyzed Ullmann cross-coupling reaction, is reductively cyclized under conventional conditions to give 6,7,8,9-tetrahydro-α-carboline that is itself readily aromatized to give α-carboline (1).

A Total Synthesis of (±)-3-O-Demethylmacronine through Rearrangement of a Precursor Embodying the Haemanthidine Alkaloid Framework.

A total synthesis of the racemic modification, (±)-2, of the tazettine-type alkaloid 3-O-demethylmacronine is described. The key steps are an intramolecular Alder-ene (IMAE) reaction and a lactam-to-lactone rearrangement of tetracycle 13, a compound that embodies the haemanthidine alkaloid framework.

Studies on the Photochemical Rearrangements of Enantiomerically Pure, Polysubstituted, and Variously Annulated Bicyclo[2.2.2]octenones.

A series of enantiomerically pure bicyclo[2.2.2]octenones, including the lactone-annulated system 26, has been prepared by engaging derivatives of an enzymatically derived and homochiral cis-1,2-dihydrocatechol in inter- or intra-molecular Diels-Alder reactions. Systems such as 26 readily participate in photochemically promoted oxa-di-π-methane rearrangement or 1,3-acyl migration processes to give products such as diquinane 34 or mixtures of cyclobutanone 36 and cyclopropane 38, respectively.

The Synthesis of Certain Derivatives and Analogues of (-)- and (+)-Galanthamine and an Assessment of their Capacities to Inhibit Acetylcholine Esterase.

Syntheses of certain di- and mono-oxygenated derivatives (e.g., 2 and 3, respectively) and analogues (e.g., 4, a D-ring monoseco-analogue of 2) of both the (-)- and (+)-enantiomeric forms of the alkaloid galanthamine [(-)-1] are reported. All have been assessed for their capacities to inhibit acetylcholine esterase but, in contrast to the predictions from docking studies, none bind strongly to this enzyme.

Total Synthesis of the Cyclic Carbonate-Containing Natural Product Aspergillusol B from d-(-)-Tartaric Acid.

A total synthesis of compound 3 from d-(-)-tartaric acid is reported, thereby establishing that the structure, including relative stereochemistry, originally assigned to the cyclic carbonate-containing natural product aspergillusol B is correct.

Active Site Desolvation and Thermostability Trade-Offs in the Evolution of Catalytically Diverse Triazine Hydrolases.

The desolvation of ionizable residues in the active sites of enzymes and the subsequent effects on catalysis and thermostability have been studied in model systems, yet little about how enzymes can naturally evolve to include active sites with highly reactive and desolvated charges is known. Variants of triazine hydrolase (TrzN) with significant differences in their active sites have been isolated from different bacterial strains: TrzN from Nocardioides sp. strain MTD22 contains a catalytic glutamate residue (Glu241) that is surrounded by hydrophobic and aromatic second-shell residues (Pro214 and Tyr215), whereas TrzN from Nocardioides sp. strain AN3 has a noncatalytic glutamine residue (Gln241) at an equivalent position, surrounded by hydrophilic residues (Thr214 and His215). To understand how and why these variants have evolved, a series of TrzN mutants were generated and characterized. These results show that desolvation by second-shell residues increases the pKa of Glu241, allowing it to act as a general acid at neutral pH. However, significant thermostability trade-offs are required to incorporate the ionizable Glu241 in the active site and to then enclose it in a hydrophobic microenvironment. Analysis of high-resolution crystal structures shows that there are almost no structural changes to the overall configuration of the active site due to these mutations, suggesting that the changes in activity and thermostability are purely based on the altered electrostatics. The natural evolution of these enzyme isoforms provides a unique system in which to study the fundamental process of charged residue desolvation in enzyme catalysis and its relative contribution to the creation and evolution of an enzyme active site.

Structure and function of an insect α-carboxylesterase (αEsterase7) associated with insecticide resistance.

Insect carboxylesterases from the αEsterase gene cluster, such as αE7 (also known as E3) from the Australian sheep blowfly Lucilia cuprina (LcαE7), play an important physiological role in lipid metabolism and are implicated in the detoxification of organophosphate (OP) insecticides. Despite the importance of OPs to agriculture and the spread of insect-borne diseases, the molecular basis for the ability of α-carboxylesterases to confer OP resistance to insects is poorly understood. In this work, we used laboratory evolution to increase the thermal stability of LcαE7, allowing its overexpression in Escherichia coli and structure determination. The crystal structure reveals a canonical α/ß-hydrolase fold that is very similar to the primary target of OPs (acetylcholinesterase) and a unique N-terminal α-helix that serves as a membrane anchor. Soaking of LcαE7 crystals in OPs led to the capture of a crystallographic snapshot of LcαE7 in its phosphorylated state, which allowed comparison with acetylcholinesterase and rationalization of its ability to protect insects against the effects of OPs. Finally, inspection of the active site of LcαE7 reveals an asymmetric and hydrophobic substrate binding cavity that is well-suited to fatty acid methyl esters, which are hydrolyzed by the enzyme with specificity constants (∼10(6) M(-1) s(-1)) indicative of a natural substrate.

Evolutionary expansion of the amidohydrolase superfamily in bacteria in response to the synthetic compounds molinate and diuron.

The amidohydrolase superfamily has remarkable functional diversity, with considerable structural and functional annotation of known sequences. In microbes, the recent evolution of several members of this family to catalyze the breakdown of environmental xenobiotics is not well understood. An evolutionary transition from binuclear to mononuclear metal ion coordination at the active sites of these enzymes could produce large functional changes such as those observed in nature, but there are few clear examples available to support this hypothesis. To investigate the role of binuclear-mononuclear active-site transitions in the evolution of new function in this superfamily, we have characterized two recently evolved enzymes that catalyze the hydrolysis of the synthetic herbicides molinate (MolA) and phenylurea (PuhB). In this work, the crystal structures, mutagenesis, metal ion analysis, and enzyme kinetics of both MolA and PuhB establish that these enzymes utilize a mononuclear active site. However, bioinformatics and structural comparisons reveal that the closest putative ancestor of these enzymes had a binuclear active site, indicating that a binuclear-mononuclear transition has occurred. These proteins may represent examples of evolution modifying the characteristics of existing catalysts to satisfy new requirements, specifically, metal ion rearrangement leading to large leaps in activity that would not otherwise be possible.

Crystal structure of the mouse interleukin-3 ß-receptor: insights into interleukin-3 binding and receptor activation.

Interleukin-3 (IL-3) is a cytokine secreted by mast cells and activated T-cells known to be an important regulator of differentiation, survival, proliferation and activation of a range of haemopoietic lineages. The effects of IL-3 on target cells are mediated by a transmembrane receptor system composed of a cytokine-specific α-subunit and a ß-subunit, the principal signalling entity. In the mouse, two ß-subunits have co-evolved: a common ß-subunit (ßc) shared between IL-3 and the related cytokines IL-5 and granulocyte/macrophage colony-stimulating factor (GM-CSF); and an IL-3-specific ß-subunit (ßIL-3). ßIL-3 differs from ßc in its specificity for IL-3 and its capacity to bind IL-3 directly in the absence of an α-subunit, and, in the absence of structural information, the basis for these properties has remained enigmatic. In the present study, we have solved the crystal structure of the ßIL-3 ectodomain at 3.45 Å (1 Å=0.1 nm) resolution. This structure provides the first evidence that ßIL-3 adopts an arch-shaped intertwined homodimer with similar topology to the paralogous ßc structure. In contrast with apo-ßc, however, the ligand-binding interface of ßIL-3 appears to pre-exist in a conformation receptive to IL-3 engagement. Molecular modelling of the IL-3-ßIL-3 interface, in conjunction with previous mutational studies, suggests that divergent evolution of both ßIL-3 and IL-3 underlies their unique capacity for direct interaction and specificity.