Folding the human proteome using BioNeMo: A fused dataset of structural models for machine learning purposes.

Human proteins are crucial players in both health and disease. Understanding their molecular landscape is a central topic in biological research. Here, we present an extensive dataset of predicted protein structures for 42,042 distinct human proteins, including splicing variants, derived from the UniProt reference proteome UP000005640. To ensure high quality and comparability, the dataset was generated by combining state-of-the-art modeling-tools AlphaFold 2, OpenFold, and ESMFold, provided within NVIDIA's BioNeMo platform, as well as homology modeling using Innophore's CavitomiX platform. Our dataset is offered in both unedited and edited formats for diverse research requirements. The unedited version contains structures as generated by the different prediction methods, whereas the edited version contains refinements, including a dataset of structures without low prediction-confidence regions and structures in complex with predicted ligands based on homologs in the PDB. We are confident that this dataset represents the most comprehensive collection of human protein structures available today, facilitating diverse applications such as structure-based drug design and the prediction of protein function and interactions.

AI-assisted structural consensus-proteome prediction of human monkeypox viruses isolated within a year after the 2022 multi-country outbreak.

IMPORTANCE: The 2022 outbreak of the monkeypox virus already involves, by April 2023, 110 countries with 86,956 confirmed cases and 119 deaths. Understanding an emerging disease on a molecular level is essential to study infection processes and eventually guide drug discovery at an early stage. To support this, we provide the so far most comprehensive structural proteome of the monkeypox virus, which includes 210 structural models, each computed with three state-of-the-art structure prediction methods. Instead of building on a single-genome sequence, we generated our models from a consensus of 3,713 high-quality genome sequences sampled from patients within 1 year of the outbreak. Therefore, we present an average structural proteome of the currently isolated viruses, including mutational analyses with a special focus on drug-binding sites. Continuing dynamic mutation monitoring within the structural proteome presented here is essential to timely predict possible physiological changes in the evolving virus.

Optimizing variant-specific therapeutic SARS-CoV-2 decoys using deep-learning-guided molecular dynamics simulations.

Treatment of COVID-19 with a soluble version of ACE2 that binds to SARS-CoV-2 virions before they enter host cells is a promising approach, however it needs to be optimized and adapted to emerging viral variants. The computational workflow presented here consists of molecular dynamics simulations for spike RBD-hACE2 binding affinity assessments of multiple spike RBD/hACE2 variants and a novel convolutional neural network architecture working on pairs of voxelized force-fields for efficient search-space reduction. We identified hACE2-Fc K31W and multi-mutation variants as high-affinity candidates, which we validated in vitro with virus neutralization assays. We evaluated binding affinities of these ACE2 variants with the RBDs of Omicron BA.3, Omicron BA.4/BA.5, and Omicron BA.2.75 in silico. In addition, candidates produced in Nicotiana benthamiana, an expression organism for potential large-scale production, showed a 4.6-fold reduction in half-maximal inhibitory concentration (IC50) compared with the same variant produced in CHO cells and an almost six-fold IC50 reduction compared with wild-type hACE2-Fc.

Structural bioinformatics analysis of SARS-CoV-2 variants reveals higher hACE2 receptor binding affinity for Omicron B.1.1.529 spike RBD compared to wild type reference.

To date, more than 263 million people have been infected with SARS-CoV-2 during the COVID-19 pandemic. In many countries, the global spread occurred in multiple pandemic waves characterized by the emergence of new SARS-CoV-2 variants. Here we report a sequence and structural-bioinformatics analysis to estimate the effects of amino acid substitutions on the affinity of the SARS-CoV-2 spike receptor binding domain (RBD) to the human receptor hACE2. This is done through qualitative electrostatics and hydrophobicity analysis as well as molecular dynamics simulations used to develop a high-precision empirical scoring function (ESF) closely related to the linear interaction energy method and calibrated on a large set of experimental binding energies. For the latest variant of concern (VOC), B.1.1.529 Omicron, our Halo difference point cloud studies reveal the largest impact on the RBD binding interface compared to all other VOC. Moreover, according to our ESF model, Omicron achieves a much higher ACE2 binding affinity than the wild type and, in particular, the highest among all VOCs except Alpha and thus requires special attention and monitoring.

Recent changes in the mutational dynamics of the SARS-CoV-2 main protease substantiate the danger of emerging resistance to antiviral drugs.

Introduction: The current coronavirus pandemic is being combated worldwide by nontherapeutic measures and massive vaccination programs. Nevertheless, therapeutic options such as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main-protease (Mpro) inhibitors are essential due to the ongoing evolution toward escape from natural or induced immunity. While antiviral strategies are vulnerable to the effects of viral mutation, the relatively conserved Mpro makes an attractive drug target: Nirmatrelvir, an antiviral targeting its active site, has been authorized for conditional or emergency use in several countries since December 2021, and a number of other inhibitors are under clinical evaluation. We analyzed recent SARS-CoV-2 genomic data, since early detection of potential resistances supports a timely counteraction in drug development and deployment, and discovered accelerated mutational dynamics of Mpro since early December 2021. Methods: We performed a comparative analysis of 10.5 million SARS-CoV-2 genome sequences available by June 2022 at GISAID to the NCBI reference genome sequence NC_045512.2. Amino-acid exchanges within high-quality regions in 69,878 unique Mpro sequences were identified and time- and in-depth sequence analyses including a structural representation of mutational dynamics were performed using in-house software. Results: The analysis showed a significant recent event of mutational dynamics in Mpro. We report a remarkable increase in mutational variability in an eight-residue long consecutive region (R188-G195) near the active site since December 2021. Discussion: The increased mutational variability in close proximity to an antiviral-drug binding site as described herein may suggest the onset of the development of antiviral resistance. This emerging diversity urgently needs to be further monitored and considered in ongoing drug development and lead optimization.

Serine 477 plays a crucial role in the interaction of the SARS-CoV-2 spike protein with the human receptor ACE2.

Since the worldwide outbreak of the infectious disease COVID-19, several studies have been published to understand the structural mechanism of the novel coronavirus SARS-CoV-2. During the infection process, the SARS-CoV-2 spike (S) protein plays a crucial role in the receptor recognition and cell membrane fusion process by interacting with the human angiotensin-converting enzyme 2 (hACE2) receptor. However, new variants of these spike proteins emerge as the virus passes through the disease reservoir. This poses a major challenge for designing a potent antigen for an effective immune response against the spike protein. Through a normal mode analysis (NMA) we identified the highly flexible region in the receptor binding domain (RBD) of SARS-CoV-2, starting from residue 475 up to residue 485. Structurally, the position S477 shows the highest flexibility among them. At the same time, S477 is hitherto the most frequently exchanged amino acid residue in the RBDs of SARS-CoV-2 mutants. Therefore, using MD simulations, we have investigated the role of S477 and its two frequent mutations (S477G and S477N) at the RBD during the binding to hACE2. We found that the amino acid exchanges S477G and S477N strengthen the binding of the SARS-COV-2 spike with the hACE2 receptor.

A Fungal Ascorbate Oxidase with Unexpected Laccase Activity.

Ascorbate oxidases are an enzyme group that has not been explored to a large extent. So far, mainly ascorbate oxidases from plants and only a few from fungi have been described. Although ascorbate oxidases belong to the well-studied enzyme family of multi-copper oxidases, their function is still unclear. In this study, Af_AO1, an enzyme from the fungus Aspergillus flavus, was characterized. Sequence analyses and copper content determination demonstrated Af_AO1 to belong to the multi-copper oxidase family. Biochemical characterization and 3D-modeling revealed a similarity to ascorbate oxidases, but also to laccases. Af_AO1 had a 10-fold higher affinity to ascorbic acid (KM = 0.16 ± 0.03 mM) than to ABTS (KM = 1.89 ± 0.12 mM). Furthermore, the best fitting 3D-model was based on the ascorbate oxidase from Cucurbita pepo var. melopepo. The laccase-like activity of Af_AO1 on ABTS (Vmax = 11.56 ± 0.15 µM/min/mg) was, however, not negligible. On the other hand, other typical laccase substrates, such as syringaldezine and guaiacol, were not oxidized by Af_AO1. According to the biochemical and structural characterization, Af_AO1 was classified as ascorbate oxidase with unusual, laccase-like activity.

Evolving the Promiscuity of Elizabethkingia meningoseptica Oleate Hydratase for the Regio- and Stereoselective Hydration of Oleic Acid Derivatives.

The addition of water to non-activated carbon-carbon double bonds catalyzed by fatty acid hydratases (FAHYs) allows for highly regio- and stereoselective oxyfunctionalization of renewable oil feedstock. So far, the applicability of FAHYs has been limited to free fatty acids, mainly owing to the requirement of a carboxylate function for substrate recognition and binding. Herein, we describe for the first time the hydration of oleic acid (OA) derivatives lacking this free carboxylate by the oleate hydratase from Elizabethkingia meningoseptica (OhyA). Molecular docking of OA to the OhyA 3D-structure and a sequence alignment uncovered conserved amino acid residues at the entrance of the substrate channel as target positions for enzyme engineering. Exchange of selected amino acids gave rise to OhyA variants which showed up to an 18-fold improved conversion of OA derivatives, while retaining the excellent regio- and stereoselectivity in the olefin hydration reaction.

Identification of Key Residues for Enzymatic Carboxylate Reduction.

Carboxylate reductases (CARs, E.C. 1.2.1.30) generate aldehydes from their corresponding carboxylic acid with high selectivity. Little is known about the structure of CARs and their catalytically important amino acid residues. The identification of key residues for carboxylate reduction provides a starting point to gain deeper understanding of enzymatic carboxylate reduction. A multiple sequence alignment of CARs with confirmed activity recently identified in our lab and from the literature revealed a fingerprint of conserved amino acids. We studied the function of conserved residues by multiple sequence alignments and mutational replacements of these residues. In this study, single-site alanine variants of Neurospora crassa CAR were investigated to determine the contribution of conserved residues to the function, expressability or stability of the enzyme. The effect of amino acid replacements was investigated by analyzing enzymatic activity of the variants in vivo and in vitro. Supported by molecular modeling, we interpreted that five of these residues are essential for catalytic activity, or substrate and co-substrate binding. We identified amino acid residues having significant impact on CAR activity. Replacement of His 237, Glu 433, Ser 595, Tyr 844, and Lys 848 by Ala abolish CAR activity, indicating their key role in acid reduction. These results may assist in the functional annotation of CAR coding genes in genomic databases. While some other conserved residues decreased activity or had no significant impact, four residues increased the specific activity of NcCAR variants when replaced by alanine. Finally, we showed that NcCAR wild-type and mutants efficiently reduce aliphatic acids.

Regioselective para-Carboxylation of Catechols with a Prenylated Flavin Dependent Decarboxylase.

The utilization of CO2 as a carbon source for organic synthesis meets the urgent demand for more sustainability in the production of chemicals. Herein, we report on the enzyme-catalyzed para-carboxylation of catechols, employing 3,4-dihydroxybenzoic acid decarboxylases (AroY) that belong to the UbiD enzyme family. Crystal structures and accompanying solution data confirmed that AroY utilizes the recently discovered prenylated FMN (prFMN) cofactor, and requires oxidative maturation to form the catalytically competent prFMNiminium species. This study reports on the in vitro reconstitution and activation of a prFMN-dependent enzyme that is capable of directly carboxylating aromatic catechol substrates under ambient conditions. A reaction mechanism for the reversible decarboxylation involving an intermediate with a single covalent bond between a quinoid adduct and cofactor is proposed, which is distinct from the mechanism of prFMN-associated 1,3-dipolar cycloadditions in related enzymes.

Exploring the Catalytic Promiscuity of Phenolic Acid Decarboxylases: Asymmetric, 1,6-Conjugate Addition of Nucleophiles Across 4-Hydroxystyrene.

The catalytic promiscuity of a ferulic acid decarboxylase from Enterobacter sp. (FDC_Es) and phenolic acid decarboxylases (PADs) for the asymmetric conjugate addition of water across the C=C bond of hydroxystyrenes was extended to the N-, C- and S-nucleophiles methoxyamine, cyanide and propanethiol to furnish the corresponding addition products in up to 91% ee. The products obtained from the biotransformation employing the most suitable enzyme/nucleophile pairs were isolated and characterized after optimizing the reaction conditions. Finally, a mechanistic rationale supported by quantum mechanical calculations for the highly (S)-selective addition of cyanide is proposed.

Small cause, large effect: Structural characterization of cutinases from Thermobifida cellulosilytica.

We have investigated the structures of two native cutinases from Thermobifida cellulosilytica, namely Thc_Cut1 and Thc_Cut2 as well as of two variants, Thc_Cut2_DM (Thc_Cut2_ Arg29Asn_Ala30Val) and Thc_Cut2_TM (Thc_Cut2_Arg19Ser_Arg29Asn_Ala30Val). The four enzymes showed different activities towards the aliphatic polyester poly(lactic acid) (PLLA). The crystal structures of the four enzymes were successfully solved and in combination with Small Angle X-Ray Scattering (SAXS) the structural features responsible for the selectivity difference were elucidated. Analysis of the crystal structures did not indicate significant conformational differences among the different cutinases. However, the distinctive SAXS scattering data collected from the enzymes in solution indicated a remarkable surface charge difference. The difference in the electrostatic and hydrophobic surface properties could explain potential alternative binding modes of the four cutinases on PLLA explaining their distinct activities. Biotechnol. Bioeng. 2017;114: 2481-2488. © 2017 Wiley Periodicals, Inc.

Polyester hydrolysis is enhanced by a truncated esterase: Less is more.

An esterase from Clostridium botulinum (Cbotu_EstA) previously reported to hydrolyze the biodegradable polyester poly(butylene adipate-co-terephthalate) was redesigned to improve the hydrolysis of synthetic polyesters. Increased activity was indeed observed for del71Cbotu_EstA variant, which performed activity on the widespread polyester polyethylene terephthalate, which was not able to be attacked by the wild-type enzyme Cbotu_EstA. Analysis of the 3D structure of the enzyme showed that removing 71 residues at the N-terminus of the enzyme exposed a hydrophobic patch on the surface and improved sorption of hydrophobic polyesters concomitantly facilitating the access of the polymer to the active site. These results show a new route for enhancing enzyme activity for hydrolysis and modification of polyesters.

PpEst is a novel PBAT degrading polyesterase identified by proteomic screening of Pseudomonas pseudoalcaligenes.

A novel esterase, PpEst, that hydrolyses the co-aromatic-aliphatic polyester poly(1,4-butylene adipate-co-terephthalate) (PBAT) was identified by proteomic screening of the Pseudomonas pseudoalcaligenes secretome. PpEst was induced by the presence of PBAT in the growth media and had predicted arylesterase (EC 3.1.1.2) activity. PpEst showed polyesterase activity on both whole and milled PBAT film releasing terephthalic acid and 4-(4-hydroxybutoxycarbonyl)benzoic acid while end product inhibition by 4-(4-hydroxybutoxycarbonyl)benzoic acid was observed. Modelling of an aromatic polyester mimicking oligomer into the PpEst active site indicated that the binding pocket could be big enough to accommodate large polymers. This is the first report of a PBAT degrading enzyme being identified by proteomic screening and shows that this approach can contribute to the discovery of new polymer hydrolysing enzymes. Moreover, these results indicate that arylesterases could be an interesting enzyme class for identifications of polyesterases.

Structures of almond hydroxynitrile lyase isoenzyme 5 provide a rationale for the lack of oxidoreductase activity in flavin dependent HNLs.

Hydroxynitrile lyases (HNLs) catalyze the asymmetric addition of HCN to aldehydes producing enantiomerically pure cyanohydrins. These enzymes can be heterologously expressed in large quantities making them interesting candidates for industrial applications. The HNLs from Rosaceae evolved from flavin dependent dehydrogenase/oxidase structures. Here we report the high resolution X-ray structure of the highly glycosylated Prunus amygdalus HNL isoenzyme5 (PaHNL5 V317A) expressed in Aspergillus niger and its complex with benzyl alcohol. A comparison with the structure of isoenzyme PaHNL1 indicates a higher accessibility to the active site and a larger cavity for PaHNL5. Additionally, the PaHNL5 complex structure with benzyl alcohol was compared with the structurally related aryl-alcohol oxidase (AAO). Even though both enzymes contain an FAD-cofactor and histidine residues at crucial positions in the active site, PaHNL5 lacks the oxidoreductase activity. The structures indicate that in PaHNLs benzyl alcohol is bound too far away from the FAD cofactor in order to be oxidized.

An Esterase from Anaerobic Clostridium hathewayi Can Hydrolyze Aliphatic-Aromatic Polyesters.

Recently, a variety of biodegradable polymers have been developed as alternatives to recalcitrant materials. Although many studies on polyester biodegradability have focused on aerobic environments, there is much less known on biodegradation of polyesters in natural and artificial anaerobic habitats. Consequently, the potential of anaerobic biogas sludge to hydrolyze the synthetic compostable polyester PBAT (poly(butylene adipate-co-butylene terephthalate) was evaluated in this study. On the basis of reverse-phase high-performance liquid chromatography (RP-HPLC) analysis, accumulation of terephthalic acid (Ta) was observed in all anaerobic batches within the first 14 days. Thereafter, a decline of Ta was observed, which occurred presumably due to consumption by the microbial population. The esterase Chath_Est1 from the anaerobic risk 1 strain Clostridium hathewayi DSM-13479 was found to hydrolyze PBAT. Detailed characterization of this esterase including elucidation of the crystal structure was performed. The crystal structure indicates that Chath_Est1 belongs to the α/ß-hydrolases family. This study gives a clear hint that also micro-organisms in anaerobic habitats can degrade manmade PBAT.

Hydrolysis of synthetic polyesters by Clostridium botulinum esterases.

Two novel esterases from the anaerobe Clostridium botulinum ATCC 3502 (Cbotu_EstA and Cbotu_EstB) were expressed in Escherichia coli BL21-Gold(DE3) and were found to hydrolyze the polyester poly(butylene adipate-co-butylene terephthalate) (PBAT). The active site residues (triad Ser, Asp, His) are present in both enzymes at the same location only with some amino acid variations near the active site at the surrounding of aspartate. Yet, Cbotu_EstA showed higher kcat values on para-nitrophenyl butyrate and para-nitrophenyl acetate and was considerably more active (sixfold) on PBAT. The entrance to the active site of the modeled Cbotu_EstB appears more narrowed compared to the crystal structure of Cbotu_EstA and the N-terminus is shorter which could explain its lower activity on PBAT. The Cbotu_EstA crystal structure consists of two regions that may act as movable cap domains and a zinc metal binding site.

Regioselective Enzymatic ß-Carboxylation of para-Hydroxy- styrene Derivatives Catalyzed by Phenolic Acid Decarboxylases.

We report on a 'green' method for the utilization of carbon dioxide as C1 unit for the regioselective synthesis of (E)-cinnamic acids via regioselective enzymatic carboxylation of para-hydroxystyrenes. Phenolic acid decarboxylases from bacterial sources catalyzed the ß-carboxylation of para-hydroxystyrene derivatives with excellent regio- and (E/Z)-stereoselectivity by exclusively acting at the ß-carbon atom of the C=C side chain to furnish the corresponding (E)-cinnamic acid derivatives in up to 40% conversion at the expense of bicarbonate as carbon dioxide source. Studies on the substrate scope of this strategy are presented and a catalytic mechanism is proposed based on molecular modelling studies supported by mutagenesis of amino acid residues in the active site.

Structure-Based Mechanism of Oleate Hydratase from Elizabethkingia meningoseptica.

Hydratases provide access to secondary and tertiary alcohols by regio- and/or stereospecifically adding water to carbon-carbon double bonds. Thereby, hydroxy groups are introduced without the need for costly cofactor recycling, and that makes this approach highly interesting on an industrial scale. Here we present the first crystal structure of a recombinant oleate hydratase originating from Elizabethkingia meningoseptica in the presence of flavin adenine dinucleotide (FAD). A structure-based mutagenesis study targeting active site residues identified E122 and Y241 as crucial for the activation of a water molecule and for protonation of the double bond, respectively. Moreover, we also observed that two-electron reduction of FAD results in a sevenfold increase in the substrate hydration rate. We propose the first reaction mechanism for this enzyme class that explains the requirement for the flavin cofactor and the involvement of conserved amino acid residues in this regio- and stereoselective hydration.

The crystal structure of D-threonine aldolase from Alcaligenes xylosoxidans provides insight into a metal ion assisted PLP-dependent mechanism.

Threonine aldolases catalyze the pyridoxal phosphate (PLP) dependent cleavage of threonine into glycine and acetaldehyde and play a major role in the degradation of this amino acid. In nature, L- as well as D-specific enzymes have been identified, but the exact physiological function of D-threonine aldolases (DTAs) is still largely unknown. Both types of enantio-complementary enzymes have a considerable potential in biocatalysis for the stereospecific synthesis of various ß-hydroxy amino acids, which are valuable building blocks for the production of pharmaceuticals. While several structures of L-threonine aldolases (LTAs) have already been determined, no structure of a DTA is available to date. Here, we report on the determination of the crystal structure of the DTA from Alcaligenes xylosoxidans (AxDTA) at 1.5 Å resolution. Our results underline the close relationship of DTAs and alanine racemases and allow the identification of a metal binding site close to the PLP-cofactor in the active site of the enzyme which is consistent with the previous observation that divalent cations are essential for DTA activity. Modeling of AxDTA substrate complexes provides a rationale for this metal dependence and indicates that binding of the ß-hydroxy group of the substrate to the metal ion very likely activates this group and facilitates its deprotonation by His193. An equivalent involvement of a metal ion has been implicated in the mechanism of a serine dehydratase, which harbors a metal ion binding site in the vicinity of the PLP cofactor at the same position as in DTA. The structure of AxDTA is completely different to available structures of LTAs. The enantio-complementarity of DTAs and LTAs can be explained by an approximate mirror symmetry of crucial active site residues relative to the PLP-cofactor.