Potent protease inhibitors of deadly lagoviruses: rabbit hemorrhagic disease virus and European brown hare syndrome virus

Rabbit hemorrhagic disease (RHD) and European brown hare syndrome (EBHS) are highly contagious diseases caused by lagoviruses in the Caliciviridae family and mainly affect rabbits and hares, respectively. These infectious diseases are associated with high mortality and a serious threat to domesticated (farmed and pet) and wild rabbits and hares, including endangered species such as Riparian brush rabbits. In the US, only isolated cases of RHD had been reported until Spring 2020. However, RHD caused by RHD type 2 virus (RHDV2) was unexpectedly reported in April 2020 in New Mexico and has subsequently spread to several US states infecting wild rabbits and hares. Since it is almost impossible to control and eradicate the virus from wild animals, it is highly likely RHD will become endemic in the US. Vaccines are available for RHD, however, there is no specific treatment for these deadly diseases. RHDV and EBHSV encode a 3C-like protease (3CLpro), which is essential for virus replication and a promising target for antiviral drug development. We have previously generated focused small molecule libraries of 3CLpro inhibitors and demonstrated the in vitro potency and in vivo efficacy of some protease inhibitors against viruses that encode 3CLpro including caliciviruses and coronaviruses. Here we established the enzyme assay and cell-based assays for these uncultivable viruses to determine the in vitro activity of 3CLpro inhibitors, including GC376, a protease inhibitor being developed for feline infectious peritonitis, and identified potent inhibitors of RHDV1 and 2 and EBHSV. In addition, structure-activity relationship study and homology modelling of the 3CLpros and inhibitors revealed that lagoviruses share similar structural requirements for 3CLpro inhibition with other caliciviruses.

Correction to: In silico structural homology modelling of EST073 motif coding protein of tea Camellia sinensis (L)

An amendment to this paper has been published and can be accessed via the original article.

Effect of Divalent Metal Ion on the Structure, Stability and Function of Klebsiella pneumoniae Nicotinate-Nucleotide Adenylyltransferase: Empirical and Computational Studies

The continuous threat of drug-resistant Klebsiella pneumoniae justifies identifying novel targets and developing effective antibacterial agents. A potential target is nicotinate nucleotide adenylyltransferase (NNAT), an indispensable enzyme in the biosynthesis of the cell-dependent metabolite, NAD+. NNAT catalyses the adenylation of nicotinamide/nicotinate mononucleotide (NMN/NaMN), using ATP to form nicotinamide/nicotinate adenine dinucleotide (NAD+/NaAD). In addition, it employs divalent cations for co-substrate binding and catalysis and has a preference for different divalent cations. Here, the biophysical structure of NNAT from K. pneumoniae (KpNNAT) and the impact of divalent cations on its activity, conformational stability and substrate-binding are described using experimental and computational approaches. The experimental study was executed using an enzyme-coupled assay, far-UV circular dichroism, extrinsic fluorescence spectroscopy, and thermal shift assays, alongside homology modelling, molecular docking, and molecular dynamic simulation. The structure of KpNNAT revealed a predominately α-helical secondary structure content and a binding site that is partially hydrophobic. Its substrates ATP and NMN share the same binding pocket with similar affinity and exhibit an energetically favourable binding. KpNNAT showed maximum activity and minimal conformational changes with Mg2+ as a cofactor compared to Zn2+, Cu2+ and Ni2+. Overall, ATP binding affects KpNNAT dynamics, and the dynamics of ATP binding depend on the presence and type of divalent cation. The data obtained from this study would serve as a basis for further evaluation towards designing structure-based inhibitors with therapeutic potential.

Molecular Informatics of Trypanothione Reductase of Leishmania major Reveals Novel Chromen-2-One Analogues as Potential Leishmanicides

Trypanothione reductase (TR), a flavoprotein oxidoreductase is an important therapeutic target for leishmaniasis. Ligand-based pharmacophore modelling and molecular docking were used to predict selective inhibitors against TR. Homology modelling was employed to generate a three-dimensional structure of Leishmania major trypanothione reductase (LmTR). A pharmacophore model used to screen a natural compound library generated 42 hits, which were docked against the LmTR protein. Compounds with lower binding energies were evaluated via in silico pharmacological profiling and bioactivity. Four compounds emerged as potential leads comprising Karatavicinol (7-[(2E,6E,10S)-10,11-dihydroxy-3,7,11-trimethyldodeca-2,6-dienoxy]chromen-2-one), Marmin (7-[(E,6R)-6,7-dihydroxy-3,7-dimethyloct-2-enoxy]chromen-2-one), Colladonin (7-[[(4aS)-6-hydroxy-5,5,8a-trimethyl-2-methylidene-3,4,4a,6,7,8-hexahydro-1H-naphthalen-1-yl]methoxy]chromen-2-one), and Pectachol (7-[(6-hydroxy-5,5,8a-trimethyl-2-methylidene-3,4,4a,6,7,8-hexahydro-1H-naphthalen-1-yl)methoxy]-6,8-dimethoxychromen-2-one) with good binding energies of −9.4, −9.3, 8.8, and −8.5 kcal/mol, respectively. These compounds bound effectively to the FAD domain of the protein with some critical residues including Asp35, Thr51, Lys61, Tyr198, and Asp327. Furthermore, molecular dynamics simulations and molecular mechanics Poisson-Boltzmann surface area (MMPBSA) computations corroborated their strong binding. The compounds were also predicted to possess anti-leishmanial activity. The molecules serves as templates for the design of potential drug candidates and can be evaluated in vitro with optimistic results in producing plausible attenuating infectivity in macrophages.

Residues R1075, D1090, R1095, and C1130 Are Critical in ADAMTS13 TSP8-Spacer Interaction Predicted by Molecular Dynamics Simulation

ADAMTS13 (A Disintegrin and Metalloprotease with Thrombospondin type 1 repeats, member 13) cleaves von Willebrand Factor (VWF) multimers to limit the prothrombotic function of VWF. The deficiency of ADAMTS13 causes a lethal thrombotic microvascular disease, thrombotic thrombocytopenic purpura (TTP). ADAMTS13 circulates in a “closed” conformation with the distal domain associating the Spacer domain to avoid off-target proteolysis or recognition by auto-antibodies. However, the interactions of the distal TSP8 domain and the Spacer domain remain elusive. Here, we constructed the TSP8-Spacer complex by a combination of homology modelling and flexible docking. Molecular dynamics simulation was applied to map the binding sites on the TSP8 or Spacer domain. The results predicted that R1075, D1090, R1095, and C1130 on the TSP8 domain were key residues that interacted with the Spacer domain. R1075 and R1095 bound exosite-4 tightly, D1090 formed multiple hydrogen bonds and salt bridges with exosite-3, and C1130 interacted with both exosite-3 and exosite-4. Specific mutations of exosite-3 (R568K/F592Y/R660K/Y661F/Y665F) or the four key residues (R1075A/D1090A/R1095A/C1130A) impaired the binding of the TSP8 domain to the Spacer domain. These results shed new light on the understanding of the auto-inhibition of ADAMTS13.

Functional modulation of PTH1R activation and signalling by RAMP2

Receptor-activity-modifying proteins (RAMPs) are ubiquitously expressed membrane proteins that associate with different G protein-coupled receptors (GPCRs) including the parathyroid hormone 1 receptor (PTH1R), a class B GPCR, and an important modulator of mineral ion homeostasis and bone metabolism. However, it is unknown whether and how RAMP proteins may affect PTH1R function. Using different optical biosensors to measure the activation of PTH1R and its downstream signalling, we describe here that RAMP2 acts as a specific allosteric modulator of PTH1R, shifting PTH1R to a unique pre-activated state that permits faster activation in a ligand-specific manner. Moreover, RAMP2 modulates PTH1R downstream signalling in an agonist-dependent manner, most notably increasing the PTH-mediated Gi3 signalling sensitivity. Additionally, RAMP2 increases both PTH- and PTHrP-triggered β-arrestin2 recruitment to PTH1R. Employing homology modelling we describe the putative structural molecular basis underlying our functional findings. These data uncover a critical role of RAMPs in the activation and signalling of a GPCR that may provide a new venue for highly specific modulation of GPCR function and advanced drug design.

A Pragmatic Pharmacophore Informatics Strategy to Discover New Potent Inhibitors Against Pim-3

Pim-3 (proviral integration site moloney murine leukemia virus-3) is an oncogene which encodes proteins belonging to serine/threonine kinase family, and PIM subfamily. It is generally over expressed in epithelial and hematological tumors. It is known to involve in numerous cellular functions such as cell growth, differentiation, survival, tumorigenesis and apoptosis. It also plays a crucial role in regulation of signal transduction cascades. Therefore it emerged as a hopefultherapeutic target for cancer treatment. In current study, indole derivatives having potent inhibitory activity against Pim 3 were taken and pharmacophore based virtual screening was carried out. A five point pharmacophore hypothesis with one hydrogen bond acceptor, one hydrogen bond donor and three aromatic rings i.e., ADRRR was developed with acceptable R2and Q2 values of 0.913 and 0.748 respectively. It was employed as a query and screening was conducted against Asinex and Otava lead library databases to screen out potent drug like candidates. The obtained compounds were subjected to SP, XP docking using 3D model of pim-3 which was constructed through comparative homology modelling and finally binding free energies were calculated for top hits. The docking and binding free energy studies revealed that six hit molecules showed higher binding energy in comparison to the best active molecule. Finally, MD simulations of the top hit with highest binding energy was carried out which indicated that the obtained hit N1 formed a stable complex with pim-3. We believe that these combined protocols will be helpful and cooperative to discover and design more potent pim-3 inhibitors in near future.

Structural Characterization of the Interaction of Hypoxia Inducible Factor-1 with Its Hypoxia Responsive Element at the −964G > A Variation Site of the HLA-G Promoter Region

Human Antigen Leukocyte-G (HLA-G) gene encodes an immune checkpoint molecule that has restricted tissue expression in physiological conditions; however, the gene may be induced in hypoxic conditions by the interaction with the hypoxia inducible factor-1 (HIF1). Hypoxia regulatory elements (HRE) located at the HLA-G promoter region and at exon 2 are the major HIF1 target sites. Since the G allele of the −964G > A transversion induces higher HLA-G expression when compared to the A allele in hypoxic conditions, here we analyzed HIF1-HRE complex interaction at the pair-atom level considering both −964G > A polymorphism alleles. Mouse HIF2 dimer crystal (Protein Data Bank ID: 4ZPK) was used as template to perform homology modelling of human HIF1 quaternary structure using MODELLER v9.14. Two 3D DNA structures were built from 5′GCRTG’3 HRE sequence containing the −964G/A alleles using x3DNA. Protein-DNA docking was performed using the HADDOCK v2.4 server, and non-covalent bonds were computed by DNAproDB server. Molecular dynamic simulation was carried out per 200 ns, using Gromacs v.2019. HIF1 binding in the HRE containing −964G allele results in more hydrogen bonds and van der Waals contact formation than HRE with −964A allele. Protein-DNA complex trajectory analysis revealed that HIF1-HRE-964G complex is more stable. In conclusion, HIF1 binds in a more stable and specific manner at the HRE with G allele.