D3Pockets: A Method and Web Server for Systematic Analysis of Protein Pocket Dynamics

The potential of Chlorella spp. as antiviral source against African swine fever virus through a virtual screening pipeline

African swine fever (ASF) causes high mortality in pigs and threatens global swine production. There is still a lack of therapeutics available, with two vaccines under scrutiny and no approved small-molecule drugs. Eleven (11) viral proteins were used to identify potential antivirals in in silico screening of secondary metabolites (127) from Chlorella spp. The metabolites were screened for affinity and binding selectivity. High-scoring compounds were assessed through in silico ADMET (Absorption, Distribution, Metabolism, Excretion, Toxicity) predictions, compared to structurally similar drugs, and checked for off-target docking with prepared swine receptors. Molecular dynamics (MD) simulations determined binding stability while binding energy was measured in Molecular Mechanics - Generalized Born Surface Area (MMGBSA) or Poisson-Boltzmann Surface Area (MMPBSA). Only six (6) compounds passed until MD analyses, of which five (5) were stable after 100 ns of MD runs. Of these five compounds, only three had binding affinities that were comparable to or stronger than controls. Specifically, phytosterols 24,25-dihydrolanosterol and CID 4206521 that interact with the RNA capping enzyme (pNP868R), and ergosterol which bound to the Erv-like thioreductase (pB119L). The compounds identified in this study can be used as a theoretical basis for in vitro screening to develop potent antiviral drugs against ASFV.

Mining the Dynamical Properties of Substrate and FAD Binding Pockets of LSD1: Hints for New Inhibitor Design Direction

Lysine-specific demethylase 1 (LSD1), a highly sophisticated epigenetic regulator, orchestrates a range of critical cellular processes, holding promising therapeutic potential for treating diverse diseases. However, the clinical research progress targeting LSD1 is very slow. After 20 years of research, only one small-molecule drug, BEA-17, targeting the degradation of LSD1 and CoREST has been approved by the U.S. Food and Drug Administration. The primary reason for this may be the lack of abundant structural data regarding its intricate functions. To gain a deeper understanding of its conformational dynamics and guide the drug design process, we conducted molecular dynamics simulations to explore the conformational states of LSD1 in the apo state and under the influence of cofactors of flavin adenine dinucleotide (FAD) and CoREST. Our results showed that, across all states, the substrate binding pocket exhibited high flexibility, whereas the FAD binding pocket remained more stable. These distinct dynamical properties are essential for LSD1's ability to bind various substrates while maintaining efficient demethylation activity. Both pockets can be enlarged by merging with adjacent pockets, although only the substrate binding pocket can shrink into smaller pockets. These new pocket shapes can inform inhibitor design, particularly for selectively FAD-competitive inhibitors of LSD1, given the presence of numerous FAD-dependent enzymes in the human body. More interestingly, in the absence of FAD binding, the united substrate and FAD binding pocket are partitioned by the conserved residue of Tyr761, offering valuable insights for the design of inhibitors that disrupt the crucial steric role of Tyr761 and the redox role of FAD. Additionally, we identified pockets that positively or negatively correlate with the substrate and FAD binding pockets, which can be exploited for the design of allosteric or concurrent inhibitors. Our results reveal the intricate dynamical properties of LSD1 as well as multiple novel conformational states, which deepen our understanding of its sophisticated functions and aid in the rational design of new inhibitors.

Subpocket-Based Analysis Approach for the Protein Pocket Dynamics

Structural and dynamic characteristics of protein pockets remarkably influence their biological functions and are also important for enzyme engineering and new drug research and development. To date, several softwares have been developed to analyze the dynamic properties of protein pockets. However, due to the complexity and diversity of the pocket information during the kinetic relaxation, further improvement and capacity expansion of current tools are required. Here, we developed a platform software AlphaTraj in which a computational strategy that divides the whole protein pocket into subpockets and examines various properties of the subpockets such as survival time, stability, and correlation was proposed and implemented. We also proposed a scoring function for the subpockets as well as the whole pocket to visualize the quality of the pocket. Furthermore, we implemented automated conformational search functions for ligand docking and ligand optimization. These functions may help us to gain a deep understanding of the dynamic properties of protein pockets and accelerate the protein engineering and the design of inhibitors and small-molecule drugs. The software is freely available at https://github.com/dooo12332/AlphaTraj.git under the GNU GPL license.

An Effective Plant Small Secretory Peptide Recognition Model Based on Feature Correction Strategy

Plant small secretory peptides (SSPs) play an important role in the regulation of biological processes in plants. Accurately predicting SSPs enables efficient exploration of their functions. Traditional experimental verification methods are very reliable and accurate, but they require expensive equipment and a lot of time. The method of machine learning speeds up the prediction process of SSPs, but the instability of feature extraction will also lead to further limitations of this type of method. Therefore, this paper proposes a new feature-correction-based model for SSP recognition in plants, abbreviated as SE-SSP. The model mainly includes the following three advantages: First, the use of transformer encoders can better reveal implicit features. Second, design a feature correction module suitable for sequences, named 2-D SENET, to adaptively adjust the features to obtain a more robust feature representation. Third, stack multiple linear modules to further dig out the deep information on the sample. At the same time, the training based on a contrastive learning strategy can alleviate the problem of sparse samples. We construct experiments on publicly available data sets, and the results verify that our model shows an excellent performance. The proposed model can be used as a convenient and effective SSP prediction tool in the future. Our data and code are publicly available at https://github.com/wrab12/SE-SSP/.

Exploring biased activation characteristics by molecular dynamics simulation and machine learning for the μ-opioid receptor

Biased ligands selectively activating specific downstream signaling pathways (termed as biased activation) exhibit significant therapeutic potential. However, the conformational characteristics revealed are very limited for the biased activation, which is not conducive to biased drug development. Motivated by the issue, we combine extensive accelerated molecular dynamics simulations and an interpretable deep learning model to probe the biased activation features for two complex systems constructed by the inactive μOR and two different biased agonists (G-protein-biased agonist TRV130 and β-arrestin-biased agonist endomorphin2). The results indicate that TRV130 binds deeper into the receptor core compared to endomorphin2, located between W2936.48 and D1142.50, and forms hydrogen bonding with D1142.50, while endomorphin2 binds above W2936.48. The G protein-biased agonist induces greater outward movements of the TM6 intracellular end, forming a typical active conformation, while the β-arrestin-biased agonist leads to a smaller extent of outward movements of TM6. Compared with TRV130, endomorphin2 causes more pronounced inward movements of the TM7 intracellular end and more complex conformational changes of H8 and ICL1. In addition, important residues determining the two different biased activation states were further identified by using an interpretable deep learning classification model, including some common biased activation residues across Class A GPCRs like some key residues on the TM2 extracellular end, ECL2, TM5 intracellular end, TM6 intracellular end, and TM7 intracellular end, and some specific important residues of ICL3 for μOR. The observations will provide valuable information for understanding the biased activation mechanism for GPCRs.

Computational Chemistry in Structure-Based Solute Carrier Transporter Drug Design: Recent Advances and Future Perspectives

Solute carrier transporters (SLCs) are a class of important transmembrane proteins that are involved in the transportation of diverse solute ions and small molecules into cells. There are approximately 450 SLCs within the human body, and more than a quarter of them are emerging as attractive therapeutic targets for multiple complex diseases, e.g., depression, cancer, and diabetes. However, only 44 unique transporters (∼9.8% of the SLC superfamily) with 3D structures and specific binding sites have been reported. To design innovative and effective drugs targeting diverse SLCs, there are a number of obstacles that need to be overcome. However, computational chemistry, including physics-based molecular modeling and machine learning- and deep learning-based artificial intelligence (AI), provides an alternative and complementary way to the classical drug discovery approach. Here, we present a comprehensive overview on recent advances and existing challenges of the computational techniques in structure-based drug design of SLCs from three main aspects: (i) characterizing multiple conformations of the proteins during the functional process of transportation, (ii) identifying druggability sites especially the cryptic allosteric ones on the transporters for substrates and drugs binding, and (iii) discovering diverse small molecules or synthetic protein binders targeting the binding sites. This work is expected to provide guidelines for a deep understanding of the structure and function of the SLC superfamily to facilitate rational design of novel modulators of the transporters with the aid of state-of-the-art computational chemistry technologies including artificial intelligence.

CETSA-MS-based target profiling of anti-aging natural compound quercetin

BACKGROUND

Quercetin is widely distributed in nature and abundant in the human diet, which exhibits diverse biological activities and potential medical benefits. However, there remains a lack of comprehensive understanding about its cellular targets, impeding its in-depth mechanistic studies and clinical applications.

PURPOSE

This study aimed to profile protein targets of quercetin at the proteome level.

METHODS

A label-free CETSA-MS proteomics technique was employed for target enrichment and identification. The R package Inflect was used for melting curve fitting and target selection. D3Pocket and LiBiSco tools were used for binding pocket prediction and binding pocket analysis. Western blotting, molecular docking, site-directed mutagenesis and pull-down assays were used for target verification and validation.

RESULTS

We curated a library of direct binding targets of quercetin in cells. This library comprises 37 proteins that show increased thermal stability upon quercetin binding and 33 proteins that display decreased thermal stability. Through Western blotting, molecular docking, site-directed mutagenesis and pull-down assays, we validated CBR1 and GSK3A from the stabilized protein group and MAPK1 from the destabilized group as direct binding targets of quercetin. Moreover, we characterized the shared chemical properties of the binding pockets of quercetin with targets.

CONCLUSION

Our findings deepen our understanding of the proteins pivotal to the bioactivity of quercetin and lay the groundwork for further exploration into its mechanisms of action and potential clinical applications.

Mechanism of Mutation-Induced Effects on the Catalytic Function of TEV Protease: A Molecular Dynamics Study

Tobacco etch virus protease (TEVp) is wildly exploited for various biotechnological applications. These applications take advantage of TEVp's ability to cleave specific substrate sequences to study protein function and interactions. A major limitation of this enzyme is its relatively slow catalytic rate. In this study, MD simulations were conducted on TEV enzymes and known highly active mutants (eTEV and uTEV3) to explore the relationship between mutation, conformation, and catalytic function. The results suggest that mutations distant from the active site can influence the substrate-binding pocket through interaction networks. MD analysis of eTEV demonstrates that, by stabilizing the orientation of the substrate at the catalytic site, mutations that appropriately enlarge the substrate-binding pocket will be beneficial for Kcat, enhancing the catalytic efficiency of the enzyme. On the contrary, mutations in uTEV3 reduced the flexibility of the active pocket and increased the hydrogen bonding between the substrate and enzyme, resulting in higher affinity. At the same time, the MD simulation demonstrates that mutations outside of the active site residues could affect the dynamic movement of the binding pocket by altering residue networks and communication pathways, thereby having a profound impact on reactivity. These findings not only provide a molecular mechanistic explanation for the excellent mutants, but also serve as a guiding framework for rational computational design.

A Structure-Based Allosteric Modulator Design Paradigm

Importance: Allosteric drugs bound to topologically distal allosteric sites hold a substantial promise in modulating therapeutic targets deemed undruggable at their orthosteric sites. Traditionally, allosteric modulator discovery has predominantly relied on serendipitous high-throughput screening. Nevertheless, the landscape has undergone a transformative shift due to recent advancements in our understanding of allosteric modulation mechanisms, coupled with a significant increase in the accessibility of allosteric structural data. These factors have extensively promoted the development of various computational methodologies, especially for machine-learning approaches, to guide the rational design of structure-based allosteric modulators. Highlights: We here presented a comprehensive structure-based allosteric modulator design paradigm encompassing 3 critical stages: drug target acquisition, allosteric binding site, and modulator discovery. The recent advances in computational methods in each stage are encapsulated. Furthermore, we delve into analyzing the successes and obstacles encountered in the rational design of allosteric modulators. Conclusion: The structure-based allosteric modulator design paradigm holds immense potential for the rational design of allosteric modulators. We hope that this review would heighten awareness of the use of structure-based computational methodologies in advancing the field of allosteric drug discovery.

Progressive computational approaches to facilitate decryption of allosteric mechanism and drug discovery

Allostery is a ubiquitous biological phenomenon where perturbation at topologically distal areas of a protein serves as a trigger to fine-tune the orthosteric site and thus regulate protein function. The investigation of allosteric regulation greatly enhances our understanding of human diseases and broadens avenue for drug discovery. For decades, owing to the difficulty in allostery characterization through serendipitous experimental screening, researchers have developed several innovative computational approaches, which proves to accelerate the elucidation of allostery. Herein, we review the state-of-the-art advance of computational methodologies for allostery study, with particular emphasis on promising trends emerging over the past two years. We expect this review will outline the latest landscape of allostery study and inspire researchers to further facilitate this field.

GPCR molecular dynamics forecasting using recurrent neural networks

G protein-coupled receptors (GPCRs) are a large superfamily of cell membrane proteins that play an important physiological role as transmitters of extracellular signals. Signal transmission through the cell membrane depends on conformational changes in the transmembrane region of the receptor, which makes the investigation of the dynamics in these regions particularly relevant. Molecular dynamics (MD) simulations provide a wealth of data about the structure, dynamics, and physiological function of biological macromolecules by modelling the interactions between their atomic constituents. In this study, a Recurrent and Convolutional Neural Network (RNN) model, namely Long Short-Term Memory (LSTM), is used to predict the dynamics of two GPCR states and three specific simulations of each one, through their activation path and focussing on specific receptor regions. Active and inactive states of the GPCRs are analysed in six scenarios involving APO, Full Agonist (BI 167107) and Partial Inverse Agonist (carazolol) of the receptor. Four Machine Learning models with increasing complexity in terms of neural network architecture are evaluated, and their results discussed. The best method achieves an overall RMSD lower than 0.139 Å and the transmembrane helices are the regions showing the minimum prediction errors and minimum relative movements of the protein.

Repurposing of rilpivirine for preventing platelet β3 integrin-dependent thrombosis by targeting c-Src active autophosphorylation

BACKGROUND

HIV-infected individuals are known to be at higher risk for thrombotic cardiovascular disease (CVD), which may also be differentially affected by components of anti-HIV drugs. To identify the effects of a series of FDA-approved anti-HIV drugs on platelet aggregation in humans, focusing on the novel pharmacological effects of rilpivirine (RPV), a reverse transcriptase inhibitor, on platelet function both in vitro and in vivo and the mechanisms involved.

METHODS AND RESULTS

In vitro studies showed that RPV was the only anti-HIV reagent that consistently and efficiently inhibited aggregation elicited by different agonists, exocytosis, morphological extension on fibrinogen, and clot retraction. Treatment of mice with RPV significantly prevented thrombus formation in FeCl3-injured mesenteric vessels, postcava with stenosis surgery, and ADP -induced pulmonary embolism models without defects in platelet viability, tail bleeding, and coagulation activities. RPV also improved cardiac performance in mice with post-ischemic reperfusion. A mechanistic study revealed that RPV preferentially attenuated fibrinogen-stimulated Tyr773 phosphorylation of β3-integrin by inhibiting Tyr419 autophosphorylation of c-Src. Molecular docking and surface plasmon resonance analyses showed that RPV can bind directly to c-Src. Further mutational analysis showed that the Phe427 residue of c-Src is critical for RPV interaction, suggesting a novel interaction site for targeting c-Src to block β3-integrin outside-in signaling.

CONCLUSION

These results demonstrated that RPV was able to prevent the progression of thrombotic CVDs by interrupting β3-integrin-mediated outside-in signaling via inhibiting c-Src activation without hemorrhagic side effects, highlighting RPV as a promising reagent for the prevention and therapy of thrombotic CVDs.

Insights into the Binding Profile of Anti-chlorpyrifos Recombinant Antibodies: From Computational Simulation to Immunoassay Validation

A novel virtual screening strategy was proposed for the profiling and discovery of active variable regions (VRs) that encode hapten-specific recombinant antibodies (rAbs). Chlorpyrifos, a hazardous organophosphorus pesticide, was selected as the target. First, a VR model-14G4 from anti-chlorpyrifos hybridoma was built via homology modeling. Its binding pattern toward seven organophosphorus analogues was assessed through virtual screening by performing molecular docking. Based on energy scoring, visual examination, and molecular interaction analysis, chlorpyrifos-methyl was also inferred as the high-affinity target for model-14G4 and was then confirmed via a non-competitive surface plasmon resonance (SPR) assay. Subsequently, we attempted to discover hapten-specific VRs by creating a collection of VR models for anonymous testing. Chlorpyrifos and model-14G4 were employed as the known hit and active VRs, respectively. After molecular docking, a novel anti-chlorpyrifos VR (model-1) was identified due to its satisfactory energy scoring and a similar binding pattern to the reference model-14G4. Expressed by HEK293(F) mammalian cells, the newly prepared full-length rAb-model-1 and rAb-14G4 exhibited high sensitivities for detecting chlorpyrifos by the indirect competitive enzyme-linked immunosorbent assay (ic-ELISA), with IC₅₀ of 3.01 ng/mL and 42.82 ng/mL, respectively. They recognized chlorpyrifos-methyl with a cross-reactivity (CR) of 2.5-17.3%. Moreover, the binding properties of rAb-model-1 for recognizing chlorpyrifos and chlorpyrifos-methyl were confirmed via a non-competitive microscale thermophoresis (MST) method. Thus, the experimental results showed good agreement with computational outputs on antibody profiling. Furthermore, the recognition diversity of rAb-model-1 for chlorpyrifos and chlorpyrifos-methyl was studied via molecular dynamics simulation. Overall, the proposed study provides a versatile and economical strategy for antibody characterization and promotes the in vitro production of rAbs for pesticide monitoring.

Whole-Genome Resequencing Reveals the Diversity of Patchouli Germplasm

As an important medicinal and aromatic plant, patchouli is distributed throughout most of Asia. However, current research on patchouli's genetic diversity is limited and lacks genome-wide studies. Here, we have collected seven representative patchouli accessions from different localities and performed whole-genome resequencing on them. In total, 402,650 single nucleotide polymorphisms (SNPs) and 153,233 insertions/deletions (INDELs) were detected. Based on these abundant genetic variants, patchouli accessions were primarily classified into the Chinese group and the Southeast Asian group. However, the accession SP (Shipai) collected from China formed a distinct subgroup within the Southeast Asian group. As SP has been used as a genuine herb in traditional Chinese medicine, its unique molecular markers have been subsequently screened and verified. For 26,144 specific SNPs and 16,289 specific INDELs in SP, 10 of them were validated using Polymerase Chain Reaction (PCR) following three different approaches. Further, we analyzed the effects of total genetic variants on genes involved in the sesquiterpene synthesis pathway, which produce the primary phytochemical compounds found in patchouli. Eight genes were ultimately investigated and a gene encoding nerolidol synthetase (PatNES) was chosen and confirmed through biochemical assay. In accession YN, genetic variants in PatNES led to a loss of synthetase activity. Our results provide valuable information for understanding the diversity of patchouli germplasm resources.

Identification Mechanism of BACE1 on Inhibitors Probed by Using Multiple Separate Molecular Dynamics Simulations and Comparative Calculations of Binding Free Energies

β-amyloid cleaving enzyme 1 (BACE1) is regarded as an important target of drug design toward the treatment of Alzheimer's disease (AD). In this study, three separate molecular dynamics (MD) simulations and calculations of binding free energies were carried out to comparatively determine the identification mechanism of BACE1 for three inhibitors, 60W, 954 and 60X. The analyses of MD trajectories indicated that the presence of three inhibitors influences the structural stability, flexibility and internal dynamics of BACE1. Binding free energies calculated by using solvated interaction energy (SIE) and molecular mechanics generalized Born surface area (MM-GBSA) methods reveal that the hydrophobic interactions provide decisive forces for inhibitor-BACE1 binding. The calculations of residue-based free energy decomposition suggest that the sidechains of residues L91, D93, S96, V130, Q134, W137, F169 and I179 play key roles in inhibitor-BACE1 binding, which provides a direction for future drug design toward the treatment of AD.

Butein inhibits cancer cell growth by rescuing the wild-type thermal stability of mutant p53

p53 is a transcription factor that activates the expression of various genes involved in the maintenance of genomic stability, and more than 50% of cancers harbor inactivating p53 mutations, which are indicative of highly aggressive cancer and poor prognosis. Pharmacological targeting of mutant p53 to restore the wild-type p53 tumor-suppressing function is a promising strategy for cancer therapy. In this study, we identified a small molecule, Butein, that reactivates mutant p53 activity in tumor cells harboring the R175H or R273H mutation. Butein restored wild-type-like conformation and DNA-binding ability in HT29 and SK-BR-3 cells harboring mutant p53-R175H and mutant p53-R273H, respectively. Moreover, Butein enabled the transactivation of p53 target genes and decreased the interactions of Hsp90 with mutant p53-R175H and mutant p53-R273H proteins, while Hsp90 overexpression reversed targeted p53 gene activation. In addition, Butein induced thermal stabilization of wild-type p53, mutant p53-R273H and mutant p53-R175H, as determined via CETSA. From docking study, we further proved that Butein binding to p53 stabilized the DNA-binding loop-sheet-helix motif of mutant p53-R175H and regulated its DNA-binding activity via an allosteric mechanism, conferring wild-type-like the DNA-binding activity of mutant p53. Collectively, the data suggest that Butein is a potential antitumor agent that restores p53 function in cancers harboring mutant p53-R273H or mutant p53-R175H. SIGNIFICANCE: Butein restores the ability of mutant p53 to bind DNA by reversing its transition to the Loop3 (L3) state, endows p53 mutants with thermal stability and re-establishes their transcriptional activity to induce cancer cell death.

The binding pocket properties were fundamental to functional diversification of the GDSL-type esterases/lipases gene family in cotton

Cotton is one of the most important crops in the world. GDSL-type esterases/lipases (GELPs) are widely present in all kingdoms and play an essential role in regulating plant growth, development, and responses to abiotic and biotic stresses. However, the molecular mechanisms underlying this functional diversity remain unclear. Here, based on the identification of the GELP gene family, we applied genetic evolution and molecular simulation techniques to explore molecular mechanisms in cotton species. A total of 1502 GELP genes were identified in 10 cotton species. Segmental duplication and differences in evolutionary rates are the leading causes of the increase in the number and diversity of GELP genes during evolution for ecological adaptation. Structural analysis revealed that the GELP family has high structural diversity. Moreover, molecular simulation studies have demonstrated significant differences in the properties of the binding pockets among cotton GELPs. In the process of adapting to the environment, GELPs not only have segmental duplication but also have different evolutionary rates, resulting in gene diversity. This diversity leads to significant differences in the 3D structure and binding pocket properties and, finally, to functional diversity. These findings provide a reference for further functional analyses of plant GELPs.

Layer-Wise Relevance Analysis for Motif Recognition in the Activation Pathway of the β2-Adrenergic GPCR Receptor

G-protein-coupled receptors (GPCRs) are cell membrane proteins of relevance as therapeutic targets, and are associated to the development of treatments for illnesses such as diabetes, Alzheimer's, or even cancer. Therefore, comprehending the underlying mechanisms of the receptor functional properties is of particular interest in pharmacoproteomics and in disease therapy at large. Their interaction with ligands elicits multiple molecular rearrangements all along their structure, inducing activation pathways that distinctly influence the cell response. In this work, we studied GPCR signaling pathways from molecular dynamics simulations as they provide rich information about the dynamic nature of the receptors. We focused on studying the molecular properties of the receptors using deep-learning-based methods. In particular, we designed and trained a one-dimensional convolution neural network and illustrated its use in a classification of conformational states: active, intermediate, or inactive, of the β2-adrenergic receptor when bound to the full agonist BI-167107. Through a novel explainability-oriented investigation of the prediction results, we were able to identify and assess the contribution of individual motifs (residues) influencing a particular activation pathway. Consequently, we contribute a methodology that assists in the elucidation of the underlying mechanisms of receptor activation-deactivation.

Integrating transcriptomics and network analysis-based multiplexed drug repurposing to screen drug candidates for M2 macrophage-associated castration-resistant prostate cancer bone metastases

Metastatic castration-resistant prostate cancer (CRPC) has long been considered to be associated with patient mortality. Among metastatic organs, bone is the most common metastatic site, with more than 90% of advanced patients developing bone metastases (BMs) before 24 months of death. Although patients were recommended to use bone-targeted drugs represented by bisphosphonates to treat BMs of CRPC, there was no significant improvement in patient survival. In addition, the use of immunotherapy and androgen deprivation therapy is limited due to the immunosuppressed state and resistance to antiandrogen agents in patients with bone metastases. Therefore, it is still essential to develop a safe and effective therapeutic schedule for CRPC patients with BMs. To this end, we propose a multiplex drug repurposing scheme targeting differences in patient immune cell composition. The identified drug candidates were ranked from the perspective of M2 macrophages by integrating transcriptome and network-based analysis. Meanwhile, computational chemistry and clinical trials were used to generate a comprehensive drug candidate list for the BMs of CRPC by drug redundancy structure filtering. In addition to docetaxel, which has been approved for clinical trials, the list includes norethindrone, testosterone, menthol and foretinib. This study provides a new scheme for BMs of CRPC from the perspective of M2 macrophages. It is undeniable that this multiplex drug repurposing scheme specifically for immune cell-related bone metastases can be used for drug screening of any immune-related disease, helping clinicians find promising therapeutic schedules more quickly, and providing reference information for drug R&D and clinical trials.

In silico prediction of the animal susceptibility and virtual screening of natural compounds against SARS-CoV-2: Molecular dynamics simulation based analysis

COVID-19 is an infectious disease caused by the SARS-CoV-2 virus. It has six open reading frames (orf1ab, orf3a, orf6, orf7a, orf8, and orf10), a spike protein, a membrane protein, an envelope small membrane protein, and a nucleocapsid protein, out of which, orf1ab is the largest ORF coding different important non-structural proteins. In this study, an effort was made to evaluate the susceptibility of different animals against SARS-CoV-2 by analyzing the interactions of Spike and ACE2 proteins of the animals and propose a list of potential natural compounds binding to orf1ab of SARS-CoV-2. Here, we analyzed structural interactions between spike proteins of SARS-CoV-2 and the ACE2 receptor of 16 different hosts. A simulation for 50 ns was performed on these complexes. Based on post-simulation analysis, Chelonia mydas was found to have a more stable complex, while Bubalus bubalis, Aquila chrysaetos chrysaetos, Crocodylus porosus, and Loxodonta africana were found to have the least stable complexes with more fluctuations than all other organisms. Apart from that, we performed domain assignment of orf1ab of SARS-CoV-2 and identified 14 distinct domains. Out of these, Domain 3 (DNA/RNA polymerases) was selected as a target, as it showed no similarities with host proteomes and was validated in silico. Then, the top 10 molecules were selected from the virtual screening of ∼1.8 lakh molecules from the ZINC database, based on binding energy, and validated for ADME and toxicological properties. Three molecules were selected and analyzed further. The structural analysis showed that these molecules were residing within the pocket of the receptor. Finally, a simulation for 200 ns was performed on complexes with three selected molecules. Based on post-simulation analysis (RMSD, RMSF, Rg, SASA, and energies), the molecule ZINC000103666966 was found as the most suitable inhibitory compound against Domain 3. As this is an in silico prediction, further experimental studies could unravel the potential of the proposed molecule against SARS-CoV-2.

Structures of the SARS-CoV-2 spike glycoprotein and applications for novel drug development

COVID-19 caused by SARS-CoV-2 has raised a health crisis worldwide. The high morbidity and mortality associated with COVID-19 and the lack of effective drugs or vaccines for SARS-CoV-2 emphasize the urgent need for standard treatment and prophylaxis of COVID-19. The receptor-binding domain (RBD) of the glycosylated spike protein (S protein) is capable of binding to human angiotensin-converting enzyme 2 (hACE2) and initiating membrane fusion and virus entry. Hence, it is rational to inhibit the RBD activity of the S protein by blocking the RBD interaction with hACE2, which makes the glycosylated S protein a potential target for designing and developing antiviral agents. In this study, the molecular features of the S protein of SARS-CoV-2 are highlighted, such as the structures, functions, and interactions of the S protein and ACE2. Additionally, computational tools developed for the treatment of COVID-19 are provided, for example, algorithms, databases, and relevant programs. Finally, recent advances in the novel development of antivirals against the S protein are summarized, including screening of natural products, drug repurposing and rational design. This study is expected to provide novel insights for the efficient discovery of promising drug candidates against the S protein and contribute to the development of broad-spectrum anti-coronavirus drugs to fight against SARS-CoV-2.

Mechanistic Insights into the Long-range Allosteric Regulation of KRAS Via Neurofibromatosis Type 1 (NF1) Scaffold Upon SPRED1 Loading

Allosteric regulation is the most direct and efficient way of regulating protein function, wherein proteins transmit the perturbations at one site to another distinct functional site. Deciphering the mechanism of allosteric regulation is of vital importance for the comprehension of both physiological and pathological events in vivo as well as the rational allosteric drug design. However, it remains challenging to elucidate dominant allosteric signal transduction pathways, especially for large and multi-component protein machineries where long-range allosteric regulation exits. One of the quintessential examples having long-range allosteric regulation is the ternary complex, SPRED1-RAS-neurofibromin type 1 (NF1, a RAS GTPase-activating protein), in which SPRED1 facilitates RAS-GTP hydrolysis by interacting with NF1 at a distal, allosteric site from the RAS binding site. To address the underlying mechanism, we performed extensive Gaussian accelerated molecular dynamics simulations and Markov state model analysis of KRAS-NF1 complex in the presence and absence of SPRED1. Our findings suggested that SPRED1 loading allosterically enhanced KRAS-NF1 binding, but hindered conformational transformation of the NF1 catalytic center for RAS hydrolysis. Moreover, we unveiled the possible allosteric pathways upon SPRED1 binding through difference contact network analysis. This study not only provided an in-depth mechanistic insight into the allosteric regulation of KRAS by SPRED1, but also shed light on the investigation of long-range allosteric regulation among complex macromolecular systems.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

Conserved protein targets for developing pan-coronavirus drugs based on sequence and 3D structure similarity analyses

There are 7 known human pathogenic coronaviruses, which are HCoV-229E, HCoV-OC43, HCoV-NL63, HCoV-HKU1, MERS-CoV, SARS-CoV and SARS-CoV-2. While SARS-CoV-2 is currently caused a severe epidemic, experts believe that new pathogenic coronavirus would emerge in the future. Therefore, developing broad-spectrum anti-coronavirus drugs is of great significance. In this study, we performed protein sequence and three-dimensional structure analyses for all the 20 virus-encoded proteins across all the 7 coronaviruses, with the purpose to identify highly conserved proteins and binding sites for developing pan-coronavirus drugs. We found that nsp5, nsp10, nsp12, nsp13, nsp14, and nsp16 are highly conserved both in protein sequences (with average identity percentage higher than 52%, average amino acid conservation scores higher than 5.2) and binding pockets (with average amino acid conservation scores higher than 5.8). We also performed the similarity comparison between these 6 proteins and all the human proteins, and found that all the 6 proteins have similarity less than 25%, indicating that the drugs targeting the 6 proteins should have little interference of human protein function. Accordingly, we suggest that nsp5, nsp10, nsp12, nsp13, nsp14, and nsp16 are potential targets for pan-coronavirus drug development.

D3AI-CoV: a deep learning platform for predicting drug targets and for virtual screening against COVID-19

Target prediction and virtual screening are two powerful tools of computer-aided drug design. Target identification is of great significance for hit discovery, lead optimization, drug repurposing and elucidation of the mechanism. Virtual screening can improve the hit rate of drug screening to shorten the cycle of drug discovery and development. Therefore, target prediction and virtual screening are of great importance for developing highly effective drugs against COVID-19. Here we present D3AI-CoV, a platform for target prediction and virtual screening for the discovery of anti-COVID-19 drugs. The platform is composed of three newly developed deep learning-based models i.e., MultiDTI, MPNNs-CNN and MPNNs-CNN-R models. To compare the predictive performance of D3AI-CoV with other methods, an external test set, named Test-78, was prepared, which consists of 39 newly published independent active compounds and 39 inactive compounds from DrugBank. For target prediction, the areas under the receiver operating characteristic curves (AUCs) of MultiDTI and MPNNs-CNN models are 0.93 and 0.91, respectively, whereas the AUCs of the other reported approaches range from 0.51 to 0.74. For virtual screening, the hit rate of D3AI-CoV is also better than other methods. D3AI-CoV is available for free as a web application at http://www.d3pharma.com/D3Targets-2019-nCoV/D3AI-CoV/index.php, which can serve as a rapid online tool for predicting potential targets for active compounds and for identifying active molecules against a specific target protein for COVID-19 treatment.

An Interpretable Convolutional Neural Network Framework for Analyzing Molecular Dynamics Trajectories: a Case Study on Functional States for G-Protein-Coupled Receptors

Molecular dynamics (MD) simulations have made great contribution to revealing structural and functional mechanisms for many biomolecular systems. However, how to identify functional states and important residues from vast conformation space generated by MD remains challenging; thus an intelligent navigation is highly desired. Despite intelligent advantages of deep learning exhibited in analyzing MD trajectory, its black-box nature limits its application. To address this problem, we explore an interpretable convolutional neural network (CNN)-based deep learning framework to automatically identify diverse active states from the MD trajectory for G-protein-coupled receptors (GPCRs), named the ICNNMD model. To avoid the information loss in representing the conformation structure, the pixel representation is introduced, and then the CNN module is constructed to efficiently extract features followed by a fully connected neural network to realize the classification task. More importantly, we design a local interpretable model-agnostic explanation interpreter for the classification result by local approximation with a linear model, through which important residues underlying distinct active states can be quickly identified. Our model showcases higher than 99% classification accuracy for three important GPCR systems with diverse active states. Notably, some important residues in regulating different biased activities are successfully identified, which are beneficial to elucidating diverse activation mechanisms for GPCRs. Our model can also serve as a general tool to analyze MD trajectory for other biomolecular systems. All source codes are freely available at https://github.com/Jane-Liu97/ICNNMD for aiding MD studies.

Identification of a Potential mRNA-based Vaccine Candidate against the SARS-CoV-2 Spike Glycoprotein: A Reverse Vaccinology Approach

The emergence of the novel coronavirus (SARS-CoV-2) in December 2019 has generated a devastating global consequence which makes the development of a rapidly deployable, effective and safe vaccine candidate an imminent global health priority. The design of most vaccine candidates has been directed at the induction of antibody responses against the trimeric spike glycoprotein of SARS-CoV-2, a class I fusion protein that aids ACE2 (angiotensin-converting enzyme 2) receptor binding. A variety of formulations and vaccinology approaches are being pursued for targeting the spike glycoprotein, including simian and human replication-defective adenoviral vaccines, subunit protein vaccines, nucleic acid vaccines and whole-inactivated SARS-CoV-2. Here, we directed a reverse vaccinology approach towards the design of a nucleic acid (mRNA-based) vaccine candidate. The "YLQPRTFLL" peptide sequence (position 269-277) which was predicted to be a B cell epitope and likewise a strong binder of the HLA*A-0201 was selected for the design of the vaccine candidate, having satisfied series of antigenicity assessments. Through the codon optimization protocol, the nucleotide sequence for the vaccine candidate design was generated and targeted at the human toll-like receptor 7 (TLR7). Bioinformatics analyses showed that the sequence "UACCUGCAGCCGCGUACCUUCCUGCUG" exhibited a strong affinity and likewise was bound to a stable cavity in the TLR7 pocket. This study is therefore expected to contribute to the research efforts directed at securing definitive preventive measures against the SARS-CoV-2 infection.

Harnessing Reversed Allosteric Communication: A Novel Strategy for Allosteric Drug Discovery

Allostery is a fundamental and extensive mechanism of intramolecular signal transmission. Allosteric drugs possess several unique pharmacological advantages over traditional orthosteric drugs, including greater selectivity, better physicochemical properties, and lower off-target toxicity. However, owing to the complexity of allosteric regulation, experimental approaches for the development of allosteric modulators are traditionally serendipitous. Recently, the reversed allosteric communication theory has been proposed, providing a feasible tool for the unbiased detection of allosteric sites. Herein, we review the latest research on the reversed allosteric communication effect using the examples of sirtuin 6, epidermal growth factor receptor, 3-phosphoinositide-dependent protein kinase 1, and Related to A and C kinases (RAC) serine/threonine protein kinase B and recapitulate the methodologies of reversed allosteric communication strategy. The novel reversed allosteric communication strategy greatly expands the horizon of allosteric site identification and allosteric mechanism exploration and is expected to accelerate an end-to-end framework for drug discovery.

Multiple regulation and targeting effects of borneol in the neurovascular unit in neurodegenerative diseases

Efficient delivery of brain-targeted drugs is highly important for the success of therapies in neurodegenerative diseases. Borneol has several biological activities, such as anti-inflammatory and cell penetration enhancing effect, and can regulate processes in the neurovascular unit (NVU), such as protein toxic stress, autophagosome/lysosomal system, oxidative stress, programmed cell death and neuroinflammation. However, the influence of borneol on NVU in neurodegenerative diseases has not been fully explained. This study searched the keywords 'borneol', 'neurovascular unit', 'endothelial cell', 'astrocyte', 'neuron', 'blood-brain barrier', 'neurodegenerative diseases' and 'brain disease', in PubMed, BioMed Central, China National Knowledge Infrastructure (CNKI), and Bing search engines to explore the influence of borneol on NVU. In addition to the principle and mechanism of penetration of borneol in the brain, this study also showed its multiple regulation effects on NVU. Borneol was able to penetrate the blood-brain barrier (BBB), affecting the signal transmission between BBB and the microenvironment of the brain, down-regulating the expression of inflammatory and oxidative stress proteins in NVU, especially in microglia and astrocytes. In summary, borneol is a potential drug delivery agent for drugs against neurodegenerative diseases.

Finding Druggable Sites in Proteins Using TACTICS

Structure-based drug discovery efforts require knowledge of where drug-binding sites are located on target proteins. To address the challenge of finding druggable sites, we developed a machine-learning algorithm called TACTICS (trajectory-based analysis of conformations to identify cryptic sites), which uses an ensemble of molecular structures (such as molecular dynamics simulation data) as input. First, TACTICS uses k-means clustering to select a small number of conformations that represent the overall conformational heterogeneity of the data. Then, TACTICS uses a random forest model to identify potentially bindable residues in each selected conformation, based on protein motion and geometry. Lastly, residues in possible binding pockets are scored using fragment docking. As proof-of-principle, TACTICS was applied to the analysis of simulations of the SARS-CoV-2 main protease and methyltransferase and the Yersinia pestis aryl carrier protein. Our approach recapitulates known small-molecule binding sites and predicts the locations of sites not previously observed in experimentally determined structures. The TACTICS code is available at https://github.com/Albert-Lau-Lab/tactics_protein_analysis.

Untangling Dual-Targeting Therapeutic Mechanism of Epidermal Growth Factor Receptor (EGFR) Based on Reversed Allosteric Communication

Dual-targeting therapeutics by coadministration of allosteric and orthosteric drugs is drawing increased attention as a revolutionary strategy for overcoming the drug-resistance problems. It was further observed that the occupation of orthosteric sites by therapeutics agents has the potential to enhance allosteric ligand binding, which leads to improved potency of allosteric drugs. Epidermal growth factor receptor (EGFR), as one of the most critical anti-cancer targets belonging to the receptor tyrosine kinase family, represents a quintessential example. It was revealed that osimertinib, an ATP-competitive covalent EGFR inhibitor, remarkably enhanced the affinity of a recently developed allosteric inhibitor JBJ-04-125-02 for EGFRL858R/T790M. Here, we utilized extensive large-scale molecular dynamics simulations and the reversed allosteric communication to untangle the detailed molecular underpinning, in which occupation of osimertinib at the orthosteric site altered the overall conformational ensemble of EGFR mutant and reshaped the allosteric site via long-distance signaling. A unique intermediate state resembling the active conformation was identified, which was further stabilized by osimertinib loading. Based on the allosteric communication pathway, we predicted a novel allosteric site positioned around K867, E868, H893, and K960 within the intermediate state. Its correlation with the orthosteric site was validated by both structural and energetic analysis, and its low sequence conservation indicated the potential for selective targeting across the human kinome. Together, these findings not only provided a mechanistic basis for future clinical application of the dual-targeting therapeutics, but also explored an innovative perception of allosteric inhibition of tyrosine kinase signaling.

D3DistalMutation: a Database to Explore the Effect of Distal Mutations on Enzyme Activity

Enzyme activity is affected by amino acid mutations, particularly mutations near the active site. Increasing evidence has shown that distal mutations more than 10 Å away from the active site may significantly affect enzyme activity. However, it is difficult to study the enzyme regulation mechanism of distal mutations due to the lack of a systematic collection of three-dimensional (3D) structures, highlighting distal mutation site and the corresponding enzyme activity change. Therefore, we constructed a distal mutation database, namely, D3DistalMutation, which relates the distal mutation to enzyme activity. As a result, we observed that approximately 80% of distal mutations could affect enzyme activity and 72.7% of distal mutations would decrease or abolish enzyme activity in D3DistalMutation. Only 6.6% of distal mutations in D3DistalMutation could increase enzyme activity, which have great potential to the industrial field. Among these mutations, the Y to F, S to D, and T to D mutations are most likely to increase enzyme activity, which sheds some light on industrial catalysis. Distal mutations decreasing enzyme activity in the allosteric pocket play an indispensable role in allosteric drug design. In addition, the pockets in the enzyme structures are provided to explore the enzyme regulation mechanism of distal mutations. D3DistalMutation is accessible free of charge at https://www.d3pharma.com/D3DistalMutation/index.php.

Ligand-based approach for predicting drug targets and for virtual screening against COVID-19

Discovering efficient drugs and identifying target proteins are still an unmet but urgent need for curing coronavirus disease 2019 (COVID-19). Protein structure-based docking is a widely applied approach for discovering active compounds against drug targets and for predicting potential targets of active compounds. However, this approach has its inherent deficiency caused by e.g. various different conformations with largely varied binding pockets adopted by proteins, or the lack of true target proteins in the database. This deficiency may result in false negative results. As a complementary approach to the protein structure-based platform for COVID-19, termed as D3Docking in our previous work, we developed in this study a ligand-based method, named D3Similarity, which is based on the molecular similarity evaluation between the submitted molecule(s) and those in an active compound database. The database is constituted by all the reported bioactive molecules against the coronaviruses, viz., severe acute respiratory syndrome coronavirus (SARS), Middle East respiratory syndrome coronavirus (MERS), severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), human betacoronavirus 2c EMC/2012 (HCoV-EMC), human CoV 229E (HCoV-229E) and feline infectious peritonitis virus (FIPV), some of which have target or mechanism information but some do not. Based on the two-dimensional (2D) and three-dimensional (3D) similarity evaluation of molecular structures, virtual screening and target prediction could be performed according to similarity ranking results. With two examples, we demonstrated the reliability and efficiency of D3Similarity by using 2D × 3D value as score for drug discovery and target prediction against COVID-19. The database, which will be updated regularly, is available free of charge at https://www.d3pharma.com/D3Targets-2019-nCoV/D3Similarity/index.php.

Next-generation computational tools and resources for coronavirus research: From detection to vaccine discovery

The COVID-19 pandemic has affected 215 countries and territories around the world with 60,187,347 coronavirus cases and 17,125,719 currently infected patients confirmed as of the November 25, 2020. Currently, many countries are working on developing new vaccines and therapeutic drugs for this novel virus strain, and a few of them are in different phases of clinical trials. The advancement in high-throughput sequence technologies, along with the application of bioinformatics, offers invaluable knowledge on genomic characterization and molecular pathogenesis of coronaviruses. Recent multi-disciplinary studies using bioinformatics methods like sequence-similarity, phylogenomic, and computational structural biology have provided an in-depth understanding of the molecular and biochemical basis of infection, atomic-level recognition of the viral-host receptor interaction, functional annotation of important viral proteins, and evolutionary divergence across different strains. Additionally, various modern immunoinformatic approaches are also being used to target the most promiscuous antigenic epitopes from the SARS-CoV-2 proteome for accelerating the vaccine development process. In this review, we summarize various important computational tools and databases available for systematic sequence-structural study on coronaviruses. The features of these public resources have been comprehensively discussed, which may help experimental biologists with predictive insights useful for ongoing research efforts to find therapeutics against the infectious COVID-19 disease.

Dynamical Correlations Reveal Allosteric Sites in G Protein-Coupled Receptors

G protein-coupled Receptors (GPCRs) play a central role in many physiological processes and, consequently, constitute important drug targets. In particular, the search for allosteric drugs has recently drawn attention, since they could be more selective and lead to fewer side effects. Accordingly, computational tools have been used to estimate the druggability of allosteric sites in these receptors. In spite of many successful results, the problem is still challenging, particularly the prediction of hydrophobic sites in the interface between the protein and the membrane. In this work, we propose a complementary approach, based on dynamical correlations. Our basic hypothesis was that allosteric sites are strongly coupled to regions of the receptor that undergo important conformational changes upon activation. Therefore, using ensembles of experimental structures, normal mode analysis and molecular dynamics simulations we calculated correlations between internal fluctuations of different sites and a collective variable describing the activation state of the receptor. Then, we ranked the sites based on the strength of their coupling to the collective dynamics. In the β2 adrenergic (β2AR), glucagon (GCGR) and M2 muscarinic receptors, this procedure allowed us to correctly identify known allosteric sites, suggesting it has predictive value. Our results indicate that this dynamics-based approach can be a complementary tool to the existing toolbox to characterize allosteric sites in GPCRs.

Computational Insights into the Conformational Accessibility and Binding Strength of SARS-CoV-2 Spike Protein to Human Angiotensin-Converting Enzyme 2

The spike protein of SARS-CoV-2 (CoV-2-S) mediates the virus entry into human cells. Experimental studies have shown the stronger binding affinity of the RBD (receptor binding domain) of CoV-2-S to angiotensin-converting enzyme 2 (ACE2) as compared to that of SARS-CoV spike (CoV-S). However, a similar or weaker binding affinity of CoV-2-S compared to that of CoV-S is observed if entire spikes are used in the bioassay. To explore the underlying mechanism, we calculated the binding affinities of the RBDs to ACE2 and simulated the transitions between ACE2-inaccessible and -accessible conformations. We found that the ACE2-accessible angle of CoV-2-S is 52.2° and that the ACE2 binding strength of CoV-2-S RBD is much stronger than that of CoV-S RBD. However, CoV-2-S has much less of an ACE2-accessible conformation and is much more difficult to shift from ACE2-inaccessible to -accessible than CoV-S, making the binding affinity of the entire protein decrease. Further analysis revealed key interactional residues for strong binding and five potential ligand-binding pockets for drug research.

Discovery of cryptic allosteric sites using reversed allosteric communication by a combined computational and experimental strategy

Allostery, which is one of the most direct and efficient methods to fine-tune protein functions, has gained increasing recognition in drug discovery. However, there are several challenges associated with the identification of allosteric sites, which is the fundamental cornerstone of drug design. Previous studies on allosteric site predictions have focused on communication signals propagating from the allosteric sites to the orthosteric sites. However, recent biochemical studies have revealed that allosteric coupling is bidirectional and that orthosteric perturbations can modulate allosteric sites through reversed allosteric communication. Here, we proposed a new framework for the prediction of allosteric sites based on reversed allosteric communication using a combination of computational and experimental strategies (molecular dynamics simulations, Markov state models, and site-directed mutagenesis). The desirable performance of our approach was demonstrated by predicting the known allosteric site of the small molecule MDL-801 in nicotinamide dinucleotide (NAD⁺)-dependent protein lysine deacetylase sirtuin 6 (Sirt6). A potential novel cryptic allosteric site located around the L116, R119, and S120 residues within the dynamic ensemble of Sirt6 was identified. The allosteric effect of the predicted site was further quantified and validated using both computational and experimental approaches. This study proposed a state-of-the-art computational pipeline for detecting allosteric sites based on reversed allosteric communication. This method enabled the identification of a previously uncharacterized potential cryptic allosteric site on Sirt6, which provides a starting point for allosteric drug design that can aid the identification of candidate pockets in other therapeutic targets.

Using reversed allosteric communication, we performed MD simulations, MSMs, and mutagenesis experiments, to discover allosteric sites. It reproduced the known allosteric site for MDL-801 on Sirt6 and uncovered a novel cryptic allosteric Pocket X.

SARS-CoV-2 transcriptome analysis and molecular cataloguing of immunodominant epitopes for multi-epitope based vaccine design

Genomics-led researches are engaged in tracing virus expression pattern, and induced immune responses in human to develop effective vaccine against COVID-19. In this study, targeted expression profiling and differential gene expression analysis of major histocompatibility complexes and innate immune system genes were performed through SARS-CoV-2 infected RNA-seq data of human cell line, and virus transcriptome was generated for T-and B-cell epitope prediction. Docking studies of epitopes with MHC and B-cell receptors were performed to identify potential T-and B-cell epitopes. Transcriptome analysis revealed the specific multiple allele expressions in cell line, genes for elicited induce immune response, and virus gene expression. Proposed T- and B-cell epitopes have high potential to elicit equivalent immune responses caused by SARS-CoV-2 infection which can be useful to provide links between elicited immune response and virus gene expression. This study will facilitate in vitro and in vivo vaccine related research studies in disease control.

•

SARS-CoV-2 transcriptome construction from RNA-seq data of infected human cell-lines.

•

T- and B-cell epitopes from SARS-CoV-2 transcriptome

•

Gene expression profiling of MHC alleles and innate immune system genes

Study on effects of co-solvents on the structure of DhaA by molecular dynamics simulation

With the increasing application of enzymes in various research fields, the choices of co-solvents in enzymatic preparations which directly related to the catalytic activity have been attracted attention. Thus, researching on the stabilization or destabilization behaviors of enzymes in different solvents is extremely essential. In this study, the structural changes of DhaA in two typical aprotic co-solvents (acetonitrile and tetrahydrofuran) were firstly investigated by molecular dynamics (MD) simulation. The simulation results revealed the strong van der Waals force between co-solvents and DhaA which could induce the structural change of enzyme. Interestingly, the differences of molecular size and the electrostatic force with enzyme of two co-solvents led to quite different influences on DhaA. As for acetonitrile, solvent molecules could penetrate into the catalytic site of DhaA which promoted by the electrostatic interaction. On the contrary, tetrahydrofuran molecules were mainly distributed around the catalytic site due to the relative weak electrostatic interaction and steric resistance effect. It can be concluded that different co-solvent can affect the key domains, substrate pathway and catalytic pocket of DhaA.Communicated by Ramaswamy H. Sarma.

D3Targets-2019-nCoV: a webserver for predicting drug targets and for multi-target and multi-site based virtual screening against COVID-19

A highly effective medicine is urgently required to cure coronavirus disease 2019 (COVID-19). For the purpose, we developed a molecular docking based webserver, namely D3Targets-2019-nCoV, with two functions, one is for predicting drug targets for drugs or active compounds observed from clinic or in vitro/in vivo studies, the other is for identifying lead compounds against potential drug targets via docking. This server has its unique features, (1) the potential target proteins and their different conformations involving in the whole process from virus infection to replication and release were included as many as possible; (2) all the potential ligand-binding sites with volume larger than 200 Å3 on a protein structure were identified for docking; (3) correlation information among some conformations or binding sites was annotated; (4) it is easy to be updated, and is accessible freely to public (https://www.d3pharma.com/D3Targets-2019-nCoV/index.php). Currently, the webserver contains 42 proteins [20 severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) encoded proteins and 22 human proteins involved in virus infection, replication and release] with 69 different conformations/structures and 557 potential ligand-binding pockets in total. With 6 examples, we demonstrated that the webserver should be useful to medicinal chemists, pharmacologists and clinicians for efficiently discovering or developing effective drugs against the SARS-CoV-2 to cure COVID-19.

A systematic strategy for the investigation of vaccines and drugs targeting bacteria

Infectious and epidemic diseases induced by bacteria have historically caused great distress to people, and have even resulted in a large number of deaths worldwide. At present, many researchers are working on the discovery of viable drug and vaccine targets for bacteria through multiple methods, including the analyses of comparative subtractive genome, core genome, replication-related proteins, transcriptomics and riboswitches, which plays a significant part in the treatment of infectious and pandemic diseases. The 3D structures of the desired target proteins, drugs and epitopes can be predicted and modeled through target analysis. Meanwhile, molecular dynamics (MD) analysis of the constructed drug/epitope-protein complexes is an important standard for testing the suitability of these screened drugs and vaccines. Currently, target discovery, target analysis and MD analysis are integrated into a systematic set of drug and vaccine analysis strategy for bacteria. We hope that this comprehensive strategy will help in the design of high-performance vaccines and drugs.