Biophysical characterization of the ETV6 PNT domain polymerization interfaces.

ETV6 is an E26 transformation specific family transcriptional repressor that self-associates by its PNT domain to facilitate cooperative DNA binding. Chromosomal translocations frequently generate constitutively active oncoproteins with the ETV6 PNT domain fused to the kinase domain of one of many protein tyrosine kinases. Although an attractive target for therapeutic intervention, the propensity of the ETV6 PNT domain to polymerize via the tight head-to-tail association of two relatively flat interfaces makes it challenging to identify suitable small molecule inhibitors of this protein-protein interaction. Herein, we provide a comprehensive biophysical characterization of the ETV6 PNT domain interaction interfaces to aid future drug discovery efforts and help define the mechanisms by which its self-association mediates transcriptional repression. Using NMR spectroscopy, X-ray crystallography, and molecular dynamics simulations, along with amide hydrogen exchange measurements, we demonstrate that monomeric PNT domain variants adopt very stable helical bundle folds that do not change in conformation upon self-association into heterodimer models of the ETV6 polymer. Surface plasmon resonance-monitored alanine scanning mutagenesis studies identified hot spot regions within the self-association interfaces. These regions include both central hydrophobic residues and flanking salt-bridging residues. Collectively, these studies indicate that small molecules targeted to these hydrophobic or charged regions within the relatively rigid interfaces could potentially serve as orthosteric inhibitors of ETV6 PNT domain polymerization.

In Silico Identification of Potential Sites for a Plastic-Degrading Enzyme by a Reverse Screening through the Protein Sequence Space and Molecular Dynamics Simulations.

The accumulation of polyethylene terephthalate (PET) seriously harms the environment because of its high resistance to degradation. The recent discovery of the bacteria-secreted biodegradation enzyme, PETase, sheds light on PET recycling; however, the degradation efficiency is far from practical use. Here, in silico alanine scanning mutagenesis (ASM) and site-saturation mutagenesis (SSM) were employed to construct the protein sequence space from binding energy of the PETase-PET interaction to identify the number and position of mutation sites and their appropriate side-chain properties that could improve the PETase-PET interaction. The binding mechanisms of the potential PETase variant were investigated through atomistic molecular dynamics simulations. The results show that up to two mutation sites of PETase are preferable for use in protein engineering to enhance the PETase activity, and the proper side chain property depends on the mutation sites. The predicted variants agree well with prior experimental studies. Particularly, the PETase variants with S238C or Q119F could be a potential candidate for improving PETase. Our combination of in silico ASM and SSM could serve as an alternative protocol for protein engineering because of its simplicity and reliability. In addition, our findings could lead to PETase improvement, offering an important contribution towards a sustainable future.

Alanine-scanning mutagenesis of protein mannosyl-transferase from Streptomyces coelicolor reveals strong activity-stability correlation.

In Actinobacteria, protein O-mannosyl transferase (Pmt)-mediated protein O-glycosylation has an important role in cell envelope physiology. In S. coelicolor, defective Pmt leads to increased susceptibility to cell wall-targeting antibiotics, including vancomycin and ß-lactams, and resistance to phage ϕC31. The aim of this study was to gain a deeper understanding of the structure and function of S. coelicolor Pmt. Sequence alignments and structural bioinformatics were used to identify target sites for an alanine-scanning mutagenesis study. Mutant alleles were introduced into pmt-deficient S. coelicolor strains using an integrative plasmid and scored for their ability to complement phage resistance and antibiotic hypersusceptibility phenotypes. Twenty-three highly conserved Pmt residues were each substituted for alanine. Six mutant alleles failed to complement the pmt▬ strains in either assay. Mapping the six corresponding residues onto a homology model of the three-dimensional structure of Pmt, indicated that five are positioned close to the predicted catalytic DE motif. Further mutagenesis to produce more conservative substitutions at these six residues produced Pmts that invariably failed to complement the DT1025 pmt▬ strain, indicating that strict residue conservation was necessary to preserve function. Cell fractionation and Western blotting of strains with the non-complementing pmt alleles revealed undetectable levels of the enzyme in either the membrane fractions or whole cell lysates. Meanwhile for all of the strains that complemented the antibiotic hypersusceptibility and phage resistance phenotypes, Pmt was readily detected in the membrane fraction. These data indicate a tight correlation between the activity of Pmt and its stability or ability to localize to the membrane.

Computational Hot-Spot Analysis of the SARS-CoV-2 Receptor Binding Domain/ACE2 Complex*.

Infection and replication of SARS CoV-2 (the virus that causes COVID-19) requires entry to the interior of host cells. In humans, a protein-protein interaction (PPI) between the SARS CoV-2 receptor-binding domain (RBD) and the extracellular peptidase domain of ACE2 on the surface of cells in the lower respiratory tract is an initial step in the entry pathway. Inhibition of the SARS CoV-2 RBD/ACE2 PPI is currently being evaluated as a target for therapeutic and/or prophylactic intervention. However, relatively little is known about the molecular underpinnings of this complex. Employing multiple computational platforms, we predicted "hot-spot" residues in a positive-control PPI (PMI/MDM2) and the CoV-2 RBD/ACE2 complex. Computational alanine scanning mutagenesis was performed to predict changes in Gibbs' free energy that are associated with mutating residues at the positive control (PMI/MDM2) or SARS RBD/ACE2 binding interface to alanine. Additionally, we used the Adaptive Poisson-Boltzmann Solver to calculate macromolecular electrostatic surfaces at the interface of the positive-control PPI and SARS CoV-2/ACE2 PPI. Finally, a comparative analysis of hot-spot residues for SARS-CoV and SARS-CoV-2, in complex with ACE2, is provided. Collectively, this study illuminates predicted hot-spot residues, and clusters, at the SARS CoV-2 RBD/ACE2 binding interface, potentially guiding the development of reagents capable of disrupting this complex and halting COVID-19.

Alanine Scanning Mutagenesis of the DRYxxI Motif and Intracellular Loop 2 of Human Melanocortin-4 Receptor.

The melanocortin-4 receptor (MC4R) is a member of the G-protein-coupled receptor (GPCR) superfamily, which has been extensively studied in obesity pathogenesis due to its critical role in regulating energy homeostasis. Both the Gs-cAMP and ERK1/2 cascades are known as important intracellular signaling pathways initiated by the MC4R. The DRYxxI motif at the end of transmembrane domain 3 and the intracellular loop 2 (ICL2) are thought to be crucial for receptor function in several GPCRs. To study the functions of this domain in MC4R, we performed alanine-scanning mutagenesis on seventeen residues. We showed that one residue was critical for receptor cell surface expression. Eight residues were important for ligand binding. Mutations of three residues impaired Gs-cAMP signaling without changing the binding properties. Investigation on constitutive activities of all the mutants in the cAMP pathway revealed that six residues were involved in constraining the receptor in inactive states and five residues were important for receptor activation in the absence of an agonist. In addition, mutations of four residues impaired the ligand-stimulated ERK1/2 signaling pathway without affecting the binding properties. We also showed that some mutants were biased to the Gs-cAMP or ERK1/2 signaling pathway. In summary, we demonstrated that the DRYxxI motif and ICL2 were important for MC4R function.

Mutational analysis of a catalytically important loop containing active site and substrate-binding site in Escherichia coli phytase AppA.

A phytase from Escherichia coli, AppA, has been the target of protein engineering to reduce the amount of undigested phosphates from livestock manure by making phosphorous from phytic acid available as a nutrient. To understand the contribution of each amino acid in the active site loop to the AppA activity, alanine and glycine scanning mutagenesis was undertaken. The results of phytase activity assay demonstrated loss of activity by mutations at charged residues within the conserved motif, supporting their importance in catalytic activity. In contrast, both conserved, non-polar residues and non-conserved residues tended to be tolerant to Ala and/or Gly mutations. Correlation analyses of chemical/structural characteristics of each mutation site against mutant activity revealed that the loop residues located closer to the substrate have greater contribution to the activity of AppA. These results may be useful in efficiently engineering AppA to improve its catalytic activity. Abbreviations: AppA: pH 2.5 acid phosphatase; CSU: contacts of structural units; HAPs: histidine acid phosphatases; SASA: solvent accessible surface area; SDS-PAGE: sodium dodecyl sulfate-polyacrylamide gel electrophoresis; SSM: site-saturation mutagenesis; WT: wild type.

Understanding the structural and energetic basis of PD-1 and monoclonal antibodies bound to PD-L1: A molecular modeling perspective.

BACKGROUND: The inhibitors blocking the interaction between programmed cell death protein 1(PD-1) and programmed death-ligand 1(PD-L1) can activate the immune response of T cell and eliminate cancer cells. The crystallographic studies have provided structural insights of the interactive interfaces between PD-L1 and its protein ligands. However, the hotspot residues on PD-L1 as well as structural and energetic basis for different protein ligands still need to be further investigated. METHODS: Molecular modeling methods including molecular dynamics simulation, per-residue free energy decomposition, virtual alanine scanning mutagenesis and residue-residue contact analysis were used to qualitatively and quantitatively analyze the interactions between PD-L1 and different protein ligands. RESULTS: The results of virtual alanine scanning mutagenesis suggest that Y56, Q66, M115, D122, Y123, R125 are the hotspot residues on PD-L1. The residue-residue contact analysis further shows that PD-1 interacts with PD-L1 mainly by F and G strands while monoclonal antibodies like avelumab and BMS-936559 mainly interact with PD-L1 by CDR2 and CDR3 loops of the heavy chain. CONCLUSIONS: A structurally similar ß-hairpin peptide with 13 or 14 residues was extracted from each protein ligand and these ß-hairpin peptides were found tightly binding to the putative hotspot residues on PD-L1. GENERAL SIGNIFICANCE: This study recognizes the hotspot residues on PD-L1 and uncovers the common structural and energetic basis of different protein ligands binding to PD-L1. These results will be valuable for the design of small molecule or peptide inhibitors targeting on PD-L1.

Whole-protein alanine-scanning mutagenesis of allostery: A large percentage of a protein can contribute to mechanism.

Many studies of allosteric mechanisms use limited numbers of mutations to test whether residues play "key" roles. However, if a large percentage of the protein contributes to allosteric function, mutating any residue would have a high probability of modifying allostery. Thus, a predicted mechanism that is dependent on only a few residues could erroneously appear to be supported. We used whole-protein alanine-scanning mutagenesis to determine which amino acid sidechains of human liver pyruvate kinase (hL-PYK; approved symbol PKLR) contribute to regulation by fructose-1,6-bisphosphate (Fru-1,6-BP; activator) and alanine (inhibitor). Each nonalanine/nonglycine residue of hL-PYK was mutated to alanine to generate 431 mutant proteins. Allosteric functions in active proteins were quantified by following substrate affinity over a concentration range of effectors. Results show that different residues contribute to the two allosteric functions. Only a small fraction of mutated residues perturbed inhibition by alanine. In contrast, a large percentage of mutated residues influenced activation by Fru-1,6-BP; inhibition by alanine is not simply the reverse of activation by Fru-1,6-BP. Moreover, the results show that Fru-1,6-BP activation would be extremely difficult to elucidate using a limited number of mutations. Additionally, this large mutational data set will be useful to train and test computational algorithms aiming to predict allosteric mechanisms.

Structure and sequence based functional annotation of Zika virus NS2b protein: Computational insights.

While Zika virus (ZIKV) outbreaks are a growing concern for global health, a deep understanding about the virus is lacking. Here we report a contribution to the basic science on the virus- a detailed computational analysis of the non structural protein NS2b. This protein acts as a cofactor for the NS3 protease (NS3Pro) domain that is important on the viral life cycle, and is an interesting target for drug development. We found that ZIKV NS2b cofactor is highly similar to other virus within the Flavivirus genus, especially to West Nile Virus, suggesting that it is completely necessary for the protease complex activity. Furthermore, the ZIKV NS2b has an important role to the function and stability of the ZIKV NS3 protease domain even when presents a low conservation score. In addition, ZIKV NS2b is mostly rigid, which could imply a non dynamic nature in substrate recognition. Finally, by performing a computational alanine scanning mutagenesis, we found that residues Gly 52 and Asp 83 in the NS2b could be important in substrate recognition.

Discovering key residues of dengue virus NS2b-NS3-protease: New binding sites for antiviral inhibitors design.

The NS2B-NS3 protease is essential for the Dengue Virus (DENV) replication process. This complex constitutes a target for efficient antiviral discovery because a drug could inhibit the viral polyprotein processing. Furthermore, since the protease is highly conserved between the four Dengue virus serotypes, it is probable that a drug would be equally effective against all of them. In this article, a strategy is reported that allowed us to identify influential residues on the function of the Dengue NS2b-NS3 Protease. Moreover, this is a strategy that could be applied to virtually any protein for the search of alternative influential residues, and for non-competitive inhibitor development. First, we incorporated several features derived from computational alanine scanning mutagenesis, sequence, structure conservation, and other structure-based characteristics. Second, these features were used as variables to obtain a multilayer perceptron model to identify defined groups (clusters) of key residues as possible candidate pockets for binding sites of new leads on the DENV protease. The identified residues included: i) amino acids close to the beta sheet-loop-beta sheet known to be important in its closed conformation for NS2b ii) residues close to the active site, iii) several residues evenly spread on the NS2b-NS3 contact surface, and iv) some inner residues most likely related to the overall stability of the protease. In addition, we found concordance on our list of residues with previously identified amino acids part of a highly conserved peptide studied for vaccine development.

Alanine scan of α-conotoxin RegIIA reveals a selective α3ß4 nicotinic acetylcholine receptor antagonist.

Activation of the α3ß4 nicotinic acetylcholine receptor (nAChR) subtype has recently been implicated in the pathophysiology of various conditions, including development and progression of lung cancer and in nicotine addiction. As selective α3ß4 nAChR antagonists, α-conotoxins are valuable tools to evaluate the functional roles of this receptor subtype. We previously reported the discovery of a new α4/7-conotoxin, RegIIA. RegIIA was isolated from Conus regius and inhibits acetylcholine (ACh)-evoked currents mediated by α3ß4, α3ß2, and α7 nAChR subtypes. The current study used alanine scanning mutagenesis to understand the selectivity profile of RegIIA at the α3ß4 nAChR subtype. [N11A] and [N12A] RegIIA analogs exhibited 3-fold more selectivity for the α3ß4 than the α3ß2 nAChR subtype. We also report synthesis of [N11A,N12A]RegIIA, a selective α3ß4 nAChR antagonist (IC50 of 370 nM) that could potentially be used in the treatment of lung cancer and nicotine addiction. Molecular dynamics simulations of RegIIA and [N11A,N12A]RegIIA bound to α3ß4 and α3ß2 suggest that destabilization of toxin contacts with residues at the principal and complementary faces of α3ß2 (α3-Tyr(92), Ser(149), Tyr(189), Cys(192), and Tyr(196); ß2-Trp(57), Arg(81), and Phe(119)) may form the molecular basis for the selectivity shift.

Active site fingerprinting and pharmacophore screening strategies for the identification of dual inhibitors of protein kinase C (ΡΚCß) and poly (ADP-ribose) polymerase-1 (PARP-1).

Current clinical studies have revealed that diabetic complications are multifactorial disorders that target two or more pathways. The majority of drugs in clinical trial target aldose reductase and protein kinase C ([Formula: see text]), while recent studies disclosed a significant role played by poly (ADP-ribose) polymerase-1 (PARP-1). In light of this, the current study was aimed to identify novel dual inhibitors of [Formula: see text] and PARP-1 using a pharmaco-informatics methodology. Pharmacophore-based 3D QSAR models for these two targets were generated using HypoGen and used to screen three commercially available chemical databases to identify dual inhibitors of [Formula: see text] and PARP-1. Overall, 18 hits were obtained from the screening process; the hits were filtered based on their drug-like properties and predicted binding affinities (docking analysis). Important amino acid residues were predicted by developing a fingerprint of the active site using alanine-scanning mutagenesis and molecular dynamics. The stability of the complexes (18 hits with both proteins) and their final binding orientations were investigated using molecular dynamics simulations. Thus, novel hits have been predicted to have good binding affinities for [Formula: see text] and PARP-1 proteins, which could be further investigated for in vitro/in vivo activity.

Glutamine Synthetase Drugability beyond Its Active Site: Exploring Oligomerization Interfaces and Pockets.

BACKGROUND: Glutamine synthetase (GS) is a crucial enzyme to the nitrogen cycle with great commercial and pharmaceutical value. Current inhibitors target the active site, affecting GS activity indiscriminately in all organisms. As the active site is located at the interface between two monomers, the protein-protein interface (PPI) of GSs gains a new role, by providing new targets for enzyme inhibition. Exploring GSs PPI could allow for the development of inhibitors selective for specific organisms. Here we map the PPI of three GSs-human (hsGS), maize (zmGS) and Mycobacterium tuberculosis (mtGS)-and unravel new drugable pockets. METHODS: The PPI binding free energy coming from key residues on three GSs from different organisms were mapped by computational alanine scan mutagenesis, applying a multiple dielectric constant MM-PBSA methodology. The most relevant residues for binding are referred as hot-spots. Drugable pockets on GS were detected with the Fpocket software. RESULTS AND CONCLUSIONS: A total of 23, 19 and 30 hot-spots were identified on hsGS, zmGS and mtGS PPI. Even possessing differences in the hot-spots, hsGS and zmGS PPI are overall very similar. On the other hand, mtGS PPI differs greatly from hsGS and zmGS PPI. A novel drugable pocket was detected on the mtGS PPI. It seems particularly promising for the development of selective anti-tuberculosis drugs given its location on a PPI region that is highly populated with hot-spots and is completely different from the hsGS and zmGS PPIs. Drugs targeting this pockets should be inactive on eukaryotic GS II enzymes.

Molecular modeling reveals binding interface of γ-tubulin with GCP4 and interactions with noscapinoids.

The initiation of microtubule assembly within cells is guided by a cone shaped multi-protein complex, γ-tubulin ring complex (γTuRC) containing γ-tubulin and atleast five other γ-tubulin-complex proteins (GCPs), i.e., GCP2, GCP3, GCP4, GCP5, and GCP6. The rim of γTuRC is a ring of γ-tubulin molecules that interacts, via one of its longitudinal interfaces, with GCP2, GCP3, or GCP4 and, via other interface, with α/ß-tubulin dimers recruited for the microtubule lattice formation. These interactions however, are not well understood in the absence of crystal structure of functional reconstitution of γTuRC subunits. In this study, we elucidate the atomic interactions between γ-tubulin and GCP4 through computational techniques. We simulated two complexes of γ-tubulin-GCP4 complex (we called dimer1 and dimer2) for 25 ns to obtain a stable complex and calculated the ensemble average of binding free energies of -158.82 and -170.19 kcal/mol for dimer1 and -79.53 and -101.50 kcal/mol for dimer2 using MM-PBSA and MM-GBSA methods, respectively. These highly favourable binding free energy values points to very robust interactions between GCP4 and γ-tubulin. From the results of the free-energy decomposition and the computational alanine scanning calculation, we identified the amino acids crucial for the interaction of γ-tubulin with GCP4, called hotspots. Furthermore, in the endeavour to identify chemical leads that might interact at the interface of γ-tubulin-GCP4 complex; we found a class of compounds based on the plant alkaloid, noscapine that binds with high affinity in a cavity close to γ-tubulin-GCP4 interface compared with previously reported compounds. All noscapinoids displayed stable interaction throughout the simulation, however, most robust interaction was observed for bromo-noscapine followed by noscapine and amino-noscapine. This offers a novel chemical scaffold for γ-tubulin binding drugs near γ-tubulin-GCP4 interface.

Structural and functional analysis of transmembrane segment IV of the salt tolerance protein Sod2.

Sod2 is the plasma membrane Na(+)/H(+) exchanger of the fission yeast Schizosaccharomyces pombe. It provides salt tolerance by removing excess intracellular sodium (or lithium) in exchange for protons. We examined the role of amino acid residues of transmembrane segment IV (TM IV) ((126)FPQINFLGSLLIAGCITSTDPVLSALI(152)) in activity by using alanine scanning mutagenesis and examining salt tolerance in sod2-deficient S. pombe. Two amino acids were critical for function. Mutations T144A and V147A resulted in defective proteins that did not confer salt tolerance when reintroduced into S. pombe. Sod2 protein with other alanine mutations in TM IV had little or no effect. T144D and T144K mutant proteins were inactive; however, a T144S protein was functional and provided lithium, but not sodium, tolerance and transport. Analysis of sensitivity to trypsin indicated that the mutations caused a conformational change in the Sod2 protein. We expressed and purified TM IV (amino acids 125-154). NMR analysis yielded a model with two helical regions (amino acids 128-142 and 147-154) separated by an unwound region (amino acids 143-146). Molecular modeling of the entire Sod2 protein suggested that TM IV has a structure similar to that deduced by NMR analysis and an overall structure similar to that of Escherichia coli NhaA. TM IV of Sod2 has similarities to TM V of the Zygosaccharomyces rouxii Na(+)/H(+) exchanger and TM VI of isoform 1 of mammalian Na(+)/H(+) exchanger. TM IV of Sod2 is critical to transport and may be involved in cation binding or conformational changes of the protein.

Identification of the ß-lactamase inhibitor protein-II (BLIP-II) interface residues essential for binding affinity and specificity for class A ß-lactamases.

The interactions between ß-lactamase inhibitory proteins (BLIPs) and ß-lactamases have been used as model systems to understand the principles of affinity and specificity in protein-protein interactions. The most extensively studied tight binding inhibitor, BLIP, has been characterized with respect to amino acid determinants of affinity and specificity for binding ß-lactamases. BLIP-II, however, shares no sequence or structural homology to BLIP and is a femtomolar to picomolar potency inhibitor, and the amino acid determinants of binding affinity and specificity are unknown. In this study, alanine scanning mutagenesis was used in combination with determinations of on and off rates for each mutant to define the contribution of residues on the BLIP-II binding surface to both affinity and specificity toward four ß-lactamases of diverse sequence. The residues making the largest contribution to binding energy are heavily biased toward aromatic amino acids near the center of the binding surface. In addition, substitutions that reduce binding energy do so by increasing off rates without impacting on rates. Also, residues with large contributions to binding energy generally exhibit low temperature factors in the structures of complexes. Finally, with the exception of D206A, BLIP-II alanine substitutions exhibit a similar trend of effect for all ß-lactamases, i.e., a substitution that reduces affinity for one ß-lactamase usually reduces affinity for all ß-lactamases tested.

Solvent-accessible surface area: How well can be applied to hot-spot detection?

A detailed comprehension of protein-based interfaces is essential for the rational drug development. One of the key features of these interfaces is their solvent accessible surface area profile. With that in mind, we tested a group of 12 SASA-based features for their ability to correlate and differentiate hot- and null-spots. These were tested in three different data sets, explicit water MD, implicit water MD, and static PDB structure. We found no discernible improvement with the use of more comprehensive data sets obtained from molecular dynamics. The features tested were shown to be capable of discerning between hot- and null-spots, while presenting low correlations. Residue standardization such as rel SASAi or rel/res SASAi , improved the features as a tool to predict ΔΔGbinding values. A new method using support machine learning algorithms was developed: SBHD (Sasa-Based Hot-spot Detection). This method presents a precision, recall, and F1 score of 0.72, 0.81, and 0.76 for the training set and 0.91, 0.73, and 0.81 for an independent test set.

The interaction of InvF-RNAP is mediated by the chaperone SicA in Salmonella sp: an in silico prediction.

In this work we carried out an in silico analysis to understand the interaction between InvF-SicA and RNAP in the bacterium Salmonella Typhimurium strain LT2. Structural analysis of InvF allowed the identification of three possible potential cavities for interaction with SicA. This interaction could occur with the structural motif known as tetratricopeptide repeat (TPR) 1 and 2 in the two cavities located in the interface of the InvF and α-CTD of RNAP. Indeed, molecular dynamics simulations showed that SicA stabilizes the Helix-turn-Helix DNA-binding motifs, i.e., maintaining their proper conformation, mainly in the DNA Binding Domain (DBD). Finally, to evaluate the role of amino acids that contribute to protein-protein affinity, an alanine scanning mutagenesis approach, indicated that R177 and R181, located in the DBD motif, caused the greatest changes in binding affinity with α-CTD, suggesting a central role in the stabilization of the complex. However, it seems that the N-terminal region also plays a key role in the protein-protein interaction, especially the amino acid R40, since we observed conformational flexibility in this region allowing it to interact with interface residues. We consider that this analysis opens the possibility to validate experimentally the amino acids involved in protein-protein interactions and explore other regulatory complexes where chaperones are involved.

Engineered ß-glycosidase from Hyperthermophilic Sulfolobus solfataricus with Improved Rd-hydrolyzing Activity for Ginsenoside Compound K Production.

Hyperthermophilic Sulfolobus solfataricus ß-glycosidase (SS-ßGly), with higher stability and activity than mesophilic enzymes, has potential for industrial ginsenosides biotransformation. However, its relatively low ginsenoside Rd-hydrolyzing activity limits the production of pharmaceutically active minor ginsenoside compound K (CK). In this study, first, we used molecular docking to predict the key enzyme residues that may hypothetically interact with ginsenoside Rd. Then, based on sequence alignment and alanine scanning mutagenesis approach, key variant sites were identified that might improve the enzyme catalytic efficiency. The enzyme catalytic efficiency (kcat/Km) and substrate affinity (Km) of the N264D variant enzyme for ginsenoside Rd increased by 60% and decreased by 17.9% compared with WT enzyme, respectively, which may be due to a decrease in the binding free energy (∆G) between the variant enzyme and substrate Rd. In addition, Markov state models (MSM) analysis during the whole 1000-ns MD simulations indicated that altering N264 to D made the variant enzyme achieve a more stable SS-ßGly conformational state than the wild-type (WT) enzyme and corresponding Rd complex. Under identical conditions, the relative activities and the CK conversion rates of the N264D enzyme were 1.7 and 1.9 folds higher than those of the WT enzyme. This study identified an excellent hyperthermophilic ß-glycosidase candidate for industrial biotransformation of ginsenosides.

Structural requirement of RARγ agonism through computational aspects.

CONTEXT: RARγ is a therapeutic target for many skin diseases and has potential in cancer treatment. In the current study, we put forward a comprehensive structure-activity relationship study of third and fourth generations of RARγ agonists, addressing multiple crystal structures of RARγ complexes and approved drugs. Adapalene and Trifarotene, through hybrid strategies including protein contacts Atlas analysis, molecular docking, dynamics simulations, MM-GBSA, ASM, and pharmacophore modeling. Our result revealed crucial amino acids Arg267, Ser278, Phe288, Phe230, Met272, Leu271, and Leu268 within the RARγ pocket, as well as pharmacophore features such as two hydrophobic groups, two aromatic rings, and negative ionic features, which are essential for the binding of RARγ agonists. Based on this study, the binding mechanism of RARγ agonists was elucidated, which will be helpful for the rational design of new RARγ agonists for skin diseases and cancer treatment. METHODS: In this study, Schrödinger suite 2021-2 with OPLS_4 force field, Discovery Studio program 3.0, LigandScout 4.3, and PyMOL are utilized in the investigation.