Performance and Its Limits in Rigid Body Protein-Protein Docking

Design and evaluation of a multi-epitope DNA vaccine against HPV16

Cervical cancer, among the deadliest cancers affecting women globally, primarily arises from persistent infection with high-risk human papillomavirus (HPV). To effectively combat persistent infection and prevent the progression of precancerous lesions into malignancy, a therapeutic HPV vaccine is under development. This study utilized an immunoinformatics approach to predict epitopes of cytotoxic T lymphocytes (CTLs) and helper T lymphocytes (HTLs) using the E6 and E7 oncoproteins of the HPV16 strain as target antigens. Subsequently, through meticulous selection of T-cell epitopes and other necessary elements, a multi-epitope vaccine was constructed, exhibiting good immunogenic, physicochemical, and structural characteristics. Furthermore, in silico simulations showed that the vaccine not only interacted well with toll-like receptors (TLR2/TLR3/TLR4), but also induced a strong innate and adaptive immune response characterized by elevated Th1-type cytokines, such as interferon-gamma (IFN-γ) and interleukin-2 (IL2). Additionally, our study investigated the effects of different immunization intervals on immune responses, aiming to optimize a time-efficient immunization program. In animal model experiments, the vaccine exhibited robust immunogenic, therapeutic, and prophylactic effects. Administered thrice, it consistently induced the expansion of specific CD4 and CD8 T cells, resulting in substantial cytokines release and increased proliferation of memory T cell subsets in splenic cells. Overall, our findings support the potential of this multi-epitope vaccine in combating HPV16 infection and signify its candidacy for future HPV vaccine development.

Super epitope dengue vaccine instigated serotype independent immune protection in-silico

Dengue becomes the most common life-threatening infectious arbovirus disease globally, with prevalence in the tropical and subtropical areas. The major clinical features include dengue haemorrhagic fever (DHF) and dengue shock syndrome (DSS), a condition of hypovolemic shock. Four different serotypes of the dengue virus, known as dengue virus serotype (DENV)- 1, 2, 3 and 4 can infect humans. Only one vaccine is available in the market, named Dengvaxia by Sanofi Pasteur, but there is no desired outcome of this treatment due the antibody dependent enhancement (ADE) of the multiple dengue serotypes. As of now, there is no cure against dengue disease. Our goal in this work was to create a subunit vaccine based on several epitopes that would be effective against every serotype of the dengue virus. Here, computational methods like- immunoinformatics and bioinformatics were implemented to find out possible dominant epitopes. A total of 21 epitopes were chosen using various in-silico techniques from the expected 133 major histocompatibility complex (MHC)- I and major histocompatibility complex (MHC)- II epitopes, along with 95 B-cell epitopes which were greatly conserved. Immune stimulant, non-allergenic and non-toxic immunodominant epitopes (super epitopes) with a suitable adjuvant (Heparin-Binding Hemagglutinin Adhesin, HBHA) were used to construct the vaccine. Following the physicochemical analysis, vaccine construct was docked with Toll-like receptors (TLRs) to predict the immune stimulation. Consequently, the optimal docked complex that demonstrated the least amount of ligand-receptor complex deformability was used to conduct the molecular dynamics analysis. By following the codon optimization, the final vaccine molecule was administered into an expressing vector to perform in-silico cloning. The robust immune responses were generated in the in-silico immune simulation analysis. Hence, this study provides a hope to control the dengue infections. For validation of the immune outcomes, in-vitro as well as in-vivo investigations are essential.

Computational investigation of the inhibitory interaction of IRF3 and SARS-CoV-2 accessory protein ORF3b

ORF3b is one of the SARS-CoV-2 accessory proteins. Previous experimental study suggested that ORF3b prevents IRF3 translocating to nucleus. However, the biophysical mechanism of ORF3b-IRF3 interaction is elusive. Here, we explored the conformation ensemble of ORF3b using all-atom replica exchange molecular dynamics simulation. Disordered ORF3b has mixed α-helix, β-turn and loop conformers. The potential ORF3b-IRF3 binding modes were searched by docking representative ORF3b conformers with IRF3, and 50 ORF3b-IRF3 complex poses were screened using molecular dynamics simulations ranging from 500 to 1000 ns. We found that ORF3b binds IRF3 predominantly on its CBP binding and phosphorylated pLxIS motifs, with CBP binding site has the highest binding affinity. The ORF3b-IRF3 binding residues are highly conserved in SARS-CoV-2. Our results provided biophysics insights into ORF3b-IRF3 interaction and explained its interferon antagonism mechanism.

In silico study to explore the mechanism of Toxoplasma-induced inflammation and target therapy based on sero and salivary Toxoplasma

We aimed to assess salivary and seroprevalence of Toxoplasma immunoglobulins in risky populations and evaluate drug docking targeting TgERP. A cross-sectional study was conducted in Alexandria University hospitals' outpatient clinics. 192 participants were enrolled from September 2022 to November 2023. Anti-Toxoplasma IgG and IgM were determined in serum and saliva by ELISA. An in-Silico study examined TgERP's protein-protein interactions (PPIs) with pro-inflammatory cytokine receptors, anti-inflammatory cytokine, cell cycle progression regulatory proteins, a proliferation marker, and nuclear envelope integrity-related protein Lamin B1. Our findings revealed that anti-T. gondii IgG were detected in serum (66.1%) and saliva (54.7%), with 2.1% of both samples were positive for IgM. Salivary IgG had 75.59% sensitivity, 86.15% specificity, 91.40% PPV, 64.40% NPP, 79.17% accuracy and fair agreement with serum IgG. On the other hand, the sensitivity, specificity, PPV, NPV, and accuracy in detecting salivary IgM were 75.0%, 99.47%, 75.0%, 99.47%, and 98.96%. AUC 0.859 indicates good discriminatory power. Examined synthetic drugs and natural products can target specific amino acids residues of TgERP that lie at the same binding interface with LB1 and Ki67, subsequently, hindering their interaction. Hence, salivary samples can be a promising diagnostic approach. The studied drugs can counteract the pro-inflammatory action of TgERP.

Toll-interacting protein is activated by the transcription factor GATA1 and Sp1 to negatively regulate NF-κB and MAPK pathways in the Japanese eel (Anguilla japonica)

Toll-interacting protein (Tollip) serves as a crucial inhibitory factor in the modulation of Toll-like receptor (TLR)-mediated innate immunological responses. The structure and function of Tollip have been well documented in mammals, yet the information in teleost remained limited. This work employed in vitro overexpression and RNA interference in vivo and in vitro to comprehensively examine the regulatory effects of AjTollip on NF-κB and MAPK signaling pathways. The levels of p65, c-Fos, c-Jun, IL-1, IL-6, and TNF-α were dramatically reduced following overexpression of AjTollip, whereas knocking down AjTollip in vivo and in vitro enhanced those genes' expression. Protein molecular docking simulations showed AjTollip interacts with AjTLR2, AjIRAK4a, and AjIRAK4b. A better understanding of the transcriptional regulation of AjTollip is crucial to elucidating the role of Tollip in fish antibacterial response. Herein, we cloned and characterized a 2.2 kb AjTollip gene promoter sequence. The transcription factors GATA1 and Sp1 were determined to be associated with the activation of AjTollip expression by using promoter truncation and targeted mutagenesis techniques. Collectively, our results indicate that AjTollip suppresses the NF-κB and MAPK signaling pathways, leading to the decreased expression of the downstream inflammatory factors, and GATA1 and Sp1 play a vital role in regulating AjTollip expression.

De novo structure prediction of meteorin and meteorin-like protein for identification of domains, functional receptor binding regions, and their high-risk missense variants

Meteorin (Metrn) and Meteorin-like (Metrnl) are homologous secreted proteins involved in neural development and metabolic regulation. In this study, we have performed de novo structure prediction and analysis of both Metrn and Metrnl using Alphafold2 (AF2) and RoseTTAfold (RF). Based on the domain and structural homology analysis of the predicted structures, we have identified that these proteins are composed of two functional domains, a CUB domain and an NTR domain, connected by a hinge/loop region. We have identified the receptor binding regions of Metrn and Metrnl using the machine-learning tools ScanNet and Masif. These were further validated by docking Metrnl with its reported KIT receptor, thus establishing the role of each domain in the receptor interaction. Also, we have studied the effect of non-synonymous SNPs on the structure and function of these proteins using an array of bioinformatics tools and selected 16 missense variants in Metrn and 10 in Metrnl that can affect the protein stability. This is the first study to comprehensively characterize the functional domains of Metrn and Metrnl at their structural level and identify the functional domains, and protein binding regions. This study also highlights the interaction mechanism of the KIT receptor and Metrnl. The predicted deleterious SNPs will allow further understanding of the role of these variants in modulating the plasma levels of these proteins in disease conditions such as diabetes.Communicated by Ramaswamy H. Sarma.

Exploring malaria parasite surface proteins to devise highly immunogenic multi-epitope subunit vaccine for Plasmodium falciparum

BACKGROUND

Malaria has remained a major health concern for decades among people living in tropical and sub-tropical countries. Plasmodium falciparum is one of the critical species that cause severe malaria and is responsible for major mortality. Moreover, the parasite has generated resistance against all WHO recommended drugs and therapies. Therefore, there is an urgent need for preventive measures in the form of reliable vaccines to achieve the target of a malaria-free world. Surface proteins are the preferable choice for subunit vaccine development because they are rapidly detected and engaged by host immune cells and vaccination-induced antibodies. Additionally, abundant surface or membrane proteins may contribute to the opsonization of pathogens by vaccine-induced antibodies.

RESULTS

In our study, we have listed all those surface proteins from the literature that could be functionally important and essential for infection and immune evasion of the malaria parasite. Eight Plasmodium surface and membrane proteins from the pre-erythrocyte and erythrocyte stages were shortlisted. Thirty-seven epitopes (B-cell, CTL, and HTL epitopes) from these proteins were predicted using immune-informatic tools and joined with suitable peptide linkers to design a vaccine construct. A TLR-4 agonist peptide adjuvant was added at the N-terminus of the multi-epitope series, followed by the PADRE sequence and EAAAK linker. The TLR-4 receptor was docked with the construct's anticipated model structure. The complex of vaccine and TLR-4, with the lowest energy -1514, was found to be stable under simulated physiological settings.

CONCLUSION

This study has provided a novel multi-epitope construct that may be exploited further for the development of an efficient vaccine for malaria.

Identification and characterisation of a novel interaction between oestrogen receptor alpha and FOXP2

Forkhead box P2 (FOXP2) regulates expression of various genes and is associated with language, speech and neural development as well as cancer. Since there may be a putative link between sex and language and because transcription factors rarely function in isolation, this study aims to investigate whether FOXP2 directly associates with oestrogen receptor α (ER1), a nuclear receptor responsible for sexual differentiation that is also associated with cancer. Isothermal titration calorimetry and fluorescence anisotropy were used to investigate the interaction between the DNA-binding forkhead domain (FHD) of FOXP2, the N-terminal region (NT) of FOXP2, and the ligand-binding domain (LBD) of ER1. ER1 LBD does not interact with FOXP2 NT but associates with apo-FOXP2 FHD in an enthalpically favourable manner. The affinity of this interaction is inversely correlated to the salt concentration. Additionally, FOXP2 FHD that is bound to ER1 LBD, has reduced ability to interact with its cognate DNA. This research identifies a novel interaction between ER1 LBD and FOXP2 FHD and shows that the interaction is regulated by salt. Moreover, FOXP2 FHD cannot bind to both ER1 LBD and DNA simultaneously, suggesting that this interaction could be involved in regulating the transcriptional pathway of FOXP2 should the interaction be found in vivo. This study could serve as a foundation for uncovering the basis of sexual dimorphism in speech and language development and related disorders and potentially offers an alternate for targeted cancer therapies.

Association mechanism between Arabidopsis immune coreceptor BAK1 and Pseudomonas syringae effector HopF2

Brassinosteroid activated kinase 1 (BAK1) is a cell-surface coreceptor which plays multiple roles in innate immunity of plants. HopF2 is an effector secreted by the bacterial pathogen Pseudomonas syringae pv. tomato DC3000 into Arabidopsis and suppresses host immune system through interaction with BAK1 as well as its downstream kinase MKK5. The association mechanism of HopF2 to BAK1 remains unclear, which prohibits our understanding and subsequent interfering of their interaction for pathogen management. Herein, we found the kinase domain of BAK1 (BAK1-KD) is sufficient for HopF2 association. With a combination of hydrogen/deuterium exchange mass spectrometry and mutational assays, we found a region of BAK1-KD N-lobe and a region of HopF2 head subdomain are critical for intermolecular interaction, which is also supported by unbiased protein-protein docking with ClusPro and kinase activity assay. Collectively, this research presents the interaction mechanism between Arabidopsis BAK1 and P. syringae HopF2, which could pave the way for bactericide development that blocking the functioning of HopF2 toward BAK1.

Targeting ARNT attenuates chemoresistance through destabilizing p38α-MAPK signaling in glioblastoma

Glioblastoma (GBM) is the most aggressive and lethal brain tumor in adults. This study aimed to investigate the functional significance of aryl hydrocarbon receptor nuclear translocator (ARNT) in the pathogenesis of GBM. Analysis of public datasets revealed ARNT is upregulated in GBM tissues compared to lower grade gliomas or normal brain tissues. Higher ARNT expression correlated with the mesenchymal subtype and poorer survival in GBM patients. Silencing ARNT using lentiviral shRNAs attenuated the proliferative, invasive, and stem-like capabilities of GBM cell lines, while ARNT overexpression enhanced these malignant phenotypes. Single-cell RNA sequencing uncovered that ARNT is highly expressed in a stem-like subpopulation and is involved in regulating glycolysis, hypoxia response, and stress pathways. Mechanistic studies found ARNT activates p38 mitogen-activated protein kinase (MAPK) signaling to promote chemoresistance in GBM cells. Disrupting the ARNT/p38α protein interaction via the ARNT PAS-A domain restored temozolomide sensitivity. Overall, this study demonstrates ARNT functions as an oncogenic driver in GBM pathogenesis and represents a promising therapeutic target.

Heliorhodopsin-mediated light-modulation of ABC transporter

Heliorhodopsins (HeRs) have been hypothesized to have widespread functions. Recently, the functions for few HeRs have been revealed; however, the hypothetical functions remain largely unknown. Herein, we investigate light-modulation of heterodimeric multidrug resistance ATP-binding cassette transporters (OmrDE) mediated by Omithinimicrobium cerasi HeR. In this study, we classifiy genes flanking the HeR-encoding genes and identify highly conservative residues for protein-protein interactions. Our results reveal that the interaction between OcHeR and OmrDE shows positive cooperatively sequential binding through thermodynamic parameters. Moreover, light-induced OcHeR upregulates OmrDE drug transportation. Hence, the binding may be crucial to drug resistance in O. cerasi as it survives in a drug-containing habitat. Overall, we unveil a function of HeR as regulatory rhodopsin for multidrug resistance. Our findings suggest potential applications in optogenetic technology.

Best practices for machine learning in antibody discovery and development

In the past 40 years, therapeutic antibody discovery and development have advanced considerably, with machine learning (ML) offering a promising way to speed up the process by reducing costs and the number of experiments required. Recent progress in ML-guided antibody design and development (D&D) has been hindered by the diversity of data sets and evaluation methods, which makes it difficult to conduct comparisons and assess utility. Establishing standards and guidelines will be crucial for the wider adoption of ML and the advancement of the field. This perspective critically reviews current practices, highlights common pitfalls and proposes method development and evaluation guidelines for various ML-based techniques in therapeutic antibody D&D. Addressing challenges across the ML process, best practices are recommended for each stage to enhance reproducibility and progress.

In Silico Analysis of Protein-Protein Interactions of Putative Endoplasmic Reticulum Metallopeptidase 1 in Schizosaccharomyces pombe

Ermp1 is a putative metalloprotease from Schizosaccharomyces pombe and a member of the Fxna peptidases. Although their function is unknown, orthologous proteins from rats and humans have been associated with the maturation of ovarian follicles and increased ER stress. This study focuses on proposing the first prediction of PPI by comparison of the interologues between humans and yeasts, as well as the molecular docking and dynamics of the M28 domain of Ermp1 with possible target proteins. As results, 45 proteins are proposed that could interact with the metalloprotease. Most of these proteins are related to the transport of Ca2+ and the metabolism of amino acids and proteins. Docking and molecular dynamics suggest that the M28 domain of Ermp1 could hydrolyze leucine and methionine residues of Amk2, Ypt5 and Pex12. These results could support future experimental investigations of other Fxna peptidases, such as human ERMP1.

Therapeutic potential of proteases in acute lung injury and respiratory distress syndrome via TLR4/Nrf2/NF-kB signaling modulation

The COVID-19 pandemic has drawn attention to acute lung injury and respiratory distress syndrome as major causes of death, underscoring the urgent need for effective treatments. Protease enzymes possess a wide range of beneficial effects, including antioxidant, anti-inflammatory, antifibrotic, and fibrinolytic effects. This study aimed to evaluate the potential therapeutic effects of bacterial protease and chymotrypsin in rats in mitigating acute lung injury induced by lipopolysaccharide. Molecular docking was employed to investigate the inhibitory effect of bacterial protease and chymotrypsin on TLR-4, the receptor for lipopolysaccharide. Bacterial protease restored TLR-4, Nrf2, p38 MAPK, NF-kB, and IKK-β levels to normal levels, while chymotrypsin normalized TLR-4, IKK-β, IL-6, and IL-17 levels. The expression of TGF-β, caspase-3, and VEGF in the bacterial protease- and chymotrypsin-treated groups was markedly reduced. Our results suggest that both therapies ameliorate LPS-induced acute lung injury and modulate the TLR4/Nrf2/NF-k signaling pathway. Each protease exhibited distinct mechanisms, with bacterial protease showing a better response to oxidative stress, edema, and fibrosis, whereas chymotrypsin provided a better response in the acute phase and innate immunity. These findings highlight the potential of each protease as a promising therapeutic option for acute lung injury and respiratory distress syndrome.

The functional and structural characterisation of the bZIP transcription factors from Myristica fragrans Houtt. associated to plant disease-resistant defence: An insight from transcriptomics and computational modelling

The bZIP transcription factors play crucial roles in various aspects of plant biology, including development, defence mechanisms, senescence, and responses to both biotic and abiotic environmental stresses. Myristica fragrans Houtt. transcriptome analysis has identified 15 bZIP transcription factors, each exhibiting major conserved domains and motifs such as BRLZ, MFMR, and DOG1. Functional characterisation of these identified MfbZIP factors indicates their predominant localisation within the nucleus. Phylogenetic analysis reveals that MfbZIP factors cluster into three subgroups alongside annotated bZIP sequences from Magnolia sinica and Arabidopsis thaliana. Moreover, gene ontology (GO) analysis highlights several key functions of MfbZIP, including involvement in defence responses, abscisic acid-induced signalling pathways, and DNA-binding transcription factor activity. Further investigation through KEGG pathway analysis reveals that the amino acid sequences of MfbZIP contain binding motifs for proteins such as TGA, implicated in plant hormone signal transduction pathways associated with disease resistance. To confirm the disease-defence-related activity of the TGA binding protein within MfbZIP, we employed amino acid sequences for 3-D ab initio modelling. Subsequently, we analysed TGA-NPR1 interactions using docking and molecular dynamics simulation analysis. These analyses shed light on the functional and structural aspects of TGA, demonstrating its stable association with NPR1 protein and its significance in the expression of PR1 protein, thus playing a pivotal role in defence responses against pathogens.

Simulation and Computational Study of RING Domain Mutants of BRCA1 and Ube2k in AD/PD Pathophysiology

Lysine-based post-translational modification (PTM) such as acylation, acetylation, deamination, methylation, SUMOylation, and ubiquitination has proven to be a major regulator of gene expression, chromatin structure, protein stability, protein-protein interaction, protein degradation, and cellular localization. However, besides all the PTMs, ubiquitination stands as the second most common PTM after phosphorylation that is involved in the etiology of neurodegenerative diseases (NDDs) namely, Alzheimer's disease (AD) and Parkinson's disease (PD). NDDs are characterized by the accumulation of misfolded protein aggregates in the brain that lead to disease-related gene mutation and irregular protein homeostasis. The ubiquitin-proteasome system (UPS) is in charge of degrading these misfolded proteins, which involve an interplay of E1, E2, E3, and deubiquitinase enzymes. Impaired UPS has been commonly observed in NDDs and E3 ligases are the key members of the UPS, thus, dysfunction of the same can accelerate the neurodegeneration process. Therefore, the aim of this study is firstly, to find E3 ligases that are common in both AD and PD through data mining. Secondly, to study the impact of mutation on its structure and function. The study deciphered 74 E3 ligases that were common in both AD and PD. Later, 10 hub genes were calculated of which protein-protein interaction, pathway enrichment, lysine site prediction, domain, and motif analysis were performed. The results predicted BRCA1, PML, and TRIM33 as the top three putative lysine-modified E3 ligases involved in AD and PD pathogenesis. However, based on structural characterization, BRCA1 was taken further to study RING domain mutation that inferred K32Y, K32L, K32C, K45V, K45Y, and K45G as potential mutants that alter the structural and functional ability of BRCA1 to interact with Ube2k, E2-conjugating enzyme. The most probable mutant observed after molecular dynamics simulation of 50 ns is K32L. Therefore, our study concludes BRCA1, a potential E3 ligase common in AD and PD, and RING domain mutation at sites K32 and K45 possibly disturbs its interaction with its E2, Ube2k.

Recognition of granulocyte-macrophage colony-stimulating factor by specific S100 proteins

Granulocyte-macrophage colony-stimulating factor (GM-CSF) is a pleiotropic myelopoietic growth factor and proinflammatory cytokine, clinically used for multiple indications and serving as a promising target for treatment of many disorders, including cancer, multiple sclerosis, rheumatoid arthritis, psoriasis, asthma, COVID-19. We have previously shown that dimeric Ca2+-bound forms of S100A6 and S100P proteins, members of the multifunctional S100 protein family, are specific to GM-CSF. To probe selectivity of these interactions, the affinity of recombinant human GM-CSF to dimeric Ca2+-loaded forms of 18 recombinant human S100 proteins was studied by surface plasmon resonance spectroscopy. Of them, only S100A4 protein specifically binds to GM-CSF with equilibrium dissociation constant, Kd, values of 0.3-2 μM, as confirmed by intrinsic fluorescence and chemical crosslinking data. Calcium removal prevents S100A4 binding to GM-CSF, whereas monomerization of S100A4/A6/P proteins disrupts S100A4/A6 interaction with GM-CSF and induces a slight decrease in S100P affinity for GM-CSF. Structural modelling indicates the presence in the GM-CSF molecule of a conserved S100A4/A6/P-binding site, consisting of the residues from its termini, helices I and III, some of which are involved in the interaction with GM-CSF receptors. The predicted involvement of the 'hinge' region and F89 residue of S100P in GM-CSF recognition was confirmed by mutagenesis. Examination of S100A4/A6/P ability to affect GM-CSF signaling showed that S100A4/A6 inhibit GM-CSF-induced suppression of viability of monocytic THP-1 cells. The ability of the S100 proteins to modulate GM-CSF activity is relevant to progression of various neoplasms and other diseases, according to bioinformatics analysis. The direct regulation of GM-CSF signaling by extracellular forms of the S100 proteins should be taken into account in the clinical use of GM-CSF and development of the therapeutic interventions targeting GM-CSF or its receptors.

Transcription factor FOXM1 specifies chromatin DNA to extracellular vesicles

Extracellular vesicle DNAs (evDNAs) hold significant diagnostic value for various diseases and facilitate transcellular transfer of genetic material. Our study identifies transcription factor FOXM1 as a mediator for directing chromatin genes or DNA fragments (termed FOXM1-chDNAs) to extracellular vesicles (EVs). FOXM1 binds to MAP1LC3/LC3 in the nucleus, and FOXM1-chDNAs, such as the DUX4 gene and telomere DNA, are designated by FOXM1 binding and translocated to the cytoplasm before being released to EVs through the secretory autophagy during lysosome inhibition (SALI) process involving LC3. Disrupting FOXM1 expression or the SALI process impairs FOXM1-chDNAs incorporation into EVs. FOXM1-chDNAs can be transmitted to recipient cells via EVs and expressed in recipient cells when they carry functional genes. This finding provides an example of how chromatin DNA fragments are specified to EVs by transcription factor FOXM1, revealing its contribution to the formation of evDNAs from nuclear chromatin. It provides a basis for further exploration of the roles of evDNAs in biological processes, such as horizontal gene transfer.Abbreviation: ATG5: autophagy related 5; CCFs: cytoplasmic chromatin fragments; ChIP: chromatin immunoprecipitation; cytoDNA: cytoplasmic DNA; CQ: chloroquine; FOXM1-DBD: FOXM1 DNA binding domain; DUX4:double homeobox 4; EVs: extracellular vesicles; evDNAs: extracellular vesicle DNAs; FOXM1: forkhead box M1; FOXM1-chDNAs: chromatin DNA fragments directed by FOXM1 to EVs; HGT: horizontal gene transfer; LC3-II: lipid modified LC3; LMNB1: lamin B1; LIR: LC3-interacting region; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; MVBs: multivesicular bodies; M1-binding DNA: a linear DNA containing 72× FOXM1 binding sites; SALI: secretory autophagy during lysosome inhibition; siRNA: small interfering RNA; TetO-DUX4: TetO array-containing DUX4 DNA; TetO: tet operator; TetR: tet repressor.

Construction of novel multi-epitope-based diagnostic biomarker HP16118P and its application in the differential diagnosis of Mycobacterium tuberculosis latent infection

Tuberculosis (TB) is an infectious disease that significantly threatens human health. However, the differential diagnosis of latent tuberculosis infection (LTBI) and active tuberculosis (ATB) remains a challenge for clinicians in early detection and preventive intervention. In this study, we developed a novel biomarker named HP16118P, utilizing 16 helper T lymphocyte (HTL) epitopes, 11 cytotoxic T lymphocyte (CTL) epitopes, and 8 B cell epitopes identified from 15 antigens associated with LTBI-RD using the IEDB database. We analyzed the physicochemical properties, spatial structure, and immunological characteristics of HP16118P using various tools, which indicated that it is a hydrophilic and relatively stable alkaline protein. Furthermore, HP16118P exhibited good antigenicity and immunogenicity, while being non-toxic and non-allergenic, with the potential to induce immune responses. We observed that HP16118P can stimulate the production of high levels of IFN-γ+ T lymphocytes in individuals with ATB, LTBI, and health controls. IL-5 induced by HP16118P demonstrated potential in distinguishing LTBI individuals and ATB patients (p=0.0372, AUC=0.8214, 95% CI [0.5843 to 1.000]) with a sensitivity of 100% and specificity of 71.43%. Furthermore, we incorporated the GM-CSF, IL-23, IL-5, and MCP-3 induced by HP16118P into 15 machine learning algorithms to construct a model. It was found that the Quadratic discriminant analysis model exhibited the best diagnostic performance for discriminating between LTBI and ATB, with a sensitivity of 1.00, specificity of 0.86, and accuracy of 0.93. In summary, HP16118P has demonstrated strong antigenicity and immunogenicity, with the induction of GM-CSF, IL-23, IL-5, and MCP-3, suggesting their potential for the differential diagnosis of LTBI and ATB.

The PKA-CREB1 axis regulates coronavirus proliferation by viral helicase nsp13 association

The COVID-19 pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has posed a worldwide threat in the past 3 years. Although it has been widely and intensively investigated, the mechanism underlying the coronavirus-host interaction requires further elucidation, which may contribute to the development of new antiviral strategies. Here, we demonstrated that the host cAMP-responsive element-binding protein (CREB1) interacts with the non-structural protein 13 (nsp13) of SARS-CoV-2, a conserved helicase for coronavirus replication, both in cells and in lung tissues subjected to SARS-CoV-2 infection. The ATPase and helicase activity of viral nsp13 were shown to be potentiated by CREB1 association, as well as by Protein kinase A (PKA)-mediated CREB1 activation. SARS-CoV-2 replication is significantly suppressed by PKA Cα, cAMP-activated protein kinase catalytic subunit alpha (PRKACA), and CREB1 knockdown or inhibition. Consistently, the CREB1 inhibitor 666-15 has shown significant antiviral effects against both the WIV04 strain and the Omicron strain of the SARS-CoV-2. Our findings indicate that the PKA-CREB1 signaling axis may serve as a novel therapeutic target against coronavirus infection.

IMPORTANCE

In this study, we provide solid evidence that host transcription factor cAMP-responsive element-binding protein (CREB1) interacts directly with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) helicase non-structural protein 13 (nsp13) and potentiate its ATPase and helicase activity. And by live SARS-CoV-2 virus infection, the inhibition of CREB1 dramatically impairs SARS-CoV-2 replication in vivo. Notably, the IC50 of CREB1 inhibitor 666-15 is comparable to that of remdesivir. These results may extend to all highly pathogenic coronaviruses due to the conserved nsp13 sequences in the virus.

Quantitative proteomics analysis identified new interacting proteins of JAL30 in Arabidopsis

Jacalin-related lectins (JALs) are a unique group of plant lectins derived from the jacalin protein family, which play important roles in plant defense responses. JAL30/PBP1 (PYK10 binding protein 1) interacts with inactive PYK10, exerting negative regulatory control over the size of the PYK10 complex, which is formed and activated upon insect or pathogen invasion. However, the precise interplay between JAL30 and other components remains elusive. In this study, we found JAL30 as a nucleocytoplasmic protein, but no obvious phenotype was observed in jal30-1 single mutant. Through immunoprecipitation (IP) enrichment combined with liquid chromatography-tandem mass spectrometry (LC-MS/MS), dozens of new JAL30 interacting proteins were found in addition to several reported ones. Gene Ontology (GO) analysis revealed that these interacting proteins were highly related to the wounding and bacterial stimuli, suggesting their potential involvement in the jasmonate (JA) response. Importantly, the expression of JAL30 was induced by MeJA treatment, further highlighting its relevance in plant defense mechanisms. A novel JAL30 interacting protein, ESM1, was identified and its interaction with JAL30 was confirmed by Co-immunoprecipitation. Moreover, ESM1 was found as an O-GlcNAcylated protein, suggesting that JAL30 may possess glycosylated protein binding ability, particularly in O-GlcNAcylated protein and peptide recognition. Overall, our study provides valuable insights into the interacting protein network and biological function of JAL30, demonstrates the interaction between JAL30 and ESM1, and uncovers the potential significance of JAL30 in plant defense system, potentially through its association with PYK10 complex or JA response. SIGNIFICANCE: The biological functions of lectin proteins, including defense responses, immunity responses, signal transduction, have been well studied. Lectin proteins were also utilized to enrich glycosylated proteins for their specific carbohydrates binding capability. Jacalin-related lectins (JALs) were found to involve in plant defense mechanism. However, it is not yet clear whether JALs could use for enrichment of glycosylated proteins. In this study, we used label-free quantification method to identify interacting proteins of JAL30. A novel interacting protein, ESM1, as an O-GlcNAcylated protein was found. ESM1 has been reported to take part in defense against insect herbivory. Therefore, our findings provided experimental evidence to confirm that JALs have potential to be developed as the bio-tools to enrich glycosylated proteins. Finally, our data not only illustrated the vital biological role of JALs in plants, but also verified unique function of JAL30 in recognizing O-GlcNAcylated proteins.

PRMT4 interacts with NCOA4 to inhibit ferritinophagy in cisplatin-induced acute kidney injury

Cisplatin-induced acute kidney injury (AKI) is commonly seen in the clinical practice, and ferroptosis, a type of non-apoptotic cell death, plays a pivotal role in it. Previous studies suggested that protein arginine methyltransferase 4 (PRMT4) was incorporated in various bioprocesses, but its role in renal injuries has not been investigated. Our present study showed that PRMT4 was highly expressed in renal proximal tubular cells, and it was downregulated in cisplatin-induced AKI. Besides, genetic disruption of PRMT4 exacerbated, while its overexpression attenuated, cisplatin-induced redox injuries in renal proximal epithelia. Mechanistically, our work showed that PRMT4 interacted with NCOA4 to inhibit ferritinophagy, a type of selective autophagy favoring lipid peroxidation to accelerate ferroptosis. Taken together, our study demonstrated that PRMT4 interacted with NCOA4 to attenuate ferroptosis in cisplatin-induced AKI, suggesting that PRMT4 might present as a new therapeutic target for cisplatin-related nephropathy.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Designing a Novel Multiepitope Vaccine from the Human Papilloma Virus E1 and E2 Proteins for Indonesia with Immunoinformatics and Molecular Dynamics Approaches

One of the deadliest malignant cancer in women globally is cervical cancer. Specifically, cervical cancer is the second most common type of cancer in Indonesia. The main infectious agent of cervical cancer is the human papilloma virus (HPV). Although licensed prophylactic vaccines are available, cervical cancer cases are on the rise. Therapy using multiepitope-based vaccines is a very promising therapy for cervical cancer. This study aimed to develop a multiepitope vaccine based on the E1 and E2 proteins of HPV 16, 18, 45, and 52 using in silico. In this study, we develop a novel multiepitope vaccine candidate using an immunoinformatic approach. We predicted the epitopes of the cytotoxic T lymphocyte (CTL) and helper T lymphocyte (HTL) and evaluated their immunogenic properties. Population coverage analysis of qualified epitopes was conducted to determine the successful use of the vaccine worldwide. The epitopes were constructed into a multiepitope vaccine by using AAY linkers between the CTL epitopes and GPGPG linkers between the HTL epitopes. The tertiary structure of the multiepitope vaccine was modeled with AlphaFold and was evaluated by Prosa-web. The results of vaccine construction were analyzed for B-cell epitope prediction, molecular docking with Toll like receptor-4 (TLR4), and molecular dynamics simulation. The results of epitope prediction obtained 4 CTL epitopes and 7 HTL epitopes that are eligible for construction of multiepitope vaccines. Prediction of the physicochemical properties of multiepitope vaccines obtained good results for recombinant protein production. The interaction showed that the interaction of the multiepitope vaccine-TLR4 complex is stable based on the binding free energy value -106.5 kcal/mol. The results of the immune response simulation show that multiepitope vaccine candidates could activate the adaptive and humoral immune systems and generate long-term B-cell memory. According to these results, the development of a multiepitope vaccine with a reverse vaccinology approach is a breakthrough to develop potential cervical cancer therapeutic vaccines.

A Conserved Intramolecular Ion-Pair Plays a Critical but Divergent Role in Regulation of Dimerization and Transport Function among the Monoamine Transporters

The monoamine transporters, including the serotonin transporter (SERT), dopamine transporter (DAT), and norepinephrine transporter (NET), are the therapeutic targets for the treatment of many neuropsychiatric disorders. Despite significant progress in characterizing the structures and transport mechanisms of these transporters, the regulation of their transport functions through dimerization or oligomerization remains to be understood. In the present study, we identified a conserved intramolecular ion-pair at the third extracellular loop (EL3) connecting TM5 and TM6 that plays a critical but divergent role in the modulation of dimerization and transport functions among the monoamine transporters. The disruption of the ion-pair interactions by mutations induced a significant spontaneous cross-linking of a cysteine mutant of SERT and an increase in cell surface expression but with an impaired specific transport activity. On the other hand, similar mutations of the corresponding ion-pair residues in both DAT and NET resulted in an opposite effect on their oxidation-induced dimerization, cell surface expression, and transport function. Reversible biotinylation experiments indicated that the ion-pair mutations slowed down the internalization of SERT but stimulated the internalization of DAT. In addition, cysteine accessibility measurements for monitoring SERT conformational changes indicated that substitution of the ion-pair residues resulted in profound effects on the rate constants for cysteine modification in both the extracellular and cytoplasmatic substrate permeation pathways. Furthermore, molecular dynamics simulations showed that the ion-pair mutations increased the interfacial interactions in a SERT dimer but decreased it in a DAT dimer. Taken together, we propose that the transport function is modulated by the equilibrium between monomers and dimers on the cell surface, which is regulated by a potential compensatory mechanism but with different molecular solutions among the monoamine transporters. The present study provided new insights into the structural elements regulating the transport function of the monoamine transporters through their dimerization.

Computational Insights into the Interaction of the Conserved Cysteine-Noose Domain of the Human Respiratory Syncytial Virus G Protein with the Canonical Fractalkine Binding site of Transmembrane Receptor CX3CR1 Isoforms

The human Respiratory Syncytial Virus (hRSV) stands as one of the most common causes of acute respiratory diseases. The infectivity of this virus is intricately linked to its membrane proteins, notably the attachment glycoprotein (G protein). The latter plays a key role in facilitating the attachment of hRSV to respiratory tract epithelial cells, thereby initiating the infection process. The present study aimed to characterize the interaction of the conserved cysteine-noose domain of hRSV G protein (cndG) with the transmembrane CX3C motif chemokine receptor 1 (CX3CR1) isoforms using computational tools of molecular modeling, docking, molecular dynamics simulations, and binding free energy calculations. From MD simulations of the molecular system embedded in the POPC lipid bilayer, we showed a stable interaction of cndG with the canonical fractalkine binding site in the N-terminal cavity of the CX3CR1 isoforms and identified that residues in the extracellular loop 2 (ECL2) region and Glu279 of this receptor are pivotal for the stabilization of CX3CR1/cndG binding, corroborating what was reported for the interaction of the chemokine fractalkine with CX3CR1 and its structure homolog US28. Therefore, the results presented here contribute by revealing key structural points for the CX3CR1/G interaction, allowing us to better understand the biology of hRSV from its attachment process and to develop new strategies to combat it.

Antibody-mediated targeting of human microglial leukocyte Ig-like receptor B4 attenuates amyloid pathology in a mouse model

Microglia help limit the progression of Alzheimer's disease (AD) by constraining amyloid-β (Aβ) pathology, effected through a balance of activating and inhibitory intracellular signals delivered by distinct cell surface receptors. Human leukocyte Ig-like receptor B4 (LILRB4) is an inhibitory receptor of the immunoglobulin (Ig) superfamily that is expressed on myeloid cells and recognizes apolipoprotein E (ApoE) among other ligands. Here, we find that LILRB4 is highly expressed in the microglia of patients with AD. Using mice that accumulate Aβ and carry a transgene encompassing a portion of the LILR region that includes LILRB4, we corroborated abundant LILRB4 expression in microglia wrapping around Aβ plaques. Systemic treatment of these mice with an anti-human LILRB4 monoclonal antibody (mAb) reduced Aβ load, mitigated some Aβ-related behavioral abnormalities, enhanced microglia activity, and attenuated expression of interferon-induced genes. In vitro binding experiments established that human LILRB4 binds both human and mouse ApoE and that anti-human LILRB4 mAb blocks such interaction. In silico modeling, biochemical, and mutagenesis analyses identified a loop between the two extracellular Ig domains of LILRB4 required for interaction with mouse ApoE and further indicated that anti-LILRB4 mAb may block LILRB4-mApoE by directly binding this loop. Thus, targeting LILRB4 may be a potential therapeutic avenue for AD.

Defining Requirements for Heme Binding in PGRMC1 and Identifying Key Elements that Influence Protein Dimerization

Progesterone receptor membrane component 1 (PGRMC1) binds heme via a surface-exposed site and displays some structural resemblance to cytochrome b5 despite their different functions. In the case of PGRMC1, it is the protein interaction with drug-metabolizing cytochrome P450s and the epidermal growth factor receptor that has garnered the most attention. These interactions are thought to result in a compromised ability to metabolize common chemotherapy agents and to enhance cancer cell proliferation. X-ray crystallography and immunoprecipitation data have suggested that heme-mediated PGRMC1 dimers are important for facilitating these interactions. However, more recent studies have called into question the requirement of heme binding for PGRMC1 dimerization. Our study employs spectroscopic and computational methods to probe and define heme binding and its impact on PGRMC1 dimerization. Fluorescence, electron paramagnetic resonance and circular dichroism spectroscopies confirm heme binding to apo-PGRMC1 and were used to demonstrate the stabilizing effect of heme on the wild-type protein. We also utilized variants (C129S and Y113F) to precisely define the contributions of disulfide bonds and direct heme coordination to PGRMC1 dimerization. Understanding the key factors involved in these processes has important implications for downstream protein-protein interactions that may influence the metabolism of chemotherapeutic agents. This work opens avenues for deeper exploration into the physiological significance of the truncated-PGRMC1 model and developing design principles for potential therapeutics to target PGRMC1 dimerization and downstream interactions.

Application of Computational Techniques in Antibody Fc-Fused Molecule Design for Therapeutics

Since the advent of hybridoma technology in the year 1975, it took a decade to witness the first approved monoclonal antibody Orthoclone OKT39 (muromonab-CD3) in the year 1986. Since then, continuous strides have been made to engineer antibodies for specific desired effects. The engineering efforts were not confined to only the variable domains of the antibody but also included the fragment crystallizable (Fc) region that influences the immune response and serum half-life. Engineering of the Fc fragment would have a profound effect on the therapeutic dose, antibody-dependent cell-mediated cytotoxicity as well as antibody-dependent cellular phagocytosis. The integration of computational techniques into antibody engineering designs has allowed for the generation of testable hypotheses and guided the rational antibody design framework prior to further experimental evaluations. In this article, we discuss the recent works in the Fc-fused molecule design that involves computational techniques. We also summarize the usefulness of in silico techniques to aid Fc-fused molecule design and analysis for the therapeutics application.

Impact of carbamazepine on SMARCA4 (BRG1) expression in colorectal cancer: modulation by KRAS mutation status

SMARCA4 is a gene traditionally considered a tumor suppressor. Recent research has however found that SMARCA4 likely promotes cancer growth and is a good target for cancer treatment. The drug carbamazepine, an autophagy inducer, was used on colorectal cancer cell lines, HCT1116 and Hke3 (KRAS mutant and wildtype). Our study finds that Carbamazepine affects SMARCA4 levels and that this effect is different depending on the KRAS mutation status. This study analyzes the effect of carbamazepine on early-stage autophagy via ULK1 as well as simulates the docking of carbamazepine on KRAS, depending on the mutation status. Our study highlights the therapeutic uses of carbamazepine on cancer, and we propose that carbamazepine in conjunction with other chemotherapies may prove useful in targeting KRAS-mutated colorectal cancer.

Identification and functional characterization of putative ligand binding domain(s) of JlpA protein of Campylobacter jejuni

Among the major Surface Exposed Colonization Proteins (SECPs) of Campylobacter jejuni (C. jejuni), Jejuni lipoprotein A (JlpA) plays a crucial role in host cell adhesion specifically by binding to the N-terminal domain of the human heat shock protein 90α (Hsp90α-NTD). Although the JlpA binding to Hsp90α activates NF-κB and p38 MAP kinase pathways, the underlying mechanism of JlpA association with the cellular receptor remains unclear. To this end, we predicted two potential receptor binding sites within the C-terminal domain of JlpA: one spanning from amino acid residues Q332-A354 and the other from S258-T295; however, the latter exhibited weaker binding. To assess the functional attributes of these predicted sequences, we generated two JlpA mutants (JlpAΔ1: S258-T295; JlpAΔ2: Q332-A354) and assessed the Hsp90α-binding affinity-kinetics by in vitro and ex vivo experiments. Our findings confirmed that the residues Q332-A354 are of greater importance in host cell adhesion with a measurable impact on cellular responses. Further, thermal denaturation by circular dichroism (CD) confirmed that the reduced binding affinity of the JlpAΔ2 to Hsp90α is not associated with protein folding or stability. Together, this study provides a possible framework for determining the molecular function of designing rational inhibitors efficiently targeting JlpA.

Sensitive detection of SARS-CoV-2 main protease 3CLpro with an engineered ribonuclease zymogen

Alongside vaccines and antiviral therapeutics, diagnostic tools are a crucial aid in combating the COVID-19 pandemic caused by the etiological agent SARS-CoV-2. All common assays for infection rely on the detection of viral sub-components, including structural proteins of the virion or fragments of the viral genome. Selective pressure imposed by human intervention of COVID-19 can, however, induce viral mutations that decrease the sensitivity of diagnostic assays based on biomolecular structure, leading to an increase in false-negative results. In comparison, mutations are unlikely to alter the function of viral proteins, and viral machinery is under less selective pressure from vaccines and therapeutics. Accordingly, diagnostic assays that rely on biomolecular function can be more robust than ones that rely on biopolymer structure. Toward this end, we used a split intein to create a circular ribonuclease zymogen that is activated by the SARS-CoV-2 main protease, 3CLpro . Zymogen activation by 3CLpro leads to a >300-fold increase in ribonucleolytic activity, which can be detected with a highly sensitive fluorogenic substrate. This coupled assay can detect low nanomolar concentrations of 3CLpro within a timeframe comparable to that of common antigen-detection protocols. More generally, the concept of detecting a protease by activating a ribonuclease could be the basis of diagnostic tools for other indications.

Phase transition of GvpU regulates gas vesicle clustering in bacteria

Gas vesicles (GVs) are microbial protein organelles that support cellular buoyancy. GV engineering has multiple applications, including reporter gene imaging, acoustic control and payload delivery. GVs often cluster into a honeycomb pattern to minimize occupancy of the cytosol. The underlying molecular mechanism and the influence on cellular physiology remain unknown. Using genetic, biochemical and imaging approaches, here we identify GvpU from Priestia megaterium as a protein that regulates GV clustering in vitro and upon expression in Escherichia coli. GvpU binds to the C-terminal tail of the core GV shell protein and undergoes a phase transition to form clusters in subsaturated solution. These properties of GvpU tune GV clustering and directly modulate bacterial fitness. GV variants can be designed with controllable sensitivity to GvpU-mediated clustering, enabling design of genetically tunable biosensors. Our findings elucidate the molecular mechanisms and functional roles of GV clustering, enabling its programmability for biomedical applications.

The contribution and therapeutic implications of IGHMBP2 mutations on IGHMBP2 biochemical activity and ABT1 association

Mutations within immunoglobulin mu DNA binding protein (IGHMBP2), an RNA-DNA helicase, result in SMA with respiratory distress type I (SMARD1) and Charcot Marie Tooth type 2S (CMT2S). The underlying biochemical mechanism of IGHMBP2 is unknown as well as the functional significance of IGHMBP2 mutations in disease severity. Here we report the biochemical mechanisms of IGHMBP2 disease-causing mutations D565N and H924Y, and their potential impact on therapeutic strategies. The IGHMBP2-D565N mutation has been identified in SMARD1 patients, while the IGHMBP2-H924Y mutation has been identified in CMT2S patients. For the first time, we demonstrate a correlation between the altered IGHMBP2 biochemical activity associated with the D565N and H924Y mutations and disease severity and pathology in patients and our Ighmbp2 mouse models. We show that IGHMBP2 mutations that alter the association with activator of basal transcription (ABT1) impact the ATPase and helicase activities of IGHMBP2 and the association with the 47S pre-rRNA 5' external transcribed spacer. We demonstrate that the D565N mutation impairs IGHMBP2 ATPase and helicase activities consistent with disease pathology. The H924Y mutation alters IGHMBP2 activity to a lesser extent while maintaining association with ABT1. In the context of the compound heterozygous patient, we demonstrate that the total biochemical activity associated with IGHMBP2-D565N and IGHMBP2-H924Y proteins is improved over IGHMBP2-D565N alone. Importantly, we demonstrate that the efficacy of therapeutic applications may vary based on the underlying IGHMBP2 mutations and the relative biochemical activity of the mutant IGHMBP2 protein.

Functional and structural characterisation of RimL from Bacillus cereus, a new Nα-acetyltransferase of ribosomal proteins that was wrongly assigned as an aminoglycosyltransferase

Enzymes of the GNAT (GCN5-relate N-acetyltransferases) superfamily are important regulators of cell growth and development. They are functionally diverse and share low amino acid sequence identity, making functional annotation difficult. In this study, we report the function and structure of a new ribosomal enzyme, Nα-acetyl transferase from Bacillus cereus (RimLBC), a protein that was previously wrongly annotated as an aminoglycosyltransferase. Firstly, extensive comparative amino acid sequence analyses suggested RimLBC belongs to a cluster of proteins mediating acetylation of the ribosomal protein L7/L12. To assess if this was the case, several well established substrates of aminoglycosyltransferases were screened. The results of these studies did not support an aminoglycoside acetylating function for RimLBC. To gain further insight into RimLBC biological role, a series of studies that included MALDI-TOF, isothermal titration calorimetry, NMR, X-ray protein crystallography, and site-directed mutagenesis confirmed RimLBC affinity for Acetyl-CoA and that the ribosomal protein L7/L12 is a substrate of RimLBC. Last, we advance a mechanistic model of RimLBC mode of recognition of its protein substrates. Taken together, our studies confirmed RimLBC as a new ribosomal Nα-acetyltransferase and provide structural and functional insights into substrate recognition by Nα-acetyltransferases and protein acetylation in bacteria.

Rational design of B-cell and T-cell multi epitope-based vaccine against Zika virus, an in silico study

The Zika virus (ZKV) is a single-stranded positive-sense, enveloped RNA virus. Zika infection during pregnancy can cause congenital microcephaly, Guillain-Barré syndrome, miscarriage, and other CNS abnormalities. The world needs safe and effective vaccinations to fight against ZIKV infection since vaccination is generally regarded as one of the most effective ways to prevent infectious diseases. In the present work, we used immunoinformatics and docking studies to construct a vaccine containing multi-epitopes using the structural and non-structural proteins of ZKV. The structural models of ZKV proteins (PrE, PrM, NS1, and NS2A) were constructed using Pyre2 and RaptorX servers. The epitopes of B-cell, T-cell (HTL and CTL), and IFN-γ were predicted, and each epitope's immunogenic nature and physiochemical properties were confirmed. As an adjuvant, the CPG-Oligodeoxynucleotide, an agonist of Toll-like receptor 9 (TLR9), is associated to cytotoxic T-lymphocytes (CTL) epitopes via PAPAP linker. To assess the binding affinity and the tendency of the designed vaccine to induce an immune response through TLR9, molecular docking was done. In the next step, molecular dynamics (MD) simulation to 100 nanoseconds (ns) was used to evaluate the stability of the interaction of the designed vaccine with TLR9. The designed vaccine is predicted to be highly antigenic, non-toxic, soluble, and stable with low flexibility in MD simulation. MD studies indicated that the finalized vaccine-TLR9 docked complex was stable during simulation time. The vaccine construct is able to stimulate both humoral and cellular immune responses. We suppose that our constructed model of the vaccine may have the ability to induce the host immune response against ZKV. Further studies, including in vitro and in vivo experimental analyses, are needed to prove the constructed vaccine's efficacy with multi-epitopes.Communicated by Ramaswamy H. Sarma.

Introduction of protein vaccine candidate based on AP65, AP33, and α-actinin proteins against Trichomonas vaginalis parasite: an immunoinformatics design

BACKGROUND

Trichomonas vaginalis is the most common nonviral sexually transmitted disease (STI) worldwide. Vaccination is generally considered to be one of the most effective methods of preventing infectious diseases. Using AP65, AP33 and α-actinin proteins, this research aims to develop a protein vaccine against Trichomonas vaginalis.

METHODS

Based on the B-cell and T-cell epitope prediction servers, the most antigenic epitopes were selected, and with the necessary evaluations, epitope-rich domains of three proteins, AP65, AP33, and α-actinin, were selected and linked. Subsequently, the ability of the vaccine to interact with toll-like receptors 2 and 4 (TLR2 and TLR4) was assessed. The stability of the interactions was also studied by molecular dynamics for a duration of 100 nanoseconds.

RESULTS

The designed protein consists of 780 amino acids with a molecular weight of 85247.31 daltons. The results of the interaction of the vaccine candidate with TLR2 and TLR4 of the immune system also showed that there are strong interactions between the vaccine candidate protein with TLR2 (-890.7 kcal mol-1) and TLR4 (-967.3 kcal mol-1). All parameters studied to evaluate the stability of the protein structure and the protein-TLR2 and protein-TLR4 complexes showed that the structure of the vaccine candidate protein is stable alone and in complex with the immune system receptors. Investigation of the ability of the designed protein to induce an immune response using the C-ImmSim web server also showed that the designed protein is capable of stimulating B- and T-cell lymphocytes to produce the necessary cytokines and antibodies against Trichomonas vaginalis.

CONCLUSIONS

Overall, our vaccine may have potential protection against Trichomonas vaginalis. However, for experimental in vivo and in vitro studies, it may be a good vaccine candidate.

Study on the structure-activity relationship of rice immunopeptides based on molecular docking

Research on food-derived immunoregulatory peptides has attracted increasing attention of scientists worldwide. However, the structure-activity relationship of rice immunopeptides was not clearly. Herein, 114 rice immunopeptides were obtained by simulating the enzymatic hydrolysis of rice proteins and were further analyzed by NetMHCIipan-4.0. Subsequently, the molecular docking was used to simulate the binding of immunoreactive peptides to major histocompatibility complex class II (MHC-II) molecules. Results show that S, R, D, E, and T amino acid could easily form hydrogen bonds with MHC-II molecules, thus enhancing innate and adaptive immunity. Finally, glucose-modified rice immunopeptides were to investigate the binding of the peptides with MHC-II molecules after glycosylation modification; this provided a theoretical basis for the targeted modification of the generated immunopeptides. All in all, the present study provides a theoretical foundation to further utilize rice processing byproducts and other food products to enhance immunity.

Integrating Reinforcement Learning and Monte Carlo Tree Search for enhanced neoantigen vaccine design

Recent advances in cancer immunotherapy have highlighted the potential of neoantigen-based vaccines. However, the design of such vaccines is hindered by the possibility of weak binding affinity between the peptides and the patient's specific human leukocyte antigen (HLA) alleles, which may not elicit a robust adaptive immune response. Triggering cross-immunity by utilizing peptide mutations that have enhanced binding affinity to target HLA molecules, while preserving their homology with the original one, can be a promising avenue for neoantigen vaccine design. In this study, we introduced UltraMutate, a novel algorithm that combines Reinforcement Learning and Monte Carlo Tree Search, which identifies peptide mutations that not only exhibit enhanced binding affinities to target HLA molecules but also retains a high degree of homology with the original neoantigen. UltraMutate outperformed existing state-of-the-art methods in identifying affinity-enhancing mutations in an independent test set consisting of 3660 peptide-HLA pairs. UltraMutate further showed its applicability in the design of peptide vaccines for Human Papillomavirus and Human Cytomegalovirus, demonstrating its potential as a promising tool in the advancement of personalized immunotherapy.

Iterative In Silico Screening for Optimizing Stable Conformation of Anti-SARS-CoV-2 Nanobodies

Nanobodies (Nbs or VHHs) are single-domain antibodies (sdAbs) derived from camelid heavy-chain antibodies. Nbs have special and unique characteristics, such as small size, good tissue penetration, and cost-effective production, making Nbs a good candidate for the diagnosis and treatment of viruses and other pathologies. Identifying effective Nbs against COVID-19 would help us control this dangerous virus or other unknown variants in the future. Herein, we introduce an in silico screening strategy for optimizing stable conformation of anti-SARS-CoV-2 Nbs. Firstly, various complexes containing nanobodies were downloaded from the RCSB database, which were identified from immunized llamas. The primary docking between Nbs and the SARS-CoV-2 spike protein receptor-binding domain was performed through the ClusPro program, with the manual screening leaving the reasonable conformation to the next step. Then, the binding distances of atoms between the antigen-antibody interfaces were measured through the NeighborSearch algorithm. Finally, filtered nanobodies were acquired according to HADDOCK scores through HADDOCK docking the COVID-19 spike protein with nanobodies under restrictions of calculated molecular distance between active residues and antigenic epitopes less than 4.5 Å. In this way, those nanobodies with more reasonable conformation and stronger neutralizing efficacy were acquired. To validate the efficacy ranking of the nanobodies we obtained, we calculated the binding affinities (∆G) and dissociation constants (Kd) of all screened nanobodies using the PRODIGY web tool and predicted the stability changes induced by all possible point mutations in nanobodies using the MAESTROWeb server. Furthermore, we examined the performance of the relationship between nanobodies' ranking and their number of mutation-sensitive sites (Spearman correlation > 0.68); the results revealed a robust correlation, indicating that the superior nanobodies identified through our screening process exhibited fewer mutation hotspots and higher stability. This correlation analysis demonstrates the validity of our screening criteria, underscoring the suitability of these nanobodies for future development and practical implementation. In conclusion, this three-step screening strategy iteratively in silico greatly improved the accuracy of screening desired nanobodies compared to using only ClusPro docking or default HADDOCK docking settings. It provides new ideas for the screening of novel antibodies and computer-aided screening methods.

Noggin promotes osteogenesis in human adipose-derived mesenchymal stem cells via FGFR2/Src/Akt and ERK signaling pathway

The balance between Noggin and bone morphogenetic proteins (BMPs) is important during early development and skeletal regenerative therapies. Noggin binds BMPs in the extracellular space, thereby preventing BMP signaling. However, Noggin may affect cell response not necessarily through the modulation of BMP signaling, raising the possibility of direct Noggin signaling through yet unspecified receptors. Here we show that in osteogenic cultures of adipose-derived stem cells (ASCs), Noggin activates fibroblast growth factor receptors (FGFRs), Src/Akt and ERK kinases, and it stabilizes TAZ proteins in the presence of dexamethasone. Overall, this leads ASCs to increased expression of osteogenic markers and robust mineral deposition. Our results also indicate that Noggin can induce osteogenic genes expression in normal human bone marrow stem cells and alkaline phosphatase activity in normal human dental pulp stem cells. Besides, Noggin can specifically activate FGFR2 in osteosarcoma cells. We believe our findings open new research avenues to further explore the involvement of Noggin in cell fate modulation by FGFR2/Src/Akt/ERK signaling and potential applications of Noggin in bone regenerative therapies.

Affinity fine-tuning anti-CAIX CAR-T cells mitigate on-target off-tumor side effects

One of the major hurdles that has hindered the success of chimeric antigen receptor (CAR) T cell therapies against solid tumors is on-target off-tumor (OTOT) toxicity due to sharing of the same epitopes on normal tissues. To elevate the safety profile of CAR-T cells, an affinity/avidity fine-tuned CAR was designed enabling CAR-T cell activation only in the presence of a highly expressed tumor associated antigen (TAA) but not when recognizing the same antigen at a physiological level on healthy cells. Using direct stochastic optical reconstruction microscopy (dSTORM) which provides single-molecule resolution, and flow cytometry, we identified high carbonic anhydrase IX (CAIX) density on clear cell renal cell carcinoma (ccRCC) patient samples and low-density expression on healthy bile duct tissues. A Tet-On doxycycline-inducible CAIX expressing cell line was established to mimic various CAIX densities, providing coverage from CAIX-high skrc-59 tumor cells to CAIX-low MMNK-1 cholangiocytes. Assessing the killing of CAR-T cells, we demonstrated that low-affinity/high-avidity fine-tuned G9 CAR-T has a wider therapeutic window compared to high-affinity/high-avidity G250 that was used in the first anti-CAIX CAR-T clinical trial but displayed serious OTOT effects. To assess the therapeutic effect of G9 on patient samples, we generated ccRCC patient derived organotypic tumor spheroid (PDOTS) ex vivo cultures and demonstrated that G9 CAR-T cells exhibited superior efficacy, migration and cytokine release in these miniature tumors. Moreover, in an RCC orthotopic mouse model, G9 CAR-T cells showed enhanced tumor control compared to G250. In summary, G9 has successfully mitigated OTOT side effects and in doing so has made CAIX a druggable immunotherapeutic target.

Investigating the interaction mechanism between gliadin and lysozyme through multispectroscopic analysis and molecular dynamic simulations

The development of stable nanocomplexes based on gliadin and other biopolymers shows potential applications as delivery vehicles in the food industry. However, there is limited study specifically targeting the gliadin-lysozyme system, and their underlying interaction mechanism remains poorly understood. Therefore, the objective of this study was to investigate the binding mechanism between gliadin and lysozyme using a combination of multispectroscopic methods and molecular dynamic simulations. Stable gliadin-lysozyme complex nanoparticles were prepared using an anti-solvent precipitation method with a gliadin-to-lysozyme mass ratio of 2:1 and pH 4.0. The characteristic changes in the UV-visible spectrum of gliadin induced by lysozyme confirmed the complex formation. The analyses of fluorescence, FT-IR spectra, and dissociation tests demonstrated the indispensability of hydrophobic, electrostatic, and hydrogen bonding interactions in the preparation of the composites. Scanning electron microscopy revealed that the surface morphology of the nanoparticles changed from smooth and spherical to rough and irregular with the addition of lysozyme. Furthermore, molecular dynamic simulations suggested that lysozyme bound to the hydrophobic region of gliadin and hydrogen bonding was crucial for the stability of the complex. These findings contribute to the advancement of gliadin-lysozyme complex nanoparticles as an efficient delivery system for encapsulating bioactive compounds in food industry.

An Aedes-Anopheles Vaccine Candidate Supplemented with BCG Epitopes Against the Aedes and Anopheles Genera to Overcome Hypersensitivity to Mosquito Bites

BACKGROUND

Skeeter syndrome is a severe local allergic response to mosquito bites that is accompanied by considerable inflammation and, in some cases, a systemic response like fever. People with the syndrome develop serious allergies, ranging from rashes to anaphylaxis or shock. The few available studies on mosquito venom immunotherapy have utilized whole-body preparations and small sample sizes. Still, owing to their little success, vaccination remains a promising alternative as well as a permanent solution for infections like Skeeter's.

METHODS

This study, therefore, illustrated the construction of an epitope-based vaccine candidate against Skeeter Syndrome using established immunoinformatic techniques. We selected three species of mosquitoes, Anopheles melas, Anopheles funestus, and Aedes aegypti, to derive salivary antigens usually found in mosquito bites. Our construct was also supplemented with bacterial epitopes known to elicit a strong TH1 response and suppress TH2 stimulation that is predicted to reduce hypersensitivity against the bites.

RESULTS

A quality factor of 98.9496, instability index of 38.55, aliphatic index of 79.42, solubility of 0.934747, and GRAVY score of -0.02 indicated the structural (tertiary and secondary) stability, thermostability, solubility, and hydrophilicity of the construct, respectively. The designed Aedes-Anopheles vaccine (AAV) candidate was predicted to be flexible and less prone to deformability with an eigenvalue of 1.5911e-9 and perfected the human immune response against Skeeter (hypersensitivity) and many mosquito-associated diseases as we noted the production of 30,000 Th1 cells per mm3 with little (insignificant production of Th2 cells. The designed vaccine also revealed stable interactions with the pattern recognition receptors of the host. The TLR2/vaccine complex interacted with a free energy of - 1069.2 kcal/mol with 26 interactions, whereas the NLRP3/vaccine complex interacted with a free energy of - 1081.2 kcal/mol with 16 molecular interactions.

CONCLUSION

Although being a pure in-silico study, the in-depth analysis performed herein speaks volumes of the potency of the designed vaccine candidate predicting that the proposition can withstand rigorous in-vitro and in-vivo clinical trials and may proceed to become the first preventative immunotherapy against mosquito bite allergy.

Decoding the genetic landscape of allergic rhinitis: a comprehensive network analysis revealing key genes and potential therapeutic targets

BACKGROUND

Allergic Rhinitis (AR), an inflammatory affliction impacting the upper respiratory tract, has been registering a substantial surge in incidence across the globe.

METHODS

We embarked on examination of differentially expressed genes (DEGs) and the Weighted Gene Co-Expression Network Analysis (WGCNA). With this armory of genes identified, we engaged the tools of Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Our study continued with the establishment of a protein-protein interaction (PPI) network and the application of LASSO regression. Finally, we leveraged a docking model to elucidate potential drug-gene interactions involving these key genes.

RESULTS

Through WGCNA and different express genes screening, PPI network was performed, identifying top 20 key genes, including CD44, CD69, CD274. LASSO regression identified three independent factors, STARD5, CST1, and CHAC1, that were significantly associated with AR. A predictive model was developed with an AUC value over 0.75. Also, 105 potential therapeutic agents were discovered, including Fluorouracil, Cyclophosphamide, Doxorubicin, and Hydrocortisone, offering promising therapeutic strategies for AR.

CONCLUSION

By fuzing DEGs with key genes derived from WGCNA, this study has illuminated a comprehensive network of gene interactions involved in the pathogenesis of AR, paving the way for future biomarker and therapeutic target discovery in AR.

Designing a T-cell epitope-based vaccine using in silico approaches against the Sal k 1 allergen of Salsola kali plant

Allergens originated from Salsola kali (Russian thistle) pollen grains are one of the most important sources of aeroallergens causing pollinosis in desert and semi-desert regions. T-cell epitope-based vaccines (TEV) are more effective among different therapeutic approaches developed to alleviate allergic diseases. The physicochemical properties, and B as well as T cell epitopes of Sal k 1 (a major allergen of S. kali) were predicted using immunoinformatic tools. A TEV was constructed using the linkers EAAAK, GPGPG and the most suitable CD4+ T cell epitopes. RS04 adjuvant was added as a TLR4 agonist to the amino (N) and carboxyl (C) terminus of the TEV protein. The secondary and tertiary structures, solubility, allergenicity, toxicity, stability, physicochemical properties, docking with immune receptors, BLASTp against the human and microbiota proteomes, and in silico cloning of the designed TEV were assessed using immunoinformatic analyses. Two CD4+ T cell epitopes of Sal k1 that had high affinity with different alleles of MHC-II were selected and used in the TEV. The molecular docking of the TEV with HLADRB1, and TLR4 showed TEV strong interactions and stable binding pose to these receptors. Moreover, the codon optimized TEV sequence was cloned between NcoI and XhoI restriction sites of pET-28a(+) expression plasmid. The designed TEV can be used as a promising candidate in allergen-specific immunotherapy against S. kali. Nonetheless, effectiveness of this vaccine should be validated through immunological bioassays.

Molecular Insights into the Inhibitory Role of α-Crystallin against γD-Crystallin Aggregation

Cataracts, a major cause of global blindness, contribute significantly to the overall prevalence of blindness. The opacification of the lens, resulting in cataract formation, primarily occurs due to the aggregation of crystallin proteins within the eye lens. Despite the high concentration of these crystallins, they remarkably maintain the lens transparency and refractive index. α-Crystallins (α-crys), acting as chaperones, play a crucial role in preventing crystallin aggregation, although the exact molecular mechanism remains uncertain. In this study, we employed a combination of molecular docking, all-atom molecular dynamics simulations, and advanced free energy calculations to investigate the interaction between γD-crystallin (γD-crys), a major structural protein of the eye lens, and α-crystallin proteins. Our findings demonstrate that α-crys exhibits an enhanced affinity for the NTD2 and CTD4 regions of γD-crys. The NTD2 and CTD4 regions form the interface between the N-terminal domain (NTD) and the C-terminal domain (CTD) of the γD-crys protein. By binding to the interface region between the NTD and CTD of the protein, α-crys effectively inhibits the formation of domain-swapped aggregates and mitigates protein aggregation. Analysis of the Markov state models using molecular dynamics trajectories confirms that minimum free energy conformations correspond to the binding of the α-crystallin domain (ACD) of α-crys to NTD2 and CTD4 that form the interdomain interface.

Computational Analysis of the Tripartite Interaction of Phasins (PhaP4 and 5)-Sigma Factor (σ24)-DNA of Azospirillum brasilense Sp7

Azospirillum brasilense Sp7 produces PHB, which is covered by granule-associated proteins (GAPs). Phasins are the main GAPs. Previous studies have shown phasins can regulate PHB synthesis. When A. brasilense grows under stress conditions, it uses sigma factors to transcribe genes for survival. One of these factors is the σ24 factor. This study determined the possible interaction between phasins and the σ24 factor or phasin-σ24 factor complex and DNA. Three-dimensional structures of phasins and σ24 factor structures were predicted using the I-TASSER and SWISS-Model servers, respectively. Subsequently, a molecular docking between phasins and the σ24 factor was performed using the ClusPro 2.0 server, followed by molecular docking between protein complexes and DNA using the HDOCK server. Evaluation of the types of ligand-receptor interactions was performed using the BIOVIA Discovery Visualizer for three-dimensional diagrams, as well as the LigPlot server to obtain bi-dimensional diagrams. The results showed the phasins (Pha4Abs7 or Pha5Abs7)-σ24 factor complex was bound near the -35 box of the promoter region of the phaC gene. However, in the individual interaction of PhaP5Abs7 and the σ24 factor, with DNA, both proteins were bound to the -35 box. This did not occur with PhaP4Abs7, which was bound to the -10 box. This change could affect the transcription level of the phaC gene and possibly affect PHB synthesis.

Fusing Peptide Epitopes for Advanced Multiplex Serological Testing for SARS-CoV-2 Antibody Detection

The tragic COVID-19 pandemic, which has seen a total of 655 million cases worldwide and a death toll of over 6.6 million seems finally tailing off. Even so, new variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continue to arise, the severity of which cannot be predicted in advance. This is concerning for the maintenance and stability of public health, since immune evasion and increased transmissibility may arise. Therefore, it is crucial to continue monitoring antibody responses to SARS-CoV-2 in the general population. As a complement to polymerase chain reaction tests, multiplex immunoassays are elegant tools that use individual protein or peptide antigens simultaneously to provide a high level of sensitivity and specificity. To further improve these aspects of SARS-CoV-2 antibody detection, as well as accuracy, we have developed an advanced serological peptide-based multiplex assay using antigen-fused peptide epitopes derived from both the spike and the nucleocapsid proteins. The significance of the epitopes selected for antibody detection has been verified by in silico molecular docking simulations between the peptide epitopes and reported SARS-CoV-2 antibodies. Peptides can be more easily and quickly modified and synthesized than full length proteins and can, therefore, be used in a more cost-effective manner. Three different fusion-epitope peptides (FEPs) were synthesized and tested by enzyme-linked immunosorbent assay (ELISA). A total of 145 blood serum samples were used, compromising 110 COVID-19 serum samples from COVID-19 patients and 35 negative control serum samples taken from COVID-19-free individuals before the outbreak. Interestingly, our data demonstrate that the sensitivity, specificity, and accuracy of the results for the FEP antigens are higher than for single peptide epitopes or mixtures of single peptide epitopes. Our FEP concept can be applied to different multiplex immunoassays testing not only for SARS-CoV-2 but also for various other pathogens. A significantly improved peptide-based serological assay may support the development of commercial point-of-care tests, such as lateral-flow-assays.

SARS-CoV-2 ORF7a Mutation Found in BF.5 and BF.7 Sublineages Impacts Its Functions

A feature of the SARS-CoV-2 Omicron subvariants BF.5 and BF.7 that recently circulated mainly in China and Japan was the high prevalence of the ORF7a: H47Y mutation, in which the 47th residue of ORF7a has been mutated from a histidine (H) to a tyrosine (Y). Here, we evaluated the effect of this mutation on the three main functions ascribed to the SARS-CoV-2 ORF7a protein. Our findings show that H47Y mutation impairs the ability of SARS-CoV-2 ORF7a to antagonize the type I interferon (IFN-I) response and to downregulate major histocompatibility complex I (MHC-I) cell surface levels, but had no effect in its anti-SERINC5 function. Overall, our results suggest that the H47Y mutation of ORF7a affects important functions of this protein, resulting in changes in virus pathogenesis.