Therapeutic p28 peptide targets essential H1N1 influenza virus proteins: insights from docking and molecular dynamics simulations.

The H1N1 influenza virus causes a severe disease that affects the human respiratory tract leading to millions of deaths every year. At present, certain vaccines and few drugs are used to control the virus during seasonal outbreaks. However, high mutation rates and genetic reassortment make it challenging to prevent and mitigate outbreaks, leading to pandemics. Thus, alternate therapies are required for its management and control. Here, we report that a bacterial protein, azurin, and its peptide derivatives p18 and p28 target critical proteins of the influenza virus in an effective manner. The molecular docking studies show that the p28 peptide could target C-PB1, NS1-ED, PB2-CBD, PB2-RBD, NP, and PA proteins. These complexes were further subjected to the simulation of molecular dynamics and binding free energy calculations. The data indicate that p28 has an unusually high affinity and forms stable complexes with the viral proteins C-PB1, PB2-CBD, PB2-RBD, and NP. We suggest that the azurin derivative p28 peptide can act as an anti-influenza agent as it can bind to multiple targets and neutralize the virus. Additional experimental studies need to be conducted to evaluate its safety and efficacy as an anti-H1N1 molecule.

A Combined Molecular Docking/Dynamics Approach to Probe the Binding Mode of Cancer Drugs with Cytochrome P450 3A4.

Cytarabine, daunorubicin, doxorubicin and vincristine are clinically used for combinatorial therapies of cancers in different combinations. However, the knowledge about the interaction of these drugs with the metabolizing enzyme cytochrome P450 is limited. Therefore, we utilized computational methods to predict and assess the drug-binding modes. In this study, we performed docking, MD simulations and free energy landscape analysis to understand the drug-enzyme interactions, protein domain motions and the most populated free energy minimum conformations of the docked protein-drug complexes, respectively. The outcome of docking and MD simulations predicted the productive, as well as the non-productive binding modes of the selected drugs. Based on these interaction studies, we observed that S119, R212 and R372 are the major drug-binding residues in CYP3A4. The molecular mechanics Poisson-Boltzmann surface area analysis revealed the dominance of hydrophobic forces in the CYP3A4-drug association. Further analyses predicted the residues that may contain favorable drug-specific interactions. The probable binding modes of the cancer drugs from this study may extend the knowledge of the protein-drug interaction and pave the way to design analogs with reduced toxicity. In addition, they also provide valuable insights into the metabolism of the cancer drugs.

Computational insights into potential marine natural products as selective inhibitors of Mycobacterium tuberculosis InhA: A structure-based virtual screening study.

Several factors are associated with the emergence of drug resistance mechanisms, such as impermeable cell walls, gene mutations, and drug efflux systems. Consequently, bacteria acquire resistance, leading to a decrease in drug efficacy. A new and innovative strategy is required to combat drug resistance in tuberculosis (TB) effectively. Therefore, targeting the mycolic acid biosynthesis pathway, which is involved in synthesising mycolic acids (MAs), essential structural components responsible for mycobacterial pathogenicity, has garnered interest in TB research and the concept of drug resistance. In this context, InhA, which plays a crucial role in the fatty acid synthase-II (FAS-II) system of the MA biosynthetic pathway, was selected as a druggable target for screening investigation. To identify potential lead molecules against InhA, diverse marine natural products (MNPs) were collected from the comprehensive marine natural products database (CMNPD). Virtual screening studies aided in selecting potential lead molecules that best fit within the substrate-binding pocket (SBP) of InhA, forming crucial hydrogen bond interaction with the catalytic residue Tyr158. Three MNPs, CMNPD30814, CMNPD1702, and CMNPD27355, were chosen as prospective alternative molecules due to their favorable pharmacokinetic properties and lack of toxicity according to ProTox-II predictions. Additionally, improved reactivity of the MNPs was observed in the results of density functional theory (DFT) studies. Furthermore, comparative molecular dynamics simulation (MDS), principal component (PC)-based free energy landscape (FEL) analysis, and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) were employed to show enhanced structural stability, increased H-bond potential, and high binding affinity toward the target InhA. Moreover, the hot spot residues that contributed to the high binding energy profile and anchored the stability of the complexes were revealed with their individual interaction energy. The computational insights from this study provide potential avenues to combat TB through the multifaceted mode of action of these marine lead molecules, which can be further explored in future experimental investigations.

Impact of nsSNPs in human AIM2 and IFI16 gene: a comprehensive in silico analysis.

AIM2 and IFI16 are the most studied members of AIM2-like receptors (ALRs) in humans and share a common N-Terminal PYD domain and C-terminal HIN domain. The HIN domain binds to dsDNA in response to the invasion of bacterial and viral DNA, and the PYD domain directs apoptosis-associated speck-like protein via protein-protein interactions. Hence, activation of AIM2 and IFI16 is crucial for protection against pathogenic assaults, and any genetic variation in these inflammasomes can dysregulate the human immune system. In this study, different computational tools were used to identify the most deleterious and disease-causing non-synonymous single nucleotide polymorphisms (nsSNPs) in AIM2 and IFI16 proteins. Molecular dynamic simulation was performed for the top damaging nsSNPs to study single amino acid substitution-induced structural alterations in AIM2 and IFI16. The observed results suggest that the variants G13V, C304R, G266R, and G266D for AIM2, and G13E and C356F are deleterious and affect structural integrity. We hope that the suggested deleterious nsSNPs and structural dynamics of AIM2 and IFI16 variants will guide future research to better understand the function of these variants with large-scale studies and may assist in fresher therapeutics focusing on these polymorphisms.Communicated by Ramaswamy H. Sarma.

Deciphering the intrinsic dynamics of unphosphorylated IRAK4 kinase bound to type I and type II inhibitors.

Interleukin-1 receptor-associated kinase 4 (IRAK4) is a vital protein involved in Toll-like and interleukin-1 receptor signal transduction. Several studies have reported regarding the crystal structure, dynamic properties, and interactions with inhibitors of the phosphorylated form of IRAK4. However, no dynamic properties of inhibitor-bound unphosphorylated IRAK4 have been previously studied. Herein, we report the intrinsic dynamics of unphosphorylated IRAK4 (uIRAK4) bound to type I and type II inhibitors. The corresponding apo and inhibitor-bound forms of uIRAK4 were subjected to three independent simulations of 500 ns (total 1.5 µs) each, and their trajectories were analyzed. The results indicated that all three systems were relatively stable, except for the type II inhibitor-bound form of uIRAK4, which exhibited less compact folding and higher solvent surface area. The intra-hydrogen bonds corroborated the structural deformation of the type-II inhibitor-bound complex, which could be attributed to the long molecular structure of the type-II inhibitor. Moreover, the type II inhibitor bound to uIRAK4 showed higher binding free energy with uIRAK4 than the type I inhibitor. The free energy landscape analysis showed a reorientation of Phe330 side chain from the DFG motif at different metastable states for all the systems. The intra-residual distance between residues Lys213, Glu233, Tyr262, and Phe330 suggests a functional interplay when the inhibitors are bound to uIRAK4, thereby hinting at their crucial role in the inhibition mechanism. Ultimately, the intrinsic dynamics study observed between type I/II inhibitor-bound forms of uIRAK4 may assist in better understanding the enzyme and designing therapeutic compounds.

Molecular dynamics of the ERRγ ligand-binding domain bound with agonist and inverse agonist.

Estrogen-related receptor gamma (ERRγ), the latest member of the ERR family, does not have any known reported natural ligands. Although the crystal structures of the apo, agonist-bound, and inverse agonist-bound ligand-binding domain (LBD) of ERRγ have been solved previously, their dynamic behavior has not been studied. Hence, to explore the intrinsic dynamics of the apo and ligand-bound forms of ERRγ, we applied long-range molecular dynamics (MD) simulations to the crystal structures of the apo and ligand-bound forms of the LBD of ERRγ. Using the MD trajectories, we performed hydrogen bond and binding free energy analysis, which suggested that the agonist displayed more hydrogen bonds with ERRγ than the inverse agonist 4-OHT. However, the binding energy of 4-OHT was higher than that of the agonist GSK4716, indicating that hydrophobic interactions are crucial for the binding of the inverse agonist. From principal component analysis, we observed that the AF-2 helix conformation at the C-terminal domain was similar to the initial structures during simulations, indicating that the AF-2 helix conformation is crucial with respect to the agonist or inverse agonist for further functional activity of ERRγ. In addition, we performed residue network analysis to understand intramolecular signal transduction within the protein. The betweenness centrality suggested that few of the amino acids are important for residue signal transduction in apo and ligand-bound forms. The results from this study may assist in designing better therapeutic compounds against ERRγ associated diseases.

Rising Influence of Nanotechnology in Addressing Oxidative Stress-Related Liver Disorders.

Reactive oxygen species (ROS) play a significant role in the survival and decline of various biological systems. In liver-related metabolic disorders such as steatohepatitis, ROS can act as both a cause and a consequence. Alcoholic steatohepatitis (ASH) and non-alcoholic steatohepatitis (NASH) are two distinct types of steatohepatitis. Recently, there has been growing interest in using medications that target ROS formation and reduce ROS levels as a therapeutic approach for oxidative stress-related liver disorders. Mammalian systems have developed various antioxidant defenses to protect against excessive ROS generation. These defenses modulate ROS through a series of reactions, limiting their potential impact. However, as the condition worsens, exogenous antioxidants become necessary to control ROS levels. Nanotechnology has emerged as a promising avenue, utilizing nanocomplex systems as efficient nano-antioxidants. These systems demonstrate enhanced delivery of antioxidants to the target site, minimizing leakage and improving targeting accuracy. Therefore, it is essential to explore the evolving field of nanotechnology as an effective means to lower ROS levels and establish efficient therapeutic interventions for oxidative stress-related liver disorders.

Therapeutic applications of nucleic acids and their analogues in Toll-like receptor signaling.

Toll-like receptors (TLRs) belong to a family of innate immune receptors that detect and clear invading microbial pathogens. Specifically intracellular TLRs such as TLR3, TLR7, TLR8 and TLR9 recognize nucleic acids such as double-stranded RNA, single-stranded RNA and CpG DNA respectively derived from microbial components. Upon infection, nucleic acid sensing TLRs signal within endosomal compartment triggering the induction of essential proinflammatory cytokines and type I interferons to initiate innate immune responses thereby leading to a critical role in the development of adaptive immune responses. Thus, stimulation of TLRs by nucleic acids is a promising area of research for the development of novel therapeutic strategies against pathogenic infection, allergies, malignant neoplasms and autoimmunity. This review summarizes the therapeutic applications of nucleic acids or nucleic acid analogues through the modulation of TLR signaling pathways.

Molecular modeling and dynamic simulation of chicken Mx protein with the S631N polymorphism.

Myxovirus resistance (Mx) proteins are antiviral GTPases induced by type I interferons (IFNs). In chickens, a single Mx protein variant, S631N, has been suggested to possess antiviral activity. However, the impact of this variant on chicken Mx (chMx) protein structure and conformation has not been investigated. Hence, in this study, we applied computational methods such as molecular modeling, molecular dynamic simulation, inter domain motion and residue networks to examine the structure and dynamic behavior of wild-type and mutant chMx. At first, we built 3-dimensional structural models for both wild-type and mutant chMx proteins, which revealed that the structural organization of chMx was similar to that of human Mx proteins. Subsequently, molecular dynamics simulations revealed that angle variation around the hinge1 region led to the different stalk domain conformations between the wild-type and mutant chMx proteins. Domain motion analysis further suggested that the conformational differences in the loop region surrounded by the mutant residue may lead to an inclined stalk domain conformation in the mutant compared to the wild-type protein. In addition, we performed betweenness centrality analysis from residue interaction networks, to identify the crucial residues for intramolecular signal flow in chMx. The results of this study provided information on the differences in structure and dynamics between wild-type and mutant chMx, which may aid in understanding the structural features of the S631N mutant, that may be associated with chMx protein antiviral activity.Communicated by Ramaswamy H. Sarma.

Computational Insights into the Structural Dynamics of MDA5 Variants Associated with Aicardi-Goutières Syndrome and Singleton-Merten Syndrome.

Melanoma differentiation-associated protein 5 (MDA5) is a crucial RIG-I-like receptor RNA helicase enzyme encoded by IFIH1 in humans. Single nucleotide polymorphisms in the IFIH1 results in fatal genetic disorders such as Aicardi-Goutières syndrome and Singleton-Merten syndrome, and in increased risk of type I diabetes in humans. In this study, we chose four different amino acid substitutions of the MDA5 protein responsible for genetic disorders: MDA5L372F, MDA5A452T, MDA5R779H, and MDA5R822Q and analyzed their structural and functional relationships using molecular dynamic simulations. Our results suggest that the mutated complexes are relatively more stable than the wild-type MDA5. The radius of gyration, interaction energies, and intra-hydrogen bond analysis indicated the stability of mutated complexes over the wild type, especially MDA5L372F and MDA5R822Q. The dominant motions exhibited by the wild-type and mutant complexes varied significantly. Moreover, the betweenness centrality of the wild-type and mutant complexes showed shared residues for intra-signal propagation. The observed results indicate that the mutations lead to a gain of function, as reported in previous studies, due to increased interaction energies and stability between RNA and MDA5 in mutated complexes. These findings are expected to deepen our understanding of MDA5 variants and may assist in the development of relevant therapeutics against the disorders.

Conformational Changes Induced by S34Y and R98C Variants in the Death Domain of Myd88.

Myeloid differentiating factor 88 (Myd88) is a universal adaptor protein that plays a critical role in innate immunity by mediating TLR downstream signaling. Myd88 death domain (DD) forms an oligomeric complex by association with other DD-containing proteins such as IRAK4. Despite its universal role, polymorphisms in Myd88 can result in several diseases. Previous studies have suggested that, out of several non-synonymous single-nucleotide polymorphisms (nsSNPs), the variants S34Y and R98C in the DD of Myd88 disrupt the formation of the Myddosome complex. Therefore, we performed molecular dynamics (MD) simulations on wild-type (Myd88WT) and mutant (Myd88S34Y, Myd88R98C) DDs to evaluate the subtle conformational changes induced by these mutations. Our results suggest that the S34Y variant induces large structural transitions compared to the R98C variant as evidenced by residual flexibility at the variable loop regions, particularly in the H1-H2 loop, and variations in the collective modes of motion observed for wild-type and mutant Myd88 DDs. The residue interaction network strongly suggests a distortion in the interaction pattern at the location of the mutated residue between the wild type and mutants. Moreover, betweenness centrality values indicate that variations in the distribution of functionally important residues may be reflected by distinct residue signal transductions in both wild-type and mutant Myd88 DDs, which may influence the interaction with other DDs in TLR downstream signaling.

Empirical Comparison and Analysis of Web-Based DNA N 4-Methylcytosine Site Prediction Tools.

DNA N 4-methylcytosine (4mC) is a crucial epigenetic modification involved in various biological processes. Accurate genome-wide identification of these sites is critical for improving our understanding of their biological functions and mechanisms. As experimental methods for 4mC identification are tedious, expensive, and labor-intensive, several machine learning-based approaches have been developed for genome-wide detection of such sites in multiple species. However, the predictions projected by these tools are difficult to quantify and compare. To date, no systematic performance comparison of 4mC tools has been reported. The aim of this study was to compare and critically evaluate 12 publicly available 4mC site prediction tools according to species specificity, based on a huge independent validation dataset. The tools 4mCCNN (Escherichia coli), DNA4mC-LIP (Arabidopsis thaliana), iDNA-MS (Fragaria vesca), DNA4mC-LIP and 4mCCNN (Drosophila melanogaster), and four tools for Caenorhabditis elegans achieved excellent overall performance compared with their counterparts. However, none of the existing methods was suitable for Geoalkalibacter subterraneus, Geobacter pickeringii, and Mus musculus, thereby limiting their practical applicability. Model transferability to five species and non-transferability to three species are also discussed. The presented evaluation will assist researchers in selecting appropriate prediction tools that best suit their purpose and provide useful guidelines for the development of improved 4mC predictors in the future.

Insights into the dynamic nature of the dsRNA-bound TLR3 complex.

Toll-like receptor 3 (TLR3), an endosomal receptor crucial for immune responses upon viral invasion. The TLR3 ectodomain (ECD) is responsible for double-stranded RNA (dsRNA) recognition and mutational analysis suggested that TLR3 ECD C-terminal dimerization is essential for dsRNA binding. Moreover, the L412F polymorphism of TLR3 is associated with human diseases. Although the mouse structure of the TLR3-dsRNA complex provides valuable insights, the structural dynamic behavior of the TLR3-dsRNA complex in humans is not completely understood. Hence, in this study, we performed molecular dynamic simulations of human wild-type and mutant TLR3 complexes. Our results suggested that apoTLR3 ECD dimers are unlikely to be stable due to the distance between the monomers are largely varied during simulations. The observed interaction energies and hydrogen bonds in dsRNA-bound TLR3 wild-type and mutant complexes indicate the presence of a weak dimer interface at the TLR3 ECD C-terminal site, which is required for effective dsRNA binding. The L412F mutant exhibited similar dominant motion compared to wild-type. Additionally, we identified the distribution of crucial residues for signal propagation in TLR3-dsRNA complex through the evaluation of residue betweenness centrality (CB). The results of this study extend our understanding of TLR3-dsRNA complex, which may assist in TLR3 therapeutics.

Polymodal Transient Receptor Potential Vanilloid Type 1 Nocisensor: Structure, Modulators, and Therapeutic Applications.

Transient receptor potential (TRP) channels belong to a superfamily of sensory-related ion channels responding to a wide variety of thermal, mechanical, or chemical stimuli. In an attempt to comprehend the piquancy and pain mechanism of the archetypal vanilloids, transient receptor potential vanilloid (TRPV) 1 was discovered. TRPV1, a well-established member of the TRP family, is implicated in a range of functions including inflammation, painful stimuli sensation, and mechanotransduction. TRPV1 channels are nonselective cation receptors that are gated by a broad array of noxious ligands. Such polymodal-sensor aspect makes the TRPV1 channel extremely versatile and important for its role in sensing burning pain. Besides ligands, TRPV1 signaling can also be modulated by lipids, secondary messengers, protein kinases, cytoskeleton, and several other proteins. Due to its central role in hyperalgesia transduction and inflammatory processes, it is considered as the primary pharmacological pain target. Moreover, understanding the structural and functional intricacies of the channel is indispensable for the therapeutic intervention of TRPV1 in pain and other pathological disorders. In this chapter, we seek to give a mechanistic outlook on the TRPV1 channel. Specifically, we will explore the TRPV1 structure, activation, modulation, ligands, and its therapeutic targeting. However, the major objective of this review is to highlight the fact that TRPV1 channel can be treated as an effective therapeutic target for treating several pain- and nonpain-related physiological and pathological states.

Role of computer-aided drug design in modern drug discovery.

Drug discovery utilizes chemical biology and computational drug design approaches for the efficient identification and optimization of lead compounds. Chemical biology is mostly involved in the elucidation of the biological function of a target and the mechanism of action of a chemical modulator. On the other hand, computer-aided drug design makes use of the structural knowledge of either the target (structure-based) or known ligands with bioactivity (ligand-based) to facilitate the determination of promising candidate drugs. Various virtual screening techniques are now being used by both pharmaceutical companies and academic research groups to reduce the cost and time required for the discovery of a potent drug. Despite the rapid advances in these methods, continuous improvements are critical for future drug discovery tools. Advantages presented by structure-based and ligand-based drug design suggest that their complementary use, as well as their integration with experimental routines, has a powerful impact on rational drug design. In this article, we give an overview of the current computational drug design and their application in integrated rational drug development to aid in the progress of drug discovery research.

Structural dynamic analysis of apo and ATP-bound IRAK4 kinase.

Interleukin-1 receptor-associated kinases (IRAKs) are Ser/Thr protein kinases that play an important role as signaling mediators in the signal transduction facilitated by the Toll-like receptor (TLR) and interleukin-1 receptor families. Among IRAK family members, IRAK4 is one of the drug targets for diseases related to the TLR and IL-1R signaling pathways. Experimental evidence suggests that the IRAK4 kinase domain is phosphorylated in its activation loop at T342, T345, and S346 in the fully activated state. However, the molecular interactions of subdomains within the active and inactive IRAK4 kinase domain are poorly understood. Hence, we employed a long-range molecular dynamics (MD) simulation to compare apo IRAK4 kinase domains (phosphorylated and unphosphorylated) and ATP-bound phosphorylated IRAK4 kinase domains. The MD results strongly suggested that lobe uncoupling occurs in apo unphosphorylated IRAK4 kinase via the disruption of the R334/T345 and R310/T345 interaction. In addition, apo unphosphorylated trajectory result in high mobility, particularly in the N lobe, activation segment, helix αG, and its adjoining loops. The Asp-Phe-Gly (DFG) and His-Arg-Asp (HRD) conserved kinase motif analysis showed the importance of these motifs in IRAK4 kinase activation. This study provides important information on the structural dynamics of IRAK4 kinase, which will aid in inhibitor development.

Multiple roles of toll-like receptor 4 in colorectal cancer.

Toll-like receptor (TLR) signaling has been implicated in the inflammatory responses in intestinal epithelial cells (IECs). Such inflammatory signals mediate complex interactions between commensal bacteria and TLRs and are required for IEC proliferation, immune response, repair, and homeostasis. The upregulation of certain TLRs in colorectal cancer (CRC) tissues suggests that TLRs may play an essential role in the prognosis of chronic and inflammatory diseases that ultimately culminate in CRC. Here, we provide a comprehensive review of the literature on the involvement of the TLR pathway in the initiation, progression, and metastasis of CRC, as well as inherited genetic variation and epigenetic regulation. The differential expression of TLRs in epithelial cells has also been discussed. In particular, we emphasize the physiological role of TLR4 in CRC development and pathogenesis, and propose novel and promising approaches for CRC therapeutics with the aid of TLR ligands.

Evolutionary, structural and functional interplay of the IκB family members.

A primary level of control for nuclear factor kappa B (NF-κB) is effected through its interactions with the inhibitor protein, inhibitor of kappa B (IκB). Several lines of evidence confirm the existence of multiple forms of IκB that appear to regulate NF-κB by distinct mechanisms. Therefore, we performed a comprehensive bioinformatics analysis to understand the evolutionary history and intrinsic functional diversity of IκB family members. Phylogenetic relationships were constructed to trace the evolution of the IκB family genes. Our phylogenetic analysis revealed 10 IκB subfamily members that clustered into 5 major clades. Since the ankyrin (ANK) domain appears to be more ancient than the Rel homology domain (RHD), our phylogenetic analysis suggests that some undefined ancestral set of ANK repeats acquired an RHD before any duplication and was later duplicated and then diverged into the different IκB subfamilies. Functional analysis identified several functionally divergent sites in the ANK repeat domains (ARDs) and revealed that this region has undergone strong purifying selection, suggesting its functional importance in IκB genes. Structural analysis showed that the major variations in the number of ANK repeats and high conformational changes in the finger loop ARD region contribute to the differing binding partner specificities, thereby leading to distinct IκB functions. In summary, our study has provided useful information about the phylogeny and structural and functional divergence of the IκB family. Additionally, we identified a number of amino acid sites that contribute to the predicted functional divergence of these proteins.

Molecular evolution and structural features of IRAK family members.

The interleukin-1 receptor-associated kinase (IRAK) family comprises critical signaling mediators of the TLR/IL-1R signaling pathways. IRAKs are Ser/Thr kinases. There are 4 members in the vertebrate genome (IRAK1, IRAK2, IRAKM, and IRAK4) and an IRAK homolog, Pelle, in insects. IRAK family members are highly conserved in vertebrates, but the evolutionary relationship between IRAKs in vertebrates and insects is not clear. To investigate the evolutionary history and functional divergence of IRAK members, we performed extensive bioinformatics analysis. The phylogenetic relationship between IRAK sequences suggests that gene duplication events occurred in the evolutionary lineage, leading to early vertebrates. A comparative phylogenetic analysis with insect homologs of IRAKs suggests that the Tube protein is a homolog of IRAK4, unlike the anticipated protein, Pelle. Furthermore, the analysis supports that an IRAK4-like kinase is an ancestral protein in the metazoan lineage of the IRAK family. Through functional analysis, several potentially diverged sites were identified in the common death domain and kinase domain. These sites have been constrained during evolution by strong purifying selection, suggesting their functional importance within IRAKs. In summary, our study highlighted the molecular evolution of the IRAK family, predicted the amino acids that contributed to functional divergence, and identified structural variations among the IRAK paralogs that may provide a starting point for further experimental investigations.