Systematic prioritization of functional hotspot in RIG-1 domains using pattern based conventional molecular dynamic simulation

Exploring okra-derived compounds as prospective aromatase inhibitors: a computational study for enhanced breast cancer therapy

Estrogen receptor-positive breast cancer represents itself as the most prevalent malignancy among postmenopausal women. One of the promising therapeutic approaches involves the use of Aromatase inhibitors, which competitively bind to Aromatase, reducing estrone and estradiol levels. While current drugs have improved survival rates, they are not without adverse effects. Consequently, this study explores the computational screening of medicinally relevant compounds derived from okra (Abelmoschus esculentus) for potential Aromatase inhibition. Molecular docking employing AMDock v1.5.2 was utilized to assess binding affinities with Aromatase (PDB:3EQM). Subsequently, in-depth molecular interactions were examined using Discovery Studio Visualizer v4.5, and the stability of docked complexes was evaluated via molecular dynamics with the GROMACS package, focusing on RMSD, RMSF, H-bond count, SASA, Free energy landscape, Principal Component Analysis and binding affinity assessment. The pharmacokinetic properties of the okra compounds were predicted using admetSAR v2.0. Our findings highlight Quercetin 3-gentiobioside as a standout candidate, demonstrating superior binding affinity (-10 kcal/mol) and an estimated Ki of 46.77 nM compared to letrozole and other okra compounds. Molecular dynamic analysis confirms the stability of Quercetin 3-gentiobioside binding in terms of H-bonds and conformational integrity. In conclusion, our computational investigation identifies Quercetin 3-gentiobioside, along with Quercetin 3-O-rutinoside and Hyperin, as promising candidates for preclinical studies in the pursuit of potential Aromatase inhibitors.Communicated by Ramaswamy H. Sarma.

Integrating Clinical Cancer and PTM Proteomics Data Identifies a Mechanism of ACK1 Kinase Activation

Beyond the most common oncogenes activated by mutation (mut-drivers), there likely exists a variety of low-frequency mut-drivers, each of which is a possible frontier for targeted therapy. To identify new and understudied mut-drivers, we developed a machine learning (ML) model that integrates curated clinical cancer data and posttranslational modification (PTM) proteomics databases. We applied the approach to 62,746 patient cancers spanning 84 cancer types and predicted 3,964 oncogenic mutations across 1,148 genes, many of which disrupt PTMs of known and unknown function. The list of putative mut-drivers includes established drivers and others with poorly understood roles in cancer. This ML model is available as a web application. As a case study, we focused the approach on nonreceptor tyrosine kinases (NRTK) and found a recurrent mutation in activated CDC42 kinase-1 (ACK1) that disrupts the Mig6 homology region (MHR) and ubiquitin-association (UBA) domains on the ACK1 C-terminus. By studying these domains in cultured cells, we found that disruption of the MHR domain helps activate the kinase while disruption of the UBA increases kinase stability by blocking its lysosomal degradation. This ACK1 mutation is analogous to lymphoma-associated mutations in its sister kinase, TNK1, which also disrupt a C-terminal inhibitory motif and UBA domain. This study establishes a mut-driver discovery tool for the research community and identifies a mechanism of ACK1 hyperactivation shared among ACK family kinases.

IMPLICATIONS

This research identifies a potentially targetable activating mutation in ACK1 and other possible oncogenic mutations, including PTM-disrupting mutations, for further study.

Extensive characterization of 28 complete chloroplast genomes of Hydrangea species: A perspective view of their organization and phylogenetic and evolutionary relationships

The tribe Hydrangeeae displays a unique, distinctive disjunct distribution encompassing East Asia, North America and Hawaii. Despite its complex trait variations and polyphyletic nature, comprehensive phylogenomic and biogeographical studies on this tribe have been lacking. To address this gap, we sequenced and characterized 28 plastomes of Hydrangeeae. Our study highlights the highly conserved nature of Hydrangeaceae chloroplast (cp) genomes in terms of gene content and arrangement. Notably, synapomorphic characteristics of tandem repeats in the conserved domain of accD were observed in the Macrophyllae, Chinenses, and Dichroa sections within the Hydrangeeae tribe. Additionally, we found lower expression of accD in these sections using structure prediction and quantitative real-time PCR analysis. Phylogenomic analyses revealed the subdivision of the Hydrangeeae tribe into two clades with robust support values. Consistent with polyphyletic relationships, sect. Broussaisia was identified as the basal group in the tribe Hydrangeeae. Our study also provides insights into the phylogenetic relationships of Hydrangea petiolaris in the Jeju and Ulleung Island populations, suggesting the need for further studies with more samples and molecular data. Divergence time estimation and biogeographical analyses suggested that the common ancestors of the tribe Hydrangeeae likely originated from North America and East Asia during the Paleocene period via the Bering Land Bridge, potentially facilitating migration within the tribe between these regions. In conclusion, this study enhances our understanding of the evolutionary history and biogeography of the tribe Hydrangeeae, shedding light on the dispersal patterns and origins of this intriguing plant group with its unique disjunct distribution.

Structural localization of pathogenic mutations in the central nucleotide-binding domain (NBD) of nucleotide-binding oligomerization domain-2 (NOD2) protein and their inference in inflammatory disorders

The human NBD domain which is centrally located in the NOD2 protein displays an essential role in oligomerization and initiates the immune response via CARD-RIPK2 interaction. The mutations associated with the NBD domain have been largely implicated in inflammatory disorders such as Blau syndrome and sarcoidosis. This study aims to determine the structural and phenotypic effect of a lethal mutation that occurs in the NBD domain which has an axiomatic impact on protein dysfunction. Initially, the most deleterious missense mutations were screened through various in silico analysis. Out of 33 variants, I-Mutant 3.0, SIFT, PolyPhen 2, Align GVGD, PHD SNP and SNP&GO have statistically identified 5 variants (R42W, D90E, E91K, G189D & W198L) as less stable, deleterious and damaging. Our predicted models have paved the way to understand the various structural properties such as physiochemical, secondary structural arrangements and stabilizing residues in folding associated with the native and mutant NBD domain especially of the functionally important regions. From the aforementioned results, R42W and G189D were found to be the more predominant among the mutants. Precisely, through molecular simulation, we have strongly justified the significant conformational disruption of R42W and G189D through the stabilization factors, folding and essential dynamics. Conclusively, these regions (α₃41-44, α₁₃185-191 and β₆133-143β₇) seem to adopt such structures that are not conducive to wild-type-like functionality. Our prediction and validation of lethal mutations based on structural stability may be useful for conducting experimental studies in detail to uncover the protein deregulation leading to inflammatory disorders.

Identification of metal binding motifs in protein frameworks to develop novel remediation strategies for Hg2+ and Cr(VI)

Amino acid sequences in metal-binding proteins with chelating properties offer exciting applications in biotechnology and medical research. To enhance their application in bioremediation studies, we explicitly aimed to identify specific metal-binding chelating motifs in protein structures for two significant pollutants, such as mercury (Hg²⁺) and chromium Cr(V1). For this purpose, we have performed an extensive coordination chemistry approach by retrieving Hg²⁺ and Cr(V1) binding protein structures from the protein database and validated using the B-factor, a term defining uncertainty of the atoms and with occupancy to obtain the best binding motifs. Our analysis revealed that acidic amino acids like aspartic acid, glutamic acid, and basic amino acids such as cysteine and histidine are predominant in coordinating with these metals. The order of preference in Hg²⁺-bound structures is predicted to be Cys > His > Asp > Glu, and for Cr(V1) is His > Asp > Glu. Examination of the atomic coordinates and their distance from each metal revealed that the sulfur atoms of cysteine showing more preference towards Hg²⁺coordination with an atomic distance ranging from 1.5 to 2.9 Å. Likewise, oxygen atoms of aspartic acid, glutamic acid and nitrogen atoms of histidine are within 2 Å of Cr(V1) coordination. Based on these observations, we obtained C-C-C, C-X(2)-C-C-(X)2-C, H-C-H motifs for Hg²⁺, and D-X(1)-D, H-X(3)-E motif for Cr(V1) to be shared within the coordination space of 3 Å. As a future scope, we propose that the identified metal-binding chelating motifs are oligopeptides and can display on the surface of microorganisms such as Escherichia coli and Saccharomyces cerevisiae for effective removal of natural Hg²⁺ and Cr(V1) through biosorption. Hence, our results will provide the basis for futuristic bioremediation.

Bioactivity of melianone against Salmonella and in silico prediction of a membrane protein target

Melianone, the protolimonoid (24, 25-epoxyflindissone), was isolated from the medicinal tree species, Swietenia mahagoni (L.) JACQ (Meliaceae). The compound isolated from petroleum ether leaf extracts (5.39%) was quantified using high-performance thin-layer chromatography (HPTLC) method. In antimicrobial assays melianone inhibited Salmonella ser. Typhi with an MIC of 0.053 µM. Induced Fit Docking (IFD) of the ligand, melianone, with proteins involved in anaerobic virulence of the pathogen, revealed that it binds with FocA (a transport protein of formate ions) at its "periplasmic opening" with a glide energy of - 51.8576 kcal mol^-1. Melianone altered the overall conformation of the protein (protomer A) by 0.347 Å RMSD. It induced a notable protein topology (Ω loop region) shift in the channel from an intermediate-open to a closed-state conformation and was supported by molecular dynamic simulations performed. FocA, a protein that contributes to its survival under anaerobic conditions, was further evaluated experimentally, after exposure of Salmonella ser. Typhi to melianone, resulting in the altered homeostasis of formate.

A molecular simulation analysis of vitamin D targets interleukin 13 (IL13) as an alternative to mometasone in asthma

Asthma, a chronic lung disease characterized by obstruction of airway passage is characterized by inflammation and hyperresponsiveness with increase in the number of eosinophils. Interleukin-13, plays a significant role in causing inflammation during an asthmatic attack by bronchial constriction. Mometasone, a glucocorticoid has been used as the first line of administration for people affected with asthma for almost a decade. However, in several cases, people treated with mometasone have faced systemic and local side effects. To reduce these side effects, we hypothesized vitamin D that can be used as a substitute to mometasone. For this purpose, we employed the use of molecular docking and simulation studies for comparative study. The docking studies revealed the binding residues of interleukin-13 which are bound to the active site. Among all, we noticed three binding residue Leu83, His84 and Arg86 common for both mometasone and vitamin D. Also, the binding energies share a significant similarity between them. The docked complexes of mometasone and vitamin D with interleukin-13 were evaluated with molecular dynamics simulation. Consistently, the MD analysis uncovered the interesting note on conformational adaptation between the complexes as well as that vitamin D has the complementary binding efficiency to interleukin-13 as compared to mometasone. The substitution of vitamin D might provide a promising gateway to reduce the side effects caused by mometasone and also reduce the cost for treatment of asthma patients.

Nearly neutral evolution in IFNL3 gene retains the immune function to detect and clear the viral infection in HCV

IFNL3 gene plays a crucial role in immune defense against viruses. It induces the interferon stimulated genes (ISGs) with antiviral properties by activating the JAK-STAT pathway. In this study, we investigated the evolutionary force involved in shaping the IFNL3 gene to perform its downstream function as a regulatory gene in HCV clearance. We have selected 25 IFNL3 coding sequences with human gene as a reference sequence and constructed a phylogeny. Furthermore, rate of variation, substitution saturation test, phylogenetic informativeness and differential selection were also analysed. The codon evolution result suggests that nearly neutral mutation is the key pattern in shaping the IFNL3 evolution. The results were validated by subjecting the human IFNL3 protein variants to that of the native through a molecular dynamics simulation study. The molecular dynamics simulation clearly depicts the negative impact on the reported variants in human IFNL3 protein. However, these detrimental mutations (R157Q and R157W) were shown to be negatively selected in the evolutionary study of the mammals. Hence, the variation revealed a mild impact on the IFNL3 function and may be removed from the population through negative selection due to its high functional constraints. In a nutshell, our study may contribute the overall evidence in phylotyping and structural transformation that takes place in the non-synonymous substitutions of IFNL3 protein. Substantially, our obtained theoretical knowledge will lay the path to extend the experimental validation in HCV clearance.