Identification of dengue viral RNA-dependent RNA polymerase inhibitor using computational fragment-based approaches and molecular dynamics study

Exploring pathogenesis and biomarkers through establishment of a rat model of male infertility with liver depression and kidney deficiency

OBJECTIVES

To establish a rat model that accurately replicates the clinical characteristics of male infertility (MI) with Liver Depression and Kidney Deficiency (LD & KD) and investigate the pathogenesis.

METHODS

After subjecting the rats to chronic restraint stress (CRS) and adenine treatment, a series of tests were conducted, including ethological assessments, evaluations of reproductive characteristics, measurements of biochemical parameters, histopathological examinations, and analyses of urinary metabolites. Additionally, bioinformatics predictions were performed for comprehensive analysis.

RESULTS

Compared to the control, the model exhibited significant manifestations of MI with LD & KD, including reduced responsiveness, diminished frequency of capturing estrous female rats, and absence of mounting behavior. Additionally, the kidney coefficient increased markedly, while the coefficients of the testis and epididymis decreased significantly. Sperm counts and viabilities decreased notably, accompanied by an increase in sperm abnormalities. Dysregulation of reproductive hormone levels in the serum was observed, accompanied by an upregulation of proinflammatory cytokines expressions in the liver and kidney, as well as exacerbated oxidative stress in the penile corpus cavernosum and testis. The seminiferous tubules in the testis exhibited a loose arrangement, loss of germ cells, and infiltration of inflammatory cells. Furthermore, utilizing urinary metabolomics and bioinformatics analysis, 5 key biomarkers and 2 crucial targets most closely linked to MI were revealed.

CONCLUSION

The study successfully established a clinically relevant animal model of MI with LD & KD. It elucidates the pathogenesis of the condition, identifies key biomarkers and targets, and provides a robust scientific foundation for the prediction, diagnosis, and treatment of MI with LD & KD.

Novel mutations in structural proteins of dengue virus genomes

BACKGROUND

Genomic characterization of the dengue virus (DENV) is useful for understanding its molecular evolution, transmission, pathogenicity and infectivity. The DENV genomic RNA encodes three structural proteins, capsid (C) envelope (E) and membrane (M) proteins mediating viral entry and assembly during host infection. The current study aims to explore the DENV serotypes and mutations in the E and M proteins.

METHODS

Twenty-three samples of DENV-positive patients were processed and selected for whole genome sequencing (WGS) from the Punjab Province of Pakistan.

RESULTS

Among the 23 WGS, 19 samples showed numerous mutations (BioProject ID PRJNA943555). DENV1 and DENV2 are the most prevalent serotypes. A total of 179 mutations were detected in the E protein, in which K203E, T88A, I114L, and I293T are novel. The I270L, T272A, S273L, and T277A were found in the "kl" β-hairpin (aa 270-279). The M protein harbors 74 mutations, of which 24 were novel. Three prominent complementary regions in the prM and E protein complex formations include R6, E46, D47, D63, and D65 on 'pr' peptide, and E84, K64, and H244, K247 on E, remain conserved except R6C. To our knowledge, it is the first comprehensive study of mutations in structural proteins.

CONCLUSION

Genomic epidemiology is critical for analyzing emerging mutations and designing new policies therapeutic efforts for future outbreaks.

Design, Synthesis, Evaluation and Molecular Dynamics Simulation of Dengue Virus NS5-RdRp Inhibitors

Dengue virus (DENV) is a major mosquito-borne human pathogen in tropical countries; however, there are currently no targeted antiviral treatments for DENV infection. Compounds 27 and 29 have been reported to be allosteric inhibitors of DENV RdRp with potent inhibitory effects. In this study, the structures of compounds 27 and 29 were optimized using computer-aided drug design (CADD) approaches. Nine novel compounds were synthesized based on rational considerations, including molecular docking scores, free energy of binding to receptor proteins, predicted Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) parameters, structural diversity, and feasibility of synthesis. Subsequently, the anti-DENV activity was assessed. In the cytopathic effect (CPE) assay conducted on BHK-21 cells using the DENV2 NGC strain, both SW-b and SW-d demonstrated comparable or superior activity against DENV2, with IC50 values of 3.58 ± 0.29 μM and 23.94 ± 1.00 μM, respectively, compared to that of compound 27 (IC50 = 19.67 ± 1.12 μM). Importantly, both SW-b and SW-d exhibited low cytotoxicity, with CC50 values of 24.65 μmol and 133.70 μmol, respectively, resulting in selectivity indices of 6.89 and 5.58, respectively. Furthermore, when compared to the positive control compound 3'-dATP (IC50 = 30.09 ± 8.26 μM), SW-b and SW-d displayed superior inhibitory activity in an enzyme inhibitory assay, with IC50 values of 11.54 ± 1.30 μM and 13.54 ± 0.32 μM, respectively. Molecular dynamics (MD) simulations elucidated the mode of action of SW-b and SW-d, highlighting their ability to enhance π-π packing interactions between benzene rings and residue W795 in the S1 fragment, compared to compounds 27 and 29. Although the transacylsulphonamide fragment reduced the interaction between T794 and NH, it augmented the interaction between R729 and T794. In summary, our study underscores the potential of SW-b and SW-d as allosteric inhibitors targeting the DENV NS5 RdRp domain. However, further in vivo studies are warranted to assess their pharmacology and toxicity profiles.

Anticorrelated position fluctuation of lipids in forming membrane water pores: molecular dynamics simulations study with dengue virus capsid protein

The traffic of molecules into or out of cells is regulated by many membrane-associated mechanisms. Membrane pores are considered as one of the major passage mechanisms, although molecular-level understanding of pore formation is still vague. The opening of a membrane pore depends on many factors, including the influence of some proteins. The ability of the cell-penetrating peptides and supercharged proteins to form membrane pores has been reported. We studied pore formation through dipalmitoylphosphatidylcholine (DPPC) lipid bilayers by supercharged dengue virus capsid (C) protein. Atomistic molecular dynamics simulations confirmed the formation of membrane pores by a combined effect of the C protein and the membrane electric field. Analyses of simulated trajectories showed highly correlated vertical position fluctuations between the C_α atom of the membrane-anchored arginine residues and the phosphorus atoms of the surrounding DPPC lipids. Certain regions of the bilayer were negatively correlated while the others were positively correlated with respect to the fluctuations of the C_α atom of the anchored arginine residues. When positively correlated lipids in one leaflet vertically aligned with the negatively correlated lipids in the other leaflet, a local anticorrelated region was generated by weakening the bilayer. The membrane pore was always formed close to this anticorrelated region. Once formed, the C protein followed the hydrated pathway provided by the water-filled pores to cross the membrane.Communicated by Ramaswamy H. Sarma.

Examples of Structural Motifs in Viral Genomes and Approaches for RNA Structure Characterization

The relationship between conserved structural motifs and their biological function in the virus replication cycle is the interest of many researchers around the world. RNA structure is closely related to RNA function. Therefore, technological progress in high-throughput approaches for RNA structure analysis and the development of new ones are very important. In this mini review, we discuss a few perspectives on the structural elements of viral genomes and some methods used for RNA structure prediction and characterization. Based on the recent literature, we describe several examples of studies concerning the viral genomes, especially severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza A virus (IAV). Herein, we emphasize that a better understanding of viral genome architecture allows for the discovery of the structure-function relationship, and as a result, the discovery of new potential antiviral therapeutics.

Molecular docking and dynamics studies on propolis sulabiroin-A as a potential inhibitor of SARS-CoV-2

Molecular docking and dynamics simulations were conducted to investigate the antiviral activity of Propolis Sulabiroin-A to inhibit the SARS-CoV-2 virus with quercetin, hesperidin, and remdesivir as control ligands. The parameters calculated were docking score and binding energy/molecular mechanics-generalized born surface area (MMGBSA), root mean square displacement (RMSD), and root mean square fluctuation (RMSF). Docking and MMGBSA scores showed that all the ligands demonstrate an excellent candidate as an inhibitor, and the order of both scores is hesperidin, remdesivir, quercetin, and sulabiroin-A. The molecular dynamics simulation showed that all the ligands are good candidates as inhibitors. Although the fluctuation of Sulabiroin-A is relatively high, it has less protein-ligand interaction time than other ligands. Overall, there is still a good possibility that sulabiroin-A can be used as an alternative inhibitor if a new structure of receptor SARS-CoV-2 is used.

Molecular level insights into the inhibition of gene expression by thiamine pyrophosphate (TPP) analogs for TPP riboswitch: A well-tempered metadynamics simulations study

Riboswitches are metabolite sensing aptamer domains present in non-coding regions in RNA and act as gene-regulating elements. Thiamine pyrophosphate (TPP) riboswitch is evolved as a new target for developing antibiotics against many pathogenic bacteria. The earlier reports suggest that the modification of the pyrophosphate group in the ligand molecule can enhance gene expression. In this work, we have examined the binding affinity and efficacy of TPP and two recently reported ligands, CH₂-TPP, and CF₂-TPP, using Well-tempered metadynamics (WT-MtD) simulations. The experimental in vitro assays show that both TPP and CH₂-TPP repress the gene expression to the same extent. The calculated binding energies correlate well with the experimental study and show the same trend of binding affinity of ligands for the TPP riboswitch. The root mean square fluctuation profiles suggest that the CH₂-TPP and TPP trigger higher fluctuations in P1 and L3 region, and such fluctuations in the P1 region is involved in the gene regulation process. The metal ion mediated contact of TPP ligand with pyrophosphate binding helix is found to be critical in the gene regulation process. The simulation results corroborate the experimental observations that the role of conformational changes occurring in different riboswitch regions upon ligand binding is essential to repress the gene expression process. This work sheds light on the subtle change in the ligand structure that can induce a more considerable impact on binding affinity and efficacy of ligands with riboswitch.

Identification of Potential Inhibitors for Targets Involved in Dengue Fever

Lethality due to dengue infection is a global threat. Nearly 400 million people are affected every year, which approximately costs 500 million dollars for surveillance and vector control itself. Many investigations on the structure-function relationship of proteins expressed by the dengue virus are being made for more than a decade and had come up with many reports on small molecule drug discovery. In this review, we present a detailed note on viral proteins and their functions as well as the inhibitors discovered/designed so far using experimental and computational methods. Further, the phytoconstituents from medicinal plants, specifically the extract of the papaya leaves, neem and bael, which combat dengue infection via dengue protease, helicase, methyl transferase and polymerase are summarized.

Computational approaches for identifying potential inhibitors on targeting protein interactions in drug discovery

In the era of big data, the interplay of artificial and human intelligence is the demanding job to address the concerns involving exchange of decisions between both sides. Drug discovery is one of the key sources of the big data, which involves synergy among various computational methods to achieve a clinical success. Rightful acquisition, mining and analysis of the data related to ligand and targets are crucial to accomplish reliable outcomes in the entire process. Novel designing and screening tactics are necessary to substantiate a potent and efficient lead compounds. Such methods are emphasized and portrayed in the current review targeting protein-ligand and protein-protein interactions involved in various diseases with potential applications.

In vitro and in silico Models to Study Mosquito-Borne Flavivirus Neuropathogenesis, Prevention, and Treatment

Mosquito-borne flaviviruses can cause disease in the nervous system, resulting in a significant burden of morbidity and mortality. Disease models are necessary to understand neuropathogenesis and identify potential therapeutics and vaccines. Non-human primates have been used extensively but present major challenges. Advances have also been made toward the development of humanized mouse models, but these models still do not fully represent human pathophysiology. Recent developments in stem cell technology and cell culture techniques have allowed the development of more physiologically relevant human cell-based models. In silico modeling has also allowed researchers to identify and predict transmission patterns and discover potential vaccine and therapeutic candidates. This review summarizes the research on in vitro and in silico models used to study three mosquito-borne flaviviruses that cause neurological disease in humans: West Nile, Dengue, and Zika. We also propose a roadmap for 21st century research on mosquito-borne flavivirus neuropathogenesis, prevention, and treatment.

Drug repositioning for dengue haemorrhagic fever by integrating multiple omics analyses

To detect drug candidates for dengue haemorrhagic fever (DHF), we employed a computational drug repositioning method to perform an integrated multiple omics analysis based on transcriptomic, proteomic, and interactomic data. We identified 3,892 significant genes, 389 proteins, and 221 human proteins by transcriptomic analysis, proteomic analysis, and human–dengue virus protein–protein interactions, respectively. The drug candidates were selected using gene expression profiles for inverse drug–disease relationships compared with DHF patients and healthy controls as well as interactomic relationships between the signature proteins and chemical compounds. Integrating the results of the multiple omics analysis, we identified eight candidates for drug repositioning to treat DHF that targeted five proteins (ACTG1, CALR, ERC1, HSPA5, SYNE2) involved in human–dengue virus protein–protein interactions, and the signature proteins in the proteomic analysis mapped to significant pathways. Interestingly, five of these drug candidates, valparoic acid, sirolimus, resveratrol, vorinostat, and Y-27632, have been reported previously as effective treatments for flavivirus-induced diseases. The computational approach using multiple omics data for drug repositioning described in this study can be used effectively to identify novel drug candidates.

Avocado (Persea americana) fruit extract (2R,4R)-1,2,4-trihydroxyheptadec-16-yne inhibits dengue virus replication via upregulation of NF-κB-dependent induction of antiviral interferon responses

Dengue virus (DENV) caused millions of infections around the world annually. Co-infection with different serotypes of DENV is associated with dengue hemorrhagic shock syndrome, leading to an estimate of 50% death rate. No approved therapies are currently available for the treatment of DENV infection. Hence, novel anti-DENV agents are urgently needed for medical therapy. Here we demonstrated that a natural product (2 R,4 R)-1,2,4-trihydroxyheptadec-16-yne (THHY), extracted from avocado (Persea americana) fruit, can inhibit DENV-2 replication in a concentration-dependent manner and efficiently suppresses replication of all DENV serotypes (1–4). We further reveal that the NF-κB-mediated interferon antiviral response contributes to the inhibitory effect of THHY on DENV replication. Using a DENV-infected ICR suckling mouse model, we found that THHY treatment caused an increased survival rate among mice infected with DENV. Collectively, these findings support THHY as a potential agent to control DENV infection.

Genomics, proteomics and evolution of dengue virus

The genome of a pathogenic organism possesses a specific order of nucleotides that contains not only information about the synthesis and expression of proteomes, which are required for its growth and survival, but also about its evolution. Inhibition of any particular protein, which is required for the survival of that pathogenic organism, can be used as a potential therapeutic target for the development of effective drugs to treat its infections. In this review, the genomics, proteomics and evolution of dengue virus have been discussed, which will be helpful in better understanding of its origin, growth, survival and evolution, and may contribute toward development of new efficient anti-dengue drugs.

Zika virus NS5 protein potential inhibitors: an enhanced in silico approach in drug discovery

The re-emerging Zika virus (ZIKV) is an arthropod-borne virus that has been described to have explosive potential as a worldwide pandemic. The initial transmission of the virus was through a mosquito vector, however, evolving modes of transmission has allowed the spread of the disease over continents. The virus has already been linked to irreversible chronic central nervous system conditions. The concerns of the scientific and clinical community are the consequences of Zika viral mutations, thus suggesting the urgent need for viral inhibitors. There have been large strides in vaccine development against the virus but there are still no FDA approved drugs available. Rapid rational drug design and discovery research is fundamental in the production of potent inhibitors against the virus that will not just mask the virus, but destroy it completely. In silico drug design allows for this prompt screening of potential leads, thus decreasing the consumption of precious time and resources. This study demonstrates an optimized and proven screening technique in the discovery of two potential small molecule inhibitors of ZIKV Methyltransferase and RNA dependent RNA polymerase. This in silico 'per-residue energy decomposition pharmacophore' virtual screening approach will be critical in aiding scientists in the discovery of not only effective inhibitors of Zika viral targets, but also a wide range of anti-viral agents.

Crystal structure of Zika virus NS5 RNA-dependent RNA polymerase

The current Zika virus (ZIKV) outbreak became a global health threat of complex epidemiology and devastating neurological impacts, therefore requiring urgent efforts towards the development of novel efficacious and safe antiviral drugs. Due to its central role in RNA viral replication, the non-structural protein 5 (NS5) RNA-dependent RNA-polymerase (RdRp) is a prime target for drug discovery. Here we describe the crystal structure of the recombinant ZIKV NS5 RdRp domain at 1.9 Å resolution as a platform for structure-based drug design strategy. The overall structure is similar to other flaviviral homologues. However, the priming loop target site, which is suitable for non-nucleoside polymerase inhibitor design, shows significant differences in comparison with the dengue virus structures, including a tighter pocket and a modified local charge distribution.

Discovery of Potent Non-Nucleoside Inhibitors of Dengue Viral RNA-Dependent RNA Polymerase from a Fragment Hit Using Structure-Based Drug Design

The discovery and optimization of non-nucleoside dengue viral RNA-dependent-RNA polymerase (RdRp) inhibitors are described. An X-ray-based fragment screen of Novartis' fragment collection resulted in the identification of a biphenyl acetic acid fragment 3, which bound in the palm subdomain of RdRp. Subsequent optimization of the fragment hit 3, relying on structure-based design, resulted in a >1000-fold improvement in potency in vitro and acquired antidengue activity against all four serotypes with low micromolar EC50 in cell-based assays. The lead candidate 27 interacts with a novel binding pocket in the palm subdomain of the RdRp and exerts a promising activity against all clinically relevant dengue serotypes.