Approaching the catalytic mechanism of protein lysine methyltransferases by biochemical and simulation techniques.

Protein lysine methyltransferases (PKMTs) transfer up to three methyl groups to the side chains of lysine residues in proteins and fulfill important regulatory functions by controlling protein stability, localization and protein/protein interactions. The methylation reactions are highly regulated, and aberrant methylation of proteins is associated with several types of diseases including neurologic disorders, cardiovascular diseases, and various types of cancer. This review describes novel insights into the catalytic machinery of various PKMTs achieved by the combined application of biochemical experiments and simulation approaches during the last years, focusing on clinically relevant and well-studied enzymes of this group like DOT1L, SMYD1-3, SET7/9, G9a/GLP, SETD2, SUV420H2, NSD1/2, different MLLs and EZH2. Biochemical experiments have unraveled many mechanistic features of PKMTs concerning their substrate and product specificity, processivity and the effects of somatic mutations observed in PKMTs in cancer cells. Structural data additionally provided information about the substrate recognition, enzyme-substrate complex formation, and allowed for simulations of the substrate peptide interaction and mechanism of PKMTs with atomistic resolution by molecular dynamics and hybrid quantum mechanics/molecular mechanics methods. These simulation technologies uncovered important mechanistic details of the PKMT reaction mechanism including the processes responsible for the deprotonation of the target lysine residue, essential conformational changes of the PKMT upon substrate binding, but also rationalized regulatory principles like PKMT autoinhibition. Further developments are discussed that could bring us closer to a mechanistic understanding of catalysis of this important class of enzymes in the near future. The results described here illustrate the power of the investigation of enzyme mechanisms by the combined application of biochemical experiments and simulation technologies.

EnzymeML: seamless data flow and modeling of enzymatic data.

The design of biocatalytic reaction systems is highly complex owing to the dependency of the estimated kinetic parameters on the enzyme, the reaction conditions, and the modeling method. Consequently, reproducibility of enzymatic experiments and reusability of enzymatic data are challenging. We developed the XML-based markup language EnzymeML to enable storage and exchange of enzymatic data such as reaction conditions, the time course of the substrate and the product, kinetic parameters and the kinetic model, thus making enzymatic data findable, accessible, interoperable and reusable (FAIR). The feasibility and usefulness of the EnzymeML toolbox is demonstrated in six scenarios, for which data and metadata of different enzymatic reactions are collected and analyzed. EnzymeML serves as a seamless communication channel between experimental platforms, electronic lab notebooks, tools for modeling of enzyme kinetics, publication platforms and enzymatic reaction databases. EnzymeML is open and transparent, and invites the community to contribute. All documents and codes are freely available at https://enzymeml.org .

Processing of a plant peptide hormone precursor facilitated by posttranslational tyrosine sulfation.

Most peptide hormones and growth factors are matured from larger inactive precursor proteins by proteolytic processing and further posttranslational modification. Whether or how posttranslational modifications contribute to peptide bioactivity is still largely unknown. We address this question here for TWS1 (Twisted Seed 1), a peptide regulator of embryonic cuticle formation in Arabidopsis thaliana. Using synthetic peptides encompassing the N- and C-terminal processing sites and the recombinant TWS1 precursor as substrates, we show that the precursor is cleaved by the subtilase SBT1.8 at both the N and the C termini of TWS1. Recognition and correct processing at the N-terminal site depended on sulfation of an adjacent tyrosine residue. Arginine 302 of SBT1.8 was found to be required for sulfotyrosine binding and for accurate processing of the TWS1 precursor. The data reveal a critical role for posttranslational modification, here tyrosine sulfation of a plant peptide hormone precursor, in mediating processing specificity and peptide maturation.

The T1150A cancer mutant of the protein lysine dimethyltransferase NSD2 can introduce H3K36 trimethylation.

Protein lysine methyltransferases (PKMTs) play essential roles in gene expression regulation and cancer development. Somatic mutations in PKMTs are frequently observed in cancer cells. In biochemical experiments, we show here that the NSD1 mutations Y1971C, R2017Q, and R2017L observed mostly in solid cancers are catalytically inactive suggesting that NSD1 acts as a tumor suppressor gene in these tumors. In contrast, the frequently observed T1150A in NSD2 and its T2029A counterpart in NSD1, both observed in leukemia, are hyperactive and introduce up to three methyl groups in H3K36 in biochemical and cellular assays, while wildtype NSD2 and NSD1 only introduce up to two methyl groups. In Molecular Dynamics simulations, we determined key mechanistic and structural features controlling the product specificity of this class of enzymes. Simulations with NSD2 revealed that H3K36me3 formation is possible due to an enlarged active site pocket of T1150A and loss of direct contacts of T1150 to critical residues which regulate the product specificity of NSD2. Bioinformatic analyses of published data suggested that the generation of H3K36me3 by NSD2 T1150A could alter gene regulation by antagonizing H3K27me3 finally leading to the upregulation of oncogenes.

Plastics degradation by hydrolytic enzymes: The plastics-active enzymes database-PAZy.

Petroleum-based plastics are durable and accumulate in all ecological niches. Knowledge on enzymatic degradation is sparse. Today, less than 50 verified plastics-active enzymes are known. First examples of enzymes acting on the polymers polyethylene terephthalate (PET) and polyurethane (PUR) have been reported together with a detailed biochemical and structural description. Furthermore, very few polyamide (PA) oligomer active enzymes are known. In this article, the current known enzymes acting on the synthetic polymers PET and PUR are briefly summarized, their published activity data were collected and integrated into a comprehensive open access database. The Plastics-Active Enzymes Database (PAZy) represents an inventory of known and experimentally verified enzymes that act on synthetic fossil fuel-based polymers. Almost 3000 homologs of PET-active enzymes were identified by profile hidden Markov models. Over 2000 homologs of PUR-active enzymes were identified by BLAST. Based on multiple sequence alignments, conservation analysis identified the most conserved amino acids, and sequence motifs for PET- and PUR-active enzymes were derived.

Expansin Engineering Database: A navigation and classification tool for expansins and homologues.

Expansins have the remarkable ability to loosen plant cell walls and cellulose material without showing catalytic activity and therefore have potential applications in biomass degradation. To support the study of sequence-structure-function relationships and the search for novel expansins, the Expansin Engineering Database (ExED, https://exed.biocatnet.de) collected sequence and structure data on expansins from Bacteria, Fungi, and Viridiplantae, and expansin-like homologues such as carbohydrate binding modules, glycoside hydrolases, loosenins, swollenins, cerato-platanins, and EXPNs. Based on global sequence alignment and protein sequence network analysis, the sequences are highly diverse. However, many similarities were found between the expansin domains. Newly created profile hidden Markov models of the two expansin domains enable standard numbering schemes, comprehensive conservation analyses, and genome annotation. Conserved key amino acids in the expansin domains were identified, a refined classification of expansins and carbohydrate binding modules was proposed, and new sequence motifs facilitate the search of novel candidate genes and the engineering of expansins.

Multicopper oxidases: modular structure, sequence space, and evolutionary relationships.

Multicopper oxidases (MCOs) use copper ions as cofactors to oxidize a variety of substrates while reducing oxygen to water. MCOs have been identified in various taxa, with notable occurrences in fungi. The role of these fungal MCOs in lignin degradation sparked an interest due to their potential for application in biofuel production and various other industries. MCOs consist of different protein domains, which led to their classification into two-, three-, and six-domain MCOs. The previously established Laccase and Multicopper Oxidase Engineering Database (https://lcced.biocatnet.de) was updated and now includes 51 058 sequences and 229 structures of MCOs. Sequences and structures of all MCOs were systematically compared. All MCOs consist of cupredoxin-like domains. Two-domain MCOs are formed by the N- and C-terminal domain (domain N and C), while three-domain MCOs have an additional domain (M) in between, homologous to domain C. The six-domain MCOs consist of alternating domains N and C, each three times. Two standard numbering schemes were developed for the copper-binding domains N and C, which facilitated the identification of conserved positions and a comparison to previously reported results from mutagenesis studies. Two sequence motifs for the copper binding sites were identified per domain. Their modularity, depending on the placement of the T1-copper binding site, was demonstrated. Protein sequence networks showed relationships between two- and three-domain MCOs, allowing for family-specific annotation and inference of evolutionary relationships.

Systematic Evaluation of Imine-Reducing Enzymes: Common Principles in Imine Reductases, ß-Hydroxy Acid Dehydrogenases, and Short-Chain Dehydrogenases/ Reductases.

The enzymatic, asymmetric reduction of imines is catalyzed by imine reductases (IREDs), members of the short-chain dehydrogenase/reductase (SDR) family, and ß-hydroxy acid dehydrogenase (ßHAD) variants. Systematic evaluation of the structures and substrate-binding sites of the three enzyme families has revealed four common principles for imine reduction: structurally conserved cofactor-binding domains; tyrosine, aspartate, or glutamate as proton donor; at least four characteristic flanking residues that adapt the donor's pKa and polarize the substrate; and a negative electrostatic potential in the substrate-binding site to stabilize the transition state. As additional catalytically relevant positions, we propose alternative proton donors in IREDs and ßHADs as well as proton relays in IREDs, ßHADs, and SDRs. The functional role of flanking residues was experimentally confirmed by alanine scanning of the imine-reducing SDR from Zephyranthes treatiae. Mutating the "gatekeeping" phenylalanine at standard position 200 resulted in a tenfold increase in imine-reducing activity.

Engineering of Thermostable ß-Hydroxyacid Dehydrogenase for the Asymmetric Reduction of Imines.

The ß-hydroxyacid dehydrogenase from Thermocrinus albus (Ta-ßHAD), which catalyzes the NADP+ -dependent oxidation of ß-hydroxyacids, was engineered to accept imines as substrates. The catalytic activity of the proton-donor variant K189D was further increased by the introduction of two nonpolar flanking residues (N192 L, N193 L). Engineering the putative alternative proton donor (D258S) and the gate-keeping residue (F250 A) led to a switched substrate specificity as compared to the single and triple variants. The two most active Ta-ßHAD variants were applied to biocatalytic asymmetric reductions of imines at elevated temperatures and enabled enhanced product formation at a reaction temperature of 50 °C.

Crossing the Border: From Keto- to Imine Reduction in Short-Chain Dehydrogenases/Reductases.

The family of NAD(P)H-dependent short-chain dehydrogenases/reductases (SDRs) comprises numerous biocatalysts capable of C=O or C=C reduction. The highly homologous noroxomaritidine reductase (NR) from Narcissus sp. aff. pseudonarcissus and Zt_SDR from Zephyranthes treatiae, however, are SDRs with an extended imine substrate scope. Comparison with a similar SDR from Asparagus officinalis (Ao_SDR) exhibiting keto-reducing activity, yet negligible imine-reducing capability, and mining the Short-Chain Dehydrogenase/Reductase Engineering Database indicated that NR and Zt_SDR possess a unique active-site composition among SDRs. Adapting the active site of Ao_SDR accordingly improved its imine-reducing capability. By applying the same strategy, an unrelated SDR from Methylobacterium sp. 77 (M77_SDR) with distinct keto-reducing activity was engineered into a promiscuous enzyme with imine-reducing activity, thereby confirming that the ability to reduce imines can be rationally introduced into members of the "classical" SDR enzyme family. Thus, members of the SDR family could be a promising starting point for protein approaches to generate new imine-reducing enzymes.

Navigating within thiamine diphosphate-dependent decarboxylases: Sequences, structures, functional positions, and binding sites.

Thiamine diphosphate-dependent decarboxylases catalyze both cleavage and formation of CC bonds in various reactions, which have been assigned to different homologous sequence families. This work compares 53 ThDP-dependent decarboxylases with known crystal structures. Both sequence and structural information were analyzed synergistically and data were analyzed for global and local properties by means of statistical approaches (principle component analysis and principal coordinate analysis) enabling complexity reduction. The different results obtained both locally and globally, that is, individual positions compared with the overall protein sequence or structure, revealed challenges in the assignment of separated homologous families. The methods applied herein support the comparison of enzyme families and the identification of functionally relevant positions. The findings for the family of ThDP-dependent decarboxylases underline that global sequence identity alone is not sufficient to distinguish enzyme function. Instead, local sequence similarity, defined by comparisons of structurally equivalent positions, allows for a better navigation within several groups of homologous enzymes. The differentiation between homologous sequences is further enhanced by taking structural information into account, such as BioGPS analysis of the active site properties or pairwise structural superimpositions. The methods applied herein are expected to be transferrable to other enzyme families, to facilitate family assignments for homologous protein sequences.

The Short-chain Dehydrogenase/Reductase Engineering Database (SDRED): A classification and analysis system for a highly diverse enzyme family.

The Short-chain Dehydrogenases/Reductases Engineering Database (SDRED) covers one of the largest known protein families (168 150 proteins). Assignment to the superfamilies of Classical and Extended SDRs was achieved by global sequence similarity and by identification of family-specific sequence motifs. Two standard numbering schemes were established for Classical and Extended SDRs that allow for the determination of conserved amino acid residues, such as cofactor specificity determining positions or superfamily specific sequence motifs. The comprehensive sequence dataset of the SDRED facilitates the refinement of family-specific sequence motifs. The glycine-rich motifs for Classical and Extended SDRs were refined to improve the precision of superfamily classification. In each superfamily, the majority of sequences formed a tightly connected sequence network and belonged to a large homologous family. Despite their different sequence motifs and their different sequence length, the two sequence networks of Classical and Extended SDRs are not separate, but connected by edges at a threshold of 40% sequence similarity, indicating that all SDRs belong to a large, connected network. The SDRED is accessible at https://sdred.biocatnet.de/.

Alteration of the Route to Menaquinone towards Isochorismate-Derived Metabolites.

Chorismate and isochorismate constitute branch-point intermediates in the biosynthesis of many aromatic metabolites in microorganisms and plants. To obtain unnatural compounds, we modified the route to menaquinone in Escherichia coli. We propose a model for the binding of isochorismate to the active site of MenD ((1R,2S, 5S,6S)-2-succinyl-5-enolpyruvyl-6-hydroxycyclohex-3-ene-1-carboxylate (SEPHCHC) synthase) that explains the outcome of the native reaction with α-ketoglutarate. We have rationally designed variants of MenD for the conversion of several isochorismate analogues. The double-variant Asn117Arg-Leu478Thr preferentially converts (5S,6S)-5,6-dihydroxycyclohexa-1,3-diene-1-carboxylate (2,3-trans-CHD), the hydrolysis product of isochorismate, with a >70-fold higher ratio than that for the wild type. The single-variant Arg107Ile uses (5S,6S)-6-amino-5-hydroxycyclohexa-1,3-diene-1-carboxylate (2,3-trans-CHA) as substrate with >6-fold conversion compared to wild-type MenD. The novel compounds have been made accessible in vivo (up to 5.3 g L-1 ). Unexpectedly, as the identified residues such as Arg107 are highly conserved (>94 %), some of the designed variations can be found in wild-type SEPHCHC synthases from other bacteria (Arg107Lys, 0.3 %). This raises the question for the possible natural occurrence of as yet unexplored branches of the shikimate pathway.

Chorismatases - the family is growing.

Chorismatases catalyse the cleavage of chorismate, yielding (dihydroxy-)benzoate derivatives, which often constitute starter units for pharmaceutically relevant secondary metabolites. Depending on their products, chorismatases have been classified into three different subfamilies. These can be assigned using a set of amino acid residues in the active site. Here, we describe five new chorismatases, two of them members of a new subfamily, which has been discovered through correlation analysis of homologous protein sequences. The enzymes from the new subfamily produce exclusively 4-hydroxybenzoate, the same compound as produced by the structurally unrelated chorismate lyases. This showcase of convergent evolution is an example of the existence of more than one pathway to central building blocks. In contrast to chorismate lyases, however, chorismatases do not suffer from product inhibition (up to 2 mM 4-HBA), while the remaining kinetic parameters are in the same range; this makes them an interesting alternative for biocatalytic applications.

Thermophysical properties of glyceline-water mixtures investigated by molecular modelling.

The effect of water content on the static and dynamic properties of the deep eutectic solvent glyceline is studied using molecular dynamics (MD) simulations. Static properties are additionally calculated using the PC-SAFT equation of state. Force fields calibrated on water-free glyceline show predictive power for density and water activity over the entire composition range. In contrast, the PC-SAFT approach using pseudo one-component or two-component modelling strategies performed better for the density or the water activity, respectively. The MD simulations show that at low water content, the hydrogen-bond network between glycerol molecules as well as between glycerol and the cholinium cation is hardly affected by the water molecules while at higher water content, glycerol-glycerol hydrogen bonds are replaced by glycerol-water hydrogen bonds indicating the formation of an aqueous solution accompanied by a strong decrease of the shear viscosity. At the same time, the thermodynamic activity of water increases such that the MD simulations are able to guide the optimal composition with respect to requirements in biocatalytic applications such as low viscosity and low water activity. The combined application of PC-SAFT to efficiently predict static properties and molecular dynamics simulations to predict static and dynamic properties offers a powerful framework in solvent design applications.

Asymmetric Enzymatic Hydration of Unactivated, Aliphatic Alkenes.

The direct enantioselective addition of water to unactivated alkenes could simplify the synthesis of chiral alcohols and solve a long-standing challenge in catalysis. Here we report that an engineered fatty acid hydratase can catalyze the asymmetric hydration of various terminal and internal alkenes. In the presence of a carboxylic acid decoy molecule for activation of the oleate hydratase from E. meningoseptica, asymmetric hydration of unactivated alkenes was achieved with up to 93 % conversion, excellent selectivity (>99 % ee, >95 % regioselectivity), and on a preparative scale.

The ω-transaminase engineering database (oTAED): A navigation tool in protein sequence and structure space.

The ω-Transaminase Engineering Database (oTAED) was established as a publicly accessible resource on sequences and structures of the biotechnologically relevant ω-transaminases (ω-TAs) from Fold types I and IV. The oTAED integrates sequence and structure data, provides a classification based on fold type and sequence similarity, and applies a standard numbering scheme to identify equivalent positions in homologous proteins. The oTAED includes 67 210 proteins (114 655 sequences) which are divided into 169 homologous families based on global sequence similarity. The 44 and 39 highly conserved positions which were identified in Fold type I and IV, respectively, include the known catalytic residues and a large fraction of glycines and prolines in loop regions, which might have a role in protein folding and stability. However, for most of the conserved positions the function is still unknown. Literature information on positions that mediate substrate specificity and stereoselectivity was systematically examined. The standard numbering schemes revealed that many positions which have been described in different enzymes are structurally equivalent. For some positions, multiple functional roles have been suggested based on experimental data in different enzymes. The proposed standard numbering schemes for Fold type I and IV ω-TAs assist with analysis of literature data, facilitate annotation of ω-TAs, support prediction of promising mutation sites, and enable navigation in ω-TA sequence space. Thus, it is a useful tool for enzyme engineering and the selection of novel ω-TA candidates with desired biochemical properties.

Terpene Cyclizations inside a Supramolecular Catalyst: Leaving-Group-Controlled Product Selectivity and Mechanistic Studies.

The tail-to-head terpene cyclization is arguably one of the most complex reactions found in nature. The hydrogen-bond-based resorcinarene capsule represents the first man-made enzyme-like catalyst that is capable of catalyzing this reaction. Based on noncovalent interactions between the capsule and the substrate, the product selectivity can be tuned by using different leaving groups. A detailed mechanistic investigation was performed to elucidate the reaction mechanism. For the cyclization of geranyl acetate, it was found that the cleavage of the leaving group is the rate-determining step. Furthermore, the studies revealed that trace amounts of acid are required as cocatalyst. A series of control experiments demonstrate that a synergistic interplay between the supramolecular capsule and the acid traces is required for catalytic activity.

Self-Assembly Nanostructures of Triglyceride-Water Interfaces Determine Functional Conformations of Candida antarctica Lipase B.

Candida antarctica lipase B (CalB) acts as a lipase when adsorbed to an acylglyceride interface and as an esterase when exposed to an aqueous environment. The effect of the molecular self-assembly nanostructure of triglyceride-water interfaces on structural conformations of adsorbed CalB and the implications to its catalytic function were studied by molecular dynamics simulations. Systems of CalB adsorbed to interfaces and solvated in water were compared. The two environments induced relative motions of helices α5 and α10 that resulted in open and closed conformations. The open conformation was stabilized by interactions between the polar and nonpolar amino acids of α5 and α10 and the nanostructure of triglyceride aggregates, which self-assembled into crystalline-like patterns of alternating polar and nonpolar lamellae. Thus, the structure of CalB has been adapted by evolution to the geometric constraints imposed by the interface nanostructure for optimized catalytic activity. Helices α5 and α10 have two functions. As mobile elements, they ensure access of bulky substrates to the active site in the open conformation. As a part of the active site pocket, they ensure binding of substrate molecules in a productive orientation near the active site. In water, access to the binding site is limited, and the smaller substrate binding site is beneficial for the binding of small, water-soluble substrates. The CalB crystal structure commonly used for protein engineering studies represents an intermediate state between open and closed, and may thus not be adequate to assess the function of CalB, neither as lipase nor as esterase.

Identification of imine reductase-specific sequence motifs.

Chiral amines are valuable building blocks for the production of a variety of pharmaceuticals, agrochemicals and other specialty chemicals. Only recently, imine reductases (IREDs) were discovered which catalyze the stereoselective reduction of imines to chiral amines. Although several IREDs were biochemically characterized in the last few years, knowledge of the reaction mechanism and the molecular basis of substrate specificity and stereoselectivity is limited. To gain further insights into the sequence-function relationships, the Imine Reductase Engineering Database (www.IRED.BioCatNet.de) was established and a systematic analysis of 530 putative IREDs was performed. A standard numbering scheme based on R-IRED-Sk was introduced to facilitate the identification and communication of structurally equivalent positions in different proteins. A conservation analysis revealed a highly conserved cofactor binding region and a predominantly hydrophobic substrate binding cleft. Two IRED-specific motifs were identified, the cofactor binding motif GLGxMGx(5 )[ATS]x(4) Gx(4) [VIL]WNR[TS]x(2) [KR] and the active site motif Gx[DE]x[GDA]x[APS]x(3){K}x[ASL]x[LMVIAG]. Our results indicate a preference toward NADPH for all IREDs and explain why, despite their sequence similarity to ß-hydroxyacid dehydrogenases (ß-HADs), no conversion of ß-hydroxyacids has been observed. Superfamily-specific conservations were investigated to explore the molecular basis of their stereopreference. Based on our analysis and previous experimental results on IRED mutants, an exclusive role of standard position 187 for stereoselectivity is excluded. Alternatively, two standard positions 139 and 194 were identified which are superfamily-specifically conserved and differ in R- and S-selective enzymes.