In silico prediction, characterization, docking studies and molecular dynamics simulation of human p97 in complex with p37 cofactor

Insights into the interactions of RWP-RK and their targets: Role of serine and its conservation across species

The RWP-RK domain is a key DNA-binding domain found in all NIN (Nodule Inception)/NLP (NIN-like proteins) and RKD (RWP-RK Domain Containing) transcription factors (TFs). The RWP-RK domain in NINs/NLPs contains a highly evolutionarily conserved sequence, RWPSRK, while in RKDs, the fourth serine (S) amino acid is substituted with either tyrosine (Y) or histidine (H). To regulate autoregulation of nodulation, the RWP-RK domain of NIN TF binds to the promoter region of CLE peptides but not RKDs. Therefore, investigating the protein-DNA interaction from a structural perspective is essential to understand the evolutionary significance of the serine (S) residue of the RWP-RK domain. Herein, we have modelled both the wild type (WT) and the variant RWP-RK domains containing substitutions like glutamic acid (E), tyrosine (Y), and histidine (H) and docked them with the modelled pCLE13 cis-element. Our docking results revealed that a helix-turn-helix (HTH) motif of the RWP-RK domain interacts with pCLE13. The WT HTH-DNA complex exhibited the most negative binding free energy, indicating a strong interaction, particularly hydrogen bonds acting between them. Simulation analysis of WT and variant models provided deeper insights into protein-DNA binding dynamics. The hydrogen bond occupancy percentage indicated that the fourth serine (S) residue is vital for maintaining a significant percentage of hydrogen bonds with DNA. The variants substituting this conserved serine (S) residue displayed energetic frustration upon binding to DNA and lost correlation among their residues. Overall, it suggested that serine (S) residue of the RWP-RK domain of all NINs/NLPs is crucial for appropriate protein-DNA interaction, which might be required for their biological relevance.

Screening of phytochemicals from Clerodendrum inerme (L.) Gaertn as potential anti-breast cancer compounds targeting EGFR: an in-silico approach

Breast cancer (BC) is the most prevalent malignancy among women around the world. The epidermal growth factor receptor (EGFR) is a tyrosine kinase receptor (RTK) of the ErbB/HER family. It is essential for triggering the cellular signaling cascades that control cell growth and survival. However, perturbations in EGFR signaling lead to cancer development and progression. Hence, EGFR is regarded as a prominent therapeutic target for breast cancer. Therefore, in the current investigation, EGFR was targeted with phytochemicals from Clerodendrum inerme (L.) Gaertn (C. inerme). A total of 121 phytochemicals identified by gas chromatography-mass spectrometry (GC-MS) analysis were screened against EGFR through molecular docking, ADMET analyses (Absorption, Distribution, Metabolism, Excretion, and Toxicity), PASS predictions, and molecular dynamics simulation, which revealed three potential hit compounds with CIDs 10586 [i.e. alpha-bisabolol (-6.4 kcal/mol)], 550281 [i.e. 2,(4,4-Trimethyl-3-hydroxymethyl-5a-(3-methyl-but-2-enyl)-cyclohexene) (-6.5 kcal/mol)], and 161271 [i.e. salvigenin (-7.4 kcal/mol)]. The FDA-approved drug gefitinib was used to compare the inhibitory effects of the phytochemicals. The top selected compounds exhibited good ADMET properties and obeyed Lipinski's rule of five (ROF). The molecular docking analysis showed that salvigenin was the best among the three compounds and formed bonds with the key residue Met 793. Furthermore, the molecular mechanics generalized born surface area (MMGBSA) calculations, molecular dynamics simulation, and normal mode analysis validated the binding affinity of the compounds and also revealed the strong stability and compactness of phytochemicals at the docked site. Additionally, DFT and DOS analyses were done to study the reactivity of the compounds and to further validate the selected phytochemicals. These results suggest that the identified phytochemicals possess high inhibitory potential against the target EGFR and can treat breast cancer. However, further in vitro and in vivo investigations are warranted towards the development of these constituents into novel anti-cancer drugs.

Amino acids characterization based on frequency and interaction analysis in human antigen-antibody complexes from Thera-SAbDab

BackgroundAntibodies are composed of light and heavy chains, both of which have constant and variable regions. The diversity, specific binding ability and therapeutic potential of antibodies are determined by hypervariable loops called complementarity-determining regions (CDRs), with the other regions being the framework regions.ObjectiveTo investigate the key amino acid patterns in various antibody regions in the human therapeutic antigen-antibody (Ag-Ab) complexes collected from the Thera-SAbDab database.MethodThe study focuses on identifying the amino acid frequency, diversity index in CDRs, paratope-epitope amino acid interactions, amino acid bond formation frequency, and bond types among selected therapeutic Ag-Ab complexes.ResultsThe results revealed that Ser is highly distributed in the overall light chain CDRs while Gly is highly distributed in the heavy chain CDRs. CDR profiling analysis indicated that the average amino acid diversity in heavy chain CDRs is 60% to 70%, while in the light chain, it is 50% to 60%. Aromatic residues such as Tyr, Trp and Phe are the top contributors to these paratope-epitope interactions in the light and heavy chains. Moreover, we examined the frequency of amino acids in light and heavy chains of Ag-Ab complexes. Importantly, the outcome of this study leverages the in depth analysis on single residues, dipeptides, and tripeptides for the therapeutic Ag-Ab complexes.ConclusionWe conclude that the amino acid frequency and interaction analysis centered on therapeutic Ag-Ab complexes will benefit antibody engineering parameters such as antibody design, optimization, affinity maturation, and overall antibody development.

In vitro and in silico study of the synergistic anticancer effect of alpinumisoflavone with gemcitabine on pancreatic ductal adenocarcinoma through suppression of ribonucleotide reductase subunit-M1

A highly aggressive neoplastic disease, pancreatic ductal adenocarcinoma (PDAC) is documented as the third chief cause of cancer-associated mortality in both sexes combined in the United States. For decades, gemcitabine-based chemotherapy has been embraced as a cornerstone drug for the treatment of PDAC. However, there have been several unsolved problems, including cytotoxicity, and chemoresistance. Gemcitabine efficacy was attributed to the attenuation of ribonucleotide reductase subunit-M1 (RRM1). Overexpression of RRM1 in PDAC is highly correlated with gemcitabine resistance and reduced gemcitabine sensitivity, resulting in a poor survival rate even after gemcitabine treatment. Moreover, the status of TP53, a tumor suppressor gene, assumes a decisive role in the response of PDAC to gemcitabine. Therefore, targeting RRM1 and P53 might be a therapeutic strategy for strengthening gemcitabine efficacy and cytotoxicity against PDAC. Alpinumisoflavone (AIF) is a prenylated isoflavone originated in Cudrania tricuspidate with versatile bioactive properties, including anticancer activity. However, there was no report whether AIF can exert anticancer effect and exhibit synergistic effect with gemcitabine against PDAC. Therefore, the anticancer properties of AIF were assessed with PANC-1 and MIA PaCa-2. In addition, synergism between AIF and gemcitabine were analyzed. Moreover, the contribution of P53 and RRM1 expression to gemcitabine resistance was assessed by comparing their protein levels in PDAC cells and normal pancreatic cells. The interactions of AIF with RRM1 protein were confirmed by molecular docking and dynamics simulation. Therefore, AIF enhances gemcitabine efficacy against PDAC through the regulation of P53 and RRM1.

Integrating protein interaction and pathway crosstalk network reveals a promising therapeutic approach for psoriasis through apoptosis induction

Psoriasis is a complex inflammatory skin disease manifested by altered proliferation and differentiation of keratinocytes with dysfunctional apoptosis. This study aimed to identify regulatory factors and comprehend the underlying mechanisms of inefficient apoptosis to open up promising therapeutic approaches. Incorporating human protein interactions, apoptosis proteins, and physical relationships of psoriasis-apoptosis proteins helped us to generate a psoriasis-apoptosis interaction (SAI) network. Subsequently, topological and functional analyses of the SAI network revealed effective proteins, functional modules, hub motifs, dysregulated pathways and transcriptional gene regulatory factors. Network pharmacology, molecular docking and molecular dynamics simulation methods identified the potential drug-target interactions. RELA, MAPK1, MAPK3, MMP9, IL1B, AKT1 and STAT1 were revealed as effective proteins. The MAPK1-MAPK3-RELA motif was identified as a hub regulator in the crosstalk between 41 pathways. Among all pathways, "lipid and atherosclerosis" was found to be the predominant pathway. Acetylcysteine, arsenic-trioxide, β-elemene, bortezomib and curcumin were identified as potential drugs to inhibit pathway crosstalk. Experimental verifications were performed using the literature search, GSE13355 and GSE14905 microarray datasets. Drug-protein-pathway interactions associated with apoptosis were deciphered. These findings highlight the role of hub motif-mediated pathway-pathway crosstalk associated with apoptosis in the complexity of psoriasis and suggest crosstalk inhibition as an effective therapeutic approach.

Cross-talk between Vimentin and autophagy regulates blood-testis barrier disruption induced by cadmium

The blood-testis barrier (BTB) plays a vital role in mammalian spermatogenesis by separating the seminiferous epithelium into an adluminal and a basal compartment. Cadmium (Cd) is a toxic heavy metal that is widely present in the environment. We observed that Cd can induce BTB disruption, leading to apoptosis of testicular cells. However, the molecular mechanisms contributing to BTB injury induced by Cd have not yet been fully clarified. Vimentin (Vim) is an important desmosome-like junction protein that mediates robust adhesion in the BTB. In this study, we investigated how Vim responds to Cd. We found that Cd treatment led to a significant decrease in Vim expression, accompanied by a marked increase in LC3-II expression and a higer number of autophagosomes. Interestingly, we also observed that Cd-induced autophagy was associated with decreased Vim activity and enhanced apoptosis of testicular cells. To further investigate the role of autophagy in Vim regulation under Cd exposure, we treated cells with an autophagy inhibitor called 3-MA. We found that 3-MA treatment enhanced Vim expression and improved the disruption of the BTB under Cd exposure. Additionally, the inhibition of Vim confirmed the role of autophagy in modulating Vim expression. These results reveal a previously unknown regulatory mechanism of Cd involving the interplay between a heavy metal and a protein.

Molecular docking and dynamics simulation studies uncover the host-pathogen protein-protein interactions in Penaeus vannamei and Vibrio parahaemolyticus

Shrimp aquaculture contributes significantly to global economic growth, and the whiteleg shrimp, Penaeus vannamei, is a leading species in this industry. However, Vibrio parahaemolyticus infection poses a major challenge in ensuring the success of P. vannamei aquaculture. Despite its significance in this industry, the biological knowledge of its pathogenesis remains unclear. Hence, this study was conducted to identify the interaction sites and binding affinity between several immune-related proteins of P. vannamei with V. parahaemolyticus proteins associated with virulence factors. Potential interaction sites and the binding affinity between host and pathogen proteins were identified using molecular docking and dynamics (MD) simulation. The P. vannamei-V. parahaemolyticus protein-protein interaction of Complex 1 (Ferritin-HrpE/YscL family type III secretion apparatus protein), Complex 2 (Protein kinase domain-containing protein-Chemotaxis CheY protein), and Complex 3 (GPCR-Chemotaxis CheY protein) was found to interact with -4319.76, -5271.39, and -4725.57 of the docked score and the formation of intermolecular bonds at several interacting residues. The docked scores of Complex 1, Complex 2, and Complex 3 were validated using MD simulation analysis, which revealed these complexes greatly contribute to the interactions between P. vannamei and V. parahaemolyticus proteins, with binding free energies of -22.50 kJ/mol, -30.20 kJ/mol, and -26.27 kJ/mol, respectively. This finding illustrates the capability of computational approaches to search for molecular binding sites between host and pathogen, which could increase the knowledge of Vibrio spp. infection on shrimps, which then can be used to assist in the development of effective treatment.

Structure-based virtual screening, molecular simulation and free energy calculations of traditional Chinese medicine, ZINC database revealed potent inhibitors of estrogen-receptor α (ERα)

Breast Cancer, a heterogeneous disease at the molecular level, is the most common cause of woman mortality worldwide. We used molecular screening and simulation approaches to target nuclear receptor protein-estrogen receptor alpha (Erα) protein to design and develop of specific and compelling drugs from traditional Chinese medicine (TCM), and ZINC database against pathophysiology of breast cancer. Using virtual screening, only six hits TCM22717, TCM23524, TCM31953, while ZINC05632920, ZINC05773243, and ZINC12780336 demonstrated better pharmacological potential than the 4-hydroxytamoxifen (OHT) taken as control. Binding mode of each of the top hit revealed that these compounds could block the main active site residues and block the function of Erα protein. Moreover, molecular simulation revealed that the identified compounds exhibit stable dynamics and may induce stronger therapeutic effects in experimental setup. All the complexes reported tighter structural packing and less flexible behaviour. We found that the average hydrogen bonds in the identified complexes remained higher than the control drug. Finally, the total binding free energy demonstrated the best hits among the all. The BF energy results revealed -30.4525 ± 3.3565 for the 4-hydroxytamoxifen (OHT)/Erα complex, for the TCM22717/Erα -57.0597 ± 3.4852 kcal/mol, for the TCM23524/Erα complex the BF energy was -56.9084 ± 3.3737 kcal/mol, for the TCM31953/Erα the BF energy was -32.4191 ± 3.8864 kcal/mol while for the ZINC05632920/Erα complex -46.3182 ± 2.7380, ZINC05773243/Erα complex -38.3690 ± 2.8240, and ZINC12780336/Erα complex the BF energy was calculated to be -35.8048 ± 4.1571 kcal/mol.Communicated by Ramaswamy H. Sarma.