Exploring the Conformational Ensembles of Protein-Protein Complex with Transformer-Based Generative Model.

Protein-protein interactions are the basis of many protein functions, and understanding the contact and conformational changes of protein-protein interactions is crucial for linking the protein structure to biological function. Although difficult to detect experimentally, molecular dynamics (MD) simulations are widely used to study the conformational ensembles and dynamics of protein-protein complexes, but there are significant limitations in sampling efficiency and computational costs. In this study, a generative neural network was trained on protein-protein complex conformations obtained from molecular simulations to directly generate novel conformations with physical realism. We demonstrated the use of a deep learning model based on the transformer architecture to explore the conformational ensembles of protein-protein complexes through MD simulations. The results showed that the learned latent space can be used to generate unsampled conformations of protein-protein complexes for obtaining new conformations complementing pre-existing ones, which can be used as an exploratory tool for the analysis and enhancement of molecular simulations of protein-protein complexes.

Leveraging the Fragment Molecular Orbital and MM-GBSA Methods in Virtual Screening for the Discovery of Novel Non-Covalent Inhibitors Targeting the TEAD Lipid Binding Pocket.

The Hippo pathway controls organ size and homeostasis and is linked to numerous diseases, including cancer. The transcriptional enhanced associate domain (TEAD) family of transcription factors acts as a receptor for downstream effectors, namely yes-associated protein (YAP) and transcriptional co-activator with PDZ-binding motif (TAZ), which binds to various transcription factors and is essential for stimulated gene transcription. YAP/TAZ-TEAD facilitates the upregulation of multiple genes involved in evolutionary cell proliferation and survival. TEAD1-4 overexpression has been observed in different cancers in various tissues, making TEAD an attractive target for drug development. The central drug-accessible pocket of TEAD is crucial because it undergoes a post-translational modification called auto-palmitoylation. Crystal structures of the C-terminal TEAD complex with small molecules are available in the Protein Data Bank, aiding structure-based drug design. In this study, we utilized the fragment molecular orbital (FMO) method, molecular dynamics (MD) simulations, shape-based screening, and molecular mechanics-generalized Born surface area (MM-GBSA) calculations for virtual screening, and we identified a novel non-covalent inhibitor-BC-001-with IC50 = 3.7 µM in a reporter assay. Subsequently, we optimized several analogs of BC-001 and found that the optimized compound BC-011 exhibited an IC50 of 72.43 nM. These findings can be used to design effective TEAD modulators with anticancer therapeutic implications.

Fragment molecular orbital-based variational quantum eigensolver for quantum chemistry in the age of quantum computing.

Quantum computers offer significant potential for complex system analysis, yet their application in large systems is hindered by limitations such as qubit availability and quantum hardware noise. While the variational quantum eigensolver (VQE) was proposed to address these issues, its scalability remains limited. Many efforts, including new ansätze and Hamiltonian modifications, have been made to overcome these challenges. In this work, we introduced the novel Fragment Molecular Orbital/Variational Quantum Eigensolver (FMO/VQE) algorithm. This method combines the fragment molecular orbital (FMO) approach with VQE and efficiently utilizes qubits for quantum chemistry simulations. Employing the UCCSD ansatz, the FMO/VQE achieved an absolute error of just 0.053 mHa with 8 qubits in a [Formula: see text] system using the STO-3G basis set, and an error of 1.376 mHa with 16 qubits in a [Formula: see text] system with the 6-31G basis set. These results indicated a significant advancement in scalability over conventional VQE, maintaining accuracy with fewer qubits. Therefore, our FMO/VQE method exemplifies how integrating fragment-based quantum chemistry with quantum algorithms can enhance scalability, facilitating more complex molecular simulations and aligning with quantum computing advancements.

Transformer-Based Molecular Generative Model for Antiviral Drug Design.

Since the Simplified Molecular Input Line Entry System (SMILES) is oriented to the atomic-level representation of molecules and is not friendly in terms of human readability and editable, however, IUPAC is the closest to natural language and is very friendly in terms of human-oriented readability and performing molecular editing, we can manipulate IUPAC to generate corresponding new molecules and produce programming-friendly molecular forms of SMILES. In addition, antiviral drug design, especially analogue-based drug design, is also more appropriate to edit and design directly from the functional group level of IUPAC than from the atomic level of SMILES, since designing analogues involves altering the R group only, which is closer to the knowledge-based molecular design of a chemist. Herein, we present a novel data-driven self-supervised pretraining generative model called "TransAntivirus" to make select-and-replace edits and convert organic molecules into the desired properties for design of antiviral candidate analogues. The results indicated that TransAntivirus is significantly superior to the control models in terms of novelty, validity, uniqueness, and diversity. TransAntivirus showed excellent performance in the design and optimization of nucleoside and non-nucleoside analogues by chemical space analysis and property prediction analysis. Furthermore, to validate the applicability of TransAntivirus in the design of antiviral drugs, we conducted two case studies on the design of nucleoside analogues and non-nucleoside analogues and screened four candidate lead compounds against anticoronavirus disease (COVID-19). Finally, we recommend this framework for accelerating antiviral drug discovery.

PhyloSophos: a high-throughput scientific name mapping algorithm augmented with explicit consideration of taxonomic science, and its application on natural product (NP) occurrence database processing.

BACKGROUND: The standardization of biological data using unique identifiers is vital for seamless data integration, comprehensive interpretation, and reproducibility of research findings, contributing to advancements in bioinformatics and systems biology. Despite being widely accepted as a universal identifier, scientific names for biological species have inherent limitations, including lack of stability, uniqueness, and convertibility, hindering their effective use as identifiers in databases, particularly in natural product (NP) occurrence databases, posing a substantial obstacle to utilizing this valuable data for large-scale research applications. RESULT: To address these challenges and facilitate high-throughput analysis of biological data involving scientific names, we developed PhyloSophos, a Python package that considers the properties of scientific names and taxonomic systems to accurately map name inputs to entries within a chosen reference database. We illustrate the importance of assessing multiple taxonomic databases and considering taxonomic syntax-based pre-processing using NP occurrence databases as an example, with the ultimate goal of integrating heterogeneous information into a single, unified dataset. CONCLUSIONS: We anticipate PhyloSophos to significantly aid in the systematic processing of poorly digitized and curated biological data, such as biodiversity information and ethnopharmacological resources, enabling full-scale bioinformatics analysis using these valuable data resources.

Leveraging the Fragment Molecular Orbital Method to Explore the PLK1 Kinase Binding Site and Polo-Box Domain for Potent Small-Molecule Drug Design.

Polo-like kinase 1 (PLK1) plays a pivotal role in cell division regulation and emerges as a promising therapeutic target for cancer treatment. Consequently, the development of small-molecule inhibitors targeting PLK1 has become a focal point in contemporary research. The adenosine triphosphate (ATP)-binding site and the polo-box domain in PLK1 present crucial interaction sites for these inhibitors, aiming to disrupt the protein's function. However, designing potent and selective small-molecule inhibitors can be challenging, requiring a deep understanding of protein-ligand interaction mechanisms at these binding sites. In this context, our study leverages the fragment molecular orbital (FMO) method to explore these site-specific interactions in depth. Using the FMO approach, we used the FMO method to elucidate the molecular mechanisms of small-molecule drugs binding to these sites to design PLK1 inhibitors that are both potent and selective. Our investigation further entailed a comparative analysis of various PLK1 inhibitors, each characterized by distinct structural attributes, helping us gain a better understanding of the relationship between molecular structure and biological activity. The FMO method was particularly effective in identifying key binding features and predicting binding modes for small-molecule ligands. Our research also highlighted specific "hot spot" residues that played a critical role in the selective and robust binding of PLK1. These findings provide valuable insights that can be used to design new and effective PLK1 inhibitors, which can have significant implications for developing anticancer therapeutics.

Targeting FLT3-TAZ signaling to suppress drug resistance in blast phase chronic myeloid leukemia.

BACKGROUND: Although the development of BCR::ABL1 tyrosine kinase inhibitors (TKIs) rendered chronic myeloid leukemia (CML) a manageable condition, acquisition of drug resistance during blast phase (BP) progression remains a critical challenge. Here, we reposition FLT3, one of the most frequently mutated drivers of acute myeloid leukemia (AML), as a prognostic marker and therapeutic target of BP-CML. METHODS: We generated FLT3 expressing BCR::ABL1 TKI-resistant CML cells and enrolled phase-specific CML patient cohort to obtain unpaired and paired serial specimens and verify the role of FLT3 signaling in BP-CML patients. We performed multi-omics approaches in animal and patient studies to demonstrate the clinical feasibility of FLT3 as a viable target of BP-CML by establishing the (1) molecular mechanisms of FLT3-driven drug resistance, (2) diagnostic methods of FLT3 protein expression and localization, (3) association between FLT3 signaling and CML prognosis, and (4) therapeutic strategies to tackle FLT3+ CML patients. RESULTS: We reposition the significance of FLT3 in the acquisition of drug resistance in BP-CML, thereby, newly classify a FLT3+ BP-CML subgroup. Mechanistically, FLT3 expression in CML cells activated the FLT3-JAK-STAT3-TAZ-TEAD-CD36 signaling pathway, which conferred resistance to a wide range of BCR::ABL1 TKIs that was independent of recurrent BCR::ABL1 mutations. Notably, FLT3+ BP-CML patients had significantly less favorable prognosis than FLT3- patients. Remarkably, we demonstrate that repurposing FLT3 inhibitors combined with BCR::ABL1 targeted therapies or the single treatment with ponatinib alone can overcome drug resistance and promote BP-CML cell death in patient-derived FLT3+ BCR::ABL1 cells and mouse xenograft models. CONCLUSION: Here, we reposition FLT3 as a critical determinant of CML progression via FLT3-JAK-STAT3-TAZ-TEAD-CD36 signaling pathway that promotes TKI resistance and predicts worse prognosis in BP-CML patients. Our findings open novel therapeutic opportunities that exploit the undescribed link between distinct types of malignancies.

Intracellular Loop in the Brain Isoforms of Anoctamin 2 Channels Regulates Calcium-dependent Activation.

Anoctamin 2 (ANO2 or TMEM16B), a calcium-activated chloride channel (CaCC), performs diverse roles in neurons throughout the central nervous system. In hippocampal neurons, ANO2 narrows action potential width and reduces postsynaptic depolarization with high sensitivity to Ca2+ at relatively fast kinetics. In other brain regions, including the thalamus, ANO2 mediates activity-dependent spike frequency adaptations with low sensitivity to Ca2+ at relatively slow kinetics. How this same channel can respond to a wide range of Ca2+ levels remains unclear. We hypothesized that splice variants of ANO2 may contribute to its distinct Ca2+ sensitivity, and thus its diverse neuronal functions. We identified two ANO2 isoforms expressed in mouse brains and examined their electrophysiological properties: isoform 1 (encoded by splice variants with exons 1a, 2, 4, and 14) was expressed in the hippocampus, while isoform 2 (encoded by splice variants with exons 1a, 2, and 4) was broadly expressed throughout the brain, including in the cortex and thalamus, and had a slower calcium-dependent activation current than isoform 1. Computational modeling revealed that the secondary structure of the first intracellular loop of isoform 1 forms an entrance cavity to the calcium-binding site from the cytosol that is relatively larger than that in isoform 2. This difference provides structural evidence that isoform 2 is involved in accommodating spike frequency, while isoform 1 is involved in shaping the duration of an action potential and decreasing postsynaptic depolarization. Our study highlights the roles and molecular mechanisms of specific ANO2 splice variants in modulating neuronal functions.

Prediction of polyreactive and nonspecific single-chain fragment variables through structural biochemical features and protein language-based descriptors.

BACKGROUND: Monoclonal antibodies (mAbs) have been used as therapeutic agents, which must overcome many developability issues after the discovery from in vitro display libraries. Especially, polyreactive mAbs can strongly bind to a specific target and weakly bind to off-target proteins, which leads to poor antibody pharmacokinetics in clinical development. Although early assessment of polyreactive mAbs is important in the early discovery stage, experimental assessments are usually time-consuming and expensive. Therefore, computational approaches for predicting the polyreactivity of single-chain fragment variables (scFvs) in the early discovery stage would be promising for reducing experimental efforts. RESULTS: Here, we made prediction models for the polyreactivity of scFvs with the known polyreactive antibody features and natural language model descriptors. We predicted 19,426 protein structures of scFvs with trRosetta to calculate the polyreactive antibody features and investigated the classifying performance of each factor for polyreactivity. In the known polyreactive features, the net charge of the CDR2 loop, the tryptophan and glycine residues in CDR-H3, and the lengths of the CDR1 and CDR2 loops, importantly contributed to the performance of the models. Additionally, the hydrodynamic features, such as partial specific volume, gyration radius, and isoelectric points of CDR loops and scFvs, were newly added to improve model performance. Finally, we made the prediction model with a robust performance ([Formula: see text]) with an ensemble learning of the top 3 best models. CONCLUSION: The prediction models for polyreactivity would help assess polyreactive scFvs in the early discovery stage and our approaches would be promising to develop machine learning models with quantitative data from high throughput assays for antibody screening.

Development of an innovative data-driven system to generate descriptive prediction equation of dielectric constant on small sample sets.

Dielectric constant (DC, ε) is a fundamental parameter in material sciences to measure polarizability of the system. In industrial processes, its value is an imperative indicator, which demonstrates the dielectric property of material and compiles information including separation information, chemical equilibrium, chemical reactivity analysis, and solubility modeling. Since, the available ε-prediction models are fairly primitive and frequently suffer from serious failures especially when deals with strong polar compounds. Therefore, we have developed a novel data-driven system to improve the efficiency and wide-range applicability of ε using in material sciences. This innovative scheme adopts the correlation distance and genetic algorithm to discriminate features' combination and avoid overfitting. Herein, the prediction output of the single ML model as a coding to estimate the target value by simulating the layer-by-layer extraction in deep learning, and enabling instant search for the optimal combination of features is recruited. Our model established an improved correlation value of 0.956 with target as compared to the previously available best traditional ML result of 0.877. Our framework established a profound improvement, especially for material systems possessing ε value >50. In terms of interpretability, we have derived a conceptual computational equation from a minimum generating tree. Our innovative data-driven system is preferentially superior over other methods due to its application for the prediction of dielectric constants as well as for the prediction of overall micro and macro-properties of any multi-components complex.

De novo molecular design with deep molecular generative models for PPI inhibitors.

We construct a protein-protein interaction (PPI) targeted drug-likeness dataset and propose a deep molecular generative framework to generate novel drug-likeness molecules from the features of the seed compounds. This framework gains inspiration from published molecular generative models, uses the key features associated with PPI inhibitors as input and develops deep molecular generative models for de novo molecular design of PPI inhibitors. For the first time, quantitative estimation index for compounds targeting PPI was applied to the evaluation of the molecular generation model for de novo design of PPI-targeted compounds. Our results estimated that the generated molecules had better PPI-targeted drug-likeness and drug-likeness. Additionally, our model also exhibits comparable performance to other several state-of-the-art molecule generation models. The generated molecules share chemical space with iPPI-DB inhibitors as demonstrated by chemical space analysis. The peptide characterization-oriented design of PPI inhibitors and the ligand-based design of PPI inhibitors are explored. Finally, we recommend that this framework will be an important step forward for the de novo design of PPI-targeted therapeutics.

Quantum computational study of chloride attack on chloromethane for chemical accuracy and quantum noise effects with UCCSD and k-UpCCGSD ansatzes.

Quantum computing is expected to play an important role in solving the problem of huge computational costs in various applications by utilizing the collective properties of quantum states, including superposition, interference, and entanglement, to perform computations. Quantum mechanical (QM) methods are candidates for various applications and can provide accurate absolute energy calculations in structure-based methods. QM methods are powerful tools for describing reaction pathways and their potential energy surfaces (PES). In this study, we applied quantum computing to describe the PES of the bimolecular nucleophilic substitution (SN2) reaction between chloromethane and chloride ions. We performed noiseless and noise simulations using quantum algorithms and compared the accuracy and noise effects of the ansatzes. In noiseless simulations, the results from UCCSD and k-UpCCGSD are similar to those of full configurational interaction (FCI) with the same active space, which indicates that quantum algorithms can describe the PES of the SN2 reaction. In noise simulations, UCCSD is more susceptible to quantum noise than k-UpCCGSD. Therefore, k-UpCCGSD can serve as an alternative to UCCSD to reduce quantum noisy effects in the noisy intermediate-scale quantum era, and k-UpCCGSD is sufficient to describe the PES of the SN2 reaction in this work. The results showed the applicability of quantum computing to the SN2 reaction pathway and provided valuable information for structure-based molecular simulations with quantum computing.

Identification of Novel Natural Product Inhibitors against Matrix Metalloproteinase 9 Using Quantum Mechanical Fragment Molecular Orbital-Based Virtual Screening Methods.

Matrix metalloproteinases (MMPs) are calcium-dependent zinc-containing endopeptidases involved in multiple cellular processes. Among the MMP isoforms, MMP-9 regulates cancer invasion, rheumatoid arthritis, and osteoarthritis by degrading extracellular matrix proteins present in the tumor microenvironment and cartilage and promoting angiogenesis. Here, we identified two potent natural product inhibitors of the non-catalytic hemopexin domain of MMP-9 using a novel quantum mechanical fragment molecular orbital (FMO)-based virtual screening workflow. The workflow integrates qualitative pharmacophore modeling, quantitative binding affinity prediction, and a raw material search of natural product inhibitors with the BMDMS-NP library. In binding affinity prediction, we made a scoring function with the FMO method and applied the function to two protein targets (acetylcholinesterase and fibroblast growth factor 1 receptor) from DUD-E benchmark sets. In the two targets, the FMO method outperformed the Glide docking score and MM/PBSA methods. By applying this workflow to MMP-9, we proposed two potent natural product inhibitors (laetanine 9 and genkwanin 10) that interact with hotspot residues of the hemopexin domain of MMP-9. Laetanine 9 and genkwanin 10 bind to MMP-9 with a dissociation constant (KD) of 21.6 and 0.614 µM, respectively. Overall, we present laetanine 9 and genkwanin 10 for MMP-9 and demonstrate that the novel FMO-based workflow with a quantum mechanical approach is promising to discover potent natural product inhibitors of MMP-9, satisfying the pharmacophore model and good binding affinity.

Evaluation of protein descriptors in computer-aided rational protein engineering tasks and its application in property prediction in SARS-CoV-2 spike glycoprotein.

The importance of protein engineering in the research and development of biopharmaceuticals and biomaterials has increased. Machine learning in computer-aided protein engineering can markedly reduce the experimental effort in identifying optimal sequences that satisfy the desired properties from a large number of possible protein sequences. To develop general protein descriptors for computer-aided protein engineering tasks, we devised new protein descriptors, one sequence-based descriptor (PCgrades), and three structure-based descriptors (PCspairs, 3D-SPIEs_5.4 Å, and 3D-SPIEs_8Å). While the PCgrades and PCspairs include general and statistical information in physicochemical properties in single and pairwise amino acids respectively, the 3D-SPIEs include specific and quantum-mechanical information with parameterized quantum mechanical calculations (FMO2-DFTB3/D/PCM). To evaluate the protein descriptors, we made prediction models with the new descriptors and previously developed descriptors for diverse protein datasets including protein expression and binding affinity change in SARS-CoV-2 spike glycoprotein. As a result, the newly devised descriptors showed a good performance in diverse datasets, in which the PCspairs showed the best performance ( R 2 = 0.783 for protein expression and R 2 = 0.711 for binding affinity). As a result, the newly devised descriptors showed a good performance in diverse datasets, in which the PCspairs showed the best performance. Similar approaches with those descriptors would be promising and useful if the prediction models are trained with sufficient quantitative experimental data from high-throughput assays for industrial enzymes or protein drugs.

Comprehensive strategies of machine-learning-based quantitative structure-activity relationship models.

Early quantitative structure-activity relationship (QSAR) technologies have unsatisfactory versatility and accuracy in fields such as drug discovery because they are based on traditional machine learning and interpretive expert features. The development of Big Data and deep learning technologies significantly improve the processing of unstructured data and unleash the great potential of QSAR. Here we discuss the integration of wet experiments (which provide experimental data and reliable verification), molecular dynamics simulation (which provides mechanistic interpretation at the atomic/molecular levels), and machine learning (including deep learning) techniques to improve QSAR models. We first review the history of traditional QSAR and point out its problems. We then propose a better QSAR model characterized by a new iterative framework to integrate machine learning with disparate data input. Finally, we discuss the application of QSAR and machine learning to many practical research fields, including drug development and clinical trials.

Hot Spot Analysis of YAP-TEAD Protein-Protein Interaction Using the Fragment Molecular Orbital Method and Its Application for Inhibitor Discovery.

The Hippo pathway is an important signaling pathway modulating growth control and cancer cell proliferation. Dysregulation of the Hippo pathway is a common feature of several types of cancer cells. The modulation of the interaction between yes-associated protein (YAP) and transcriptional enhancer associated domain (TEAD) in the Hippo pathway is considered an attractive target for cancer therapeutic development, although the inhibition of PPI is a challenging task. In order to investigate the hot spots of the YAP and TEAD1 interacting complex, an ab initio Fragment Molecular Orbital (FMO) method was introduced. With the hot spots, pharmacophores for the inhibitor design were constructed, then virtual screening was performed to an in-house library. Next, we performed molecular docking simulations and FMO calculations for screening results to study the binding modes and affinities between PPI inhibitors and TEAD1. As a result of the virtual screening, three compounds were selected as virtual hit compounds. In order to confirm their biological activities, cellular (luciferase activity, proximity ligation assay and wound healing assay in A375 cells, qRT-PCR in HEK 293T cells) and biophysical assays (surface plasmon resonance assays) were performed. Based on the findings of the study, we propose a novel PPI inhibitor BY03 and demonstrate a profitable strategy to analyze YAP-TEAD PPI and discover novel PPI inhibitors.

Development of an Infinite Dilution Activity Coefficient Prediction Model for Organic Solutes in Ionic Liquids with Modified Partial Equalization Orbital Electronegativity Method Derived Descriptors.

The objective of this study was to develop a robust prediction model for the infinite dilution activity coefficients (γ ∞) of organic molecules in diverse ionic liquid (IL) solvents. Electrostatic, hydrogen bond, polarizability, molecular structure, and temperature terms were used in model development. A feed-forward model based on artificial neural networks was developed with 34,754 experimental activity coefficients, a combination of 195 IL solvents (88 cations and 38 anions), and 147 organic solutes at a temperature range of 298 to 408 K. The root mean squared error (RMSE) of the training set and test set was 0.219 and 0.235, respectively. The R 2 of the training set and the test set was 0.984 and 0.981, respectively. The applicability domain was determined through a Williams plot, which implied that water and halogenated compounds were outside of the applicability domain. The robustness test shows that the developed model is robust. The web server supports using the developed prediction model and is freely available at https://preadmet.bmdrc.kr/activitycoefficient_mainpage/prediction/.

Modeling lipid layers of atopic skin and observation of changes in lipid layer properties with changes in ceramide content.

BACKGROUND: Studies have shown that there is a high correlation between atopic dermatitis and decrease in ceramide content in the lipid bilayer of skin. Moreover, it has been shown that the reduction in ceramide content in the stratum corneum is unique to atopic dermatitis, indicating that there are particular structural differences between the lipid bilayers of normal and atopic skin. AIM: This study aimed to compare the lipid bilayer of the atopic skin with that of the healthy skin and to establish a structural model of the lipid bilayer for atopy. METHODS: Molecular dynamics simulations were performed using NAMD 2.8. Models of lipid bilayers of normal skin and atopic skin, and a model of lipid bilayer containing only ceramide were built with CHARMM-GUI. The thickness, area occupied per lipid, and alignment of lipids were compared among the three models. Potential mean force (PMF) of the sodium laureth sulfate (SLES) on lipid bilayers was calculated to predict the affinity between SLES and lipid bilayers. RESULTS: Potential mean force calculations showed that the lipid bilayer of atopic skin was able to absorb the surfactant more easily than that of normal skin. CONCLUSIONS: When the ceramide ratio is low, the thickness of lipid bilayer is reduced and its structure is weakened. Other structural differences between the lipid layers of normal and atopic skin included increased area per lipid and poor alignment of lipids. Further, the atopy lipid bilayer model was found to absorb more SLES than the normal skin lipid bilayer model.

Isolindleyin exerts anti-melanogenic effects in human epidermal melanocytes via direct binding to tyrosinase.

To overcome dermatological concerns causing abnormally excessive melanin synthesis, highly effective and safe skin depigmentation compounds have been identified in the cosmetic and pharmaceutical industries. Among several methods used to achieve skin depigmentation, inhibition of tyrosinase is one of the most effective, since tyrosinase is a crucial enzyme in melanogenesis. Herein, isolindleyin, a novel inhibitor of human tyrosinase, was introduced and evaluated for its anti-melanogenic effects in human epidermal melanocytes. The results revealed that isolindleyin was directly bound to tyrosinase and it suppressed melanin synthesis. The binding mode between isolindleyin and the active sites of human tyrosinase was investigated using computational molecular docking at the atomic level. Isolindleyin binding was found to be stabilized by hydrophobic interactions between His 367 and Val 377 and by hydrogen bonds between Ser 380 and Asn 364. The results of this study revealed the anti-melanogenic effects of isolindleyin that could contribute toward overcoming dermatological concerns that cause abnormally excessive melanin synthesis.

Hot spot profiles of SARS-CoV-2 and human ACE2 receptor protein protein interaction obtained by density functional tight binding fragment molecular orbital method.

The prevalence of a novel ß-coronavirus (SARS-CoV-2) was declared as a public health emergency of international concern on 30 January 2020 and a global pandemic on 11 March 2020 by WHO. The spike glycoprotein of SARS-CoV-2 is regarded as a key target for the development of vaccines and therapeutic antibodies. In order to develop anti-viral therapeutics for SARS-CoV-2, it is crucial to find amino acid pairs that strongly attract each other at the interface of the spike glycoprotein and the human angiotensin-converting enzyme 2 (hACE2) complex. In order to find hot spot residues, the strongly attracting amino acid pairs at the protein-protein interaction (PPI) interface, we introduce a reliable inter-residue interaction energy calculation method, FMO-DFTB3/D/PCM/3D-SPIEs. In addition to the SARS-CoV-2 spike glycoprotein/hACE2 complex, the hot spot residues of SARS-CoV-1 spike glycoprotein/hACE2 complex, SARS-CoV-1 spike glycoprotein/antibody complex, and HCoV-NL63 spike glycoprotein/hACE2 complex were obtained using the same FMO method. Following this, a 3D-SPIEs-based interaction map was constructed with hot spot residues for the hACE2/SARS-CoV-1 spike glycoprotein, hACE2/HCoV-NL63 spike glycoprotein, and hACE2/SARS-CoV-2 spike glycoprotein complexes. Finally, the three 3D-SPIEs-based interaction maps were combined and analyzed to find the consensus hot spots among the three complexes. As a result of the analysis, two hot spots were identified between hACE2 and the three spike proteins. In particular, E37, K353, G354, and D355 of the hACE2 receptor strongly interact with the spike proteins of coronaviruses. The 3D-SPIEs-based map would provide valuable information to develop anti-viral therapeutics that inhibit PPIs between the spike protein of SARS-CoV-2 and hACE2.