Exploring the immune evasion of SARS-CoV-2 variant harboring E484K by molecular dynamics simulations.

Although the current coronavirus disease 2019 (COVID-19) vaccines have been used worldwide to halt spread of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the emergence of new SARS-CoV-2 variants with E484K mutation shows significant resistance to the neutralization of vaccine sera. To better understand the resistant mechanism, we calculated the binding affinities of 26 antibodies to wild-type (WT) spike protein and to the protein harboring E484K mutation, respectively. The results showed that most antibodies (~85%) have weaker binding affinities to the E484K mutated spike protein than to the WT, indicating the high risk of immune evasion of the mutated virus from most of current antibodies. Binding free energy decomposition revealed that the residue E484 forms attraction with most antibodies, while the K484 has repulsion from most antibodies, which should be the main reason of the weaker binding affinities of E484K mutant to most antibodies. Impressively, a monoclonal antibody (mAb) combination was found to have much stronger binding affinity with E484K mutant than WT, which may work well against the mutated virus. Based on binding free energy decomposition, we predicted that the mutation of four more residues on receptor-binding domain (RBD) of spike protein, viz., F490, V483, G485 and S494, may have high risk of immune evasion, which we should pay close attention on during the development of new mAb therapeutics.

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.

D3Rings: A Fast and Accurate Method for Ring System Identification and Deep Generation of Drug-Like Cyclic Compounds.

Continuous exploration of the chemical space of molecules to find ligands with high affinity and specificity for specific targets is an important topic in drug discovery. A focus on cyclic compounds, particularly natural compounds with diverse scaffolds, provides important insights into novel molecular structures for drug design. However, the complexity of their ring structures has hindered the applicability of widely accepted methods and software for the systematic identification and classification of cyclic compounds. Herein, we successfully developed a new method, D3Rings, to identify acyclic, monocyclic, spiro ring, fused and bridged ring, and cage ring compounds, as well as macrocyclic compounds. By using D3Rings, we completed the statistics of cyclic compounds in three different databases, e.g., ChEMBL, DrugBank, and COCONUT. The results demonstrated the richness of ring structures in natural products, especially spiro, macrocycles, and fused and bridged rings. Based on this, three deep generative models, namely, VAE, AAE, and CharRNN, were trained and used to construct two data sets similar to DrugBank and COCONUT but 10 times larger than them. The enlarged data sets were then used to explore the molecular chemical space, focusing on complex ring structures, for novel drug discovery and development. Docking experiments with the newly generated COCONUT-like data set against three SARS-CoV-2 target proteins revealed that an expanded compound database improves molecular docking results. Cyclic structures exhibited the best docking scores among the top-ranked docking molecules. These results suggest the importance of exploring the chemical space of structurally novel cyclic compounds and continuous expansion of the library of drug-like compounds to facilitate the discovery of potent ligands with high binding affinity to specific targets. D3Rings is now freely available at http://www.d3pharma.com/D3Rings/.

TRIP13 - a potential drug target in cancer pharmacotherapy.

ATPases Associated with Diverse Cellular Activity (AAA+ATPases) are important enzymatic functional proteins in human cells. Thyroid Hormone Receptor Interacting Protein-13 (TRIP13) is a member of this protein superfamily, that partly regulates DNA repair pathways and spindle assembly checkpoints during mitosis. TRIP13 is reported as an oncogene involving multiple pathways in many human malignancies, including multiple myeloma, brain tumors, etc. The structure of TRIP13 reveals the mechanisms for ATP binding and how TRIP13 recognizes the Mitotic Arrest Deficiency-2 (MAD2) protein, with p31comet acting as an adapter protein. DCZ0415, TI17, DCZ5417, and DCZ5418 are the reported small-molecule inhibitors of TRIP13, which have been demonstrated to inhibit TRIP13's biological functions significantly and effective in suppressing various types of malignant cells, indicating that TRIP13 is a significant anticancer drug target. Currently, no systematic reviews are cutting across the functions, structure, and novel inhibitors of TRIP13. This review provides a comprehensive overview of TRIP13's biological functions, its roles in eighteen different cancers, four small molecule inhibitors, different underlying molecular mechanisms, and its functionality as a potential anticancer drug target.

Preferred microenvironments of halogen bonds and hydrogen bonds revealed using statistics and QM/MM calculation studies.

Hydrogen bonds (HBs) and halogen bonds (XBs) are two essential non-covalent interactions for molecular recognition and drug design. As proteins are heterogeneous in structure, the microenvironments of protein structures should have effects on the formation of HBs and XBs with ligands. However, there are no systematic studies reported on this effect to date. For quantitatively describing protein microenvironments, we defined the local hydrophobicities (LHs) and local dielectric constants (LDCs) in this study. With the defined parameters, we conducted an elaborate database survey on the basis of 22 011 ligand-protein structures to explore the microenvironmental preference of HBs (91 966 in total) and XBs (1436 in total). The statistics show that XBs prefer hydrophobic microenvironments compared to HBs. The polar residues like ASP are more likely to form HBs with ligands, while nonpolar residues such as PHE and MET prefer XBs. Both the LHs and LDCs (10.69 ± 4.36 for HBs; 8.86 ± 4.00 for XBs) demonstrate that XBs are prone to hydrophobic microenvironments compared with HBs with significant differences (p < 0.001), indicating that evaluating their strengths in the corresponding environments should be necessary. Quantum Mechanics-Molecular Mechanics (QM/MM) calculations reveal that in comparison with vacuum environments, the interaction energies of HBs and XBs are decreased to varying degrees given different microenvironments. In addition, the strengths of HBs are impaired more than those of XBs when the local dielectric constant's difference between the XB microenvironments and the HB microenvironments is large.

Exploring the Binding Affinity and Mechanism between ACE2 and the Trimers of Delta and Omicron Spike Proteins by Molecular Dynamics Simulation and Bioassay.

Five major variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have emerged and posed challenges in controlling the pandemic. Among them, the current dominant variant, viz., Omicron, has raised serious concerns about its infectiousness and antibody neutralization. However, few studies pay attention to the effect of the mutations on the dynamic interaction network of Omicron S protein trimers binding to the host angiotensin-converting enzyme 2 (ACE2). In this study, we conducted molecular dynamics (MD) simulations and enzyme linked immunosorbent assay (ELISA) to explore the binding strength and mechanism of wild type (WT), Delta, and Omicron S protein trimers to ACE2. The results showed that the binding capacities of both the two variants' S protein trimers to ACE2 are enhanced in varying degrees, indicating possibly higher cell infectiousness. Energy decomposition and protein-protein interaction network analysis suggested that both the mutational and conserved sites make effects on the increase in the overall affinity through a variety of interactions. The experimentally determined KD values by biolayer interferometry (BLI) and the predicted binding free energies of the RBDs of Delta and Omicron to mAb HLX70 revealed that the two variants may have the high risk of immune evasion from the mAb. These results are not only helpful in understanding the binding strength and mechanism of S protein trimer-ACE2 but also beneficial for drug, especially for antibody development.

Predicting spike protein NTD mutations of SARS-CoV-2 causing immune evasion by molecular dynamics simulations.

The coronavirus disease 2019 (COVID-19) pandemic was caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Among all the potential targets studied for developing drugs and antibodies, the spike (S) protein is the most striking one, which is on the surface of the virus. In contrast with the intensively investigated immunodominant receptor-binding domain (RBD) of the protein, little is known about the neutralizing antibody binding mechanisms of the N-terminal domain (NTD), let alone the effects of NTD mutations on antibody binding and thereby the risk of immune evasion. Based on 400 ns molecular dynamics simulation for 11 NTD-antibody complexes together with other computational approaches in this study, we investigated critical residues for NTD-antibody binding and their detailed mechanisms. The results show that 36 residues on the NTD including R246, Y144, K147, Y248, L249 and P251 are critically involved in the direct interaction of the NTD with many monoclonal antibodies (mAbs), indicating that the viruses harboring these residue mutations may have a high risk of immune evasion. Binding free energy calculations and an interaction mechanism study reveal that R246I, which is present in the Beta (B.1.351/501Y.V2) variant, may have various impacts on current NTD antibodies through abolishing the hydrogen bonds and electrostatic interaction with the antibodies or affecting other interface residues. Therefore, special attention should be paid to the mutations of these key residues in future antibody and vaccine design and development.

Impact of Halogen Bonds on Protein-Peptide Binding and Protein Structural Stability Revealed by Computational Approaches.

Halogen bonds (XBs) are essential noncovalent interactions in molecular recognition and drug design. Current studies on XBs in drug design mainly focus on the interactions between halogenated ligands and target proteins, lacking a systematic study of naturally existing and artificially prepared halogenated residue XBs (hr_XBs) and their characteristics. Here, we conducted a computational study on the potential hr_XBs in proteins/peptides using database searching, quantum mechanics calculations, and molecular dynamics simulations. XBs at the protein-peptide interaction interfaces are found to enhance their binding affinity. Additionally, the formation of intramolecular XBs (intra_XBs) within proteins may significantly contribute to the structural stability of structurally flexible proteins while having a minor impact on proteins with inherently high structural rigidity. Impressively, introducing halogens without the formation of intra_XBs may lead to a decrease in the protein structural stability. This study enriches our understanding of the roles and effects of halogenated residue XBs in biological systems.

Unexpected effect of halogenation on the water solubility of small organic compounds.

Halogenation is an indispensable method in the structural modification of lead compounds. It is known to increase lipophilicity and is hence used to improve membrane permeability and thus bioavailability. In this study, we compare the water solubility (logS) of organohalogen compounds and their non-halogenated parent compounds using the molecular matched pair (MMP) analysis method. Unexpectedly, 19.9% of the compounds increased their water solubility upon halogenation. Iodination was observed to have the greatest effect on solubility, followed by chlorination, bromination, and fluorination. Introducing amino, hydroxyl and carboxyl groups into organohalogens improves their aqueous solubilities, whereas introducing a trifluoromethyl group has the opposite effect. According to our quantum chemical calculations, the increased water solubility upon halogenation is, at least partially, attributed to an increased polarity and polarizability. These results improve our understanding of the influence of halogenation on bioactivity.

Global dynamics and pattern formation for predator-prey system with density-dependent motion.

In this paper, we concern with the predator-prey system with generalist predator and density-dependent prey-taxis in two-dimensional bounded domains. We derive the existence of classical solutions with uniform-in-time bound and global stability for steady states under suitable conditions through the Lyapunov functionals. In addition, by linear instability analysis and numerical simulations, we conclude that the prey density-dependent motility function can trigger the periodic pattern formation when it is monotone increasing.

D3CARP: a comprehensive platform with multiple-conformation based docking, ligand similarity search and deep learning approaches for target prediction and virtual screening.

Resource- and time-consuming biological experiments are unavoidable in traditional drug discovery, which have directly driven the evolution of various computational algorithms and tools for drug-target interaction (DTI) prediction. For improving the prediction reliability, a comprehensive platform is highly expected as some previously reported webservers are small in scale, single-method, or even out of service. In this study, we integrated the multiple-conformation based docking, 2D/3D ligand similarity search and deep learning approaches to construct a comprehensive webserver, namely D3CARP, for target prediction and virtual screening. Specifically, 9352 conformations with positive control of 1970 targets were used for molecular docking, and approximately 2 million target-ligand pairs were used for 2D/3D ligand similarity search and deep learning. Besides, the positive compounds were added as references, and related diseases of therapeutic targets were annotated for further disease-based DTI study. The accuracies of the molecular docking and deep learning approaches were 0.44 and 0.89, respectively. And the average accuracy of five ligand similarity searches was 0.94. The strengths of D3CARP encompass the support for multiple computational methods, ensemble docking, utilization of positive controls as references, cross-validation of predicted outcomes, diverse disease types, and broad applicability in drug discovery. The D3CARP is freely accessible at https://www.d3pharma.com/D3CARP/index.php.

Molecular dynamics study on the behavior and binding mechanism of target protein Transgelin-2 with its agonist TSG12 for anti-asthma drug discovery.

Transgelin-2 (TG2) is a novel promising therapeutic target for the treatment of asthma as it plays an important role in relaxing airway smooth muscles and reducing pulmonary resistance in asthma. The compound TSG12 is the only reported TG2 agonist with in vivo anti-asthma activity. However, the dynamic behavior and ligand binding sites of TG2 and its binding mechanism with TSG12 remain unclear. In this study, we performed 12.6 µs molecular dynamics (MD) simulations for apo-TG2 and TG2-TSG12 complex, respectively. The results suggested that the apo-TG2 has 4 most populated conformations, and that its binding of the agonist could expand the conformation distribution space of the protein. The simulations revealed 3 potential binding sites in 3 most populated conformations, one of which is induced by the agonist binding. Free energy decomposition uncovered 8 important residues with contributions stronger than -1 kcal/mol. Computational alanine scanning for the important residues by 100 ns conventional MD simulation for each mutated TG2-TSG12 complexes demonstrated that E27, R49 and F52 are essential residues for the agonist binding. These results should be helpful to understand the dynamic behavior of TG2 and its binding mechanism with the agonist TSG12, which could provide some structural insights into the novel mechanism for anti-asthma drug development.

Discovery of zolinium TSG1180 as a novel agonist of transgelin-2 for treating asthma.

Asthma is a complex and heterogeneous respiratory disease that causes serious social and economic burdens. Current drugs such as ß2-agonists cannot fully control asthma. Our previous study found that Transgelin-2 is a potential target for treating asthmatic pulmonary resistance. Herein, we discovered a zolinium compound, TSG1180, that showed a strong interaction with Transgelin-2. The equilibrium dissociation constants (KD) of TSG1180 to Transgelin-2 were determined to be 5.363 × 10-6 and 9.81 × 10-6 M by surface plasmon resonance (SPR) and isothermal titration calorimetry (ITC). Cellular thermal shift assay (CETSA) results showed that the thermal stability of Transgelin-2 increased after coincubation of TSG1180 with lysates of airway smooth muscle cells (ASMCs). Molecular docking showed that Arg39 may be the key residue for the binding. Then, the SPR result showed that the binding affinity of TSG1180 to Transgelin-2 mutant (R39E) was decreased by 1.69-fold. Real time cell analysis (RTCA) showed that TSG1180 treatment could relax ASMCs by 19 % (P < 0.05). Once Transgelin-2 was inhibited, TSG1180 cannot induce a relaxation effect, suggesting that the relaxation effect was specifically mediated by Transgelin-2. In vivo study showed TSG1180 effectively reduced pulmonary resistance by 64 % in methacholine-induced mice model (P < 0.05). Furthermore, the phosphorylation of Ezrin at T567 was increased by 8.06-fold, the phosphorylation of ROCK at Y722 was reduced by 38 % and the phosphorylation of RhoA at S188 was increased by 52 % after TSG1180 treatment. These results suggested that TSG1180 could be a Transgelin-2 agonist for further optimization and development as an anti-asthma drug.

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.