Evolution of Protein Quaternary Structure in Response to Selective Pressure for Increased Thermostability

Site-specific crosslinking and assembly of tetrameric β-glucuronidase improve glycyrrhizin hydrolysis

In this study, eight nonconserved residues with exposed surfaces and flexible conformations of the homotetrameric PGUS (β-glucuronidase from Aspergillus oryzae Li-3) were identified. Single-point mutation into cysteine enabled the thiol-maleimide reaction and site-specific protein assembly using a two-arm polyethylene glycol (PEG)-maleimide crosslinker (Mal2 ). The Mal2 (1k) (with 1 kDa PEG spacer)-crosslinked PGUS assemblies showed low crosslinking efficiency and unimproved thermostability except for G194C-Mal2 (1k). To improve the crosslinking efficiency, a lengthened crosslinker Mal2 (2k) (with 2 kDa PEG spacer) was used to produce PGUS assembly and a highly improved thermostability was achieved with a half-life of 47.2-169.2 min at 70°C, which is 1.04-3.74 times that of wild type PGUS. It is found that the thermostability of PGUS assembly was closely associated with the formation of inter-tetramer assembly and intratetramer crosslinking, rather than the PEGylation of the enzyme. Therefore, the four-arm PEG-maleimide crosslinker Mal4 (2k) (with 2 kDa PEG spacer) was employed to simultaneously increase the inter-tetramer assembly and intratetramer crosslinking, and the resulting PGUS assemblies showed further improved thermostabilities compared with Mal2 (2k)-crosslinked assemblies. Finally, the application of PGUS assemblies with significantly improved thermostability to the bioconversion of GL proved that the PGUS assembly is a strong catalyst for glycyrrhizin (GL) hydrolysis in industrial applications.

Structural and functional analyses of Pcal_0917, an α-glucosidase from hyperthermophilic archaeon Pyrobaculum calidifontis

Genome analysis of Pyrobaculum calidifontis revealed the presence of α-glucosidase (Pcal_0917) gene. Structural analysis affirmed the presence of signature sequences of Type II α-glucosidases in Pcal_0917. We have heterologously expressed the gene and produced recombinant Pcal_0917 in Escherichia coli. Biochemical characteristics of the recombinant enzyme resembled to that of Type I α-glucosidases, instead of Type II. Recombinant Pcal_0917 existed in a tetrameric form in solution and displayed highest activity at 95 °C and pH 6.0, independent of any metal ions. A short heat-treatment at 90 °C resulted in a 35 % increase in enzyme activity. A slight structural shift was observed by CD spectrometry at this temperature. Half-life of the enzyme was >7 h at 90 °C. Pcal_0917 exhibited apparent V_max values of 1190 ± 5 and 3.9 ± 0.1 U/mg against p-nitrophenyl α-D-glucopyranoside and maltose, respectively. To the best of our knowledge, Pcal_0917 displayed the highest ever reported p-nitrophenyl α-D-glucopyranosidase activity among the characterized counterparts. Moreover, Pcal_0917 displayed transglycosylation activity in addition to α-glucosidase activity. Furthermore, in combination with α-amylase, Pcal_0917 was capable of producing glucose syrup from starch with >40 % glucose content. These properties make Pcal_0917 a potential candidate for starch hydrolyzing industry.

Stability of enzyme immobilized on the nanofluidic channel surface

The lifetime of an enzyme is critical to prevent system failure and optimize maintenance schedules in biological and analytical chemistry. The lifetime metrics of an enzyme can be evaluated from enzyme activity in terms of catalytic cycles per enzyme at various storage times. Trypsin, which is a gold-standard enzyme in proteomics, has been known to decrease activity due to self-digestion. To improve the activity of trypsin, enzyme reactors have developed by immobilizing in micro and nanospace. However, an evaluation method for the catalytic cycle has not been established due to major issues such as nonuniform space, unstable liquid transport, and self-digestion during immobilization in conventional work. To solve these issues, we have previously developed an ultra-fast enzyme reactor with a well-defined nanofabrication method, stable liquid transport, and partial enzyme modification. Here, we aimed to investigate catalytic cycles in a nanochannel. To extend enzyme lifetime efficiently, we have evaluated the optimal immobilization process and catalytic cycles of trypsin. As a result, immobilized enzyme densities by the trypsinogen immobilization process were increased at all concentrations compared to the trypsin immobilization process. To evaluate the lifetime of trypsin, the immobilized enzyme densities and activities were almost the same before and after 72 h of enzyme storage, and the calculated catalytic cycles were 1740. These results indicated that self-digestion of the immobilized enzyme was highly suppressed. Consequently, the reaction efficiency has been evaluated depending on the catalytic cycles from the substrate for the first time, while preventing self-digestion by trypsin.

Tetramerization of GH2 β-Glucuronidases is Essential for Catalyzing the Hydrolysis of the Large Substrate Glycyrrhizin

In this study, structural analysis was employed to identify three hotspot residues that contribute most to the tetramer formation of a glycoside hydrolase family 2 (GH2) β-glucuronidase (GUS) from Aspergillus oryzae Li-3. Single-point mutation at these sites completely disrupted the tetramer structure and abolished the glycyrrhizin (GL)-hydrolyzing activity. Then, the W522A dimer was refactored into a tetramer by disulfide bonding, and partial GL activity was restored. Further saturated mutation showed a strong correlation between the GL activity of the mutants and their tetramer ratios. Molecular simulations were employed to illustrate the critical role of the tetramer interface in maintaining a functional active-site structure. The three highly conserved tetramer-forming residues were finally applied to two other GH2 GUSs for tetramer dissociation and demonstrated the significance of the homotetramerization for GL-hydrolyzing activity of GH2 GUSs. This study lays foundation for engineering GL-hydrolyzing GUSs at the quaternary structure level for function regulations.

Simple mechanisms for the evolution of protein complexity

Proteins are tiny models of biological complexity: specific interactions among their many amino acids cause proteins to fold into elaborate structures, assemble with other proteins into higher-order complexes, and change their functions and structures upon binding other molecules. These complex features are classically thought to evolve via long and gradual trajectories driven by persistent natural selection. But a growing body of evidence from biochemistry, protein engineering, and molecular evolution shows that naturally occurring proteins often exist at or near the genetic edge of multimerization, allostery, and even new folds, so just one or a few mutations can trigger acquisition of these properties. These sudden transitions can occur because many of the physical properties that underlie these features are present in simpler proteins as fortuitous by-products of their architecture. Moreover, complex features of proteins can be encoded by huge arrays of sequences, so they are accessible from many different starting points via many possible paths. Because the bridges to these features are both short and numerous, random chance can join selection as a key factor in explaining the evolution of molecular complexity.

Role of C-terminal domain in a manganese-catalase from Geobacillus thermopakistaniensis

Catalases catalyze the decomposition of hydrogen peroxide into water and oxygen. We have characterized two manganese-catalases from Geobacillus thermopakistaniensis, Cat_Gt and Cat-II_Gt, which exhibited significant variation in their sequence, structure and properties. There was only 23% sequence identity between the two. The striking structural difference was the presence of an extended C-terminal domain in Cat_Gt. Molecular modelling and docking studies revealed that deletion of the C-terminal domain removes non-specific binding, which results in increased substrate affinity. To verify experimentally, a C-terminal truncated version of Cat_Gt, named as Cat_Gt-ΔC, was produced in Escherichia coli and effects of deletion were analyzed. There was no significant difference in optimal pH, optimal temperature and substrate specificity of Cat_Gt and Cat_Gt-ΔC. However, K_m value was reduced from 259 to 157 mM and Cat_Gt-ΔC exhibited ∼1.5-fold higher catalytic efficiency as compared to Cat_Gt. Furthermore, removal of the C-terminal domain converted the tetrameric nature to monomeric, and reduced the thermostability of the truncated protein. These results demonstrate that C-terminal domain of Cat_Gt might have little role in maintaining enzyme function but provides additional structural stability to the protein, which is a desired property for industrial applications.

Molecular cloning, expression in Escherichia coli and structural-functional analysis of a pyruvate kinase from Pyrobaculum calidifontis

This manuscript describes recombinant production, characterization and structural analysis of wild-type and mutant Pcal_0029, a pyruvate kinase from Pyrobaculum calidifontis. Recombinant Pcal_0029 was produced in soluble and highly active form in Escherichia coli. Purified protein exhibited divalent metal-dependent activity which increased with the increase in temperature till 85 °C. Recombinant Pcal_0029 was highly thermostable with no significant loss in activity even after an incubation of 120 min at 100 °C. The enzyme exhibited apparent S_0.5 and V_max values of 0.44 ± 0.05 mM and 840 ± 39 units, respectively, towards phosphoenolpyruvate. These values towards adenosine-5'-diphosphate were 0.5 ± 0.07 mM and 870 ± 26 units, respectively. In silico structural analysis and comparison with the characterized enzymes revealed the presence of eight conserved regions. Two substitutions, K130E and S155G, resulted in a 10-fold decrease in activity. Secondary structure analysis indicated similar structures for the wild-type and the mutant enzymes. Bioinformatics analysis revealed disruption of interatomic interactions and hydrogen bond formation, leading to a decreased flexibility and solvent accessibility, which may have led to decrease in activity. To the best of our knowledge, Pcal_0029 is the most thermostable pyruvate kinase reported so far. Moreover, this is the first study on the role of non-catalytic residues in a pyruvate kinase.

The activity and stability of a cold-active acylaminoacyl peptidase rely on its dimerization by domain swapping

The study of enzymes from extremophiles arouses interest in Protein Science because of the amazing solutions these proteins adopt to cope with extreme conditions. Recently solved, the structure of the psychrophilic acyl aminoacyl peptidase from Sporosarcina psychrophila (SpAAP) pinpoints a mechanism of dimerization unusual for this class of enzymes. The quaternary structure of SpAAP relies on a domain-swapping mechanism involving the N-terminal A1 helix. The A1 helix is conserved among homologous mesophilic and psychrophilic proteins and its deletion causes the formation of a monomeric enzyme, which is inactive and prone to aggregate. Here, we investigate the dimerization mechanism of SpAAP through the analysis of chimeric heterodimers where a protomer lacking the A1 helix combines with a protomer carrying the inactivated catalytic site. Our results indicate that the two active sites are independent, and that a single A1 helix is sufficient to partially recover the quaternary structure and the activity of chimeric heterodimers. Since catalytically competent protomers are unstable and inactive unless they dimerize, SpAAP reveals as an "obligomer" for both structural and functional reasons.

Structural and Biochemical Characterization of a Cold-Active PMGL3 Esterase with Unusual Oligomeric Structure

The gene coding for a novel cold-active esterase PMGL3 was previously obtained from a Siberian permafrost metagenomic DNA library and expressed in Escherichia coli. We elucidated the 3D structure of the enzyme which belongs to the hormone-sensitive lipase (HSL) family. Similar to other bacterial HSLs, PMGL3 shares a canonical α/β hydrolase fold and is presumably a dimer in solution but, in addition to the dimer, it forms a tetrameric structure in a crystal and upon prolonged incubation at 4 °C. Detailed analysis demonstrated that the crystal tetramer of PMGL3 has a unique architecture compared to other known tetramers of the bacterial HSLs. To study the role of the specific residues comprising the tetramerization interface of PMGL3, several mutant variants were constructed. Size exclusion chromatography (SEC) analysis of D7N, E47Q, and K67A mutants demonstrated that they still contained a portion of tetrameric form after heat treatment, although its amount was significantly lower in D7N and K67A compared to the wild type. Moreover, the D7N and K67A mutants demonstrated a 40 and 60% increase in the half-life at 40 °C in comparison with the wild type protein. K_m values of these mutants were similar to that of the wt PMGL3. However, the catalytic constants of the E47Q and K67A mutants were reduced by ~40%.

Altered conformational sampling along an evolutionary trajectory changes the catalytic activity of an enzyme

Several enzymes are known to have evolved from non-catalytic proteins such as solute-binding proteins (SBPs). Although attention has been focused on how a binding site can evolve to become catalytic, an equally important question is: how do the structural dynamics of a binding protein change as it becomes an efficient enzyme? Here we performed a variety of experiments, including propargyl-DO3A-Gd(III) tagging and double electron-electron resonance (DEER) to study the rigid body protein dynamics of reconstructed evolutionary intermediates to determine how the conformational sampling of a protein changes along an evolutionary trajectory linking an arginine SBP to a cyclohexadienyl dehydratase (CDT). We observed that primitive dehydratases predominantly populate catalytically unproductive conformations that are vestiges of their ancestral SBP function. Non-productive conformational states, including a wide-open state, are frozen out of the conformational landscape via remote mutations, eventually leading to extant CDT that exclusively samples catalytically relevant compact states. These results show that remote mutations can reshape the global conformational landscape of an enzyme as a mechanism for increasing catalytic activity.

A conserved π-helix plays a key role in thermoadaptation of catalysis in the glycoside hydrolase family 4

Here, we characterize the role of a π-helix in the molecular mechanisms underlying thermoadaptation in the glycoside hydrolase family 4 (GH4). The interspersed π-helix present in a subgroup is evolutionarily related to a conserved α-helix in other orthologs by a single residue insertion/deletion event. The insertional residue, Phe407, in a hyperthermophilic α-glucuronidase, makes specific interactions across the inter-subunit interface. In order to establish the sequence-structure-stability implications of the π-helix, the wild-type and the deletion variant (Δ407) were characterized. The variant showed a significant lowering of melting temperature and optimum temperature for the highest activity. Crystal structures of the proteins show a transformation of the π-helix to a continuous α-helix in the variant, identical to that in orthologs lacking this insertion. Thermodynamic parameters were determined from stability curves representing the temperature dependence of unfolding free energy. Though the proteins display maximum stabilities at similar temperatures, a higher melting temperature in the wild-type is achieved by a combination of higher enthalpy and lower heat capacity of unfolding. Comparisons of the structural changes, and the activity and thermodynamic profiles allow us to infer that specific non-covalent interactions, and the existence of residual structure in the unfolded state, are crucial determinants of its thermostability. These features permit the enzyme to balance the preservation of structure at a higher temperature with the thermodynamic stability required for optimum catalysis.

Shotgun proteomics reveals putative polyesterases in the secretome of the rock-inhabiting fungus Knufia chersonesos

Knufia chersonesos is an ascomycotal representative of black fungi, a morphological group of polyextremotolerant melanotic fungi, whose ability to resort to recalcitrant carbon sources makes it an interesting candidate for degradation purposes. A secretome screening towards polyesterases was carried out for the fungus and its non-melanized mutant, grown in presence of the synthetic copolyester Polybutylene adipate terephthalate (PBAT) as additional or sole carbon source, and resulted in the identification of 37 esterolytic and lipolytic enzymes across the established cultivation conditions. Quantitative proteomics allowed to unveil 9 proteins being constitutively expressed at all conditions and 7 which were instead detected as up-regulated by PBAT exposure. Protein functional analysis and structure prediction indicated similarity of these enzymes to microbial polyesterases of known biotechnological use such as MHETase from Ideonella sakaiensis and CalA from Candida albicans. For both strains, PBAT hydrolysis was recorded at all cultivation conditions and primarily the corresponding monomers were released, which suggests degradation to the polymer's smallest building block. The work presented here aims to demonstrate how investigations of the secretome can provide new insights into the eco-physiology of polymer degrading fungi and ultimately aid the identification of novel enzymes with potential application in polymer processing, recycling and degradation.

Two Highly Similar Chitinases from Marine Vibrio Species have Different Enzymatic Properties

Chitinase, as one of the most important extracellular enzymes in the marine environment, has great ecological and applied values. In this study, two chitinases (Chi1557 and Chi4668) with 97.33% amino acid sequences identity were individually found in Vibrio rotiferianus and Vibrio harveyi. They both were encoding by 561 amino acids, but differed in 15 amino acids and showed different enzymatic properties. The optimal temperature and pH ranges were 45-50 °C and pH 5.0-7.0 for Chi1557, while ~50 °C and pH 3.0-6.0 for Chi4668. K+, Mg2+, and EDTA increased the enzymatic activity of Chi4668 significantly, yet these factors were inhibitory to Chi1557. Moreover, Chi1557 degraded colloidal chitin to produce (GlcNAc)2 and minor GlcNAc, whereas Chi4668 produce (GlcNAc)2 with minor (GlcNAc)3 and (GlcNAc)4. The Kcat/Km of Chi4668 was ~4.7 times higher than that of Chi1557, indicating that Chi4668 had stronger catalytic activity than Chi1557. Furthermore, site-directed mutagenesis was performed on Chi1557 focusing on seven conserved amino acid residues of family GH18 chitinases. Chi1557 was almost completely inactive after Glu154, Gln219, Tyr221, or Trp312 was individually mutated, retained ~50% activity after Tyr37 was mutated, and increased two times activity after Asp152 was mutated, indicating that these six amino acids were key sites for Chi1557.

The interdimeric interface controls function and stability of Ureaplasma urealiticum methionine S-adenosyltransferase

Methionine S-adenosyltransferases (MATs) are predominantly homotetramers, comprised of dimers of dimers. The larger, highly conserved intradimeric interface harbors two active sites, making the dimer the obligatory functional unit. However, functionality of the smaller, more diverged, and recently evolved interdimeric interface is largely unknown. Here, we show that the interdimeric interface of Ureaplasmaurealiticum MAT has evolved to control the catalytic activity and structural integrity of the homotetramer in response to product accumulation. When all four active sites are occupied with the product, S-adenosylmethionine (SAM), binding of four additional SAM molecules to the interdimeric interface prompts a ∼45° shift in the dimer orientation and a concomitant ∼60% increase in the interface area. This rearrangement inhibits the enzymatic activity by locking the flexible active site loops in a closed state and renders the tetramer resistant to proteolytic degradation. Our findings suggest that the interdimeric interface of MATs is subject to rapid evolutionary changes that tailor the molecular properties of the entire homotetramer to the specific needs of the organism.

Overcoming insecticide resistance through computational inhibitor design

Insecticides allow control of agricultural pests and disease vectors and are vital for global food security and health. The evolution of resistance to insecticides, such as organophosphates (OPs), is a serious and growing concern. OP resistance often involves sequestration or hydrolysis of OPs by carboxylesterases. Inhibiting carboxylesterases could, therefore, restore the effectiveness of OPs for which resistance has evolved. Here, we use covalent virtual screening to produce nano-/picomolar boronic acid inhibitors of the carboxylesterase αE7 from the agricultural pest Lucilia cuprina as well as a common Gly137Asp αE7 mutant that confers OP resistance. These inhibitors, with high selectivity against human acetylcholinesterase and low to no toxicity in human cells and in mice, act synergistically with the OPs diazinon and malathion to reduce the amount of OP required to kill L. cuprina by up to 16-fold and abolish resistance. The compounds exhibit broad utility in significantly potentiating another OP, chlorpyrifos, against the common pest, the peach-potato aphid (Myzus persicae). These compounds represent a solution to OP resistance as well as to environmental concerns regarding overuse of OPs, allowing significant reduction of use without compromising efficacy.

Structure of a hyperthermostable dimeric archaeal Rubisco fromHyperthermus butylicus

The crystal structure of ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) from the hyperthermophilic archaeonHyperthermus butylicusis presented at 1.8 Å resolution. Previous structures of archaeal Rubisco have been found to assemble into decamers, and this oligomerization was thought to be required for a highly thermally stable enzyme. In the current study,H. butylicusRubisco is shown to exist as a dimer in solution, yet has a thermal denaturation midpoint of 114°C, suggesting that high thermal stability can be achieved without an increased oligomeric state. This increased thermal stability appears to be due to an increased number of electrostatic interactions within the monomeric subunit. As such,H. butylicusRubisco presents a well characterized system in which to investigate the role of assembly and thermal stability in enzyme function.

Structural Insights into the Molecular Evolution of the Archaeal Exo-β-d-Glucosaminidase

The archaeal exo-β-d-glucosaminidase (GlmA), a thermostable enzyme belonging to the glycosidase hydrolase (GH) 35 family, hydrolyzes chitosan oligosaccharides into monomer glucosamines. GlmA is a novel enzyme in terms of its primary structure, as it is homologous to both GH35 and GH42 β-galactosidases. The catalytic mechanism of GlmA is not known. Here, we summarize the recent reports on the crystallographic analysis of GlmA. GlmA is a homodimer, with each subunit comprising three distinct domains: a catalytic TIM-barrel domain, an α/β domain, and a β1 domain. Surprisingly, the structure of GlmA presents features common to GH35 and GH42 β-galactosidases, with the domain organization resembling that of GH42 β-galactosidases and the active-site architecture resembling that of GH35 β-galactosidases. Additionally, the GlmA structure also provides critical information about its catalytic mechanism, in particular, on how the enzyme can recognize glucosamine. Finally, we postulate an evolutionary pathway based on the structure of an ancestor GlmA to extant GH35 and GH42 β-galactosidases.

N-terminal domain replacement changes an archaeal monoacylglycerol lipase into a triacylglycerol lipase

BACKGROUND

Lipolytic enzymes of hyperthermophilic archaea generally prefer small carbon chain fatty acid esters (C₂-C₁₂) and are categorized as esterases. However, a few have shown activity with long-chain fatty acid esters, but none of them have been classified as a true lipase except a lipolytic enzyme AFL from Archaeglobus fulgidus. Thus, our main objective is to engineer an archaeal esterase into a true thermostable lipase for industrial applications. Lipases which hydrolyze long-chain fatty acid esters display an interfacial activation mediated by the lid domain which lies over active site and switches to open conformation at the oil-water interface. Lid domains modulate enzyme activities, substrate specificities, and stabilities which have been shown by protein engineering and mutational analyses. Here, we report engineering of an uncharacterized monoacylglycerol lipase (TON-LPL) from an archaeon Thermococcus onnurineus (strain NA1) into a triacylglycerol lipase (rc-TGL) by replacing its 61 N-terminus amino acid residues with 118 residues carrying lid domain of a thermophilic fungal lipase-Thermomyces lanuginosus (TLIP).

RESULTS

TON-LPL and rc-TGL were cloned and overexpressed in E. coli, and the proteins were purified by Ni-NTA affinity chromatography for biochemical studies. Both enzymes were capable of hydrolyzing various monoglycerides and shared the same optimum pH of 7.0. However, rc-TGL showed a significant decrease of 10 °C in its optimum temperature (T_opt). The far UV-CD spectrums were consistent with a well-folded α/β-hydrolase fold for both proteins, but gel filtration chromatography revealed a change in quaternary structure from trimer (TON-LPL) to monomer (rc-TGL). Seemingly, the difference in the oligomeric state of rc-TGL may be linked to a decrease in temperature optimum. Nonetheless, rc-TGL hydrolyzed triglycerides and castor oil, while TON-LPL was not active with these substrates.

CONCLUSIONS

Here, we have confirmed the predicted esterase activity of TON-LPL and also performed the lid engineering on TON-LPL which effectively expanded its substrate specificity from monoglycerides to triglycerides. This approach provides a way to engineer other hyperthermophilic esterases into industrially suitable lipases by employing N-terminal domain replacement. The immobilized preparation of rc-TGL has shown significant activity with castor oil and has a potential application in castor oil biorefinery to obtain value-added chemicals.

Contribution of the Oligomeric State to the Thermostability of Isoenzyme 3 from Candida rugosa

Thermophilic proteins have evolved different strategies to maintain structure and function at high temperatures; they have large, hydrophobic cores, and feature increased electrostatic interactions, with disulfide bonds, salt-bridging, and surface charges. Oligomerization is also recognized as a mechanism for protein stabilization to confer a thermophilic adaptation. Mesophilic proteins are less thermostable than their thermophilic homologs, but oligomerization plays an important role in biological processes on a wide variety of mesophilic enzymes, including thermostabilization. The mesophilic yeast Candida rugosa contains a complex family of highly related lipase isoenzymes. Lip3 has been purified and characterized in two oligomeric states, monomer (mLip3) and dimer (dLip3), and crystallized in a dimeric conformation, providing a perfect model for studying the effects of homodimerization on mesophilic enzymes. We studied kinetics and stability at different pHs and temperatures, using the response surface methodology to compare both forms. At the kinetic level, homodimerization expanded Lip3 specificity (serving as a better catalyst on soluble substrates). Indeed, dimerization increased its thermostability by more than 15 °C (maximum temperature for dLip3 was out of the experimental range; >50 °C), and increased the pH stability by nearly one pH unit, demonstrating that oligomerization is a viable strategy for the stabilization of mesophilic enzymes.

Multiple mutations and overexpression of the MdaE7 carboxylesterase gene associated with male-linked malathion resistance in housefly, Musca domestica (Diptera: Muscidae)

Two unique housefly strains, MSS and N-MRS, were selected and used to clarify mechanisms of sex-associated malathion resistance in the housefly, Musca domestica. Compared with the lab-susceptible CSS strain, susceptible females and resistant males were observed in the malathion-susceptible MSS strain, while the malathion-resistant near-isogenic line, N-MRS, achieved similar resistance level between genders. Significant synergistic effect of the esterase-inhibitor DEF on resistant houseflies pointed to the important involvement of esterase in this specific malathion resistance. Examination of the carboxylesterase gene MdαE7 in malathion resistant housefly populations found seven, non-synonymous SNP mutations (Ser²⁵⁰-Thr, Trp²⁵¹-Ser, Met³⁰³-Ile, Leu³⁵⁴-Phe, Ser³⁵⁷-Leu, Trp³⁷⁸-Arg and Ser³⁸³-Thr), not found in susceptible houseflies, revealing a strong correlation between these mutations and the development of malathion resistance. Further genetic analysis conducted with bioassays by topical application and nucleotide polymorphism detection provided a first line of molecular evidence for a linkage between a male-determining factor and MdαE7 gene in the MSS and N-MRS males. This linkage results in a much higher level of malathion resistance for males than females in the MSS strain. Lastly, quantitative real-time PCR showed that MdαE7 was over expressed in the resistant strain due to the increased transcription level of mRNA rather than gene duplication.

Proteome Evolution of Deep-Sea Hydrothermal Vent Alvinellid Polychaetes Supports the Ancestry of Thermophily and Subsequent Adaptation to Cold in Some Lineages

Temperature, perhaps more than any other environmental factor, is likely to influence the evolution of all organisms. It is also a very interesting factor to understand how genomes are shaped by selection over evolutionary timescales, as it potentially affects the whole genome. Among thermophilic prokaryotes, temperature affects both codon usage and protein composition to increase the stability of the transcriptional/translational machinery, and the resulting proteins need to be functional at high temperatures. Among eukaryotes less is known about genome evolution, and the tube-dwelling worms of the family Alvinellidae represent an excellent opportunity to test hypotheses about the emergence of thermophily in ectothermic metazoans. The Alvinellidae are a group of worms that experience varying thermal regimes, presumably having evolved into these niches over evolutionary times. Here we analyzed 423 putative orthologous loci derived from 6 alvinellid species including the thermophilic Alvinella pompejana and Paralvinella sulfincola. This comparative approach allowed us to assess amino acid composition, codon usage, divergence, direction of residue changes and the strength of selection along the alvinellid phylogeny, and to design a new eukaryotic thermophilic criterion based on significant differences in the residue composition of proteins. Contrary to expectations, the alvinellid ancestor of all present-day species seems to have been thermophilic, a trait subsequently maintained by purifying selection in lineages that still inhabit higher temperature environments. In contrast, lineages currently living in colder habitats likely evolved under selective relaxation, with some degree of positive selection for low-temperature adaptation at the protein level.

Functional determinants of protein assembly into homomeric complexes

Approximately half of proteins with experimentally determined structures can interact with other copies of themselves and assemble into homomeric complexes, the overwhelming majority of which (>96%) are symmetric. Although homomerisation is often assumed to a functionally beneficial result of evolutionary selection, there has been little systematic analysis of the relationship between homomer structure and function. Here, utilizing the large numbers of structures and functional annotations now available, we have investigated how proteins that assemble into different types of homomers are associated with different biological functions. We observe that homomers from different symmetry groups are significantly enriched in distinct functions, and can often provide simple physical and geometrical explanations for these associations in regards to substrate recognition or physical environment. One of the strongest associations is the tendency for metabolic enzymes to form dihedral complexes, which we suggest is closely related to allosteric regulation. We provide a physical explanation for why allostery is related to dihedral complexes: it allows for efficient propagation of conformational changes across isologous (i.e. symmetric) interfaces. Overall we demonstrate a clear relationship between protein function and homomer symmetry that has important implications for understanding protein evolution, as well as for predicting protein function and quaternary structure.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.