AI-Powered Mining of Highly Customized and Superior ESIPT-Based Fluorescent Probes.

Excited-state intramolecular proton transfer (ESIPT) has attracted great attention in fluorescent sensors and luminescent materials due to its unique photobiological and photochemical features. However, the current structures are far from meeting the specific demands for ESIPT molecules in different scenarios; the try-and-error development method is labor-intensive and costly. Therefore, it is imperative to devise novel approaches for the exploration of promising ESIPT fluorophores. This research proposes an artificial intelligence approach aiming at exploring ESIPT molecules efficiently. The first high-quality ESIPT dataset and a multi-level prediction system are constructed that realized accurate identification of ESIPT molecules from a large number of compounds under a stepwise distinguishing from conventional molecules to fluorescent molecules and then to ESIPT molecules. Furthermore, key structural features that contributed to ESIPT are revealed by using the SHapley Additive exPlanations (SHAP) method. Then three strategies are proposed to ensure the ESIPT process while keeping good safety, pharmacokinetic properties, and novel structures. With these strategies, >700 previously unreported ESIPT molecules are screened from a large pool of 570 000 compounds. The ESIPT process and biosafety of optimal molecules are successfully validated by quantitative calculation and experiment. This novel approach is expected to bring a new paradigm for exploring ideal ESIPT molecules.

Rational Design of Organelle-Targeted Fluorescent Probes: Insights from Artificial Intelligence.

Monitoring the physiological changes of organelles is essential for understanding the local biological information of cells and for improving the diagnosis and therapy of diseases. Currently, fluorescent probes are considered as the most powerful tools for imaging and have been widely applied in biomedical fields. However, the expected targeting effects of these probes are often inconsistent with the real experiments. The design of fluorescent probes mainly depends on the empirical knowledge of researchers, which was inhibited by limited chemical space and low efficiency. Herein, we proposed a novel multilevel framework for the prediction of organelle-targeted fluorescent probes by employing advanced artificial intelligence algorithms. In this way, not only the targeting mechanism could be interpreted beyond intuitions but also a quick evaluation method could be established for the rational design. Furthermore, the targeting and imaging powers of the optimized and synthesized probes based on this methodology were verified by quantitative calculation and experiments.

Dynamic Interactions Between Brilliant Green and MscL Investigated by Solid-State NMR Spectroscopy and Molecular Dynamics Simulations.

The mechanosensitive ion channel of large conductance (MscL) is a promising template for the development of new antibiotics due to its high conservation and uniqueness to microbes. Brilliant green (BG), a triarylmethane dye, has been identified as a new antibiotic targeted MscL. However, the detailed binding sites to MscL and the dynamic pathway of BG through the MscL channel remain unknown. Here, the dynamic interactions between BG and MscL were investigated using solid-state NMR spectroscopy and molecule dynamics (MD) simulations. Residue site-specific binding sites of BG to the MscL channel were identified by solid-state NMR. In addition, MD simulations revealed that BG conducts through the MscL channel via residues along the inner surface of the pore sequentially, in which the strong hydrophobic interactions between BG and hydrophobic residues F23 and I27 in the hydrophobic gate region of the MscL channel are major restrictions. Particularly, it was demonstrated that BG activates the MscL channel by reducing the hydrophobicity of the F23 in the gate region by water molecules that are bound to BG. Taken together, these simulations and experimental data provide novel insights into the dynamic interactions between BG and MscL, based on which new hydrophobic antibiotics and adjuvants targeting MscL can be developed.

Molecular Mechanisms of Mercury-Sensitive Aquaporins.

Aquaporins are transmembrane channels that allow for the passive permeation of water and other small molecules across biological membranes. Their channel activities are sensitive to mercury ions. Intriguingly, while most aquaporins are inhibited by mercury ions, several aquaporins are activated by mercury ions. The molecular basis of the opposing aquaporin regulation by mercury remains poorly understood. Herein, we investigated AqpZ inhibition and AQP6 activation upon binding of mercury ions using solid-state NMR (ssNMR) and molecular dynamics (MD) simulations. Based on the structure of the Hg-AqpZ complex constructed by MD simulations and ssNMR, we identified that the pore closure was caused by mercury-induced conformational changes of the key residue R189 in the selectivity filter region, while pore opening was caused by conformational changes of residues H181 and R196 in the selectivity filter region in AQP6. Both conformational changes were caused by the disruption of the H-bond network of R189/R196 by mercury. The molecular details provided a structural basis for mercury-mediated functional changes in aquaporins.

Mechanism of unusual AQP6 activation by mercury binding to a pore-external residue C155.

Aquaporins (AQPs) transport water molecules across cell membranes. Although most aquaporins are inhibited by mercury ions, AQP6 was reported to be activated by binding mercury ions to residues C155 and C190. Different from C190 and the other pore-line cysteine residues, C155 is located outside the pore, thus not directly affecting the internal pathway by mercury binding to it. The molecular mechanism of unusual water channel activation by mercury ion binding to the C155 site remains unknown. Here, we investigate the activation of AQP6 by mercury ions binding to C155 by molecular dynamics (MD) simulations. The MD simulation results show that the mercury-induced water permeation activation is derived from the conformational change of a pore-line residue M160, from a point-to-pore conformation before mercury binding to an away-pore conformation after mercury binding. The conformation change of M160 is derived from the reduction of the hydrogen bonding between C155 and S159 in the α-helix with the coordination of C155 to mercury ion altering their conformation significantly. This study reveals the complex mechanism of water channel activation by mercury ion binding to pore-external residues in water channels.

Versatile ginsenoside Rg3 liposomes inhibit tumor metastasis by capturing circulating tumor cells and destroying metastatic niches.

Limited circulating tumor cells (CTCs) capturing efficiency and lack of regulation capability on CTC-supportive metastatic niches (MNs) are two main obstacles hampering the clinical translation of conventional liposomes for the treatment of metastatic breast cancers. Traditional delivery strategies, such as ligand modification and immune modulator co-encapsulation for nanocarriers, are inefficient and laborious. Here, a multifunctional Rg3 liposome loading with docetaxel (Rg3-Lp/DTX) was developed, in which Rg3 was proved to intersperse in the phospholipid bilayer and exposed its glycosyl on the liposome surface. Therefore, it exhibited much higher CTC-capturing efficiency via interaction with glucose transporter 1 (Glut1) overexpressed on CTCs. After reaching the lungs with CTCs, Rg3 inhibited the formation of MNs by reversing the immunosuppressive microenvironment. Together, Rg3-Lp/DTX exhibited excellent metastasis inhibition capacity by CTC ("seeds") neutralization and MN ("soil") inhibition. The strategy has great clinical translation prospects for antimetastasis treatment with enhanced therapeutic efficacy and simple preparation process.

Molecular docking-based computational platform for high-throughput virtual screening.

Structure-based virtual screening is a key, routine computational method in computer-aided drug design. Such screening can be used to identify potentially highly active compounds, to speed up the progress of novel drug design. Molecular docking-based virtual screening can help find active compounds from large ligand databases by identifying the binding affinities between receptors and ligands. In this study, we analyzed the challenges of virtual screening, with the aim of identifying highly active compounds faster and more easily than is generally possible. We discuss the accuracy and speed of molecular docking software and the strategy of high-throughput molecular docking calculation, and we focus on current challenges and our solutions to these challenges of ultra-large-scale virtual screening. The development of Web services helps lower the barrier to drug virtual screening. We introduced some related web sites for docking and virtual screening, focusing on the development of pre- and post-processing interactive visualization and large-scale computing.

Glycan Epitopes and Potential Glycoside Antagonists of DC-SIGN Involved in COVID-19: In Silico Study.

Glycosylation is an important post-translational modification that affects a wide variety of physiological functions. DC-SIGN (Dendritic Cell-Specific Intercellular adhesion molecule-3-Grabbing Non-integrin) is a protein expressed in antigen-presenting cells that recognizes a variety of glycan epitopes. Until now, the binding of DC-SIGN to SARS-CoV-2 Spike glycoprotein has been reported in various articles and is regarded to be a factor in systemic infection and cytokine storm. The mechanism of DC-SIGN recognition offers an alternative method for discovering new medication for COVID-19 treatment. Here, we discovered three potential pockets that hold different glycan epitopes by performing molecular dynamics simulations of previously reported oligosaccharides. The "EPN" motif, "NDD" motif, and Glu354 form the most critical pocket, which is known as the Core site. We proposed that the type of glycan epitopes, rather than the precise amino acid sequence, determines the recognition. Furthermore, we deduced that oligosaccharides could occupy an additional site, which adds to their higher affinity than monosaccharides. Based on our findings and previously described glycoforms on the SARS-CoV-2 Spike, we predicted the potential glycan epitopes for DC-SIGN. It suggested that glycan epitopes could be recognized at multiple sites, not just Asn234, Asn149 and Asn343. Subsequently, we found that Saikosaponin A and Liquiritin, two plant glycosides, were promising DC-SIGN antagonists in silico.

Molecular simulation studies of the interactions between the human/pangolin/cat/bat ACE2 and the receptor binding domain of the SARS-CoV-2 spike protein.

The recent outbreak of SARS-CoV-2 has had a profound effect on the world. Similar to that in SARS-CoV, the entry receptor of SARS-CoV-2 is ACE2. The binding of SARS-CoV-2 spike protein to ACE2 is the critical to the virus infection. Recently multiple species (human, Chinese chrysanthemum, Malay pangolin and cat) have been reported to be susceptible to the virus infection. However, the binding capacity and the detailed binding mechanism of SARS-CoV-2 spike protein to ACE2 of these species remains unexplored. Herein free energy calculations with MM-GBSA and Potential of Mean Forces together reveal that the Human-SARS-CoV-2 has a higher stability tendency than Human-SARS-CoV. Meanwhile, we uncover that SARS-CoV-2 has an enhanced ability to bind with the ACE2 in humans, pangolins and cats compared to that in bats. Analysis of key residues with energy decomposition and residue contact maps reveal several important consensus sites in ACE2s among the studied species, and determined the more favorable specified residues among the different types of amino acids. These results provide important implications for understanding SARS-CoV-2 host range which will make it possible to control the spread of the virus and use of animal models, targeted drug screening and vaccine candidates against SARS-CoV-2.

Perspective of drug design with high-performance computing.

The representative applications, recent advances and possible future directions of computational drug design were summarized, aiming to accelerate the drug discovery with the assistance of the fast-developing high-performance computing.

Identify potent SARS-CoV-2 main protease inhibitors via accelerated free energy perturbation-based virtual screening of existing drugs.

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has become a global crisis. There is no therapeutic treatment specific for COVID-19. It is highly desirable to identify potential antiviral agents against SARS-CoV-2 from existing drugs available for other diseases and thus repurpose them for treatment of COVID-19. In general, a drug repurposing effort for treatment of a new disease, such as COVID-19, usually starts from a virtual screening of existing drugs, followed by experimental validation, but the actual hit rate is generally rather low with traditional computational methods. Here we report a virtual screening approach with accelerated free energy perturbation-based absolute binding free energy (FEP-ABFE) predictions and its use in identifying drugs targeting SARS-CoV-2 main protease (Mpro). The accurate FEP-ABFE predictions were based on the use of a restraint energy distribution (RED) function, making the practical FEP-ABFE-based virtual screening of the existing drug library possible. As a result, out of 25 drugs predicted, 15 were confirmed as potent inhibitors of SARS-CoV-2 Mpro The most potent one is dipyridamole (inhibitory constant Ki = 0.04 µM) which has shown promising therapeutic effects in subsequently conducted clinical studies for treatment of patients with COVID-19. Additionally, hydroxychloroquine (Ki = 0.36 µM) and chloroquine (Ki = 0.56 µM) were also found to potently inhibit SARS-CoV-2 Mpro We anticipate that the FEP-ABFE prediction-based virtual screening approach will be useful in many other drug repurposing or discovery efforts.

Anti-SARS-CoV-2 activities in vitro of Shuanghuanglian preparations and bioactive ingredients.

Human infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19) and there is no cure currently. The 3CL protease (3CLpro) is a highly conserved protease which is indispensable for CoVs replication, and is a promising target for development of broad-spectrum antiviral drugs. In this study we investigated the anti-SARS-CoV-2 potential of Shuanghuanglian preparation, a Chinese traditional patent medicine with a long history for treating respiratory tract infection in China. We showed that either the oral liquid of Shuanghuanglian, the lyophilized powder of Shuanghuanglian for injection or their bioactive components dose-dependently inhibited SARS-CoV-2 3CLpro as well as the replication of SARS-CoV-2 in Vero E6 cells. Baicalin and baicalein, two ingredients of Shuanghuanglian, were characterized as the first noncovalent, nonpeptidomimetic inhibitors of SARS-CoV-2 3CLpro and exhibited potent antiviral activities in a cell-based system. Remarkably, the binding mode of baicalein with SARS-CoV-2 3CLpro determined by X-ray protein crystallography was distinctly different from those of known 3CLpro inhibitors. Baicalein was productively ensconced in the core of the substrate-binding pocket by interacting with two catalytic residues, the crucial S1/S2 subsites and the oxyanion loop, acting as a "shield" in front of the catalytic dyad to effectively prevent substrate access to the catalytic dyad within the active site. Overall, this study provides an example for exploring the in vitro potency of Chinese traditional patent medicines and effectively identifying bioactive ingredients toward a specific target, and gains evidence supporting the in vivo studies of Shuanghuanglian oral liquid as well as two natural products for COVID-19 treatment.

Structural Mechanism of Barriers to Interspecies Seeding Transmissibility of Full-Length Prion Protein Amyloid.

A puzzling feature of prion diseases is the cross-species barriers. The detailed molecular mechanisms underlying these interspecies barriers remain poorly understood because of a lack of high-resolution structural information on the scrapie isoform of the prion protein (PrPSc ). In this study we identified the critical role of the residues 165/167 in the barrier to seeding mouse PrP (mPrP) fibril seeds to human cellular prion protein (PrPC ). Solid-state NMR revealed a C-terminal ß-sheet core spanning residues 165-230 and the packing arrangement of mPrP fibrils. Residues 165/167 are located on one end of the fibril core. Molecular dynamics simulations demonstrated that the stabilities of the seeding-induced ß-strand structures are significantly impacted by hydrogen bonds involving the side chain of residue 167 and steric resistance involving residue 165. These findings suggest that the α2-ß2 loop containing residues 165/167 could be the initial site of seed-template conformational conversion.

Structural Insight into the Mechanism of Staphylococcus aureus Stp1 Phosphatase.

Staphylococcus aureus Stp1, which belongs to the bacterial metal-dependent protein phosphatase (PPM) family, is a promising candidate for antivirulence targeting. How Stp1 recognizes the phosphorylated peptide remains unclear, however. In order to investigate the recognition mechanism of Stp1 in depth, we have determined a series of crystal structures of S. aureus Stp1 in different states and the structural complex of Stp1 bound with a phosphorylated peptide His12. Different phosphorylated peptides, including MgrA- and GraR-derived phosphopeptides, are substrates of Stp1, which supports the function of Stp1 as a selective Ser/Thr phosphatase. In addition, interestingly, the crystal structures of R161-Stp1 variants combined with the biochemical activity validations have uncovered that R161 residue plays a key role to control the conformation switches of the flap domain in order to facilitate substrate binding and the dephosphorylation process. Our findings provide crucial structural insight into the molecular mechanism of S. aureus Stp1 phosphatase and reveal the phosphorylated peptides for biochemistry study and inhibitor screening of Stp1.

Sensitive Water-Soluble Fluorescent Probe Based on Umpolung and Aggregation-Induced Emission Strategies for Selective Detection of Hg2+ in Living Cells and Zebrafish.

Using Hg2+-induced umpolung reaction and aggregation induced emission (AIE), we have rationally developed a water-soluble fluorescent probe 2,2'-(((4-(4,5-bis(4-methoxyphenyl)-1-phenyl-1H-imidazol-2-yl)phenyl)methylene)bis(sulfanediyl))diethanol (MPIPBS) for Hg2+ detection. MPIPBS was found to have high selectivity and sensitivity toward Hg2+ detection. The mechanism of MPIPBS response to Hg2+ was verified by 1H NMR titration, HPLC, and HRMS spectroscopy. The detection limit was examined to be 1.45 nM, which is lower than most reported probes for Hg2+. Taking advantage of excellent optical properties of AIEgen, a paper based sensor for Hg2+ detection was fabricated by immobilizing the MPIPBS on Waterman test paper. Meanwhile, MPIPBS showed satisfactory analytical performance in real water and urine samples. Further, thanks to the high water solubility, cell membrane permeability and low cytotoxicity, MPIPBS was further used to detect Hg2+ both in living cells and zebrafish. We anticipate that the prepared probe was available to detect Hg2+ in environment and biosamples.

The inhibitory mechanism of aurintricarboxylic acid targeting serine/threonine phosphatase Stp1 in Staphylococcus aureus: insights from molecular dynamics simulations.

Serine/threonine phosphatase (Stp1) is a member of the bacterial Mg2+- or Mn2+- dependent protein phosphatase/protein phosphatase 2C family, which is involved in the regulation of Staphylococcus aureus virulence. Aurintricarboxylic acid (ATA) is a known Stp1 inhibitor with an IC50 of 1.03 µM, but its inhibitory mechanism has not been elucidated in detail because the Stp1-ATA cocrystal structure has not been determined thus far. In this study, we performed 400 ns molecular dynamics (MD) simulations of the apo-Stp1 and Stp1-ATA complex models. During MD simulations, the flap subdomain of the Stp1-ATA complex experienced a clear conformational transition from an open state to a closed state, whereas the flap domain of apo-Stp1 changed from an open state to a semi-open state. In the Stp1-ATA complex model, the hydrogen bond (H-bond) between D137 and N142 disappeared, whereas critical H-bond interactions were formed between Q160 and H13, Q160/R161 and ATA, as well as N162 and D198. Finally, four residues (D137, N142, Q160, and R161) in Stp1 were mutated to alanine and the mutant enzymes were assessed using phosphate enzyme activity assays, which confirmed their important roles in maintaining Stp1 activity. This study indicated the inhibitory mechanism of ATA targeting Stp1 using MD simulations and sheds light on the future design of allosteric Stp1 inhibitors.

Protein encapsulation: a new approach for improving the capability of small-molecule fluorogenic probes.

Herein, we report a protein-based hybridization strategy that exploits the host-guest chemistry of HSA (human serum albumin) to solubilize the otherwise cell impermeable ONOO- fluorescent probe Pinkment-OAc. Formation of a HSA/Pinkment-OAc supramolecular hybrid was confirmed by SAXS and solution-state analyses. This HSA/Pinkment-OAc hybrid provided an enhanced fluorescence response towards ONOO- versus Pinkment-OAc alone, as determined by in vitro experiments. The HSA/Pinkment-OAc hybrid was also evaluated in RAW 264.7 macrophages and HeLa cancer cell lines, which displayed an enhanced cell permeability enabling the detection of SIN-1 and LPS generated ONOO- and the in vivo imaging of acute inflammation in LPS-treated mice. A remarkable 5.6 fold (RAW 264.7), 8.7-fold (HeLa) and 2.7-fold increased response was seen relative to Pinkment-OAc alone at the cellular level and in vivo, respectively. We anticipate that HSA/fluorescent probe hybrids will soon become ubiquitous and routinely applied to overcome solubility issues associated with hydrophobic fluorescent imaging agents designed to detect disease related biomarkers.

Biochemical Studies and Molecular Dynamic Simulations Reveal the Molecular Basis of Conformational Changes in DNA Methyltransferase-1.

DNA methyltransferase-1 (DNMT1) plays a crucial role in the maintenance of genomic methylation patterns. The crystal structure of DNMT1 was determined in two different states in which the helix that follows the catalytic loop was either kinked (designated helix-kinked) or well folded (designated helix-straight state). Here, we show that the proper structural transition between these two states is required for DNMT1 activity. The mutations of N1248A and R1279D, which did not affect interactions between DNMT1 and substrates or cofactors, allosterically reduced enzymatic activities in vitro by decreasing kcat/ Km for AdoMet. The crystallographic data combined with molecular dynamic (MD) simulations indicated that the N1248A and R1279D mutants bias the catalytic helix to either the kinked or straight conformation. In addition, genetic complementation assays for the two mutants suggested that disturbing the conformational transition reduced DNMT1 activity in cells, which could act additively with existing DNMT inhibitors to decrease DNA methylation. Collectively, our studies provide molecular insights into conformational changes of the catalytic helix, which is essential for DNMT1 catalytic activity, and thus aid in better understanding the relationship between DNMT1 dynamic switching and enzymatic activity.

Structural insights into the potency of SK channel positive modulators.

Small-conductance Ca2+-activated K+ (SK) channels play essential roles in the regulation of cellular excitability and have been implicated in neurological and cardiovascular diseases through both animal model studies and human genetic association studies. Over the past two decades, positive modulators of SK channels such as NS309 and 1-EBIO have been developed. Our previous structural studies have identified the binding pocket of 1-EBIO and NS309 that is located at the interface between the channel and calmodulin. In this study, we took advantage of four compounds with potencies varying over three orders of magnitude, including 1-EBIO, NS309, SKS-11 (6-bromo-5-methyl-1H-indole-2,3-dione-3-oxime) and SKS-14 (7-fluoro-3-(hydroxyimino)indolin-2-one). A combination of x-ray crystallographic, computational and electrophysiological approaches was utilized to investigate the interactions between the positive modulators and their binding pocket. A strong trend exists between the interaction energy of the compounds within their binding site calculated from the crystal structures, and the potency of these compounds in potentiating the SK2 channel current determined by electrophysiological recordings. Our results further reveal that the difference in potency of the positive modulators in potentiating SK2 channel activity may be attributed primarily to specific electrostatic interactions between the modulators and their binding pocket.

Conformation and dynamics of the C-terminal region in human phosphoglycerate mutase 1.

Phosphoglycerate mutase 1 (PGAM1), an important enzyme in glycolysis, is overexpressed in a number of human cancers, thus has been proposed as a promising metabolic target for cancer treatments. The C-terminal portion of the available crystal structures of PGAM1 and its homologous proteins is partially disordered, as evidenced by weak electron density. In this study, we identified the conformational behavior of the C-terminal region of PGAM1 as well as its role during the catalytic cycle. Using the PONDR-FIT server, we demonstrated that the C-terminal region was intrinsically disordered. We applied the Monte Carlo (MC) method to explore the conformational space of the C-terminus and conducted a series of explicit-solvent molecular dynamics (MD) simulations, and revealed that the C-terminal region is inherently dynamic; large-scale conformational changes in the C-terminal segment led to the structural transition of PGAM1 from the closed state to the open state. Furthermore, the C-terminal segment influenced 2,3-bisphosphoglycerate (2,3-BPG) binding. The proposed swing model illustrated a critical role of the C-terminus in the catalytic cycle through the conformational changes. In conclusion, the C-terminal region induces large movements of PGAM1 from the closed state to the open state and influences cofactor binding during the catalytic cycle. This report describes the dynamic features of the C-terminal region in detail and should aid in design of novel and efficient inhibitors of PGAM1. A swing mechanism of the C-terminal region is proposed, to facilitate further studies of the catalytic mechanism and the physiological functions of its homologues.