Computational prediction of B and T-cell epitopes of Kyasanur Forest Disease virus marker proteins towards the development of precise diagnosis and potent subunit vaccine

Kyasanur Forest Disease (KFD), also known as 'monkey fever', caused by KFD Virus (KFDV), is a highly neglected tropical disease endemic to Western Ghat region of Karnataka, India. Recently, KFD, which is fatal for both monkeys and humans with a mortality rate of 2-10% has been found to spread from its epicenter to neighboring districts and states also. The current ELISA based KFD detection method is very non-specific due to cross-reactivity with other flaviviruses. Further, presently available formalin-inactivated vaccine has been found to be less effective leading to disease susceptibility and severity. To address these, the present study was aimed at predicting the potent specific B and T-cell epitopes of KFDV immunogenic marker proteins using diverse computational tools aiming at developing precise diagnostic method and an effective subunit vaccine. Here, we have chosen E, NS1 and NS5 proteins as markers of KFDV by taking into account of their differential and non-overlapping sequences with selected arboviruses. Based on the linear and nonlinear epitope prediction tools and important biophysical parameters, we identified three potential linear and ten nonlinear B-cell epitopes. We also predicted T-cell epitope peptides which binds to MHC class-I and class-II receptors for the effective T-cell activation. Thus, our molecular docking and molecular dynamics simulation analysis has identified six different TH-cell epitopes based on the distribution frequency of MHC-II haplotypes in the human population and one TC-cell epitope from NS5 protein that has maximum interaction with class-I MHC. Overall, we have successfully identified potential B and T-cell epitope marker peptides present in the envelope and two non-structural proteins.Communicated by Ramaswamy H. Sarma.

Exploring the potential and identifying Withania somnifera alkaloids as novel dihydrofolate reductase (DHFR) inhibitors by the AlteQ method

There is an urgent need to discover and develop novel drugs to combat Mycobacterium tuberculosis, the causative agent of tuberculosis (TB) in humans. Alkaloids have been shown to have wide-ranging therapeutic application and could be ideal candidates for drug development, and research is underway to develop new anti-tubercular drugs from natural sources. In this regard, the current research deals with finding novel lead compounds from the Withania somnifera (WS) plant. Broad health benefits of WS are due to the presence of diverse chemical constituents which include anaferine and anahygrine and which belong to the alkaloid family. In the present study, these two compounds have been theoretically studied to understand their electronic properties using the density functional theory (DFT) at the B3LYP/6-311 + G (d,p) level. HOMO and LUMO properties and molecular electrostatic potential (MEP) surface were calculated. Further, to understand the mechanism of action of these compounds and to identify their putative drug target, molecular docking and dynamics studies were employed against Mycobacterium tuberculosis dihydrofolate reductase (DHFR). It was determined that NADP+ affects stability of the complexes by reducing fluctuations of residues 14-23 and 117-126. It was also found that Ile5 and Gln28 play an important role in complexation. Electron density analysis (using the AlteQ method) of the intermolecular region, analyzing both the anaferin-NADP+ and anahygrin-NADP+ complexes showed that anaferin and anahygrin complexes are more stable in the presence of NADP+. It has been established that in most intermolecular contacts the contribution of the ligand to the electron density is greater than that of NADP+. The present study thus provides an excellent way to analyze the effect of anaferine and anahygrine in essential processes of M. tuberculosis.Communicated by Ramaswamy H. Sarma.

Proposition of a new allosteric binding site for potential SARS-CoV-2 3CL protease inhibitors by utilizing molecular dynamics simulations and ensemble docking

The SARS-CoV-2 3CL protease (3CLpro) shows a high similarity with 3CL proteases of other beta-coronaviruses, such as SARS and MERS. It is the main enzyme involved in generating various non-structural proteins that are important for viral replication and is one of the most important proteins responsible for SARS-CoV-2 virulence. In this study, we have conducted an ensemble docking of molecules from the DrugBank database using both the crystallographic structure of the SARS-CoV-2 3CLpro, as well as five conformations obtained after performing a cluster analysis of a 300 ns molecular dynamics (MD) simulation. This procedure elucidated the inappropriateness of the active site for non-covalent inhibitors, but it has also shown that there exists an additional, more favorable, allosteric binding site, which could be a better target for non-covalent inhibitors, as it could prevent dimerization and activation of SARS-CoV-2 3CLpro. Two such examples are radotinib and nilotinib, tyrosine kinase inhibitors already in use for treatment of leukemia and which binding to the newly found allosteric binding site was also confirmed using MD simulations. Communicated by Ramaswamy H. Sarma.

Hybrid Quinazoline 1,3,5-Triazines as Epidermal Growth Factor Receptor (EGFR) Inhibitors with Anticancer Activity: Design, Synthesis, and Computational Study

We report a series of hybrid quinazoline-1,3,5-triazine derivatives as EGFR inhibitors, which were synthesised and tested by using a variety of in vitro, in silico, and in vivo techniques. The derivatives were found to be active against different cancer cell lines and nontoxic against normal ones, with compounds 7 c, 7 d, 7 e, and 7 j being the most potent ones. The derivatives were also evaluated for angiogenesis inhibition potency in chicken eggs, and molecular docking and dynamics simulation studies were carried out to elucidate the fundamental substituent groups essential for their bioactivity. Additionally, a SAR study of the derivatives was performed for future compound optimisation. These studies suggested that the derivatives have a high affinity towards EGFR with favourable pharmacological properties. The most active compound (7 e) was further evaluated for in vivo anticancer activity against DMBA-induced tumours in female Sprague-Dawley rats as well as its effects on plasma antioxidant status, biotransformation enzymes, and lipid profile. The study suggested that 7 e has lead properties against breast cancer and can serve as a starting compound for further development of anti-EGFR compounds.

Pathway enrichment analysis of virus-host interactome and prioritization of novel compounds targeting the spike glycoprotein receptor binding domain-human angiotensin-converting enzyme 2 interface to combat SARS-CoV-2

SARS-CoV-2 has become a pandemic causing a serious global health concern. The absence of effective drugs for treatment of the disease has caused its rapid spread on a global scale. Similarly to the SARS-CoV, the SARS-CoV-2 is also involved in a complex interplay with the host cells. This infection is characterized by a diffused alveolar damage consistent with the Acute Respiratory Disease Syndrome (ARDS). To explore the complex mechanisms of the disease at the system level, we used a network medicine tools approach. The protein-protein interactions (PPIs) between the SARS-CoV and the associated human cell proteins are crucial for the viral pathogenesis. Since the cellular entry of SARS-CoV-2 is accomplished by binding of the spike glycoprotein binding domain (RBD) to the human angiotensin-converting enzyme 2 (hACE2), a molecule that can bind to the spike RDB-hACE2 interface could block the virus entry. Here, we performed a virtual screening of 55 compounds to identify potential molecules that can bind to the spike glycoprotein and spike-ACE2 complex interface. It was found that the compound ethyl 1-{3-[(2,4-dichlorobenzyl) carbamoyl]-1-ethyl-6-fluoro-4-oxo-1,4-dihydro-7-quinolinyl}-4-piperidine carboxylate (the S54 ligand) and ethyl 1-{3-[(2,4-dichlorobenzyl) carbamoyl]-1-ethyl-6-fluoro-4-oxo-1,4-dihydro-7-quinolinyl}-4 piperazine carboxylate (the S55 ligand) forms hydrophobic interactions with Tyr41A, Tyr505B and Tyr553B, Leu29A, Phe495B, respectively of the spike glycoprotein, the hotspot residues in the spike glycoprotein RBD-hACE2 binding interface. Furthermore, molecular dynamics simulations and free energy calculations using the MM-GBSA method showed that the S54 ligand is a stronger binder than a known SARS-CoV spike inhibitor SSAA09E3 (N-(9,10-dioxo-9, 10-dihydroanthracen-2-yl) benzamide).

Communicated by Ramaswamy H. Sarma

What a Difference a Methyl Group Makes: The Selectivity of Monoamine Oxidase B Towards Histamine and N-Methylhistamine

Monoamine oxidase (MAO) enzymes catalyze the degradation of a very broad range of biogenic and dietary amines including many neurotransmitters in the brain, whose imbalance is extensively linked with the biochemical pathology of various neurological disorders. Although sharing around 70 % sequence identity, both MAO A and B isoforms differ in substrate affinities and inhibitor sensitivities. Inhibitors that act on MAO A are used to treat depression, due to their ability to raise serotonin concentrations, whereas MAO B inhibitors decrease dopamine degradation and improve motor control in patients with Parkinson disease. Despite this functional importance, the factors affecting MAO selectivity are poorly understood. Here, we used a combination of molecular dynamics (MD) simulations, molecular mechanics with Poisson-Boltzmann and surface area solvation (MM-PBSA) binding free energy evaluations, and quantum mechanical (QM) cluster calculations to address the unexpected, yet challenging MAO B selectivity for N-methylhistamine (NMH) over histamine (HIS), differing only in a single methyl group distant from the reactive ethylamino center. This study shows that a dominant selectivity contribution is offered by a lower activation free energy for NMH by 2.6 kcal mol^-1 , in excellent agreement with the experimental ΔΔG^≠_EXP =1.4 kcal mol^-1 , together with a more favorable reaction exergonicity and active-site binding. This study also confirms the hydrophobic nature of the MAO B active site and underlines the important role of Ile199, Leu171, and Leu328 in properly orienting substrates for the reaction.