Multidomain Protein Research Articles

EXTH-21. FKBP38 REGULATES SELF-RENEWAL AND SURVIVAL OF GBM NEUROSPHERES

Abstract Glioblastoma (GBM) is the most common malignant primary brain tumor. The outcome is dismal, despite the multimodal therapeutic approach that includes surgical resection followed by radiation and chemotherapy. The quest for novel therapeutic targets to treat GBM is underway. FKBP38, a member of the immunophilin family of proteins, is a multidomain protein that plays an important role in the regulation of cellular functions including apoptosis and autophagy. In this study, we tested the role of FKBP38 in GBM tumor biology. Expression of FKBP38 is upregulated in the patient-derived primary GBM neurospheres (GBMNS) compared to normal human astrocytes. Attenuation of FKBP38 expression decreased the viability of GBMNS and increased the caspase 3/7 activity indicating that FKBP38 is required for the survival of GBMNS. Further, the depletion of FKBP38 significantly reduced the number of neurospheres that were formed implying that FKBP38 regulates the self-renewal of GBMNS. Also, the transient knockdown of FKBP38 increased the LC3-II/I ratio suggesting the induction of autophagy with the depletion of FKBP38. Further investigation showed that the negative regulation of autophagy by FKBP38 in GBMNS is mediated through the JNK/C-Jun-PTEN-AKT pathway. In vivo, FKBP38 depletion significantly extended the survival of tumor-bearing mice. Overall, our results suggest targeting FKBP38 imparts an anti-GBM effect by inducing apoptosis and autophagy and thus can be a potential therapeutic target for GBM therapy.

FKBP38 Regulates Self-Renewal and Survival of GBM Neurospheres.

Glioblastoma is the most common malignant primary brain tumor. The outcome is dismal, despite the multimodal therapeutic approach that includes surgical resection, followed by radiation and chemotherapy. The quest for novel therapeutic targets to treat glioblastoma is underway. FKBP38, a member of the immunophilin family of proteins, is a multidomain protein that plays an important role in the regulation of cellular functions, including apoptosis and autophagy. In this study, we tested the role of FKBP38 in glioblastoma tumor biology. Expression of FKBP38 was upregulated in the patient-derived primary glioblastoma neurospheres (GBMNS), compared to normal human astrocytes. Attenuation of FKBP38 expression decreased the viability of GBMNSs and increased the caspase 3/7 activity, indicating that FKBP38 is required for the survival of GBMNSs. Further, the depletion of FKBP38 significantly reduced the number of neurospheres that were formed, implying that FKBP38 regulates the self-renewal of GBMNSs. Additionally, the transient knockdown of FKBP38 increased the LC3-II/I ratio, suggesting the induction of autophagy with the depletion of FKBP38. Further investigation showed that the negative regulation of autophagy by FKBP38 in GBMNSs is mediated through the JNK/C-Jun-PTEN-AKT pathway. In vivo, FKBP38 depletion significantly extended the survival of tumor-bearing mice. Overall, our results suggest that targeting FKBP38 imparts an anti-glioblastoma effect by inducing apoptosis and autophagy and thus can be a potential therapeutic target for glioblastoma therapy.

Exoribonuclease RNase R protects Antarctic Pseudomonas syringae Lz4W from DNA damage and oxidative stress.

Bacterial exoribonucleases play a crucial role in RNA maturation, degradation, quality control, and turnover. In this study, we have uncovered a previously unknown role of 3'-5' exoribonuclease RNase R of Pseudomonas syringae Lz4W in DNA damage and oxidative stress response. Here, we show that neither the exoribonuclease function of RNase R nor its association with the RNA degradosome complex is essential for this function. Interestingly, in P. syringae Lz4W, hydrolytic RNase R exhibits physiological roles similar to phosphorolytic 3'-5' exoribonuclease PNPase of E. coli. Our data suggest that during the course of evolution, mesophilic E. coli and psychrotrophic P. syringae have apparently swapped these exoribonucleases to adapt to their respective environmental growth conditions.

Retracing the evolution of a modern periplasmic binding protein.

Investigating the evolution of structural features in modern multidomain proteins helps to understand their immense diversity and functional versatility. The class of periplasmic binding proteins (PBPs) offers an opportunity to interrogate one of the main processes driving diversification: the duplication and fusion of protein sequences to generate new architectures. The symmetry of their two-lobed topology, their mechanism of binding, and the organization of their operon structure led to the hypothesis that PBPs arose through a duplication and fusion event of a single common ancestor. To investigate this claim, we set out to reverse the evolutionary process and recreate the structural equivalent of a single-lobed progenitor using ribose-binding protein (RBP) as our model. We found that this modern PBP can be deconstructed into its lobes, producing two proteins that represent possible progenitor halves. The isolated halves of RBP are well folded and monomeric proteins, albeit with a lower thermostability, and do not retain the original binding function. However, the two entities readily form a heterodimer in vitro and in-cell. The x-ray structure of the heterodimer closely resembles the parental protein. Moreover, the binding function is fully regained upon formation of the heterodimer with a ligand affinity similar to that observed in the modern RBP. This highlights how a duplication event could have given rise to a stable and functional PBP-like fold and provides insights into how more complex functional structures can evolve from simpler molecular components.

Accurate prediction of protein folding mechanisms by simple structure-based statistical mechanical models

Recent breakthroughs in highly accurate protein structure prediction using deep neural networks have made considerable progress in solving the structure prediction component of the ‘protein folding problem’. However, predicting detailed mechanisms of how proteins fold into specific native structures remains challenging, especially for multidomain proteins constituting most of the proteomes. Here, we develop a simple structure-based statistical mechanical model that introduces nonlocal interactions driving the folding of multidomain proteins. Our model successfully predicts protein folding processes consistent with experiments, without the limitations of protein size and shape. Furthermore, slight modifications of the model allow prediction of disulfide-oxidative and disulfide-intact protein folding. These predictions depict details of the folding processes beyond reproducing experimental results and provide a rationale for the folding mechanisms. Thus, our physics-based models enable accurate prediction of protein folding mechanisms with low computational complexity, paving the way for solving the folding process component of the ‘protein folding problem’.

Comparative analysis of permanent and transient domain-domain interactions in multi-domain proteins.

Protein domains are structural, functional, and evolutionary units. These domains bring out the diversity of functionality by means of interactions with other co-existing domains and provide stability. Hence, it is important to study intra-protein inter-domain interactions from the perspective of types of interactions. Domains within a chain could interact over short timeframes or permanently, rather like protein-protein interactions (PPIs). However, no systematic study has been carried out between two classes, namely permanent and transient domain-domain interactions. In this work, we studied 263 two-domain proteins, belonging to either of these classes and their interfaces on the basis of several factors, such as interface area and details of interactions (number, strength, and types of interactions). We also characterized them based on residue conservation at the interface, correlation of residue motions across domains, its involvement in repeat formation, and their involvement in particular molecular processes. Finally, we could analyze the interactions arising from domains in two-domain monomeric proteins, and we observed significant differences between these two classes of domain interactions and a few similarities. This study will help to obtain a better understanding of structure-function and folding principles of multi-domain proteins.

Unveiling Metastable Ensembles of GRB2 and the Relevance of Interdomain Communication during Folding.

The folding process of multidomain proteins is a highly intricate phenomenon involving the assembly of distinct domains into a functional three-dimensional structure. During this process, each domain may fold independently while interacting with others. The folding of multidomain proteins can be influenced by various factors, including their composition, the structure of each domain, or the presence of disordered regions, as well as the surrounding environment. Misfolding of multidomain proteins can lead to the formation of nonfunctional structures associated with a range of diseases, including cancers or neurodegenerative disorders. Understanding this process is an important step for many biophysical analyses such as stability, interaction, malfunctioning, and rational drug design. One such multidomain protein is growth factor receptor-bound protein 2 (GRB2), an adaptor protein that is essential in regulating cell survival. GRB2 consists of one central Src homology 2 (SH2) domain flanked by two Src homology 3 (SH3) domains. The SH2 domain interacts with phosphotyrosine regions in other proteins, while the SH3 domains recognize proline-rich regions on protein partners during cell signaling. Here, we combined computational and experimental techniques to investigate the folding process of GRB2. Through computational simulations, we sampled the conformational space and mapped the mechanisms involved by the free energy profiles, which may indicate possible intermediate states. From the molecular dynamics trajectories, we used the energy landscape visualization method (ELViM), which allowed us to visualize a three-dimensional (3D) representation of the overall energy surface. We identified two possible parallel folding routes that cannot be seen in a one-dimensional analysis, with one occurring more frequently during folding. Supporting these results, we used differential scanning calorimetry (DSC) and fluorescence spectroscopy techniques to confirm these intermediate states in vitro. Finally, we analyzed the deletion of domains to compare our model outputs to previously published results, supporting the presence of interdomain modulation. Overall, our study highlights the significance of interdomain communication within the GRB2 protein and its impact on the formation, stability, and structural plasticity of the protein, which are crucial for its interaction with other proteins in key signaling pathways.

Structure of full-length cobalamin-dependent methionine synthase and cofactor loading captured in crystallo

Cobalamin-dependent methionine synthase (MS) is a key enzyme in methionine and folate one-carbon metabolism. MS is a large multi-domain protein capable of binding and activating three substrates: homocysteine, folate, and S-adenosylmethionine for methylation. Achieving three chemically distinct methylations necessitates significant domain rearrangements to facilitate substrate access to the cobalamin cofactor at the right time. The distinct conformations required for each reaction have eluded structural characterization as its inherently dynamic nature renders structural studies difficult. Here, we use a thermophilic MS homolog (tMS) as a functional MS model. Its exceptional stability enabled characterization of MS in the absence of cobalamin, marking the only studies of a cobalamin-binding protein in its apoenzyme state. More importantly, we report the high-resolution full-length MS structure, ending a multi-decade quest. We also capture cobalamin loading in crystallo, providing structural insights into holoenzyme formation. Our work paves the way for unraveling how MS orchestrates large-scale domain rearrangements crucial for achieving challenging chemistries.

Computational designing of the ligands of Protein L affinity chromatography based on molecular docking and molecular dynamics simulations

Protein L is a multidomain protein from Peptostreptococcus magnus with binding affinity to kappa light chain of human immunoglobulin (Ig) which is used for the purification of antibody fragments by affinity chromatography. The advances in protein engineering and computational biology approaches lead to the development of engineered affinity ligands with improved properties including binding affinity. In this study, molecular dynamics simulations (MDs) and Osprey software were used to design single B domains of the Protein L with higher affinity to antibody fragments. The modified B domains were then polymerized to ligand with six B domains by homology modeling methods. The results showed that single B domain mutants of MB1 (Thr865Trp) and MB2 (Thr847Met-Thr865Trp) had higher binding affinity to Fab compared to the wild single B domain. Also, MDs and molecular docking results showed that the polymerized Proteins L including the wild and mutated six B domains (6B0, 6B1, and 6B2) were stable during MDs and the two mutants of 6B1 and 6B2 showed higher binding affinity to Fab relative to the wild type. Communicated by Ramaswamy H. Sarma

THU393 An Internal NPT2A Cryptic PDZ Motif Is A Novel Regulator of PTH- And FGF23-sensitive Phosphate Transport

Abstract Disclosure: B.W. Sneddon: None. P.A. Friedman: None. T. Mamonova: None. The Na+-dependent NPT2A (SLC34A1) phosphate co-transporter is a major regulator of phosphate homeostasis. Its most notable structural element is a carboxy-terminal PDZ ligand that binds Na+/H+ Exchanger Regulatory Factor-1 (NHERF1, SLC9A3R1), a multidomain PDZ protein that determines apical NPT2A localization and is required for phosphate transport. NPT2A also possesses an uncharacterized internal PDZ ligand. Two recent clinical reports describe congenital hypophosphatemia in children harboring Arg495His or Arg495Cys variants within the internal PDZ motif. Wild-type internal 494TRL496 PDZ binds NHERF1 PDZ2, which we consider a regulatory domain. Ablating the internal PDZ ligand with a 494AAA496 variant blocked hormone-sensitive phosphate transport. Complementary approaches, including CRISPR/Cas9 technology, side-directed mutagenesis, confocal microscopy, and modeling, showed that NPT2A Arg495 His or Arg495Cys variants do not support PTH or FGF23 action on phosphate transport. CoIP experiments indicate that both variants bound NHERF1 similarly to WT NPT2A. However, in contrast to WT NPT2A, NPT2A Arg495His or Arg495Cys variants remain at the apical membrane and are not internalized in response to PTH. We predict that His or Cys substitution of the charged Arg495 changes the electrostatics, preventing phosphorylation of the upstream Thr494, interfering with phosphate uptake in response to hormone action, and inhibiting NPT2A trafficking. We advance a model where the carboxy-terminal PDZ ligand defines apical localization of NPT2A, while the internal PDZ ligand is essential for hormone-triggered phosphate transport. Presentation: Thursday, June 15, 2023

Droplet digital PCR: A comprehensive tool for genetic analysis and prediction of bispecific antibody assembly during cell line development

Molecular biological methods have emerged as inevitable tools to accompany the process of cell line development for the generation of stable and highly productive manufacturing cell lines in the biopharmaceutical industry. PCR-based methods are especially useful for screening and characterization of cell lines due to their low cost, scalability, precision and propensity for multidimensional read-outs. In this study, the diverse applications of droplet digital PCR (ddPCR) as a molecular biological tool for cell line development are demonstrated. Specifically, it is shown that ddPCR can be used to enable precise, sensitive and reproducible absolute quantification of genomically integrated transgene copies during cell line development and cell bank characterization. Additionally, an amplitude multiplexing approach is applied to simultaneously run multiple assays on different genetic targets in a single reaction and advance clonal screening by measuring gene expression profiles to predict the assembly and homogeneity of difficult-to-express (DTE) proteins. The implementation of ddPCR-based assays during cell line development allows for early screening at a transcriptional level, particularly for complex, multidomain proteins, where balanced polypeptide chain ratios are of primary importance. Moreover, it is demonstrated that ddPCR-based genomic characterization improves the robustness, efficiency and comparability of absolute transgene copy number quantification, an essential genetic parameter that must be demonstrated to regulatory authorities during clinical trial and market authorization application submissions to support genetic stability and consistency of the selected cell substrate.

E2EDA: Protein Domain Assembly Based on End-to-End Deep Learning.

With the development of deep learning, almost all single-domain proteins can be predicted at experimental resolution. However, the structure prediction of multi-domain proteins remains a challenge. Achieving end-to-end protein domain assembly and further improving the accuracy of the full-chain modeling by accurately predicting inter-domain orientation while improving the assembly efficiency will provide significant insights into structure-based drug discovery. In this work, we propose an End-to-End Domain Assembly method based on deep learning, named E2EDA. We first develop RMNet, an EfficientNetV2-based deep learning model that fuses multiple features using an attention mechanism to predict inter-domain rigid motion. Then, the predicted rigid motions are transformed into inter-domain spatial transformations to directly assemble the full-chain model. Finally, the scoring strategy RMscore is designed to select the best model from multiple assembled models. The experimental results show that the average TM-score of the model assembled by E2EDA on the benchmark set (282) is 0.827, which is better than those of other domain assembly methods SADA (0.792) and DEMO (0.730). Meanwhile, on our constructed multi-domain data set from AlphaFold DB, the model reassembled by E2EDA is 7.0% higher in TM-score compared to the full-chain model predicted by AlphaFold2, indicating that E2EDA can capture more accurate inter-domain orientations to improve the quality of the model predicted by AlphaFold2. Furthermore, compared to SADA and AlphaFold2, E2EDA reduced the average runtime on the benchmark by 64.7% and 19.2%, respectively, indicating that E2EDA can significantly improve assembly efficiency through an end-to-end approach. The online server is available at http://zhanglab-bioinf.com/E2EDA.

Immunohistochemical expression of paxillin in ameloblastoma and odontogenic keratocyst: An observational study.

Cell adhesion molecules (CAMs) are found on the surface of all cells, where they allow dynamic processes to take place. These include cadherins, integrins, selectins and Immunoglobulin superfamily. Directly associated with β-integrin tails is a multidomain protein known as paxillin. However, CAMs participate in cell-cell and extracellular matrix-cell interactions during histomorphogenesis in the various phases of odontogenesis. Some tumours or cysts like ameloblastoma (AB) or odontogenic keratocyst (OKC) having odontogenic origin show disturbance in the interaction of these CAMs. Hence, the assessment of paxillin expression in AB and OKC was carried out. The present observational study comprised 30 clinically and histologically confirmed cases of AB and OKC. All the slides were stained immunohistochemically using a paxillin antibody. Upon comparison of staining intensity of paxillin among AB and OKC showed statistically significant result, whereas quantitative staining and final summation showed non-significant result. Gender-wise comparison of paxillin staining intensity, quantitative staining and final summation among OKC showed significant result; however, in AB, staining intensity showed non-significant result, whereas quantitative staining and final summation showed significant result. Paxillin has the greatest influence on tissue morphogenesis and development. The regulation of cell mobility is aided by the multiple roles that paxillin plays in a range of cells and tissues. However, further studies using a large sample size, along with other molecular analytical methods, may be essential to draw a definite conclusion about the association of paxillin and its exact function in OKC and AB.

Molecular Mechanism of STIL Coiled-Coil Domain Oligomerization.

Coiled-coil domains (CCDs) play key roles in regulating both healthy cellular processes and the pathogenesis of various diseases by controlling protein self-association and protein-protein interactions. Here, we probe the mechanism of oligomerization of a peptide representing the CCD of the STIL protein, a tetrameric multi-domain protein that is over-expressed in several cancers and associated with metastatic spread. STIL tetramerization is mediated both by an intrinsically disordered domain (STIL400-700) and a structured CCD (STIL CCD718-749). Disrupting STIL oligomerization via the CCD inhibits its activity in vivo. We describe a comprehensive biophysical and structural characterization of the concentration-dependent oligomerization of STIL CCD peptide. We combine analytical ultracentrifugation, fluorescence and circular dichroism spectroscopy to probe the STIL CCD peptide assembly in solution and determine dissociation constants of both the dimerization, (KD = 8 ± 2 µM) and tetramerization (KD = 68 ± 2 µM) of the WT STIL CCD peptide. The higher-order oligomers result in increased thermal stability and cooperativity of association. We suggest that this complex oligomerization mechanism regulates the activated levels of STIL in the cell and during centriole duplication. In addition, we present X-ray crystal structures for the CCD containing destabilising (L736E) and stabilising (Q729L) mutations, which reveal dimeric and tetrameric antiparallel coiled-coil structures, respectively. Overall, this study offers a basis for understanding the structural molecular biology of the STIL protein, and how it might be targeted to discover anti-cancer reagents.

The Dictyostelium discoideum FimA protein, unlike yeast and plant fimbrins, is regulated by calcium similar to mammalian plastins

Plastins, also known as fimbrins, are highly conserved eukaryotic multidomain proteins that are involved in actin-bundling. They all contain four independently folded Calponin Homology-domains and an N-terminal headpiece that is comprised of two calcium-binding EF-hand motifs. Since calcium-binding has been shown to be integral to regulating the activity of the three mammalian plastin proteins, we decided to study the properties of the headpiece regions of fimbrins from the model plant Arabidopsis thaliana, the yeasts Saccharomyces cerevisiae and Schizosaccharomyces pombe and the amoeba Dictyostelium discoideum. Of these protein domains only the FimA headpiece from the amoeba protein possesses calcium binding properties. Structural characterization of this protein domain by multidimensional NMR and site-directed mutagenesis studies indicates that this EF-hand region of FimA also contains a regulatory ‘switch helix’ that is essential to regulating the activity of the human L-plastin protein. Interestingly this regulatory helical region seems to be lacking in the plant and yeast proteins and in fimbrins from all other nonmotile systems. Typical calmodulin antagonists can displace the switch-helix from the FimA headpiece, suggesting that such drugs can deregulate the Ca2+-regulation of the actin-bunding in the amoeba, thereby making it a useful organism for drug screening against mammalian plastins.

Time-resolved study on signaling pathway of photoactivated adenylate cyclase and its nonlinear optical response

Photoactivated adenylate cyclases (PACs) are multidomain BLUF proteins that regulate the cellular levels of cAMP in a light-dependent manner. The signaling route and dynamics of PAC from Oscillatoria acuminata (OaPAC), which consists of a light sensor BLUF domain, an adenylate cyclase domain, and a connector helix (α3-helix), were studied by detecting conformational changes in the protein moiety. Although circular dichroism and small-angle X-ray scattering measurements did not show significant changes upon light illumination, the transient grating method successfully detected light-induced changes in the diffusion coefficient (diffusion-sensitive conformational change (DSCC)) of full-length OaPAC and the BLUF domain with the α3-helix. DSCC of full-length OaPAC was observed only when both protomers in a dimer were photoconverted. This light intensity dependence suggests that OaPAC is a cyclase with a nonlinear light intensity response. The enzymatic activity indeed nonlinearly depends on light intensity, that is, OaPAC is activated under strong light conditions. It was also found that both DSCC and enzymatic activity were suppressed by a mutation in the W90 residue, indicating the importance of the highly conserved Trp in many BLUF domains for the function. Based on these findings, a reaction scheme was proposed together with the reaction dynamics.

An integrated approach of NMR experiments and MD simulations visualizes structural dynamics of a cyclic multi-domain protein.

Cyclization can stabilize the structure of proteins, as previously demonstrated in single-domain proteins. Although Lys48-linked polyubiquitin, a multi-domain protein, is also known to be cyclized in human cells, the structural effects of cyclization remain unclear. Here, we examined the impact of cyclization on the structural stability and dynamics of cyclic Lys48-linked diubiquitin (Ub2 ). As expected, cyclization increased the thermal stability of Ub2 and its resistance to proteolytic digestion, indicating that cyclization stabilized the structure of Ub2 . Furthermore, cyclization repressed the interdomain motion in Ub2 , but cyclic Ub2 still exhibited microsecond conformational exchange in NMR relaxation dispersion experiments. A series of long coarse-grained (CG) MD simulations visualized how cyclization slowed down the intrinsic nanosecond open-closed domain motion of Ub2 to microseconds. Thus, CG-MD analysis helped to explain the unexpected NMR relaxation results, thereby facilitating characterization of the structural stabilization of cyclic Ub2 .

Structure/function studies of the NAD+-dependent DNA ligase from the poly-extremophile Deinococcus radiodurans reveal importance of the BRCT domain for DNA binding

Bacterial NAD+-dependent DNA ligases (LigAs) are enzymes involved in replication, recombination, and DNA-repair processes by catalyzing the formation of phosphodiester bonds in the backbone of DNA. These multidomain proteins exhibit four modular domains, that are highly conserved across species, with the BRCT (breast cancer type 1 C-terminus) domain on the C-terminus of the enzyme. In this study, we expressed and purified both recombinant full-length and a C-terminally truncated LigA from Deinococcus radiodurans (DrLigA and DrLigA∆BRCT) and characterized them using biochemical and X-ray crystallography techniques. Using seeds of DrLigA spherulites, we obtained ≤ 100 µm plate crystals of DrLigA∆BRCT. The crystal structure of the truncated protein was obtained at 3.4 Å resolution, revealing DrLigA∆BRCT in a non-adenylated state. Using molecular beacon-based activity assays, we demonstrated that DNA ligation via nick sealing remains unaffected in the truncated DrLigA∆BRCT. However, DNA-binding assays revealed a reduction in the affinity of DrLigA∆BRCT for dsDNA. Thus, we conclude that the flexible BRCT domain, while not critical for DNA nick-joining, plays a role in the DNA binding process, which may be a conserved function of the BRCT domain in LigA-type DNA ligases.

Klebsiella phage KP34gp57 capsular depolymerase structure and function: from a serendipitous finding to the design of active mini-enzymes against K. pneumoniae.

In this work, we determined the structure of Klebsiella phage KP34p57 capsular depolymerase and dissected the role of individual domains in trimerization and functional activity. The crystal structure serendipitously revealed that the enzyme can exist in a monomeric state once deprived of its C-terminal domain. Based on the crystal structure and site-directed mutagenesis, we localized the key catalytic residues in an intra-subunit deep groove. Consistently, we show that C-terminally trimmed KP34p57 variants are monomeric, stable, and fully active. The elaboration of monomeric, fully active phage depolymerases is innovative in the field, as no previous example exists. Indeed, mini phage depolymerases can be combined in chimeric enzymes to extend their activity ranges, allowing their use against multiple serotypes.

Targeting SARS-CoV-2 Macrodomain-1 to Restore the Innate Immune Response Using In Silico Screening of Medicinal Compounds and Free Energy Calculation Approaches.

Among the different drug targets of SARS-CoV-2, a multi-domain protein known as NSP3 is a critical element of the translational and replication machinery. The macrodomain-I, in particular, has been reported to have an essential role in the viral attack on the innate immune response. In this study, we explore natural medicinal compounds and identify potential inhibitors to target the SARS-CoV-2-NSP3 macrodomain-I. Computational modeling and simulation tools were utilized to investigate the structural-dynamic properties using triplicates of 100 ns MD simulations. In addition, the MM/GBSA method was used to calculate the total binding free energy of each inhibitor bound to macrodomain-I. Two significant hits were identified: 3,5,7,4'-tetrahydroxyflavanone 3'-(4-hydroxybenzoic acid) and 2-hydroxy-3-O-beta-glucopyranosyl-benzoic acid. The structural-dynamic investigation of both compounds with macrodomain-I revealed stable dynamics and compact behavior. In addition, the total binding free energy for each complex demonstrated a robust binding affinity, of ΔG -61.98 ± 0.9 kcal/mol for Compound A, while for Compound B, the ΔG was -45.125 ± 2.8 kcal/mol, indicating the inhibitory potential of these compounds. In silico bioactivity and dissociation constant (KD) determination for both complexes further validated the inhibitory potency of each compound. In conclusion, the aforementioned natural products have the potential to inhibit NSP3, to directly rescue the host immune response. The current study provides the basis for novel drug development against SARS-CoV-2 and its variants.