Caver Web 1.0: identification of tunnels and channels in proteins and analysis of ligand transport

Comparison of two peroxidases with high potential for biotechnology applications - HRP vs. APEX2

Peroxidases are essential elements in many biotechnological applications. An especially interesting concept involves split enzymes, where the enzyme is separated into two smaller and inactive proteins that can dimerize into a fully active enzyme. Such split forms were developed for the horseradish peroxidase (HRP) and ascorbate peroxidase (APX) already. Both peroxidases have a high potential for biotechnology applications. In the present study, we performed biophysical comparisons of these two peroxidases and their split analogues. The active site availability is similar for all four structures. The split enzymes are comparable in stability with their native analogues, meaning that they can be used for further biotechnology applications. Also, the tertiary structures of the two peroxidases are similar. However, differences that might help in choosing one system over another for biotechnology applications were noticed. The main difference between the two systems is glycosylation which is not present in the case of APX/sAPEX2, while it has a high impact on the HRP/sHRP stability. Further differences are calcium ions and cysteine bridges that are present only in the case of HRP/sHRP. Finally, computational results identified sAPEX2 as the systems with the smallest structural variations during molecular dynamics simulations showing its dominant stability comparing to other simulated proteins. Taken all together, the sAPEX2 system has a high potential for biotechnological applications due to the lack of glycans and cysteines, as well as due to high stability.

Discovery and biocatalytic characterization of opine dehydrogenases by metagenome mining

Enzymatic processes play an increasing role in synthetic organic chemistry which requires the access to a broad and diverse set of enzymes. Metagenome mining is a valuable and efficient way to discover novel enzymes with unique properties for biotechnological applications. Here, we report the discovery and biocatalytic characterization of six novel metagenomic opine dehydrogenases from a hot spring environment (mODHs) (EC 1.5.1.X). These enzymes catalyze the asymmetric reductive amination between an amino acid and a keto acid resulting in opines which have defined biochemical roles and represent promising building blocks for pharmaceutical applications. The newly identified enzymes exhibit unique substrate specificity and higher thermostability compared to known examples. The feature that they preferably utilize negatively charged polar amino acids is so far unprecedented for opine dehydrogenases. We have identified two spatially correlated positions in their active sites that govern this substrate specificity and demonstrated a switch of substrate preference by site-directed mutagenesis. While they still suffer from a relatively narrow substrate scope, their enhanced thermostability and the orthogonality of their substrate preference make them a valuable addition to the toolbox of enzymes for reductive aminations. Importantly, enzymatic reductive aminations with highly polar amines are very rare in the literature. Thus, the preparative-scale enzymatic production, purification, and characterization of three highly functionalized chiral secondary amines lend a special significance to our work in filling this gap. KEY POINTS: • Six new opine dehydrogenases have been discovered from a hot spring metagenome • The newly identified enzymes display a unique substrate scope • Substrate specificity is governed by two correlated active-site residues.

Unlocking the potential of enzyme engineering via rational computational design strategies

Enzymes play a pivotal role in various industries by enabling efficient, eco-friendly, and sustainable chemical processes. However, the low turnover rates and poor substrate selectivity of enzymes limit their large-scale applications. Rational computational enzyme design, facilitated by computational algorithms, offers a more targeted and less labor-intensive approach. There has been notable advancement in employing rational computational protein engineering strategies to overcome these issues, it has not been comprehensively reviewed so far. This article reviews recent developments in rational computational enzyme design, categorizing them into three types: structure-based, sequence-based, and data-driven machine learning computational design. Case studies are presented to demonstrate successful enhancements in catalytic activity, stability, and substrate selectivity. Lastly, the article provides a thorough analysis of these approaches, highlights existing challenges and potential solutions, and offers insights into future development directions.

A contribution to lipid digestion of Odobenidae family: Computational analysis of gastric and pancreatic lipases from walruses (Odobenus rosmarus divergens)

Triacylglycerols (TAGs) are a primary energy source for marine mammals during lipid digestion. Walruses (Odobenus rosmarus divergens) consume prey with a high content of long-chain polyunsaturated fatty acids; however, their digestive physiology and lipid digestion remain poorly studied. The present study aims to model and characterize the gastric (PWGL) and pancreatic (PWPL) lipases of Pacific walruses using an in-silico approach. The confident 3D models of PWGL and PWPL were obtained via homology modeling and protein threading and displayed the structural features of lipases. Molecular docking analysis demonstrated substrate selectivity for long-chain TAG (Trieicosapentaenoin; TC20:5n-3) in PWGL and short-chain TAG (Trioctanoin; TC8:0) in PWPL. Molecular dynamics simulations demonstrate that PWGL maintains structural stability at salinity conditions, with no significant conformational changes observed. In the simulations of PWGL bound to tridocosahexaenoin (TC22:6n-3), the protein is considerably stable at all three salinity conditions, but fluctuations are observed in the regions associated with catalytic sites and the lid, indicating the hydrolysis of the substrate. This is the first study to report on the digestion of TAGs in walruses, including modeling and lipases characterization and proposing a digestive tract for pinnipeds.

Imidazole[1,5-a]pyridine derivatives as EGFR tyrosine kinase inhibitors unraveled by umbrella sampling and steered molecular dynamics simulations

Although the use of the tyrosine kinase inhibitors (TKIs) has been proved that it can save live in a cancer treatment, the currently used drugs bring in many undesirable side-effects. Therefore, the search for new drugs and an evaluation of their efficiency are intensively carried out. Recently, a series of eighteen imidazole[1,5-a]pyridine derivatives were synthetized by us, and preliminary analyses pointed out their potential to be an important platform for pharmaceutical development owing to their promising actions as anticancer agents and enzyme (kinase, HIV-protease,…) inhibitors. In the present theoretical study, we further analyzed their efficiency in using a realistic scenario of computational drug design. Our protocol has been developed to not only observe the atomistic interaction between the EGFR protein and our 18 novel compounds using both umbrella sampling and steered molecular dynamics simulations, but also determine their absolute binding free energies. Calculated properties of the 18 novel compounds were in detail compared with those of two known drugs, erlotinib and osimertinib, currently used in cancer treatment. Inspiringly the simulation results promote three imidazole[1,5-a]pyridine derivatives as promising inhibitors into a further step of clinical trials.

Treatment of flexibility of protein backbone in simulations of protein-ligand interactions using steered molecular dynamics

To ensure that an external force can break the interaction between a protein and a ligand, the steered molecular dynamics simulation requires a harmonic restrained potential applied to the protein backbone. A usual practice is that all or a certain number of protein's heavy atoms or Cα atoms are fixed, being restrained by a small force. This present study reveals that while fixing both either all heavy atoms and or all Cα atoms is not a good approach, while fixing a too small number of few atoms sometimes cannot prevent the protein from rotating under the influence of the bulk water layer, and the pulled molecule may smack into the wall of the active site. We found that restraining the Cα atoms under certain conditions is more relevant. Thus, we would propose an alternative solution in which only the Cα atoms of the protein at a distance larger than 1.2 nm from the ligand are restrained. A more flexible, but not too flexible, protein will be expected to lead to a more natural release of the ligand.

Functional and structural characterization of an IclR family transcription factor for the development of dicarboxylic acid biosensors

Prokaryotic transcription factors (TFs) regulate gene expression in response to small molecules, thus representing promising candidates as versatile small molecule-detecting biosensors valuable for synthetic biology applications. The engineering of such biosensors requires thorough in vitro and in vivo characterization of TF ligand response as well as detailed molecular structure information. In this work, we functionally and structurally characterize the Pca regulon regulatory protein (PcaR) transcription factor belonging to the IclR transcription factor family. Here, we present in vitro functional analysis of the ligand profile of PcaR and the construction of genetic circuits for the characterization of PcaR as an in vivo biosensor in the model eukaryote Saccharomyces cerevisiae. We report the crystal structures of PcaR in the apo state and in complex with one of its ligands, succinate, which suggests the mechanism of dicarboxylic acid recognition by this transcription factor. This work contributes key structural and functional insights enabling the engineering of PcaR for dicarboxylic acid biosensors, in addition to providing more insights into the IclR family of regulators.

Structure of adenylyl cyclase 5 in complex with Gβγ offers insights into ADCY5-related dyskinesia

The nine different membrane-anchored adenylyl cyclase isoforms (AC1-9) in mammals are stimulated by the heterotrimeric G protein, Gαs, but their response to Gβγ regulation is isoform specific. In the present study, we report cryo-electron microscope structures of ligand-free AC5 in complex with Gβγ and a dimeric form of AC5 that could be involved in its regulation. Gβγ binds to a coiled-coil domain that links the AC transmembrane region to its catalytic core as well as to a region (C1b) that is known to be a hub for isoform-specific regulation. We confirmed the Gβγ interaction with both purified proteins and cell-based assays. Gain-of-function mutations in AC5 associated with human familial dyskinesia are located at the interface of AC5 with Gβγ and show reduced conditional activation by Gβγ, emphasizing the importance of the observed interaction for motor function in humans. We propose a molecular mechanism wherein Gβγ either prevents dimerization of AC5 or allosterically modulates the coiled-coil domain, and hence the catalytic core. As our mechanistic understanding of how individual AC isoforms are uniquely regulated is limited, studies such as this may provide new avenues for isoform-specific drug development.

Opportunities and challenges in design and optimization of protein function

The field of protein design has made remarkable progress over the past decade. Historically, the low reliability of purely structure-based design methods limited their application, but recent strategies that combine structure-based and sequence-based calculations, as well as machine learning tools, have dramatically improved protein engineering and design. In this Review, we discuss how these methods have enabled the design of increasingly complex structures and therapeutically relevant activities. Additionally, protein optimization methods have improved the stability and activity of complex eukaryotic proteins. Thanks to their increased reliability, computational design methods have been applied to improve therapeutics and enzymes for green chemistry and have generated vaccine antigens, antivirals and drug-delivery nano-vehicles. Moreover, the high success of design methods reflects an increased understanding of basic rules that govern the relationships among protein sequence, structure and function. However, de novo design is still limited mostly to α-helix bundles, restricting its potential to generate sophisticated enzymes and diverse protein and small-molecule binders. Designing complex protein structures is a challenging but necessary next step if we are to realize our objective of generating new-to-nature activities.

Coupling and regulation mechanisms of the flavin-dependent halogenase PyrH observed by infrared difference spectroscopy

Flavin-dependent halogenases are central enzymes in the production of halogenated secondary metabolites in various organisms and they constitute highly promising biocatalysts for regioselective halogenation. The mechanism of these monooxygenases includes formation of hypohalous acid from a reaction of fully reduced flavin with oxygen and halide. The hypohalous acid then diffuses via a tunnel to the substrate-binding site for halogenation of tryptophan and other substrates. Oxidized flavin needs to be reduced for regeneration of the enzyme, which can be performed in vitro by a photoreduction with blue light. Here, we employed this photoreduction to study characteristic structural changes associated with the transition from oxidized to fully reduced flavin in PyrH from Streptomyces rugosporus as a model for tryptophan-5-halogenases. The effect of the presence of bromide and chloride or the absence of any halides on the UV-vis spectrum of the enzyme demonstrated a halide-dependent structure of the flavin-binding pocket. Light-induced FTIR difference spectroscopy was applied and the signals assigned by selective isotope labeling of the protein moiety. The identified structural changes in α-helix and β-sheet elements were strongly dependent on the presence of bromide, chloride, the substrate tryptophan, and the product 5-chloro-tryptophan, respectively. We identified a clear allosteric coupling in solution at ambient conditions between cofactor-binding site and substrate-binding site that is active in both directions, despite their separation by a tunnel. We suggest that this coupling constitutes a fine-tuned mechanism for the promotion of the enzymatic reaction of flavin-dependent halogenases in dependence of halide and substrate availability.

Structure-Guided Evolution Modulate Alcohol Oxidase to Improve Ethanol Oxidation Performance

A high ethanol usage of alcohol oxidase (AOX) was required in industry. In this study, a "expand substrate pocket" strategy achieved a high activity AOX from Hansenula polymorpha (H. polymorpha) by Phe to Val residue (F/V) site-directed mutation to enlarge ethanol channel. Although H. Polymorpha AOX (HpAOX) possessed respectively 71.3% and 76.1% similarity with AOX (PpAOX) from Pichia pastoris (P. pastoris) in DNA and protein sequences, their active site structures including catalytic site and substrate channel were similar according to computer-aided analysis. After 3D structure analysis, Phe99 residue of their substrate channels was the most important residue to impact enzyme activity because of its large aromatic side chains. F99V mutation of HpAOX (HpAOXF99V) was designed and executed based on the enzyme catalytic mechanism and molecular computation in order to allow more larger size ethanol into active site. The highest enzyme activity of the fourth strains of HpAOXF99V mutant strain exhibited 12.06-folds increase than that of the host GS115 strain. Furthermore, kinetic studies indicated that the HpAOXF99V significantly promoted catalytic efficiency of ethanol than HpAOX, including Km, Vmax, kcat and kcat/Km. We also provided a new insight that the cofactor FAD irritated both active AOX octamer biosynthesis production and enzyme-catalysed ability due to help enzyme assembly and redox potential.

Ionic liquid modification reshapes the substrate pockets of lipase to boost its stability and activity in vitamin E succinate synthesis

BACKGROUND

The relative low stability, reusability and activity of enzymes made the industrial production of vitamin E succinate (VES) can only be performed with complex processes and high cost using chemical methods. To address these issues, in the present study, an ionic liquids (ILs) modification strategy was developed to improve the activity and stability of lipases in VES synthesis.

RESULTS

The results showed that the [1-butyl-3-methyl imidazole] [N-acetyl-l-proline] ILs modified Candida rugosa lipase (CRL) has the highest modification degree (48.28%), activity (774 U g-1 ), thermostability and solvent tolerance in three selected modifiers. Additionally, after reaction condition optimization, the highest yield of VES can be improved to 95.18% at 45 °C for 15 h, which was significantly improved compared to some previous studies.

CONCLUSION

In the present study, a high-efficiency VES synthesis strategy was successfully developed via modification of lipase. Moreover, the mechanism by which ILs modification can enhance the activity and stability of lipase was investigated via both experimental and computational-aided methods. Molecular dynamics simulation suggested that ILs modification changed the geometry of Phe344 from flat to upright, which significantly reshaped and enhanced the size of substrate binding pocket of CRL. It is also agreement with our circular dichroism and fluorescence spectroscopy results, which suggested that the modification changed the secondary structure of CRL to a certain extent. The larger pocket also endowed the suitable binding pose of succinate, which made the hydrogen bonds between succinate and active site Ser209 become stronger, and thus improving the yield of VES. © 2023 Society of Chemical Industry.

Achieving cinnamic acid amides in water by a variant of acyltransferase from Mycobacterium smegmatis and its immobilized form using Ni-NTA modified aspen powder as a carrier

An aqueous N-acylation reaction for preparing cinnamic acid amides was realized by using a variant of acyltransferase from Mycobacterium smegmatis (MsAcT-L12A), whereas the wild-type MsAcT showed no activity. MsAcT-L12A exhibited broad substrate adaptability, and preferred the substrates with electron-donating group. When the vinyl cinnamate (1a, 40 mM) and p-methoxyaniline (2a, 4 mM) were involved in the reaction, the excellent yield reached to 86.7 % ± 2.1 % within 3 h by MsAcT-L12A (1 mgpro./mL) in a PBS buffer (100 mM, pH 8.0) at 25 °C. The aqueous N-acylation reaction could be further improved by using an immobilized MsAcT-L12A. The biomass aspen powder (AP) as a carrier provided a low-cost, green, and environmental-friendly immobilization strategy. After it was modified by Ni-NTA, the obtained Ni-NAP could realize one-step purification and immobilization of MsAcT-L12A. The accomplished MsAcT-L12A-Ni-NAP exhibited excellent stability and recyclability, and retained its relative yield as 83.3 % ± 2.2 % even after the 7th cycle of reuse. Using only PBS buffer as a reaction medium, the operation for MsAcT-L12A-catalyzed acyl transfer was greatly simplified, and the improved stabilities of MsAcT-L12A-Ni-NAP could enhance its application potential.

Dirigent isoflavene-forming PsPTS2: 3D structure, stereochemical, and kinetic characterization comparison with pterocarpan-forming PsPTS1 homolog in pea

Pea phytoalexins (-)-maackiain and (+)-pisatin have opposite C6a/C11a configurations, but biosynthetically how this occurs is unknown. Pea dirigent-protein (DP) PsPTS2 generates 7,2'-dihydroxy-4',5'-methylenedioxyisoflav-3-ene (DMDIF), and stereoselectivity toward four possible 7,2'-dihydroxy-4',5'-methylenedioxyisoflavan-4-ol (DMDI) stereoisomers was investigated. Stereoisomer configurations were determined using NMR spectroscopy, electronic circular dichroism, and molecular orbital analyses. PsPTS2 efficiently converted cis-(3R,4R)-DMDI into DMDIF 20-fold faster than the trans-(3R,4S)-isomer. The 4R-configured substrate's near β-axial OH orientation significantly enhanced its leaving group abilities in generating A-ring mono-quinone methide (QM), whereas 4S-isomer's α-equatorial-OH was a poorer leaving group. Docking simulations indicated that the 4R-configured β-axial OH was closest to Asp51, whereas 4S-isomer's α-equatorial OH was further away. Neither cis-(3S,4S)- nor trans-(3S,4R)-DMDIs were substrates, even with the former having C3/C4 stereochemistry as in (+)-pisatin. PsPTS2 used cis-(3R,4R)-7,2'-dihydroxy-4'-methoxyisoflavan-4-ol [cis-(3R,4R)-DMI] and C3/C4 stereoisomers to give 2',7-dihydroxy-4'-methoxyisoflav-3-ene (DMIF). DP homologs may exist in licorice (Glycyrrhiza pallidiflora) and tree legume Bolusanthus speciosus, as DMIF occurs in both species. PsPTS1 utilized cis-(3R,4R)-DMDI to give (-)-maackiain 2200-fold more efficiently than with cis-(3R,4R)-DMI to give (-)-medicarpin. PsPTS1 also slowly converted trans-(3S,4R)-DMDI into (+)-maackiain, reflecting the better 4R configured OH leaving group. PsPTS2 and PsPTS1 provisionally provide the means to enable differing C6a and C11a configurations in (+)-pisatin and (-)-maackiain, via identical DP-engendered mono-QM bound intermediate generation, which PsPTS2 either re-aromatizes to give DMDIF or PsPTS1 intramolecularly cyclizes to afford (-)-maackiain. Substrate docking simulations using PsPTS2 and PsPTS1 indicate cis-(3R,4R)-DMDI binds in the anti-configuration in PsPTS2 to afford DMDIF, and the syn-configuration in PsPTS1 to give maackiain.

Enzymatic Synthesis of Structured Lipids Enriched with Medium- and Long-Chain Triacylglycerols via Pickering Emulsion-Assisted Interfacial Catalysis: A Preliminary Exploration

Medium- and long-chain triacylglycerol (MLCT), as a novel functional lipid, is valuable due to its special nutritional properties. Its low content in natural resources and inefficient synthesis during preparation have limited its practical applications. In this study, we developed an effective Pickering emulsion interfacial catalysis system (PE system) for the enzymatic synthesis of MLCT by trans-esterification. Lipase NS 40086 served simultaneously as a catalyst and a solid emulsifier to stabilize the Pickering emulsion. Benefitting from the sufficient oil-water interface, the obtained PE system exhibited outstanding catalytic efficiency, achieving 77.5% of MLCT content within 30 min, 26% higher than that of a water-free system. The Km value (0.259 mM) and activation energy (14.45 kJ mol-1) were 6.8-fold and 1.6-fold lower than those of the water-free system, respectively. The kinetic parameters as well as the molecular dynamics simulation and the tunnel analysis implied that the oil-water interface enhanced the binding between substrate and lipase and thus boosted catalytic efficiency. The conformational changes in the lipase were further explored by FT-IR. This method could give a novel strategy for enhancing lipase activity and the design of efficient catalytic systems to produce added-value lipids. This work will open a new methodology for the enzymatic synthesis of structured lipids.

Hormonal steroids induce multidrug resistance and stress response genes in Neisseria gonorrhoeae by binding to MtrR

Transcriptional regulator MtrR inhibits the expression of the multidrug efflux pump operon mtrCDE in the pathogenic bacterium Neisseria gonorrhoeae. Here, we show that MtrR binds the hormonal steroids progesterone, β-estradiol, and testosterone, which are present at urogenital infection sites, as well as ethinyl estrogen, a component of some hormonal contraceptives. Steroid binding leads to the decreased affinity of MtrR for cognate DNA, increased mtrCDE expression, and enhanced antimicrobial resistance. Furthermore, we solve crystal structures of MtrR bound to each steroid, thus revealing their binding mechanisms and the conformational changes that induce MtrR.

In silico investigation of cholesterol-lowering drugs to find potential inhibitors of dehydrosqualene synthase in Staphylococcus aureus

Staphylococcus aureus is a lethal pathogen that can cause various bacterial infections. This study targets the CrtM enzyme of S. aureus, which is crucial for synthesizing golden carotenoid pigment: staphyloxanthin, which provides anti-oxidant activity to this bacterium for combating antimicrobial resistance inside the host cell. The present investigation quests for human SQS inhibitors against the CrtM enzyme by employing structure-based drug design approaches including induced fit docking (IFD), molecular dynamic (MD) simulations, and binding free energy calculations. Depending upon the docking scores, two compounds, lapaquistat acetate and squalestatin analog 20, were identified as the lead molecules exhibit higher affinity toward the CrtM enzyme. These docked complexes were further subjected to 100 ns MD simulation and several thermodynamics parameters were analyzed. Further, the binding free energies (ΔG) were calculated for each simulated protein-ligand complex to study the stability of molecular contacts using the MM-GBSA approach. Pre-ADMET analysis was conducted for systematic evaluation of physicochemical and medicinal chemistry properties of these compounds. The above study suggested that lapaquistat acetate and squalestatin analog 20 can be selected as potential lead candidates with promising binding affinity for the S. aureus CrtM enzyme. This study might provide insights into the discovery of potential drug candidates for S. aureus with a high therapeutic index.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03862-y.

Tunnel engineering of gas-converting enzymes for inhibitor retardation and substrate acceleration

Carbon monoxide dehydrogenase (CODH), formate dehydrogenase (FDH), hydrogenase (H2ase), and nitrogenase (N2ase) are crucial enzymatic catalysts that facilitate the conversion of industrially significant gases such as CO, CO2, H2, and N2. The tunnels in the gas-converting enzymes serve as conduits for these low molecular weight gases to access deeply buried catalytic sites. The identification of the substrate tunnels is imperative for comprehending the substrate selectivity mechanism underlying these gas-converting enzymes. This knowledge also holds substantial value for industrial applications, particularly in addressing the challenges associated with separation and utilization of byproduct gases. In this comprehensive review, we delve into the emerging field of tunnel engineering, presenting a range of approaches and analyses. Additionally, we propose methodologies for the systematic design of enzymes, with the ultimate goal of advancing protein engineering strategies.

Emerging Strategies for Modifying Cytochrome P450 Monooxygenases into Peroxizymes

ConspectusCytochrome P450 monooxygenase is a versatile oxidizing enzyme with great potential in synthetic chemistry and biology. However, the dependence of its catalytic function on the nicotinamide cofactor NAD(P)H and redox partner proteins limits the practical catalytic application of P450 in vitro. An alternative to expensive cofactors is low-cost H2O2, which can be used directly to exploit the catalytic potential of P450s. However, the peroxide shunt pathway is generally inefficient at driving P450 catalysis compared to normal NAD(P)H-dependent activity. Over the last few decades, the scientific community has made continuous efforts to use directed evolution or site-directed mutagenesis to modify P450 monooxygenases into their peroxizyme modes─peroxygenase and peroxidase. Despite significant progress, obtaining efficient P450 peroxizymes remains a huge challenge. Here, we summarize our efforts to modulate peroxizyme activity in P450 monooxygenases and exploit their catalytic applications in challenging selective C-H oxidation, oxygenation, and oxyfunctionalization over the past seven years. We first developed a dual-functional small molecule (DFSM) strategy for transforming P450BM3 monooxygenase into peroxygenase. In this strategy, the typical DFSM, such as N-(ω-imidazolyl)-hexanoyl-l-phenylalanine (Im-C6-Phe), binds to the P450BM3 protein with an anchoring group at one end and plays a general acid-base catalytic role in the activation of H2O2 with an imidazolyl group at the other end. Compared with the O-O homolysis mechanism in the absence of DFSM, the addition of DFSM efficiently enables the heterolytic O-O cleavage of the adduct Fe-O-OH, thus being favored for the formation of active species compound I, which has been demonstrated by combining crystallographic and theoretical calculations. Furthermore, protein engineering showed the unique catalytic performance of DFSM-facilitated P450 peroxygenase for the highly difficult selective oxidation of C-H bonds. This catalytic performance was demonstrated during the chemoselective hydroxylation of gaseous alkanes, regioselective O-demethylation of aryl ethers, highly (R)-enantioselective epoxidation of styrene, and regio- and enantiomerically diverse hydroxylation of alkylbenzenes. Second, we demonstrated that DFSM-facilitated P450BM3 peroxygenase could be effectively switched to an efficient peroxidase mode through mechanism-guided protein engineering of redox-sensitive residues. Utilizing the peroxidase function of P450 enabled the direct nitration of unsaturated hydrocarbons including phenols, aromatic amines, and styrene derivatives, which was not only the P450-catalyzed direct nitration of phenols and aromatic amines for the first time but also the first example of the direct biological nitration of olefins. Finally, we report an H2O2 tunnel engineering strategy to enable peroxygenase activity in several different P450 monooxygenases for the first time, providing a general approach for accessing engineered P450 peroxygenases. In this Account, we highlight the emerging strategies we have developed for producing practical P450 peroxizyme biocatalysts. Although the DFSM strategy is primarily applied to P450BM3 to date, both strategies of redox-sensitive residue engineering and H2O2 tunnel engineering show great potential to extend to other P450s. These strategies have expanded the scope of applications of P450 chemistry and catalysis. Additionally, they provide a unique solution to the challenging selective oxidation of inert C-H bonds in synthetic chemistry.

The shortest path method (SPM) webserver for computational enzyme design

SPMweb is the online webserver of the Shortest Path Map (SPM) tool for identifying the key conformationally-relevant positions of a given enzyme structure and dynamics. The server is built on top of the DynaComm.py code and enables the calculation and visualization of the SPM pathways. SPMweb is easy-to-use as it only requires three input files: the three-dimensional structure of the protein of interest, and the two matrices (distance and correlation) previously computed from a Molecular Dynamics simulation. We provide in this publication information on how to generate the files for SPM construction even for non-expert users and discuss the most relevant parameters that can be modified. The tool is extremely fast (it takes less than one minute per job), thus allowing the rapid identification of distal positions connected to the active site pocket of the enzyme. SPM applications expand from computational enzyme design, especially if combined with other tools to identify the preferred substitution at the identified position, but also to rationalizing allosteric regulation, and even cryptic pocket identification for drug discovery. The simple user interface and setup make the SPM tool accessible to the whole scientific community. SPMweb is freely available for academia at http://spmosuna.com/.

ChannelsDB 2.0: a comprehensive database of protein tunnels and pores in AlphaFold era

ChannelsDB 2.0 is an updated database providing structural information about the position, geometry and physicochemical properties of protein channels-tunnels and pores-within deposited biomacromolecular structures from PDB and AlphaFoldDB databases. The newly deposited information originated from several sources. Firstly, we included data calculated using a popular CAVER tool to complement the data obtained using original MOLE tool for detection and analysis of protein tunnels and pores. Secondly, we added tunnels starting from cofactors within the AlphaFill database to enlarge the scope of the database to protein models based on Uniprot. This has enlarged available channel annotations ∼4.6 times as of 1 September 2023. The database stores information about geometrical features, e.g. length and radius, and physico-chemical properties based on channel-lining amino acids. The stored data are interlinked with the available UniProt mutation annotation data. ChannelsDB 2.0 provides an excellent resource for deep analysis of the role of biomacromolecular tunnels and pores. The database is available free of charge: https://channelsdb2.biodata.ceitec.cz.

Structure and function of Campylobacter jejuni polynucleotide phosphorylase (PNPase): Insights into the role of this RNase in pathogenicity

Ribonucleases are in charge of the processing, degradation and quality control of all cellular transcripts, which makes them crucial factors in RNA regulation. This post-transcriptional regulation allows bacteria to promptly react to different stress conditions and growth phase transitions, and also to produce the required virulence factors in pathogenic bacteria. Campylobacter jejuni is the main responsible for human gastroenteritis in the world. In this foodborne pathogen, exoribonuclease PNPase (CjPNP) is essential for low-temperature cell survival, affects the synthesis of proteins involved in virulence and has an important role in swimming, cell adhesion/invasion ability, and chick colonization. Here we report the crystallographic structure of CjPNP, complemented with SAXS, which confirms the characteristic doughnut-shaped trimeric arrangement and evaluates domain arrangement and flexibility. Mutations in highly conserved residues were constructed to access their role in RNA degradation and polymerization. Surprisingly, we found two mutations that altered CjPNP into a protein that is only capable of degrading RNA even in conditions that favour polymerization. These findings will be important to develop new strategies to combat C. jejuni infections.

Highly efficient site-specific protein modification using tyrosinase from Streptomyces avermitilis: Structural insight

Tyrosinase-mediated protein conjugation has recently drawn attention as a site-specific protein modification tool under mild conditions. However, the tyrosinases reported to date act only on extremely exposed tyrosine residues, which limits where the target tyrosine can be located. Herein, we report a tyrosinase from Streptomyces avermitilis (SaTYR), that exhibits a much higher activity against tyrosine residues on the protein surface than other tyrosinases. We determined the crystal structure of SaTYR and revealed that the enzyme has a relatively flat and shallow substrate-binding pocket to accommodate a protein substrate. We demonstrated SaTYR-mediated fluorescence dye tagging and PEGylation of a surface tyrosine residue that was unreacted by other tyrosinases with an approximately 95.2 % conjugation yield in 1 h. We also present a structural rationale that considers the steric hindrance from adjacent residues and surrounding structures along with the extent of solvent exposure of residues, as necessary when determining the optimal positions for introducing target tyrosine residues in SaTYR-mediated protein modification. The study demonstrated that the novel tyrosinase, SaTYR, extends the scope of tyrosinase-mediated protein modification, and we propose that site-specific tyrosine conjugation using SaTYR is a promising strategy for protein bioconjugation in various applications.

A cytochrome P450 insecticide detoxification mechanism is not conserved across the Megachilidae family of bees

Recent work has demonstrated that many bee species have specific cytochrome P450 enzymes (P450s) that can efficiently detoxify certain insecticides. The presence of these P450s, belonging or closely related to the CYP9Q subfamily (CYP9Q-related), is generally well conserved across the diversity of bees. However, the alfalfa leafcutter bee, Megachile rotundata, lacks CYP9Q-related P450s and is 170-2500 times more sensitive to certain insecticides than bee pollinators with these P450s. The extent to which these findings apply to other Megachilidae bee species remains uncertain. To address this knowledge gap, we sequenced the transcriptomes of four Megachile species and leveraged the data obtained, in combination with publicly available genomic data, to investigate the evolution and function of P450s in the Megachilidae. Our analyses reveal that several Megachilidae species, belonging to the Lithurgini, Megachilini and Anthidini tribes, including all species of the Megachile genus investigated, lack CYP9Q-related genes. In place of these genes Megachile species have evolved phylogenetically distinct CYP9 genes, the CYP9DM lineage. Functional expression of these P450s from M. rotundata reveal they lack the capacity to metabolize the neonicotinoid insecticides thiacloprid and imidacloprid. In contrast, species from the Osmiini and Dioxyini tribes of Megachilidae have CYP9Q-related P450s belonging to the CYP9BU subfamily that are able to detoxify thiacloprid. These findings provide new insight into the evolution of P450s that act as key determinants of insecticide sensitivity in bees and have important applied implications for pesticide risk assessment.

Recent advances in structure-based enzyme engineering for functional reconstruction

Structural information can help engineer enzymes. Usually, specific amino acids in particular regions are targeted for functional reconstruction to enhance the catalytic performance, including activity, stereoselectivity, and thermostability. Appropriate selection of target sites is the key to structure-based design, which requires elucidation of the structure-function relationships. Here, we summarize the mutations of residues in different specific regions, including active center, access tunnels, and flexible loops, on fine-tuning the catalytic performance of enzymes, and discuss the effects of altering the local structural environment on the functions. In addition, we keep up with the recent progress of structure-based approaches for enzyme engineering, aiming to provide some guidance on how to take advantage of the structural information.

Illuminating the mechanism and allosteric behavior of NanoLuc luciferase

NanoLuc, a superior β-barrel fold luciferase, was engineered 10 years ago but the nature of its catalysis remains puzzling. Here experimental and computational techniques are combined, revealing that imidazopyrazinone luciferins bind to an intra-barrel catalytic site but also to an allosteric site shaped on the enzyme surface. Structurally, binding to the allosteric site prevents simultaneous binding to the catalytic site, and vice versa, through concerted conformational changes. We demonstrate that restructuration of the allosteric site can boost the luminescent reaction in the remote active site. Mechanistically, an intra-barrel arginine coordinates the imidazopyrazinone component of luciferin, which reacts with O2 via a radical charge-transfer mechanism, and then it also protonates the resulting excited amide product to form a light-emitting neutral species. Concomitantly, an aspartate, supported by two tyrosines, fine-tunes the blue color emitter to secure a high emission intensity. This information is critical to engineering the next-generation of ultrasensitive bioluminescent reporters.

PredictONCO: a web tool supporting decision-making in precision oncology by extending the bioinformatics predictions with advanced computing and machine learning

PredictONCO 1.0 is a unique web server that analyzes effects of mutations on proteins frequently altered in various cancer types. The server can assess the impact of mutations on the protein sequential and structural properties and apply a virtual screening to identify potential inhibitors that could be used as a highly individualized therapeutic approach, possibly based on the drug repurposing. PredictONCO integrates predictive algorithms and state-of-the-art computational tools combined with information from established databases. The user interface was carefully designed for the target specialists in precision oncology, molecular pathology, clinical genetics and clinical sciences. The tool summarizes the effect of the mutation on protein stability and function and currently covers 44 common oncological targets. The binding affinities of Food and Drug Administration/ European Medicines Agency -approved drugs with the wild-type and mutant proteins are calculated to facilitate treatment decisions. The reliability of predictions was confirmed against 108 clinically validated mutations. The server provides a fast and compact output, ideal for the often time-sensitive decision-making process in oncology. Three use cases of missense mutations, (i) K22A in cyclin-dependent kinase 4 identified in melanoma, (ii) E1197K mutation in anaplastic lymphoma kinase 4 identified in lung carcinoma and (iii) V765A mutation in epidermal growth factor receptor in a patient with congenital mismatch repair deficiency highlight how the tool can increase levels of confidence regarding the pathogenicity of the variants and identify the most effective inhibitors. The server is available at https://loschmidt.chemi.muni.cz/predictonco.

Unveiling the Catalytic Mechanism of a Processive Metalloaminopeptidase

Intracellular leucine aminopeptidases (PepA) are metalloproteases from the family M17. These enzymes catalyze peptide bond cleavage, removing N-terminal residues from peptide and protein substrates, with consequences for protein homeostasis and quality control. While general mechanistic studies using model substrates have been conducted on PepA enzymes from various organisms, specific information about their substrate preferences and promiscuity, choice of metal, activation mechanisms, and the steps that limit steady-state turnover remain unexplored. Here, we dissected the catalytic and chemical mechanisms of PaPepA: a leucine aminopeptidase from Pseudomonas aeruginosa. Cleavage assays using peptides and small-molecule substrate mimics allowed us to propose a mechanism for catalysis. Steady-state and pre-steady-state kinetics, pH rate profiles, solvent kinetic isotope effects, and biophysical techniques were used to evaluate metal binding and activation. This revealed that metal binding to a tight affinity site is insufficient for enzyme activity; binding to a weaker affinity site is essential for catalysis. Progress curves for peptide hydrolysis and crystal structures of free and inhibitor-bound PaPepA revealed that PaPepA cleaves peptide substrates in a processive manner. We propose three distinct modes for activity regulation: tight packing of PaPepA in a hexameric assembly controls substrate length and reaction processivity; the product leucine acts as an inhibitor, and the high concentration of metal ions required for activation limits catalytic turnover. Our work uncovers catalysis by a metalloaminopeptidase, revealing the intricacies of metal activation and substrate selection. This will pave the way for a deeper understanding of metalloenzymes and processive peptidases/proteases.

Analysis of Zinc and Stromal Immunity in Disuse Osteoporosis: Mendelian Randomization and Transcriptomic Analysis

OBJECTIVE

Disuse osteoporosis is known to be primarily caused by a lack of exercise. However, the causal relationships between zinc and immunity and disuse osteoporosis remain unknown. This study investigated these relationships and their potential mechanisms.

METHODS

This study was an integrative study combining genome-wide association studies and transcriptomics. Two-sample Mendelian randomization analysis (MR) was used to analyze the causal relationships between exposures (zinc, immunity, physical activity) and the outcome (osteoporosis) with the aid of single-nucleotide polymorphisms (SNPs) as instrumental variables (IVs). Four models, MR-Egger, inverse variance weighted, weighted median and MR-Pleiotrophy RESidual Sum and Outlier (MRPRESSO), were used to calculate odds ratio values. Sensitivity and heterogeneity analyses were also performed using MRPRESSO and MR-Egger methods. The mRNA transcriptomic analysis was subsequently conducted. Zinc metabolism scores were acquired through single-sample Gene Set Enrichment Analysis algorithms. Stromal scores were obtained using the R Package "estimate" algorithms. Important Kyoto Encyclopedia of Genes and Genomes and Gene Ontology pathways were also derived through gene set variation analysis. Cytoscape software helped construct the transcription factor (TF)-mRNA-microRNA (miRNA) network. Virtual screening and molecular docking were performed. Polymerase chain reaction validation was also carried out in vivo.

RESULTS

Causal relationships were demonstrated between zinc and exercise (95% confidence interval [CI] = 1.30-2.95, p = 0.001), exercise and immunity (95% CI = 0.36-0.80, p = 0.002), exercise and osteoporosis (95% CI = 0.97-0.99, p = 0.0007), and immunity disorder and osteoporosis (95% CI = 1.30-2.03, p = 0.00002). One hundred and seventy-nine mRNAs in important modules were screened. Combining the differential expressional genes (DEGs) and the Boruta selection, six DEGs were screened (AHNAK, CSF2, ADAMTS12, SRA1, RUNX2, and SLC39A14). TF HOXC10 and miRNA hsa-miR-204 were predicted. Then, the TF-mRNA-miRNA network was successfully constructed. RUNX2 and SLC39A14 were identified as hub mRNAs in the TF-mRNA-miRNA network. Eventually, the novel small drug C6O4NH5 was designed according to the pharmacophore structure of SLC39A14. The docking energy for the novel drug was -5.83 kcal/mol. SLC39A14 and RUNX2 were downregulated (of statistical significance p-value < 0.05) in our animal experiment.

CONCLUSION

This study revealed that zinc had a protective causal relationship with disuse osteoporosis by promoting exercise and immunity. SLC39A14 and RUNX2 mRNA participated in this zinc-related mechanism.

CLICK-chemoproteomics and molecular dynamics simulation reveals pregnenolone targets and their binding conformations in Th2 cells

Pregnenolone (P5) is synthesized as the first bioactive steroid in the mitochondria from cholesterol. Clusters of differentiation 4 (CD4+) and Clusters of differentiation 8 (CD8+) immune cells synthesize P5 de novo; P5, in turn, play important role in immune homeostasis and regulation. However, P5's biochemical mode of action in immune cells is still emerging. We envisage that revealing the complete spectrum of P5 target proteins in immune cells would have multifold applications, not only in basic understanding of steroids biochemistry in immune cells but also in developing new therapeutic applications. We employed a CLICK-enabled probe to capture P5-binding proteins in live T helper cell type 2 (Th2) cells. Subsequently, using high-throughput quantitative proteomics, we identified the P5 interactome in CD4+ Th2 cells. Our study revealed P5's mode of action in CD4+ immune cells. We identified novel proteins from mitochondrial and endoplasmic reticulum membranes to be the primary mediators of P5's biochemistry in CD4+ and to concur with our earlier finding in CD8+ immune cells. Applying advanced computational algorithms and molecular simulations, we were able to generate near-native maps of P5-protein key molecular interactions. We showed bonds and interactions between key amino acids and P5, which revealed the importance of ionic bond, hydrophobic interactions, and water channels. We point out that our results can lead to designing of novel molecular therapeutics strategies.

Awakening the natural capability of psicose production in Escherichia coli

Due to the rampant rise in obesity and diabetes, consumers are desperately seeking for ways to reduce their sugar intake, but to date there are no options that are both accessible and without sacrifice of palatability. One of the most promising new ingredients in the food system as a non-nutritive sugar substitute with near perfect palatability is D-psicose. D-psicose is currently produced using an in vitro enzymatic isomerization of D-fructose, resulting in low yield and purity, and therefore requiring substantial downstream processing to obtain a high purity product. This has made adoption of D-psicose into products limited and results in significantly higher per unit costs, reducing accessibility to those most in need. Here, we found that Escherichia coli natively possesses a thermodynamically favorable pathway to produce D-psicose from D-glucose through a series of phosphorylation-epimerization-dephosphorylation steps. To increase carbon flux towards D-psicose production, we introduced a series of genetic modifications to pathway enzymes, central carbon metabolism, and competing metabolic pathways. In an attempt to maximize both cellular viability and D-psicose production, we implemented methods for the dynamic regulation of key genes including clustered regularly interspaced short palindromic repeats inhibition (CRISPRi) and stationary-phase promoters. The engineered strains achieved complete consumption of D-glucose and production of D-psicose, at a titer of 15.3 g L-1, productivity of 2 g L-1 h-1, and yield of 62% under test tube conditions. These results demonstrate the viability of whole-cell catalysis as a sustainable alternative to in vitro enzymatic synthesis for the accessible production of D-psicose.

EzMechanism: an automated tool to propose catalytic mechanisms of enzyme reactions

Over the years, hundreds of enzyme reaction mechanisms have been studied using experimental and simulation methods. This rich literature on biological catalysis is now ripe for use as the foundation of new knowledge-based approaches to investigate enzyme mechanisms. Here, we present a tool able to automatically infer mechanistic paths for a given three-dimensional active site and enzyme reaction, based on a set of catalytic rules compiled from the Mechanism and Catalytic Site Atlas, a database of enzyme mechanisms. EzMechanism (pronounced as 'Easy' Mechanism) is available to everyone through a web user interface. When studying a mechanism, EzMechanism facilitates and improves the generation of hypotheses, by making sure that relevant information is considered, as derived from the literature on both related and unrelated enzymes. We validated EzMechanism on a set of 62 enzymes and have identified paths for further improvement, including the need for additional and more generic catalytic rules.

Computational insights of the molecular recognition between volatile molecules and odorant binding proteins from the red palm weevil Rhynchophorus ferrugineus

The red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae) is one of the most harmful pests for palm trees, causing serious economic damage worldwide. The present work aims to model and study the 3D structures of highly expressed odorant binding proteins from R. ferrugineus (RferOBPs) and identify possible binding modes and ligand release mechanism by docking and molecular dynamics. Highly confident 3D structures of a total of 11 odorant binding proteins (OBPs) were obtained with AlphaFold2. All 3D RferOBPs modeled structures displayed six characteristic α-helices, except for RfeOBP7 and RfeOBP10, which had an extra terminal α-helix. Among the eleven modeled RferOBPs, RferOBP4 was highly expressed in the antennae and subsequently selected for further analyses. Molecular docking analyses demonstrated that ferruginol, α-pinene, DEET, and picaridin can favorably bind the RferOBP4 cavity with low affinity energies. Molecular dynamic simulations of RferOBP4 bound to ferruginol at different pH values showed that low pH environments dictate a structural change into an apo-state that modifies the number of tunnels where the ligand can coexist, further triggering ligand release by a pH-dependent mechanism. This is the first report concerning the modelling and study of ligand binding modes and release mechanism of R. ferrugineus OBPs.Communicated by Ramaswamy H. Sarma.

Fishing for Catalysis: Experimental Approaches to Narrowing Search Space in Directed Evolution of Enzymes

Directed evolution has transformed protein engineering offering a path to rapid improvement of protein properties. Yet, in practice it is limited by the hyper-astronomic protein sequence search space, and approaches to identify mutagenic hot spots, i.e., locations where mutations are most likely to have a productive impact, are needed. In this perspective, we categorize and discuss recent progress in the experimental approaches (broadly defined as structural, bioinformatic, and dynamic) to hot spot identification. Recent successes in harnessing protein dynamics and machine learning approaches provide new opportunities for the field and will undoubtedly help directed evolution reach its full potential.

Structure-Based Design of Small Imine Reductase Panels for Target Substrates

Biocatalysis is important in the discovery, development, and manufacture of pharmaceuticals. However, the identification of enzymes for target transformations of interest requires major screening efforts. Here, we report a structure-based computational workflow to prioritize protein sequences by a score based on predicted activities on substrates, thereby reducing a resource-intensive laboratory-based biocatalyst screening. We selected imine reductases (IREDs) as a class of biocatalysts to illustrate the application of the computational workflow termed IREDFisher. Validation by using published data showed that IREDFisher can retrieve the best enzymes and increase the hit rate by identifying the top 20 ranked sequences. The power of IREDFisher is confirmed by computationally screening 1400 sequences for chosen reductive amination reactions with different levels of complexity. Highly active IREDs were identified by only testing 20 samples in vitro. Our speed test shows that it only takes 90 min to rank 85 sequences from user input and 30 min for the established IREDFisher database containing 591 IRED sequences. IREDFisher is available as a user-friendly web interface (https://enzymeevolver.com/IREDFisher). IREDFisher enables the rapid discovery of IREDs for applications in synthesis and directed evolution studies, with minimal time and resource expenditure. Future use of the workflow with other enzyme families could be implemented following the modification of the workflow scoring function.

Worth the Weight: Sub-Pocket EXplorer (SubPEx), a Weighted Ensemble Method to Enhance Binding-Pocket Conformational Sampling

Structure-based virtual screening (VS) is an effective method for identifying potential small-molecule ligands, but traditional VS approaches consider only a single binding-pocket conformation. Consequently, they struggle to identify ligands that bind to alternate conformations. Ensemble docking helps address this issue by incorporating multiple conformations into the docking process, but it depends on methods that can thoroughly explore pocket flexibility. We here introduce Sub-Pocket EXplorer (SubPEx), an approach that uses weighted ensemble (WE) path sampling to accelerate binding-pocket sampling. As proof of principle, we apply SubPEx to three proteins relevant to drug discovery: heat shock protein 90, influenza neuraminidase, and yeast hexokinase 2. SubPEx is available free of charge without registration under the terms of the open-source MIT license: http://durrantlab.com/subpex/.

Immobilization of Candida antarctica lipase B on ILs modified CNTs with different chain lengths: Regulation of substrate tunnel "Leucine gating"

Ionic liquids (ILs) have been widely used as chemical modifiers to modify the carriers and thus improve the efficiency, activity and stability of the enzymes. However, as thousands of ILs have been found up to date, it's a huge work for screening and designing suitable ILs for immobilization of enzymes. Moreover, the mechanism of improving enzymes catalytic performance is still remain ambiguous. Thus, this study investigated the impact of ILs with different chain lengths on the enzymatic properties of Candida antarctica lipase B (CALB). Molecular dynamics simulations were employed to examine the interaction between ILs modified CNTs and CALB, as well as their effects on CALB's structure. The results revealed that ILs with different chain lengths significantly influenced the absorption orientation of CALB. Tunnel analysis identified a key role for Leu278 in regulating the open or closed state of Tunnel 2 during CALB's catalytic cycle. The weak interaction analysis demonstrated that ILs with suitable chain lengths provided spatial freedom and formed strong interactions with CNTs and ILs (vdW and hbond). This led to a conformational flip of Leu278, stabilizing the open state of Tunnel 2 and improving the activity and stability of immobilized CALB. This study provides novel insights into the design of new green modifiers to modulate carrier performance and obtain immobilized enzymes with better performance, and establishes a theoretical basis for the design and selection of modifiers for ILs in future work.

Induction of AHR Signaling in Response to the Indolimine Class of Microbial Stress Metabolites

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor that plays an important role in gastrointestinal barrier function, tumorigenesis, and is an emerging drug target. The resident microbiota is capable of metabolizing tryptophan to metabolites that are AHR ligands (e.g., indole-3-acetate). Recently, a novel set of mutagenic tryptophan metabolites named indolimines have been identified that are produced by M. morganii in the gastrointestinal tract. Here, we determined that indolimine-200, -214, and -248 are direct AHR ligands that can induce Cyp1a1 transcription and subsequent CYP1A1 enzymatic activity capable of metabolizing the carcinogen benzo(a)pyrene in microsomal assays. In addition, indolimines enhance IL6 expression in a colonic tumor cell line in combination with cytokine treatment. The concentration of indolimine-248 that induces AHR transcriptional activity failed to increase DNA damage. These observations reveal an additional aspect of how indolimines may alter colonic tumorigenesis beyond mutagenic activity.

Heterologous Expression and Biochemical Characterization of a New Chloroperoxidase Isolated from the Deep-Sea Hydrothermal Vent Black Yeast Hortaea werneckii UBOCC-A-208029

The initiation of this study relies on a targeted genome-mining approach to highlight the presence of a putative vanadium-dependent haloperoxidase-encoding gene in the deep-sea hydrothermal vent fungus Hortaea werneckii UBOCC-A-208029. To date, only three fungal vanadium-dependent haloperoxidases have been described, one from the terrestrial species Curvularia inaequalis, one from the fungal plant pathogen Botrytis cinerea, and one from a marine derived isolate identified as Alternaria didymospora. In this study, we describe a new vanadium chloroperoxidase from the black yeast H. werneckii, successfully cloned and overexpressed in a bacterial host, which possesses higher affinity for bromide (Km = 26 µM) than chloride (Km = 237 mM). The enzyme was biochemically characterized, and we have evaluated its potential for biocatalysis by determining its stability and tolerance in organic solvents. We also describe its potential three-dimensional structure by building a model using the AlphaFold 2 artificial intelligence tool. This model shows some conservation of the 3D structure of the active site compared to the vanadium chloroperoxidase from C. inaequalis but it also highlights some differences in the active site entrance and the volume of the active site pocket, underlining its originality.

Investigation of molecular details of a bacterial cationic amino acid transporter (GkApcT) during arginine transportation using molecular dynamics simulation and umbrella sampling techniques

CONTEXT

Cationic amino acid transporters (CATs) facilitate arginine transport across membranes and maintain its levels in various tissues and organs, but their overexpression has been associated with severe cancers. A recent study identified the alternating access mechanism and critical residues involved in arginine transportation in a cationic amino acid transporter from Geobacillus kaustophilus (GkApcT). Here, we used molecular dynamics (MD) simulation methods to investigate the transportation mechanism of arginine (Arg) through GkApcT. The results revealed that arginine strongly interacts with specific binding site residues (Thr43, Asp111, Glu115, Lys191, Phe231, Ile234, and Asp237). Based on the umbrella sampling, the main driving force for arginine transport is the polar interactions of the arginine with channel-lining residues. An in-depth description of the dissociation mechanism and binding energy analysis brings valuable insight into the interactions between arginine and transporter residues, facilitating the design of effective CAT inhibitors in cancer cells.

METHODS

The membrane-protein system was constructed by uploading the prokaryotic CAT (PDB ID: 6F34) to the CHARMM-GUI web server. Molecular dynamics simulations were done using the GROMACS package, version 5.1.4, with the CHARMM36 force field and TIP3P water model. The MM-PBSA approach was performed for determining the arginine binding free energy. Furthermore, the hotspot residues were identified through per-residue decomposition analysis. The characteristics of the channel such as bottleneck radius and channel length were analyzed using the CaverWeb 1.1 web server. The proton wire inside the transporter was investigated based on the classic Grotthuss mechanism. We also investigated the atomistic details of arginine transportation using the path-based free energy umbrella sampling technique (US).

Identification of key amino acid residues in AtUMAMIT29 for transport of glucosinolates

Glucosinolates are key defense compounds of plants in Brassicales order, and their accumulation in seeds is essential for the protection of the next generation. Recently, members of the Usually Multiple Amino acids Move In and Out Transporter (UMAMIT) family were shown to be essential for facilitating transport of seed-bound glucosinolates from site of synthesis within the reproductive organ to seeds. Here, we set out to identify amino acid residues responsible for glucosinolate transport activity of the main seed glucosinolate exporter UMAMIT29 in Arabidopsis thaliana. Based on a predicted model of UMAMIT29, we propose that the substrate transporting cavity consists of 51 residues, of which four are highly conserved residues across all the analyzed homologs of UMAMIT29. A comparison of the putative substrate binding site of homologs within the brassicaceous-specific, glucosinolate-transporting clade with the non-brassicaceous-specific, non-glucosinolate-transporting UMAMIT32 clade identified 11 differentially conserved sites. When each of the 11 residues of UMAMIT29 was individually mutated into the corresponding residue in UMAMIT32, five mutant variants (UMAMIT29#V27F, UMAMIT29#M86V, UMAMIT29#L109V, UMAMIT29#Q263S, and UMAMIT29#T267Y) reduced glucosinolate transport activity over 75% compared to wild-type UMAMIT29. This suggests that these residues are key for UMAMIT29-mediated glucosinolate transport activity and thus potential targets for blocking the transport of glucosinolates to the seeds.

Fungal P450 Deconstructs the 2,5-Diazabicyclo[2.2.2]octane Ring En Route to the Complete Biosynthesis of 21R-Citrinadin A

Prenylated indole alkaloids (PIAs) possess great structural diversity and show biological activities. Despite significant efforts in investigating the biosynthetic mechanism, the key step in the transformation of 2,5-diazabicyclo[2.2.2]octane-containing PIAs into a distinct class of pentacyclic compounds remains unknown. Here, using a combination of gene deletion, heterologous expression, and biochemical characterization, we show that a unique fungal P450 enzyme CtdY catalyzes the cleavage of the amide bond in the 2,5-diazabicyclo[2.2.2]octane system, followed by a decarboxylation step to form the 6/5/5/6/6 pentacyclic ring in 21R-citrinadin A. We also demonstrate the function of a subsequent cascade of stereospecific oxygenases to further modify the 6/5/5/6/6 pentacyclic intermediate en route to the complete 21R-citrinadin A biosynthesis. Our findings reveal a key enzyme CtdY for the pathway divergence in the biosynthesis of PIAs and uncover the complex late-stage post-translational modifications in 21R-citrinadin A biosynthesis.

Comparison of Fungal Thermophilic and Mesophilic Catalase-Peroxidases for Their Antioxidative Properties

Catalase-peroxidases (KatGs) are unique bifunctional oxidoreductases that contain heme in their active centers allowing both the peroxidatic and catalatic reaction modes. These originally bacterial enzymes are broadly distributed among various fungi allowing them to cope with reactive oxygen species present in the environment or inside the cells. We used various biophysical, biochemical, and bioinformatics methods to investigate differences between catalase-peroxidases originating in thermophilic and mesophilic fungi from different habitats. Our results indicate that the architecture of the active center with a specific post-translational modification is highly similar in mesophilic and thermophilic KatG and also the peroxidatic acitivity with ABTS, guaiacol, and L-DOPA. However, only the thermophilic variant CthedisKatG reveals increased manganese peroxidase activity at elevated temperatures. The catalatic activity releasing molecular oxygen is comparable between CthedisKatG and mesophilic MagKatG1 over a broad temperature range. Two constructed point mutations in the active center were performed selectively blocking the formation of described post-translational modification in the active center. They exhibited a total loss of catalatic activity and changes in the peroxidatic activity. Our results indicate the capacity of bifunctional heme enzymes in the variable reactivity for potential biotech applications.

Siderophore Synthetase DesD Catalyzes N-to-C Condensation in Desferrioxamine Biosynthesis

Desferrioxamine siderophores are assembled by the nonribosomal-peptide-synthetase-independent siderophore (NIS) synthetase enzyme DesD via ATP-dependent iterative condensation of three N1-hydroxy-N1-succinyl-cadaverine (HSC) units. Current knowledge of NIS enzymology and the desferrioxamine biosynthetic pathway does not account for the existence of most known members of this natural product family, which differ in substitution patterns of the N- and C-termini. The directionality of desferrioxamine biosynthetic assembly, N-to-C versus C-to-N, is a longstanding knowledge gap that is limiting further progress in understanding the origins of natural products in this structural family. Here, we establish the directionality of desferrioxamine biosynthesis using a chemoenzymatic approach with stable isotope incorporation and dimeric substrates. We propose a mechanism where DesD catalyzes the N-to-C condensation of HSC units to establish a unifying biosynthetic paradigm for desferrioxamine natural products in Streptomyces.

Hormonal steroids bind the Neisseria gonorrhoeae multidrug resistance regulator, MtrR, to induce a multidrug binding efflux pump and stress-response sigma factor

Overexpression of the multidrug efflux pump MtrCDE, a critical factor of multidrug-resistance in Neisseria gonorrhoeae , the causative agent of gonorrheae, is repressed by the transcriptional regulator, MtrR (multiple transferable resistance repressor). Here, we report the results from a series of in vitro experiments to identify innate, human inducers of MtrR and to understand the biochemical and structural mechanisms of the gene regulatory function of MtrR. Isothermal titration calorimetry experiments reveal that MtrR binds the hormonal steroids progesterone, β-estradiol, and testosterone, all of which are present at significant concentrations at urogenital infection sites as well as ethinyl estrogen, a component of some birth control pills. Binding of these steroids results in decreased affinity of MtrR for cognate DNA, as demonstrated by fluorescence polarization-based assays. The crystal structures of MtrR bound to each steroid provided insight into the flexibility of the binding pocket, elucidated specific residue-ligand interactions, and revealed the conformational consequences of the induction mechanism of MtrR. Three residues, D171, W136 and R176 are key to the specific binding of these gonadal steroids. These studies provide a molecular understanding of the transcriptional regulation by MtrR that promotes N. gonorrhoeae survival in its human host.

KVFinder-web: a web-based application for detecting and characterizing biomolecular cavities

Molecular interactions that modulate catalytic processes occur mainly in cavities throughout the molecular surface. Such interactions occur with specific small molecules due to geometric and physicochemical complementarity with the receptor. In this scenario, we present KVFinder-web, an open-source web-based application of parKVFinder software for cavity detection and characterization of biomolecular structures. The KVFinder-web has two independent components: a RESTful web service and a web graphical portal. Our web service, KVFinder-web service, handles client requests, manages accepted jobs, and performs cavity detection and characterization on accepted jobs. Our graphical web portal, KVFinder-web portal, provides a simple and straightforward page for cavity analysis, which customizes detection parameters, submits jobs to the web service component, and displays cavities and characterizations. We provide a publicly available KVFinder-web at https://kvfinder-web.cnpem.br, running in a cloud environment as docker containers. Further, this deployment type allows KVFinder-web components to be configured locally and customized according to user demand. Hence, users may run jobs on a locally configured service or our public KVFinder-web.

Isoform Specific Regulation of Adenylyl Cyclase 5 by Gβγ

The nine different membrane-anchored adenylyl cyclase isoforms (AC1-9) in mammals are stimulated by the heterotrimeric G protein Gαs, but their response to Gβγ regulation is isoform-specific. For example, AC5 is conditionally activated by Gβγ. Here, we report cryo-EM structures of ligand-free AC5 in complex with Gβγ and of a dimeric form of AC5 that could be involved in its regulation. Gβγ binds to a coiled-coil domain that links the AC transmembrane region to its catalytic core as well as to a region (C1b) that is known to be a hub for isoform-specific regulation. We confirmed the Gβγ interaction with both purified proteins and cell-based assays. The interface with Gβγ involves AC5 residues that are subject to gain-of-function mutations in humans with familial dyskinesia, indicating that the observed interaction is important for motor function. A molecular mechanism wherein Gβγ either prevents dimerization of AC5 or allosterically modulates the coiled-coil domain, and hence the catalytic core, is proposed. Because our mechanistic understanding of how individual AC isoforms are uniquely regulated is limited, studies such as this may provide new avenues for isoform-specific drug development.

Water potential governs the effector specificity of the transcriptional regulator XylR of Pseudomonas putida

The biodegradative capacity of bacteria in their natural habitats is affected by water availability. In this work, we have examined the activity and effector specificity of the transcriptional regulator XylR of the TOL plasmid pWW0 of Pseudomonas putida mt-2 for biodegradation of m-xylene when external water potential was manipulated with polyethylene glycol PEG8000. By using non-disruptive luxCDEAB reporter technology, we noticed that the promoter activated by XylR (Pu) restricted its activity and the regulator became more effector-specific towards head TOL substrates when cells were grown under water subsaturation. Such a tight specificity brought about by water limitation was relaxed when intracellular osmotic stress was counteracted by the external addition of the compatible solute glycine betaine. With these facts in hand, XylR variants isolated earlier as effector-specificity responders to the non-substrate 1,2,4-trichlorobenzene under high matric stress were re-examined and found to be unaffected by water potential in vivo. All these phenomena could be ultimately explained as the result of water potential-dependent conformational changes in the A domain of XylR and its effector-binding pocket, as suggested by AlphaFold prediction of protein structures. The consequences of this scenario for the evolution of specificities in regulators and the emergence of catabolic pathways are discussed.

Soybean (Glycine max L.) Lipoxygenase 1 (LOX 1) Is Modulated by Nitric Oxide and Hydrogen Sulfide: An In Vitro Approach

Hydrogen sulfide (H₂S) and nitric oxide (NO) are two relevant signal molecules that can affect protein function throughout post-translational modifications (PTMs) such as persulfidation, S-nitrosation, metal-nitrosylation, and nitration. Lipoxygenases (LOXs) are a group of non-heme iron enzymes involved in a wide range of plant physiological functions including seed germination, plant growth and development, and fruit ripening and senescence. Likewise, LOXs are also involved in the mechanisms of response to diverse environmental stresses. Using purified soybean (Glycine max L.) lipoxygenase type 1 (LOX 1) and nitrosocysteine (CysNO) and sodium hydrosulfide (NaHS) as NO and H₂S donors, respectively, the present study reveals that both compounds negatively affect LOX activity, suggesting that S-nitrosation and persulfidation are involved. Mass spectrometric analysis of nitrated soybean LOX 1 using a peroxynitrite (ONOO^-) donor enabled us to identify that, among the thirty-five tyrosine residues present in this enzyme, only Y214 was exclusively nitrated by ONOO^-. The nitration of Y214 seems to affect its interaction with W500, a residue involved in the substrate binding site. The analysis of the structure 3PZW demonstrates the existence of several tunnels that directly communicate the surface of the protein with different internal cysteines, thus making feasible their potential persulfidation, especially C429 and C127. On the other hand, the CysNO molecule, which is hydrophilic and bulkier than H₂S, can somehow be accommodated throughout the tunnel until it reaches C127, thus facilitating its nitrosation. Overall, a large number of potential persulfidation targets and the ease by which H₂S can reach them through the diffuse tunneling network could be behind their efficient inhibition.

Modulating Substrate Specificity of Rhizobium sp. Histamine Dehydrogenase through Protein Engineering for Food Quality Applications

Histamine is a biogenic amine found in fish-derived and fermented food products with physiological relevance since its concentration is proportional to food spoilage and health risk for sensitive consumers. There are various analytical methods for histamine quantification from food samples; however, a simple and quick enzymatic detection and quantification method is highly desirable. Histamine dehydrogenase (HDH) is a candidate for enzymatic histamine detection; however, other biogenic amines can change its activity or produce false positive results with an observed substrate inhibition at higher concentrations. In this work, we studied the effect of site saturation mutagenesis in Rhizobium sp. Histamine Dehydrogenase (Rsp HDH) in nine amino acid positions selected through structural alignment analysis, substrate docking, and proximity to the proposed histamine-binding site. The resulting libraries were screened for histamine and agmatine activity. Variants from two libraries (positions 72 and 110) showed improved histamine/agmatine activity ratio, decreased substrate inhibition, and maintained thermal resistance. In addition, activity characterization of the identified Phe72Thr and Asn110Val HDH variants showed a clear substrate inhibition curve for histamine and modified kinetic parameters. The observed maximum velocity (V_max) increased for variant Phe72Thr at the cost of an increased value for the Michaelis-Menten constant (K_m) for histamine. The increased K_m value, decreased substrate inhibition, and biogenic amine interference observed for variant Phe72Thr support a tradeoff between substrate affinity and substrate inhibition in the catalytic mechanism of HDHs. Considering this tradeoff for future enzyme engineering of HDH could lead to breakthroughs in performance increases and understanding of this enzyme class.