Nat Communications | Structure of human NaV1.6 channel reveals Na+ selectivity and pore blockade by 4,9-anhydro-tetrodotoxin

The sodium channel NaV1.6 is widely expressed in neurons of the central and peripheral nervous systems, which plays a critical role in regulating neuronal excitability. Dysfunction of NaV1.6 has been linked to epileptic encephalopathy, intellectual disability and movement disorders. Here we present cryo-EM structures of human NaV1.6/β1/β2 alone and complexed with a guanidinium neurotoxin 4,9-anhydro-tetrodotoxin (4,9-ah-TTX), revealing molecular mechanism of NaV1.6 inhibition by the blocker. The apo-form structure reveals two potential Na+ binding sites within the selectivity filter, suggesting a possible mechanism for Na+ selectivity and conductance. In the 4,9-ah-TTX bound structure, 4,9-ah-TTX binds to a pocket similar to the tetrodotoxin (TTX) binding site, which occupies the Na+ binding sites and completely blocks the channel. Molecular dynamics simulation results show that subtle conformational differences in the selectivity filter affect the affinity of TTX analogues. Taken together, our results provide important insights into NaV1.6 structure, ion conductance, and inhibition.

JCIM Cover | Observing Noncovalent Interactions in Experimental Electron Density for Macromolecular Systems: A Novel Perspective for Protein–Ligand Interaction Research

We report for the first time the use of experimental electron density (ED) in the Protein Data Bank for modeling of noncovalent interactions (NCIs) for protein–ligand complexes. Our methodology is based on reduced electron density gradient (RDG) theory describing intermolecular NCIs by ED and its first derivative. We established a database named Experimental NCI Database (ExptNCI; http://ncidatabase.stonewise.cn/#/nci) containing ED saddle points, indicating ∼200,000 NCIs from over 12,000 protein–ligand complexes. We also demonstrated the usage of the database in the case of depicting amide−π interactions in protein–ligand binding systems. In summary, the database provides details on experimentally observed NCIs for protein–ligand complexes and can support future studies including studies on rarely documented NCIs and the development of artificial intelligence models for protein–ligand binding prediction.

Cell | A mechanism for SARS-CoV-2 RNA capping and its inhibition by nucleotide analogue inhibitors

Decoration of cap on viral RNA plays essential roles in SARS-CoV-2 proliferation. Here we report a mechanism for SARS-CoV-2 RNA capping and document structural details at atomic resolution. The NiRAN domain in polymerase catalyzes the covalent link of RNA 5’ end to the first residue of nsp9 (termed as RNAylation), thus being an intermediate to form cap core (GpppA) with GTP catalyzed again by NiRAN. We also reveal that triphosphorylated nucleotide analogue inhibitors can be bonded to nsp9 and fit into a previously unknown ‘Nuc-pocket’ in NiRAN, thus inhibiting nsp9 RNAylation and formation of GpppA. S-loop (residues 50-KTN-52) in NiRAN presents a remarkable conformational shift observed in RTC bound with sofosbuvir monophosphate, reasoning an ‘induce-and-lock’ mechanism to design inhibitors. These findings not only improve the understanding of SARS-CoV-2 RNA capping and the mode of action of NAIs, but also provide a strategy to design antiviral drugs.

JCIM | Modified Electrostatic Complementary Score Function and Its Application Boundary Exploration in Drug Design

roperties with accuracy comparable to high-level quantum chemistry methods. One such example is the calculation of electrostatic potential (ESP). Different ESP prediction ML models were proposed to generate surface molecular charge distribution. Electrostatic complementarity (EC) can apply ESP data to quantify the complementarity between a ligand and its binding pocket, leading to the potential to increase the efficiency of drug design. However, there is not much research discussing EC score functions and their applicability domain. We propose a new EC score function modified from the one originally developed by Bauer and Mackey, and confirm its effectiveness against the available Pearson's R correlation coefficient. Additionally, the applicability domain of the EC score and two indices used to define the EC score application scope will be discussed.

Scientific Reports | A pocket-based 3D molecule generative model fueled by experimental electron density

We report for the first time the use of experimental electron density (ED) as training data for the generation of drug-like three-dimensional molecules based on the structure of a target protein pocket. Similar to a structural biologist building molecules based on their ED, our model functions with two main components: a generative adversarial network (GAN) to generate the ligand ED in the input pocket and an ED interpretation module for molecule generation. The model was tested on three targets: a kinase (hematopoietic progenitor kinase 1), protease (SARS-CoV-2 main protease), and nuclear receptor (vitamin D receptor), and evaluated with a reference dataset composed of over 8000 compounds that have their activities reported in the literature. The evaluation considered the chemical validity, chemical space distribution-based diversity, and similarity with reference active compounds concerning the molecular structure and pocket-binding mode. Our model can generate molecules with similar structures to classical active compounds and novel compounds sharing similar binding modes with active compounds, making it a promising tool for library generation supporting high-throughput virtual screening. The ligand ED generated can also be used to support fragment-based drug design. Our model is available as an online service to academic users via https://edmg.stonewise.cn/#/create .

Nat Commun. | N-type fast inactivation of a eukaryotic voltage-gated sodium channel

Voltage-gated sodium (NaV) channels initiate action potentials. Fast inactivation of NaV channels, mediated by an Ile-Phe-Met motif, is crucial for preventing hyperexcitability and regulating firing frequency. Here we present cryo-electron microscopy structure of NaVEh from the coccolithophore Emiliania huxleyi, which reveals an unexpected molecular gating mechanism for NaV channel fast inactivation independent of the Ile-Phe-Met motif. An N-terminal helix of NaVEh plugs into the open activation gate and blocks it. The binding pose of the helix is stabilized by multiple electrostatic interactions. Deletion of the helix or mutations blocking the electrostatic interactions completely abolished the fast inactivation. These strong interactions enable rapid inactivation, but also delay recovery from fast inactivation, which is ~160-fold slower than human NaV channels. Together, our results provide mechanistic insights into fast inactivation of NaVEh that fundamentally differs from the conventional local allosteric inhibition, revealing both surprising structural diversity and functional conservation of ion channel inactivation.

Nat Commun | Structure, gating, and pharmacology of human Cav3.3 channel

The low-voltage activated T-type calcium channels regulate cellular excitability and oscillatory behavior of resting membrane potential which trigger many physiological events and have been implicated with many diseases. Here, we determine structures of the human T-type CaV3.3 channel, in the absence and presence of antihypertensive drug mibefradil, antispasmodic drug otilonium bromide and antipsychotic drug pimozide. CaV3.3 contains a long bended S6 helix from domain III, with a positive charged region protruding into the cytosol, which is critical for T-type CaV channel activation at low voltage. The drug-bound structures clearly illustrate how these structurally different compounds bind to the same central cavity inside the CaV3.3 channel, but are mediated by significantly distinct interactions between drugs and their surrounding residues. Phospholipid molecules penetrate into the central cavity in various extent to shape the binding pocket and play important roles in stabilizing the inhibitor. These structures elucidate mechanisms of channel gating, drug recognition, and actions, thus pointing the way to developing potent and subtype-specific drug for therapeutic treatments of related disorders.

IEEE | Doc-to-Doc Recommender for Medical Literature with Similarity of Molecule Graphs

The increasing literature leads to formidable pressure for medical researchers. Most existing recommender approaches mainly depend on text-based information. How to extract and utilize the heterogeneous information, especially the graphic ones, to improve the recommender is worthy of further exploring. To this end, we establish a document-to-document recommender system for medical literature (D2D-MR). Specifically, we proposed HB-GED, the Half-branch GED algorithm, and the bipartite-graph-based algorithm for solving the molecule similarity and the paper similarity, respectively. Experimental results on real-world datasets demonstrate the effectiveness of the proposed recommender system.

JCIM | AIScaffold: A Web-Based Tool for Scaffold Diversification Using Deep Learning

AIScaffold , a web-based tool for scaffold diversification using the deep generative model, can perform large-scale (up to 500,000 molecules) diversification in several minutes and recommend the top 500 (top 0.1%) molecules. Features such as site-specific diversification are also supported. This tool can facilitate the scaffold diversification process for medicinal chemists, thereby accelerating drug design.

JCIM | DeepScaffold: A Comprehensive Tool for Scaffold-Based De Novo Drug Discovery Using Deep Learning

The ultimate goal of drug design is to find novel compounds with desirable pharmacological properties. Designing molecules retaining particular scaffolds as their core structures is an efficient way to obtain potential drug candidates. We propose a scaffold-based molecular generative model for drug discovery, which performs molecule generation based on a wide spectrum of scaffold definitions, including Bemis-Murcko scaffolds, cyclic skeletons, and scaffolds with specifications on side-chain properties. The model can generalize the learned chemical rules of adding atoms and bonds to a given scaffold. The generated compounds were evaluated by molecular docking in DRD2 targets, and the results demonstrated that this approach can be effectively applied to solve several drug design problems, including the generation of compounds containing a given scaffold and de novo drug design of potential drug candidates with specific docking scores.