Structural basis for ion selectivity in potassium-selective channelrhodopsins.

KCR channelrhodopsins (K+-selective light-gated ion channels) have received attention as potential inhibitory optogenetic tools but more broadly pose a fundamental mystery regarding how their K+ selectivity is achieved. Here, we present 2.5-2.7 Å cryo-electron microscopy structures of HcKCR1 and HcKCR2 and of a structure-guided mutant with enhanced K+ selectivity. Structural, electrophysiological, computational, spectroscopic, and biochemical analyses reveal a distinctive mechanism for K+ selectivity; rather than forming the symmetrical filter of canonical K+ channels achieving both selectivity and dehydration, instead, three extracellular-vestibule residues within each monomer form a flexible asymmetric selectivity gate, while a distinct dehydration pathway extends intracellularly. Structural comparisons reveal a retinal-binding pocket that induces retinal rotation (accounting for HcKCR1/HcKCR2 spectral differences), and design of corresponding KCR variants with increased K+ selectivity (KALI-1/KALI-2) provides key advantages for optogenetic inhibition in vitro and in vivo. Thus, discovery of a mechanism for ion-channel K+ selectivity also provides a framework for next-generation optogenetics.

Mechanism of Substrate Translocation in an Alternating Access Transporter.

Transporters shuttle molecules across cell membranes by alternating among distinct conformational states. Fundamental questions remain about how transporters transition between states and how such structural rearrangements regulate substrate translocation. Here, we capture the translocation process by crystallography and unguided molecular dynamics simulations, providing an atomic-level description of alternating access transport. Simulations of a SWEET-family transporter initiated from an outward-open, glucose-bound structure reported here spontaneously adopt occluded and inward-open conformations. Strikingly, these conformations match crystal structures, including our inward-open structure. Mutagenesis experiments further validate simulation predictions. Our results reveal that state transitions are driven by favorable interactions formed upon closure of extracellular and intracellular "gates" and by an unfavorable transmembrane helix configuration when both gates are closed. This mechanism leads to tight allosteric coupling between gates, preventing them from opening simultaneously. Interestingly, the substrate appears to take a "free ride" across the membrane without causing major structural rearrangements in the transporter.

Cryo-EM reveals mechanistic insights into lipid-facilitated polyamine export by human ATP13A2.

The cytoplasmic polyamine maintains cellular homeostasis by chelating toxic metal cations, regulating transcriptional activity, and protecting DNA. ATP13A2 was identified as a lysosomal polyamine exporter responsible for polyamine release into the cytosol, and its dysfunction is associated with Alzheimer's disease and other neural degradation diseases. ATP13A2 belongs to the P5 subfamily of the P-type ATPase family, but its mechanisms remain unknown. Here, we report the cryoelectron microscopy (cryo-EM) structures of human ATP13A2 under four different conditions, revealing the structural coupling between the polyamine binding and the dephosphorylation. Polyamine is bound at the luminal tunnel and recognized through numerous electrostatic and π-cation interactions, explaining its broad specificity. The unique N-terminal domain is anchored to the lipid membrane to stabilize the E2P conformation, thereby accelerating the E1P-to-E2P transition. These findings reveal the distinct mechanism of P5B ATPases, thereby paving the way for neuroprotective therapy by activating ATP13A2.

A Zn-dependent structural transition of SOD1 modulates its ability to undergo phase separation.

The misfolding and mutation of Cu/Zn superoxide dismutase (SOD1) is commonly associated with amyotrophic lateral sclerosis (ALS). SOD1 can accumulate within stress granules (SGs), a type of membraneless organelle, which is believed to form via liquid-liquid phase separation (LLPS). Using wild-type, metal-deficient, and different ALS disease mutants of SOD1 and computer simulations, we report here that the absence of Zn leads to structural disorder within two loop regions of SOD1, triggering SOD1 LLPS and amyloid formation. The addition of exogenous Zn to either metal-free SOD1 or to the severe ALS mutation I113T leads to the stabilization of the loops and impairs SOD1 LLPS and aggregation. Moreover, partial Zn-mediated inhibition of LLPS was observed for another severe ALS mutant, G85R, which shows perturbed Zn-binding. By contrast, the ALS mutant G37R, which shows reduced Cu-binding, does not undergo LLPS. In addition, SOD1 condensates induced by Zn-depletion exhibit greater cellular toxicity than aggregates formed by prolonged incubation under aggregating conditions. Overall, our work establishes a role for Zn-dependent modulation of SOD1 conformation and LLPS properties that may contribute to amyloid formation.

Specific residues and conformational plasticity define the substrate specificity of short-chain dehydrogenases/reductases.

Short-chain dehydrogenases/reductases (SDRs) are one of the most prevalent enzyme families distributed among the sequenced microorganisms. Despite the presence of a conserved catalytic tetrad and high structural similarity, these enzymes exhibit different substrate specificities. The insufficient knowledge regarding the amino acids underlying substrate specificity hinders the understanding of the SDRs' roles in diverse and significant biological processes. Here, we performed bioinformatic analysis, molecular modeling, and mutagenesis studies to identify the key residues that regulate the substrate specificities of two homologous microbial SDRs (i.e., DesE and KduD). Further, we investigated the impact of altering the physicochemical properties of these amino acids on enzyme activity. Interestingly, molecular dynamics simulations also suggest a critical role of enzyme conformational flexibility in substrate recognition and catalysis. Overall, our findings improve the understanding of microbial SDR substrate specificity and shed light on future rational design of more efficient and effective biocatalysts.

High-Strength Polycrystalline Covalent Organic Framework with Abnormal Thermal Transport Insensitive to Grain Boundary.

Grain boundaries (GBs) in two-dimensional (2D) covalent organic frameworks (COFs) unavoidably form during the fabrication process, playing pivotal roles in the physical characteristics of COFs. Herein, molecular dynamics simulations were employed to elucidate the fracture failure and thermal transport mechanisms of polycrystalline COFs (p-COFs). The results revealed that the tilt angle of GBs significantly influences out-of-plane wrinkles and residual stress in monolayer p-COFs. The tensile strength of p-COFs can be enhanced and weakened with the tilt angle, which exhibits an inverse relationship with the defect density. The crack always originates from weaker heptagon rings during uniaxial tension. Notably, the thermal transport in p-COFs is insensitive to the GBs due to the variation of minor polymer chain length at defects, which is abnormal for other 2D crystalline materials. This study contributes insights into the impact of GBs in p-COFs and offers theoretical guidance for structural design and practical applications of advanced COFs.

Time-resolved fluorescence of tryptophan characterizes membrane perturbation by cyclic lipopeptides.

Viscosin is a membrane-permeabilizing, cyclic lipopeptide (CLiP) produced by Pseudomonas species. Here, we have studied four synthetic analogs (L1W, V4W, L5W, L7W), each with one leucine (Leu; L) or valine residue exchanged for tryptophan (Trp; W) by means of time-resolved fluorescence spectroscopy of Trp. To this end, we recorded the average fluorescence lifetime, rotational correlation time and limiting anisotropy, dipolar relaxation time and limiting extent of relaxation, rate constant of acrylamide quenching, effect of H2O-D2O exchange, and time-resolved halfwidth of the spectrum in the absence and presence of POPC liposomes. Structure, localization, and hydration of the peptides were described by molecular dynamics simulations. The combination of the parameters provides a good description of the molecular environments of the Trp positions and the behavior of viscosin as a whole. Of particular value for characterizing the impact of viscosin on the membrane is the dipolar relaxation of Trp4 in V4W, which is deeply embedded in the hydrophobic core. The limiting relaxation level represents the membrane perturbation - unlike typical membrane probes - at the site of the perturbant. Fractions of Trp4 relax at different rates; the one not in contact with water upon excitation relaxes via recruitment of a water molecule on the 10 ns time scale. This rate is sensitive to the concerted membrane perturbation by more than one lipopeptide, which appears at high lipopeptide concentration and is assumed a prerequisite for the final formation of a membrane-permeabilizing defect. Trp7 relaxes primarily with respect to neighboring Ser residues. Trp5 flips between a membrane-inserted and surface-exposed orientation.

Different conformational dynamics of SNARE protein Ykt6 among yeast and mammals.

Ykt6 is one of the most conserved SNARE (N-ethylmaleimide-sensitive factor attachment protein receptor) proteins involved in multiple intracellular membrane trafficking processes. The membrane-anchoring function of Ykt6 has been elucidated to result from its conformational transition from a closed state to an open state. Two ways of regulating the conformational transition were proposed: the C-terminal lipidation and the phosphorylation at the SNARE core. Despite many aspects of common properties, Ykt6 displays differential cellular localizations and functional behaviors in different species, such as yeast, mammals, and worms. The structure-function relationship underlying these differences remains elusive. Here, we combined biochemical characterization, single-molecule FRET measurement, and molecular dynamics simulation to compare the conformational dynamics of yeast and rat Ykt6. Compared to rat Ykt6 (rYkt6), yeast Ykt6 (yYkt6) has more open conformations and could not bind dodecylphosphocholine that inhibits rYkt6 in the closed state. A point mutation T46L/Q57A was shown to be able to convert yYkt6 to a more closed and dodecylphosphocholine-bound state, where Leu46 contributes key hydrophobic interactions for the closed state. We also demonstrated that the phospho-mutation S174D could shift the conformation of rYkt6 to a more open state, but the corresponding mutation S176D in yYkt6 leads to a slightly more closed conformation. These observations shed light on the regulatory mechanism underlying the variations of Ykt6 functions across species.

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

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

Structural and thermodynamic properties of conserved water molecules in Mpro native: A combined approach by MD simulation and Grid Inhomogeneous Solvation Theory.

The new viral strains of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) are continuously rising, becoming more virulent, and transmissible. Therefore, the development of new antiviral drugs is essential. Due to its significant role in the viral life cycle of SARS-CoV-2, the main protease (Mpro) enzyme is a leading target for antiviral drug design. The Mpro monomer consists of domain DI, DII, and DI-DII interface. Twenty-one conserved water molecules (W4-W24) are occupied at these domains according to multiple crystal structure analyses. The crystal and MD structures reveal the presence of eight conserved water sites in domain DI, DII and remaining in the DI-DII interface. Grid-based inhomogeneous fluid solvation theory (GIST) was employed on MD structures of Mpro native to predict structural and thermodynamic properties of each conserved water site for focusing to identify the specific conserved water molecules that can easily be displaced by proposed ligands. Finally, MD water W13 is emerged as a promising candidate for water mimic drug design due to its low mean interaction energy, loose binding character with the protein, and its involvement in a water-mediated H-bond with catalytic His41 via the interaction Thr25(OG)---W13---W---His41(NE2). In this context, water occupancy, relative interaction energy, entropy, and topologies of W13 are thermodynamically acceptable for the water displacement method. Therefore, the strategic use of W13's geometrical position in the DI domain may be implemented for drug discovery against COVID disease by designing new ligands with appropriately oriented chemical groups to mimic its structural, electronic, and thermodynamic properties.

In silico and structure-based evaluation of deleterious mutations identified in human Chk1, Chk2, and Wee1 protein kinase.

Checkpoint kinases Chk1, Chk2, Wee1 are playing a key role in DNA damage response and genomic integrity. Cancer-associated mutations identified in human Chk1, Chk2, and Wee1 were retrieved to understand the function associated with the mutation and also alterations in the folding pattern. Therefore, an attempt has been made to identify deleterious effect of variants using in silico and structure-based approach. Variants of uncertain significance for Chk1, Chk2, and Wee1 were retrieved from different databases and four prediction servers were employed to predict pathogenicity of mutations. Further, Interpro, I-Mutant 3.0, Consurf, TM-align, and have (y)our protein explained were used for comprehensive study of the deleterious effects of variants. The sequences of Chk1, Chk2, and Wee1 were analyzed using Clustal Omega, and the three-dimensional structures of the proteins were aligned using TM-align. The molecular dynamics simulations were performed to explore the differences in folding pattern between Chk1, Chk2, Wee1 wild-type, and mutant protein and also to evaluate the structural integrity. Thirty-six variants in Chk1, 250 Variants in Chk2, and 29 in Wee1 were categorized as pathogenic using in silico prediction tools. Furthermore, 25 mutations in Chk1, 189 in Chk2, and 14 in Wee1 were highly conserved, possessing deleterious effect and also influencing the protein structure and function. These identified mutations may provide underlying genetic intricacies to serve as potential targets for therapeutic inventions and clinical management.

In Silico Prediction of Potential Inhibitors for Targeting RNA CAG Repeats via Molecular Docking and Dynamics Simulation: A Drug Discovery Approach.

Spinocerebellar ataxia (SCA) is a rare neurological illness inherited dominantly that causes severe impairment and premature mortality. While each rare disease may affect individuals infrequently, collectively they pose a significant healthcare challenge. It is mainly carried out due to the expansion of RNA triplet (CAG) repeats, although missense or point mutations can also be induced. Unfortunately, there is no cure; only symptomatic treatments are available. To date, SCA has about 48 subtypes, the most common of these being SCA 1, 2, 3, 6, 7, 12, and 17 having CAG repeats. Using molecular docking and molecular dynamics (MD) simulation, this study seeks to investigate effective natural herbal neuroprotective compounds against CAG repeats, which are therapeutically significant in treating SCA. Initially, virtual screening followed by molecular docking was used to estimate the binding affinity of neuroprotective natural compounds toward CAG repeats. The compound with the highest binding affinity, somniferine, was then chosen for MD simulation. The structural stability, interaction mechanism, and conformational dynamics of CAG repeats and somniferine were investigated via MD simulation. The MD study revealed that during the simulation period, the interaction between CAG repeats and somniferine stabilizes and results in fewer conformational variations. This in silico study suggests that Somniferine can be used as a therapeutic medication against RNA CAG repeats in SCA.

Concurrent inhibition of IR, ITGB1, and CD36 perturbated the interconnected network of energy metabolism and epithelial-to-mesenchymal transition in breast cancer cells.

Altered energy metabolism is an emerging hallmark of cancer and plays a pivotal in cell survival, proliferation, and biosynthesis. In a rapidly proliferating cancer, energy metabolism acts in synergism with epithelial-to-mesenchymal transition (EMT), enabling cancer stemness, dissemination, and metastasis. In this study, an interconnected functional network governing energy metabolism and EMT signaling pathways was targeted through the concurrent inhibition of IR, ITGB1, and CD36 activity. A novel multicomponent MD simulation approach was employed to portray the simultaneous inhibition of IR, ITGB1, and CD36 by a 2:1 combination of Pimozide and Ponatinib. Further, in-vitro studies revealed the synergistic anticancer efficacy of drugs against monolayer as well as tumor spheroids of breast cancer cell lines (MCF-7 and MDA-MB-231). In addition, the combination therapy exerted approximately 40% of the apoptotic population and more than 1.5- to 3-fold reduction in the expression of ITGB1, IR, p-IR, IRS-1, and p-AKT in MCF-7 and MDA-MB-231 cell lines. Moreover, the reduction in fatty acid uptake, lipid droplet accumulation, cancer stemness, and migration properties were also observed. Thus, targeting IR, ITGB1, and CD36 in the interconnected network with the combination of Pimozide and Ponatinib represents a promising therapeutic approach for breast cancer.

Validating Fractalkine receptor as a target and identifying candidates for drug discovery against type 2 diabetes.

Type 2 diabetes mellitus (T2DM) is one of the most common chronic diseases employing abnormal levels of insulin. Enhancing the insulin production is greatly aided by the regulatory mechanisms of the Fractalkine receptor (CX3CR1) system in islet ß-cell function. However, elements including a high-fat diet, obesity, and ageing negatively impact the expression of CX3CR1 in islets. CX3CL1/CX3CR1 receptor-ligand complex is now recognized as a novel therapeutic target. It suggests that T2DM-related ß-cell dysfunction may result from lower amount of these proteins. We analyzed the differential expression of CX3CR1 gene samples taken from persons with T2DM using data obtained from the Gene Expression Omnibus database. Homology modeling enabled us to generate the three-dimensional structure of CX3CR1 and a possible binding pocket. The optimized CX3CR1 structure was subjected to rigorous screening against a massive library of 693 million drug-like molecules from the ZINC15 database. This screening process led to the identification of three compounds with strong binding affinity at the identified binding pocket of CX3CR1. To further evaluate the potential of these compounds, molecular dynamics simulations were conducted over a 50 ns time scale to assess the stability of the protein-ligand complexes. These simulations revealed that ZINC000032506419 emerged as the most promising drug-like compound among the three potent molecules. The discovery of ZINC000032506419 holds exciting promise as a potential therapeutic agent for T2D and other related metabolic disorders. These findings pave the way for the development of effective medications to address the complexities of T2DM and its associated metabolic diseases.

In Silico and In Vitro Study of Isoquercitrin against Kidney Cancer and Inflammation by Triggering Potential Gene Targets.

Kidney cancer has emerged as a major medical problem in recent times. Multiple compounds are used to treat kidney cancer by triggering cancer-causing gene targets. For instance, isoquercitrin (quercetin-3-O-ß-d-glucopyranoside) is frequently present in fruits, vegetables, medicinal herbs, and foods and drinks made from plants. Our previous study predicted using protein-protein interaction (PPI) and molecular docking analysis that the isoquercitrin compound can control kidney cancer and inflammation by triggering potential gene targets of IGF1R, PIK3CA, IL6, and PTGS2. So, the present study is about further in silico and in vitro validation. We performed molecular dynamic (MD) simulation, gene ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, cytotoxicity assay, and RT-PCR and qRT-PCR validation. According to the MD simulation (250 ns), we found that IGF1R, PIK3CA, and PTGS2, except for IL6 gene targets, show stable binding energy with a stable complex with isoquercitrin. We also performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of the final targets to determine their regulatory functions and signaling pathways. Furthermore, we checked the cytotoxicity effect of isoquercitrin (IQ) and found that 5 µg/mL and 10 µg/mL doses showed higher cell viability in a normal kidney cell line (HEK 293) and also inversely showed an inhibition of cell growth at 35% and 45%, respectively, in the kidney cancer cell line (A498). Lastly, the RT-PCR and qRT-PCR findings showed a significant decrease in PTGS2, PIK3CA, and IGF1R gene expression, except for IL6 expression, following dose-dependent treatments with IQ. Thus, we can conclude that isoquercitrin inhibits the expression of PTGS2, PIK3CA, and IGF1R gene targets, which in turn controls kidney cancer and inflammation.

Structural basis of pri-let-7 recognition by human pseudouridine synthase TruB1.

Let-7 was one of the first microRNAs (miRNAs) to be discovered and its expression promotes differentiation during development and function as tumor suppressors in various cancers. The maturation process of let-7 miRNA is tightly regulated by multiple RNA-binding proteins. For example, LIN28 binds to the terminal loops of the precursors of let-7 family and block their processing into mature miRNAs. Trim25 promotes the uridylation-mediated degradation of pre-let-7 modified by LIN28/TUT4. Recently, human pseudouridine synthase TruB1 has been reported to facilitate let-7 maturation by directly binding to pri-let-7 and recruiting Drosha-DGCR8 microprocessor. Through biochemical assay and structural investigation, we show that human TruB1 binds specifically the terminal loop of pri-let-7a1 at nucleotides 31-41, which folds as a small stem-loop architecture. Although TruB1 recognizes the terminal loop of pri-let-7a1 in a way similar to how E. coli TruB interacts with tRNA, a conserved KRKK motif in human and other higher eukaryotes adds an extra binding interface and strengthens the recognition of TruB1 for pri-let-7a1 through electrostatic interactions. These findings reveal the structural basis of TruB1-pri-let-7 interaction which may assists the elucidation of precise role of TruB1 in biogenesis of let-7.

Structural insights in to the atypical type-I ABC Glucose-6-phosphate importer VCA0625-27 of Vibrio cholerae.

Sugar phosphates are potential sources of carbon and phosphate for bacteria. Despite that the process of internalization of Glucose-6-Phosphate (G6P) through plasma membrane remained elusive in several bacteria. VCA0625-27, made of periplasmic ligand binding protein (PLBP) VCA0625, an atypical monomeric permease VCA0626, and a cytosolic ATPase VCA0627, recently emerged as hexose-6-phosphate uptake system of Vibrio cholerae. Here we report high resolution crystal structure of VCA0625 in G6P bound state that largely resembles AfuA of Actinobacillus pleuropneumoniae. MD simulations on VCA0625 in apo and G6P bound states unraveled an 'open to close' and swinging bi-lobal motions, which are diminished upon G6P binding. Mutagenesis followed by biochemical assays on VCA0625 underscored that R34 works as gateway to bind G6P. Although VCA0627 binds ATP, it is ATPase deficient in the absence of VCA0625 and VCA0626, which is a signature phenomenon of type-I ABC importer. Further, modeling, docking and systematic sequence analysis allowed us to envisage the existence of similar atypical type-I G6P importer with fused monomeric permease in 27 other gram-negative bacteria.

Probing the range of applicability of structure- and energy-adjusted QM/MM link bonds III: QM/MM MD simulations of solid-state systems at the example of layered carbon structures.

The previously introduced workflow to achieve an energetically and structurally optimized description of frontier bonds in quantum mechanical/molecular mechanics (QM/MM)-type applications was extended into the regime of computational material sciences at the example of a layered carbon model systems. Optimized QM/MM link bond parameters at HSEsol/6-311G(d,p) and self-consistent density functional tight binding (SCC-DFTB) were derived for graphitic systems, enabling detailed investigation of specific structure motifs occurring in graphene-derived structures v i a $$ via $$ quantum-chemical calculations. Exemplary molecular dynamics (MD) simulations in the isochoric-isothermic (NVT) ensemble were carried out to study the intercalation of lithium and the properties of the Stone-Thrower-Wales defect. The diffusivity of lithium as well as hydrogen and proton adsorption on a defective graphene surface served as additional example. The results of the QM/MM MD simulations provide detailed insight into the applicability of the employed link-bond strategy when studying intercalation and adsorption properties of graphitic materials.

Characterizing the stabilization effects of stabilizers in protein-protein systems with end-point binding free energy calculations.

Drug design targeting protein-protein interactions (PPIs) associated with the development of diseases has been one of the most important therapeutic strategies. Besides interrupting the PPIs with PPI inhibitors/blockers, increasing evidence shows that stabilizing the interaction between two interacting proteins may also benefit the therapy, such as the development of various types of molecular glues/stabilizers that mostly work by stabilizing the two interacting proteins to regulate the downstream biological effects. However, characterizing the stabilization effect of a stabilizer is usually hard or too complicated for traditional experiments since it involves ternary interactions [protein-protein-stabilizer (PPS) interaction]. Thus, developing reliable computational strategies will facilitate the discovery/design of molecular glues or PPI stabilizers. Here, by fully analyzing the energetic features of the binary interactions in the PPS ternary complex, we systematically investigated the performance of molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) and molecular mechanics generalized Born surface area (MM/GBSA) methods on characterizing the stabilization effects of stabilizers in 14-3-3 systems. The results show that both MM/PBSA and MM/GBSA are powerful tools in distinguishing the stabilizers from the decoys (with area under the curves of 0.90-0.93 for all tested cases) and are reasonable for ranking protein-peptide interactions in the presence or absence of stabilizers as well (with the average Pearson correlation coefficient of ~0.6 at a relatively high dielectric constant for both methods). Moreover, to give a detailed picture of the stabilization effects, the stabilization mechanism is also analyzed from the structural and energetic points of view for individual systems containing strong or weak stabilizers. This study demonstrates a potential strategy to accelerate the discovery of PPI stabilizers.

Structure-Guided Evolution of Cyclohexanone Monooxygenase Toward Bulky Omeprazole Sulfide: Substrate Migration and Stereoselectivity Inversion.

Structure-guided engineering of a CHMO from Amycolatopsis methanolica (AmCHMO) was performed for asymmetric sulfoxidation activity and stereoselectivity toward omeprazole sulfide. Initially, combinatorial active-site saturation test (CASTing) and iteratively saturation mutagenesis (ISM) were performed on 5 residues at the "bottleneck" of substrate tunnel, and MT3 was successfully obtained with a specific activity of 46.19 U/g and R-stereoselectivity of 99 % toward OPS. Then, 4 key mutations affecting the stereoselectivity were identified through multiple rounds of ISM on residues at the substrate binding pocket region, resulting MT8 with an inversed stereoselectivity from 99 % (R) to 97 % (S). MT8 has a greatly compromised specific activity of 0.08 U/g. By introducing additional beneficial mutations, MT11 was constructed with significantly increased specific activity of 2.29 U/g and stereoselectivity of 97 % (S). Enlarged substrate tunnel is critical to the expanded substrate spectrum of AmCHMO, while reshaping of substrate binding pocket is important for stereoselective inversion. Based on MD simulation, pre-reaction states of MT3-OPSproR, MT8-OPSproS, and MT11-OPSproS were calculated to be 45.56 %, 17.94 %, and 28.65 % respectively, which further confirm the experimental data on activity and stereoselectivity. Our results pave the way for engineering distinct activity and stereoselectivity of BVMOs toward bulky prazole thioethers.