CAVER 3.0: a tool for the analysis of transport pathways in dynamic protein structures

Enhancing cryo-enzymatic efficiency in cold-adapted lipase from Psychrobacter sp. C18 via site-directed mutagenesis

As industrial demands for cold-active enzymes have been increased, psychrophilic lipases present a promising solution with potential for innovation and growth in food, pharmaceutical, and detergent industries. Cold-adapted enzymes achieve high catalytic efficiency at low temperatures through their structural flexibility and conformational adaptability. Therefore, in this study, the lipase gene from Psychrobacter sp. C18 was cloned and subjected to site-directed mutagenesis based on computer aided predictions to enhance the enzyme's cold-adapted properties and flexibility. Mutations were strategically selected in loops of the active site to improve the enzyme's accessibility to the substrate under cold conditions. The P163G, L186G, and Q239W mutations were selected for further analysis. Enzyme activity, along with its stability and structural flexibility, was assessed using techniques including UV-Vis spectroscopy, fluorescence, and circular dichroism (CD) spectroscopy. The obtained data revealed that the optimal temperature for the wild-type lipase was 30 °C, which shifted to lower temperatures in the mutants: 15 °C for P163G and L186G, and 20 °C for Q239W. Additionally, the optimal pH of the mutant lipases shifted to more alkaline conditions compared to the wild-type enzyme. While the thermal and pH stability of the mutant enzymes slightly decreased, these findings can be attributed to their enhanced flexibility. Far-UV CD spectroscopy revealed a reduction in α-helical content of the mutant enzymes. Molecular dynamics simulations corroborated these findings, confirming increased structural flexibility in all three mutants compared to the wild-type enzyme. This research underlines the importance of applying engineered cold-adapted enzymes for industrial application.

Crystal structures of the putative endoribonuclease L-PSP from Entamoeba histolytica

Entamoeba histolytica causes amebiasis, a neglected disease that kills ∼100 000 people globally each year. Due to emerging drug resistance, E. histolytica is one of the target organisms for structure-based drug discovery by the Seattle Structural Genomics Center for Infectious Disease (SSGCID). Purification, crystallization and three structures of the putative drug target endoribonuclease L-PSP from E. histolytica (EhL-PSP) are presented. EhL-PSP has a two-layer α/β-sandwich with structural homology to endoribonuclease L-PSP. All three structures reveal the prototypical YjgF/YER057c/UK114 family trimer topology with accessible allosteric active sites. Citrate molecules from the crystallization solution are bound to the allosteric site in two of the three reported structures. The large allosteric site of EhL-PSP is well conserved with bacterial YjgF/YER057c/UK114 family members and could be targeted for inhibition, drug discovery or repurposing.

Sucrose and Gibberellic Acid Binding Stabilize the Inward-Open Conformation of AtSWEET13: A Molecular Dynamics Study

In plants, sugar will eventually be exported transporters (SWEETs) facilitate the translocation of mono- and disaccharides across membranes and play a critical role in modulating responses to gibberellin (GA3), a key growth hormone. However, the dynamic mechanisms underlying sucrose and GA3 binding and transport remain elusive. Here, we employed microsecond-scale molecular dynamics (MD) simulations to investigate the influence of sucrose and GA3 binding on SWEET13 transporter motions. While sucrose exhibits high flexibility within the binding pocket, GA3 remains firmly anchored in the central cavity. Binding of both ligands increases the average channel radius along the transporter's principal axis. In contrast to the apo form, which retains its initial conformation throughout the simulation, ligand-bound complexes undergo a significant conformational transition characterized by further opening of the intracellular gate relative to the inward-open crystal structure (5XPD). This opening is driven by ligand-induced bending of helix V, stabilizing the inward-open state. Sucrose binding notably enhances the flexibility of the intracellular gate and amplifies anticorrelated motions between the N- and C-terminal domains at the intracellular side, suggesting an opening-closing motion of these domains. Principal component analysis revealed that this gating motion is most pronounced in the sucrose complex and minimal in the apo form, highlighting sucrose's ability to induce high-amplitude gating. Our binding free energy calculations indicate that SWEET13 has lower binding affinity for sucrose compared to GA3, consistent with its role in sugar transport. These results provide insight into key residues involved in sucrose and GA3 binding and transport, advancing our understanding of SWEET13 dynamics.

An Effective Computational Strategy for UGTs Catalytic Function Prediction

The GT-B type glycosyltransferases play a crucial post-modification role in synthesizing natural products, such as triterpenoid and steroidal saponins, renowned for their diverse pharmacological activities. Despite phylogenetic analysis aiding in enzyme family classification, distinguishing substrate specificity between triterpenoid and steroidal saponins, with their highly similar cyclic scaffolds, remains a formidable challenge. Our studies unveil the potential transport tunnels for the glycosyl donor and acceptor in PpUGT73CR1, by molecular dynamics simulations. This revelation leads to a plausible substrate transport mechanism, highlighting the regulatory role of the N-terminal domain (NTD) in glycosyl acceptor binding and transport. Inspired by these structural and mechanistic insights, we further analyze the binding pockets of 44 plant-derived UGTs known to glycosylate triterpenes and sterols. Notably, sterol UGTs are found to harbor aromatic and hydrophobic residues with polar residues typically present at the bottom of the active pocket. Drawing inspiration from the substrate binding and product release mechanism revealed through structure-based molecular modeling, we devised a fast sequence-based method for classifying UGTs using the pre-trained ESM2 protein model. This method involved extracting the NTD features of UGTs and performing PCA clustering analysis, enabling accurate identification of enzyme function, and even differentiation of substrate specificity/promiscuity between structurally similar triterpenoid and steroidal substrates, which is further validated by experiments. This work not only deepens our understanding of substrate binding mechanisms but also provides an effective computational protocol for predicting the catalytic function of unknown UGTs.

The identification of XPR1 as a voltage- and phosphate-activated phosphate-permeable ion channel

Maintaining a balance of inorganic phosphate (Pi) is vital for cellular functionality. Proper phosphate levels are managed through Pi import and export; and the processes governing Pi export remain the least understood. Xenotropic and Polytropic retrovirus Receptor 1 (XPR1) has been identified as the only known Pi export protein in mammals. In this study, we introduce the cryogenic electron microscopy structure of human XPR1 (hXPR1), unveiling a structural arrangement distinct from that of any known ion transporter. Our structural results suggest that hXPR1 may operate as an ion channel, a hypothesis supported by patch clamp recordings revealing hXPR1's voltage- and Pi-dependent activity and large unitary conductance. Further analyses, including the structure of hXPR1 in presence of Pi, and mutagenesis studies at one of the putative Pi binding sites, lead us to propose a plausible ion permeation pathway. Together, our results provide novel perspectives on the Pi transport mechanism of XPR1.

Biochemical mapping reveals a conserved heme transport mechanism via CcmCD in System I bacterial cytochrome c biogenesis

Heme is a redox-active cofactor for essential processes across all domains of life. Heme's redox capabilities are responsible for its biological significance but also make it highly cytotoxic, requiring tight intracellular regulation. Thus, the mechanisms of heme trafficking are still not well understood. To address this, the bacterial cytochrome c biogenesis pathways are being developed into model systems to elucidate mechanisms of heme trafficking. These pathways function to attach heme to apocytochrome c, which requires the transport of heme from inside to outside of the cell. Here, we focus on the System I pathway (CcmABCDEFGH) which is proposed to function in two steps: CcmABCD transports heme across the membrane and attaches it to CcmE. HoloCcmE then transports heme to the holocytochrome c synthase, CcmFH, for attachment to apocytochrome c. To interrogate heme transport across the membrane, we focus on CcmCD, which can form holoCcmE independent of CcmAB, leading to the hypothesis that CcmCD is a heme transporter. A structure-function analysis via cysteine/heme crosslinking identified a heme acceptance domain and heme transport channel in CcmCD. Bioinformatic analysis and structural predictions across prokaryotic organisms determined that the heme acceptance domains are structurally variable, potentially to interact with diverse heme delivery proteins. In contrast, the CcmC transmembrane heme channel is structurally conserved, indicating a common mechanism for transmembrane heme transport. We provide direct biochemical evidence mapping the CcmCD heme channel and providing insights into general mechanisms of heme trafficking by other putative heme transporters.

IMPORTANCE

Heme is a biologically important cofactor for proteins involved with essential cellular functions, such as oxygen transport and energy production. Heme can also be toxic to cells and thus requires tight regulation and specific trafficking pathways. As a result, much effort has been devoted to understanding how this important, yet cytotoxic, molecule is transported. While several heme transporters/importers/exporters have been identified, the biochemical mechanisms of transport are not well understood, representing a major knowledge gap. Here, the bacterial cytochrome c biogenesis pathway, System I (CcmABCDEFGH), is used to elucidate the transmembrane transport of heme via CcmCD. We utilize a cysteine/heme crosslinking approach, which can trap endogenous heme in specific domains, to biochemically map the heme transport channel in CcmCD, demonstrating that CcmCD is a heme transporter. These results suggest a model for heme trafficking by other heme transporters in both prokaryotes and eukaryotes.

Surveying Enzyme Crystal Structures Reveals the Commonality of Active-Site Solvent Accessibility and Enzymatic Water Networks

Despite the demonstrable dependence of enzyme functionality on solvation, the notion of water being directly chemically required for catalysis inside active sites remains unexplored. Here we report that over 99% of 1013 enzyme crystals obtained by X-ray crystallography with high resolution (<1.5 Å) contain continuous chains of water linking residues within the active site to bulk water. Also reported are the findings which inspired this study-that electric fields experienced by water hydrogen atoms are on average twice as strong in the active sites of both chains of bacterial polynucleotide kinase (PDB 4QM6) structures compared to those in bulk water. These results point to the possibility that water molecules within active sites may be paramount to the immense catalytic power of enzymes, especially for mechanisms requiring hydronium or hydroxide ions.

Evolutionary insights into the selectivity of sterol oxidising cytochrome P450 enzymes based on ancestral sequence reconstruction

The cytochrome P450 (CYP) enzyme CYP125A1 is a crucial enzyme for the long-term survival and pathogenicity of Mycobacterium tuberculosis. CYP125 genes are found not only in pathogenic mycobacteria but are also widely dispersed within the Actinobacteria phylum, with many species possessing multiple copies of CYP125 encoding genes. Their primary function is the catalytic hydroxylation of the terminal methyl group of cholesterol and phytosterols. We have previously shown that CYP125 enzymes from distinct mycobacteria have substrate selectivity preferences for animal versus plant steroid oxidation. An evolutionary understanding of this selectivity is not known. Here, we use Ancestral Sequence Reconstruction (ASR), to support the hypothesis that some CYP125 enzymes evolved in a manner reflective of their adaptation to a pathogenic niche. We constructed a maximum-likelihood, most-recent common ancestor of the CYP125 clade (CYP125MRCA). We were then able to produce and characterise this enzyme both functionally and structurally. We found that CYP125MRCA was able to catalyse the terminal hydroxylation of cholesterol, phytosterols, and vitamin D3 (cholecalciferol); the latter was hydroxylated at both C-25 and C-26. This is the first example to date of vitamin D3 oxidation by a CYP125 enzyme, thereby demonstrating an increased substrate range of CYP125MRCA relative to its characterised extant relatives. The X-ray crystal structures of CYP125MRCA bound with sitosterol and vitamin D3 were determined, providing important insight into the changes that enable the expanded substrate range.

Structural insights into auxin influx mediated by the Arabidopsis AUX1

Auxin is crucial in orchestrating diverse aspects of plant growth and development and modulating responses to environmental signals. The asymmetric spatiotemporal distribution of auxin generates local gradient patterns, which are regulated by both cellular auxin influx and efflux. The AUXIN1/LIKE-AUX1 (AUX1/LAX) family transporters have been identified as major auxin influx carriers. Here, we characterize the auxin uptake mediated by AUX1 from Arabidopsis thaliana. Using cryoelectron microscopy (cryo-EM), we determine its structure in three states: the auxin-unbound, the auxin-bound, and the competitive inhibitor, 3-chloro-4-hydroxyphenylacetic acid (CHPAA)-bound state. All structures adopt an inward-facing conformation. In the auxin-bound structure, indole-3-acetic acid (IAA) is coordinated to AUX1 primarily through hydrogen bonds with its carboxyl group. The functional roles of key residues in IAA binding are validated by in vitro and in planta analyses. CHPAA binds to the same site as IAA. These findings advance our understanding of auxin transport in plants.

Protein Engineering for Enhancing Electron Transfer in P450-Mediated Catalysis

Cytochrome P450 enzymes (P450s) are versatile biocatalysts with applications spanning pharmaceutical development and natural product biosynthesis. A critical bottleneck in P450-mediated reactions is the electron transfer process, which often limits catalytic efficiency and promotes uncoupling events leading to reactive oxygen species (ROS) formation. This review comprehensively examines recent protein engineering strategies aimed at enhancing electron transfer efficiency in P450 systems. We explore the design and application of different fusion constructs, which improve proximity between the P450 enzyme and its redox partners (RPs), as well as scaffold-mediated protein assembly, enabling precise spatial organization of P450s and RPs. Furthermore, we discuss targeted modifications at the P450-RP interaction interface and optimization of electron transfer pathways through site-directed mutagenesis and directed evolution. These strategies enhance catalytic activity, improve coupling efficiency, and reduce ROS formation. Finally, we address the remaining challenges in understanding and engineering P450 electron transfer, and discuss the future directions, emphasizing the need for advanced computational modeling, structural characterization, and integration of synthetic and systems biology approaches.

Caver Web 2.0: analysis of tunnels and ligand transport in dynamic ensembles of proteins

Enzymes with buried active sites utilize molecular tunnels to exchange substrates, products, and solvent molecules with the surface. These transport mechanisms are crucial for protein function and influence various properties. As proteins are inherently dynamic, their tunnels also vary structurally. Understanding these dynamics is essential for elucidating structure-function relationships, drug discovery, and bioengineering. Caver Web 2.0 is a user-friendly web server that retains all Caver Web 1.0 functionalities while introducing key improvements: (i) generation of dynamic ensembles via automated molecular dynamics with YASARA, (ii) analysis of dynamic tunnels with CAVER 3.0, (iii) prediction of ligand trajectories in multiple snapshots with CaverDock 1.2, and (iv) customizable ligand libraries for virtual screening. Users can assess protein flexibility, identify and characterize tunnels, and predict ligand trajectories and energy profiles in both static and dynamic structures. Additionally, the platform supports virtual screening with FDA/EMA-approved drugs and user-defined datasets. Caver Web 2.0 is a versatile tool for biological research, protein engineering, and drug discovery, aiding the identification of strong inhibitors or new substrates to bind to the active sites or tunnels, and supporting drug repurposing efforts. The server is freely accessible at https://loschmidt.chemi.muni.cz/caverweb.

Biochemical, Structural, and Conformational Characterization of a Fungal Ethylene-Forming Enzyme

The ethylene-forming enzyme (EFE) from the fungus Penicillium digitatum strain Pd1 was heterologously produced in Escherichia coli and its properties were compared to the extensively characterized bacterial enzyme from Pseudomonas savastanoi strain PK2. Both enzymes catalyze four reactions: the conversion of 2-oxoglutarate (2OG) to ethylene and CO2, oxidative decarboxylation of 2OG coupled to l-arginine (l-Arg) hydroxylation, uncoupled oxidative decarboxylation of 2OG, and the production of 3-hydroxypropionate (3-HP) from 2OG. The strain Pd1 enzyme exhibited a greater ratio of ethylene production over l-Arg hydroxylation than the PK2 strain EFE. The uncoupled decarboxylation of 2OG and 3-HP production are minor reactions in both cases, but they occur to a greater extent using the fungal enzyme. Additional distinctions of the fungal versus bacterial enzyme are noted in the absorbance maxima and l-Arg dependence of their anaerobic electronic spectra. The structures of the Pd1 EFE apoprotein and the EFE·Mn(II)·2OG complex resembled the corresponding structures of the PK2 enzyme, but notable structural differences were observed in the computationally predicted Pd1 EFE·Fe(II)·2OG·l-Arg complex versus the PK2 EFE·Mn(II)·2OG·l-Arg crystal structure. These studies extend our biochemical understanding and represent the first structural and conformational characterization of a eukaryotic EFE.

Molecular cloning, overexpression, characterization, and mechanism explanation of an esterase RasEst3 for ester synthesis under aqueous phase

Fatty acid esters are widely used in fragrance compounds, solvents, lubricants, and biofuels. Enzymatic synthesis of these esters in aqueous phase is an environmentally friendly approach. In this study, an esterase RasEst3 from Rasamsonia emersonii was identified for fatty acid ester synthesis through sequence alignment. The gene encoding RasEst3 was heterologously expressed in Escherichia coli BL21(DE3), and its enzymatic properties were analyzed. The enzyme exhibited optimal activity at pH 3.5 and 30 °C, with a preference for medium-chain substrates. Structurally, RasEst3 contains a lid domain and a catalytic domain, with a catalytic triad composed of Ser146-His227-Asp214. The smaller pocket spatial site resistance and the hydrophobicity of the substrate channel facilitate effective substrate binding to the active center. Site-directed mutagenesis and molecular dynamics simulations revealed that the oxygen anion holes formed by Gly69 and Thr70, along with the π-bond stacking formed by Tyr112 and Tyr145, play crucial roles in catalysis. After removing a loop region from RasEst3, its ethyl octanoate synthesis activity increased by 253.22 % compared to the wild-type enzyme. This study not only clarifies the structure-function relationship of RasEst3 but also provides valuable insights for developing novel biocatalysts in green chemistry.

A leigh syndrome mutation perturbs long-range energy coupling in respiratory complex I

Respiratory complex I is a central enzyme of cellular energy metabolism that couples electron transfer with proton translocation across a biological membrane. In doing so, it powers oxidative phosphorylation that drives energy consuming processes. Mutations in complex I lead to severe neurodegenerative diseases in humans. However, the biochemical consequences of these mutations remain largely unknown. Here, we use the Escherichia coli complex I as a model to biochemically characterize the F124LMT-ND5 mutation found in patients suffering from Leigh syndrome. We show that the mutation drastically perturbs proton translocation and electron transfer activities to the same extent, despite the remarkable 140 Å distance between the mutated position and the electron transfer domain. Our molecular dynamics simulations suggest that the disease-causing mutation induces conformational changes that hamper the propagation of an electric wave through an ion-paired network essential for proton translocation. Our findings imply that malfunction of the proton translocation domain is entirely transmitted to the electron transfer domain underlining the action-at-a-distance coupling in the proton-coupled electron transfer of respiratory complex I.

Structural Basis for the Catalytic Mechanism of ATP-Dependent Diazotase CmaA6

Although several diazotases have been recently reported, the details of the reaction mechanism are not yet understood. In this study, we investigated the mechanism of CmaA6, an ATP-dependent diazotase, which catalyzes the diazotization of 3-aminocoumaric acid using nitrous acid. X-ray crystallography and cryogenic electron microscopy-single particle analysis revealed CmaA6 structures in the substrate-free and AMP-binding states. Kinetic analysis suggested that CmaA6 catalyzes diazotization via a sequential reaction mechanism in which three substrates (nitrous acid, ATP, and 3-aminocoumaric acid) are simultaneously bound in the reaction pocket. The nitrous acid and 3-aminocoumaric acid binding sites were predicted based on the AMP-binding state and confirmed by site-directed mutagenesis. In addition, computational analysis revealed a tunnel for 3-aminocoumaric acid to enter the reaction pocket, which was advantageous for the sequential reaction mechanism. This study provides important insights into the catalytic mechanism of diazotization in natural product biosynthesis.

Probing substrate water access through the O1 channel of Photosystem II by single site mutations and membrane inlet mass spectrometry

Light-driven water oxidation by photosystem II sustains life on Earth by providing the electrons and protons for the reduction of CO2 to carbohydrates and the molecular oxygen we breathe. The inorganic core of the oxygen evolving complex is made of the earth-abundant elements manganese, calcium and oxygen (Mn4CaO5 cluster), and is situated in a binding pocket that is connected to the aqueous surrounding via water-filled channels that allow water intake and proton egress. Recent serial crystallography and infrared spectroscopy studies performed with PSII isolated from Thermosynechococcus vestitus (T. vestitus) support that one of these channels, the O1 channel, facilitates water access to the Mn4CaO5 cluster during its S2→S3 and S3→S4→S0 state transitions, while a subsequent CryoEM study concluded that this channel is blocked in the cyanobacterium Synechocystis sp. PCC 6803, questioning the role of the O1 channel in water delivery. Employing site-directed mutagenesis we modified the two O1 channel bottleneck residues D1-E329 and CP43-V410 (T. vestitus numbering) and probed water access and substrate exchange via time resolved membrane inlet mass spectrometry. Our data demonstrates that water reaches the Mn4CaO5 cluster via the O1 channel in both wildtype and mutant PSII. In addition, the detailed analysis provides functional insight into the intricate protein-water-cofactor network near the Mn4CaO5 cluster that includes the pentameric, near planar 'water wheel' of the O1 channel.

Decoding Allosteric Effects of Missense Variations in Drug Metabolism: Afrocentric CYP3A4 Alleles Explored

There is growing research on the allosteric behaviour of proteins, including studies on allosteric mutations that contribute to human diseases and the development of allosteric drugs. Allostery also plays a key role in drug metabolism, an essential factor in drug development. However, population specific variations, particularly in 3D protein structures, remain understudied. This study focuses on CYP3A4, a key enzyme responsible for metabolizing over 50% of FDA-approved drugs and often linked to adverse drug reactions. Given the vast genetic diversity of Africa, we investigated 13 CYP3A4 alleles from African populations using post-molecular dynamics analyses, with 12 being single variations and one containing a double variation. Except for one, all allele variations were located away from the active site, suggesting allosteric effects. Our comparative analyses of reference and variant structures, through hydrogen bond interactions, dynamic residue network analysis and substrate channel dynamics, revealed notable differences at both global and residue levels. The *32-I335T variant showed the largest changes compared to the reference structure, while *3-M445T (near normal metabolizer) exhibited the least change, with other variants falling in between. The *32-I335T variant showed a distorted conformation in the radius of gyration, a distinct kink in the I helix with specific hydrogen bonds and altered channel patterns. The *12-L373F variant, associated with reduced metabolism of midazolam and quinine, showed increased rigidity in its vicinity, potentially interfering with catalytic activity. Our findings align with clinical and wet lab data, suggesting that our approaches could be applied to analyse variants without clinical evidence.

Mechanistic Study of Protein Interaction with Natto Inhibitory Peptides Targeting Xanthine Oxidase: Insights from Machine Learning and Molecular Dynamics Simulations

Bioactive peptides from food sources offer a safe and biocompatible approach to enzyme inhibition, with potential applications in managing metabolic disorders such as hyperuricemia and gout, conditions linked to excessive xanthine oxidase activity. Using a machine learning-based screening approach inspired by the bioactivity of natto, two peptides, ECFK and FECK, were identified from the Bacillus subtilis proteome and validated as xanthine oxidase inhibitors with IC50 values of 37.36 and 71.57 mM, respectively. Further experiments confirmed their safety through cytotoxicity assays, and electronic tongue analysis demonstrated their mild sensory properties, supporting their edibility. Molecular dynamics simulations revealed that these peptides stabilize critical enzyme regions, with ECFK showing a higher dissociation energy barrier (52.08 kcal/mol) than FECK (46.39 kcal/mol), indicating strong, stable interactions. This study highlights food-derived peptides as safe and natural inhibitors of xanthine oxidase, offering promising therapeutic potential for metabolic disorder management.

Discovering cryptic pocket opening and binding of a stimulant derivative in a vestibular site of the 5-HT3A receptor

A diverse set of modulators, including stimulants and anesthetics, regulates ion channel function in our nervous system. However, structures of ligand-bound complexes can be difficult to capture by experimental methods, particularly when binding is dynamic. Here, we used computational methods and electrophysiology to identify a possible bound state of a modulatory stimulant derivative in a cryptic vestibular pocket of a mammalian serotonin-3 receptor. We first applied a molecular dynamics simulation-based goal-oriented adaptive sampling method to identify possible open-pocket conformations, followed by Boltzmann docking that combines traditional docking with Markov state modeling. Clustering and analysis of stability and accessibility of docked poses supported a preferred binding site; we further validated this site by mutagenesis and electrophysiology, suggesting a mechanism of potentiation by stabilizing intersubunit contacts. Given the pharmaceutical relevance of serotonin-3 receptors in emesis, psychiatric, and gastrointestinal diseases, characterizing relatively unexplored modulatory sites such as these could open valuable avenues to understanding conformational cycling and designing state-dependent drugs.

Polyethylenimine-Assisted Interfacial Modulation Based on Electrostatic Balancing and Hierarchical Channel to the Cytidine 5'-Monophosphate Conversion Performance of Uridine-Cytosine Kinase

Interfacial modulation of protein microenvironments plays a pivotal role in enhancing enzyme immobilization and catalysis. In this study, we proposed a polyethylenimine (PEI)-assisted strategy that combines electrostatic and affinity interactions to improve the performance of Cytidine 5'-Monophosphate (CMP) conversion by uridine-cytidine kinase (UCK). The PEI-modified interface creates an optimal local microenvironment that maintains a balanced charge distribution, stabilizes UCK's conformation, and prevents denaturation. Electrostatic interactions promote product adsorption, enhance diffusion, and reduce substrate accumulation, boosting reaction efficiency. The kinetic assays revealed an increase in the maximum reaction rate from 16.8 to 113.2 μM·min-1 with a remarkable increase in substrate affinity and enzyme activity. The relative enzyme activity at the optimal substrate concentration increased from 70.4 to 106.9%, and by 113.3% under conditions of substrate inhibition. This study provides theoretical and technical support for the efficient production of CMP with promising applications in food, feed, and medical fields.

Generative Landscapes and Dynamics to Design Multidomain Artificial Transmembrane Transporters

Protein design is challenging as it requires simultaneous consideration of interconnected factors, such as fold, dynamics, and function. These evolutionary constraints are encoded in protein sequences and can be learned through the latent generative landscape (LGL) framework to predict functional sequences by leveraging evolutionary patterns, enabling exploration of uncharted sequence space. By simulating designed proteins through molecular dynamics (MD), we gain deeper insights into the interdependencies governing structure and dynamics. We present a synergized workflow combining LGL with MD and biochemical characterization, allowing us to explore the sequence space effectively. This approach has been applied to design and characterize two artificial multidomain ATP-driven transmembrane copper transporters, with native-like functionality. This integrative approach proved effective in unraveling the intricate relationships between sequence, structure, and function.

Ancestral Sequence Reconstruction and Semirational Engineering of Glycosyltransferase for Efficient Synthesis of Rare Ginsenoside Rh1

Rare ginsenoside Rh1, exhibiting great potential in the food industry, is limited by its natural scarcity. This constraint has driven the development of biocatalytic synthesis approaches, yet robust enzymes capable of efficient production remain elusive. Here, we employed the ancestral sequence reconstruction (ASR) approach to create a thermostable UDP-dependent glycosyltransferase (UGT227) for Rh1 synthesis from 20(S)-protopanaxatriol (PPT). UGT227 exhibited enhanced thermostability (t1/2 = 44.2 h at 60 °C) but initially yielded only 15% Rh1. Semirational engineering generated the I83A/F285 M variant, increasing the yield to 92%. For economic viability, the I83A/F285 M variant was coexpressed with Arabidopsis thaliana sucrose synthase (AtSUS1), enabling the use of cost-effective sucrose for UDP-glucose regeneration. This integration achieved a 99.9% yield at a 1 mM PPT. Molecular dynamics simulations revealed that the enlarged binding pocket entrance of I83A/F285 M contributed to the enhanced Rh1 yield. Our findings offer strategies for efficient biosynthesis of Rh1 and pave the way for economically feasible production.

Peroxisomal H2O2-generating sulfite oxidase (SOX) from pepper fruits is negatively modulated by NO and H2S

Nitric oxide and hydrogen sulfide are signal molecules that can exert regulatory functions in diverse plant processes including fruit ripening. Sulfite oxidase (SOX) is a peroxisomal enzyme that catalyzes the oxidation of sulfite (SO32-) to sulfate (SO42-) with the concomitant generation of H2O2. SOX requires the molybdenum cofactor (Moco) and it has been proposed that SOX functions as a mechanism of protection against sulfite toxicity. Based on the analysis of the pepper genome and fruit transcriptome (RNA-seq), a single gene encoding for a SOX, was identified in chromosome 2. The CaSOX gene expression analysis during fruit ripening, from green immature (G) to red ripe (R) indicates that its expression increased. In-gel analysis using non-denaturing PAGE of a 50-75% (NH4)2SO4 protein fraction allowed the detection of its SOX activity in green pepper fruits. In vitro assay of the SOX from pepper fruits showed that the SOX activity is differently regulated by NO and H2S. Mass spectrometric analysis of the nitrated recombinant pepper SOX enables us to corroborate that this enzyme undergoes inhibition by nitration in Tyr10. Protein modeling analysis also reveals that Cys70 and Cys163 are susceptible targets for S-nitrosation and persulfidation. These findings suggest that NO and H2S could function upstream of the peroxisomal H2O2-generating SOX, highlighting the intricate network of signaling molecules within this subcellular compartment.

A Versatile Protein Scaffold Engineered for the Hierarchical Assembly of Robust and Highly Active Enzymes

Scaffold proteins play immense roles in bringing enzymes together to enhance their properties. However, the direct fusion of scaffold with bulky guest enzymes may disrupt the assembly process or diminish catalytic efficiency. Most self-assembling protein scaffolds are engineered to form structures beforehand, and then carry guest proteins via different conjugation strategies in vitro. Here, a robust self-assembling scaffold is presented, engineered from Methanococcus jannaschii using disulfide bonds, which efficiently assembles bulky enzymes into higher-order helices without additional chemistry or bio-conjugation in vitro. When fused directly with monomeric Endo-1,4-beta-xylanase A, the catalytic efficiency of the guest enzyme increased by 2.5 times with enhanced thermostability. Additionally, integrating the scaffold with the multimeric metalloenzyme nitrile hydratase overcame the typical stability-activity trade-off of such industrial enzyme, yielding three-fold higher activity and 28-fold higher thermostability. Structural analyses suggest that the artificially made helical twist structures create new interface interactions and provide a concentration of active sites of guest enzymes. Further fusion of fluorescent protein pairs with the scaffold exhibited a 12-fold higher FRET efficiency, suggesting its potential for dual-enzyme cascade applications. Overall, this study showcases a simple yet powerful protein scaffold that organizes guest enzymes into hierarchical structures with enhanced catalytic performance.

Conformational dynamics and activation of membrane-associated human Group IVA cytosolic phospholipase A 2 (cPLA 2 )

Cytosolic phospholipase A 2 (cPLA 2 ) associates with membranes where it hydrolyzes phospholipids containing arachidonic acid to initiate an inflammatory cascade. All-atom molecular dynamics simulations were employed to understand the activation process when cPLA 2 associates with the endoplasmic reticulum (ER) membrane of macrophages where it acts. We found that membrane association causes the lid region of cPLA 2 to undergo a closed-to-open state transition that is accompanied by the sideways movement of loop 495-540, allowing the exposure of a cluster of lysine residues (K488, K541, K543, and K544), which binds the allosteric activator PIP 2 in the membrane. The active site of the open form of cPLA 2 , containing the catalytic dyad residues S228 and D549, exhibited a three-fold larger cavity than the closed form of cPLA 2 in aqueous solution. These findings provide mechanistic insight as to how cPLA 2 ER membrane association promotes major transitions between conformational states critical to allosteric activation and enzymatic phospholipid hydrolysis.

Detection of Putative Ligand Dissociation Pathways in Proteins Using Site-Identification by Ligand Competitive Saturation

Drug efficacy often correlates better with dissociation kinetics than binding affinity alone. To study binding kinetics computationally, it is necessary to identify all of the possible ligand dissociation pathways. The site identification by ligand competitive saturation (SILCS) method involves the precomputation of a set of maps (FragMaps), which describe the free energy landscapes of typical chemical functionalities in and around a target protein or RNA. In the current work, we present and implement a method to use SILCS to identify ligand dissociation pathways, termed "SILCS-Pathway." The A* pathfinding algorithm is utilized to enumerate ligand dissociation pathways between the ligand binding site and the surrounding bulk solvent environment defined on evenly spaced points around the protein based on a Fibonacci lattice. The cost function for the A* algorithm is calculated using the SILCS exclusion maps and the SILCS grid free energy scores, thereby identifying paths that account for local protein flexibility and potential favorable interactions with the ligand. By traversing all evenly distributed bulk solvent points around the protein, we located all possible dissociation pathways and clustered them to identify general ligand unbinding pathways. The procedure is verified by using proteins studied previously with enhanced sampling molecular dynamics (MD) techniques and is shown to be capable of capturing important ligand dissociation routes in a highly computationally efficient manner. The identified pathways will act as the foundation for determining ligand dissociation kinetics using SILCS free energy profiles, which will be described in a subsequent article.

Crystal Structure, Modeling, and Identification of Key Residues Provide Insights into the Mechanism of the Key Toxoflavin Biosynthesis Protein ToxD

Toxoflavin, a toxic secondary metabolite produced by a variety of bacteria, has been implicated as a causative agent in food poisoning and a virulence factor in phytopathogenic bacteria. This toxin is produced by genes encoded in the tox operon in Burkholderia glumae, in which the encoded protein, ToxD, was previously characterized as essential for toxoflavin production. To better understand the function of ToxD in toxoflavin biosynthesis and provide a basis for future work to develop inhibitors of ToxD, we undertook the identification of structurally and catalytically important amino acid residues through a combination of X-ray crystallography and site directed mutagenesis. We solved the structure of BgToxD, which crystallized as a dimer, to 1.8 Å resolution. We identified a citrate molecule in the putative active site. To investigate the role of individual residues, we used Pseudomonas protegens Pf-5, a BL1 plant protective bacterium known to produce toxoflavin, and created mutants in the ToxD-homologue PFL1035. Using a multiple sequence alignment and the BgToxD structure, we identified and explored the functional importance of 12 conserved residues in the putative active site. Eight variants of PFL1035 resulted in no observable production of toxoflavin. In contrast, four ToxD variants resulted in reduced but detectable toxoflavin production suggesting a nonessential role. The crystal structure and structural models of the substrate and intermediate bound enzyme provide a molecular interpretation for the mutagenesis data.

From closed to open: three dynamic states of membrane-bound cytochrome P450 3A4

Cytochrome P450 3A4 (CYP3A4) is a membrane bound monooxygenase. It metabolizes the largest proportion of all orally ingested drugs. Ligands can enter and exit the enzyme through flexible tunnels, which co-determine CYP3A4's ligand promiscuity. The flexibility can be represented by distinct conformational states of the enzyme. However, previous state definitions relied solely on crystal structures. We employed conventional molecular dynamics (cMD) simulations to sample the conformational space of CYP3A4. Five conformationally different crystal structures embedded in a membrane were simulated for 1 µs each. A Markov state model (MSM) coupled with spectral clustering (Robust Perron Cluster Analysis PCCA +) resulted in three distinct states: Two open conformations and an intermediate conformation. The tunnels inside CYP3A4 were calculated with CAVER3.0. Notably, we observed variations in bottleneck radii compared to those derived from crystallographic data. We want to point out the importance of simulations to characterize the dynamic behaviour. Moreover, we identified a mechanism, in which the membrane supports the opening of a tunnel. Therefore, CYP3A4 must be investigated in its membrane-bound state.

Activity Regulation of a Glutamine Amidotransferase Bienzyme Complex by Substrate-Induced Subunit Interface Expansion

Glutamine amidotransferases are multienzyme machineries in which reactive ammonia is generated by a glutaminase and then transferred through a sequestered protein tunnel to a synthase active site for incorporation into diverse metabolites. To avoid wasteful metabolite consumption, there is a requirement for synchronized catalysis, but any generally applicable mechanistic insight is still lacking. As synthase activity depends on glutamine turnover, we investigated possible mechanisms controlling glutaminase catalysis using aminodeoxychorismate synthase involved in folate biosynthesis as a model. By analyzing this system in distinct states of catalysis, we found that incubation with glutamine leads to a subunit interface expansion by one-third of its original area. These changes completely enclose the glutaminase active site for sequestered catalysis and the subsequent transport of volatile ammonia to the synthase active site. In view of similar rearrangements in other glutamine amidotransferases, our observations may provide a general mechanism for the catalysis synchronization of this multienzyme family.

Rational multienzyme architecture design with iMARS

Biocatalytic cascades with spatial proximity can orchestrate multistep pathways to form metabolic highways, which enhance the overall catalytic efficiency. However, the effect of spatial organization on catalytic activity is poorly understood, and multienzyme architectural engineering with predictable performance remains unrealized. Here, we developed a standardized framework, called iMARS, to rapidly design the optimal multienzyme architecture by integrating high-throughput activity tests and structural analysis. The approach showed potential for industrial-scale applications, with artificial fusion enzymes designed by iMARS significantly improving the production of resveratrol by 45.1-fold and raspberry ketone by 11.3-fold in vivo, as well as enhancing ergothioneine synthesis in fed-batch fermentation. In addition, iMARS greatly enhanced the in vitro catalytic efficiency of the multienzyme complexes for PET plastic depolymerization and vanillin biosynthesis. As a generalizable and flexible strategy at molecular level, iMARS could greatly facilitate green chemistry, synthetic biology, and biomanufacturing.

Extremophilic hemoglobins: The structure of Shewanella benthica truncated hemoglobin N

Truncated hemoglobins (TrHbs) have an ancient origin and are widely distributed in microorganisms where they often serve roles other than dioxygen transport and storage. In extremophiles, these small heme proteins must have features that secure function under challenging conditions: at minimum, they must be folded, retain the heme group, allow substrates to access the heme cavity, and maintain their quaternary structure if present and essential. The genome of the obligate psychropiezophile Shewanella benthica strain KT99 harbors a gene for a TrHb belonging to a little-studied clade of globins (subgroup 2 of group N). In the present work, we characterized the structure of this protein (SbHbN) with electronic absorption spectroscopy and X-ray crystallography and inspected its structural integrity under hydrostatic pressure with NMR spectroscopy and small-angle X-ray scattering. We found that SbHbN self-associates weakly in solution and contains an extensive network of hydrophobic tunnels connecting the active site to the surface. Amino acid replacements at the dimeric interface formed by helices G and H in the crystal confirmed this region to be the site of intermolecular interactions. High hydrostatic pressure dissociated the assemblies while the porous subunits resisted unfolding and heme loss. Preservation of structural integrity under pressure is also observed in nonpiezophilic TrHbs, which suggests that this ancient property is derived from functional requirements. Added to the inability of SbHbN to combine reversibly with dioxygen and a propensity to form heme d, the study broadens our perception of the TrHb lineage and the resistance of globins to extreme environmental conditions.

Molecular basis of H2O2/O2.-/.OH discrimination during electrochemical activation of DyP peroxidases: The critical role of the distal residues

Here, we show that the replacement of the distal residues Asp and/or Arg of the DyP peroxidases from Bacillus subtilis and Pseudomonas putida results in functional enzymes, albeit with spectroscopically perturbed active sites. All the enzymes can be activated either by the addition of exogenous H2O2 or by in situ electrochemical generation of the reactive oxygen species (ROS) •OH, O2•- and H2O2. The latter method leads to broader and upshifted pH-activity profiles. Both WT enzymes exhibit a differential predominance of ROS involved in their electrochemical activation, which follows the order •OH > O2•- > H2O2 for BsDyP and O2•- > H2O2 > •OH for PpDyP. This ROS selectivity is preserved in mutants with unperturbed sites but is blurred out for distorted sites. The underlying molecular basis of the selectivity mechanisms is analysed through molecular dynamics simulations, which reveal distorted hydrogen bonding networks and higher throughput of the access tunnels in the variants exhibiting no selectivity. The electrochemical activation method provides superior performance for protein variants with a high prevalence of the alternative •OH and O2•- species. These results constitute a promising advance towards engineering DyPs for electrocatalytic applications.

Distal mutations enhance catalysis in designed enzymes by facilitating substrate binding and product release

The role of amino-acid residues distant from an enzyme's active site in facilitating the complete catalytic cycle-including substrate binding, chemical transformation, and product release-remains poorly understood. Here, we investigate how distal mutations promote the catalytic cycle by engineering mutants of three de novo Kemp eliminases containing either active-site or distal mutations identified through directed evolution. Kinetic analyses, X-ray crystallography, and molecular dynamics simulations reveal that while active-site mutations create preorganized catalytic sites for efficient chemical transformation, distal mutations enhance catalysis by facilitating substrate binding and product release through tuning structural dynamics to widen the active-site entrance and reorganize surface loops. These distinct contributions work synergistically to improve overall activity, demonstrating that a well-organized active site, though necessary, is not sufficient for optimal catalysis. Our findings reveal critical roles that distal residues play in shaping the catalytic cycle to enhance efficiency, yielding valuable insights for enzyme design.

Structure of the human TWIK-2 potassium channel and its inhibition by pimozide

The potassium channel TWIK-2 is crucial for ATP-induced activation of the NLRP3 inflammasome in macrophages. The channel is a member of the two-pore domain potassium (K2P) channel superfamily and an emerging therapeutic target to mitigate severe inflammatory injury involving NLRP3 activation. We report the cryo-EM structure of human TWIK-2. In comparison to other K2P channels, the structure reveals a unique 'up' conformation of Tyr111 in the selectivity filter and a SF1-P1 pocket behind the filter that could serve as a binding site for channel modulators. Density for acyl chains is present in fenestrations within the transmembrane region that connect the central cavity of the pore to the lipid membrane. Limited pharmacological tools are available for TWIK-2 despite its importance as a drug target. We show that the small molecule pimozide inhibits TWIK-2 and determine a structure of the channel with pimozide. Pimozide displaces the acyl chains and binds below the selectivity filter to block ion conduction. The drug may access its binding site via the membrane, suggesting that other hydrophobic small molecules could have utility for inhibiting TWIK-2. The work defines the structure of TWIK-2 and provides a structural foundation for development of specific inhibitors with potential utility as anti-inflammatory drugs.

Significance Statement

The TWIK-2 potassium channel is a member of the two-pore domain potassium (K2P) channel superfamily and a potential therapeutic target to control severe inflammatory injury involving the NLRP3 inflammasome. We report the cryo-EM structure of the human TWIK-2 channel at 2.85 Å resolution, revealing differences in comparison to other K2P channels. We identify that pimozide, an FDA-approved drug for Tourette syndrome, inhibits TWIK-2. A cryo-EM structure of TWIK-2 in complex with pimozide identifies its binding location and mechanism of inhibition. The work provides a structural foundation for development of specific TWIK-2 inhibitors that have potential therapeutic utility for inflammatory diseases involving NLRP3 activation.

Improving the hydrophilic microenvironment surrounding the catalytic site of fructosyltransferase enhances its catalytic ability

The hydrophilic microenvironment surrounding an enzyme's active site can influence its catalytic activity. This study examines the effect of enhancing this environment in the Aspergillus niger fructosyltransferase, SucC. Bioinformatics analysis identified a cysteine residue (C66) near the catalytic triad (D64, D194, E271) as vital for maintaining the active site's structure and facilitating substrate transport. Simulated mutagenesis suggested that mutating cysteine to serine (C66S) could increase hydrophilicity without altering the structure significantly. This mutation was predicted to enhance substrate affinity, with binding energy changing from -3.65 to -4.14 kcal mol-1. The C66S mutant, expressed in Pichia pastoris GS115, showed a 61.3% increase in specific activity, a 13.5% decrease in Km (82.20/71.14 mM), and a 21.6% increase in kcat (112.23/136.48 min-1), resulting in a 40.1% increase in catalytic efficiency (1.37/1.92 min-1 mM-1). For fructooligosaccharides (FOS) production, C66S demonstrated enhanced transfructosylation, particularly in the initial stages of the reaction, achieving higher overall FOS yields. These findings highlight that modifying the active site hydrophilicity, without causing major structural changes, is a promising strategy for improving an enzyme's catalytic efficiency.

HydE Catalytic Mechanism Is Powered by a Radical Relay with Redox-Active Fe(I)-Containing Intermediates

[FeFe]-hydrogenases are enzymes that catalyze the redox interconversion of H+ and H2 using a six-iron active site, known as the H-cluster, which consists of a structurally unique [2Fe]H subcluster linked to a [4Fe-4S]H subcluster. A set of enzymes, HydG, HydE, and HydF, are responsible for the biosynthesis of the [2Fe]H subcluster. Among them, it is well established that HydG cleaves tyrosine into CO and CN- and forms a mononuclear [Fe(II)(Cys)(CO)2(CN)] complex. Recent work using EPR spectroscopy and X-ray crystallography show that HydE uses this organometallic Fe complex as its native substrate. The low spin Fe(II) center is reduced into an adenosylated Fe(I) species, which is proposed to form an Fe(I)Fe(I) dimer within HydE. The highly unusual transformation catalyzed by HydE draws interest in both biochemistry and organometallic chemistry. Due to the instability of the substrate, the intermediates, and the proposed product, experimental characterization of the detailed HydE mechanism and its final product is challenging. Herein, the catalytic mechanism of HydE is studied using hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics simulations. A radical relay mechanism was found for the cleavage of the cysteine S-Cβ bond that is energetically favored with respect to a closed-shell mechanism involving unconventional proton transfer. In addition, we propose a pathway for the dimerization of two Fe(I) complexes within the HydE hydrophobic cavity, which is consistent with the recent experimental result that HydF can perform [FeFe]-hydrogenase maturation with a synthetic dimer complex as the substrate. These simulation results take us further down the path to a more complete understanding of these enzymes that synthesize one of Nature's most efficient energy conversion catalysts.

A hybrid machine learning framework for functional annotation of mitochondrial glutathione transport and metabolism proteins in cancers

BACKGROUND

Alterations of metabolism, including changes in mitochondrial metabolism as well as glutathione (GSH) metabolism are a well appreciated hallmark of many cancers. Mitochondrial GSH (mGSH) transport is a poorly characterized aspect of GSH metabolism, which we investigate in the context of cancer. Existing functional annotation approaches from machine (ML) or deep learning (DL) models based only on protein sequences, were unable to annotate functions in biological contexts.

RESULTS

We develop a flexible ML framework for functional annotation from diverse feature data. This hybrid ML framework leverages cancer cell line multi-omics data and other biological knowledge data as features, to uncover potential genes involved in mGSH metabolism and membrane transport in cancers. This framework achieves strong performance across functional annotation tasks and several cell line and primary tumor cancer samples. For our application, classification models predict the known mGSH transporter SLC25A39 but not SLC25A40 as being highly probably related to mGSH metabolism in cancers. SLC25A10, SLC25A50, and orphan SLC25A24, SLC25A43 are predicted to be associated with mGSH metabolism in multiple biological contexts and structural analysis of these proteins reveal similarities in potential substrate binding regions to the binding residues of SLC25A39.

CONCLUSION

These findings have implications for a better understanding of cancer cell metabolism and novel therapeutic targets with respect to GSH metabolism through potential novel functional annotations of genes. The hybrid ML framework proposed here can be applied to other biological function classifications or multi-omics datasets to generate hypotheses in various biological contexts. Code and a tutorial for generating models and predictions in this framework are available at: https://github.com/lkenn012/mGSH_cancerClassifiers .

Characterization of Delta-7 Alkenone Desaturase in Haptophyte Gephyrocapsa huxleyi Through Heterologous Expression in Tisochrysis lutea

The marine haptophyte Gephyrocapsa huxleyi is an ecologically and geochemically important phytoplankton due to its contribution to the global carbon cycle and its ability to biosynthesize certain alkenones. These alkenones are long-chain alkyl ketones with two to four trans-type double bonds. The genes encoding alkenone desaturase in G. huxleyi have not been experimentally characterized so far, partly due to the difficulty of inducing genetic transformation in G. huxleyi. Therefore, we introduced the putative alkenone delta-7 desaturase of G. huxleyi (designated "DesT") to the transformable and alkenone-producing haptophyte Tisochrysis lutea. We found two types of coding sequences for DesT, which are probably derived from the expression products of different alleles, and designated them "DesT-1" and "DesT-2." The ratio of C37:3 to C37:2 methyl alkenone in the DesT-1 transformant was significantly higher than that in the mock strain that expressed only the hygromycin resistance gene, suggesting that DesT-1 was an alkenone delta-7 desaturase in G. huxleyi. In the protein structure, a tunnel where a substrate alkenone penetrates was predicted to be located around the histidine box of DesT, and hydrophilic and hydrophobic amino acids were respectively located at the proximal (near side to the histidine box) and distal ends of the tunnel. This is the first study to conduct experimental characterization of the alkenone metabolism-related gene in G. huxleyi. The heterologous expression system using T. lutea paves the way for further characterization of the alkenone metabolism-related genes in less transformable haptophytes.

Packaging "vegetable oils": Insights into plant lipid droplet proteins

Plant neutral lipids, also known as "vegetable oils", are synthesized within the endoplasmic reticulum (ER) membrane and packaged into subcellular compartments called lipid droplets (LDs) for stable storage in the cytoplasm. The biogenesis, modulation, and degradation of cytoplasmic LDs in plant cells are orchestrated by a variety of proteins localized to the ER, LDs, and peroxisomes. Recent studies of these LD-related proteins have greatly advanced our understanding of LDs not only as steady oil depots in seeds but also as dynamic cell organelles involved in numerous physiological processes in different tissues and developmental stages of plants. In the past 2 decades, technology advances in proteomics, transcriptomics, genome sequencing, cellular imaging and protein structural modeling have markedly expanded the inventory of LD-related proteins, provided unprecedented structural and functional insights into the protein machinery modulating LDs in plant cells, and shed new light on the functions of LDs in nonseed plant tissues as well as in unicellular algae. Here, we review critical advances in revealing new LD proteins in various plant tissues, point out structural and mechanistic insights into key proteins in LD biogenesis and dynamic modulation, and discuss future perspectives on bridging our knowledge gaps in plant LD biology.

Directed Evolution and Unusual Protonation Mechanism of Pyridoxal Radical C-C Coupling Enzymes for the Enantiodivergent Photobiocatalytic Synthesis of Noncanonical Amino Acids

Visible light-driven pyridoxal radical biocatalysis has emerged as a new strategy for the stereoselective synthesis of valuable noncanonical amino acids in a protecting-group-free fashion. In our previously developed dehydroxylative C-C coupling using engineered PLP-dependent tryptophan synthases, an enzyme-controlled unusual α-stereochemistry reversal and pH-controlled enantiopreference were observed. Herein, through high-throughput photobiocatalysis, we evolved a set of stereochemically complementary PLP radical enzymes, allowing the synthesis of both l- and d-amino acids with enhanced enantiocontrol across a broad pH window. These newly engineered l- and d-amino acid synthases permitted the use of a broad range of organoboron substrates, including boronates, trifluoroborates, and boronic acids, with excellent efficiency. Mechanistic studies unveiled unexpected PLP racemase activity with our earlier PLP enzyme variants. This promiscuous racemase activity was abolished in our evolved amino acid synthases, shedding light on the origin of enhanced enantiocontrol. Further mechanistic investigations suggest a switch of proton donor to account for the stereoinvertive formation of d-amino acids, highlighting an unusual stereoinversion mechanism that is rare in conventional two-electron PLP enzymology.

Assessing the role of deep eutectic solvents in Yarrowia lipolytica inhibition

Yarrowia lipolytica has gained recognition as a microorganism with biological relevance and extensive biotechnological applications. Some of its features include a high enzyme secretion capacity and a high cell-density fermentation mode. Hexokinase (YlHxk) is a vital enzyme in Y. lipolytica growth since it catalyzes glucose metabolism through phosphorylation in the glycolytic pathway. Given the potential application of deep eutectic solvents (DES) as novel solvents in biotechnological processes, this study evaluated the influence of eighteen DES on the growth of Y. lipolytica. Furthermore, this work examined the effects of individual ions on the YlHxk enzyme by analyzing its enzymatic tunnel structure, molecule transport, and molecular docking. The results revealed a significant reduction in yeast growth in the presence of most DES compared to the control (medium without DES), with the exception of the [N8881]Cl: hexanoic acid (1:1) DES. The growth varied between 11.95 ± 0.60 and 0.68 ± 0.17 g dry cell weight L-1. According to the enzymatic tunnel analysis, DES components associated with the lowest microbial growth values were transported through tunnel 1. On the other hand, DES components had their pathway facilitated through tunnel 2 ([N8881]+ and hexanoic acid) and showed growth values close to the control. Molecular docking analysis identified a similarity between all the ligands in this tunnel (including substrate and product), presenting binding interactions with the ASN273 amino acid of the YlHxk active site. Combining experimental results with computational tools provided promising insights at the molecular level, while also potentially reducing analysis costs and time, paving the way for similar approaches in broad biocatalytic reactions.

The NADH-dependent flavin reductase ThdF follows an ordered sequential mechanism though crystal structures reveal two FAD molecules in the active site

Two-component flavin-dependent monooxygenases are of great interest as biocatalysts for the production of pharmaceuticals and other relevant molecules, as they catalyze chemically important reactions such as hydroxylation, epoxidation, and halogenation. The monooxygenase components require a separate flavin reductase which provides the necessary reduced flavin cofactor. The tryptophan halogenase Thal from Streptomyces albogriseolus is a well-characterized two-component flavin-dependent halogenase. Thal exhibits some limitations in terms of halogenation efficiency, also caused by unproductive enzyme-substrate complexes with reduced flavin adenine dinucleotide (FAD). Since the reductase components have an important regulatory function for the activity and efficiency of the monooxygenase by controlling the supply of reduced flavin, here, we studied the so far uncharacterized flavin reductase ThdF from the same gene cluster in S. albogriseolus, which potentially cooperates with Thal. A crystal structure of ThdF in complex with both substrates, FAD and NADH, revealed their orientation for hydride transfer. We obtained two further ThdF structures with two FAD molecules bound to the active site, suggesting a ping-pong bi-bi mechanism. In contrast, steady-state enzyme kinetics clearly showed that ThdF catalyzes flavin reduction via an ordered sequential mechanism, with FAD being bound first and FADH2 released last. Compared to related flavin reductases, ThdF has a low kcat and low KM value. The inhibition of ThdF by NAD+ might limit Thal's halogenation activity when the cellular NADH level is low. These results provide first insights into how the efficiency of Thal could be controlled by flavin reduction at the reductase ThdF.

PathInHydro, a Set of Machine Learning Models to Identify Unbinding Pathways of Gas Molecules in [NiFe] Hydrogenases

Machine learning (ML) is a powerful tool for the automated data analysis of molecular dynamics (MD) simulations. Recent studies showed that ML models can be used to identify protein-ligand unbinding pathways and understand the underlying mechanism. To expedite the examination of MD simulations, we constructed PathInHydro, a set of supervised ML models capable of automatically assigning unbinding pathways for the dissociation of gas molecules from [NiFe] hydrogenases, using the unbinding trajectories of CO and H2 fromDesulfovibrio fructosovorans [NiFe] hydrogenase as a training set. [NiFe] hydrogenases are receiving increasing attention in biotechnology due to their high efficiency in the generation of H2, which is considered by many to be the fuel of the future. However, some of these enzymes are sensitive to O2 and CO. Many efforts have been made to rectify this problem and generate air-stable enzymes by introducing mutations that selectively regulate the access of specific gas molecules to the catalytic site. Herein, we showcase the performance of PathInHydro for the identification of unbinding paths in different test sets, including another gas molecule and a different [NiFe] hydrogenase, which demonstrates its feasibility for the trajectory analysis of a diversity of gas molecules along enzymes with mutations and sequence differences. PathInHydro allows the user to skip time-consuming manual analysis and visual inspection, facilitating data analysis for MD simulations of ligand unbinding from [NiFe] hydrogenases. The codes and data sets are available online: https://github.com/FarzinSohraby/PathInHydro.

Dissecting Methionine Oxidation by Hydrogen Peroxide in Proteins: Thermodynamics, Kinetics, and Susceptibility Descriptors

The oxidation of Met residues in proteins is a complex process, where protein-specific structural and dynamical features play a relevant role in determining the reaction kinetics. Aiming to a full-side perspective, we report here a comprehensive characterization of Met oxidation kinetics by hydrogen peroxide in a leptin protein case study. To do that, we estimated the reaction-free energy profile of the Met oxidation via a QM/MM approach, while the kinetics of the formation of the reactive species were calculated using classical molecular dynamics (MD) simulations. Our data, validated against the available experimental data on the Met oxidation in this protein, indicated that the protein's local and global motion represent the primary discriminating factor among residues' oxidation rates. Moreover, assuming that the free energy profile is independent of the specific protein system, the different reactivities of Met residues within five proteins (hGCSF, IL-1ra, leptin, somatotropin, and RNase) were qualitatively analyzed in terms of well-known structural/dynamic features, which can affect the kinetics of the whole process. The comprehensive analysis of the reaction thermodynamics and kinetics fingerprint enabled the identification of additional descriptors, helpful in assessing the susceptibility of protein-bound Met residues to oxidation.

Semi-rational design of an aromatic dioxygenase by substrate tunnel redirection

Lignin valorization is crucial for achieving economic and sustainable biorefinery processes. However, the enzyme substrate preferences involved in lignin degradation remain poorly understood, and low activity toward specific substrates presents a significant challenge to the efficient utilization of lignin. In this study, we investigated the substrate promiscuity of ThAdo, a key enzyme involved in lignin valorization. Pre-reaction state analysis revealed that a hydrogen bond network is critical in determining substrate selectivity. By performing targeted saturation mutagenesis on residues surrounding the substrate tunnels, we identified the Y205W and Y205Q mutants, which demonstrated 0.73-fold and 0.72-fold enhancements in activity, respectively. Structural analysis indicated that the redirection of the original substrate tunnel may be responsible for the improved activity. Our study provides essential insights into the substrate preference mechanisms of lignin degrading enzymes and suggests that this tunnel-redirection strategy can be extended to other promiscuous enzymes.

Structural basis of phosphate export by human XPR1

Phosphorus in crucial for all living organisms. In vertebrate, cellular phosphate homeostasis is partly controlled by XPR1, a poorly characterized inositol pyrophosphate-dependent phosphate exporter. Here, we report the cryo-EM structure of human XPR1, which forms a loose dimer with 10 transmembrane helices (TM) in each protomer. The structure consists of a scaffold domain (TM1, TM3-4) and a core domain (TM2, TM5-10) structurally related to ion-translocating rhodopsins. Bound phosphate is observed in a tunnel within the core domain at a narrow point that separates the tunnel into intracellular and extracellular vestibules. This site contains a cluster of basic residues that coordinate phosphate and a conserved W573 essential for export function. Loss of inositol pyrophosphate binding is accompanied by structural movements in TM9 and the W573 sidechain, closing the extracellular vestibule and blocking phosphate export. These findings provide insight into XPR1 mechanism and pave the way for further in-depth XPR1 studies.

Water Migration through Enzyme Tunnels Is Sensitive to the Choice of Explicit Water Model

The utilization of tunnels and water transport within enzymes is crucial for their catalytic function as water molecules can stabilize bound substrates and help with unbinding processes of products and inhibitors. Since the choice of water models for molecular dynamics simulations was shown to determine the accuracy of various calculated properties of the bulk solvent and solvated proteins, we have investigated if and to what extent water transport through the enzyme tunnels depends on the selection of the water model. Here, we focused on simulating enzymes with various well-defined tunnel geometries. In a systematic investigation using haloalkane dehalogenase as a model system, we focused on the well-established TIP3P, OPC, and TIP4P-Ew water models to explore their impact on the use of tunnels for water molecule transport. The TIP3P water model showed significantly faster migration, resulting in the transport of approximately 2.5 times more water molecules compared to that of the OPC and 1.7 times greater than that of the TIP4P-Ew. Finally, the transport was 1.4-fold more pronounced in TIP4P-Ew than in OPC. The increase in migration of TIP3P water molecules was mainly due to faster transit times through dehalogenase tunnels. We observed similar behavior in two different enzymes with buried active sites and different tunnel network topologies, i.e., alditol oxidase and cytochrome P450, indicating that our findings are likely not restricted to a particular enzyme family. Overall, this study showcases the critical importance of water models in comprehending the use of enzyme tunnels for small molecule transport. Given the significant role of water availability in various stages of the catalytic cycle and the solvation of substrates, products, and drugs, choosing an appropriate water model may be crucial for accurate simulations of complex enzymatic reactions, rational enzyme design, and predicting drug residence times.

Product release and substrate entry of aldehyde deformylating oxygenase revealed by molecular dynamics simulations

Cyanobacteria can produce alkanes equivalent to diesel fuels through a two-step enzymatic process involving acyl-(acyl carrier protein) reductase (AAR) and aldehyde deformylating oxygenase (ADO), providing a potential renewable biofuel source. AAR binds to ADO for efficient delivery of an aldehyde substrate and they have been proposed to dissociate when the alkane product is released from the same site as the substrate entrance of ADO. However, the dynamics of the substrate and product in ADO during substrate entry and product release are poorly understood. Here, we performed molecular dynamics (MD) simulations of ADO in the presence of substrate or product. We found that while the aldehyde substrate remains close to the active center of ADO before catalysis, the alkane product can dynamically rotate within the hydrophobic tunnel inside ADO toward the product exit after catalysis. Furthermore, the parallel cascade selection (PaCS)-MD simulations of ADO and the AAR/ADO complex identified the locations of the substrate entrance and the multiple exits for product release on ADO. Strikingly, the PaCS-MD simulations revealed that the alkane product can be released from the exit different from the substrate entrance without dissociation of AAR. Based on these results, we propose a reaction model for efficient alkane production by the AAR/ADO complex in which aldehydes and alkanes are synthesized simultaneously while AAR and ADO remain bound, and the aldehyde substrate can be delivered to ADO immediately after alkane release. Our study will be useful in improving the efficiency of bioalkane production using AAR and ADO.

Characterization and functional analysis of type III polyketide synthases in Selaginella moellendorffii

MAIN CONCLUSION

The evolutionary conservation of type III polyketide synthases (PKS) in Selaginella has been elucidated, and the critical amino acid residues of the anther-specific chalcone synthase-like enzyme (SmASCL) have been identified. Selaginella species are the oldest known vascular plants and a valuable resource for the study of metabolic evolution in land plants. Polyketides, especially flavonoids and sporopollenin precursors, are essential prerequisites for plant land colonization. Although type III polyketide synthases (PKS) are widely studied in seed plants, the related enzymes in Selaginella remain poorly characterized. Here, eight type III PKSs were identified in the Selaginella moellendorffii genome and classified into three clusters. Two PKSs were selected for further research based on their phylogenetic relationships and protein sequence similarity. Functional studies revealed that they were chalcone synthase (SmCHS) and anther-specific CHS-like enzyme (SmASCL). These enzymes are involved in the biosynthesis of flavonoids and sporopollenin, respectively. Their sequence information and enzymatic activity are similar to the orthologs in other plants. Phylogenetic analysis revealed that the ASCL and CHS enzymes were separated into two clades from the Bryophyta. These results suggest that CHS and ASCL emerged in the first land plants and then remained conserved during plant evolution. To study the structural basis of the enzymatic function of SmASCL, a series of mutants were constructed. The number of condensation reactions catalyzed by the P210L/Y211D and I200V/G201T double mutants exceeds that of the wild-type enzyme. Our study provides insight into the characteristics and functions of type III PKSs in S. moellendorffii. It also offers clues for a deeper understanding of the relationship between active sites and the enzymatic function of ASCLs.