CABS-flex 2.0: a web server for fast simulations of flexibility of protein structures

Interaction of DisBa01 peptide from Bothrops alternatus venom with BRAF melanoma receptors: Modeling and molecular docking

Metastatic melanoma is highly aggressive and challenging, often leading to a grim prognosis. Its progression is swift, especially when mutations like BRAFV600E continuously activate pathways vital for cell growth and survival. Although several treatments target this mutation, resistance typically emerges over time. In recent decades, research has underscored the potential of snake venoms and peptides as bioactive substances for innovative drugs, including anti-coagulants, anti-microbial, and anti-cancer agents. Leveraging this knowledge, we propose employing a bioinformatics simulation approach to: a) Predict how well a peptide (DisBa01) from Bothrops alternatus snake venom binds to the melanoma receptor BRAFV600E via Molecular Docking. b) Identify the specific peptide binding sites on receptors and analyze their proximity to active receptor sites. c) Evaluate the behavior of resulting complexes through molecular dynamics simulations. d) Assess whether this peptide qualifies as a candidate for anti-melanoma therapy. Our findings reveal that DisBa01 enhances stability in the BRAFV600E melanoma receptor structure by binding to its RGD motif, an interaction absent in the BRAF WT model. Consequently, both docking and molecular dynamics simulations suggest that DisBa01 shows promise as a BRAFV600E inhibitor.

Humanization of the antigen-recognition domain does not impinge on the antigen-binding, cytokine secretion, and antitumor reactivity of humanized nanobody-based CD19-redirected CAR-T cells

BACKGROUND

The immunogenicity of the antigen-recognition domains of chimeric antigen receptor (CAR)-T cells leads to immune responses that may compromise the antitumor effects of the adoptively transferred T cells. Herein, we attempt to humanize a CD19-specific VHH (named H85) using in silico techniques and investigate the impact of antigen-recognition domain humanization on CAR expression and density, cytokine secretion, and cytolytic reactivity of CAR-T cells based on the humanized VHH.

METHODS

H85 was humanized (named HuH85), and then HuH85 was compared with H85 in terms of conformational structure, physicochemical properties, antigenicity and immunogenicity, solubility, flexibility, stability, and CD19-binding capacity using in silico techniques. Next, H85CAR-T cells and HuH85CAR-T cells were developed and CAR expression and surface density were assessed via flow cytometry. Ultimately, the antitumor reactivity and secreted levels of IFN-γ, IL-2, and TNF-α were assessed following the co-cultivation of the CAR-T cells with Ramos, Namalwa, and K562 cells.

RESULTS

In silico findings demonstrated no negative impacts on HuH85 as a result of humanization. Ultimately, H85CAR and HuH85CAR could be surface-expressed on transduced T cells at comparable levels as assessed via mean fluorescence intensity. Moreover, H85CAR-T cells and HuH85CAR-T cells mediated comparable antitumor effects against Ramos and Namalwa cells and secreted comparable levels of IFN-γ, IL-2, and TNF-α following co-cultivation.

CONCLUSION

HuH85 can be used to develop immunotherapeutics against CD19-associated hematologic malignancies. Moreover, HuH85CAR-T cells must be further investigated in vitro and in preclinical xenograft models of CD19+ leukemias and lymphomas before advancing into clinical trials.

Exploring the nuclear proteins, viral capsid protein, and early antigen protein using immunoinformatic and molecular modeling approaches to design a vaccine candidate against Epstein Barr virus

The available Epstein Barr virus vaccine has tirelessly harnessed the gp350 glycoprotein as its target epitope, but the result has not been preventive. Right here, we designed a global multi-epitope vaccine for EBV; with special attention to making sure all strains and preventive antigens are covered. Using a robust computational vaccine design approach, our proposed vaccine is armed with 6-16 mers linear B-cell epitopes, 4-9 mer CTL epitopes, and 8-15 mer HTL epitopes which are verified to induce interleukin 4, 10 & IFN-gamma. We employed deep computational mining coupled with expert intelligence in designing the vaccine, using human Beta defensin-3-which has been reported to induce the same TLRs as EBV-as the adjuvant. The tendency of the vaccine to cause autoimmune disorder is quenched by the assurance that the construct contains no EBNA-1 homolog. The protein vaccine construct exhibited excellent physicochemical attributes such as Aliphatic index 59.55 and GRAVY - 0.710; and a ProsaWeb Z score of - 3.04. Further computational analysis revealed the vaccine docked favorably with EBV indicted TLR 1, 2, 4 & 9 with satisfactory interaction patterns. With global coverage of 85.75% and the stable molecular dynamics result obtained for the best two interactions, we are optimistic that our nontoxic, non-allergenic multi-epitope vaccine will help to ameliorate the EBV-associated diseases-which include various malignancies, tumors, and cancers-preventively.

Aggrescan4D: structure-informed analysis of pH-dependent protein aggregation

Protein aggregation is behind the genesis of incurable diseases and imposes constraints on drug discovery and the industrial production and formulation of proteins. Over the years, we have been advancing the Aggresscan3D (A3D) method, aiming to deepen our comprehension of protein aggregation and assist the engineering of protein solubility. Since its inception, A3D has become one of the most popular structure-based aggregation predictors because of its performance, modular functionalities, RESTful service for extensive screenings, and intuitive user interface. Building on this foundation, we introduce Aggrescan4D (A4D), significantly extending A3D's functionality. A4D is aimed at predicting the pH-dependent aggregation of protein structures, and features an evolutionary-informed automatic mutation protocol to engineer protein solubility without compromising structure and stability. It also integrates precalculated results for the nearly 500,000 jobs in the A3D Model Organisms Database and structure retrieval from the AlphaFold database. Globally, A4D constitutes a comprehensive tool for understanding, predicting, and designing solutions for specific protein aggregation challenges. The A4D web server and extensive documentation are available at https://biocomp.chem.uw.edu.pl/a4d/. This website is free and open to all users without a login requirement.

Multireceptor Analysis for Evaluating the Antidiabetic Efficacy of Karanjin: A Computational Approach

BACKGROUND

Diabetes mellitus, notably type 2, is a rising global health challenge, prompting the need for effective management strategies. Common medications such as metformin, insulin, repaglinide and sitagliptin can induce side effects like gastrointestinal disturbances, hypoglycemia, weight gain and specific organ risks. Plant-derived therapies like Karanjin from Pongamia pinnata present promising alternatives due to their historical use, holistic health benefits and potentially fewer adverse effects. This study employs in silico analysis to explore Karanjin's interactions with diabetes-associated receptors, aiming to unveil its therapeutic potential while addressing the limitations and side effects associated with conventional medications.

METHODOLOGY

The research encompassed the selection of proteins from the Protein Data Bank (PDB), followed by structural refinement processes and optimization. Ligands such as Karanjin and standard drugs were retrieved from PubChem, followed by a comprehensive analysis of their ADMET profiling and pharmacokinetic properties. Protein-ligand interactions were evaluated through molecular docking using AutoDockTools 1.5.7, followed by the analysis of structural stability using coarse-grained simulations with CABS Flex 2.0. Molecular dynamics simulations were performed using Desmond 7.2 and the OPLS4 force field to explore how Karanjin interacts with proteins over 100 nanoseconds, focusing on the dynamics and structural stability.

RESULTS

Karanjin, a phytochemical from Pongamia pinnata, shows superior drug candidate potential compared to common medications, offering advantages in efficacy and reduced side effects. It adheres to drug-likeness criteria and exhibits optimal ADMET properties, including moderate solubility, high gastrointestinal absorption and blood-brain barrier penetration. Molecular docking revealed Karanjin's highest binding energy against receptor 3L2M (Pig pancreatic alpha-amylase) at -9.1 kcal/mol, indicating strong efficacy potential. Molecular dynamics simulations confirmed stable ligand-protein complexes with minor fluctuations in RMSD and RMSF, suggesting robust interactions with receptors 3L2M.

CONCLUSION

Karanjin demonstrates potential in pharmaceutical expansion for treating metabolic disorders such as diabetes, as supported by computational analysis. Prospects for Karanjin in pharmaceutical development include structural modifications for enhanced efficacy and safety. Nanoencapsulation may improve bioavailability and targeted delivery to pancreatic cells, while combination therapies could optimize treatment outcomes in diabetes management. Clinical trials and experimental studies are crucial to validate its potential as a novel therapeutic agent.

Integrating pan-genome and reverse vaccinology to design multi-epitope vaccine against Herpes simplex virus type-1

Herpes simplex virus type-1 (HSV-1), the etiological agent of sporadic encephalitis and recurring oral (sometimes genital) infections in humans, affects millions each year. The evolving viral genome reduces susceptibility to existing antivirals and, thus, necessitates new therapeutic strategies. Immunoinformatics strategies have shown promise in designing novel vaccine candidates in the absence of a clinically licensed vaccine to prevent HSV-1. However, to encourage clinical translation, the HSV-1 pan-genome was integrated with the reverse-vaccinology pipeline for rigorous screening of universal vaccine candidates. Viral targets were screened from 104 available complete genomes. Among 364 proteins, envelope glycoprotein D being an outer membrane protein with a high antigenicity score (> 0.4) and solubility (> 0.6) was selected for epitope screening. A total of 17 T-cell and 4 B-cell epitopes with highly antigenic, immunogenic, non-toxic properties and high global population coverage were identified. Furthermore, 8 vaccine constructs were designed using different combinations of epitopes and suitable linkers. VC-8 was identified as the most potential vaccine candidate regarding chemical and structural stability. Molecular docking revealed high interactive affinity (low binding energy: - 56.25 kcal/mol) of VC-8 with the target elicited by firm intermolecular H-bonds, salt-bridges, and hydrophobic interactions, which was validated with simulations. Compatibility of the vaccine candidate to be expressed in pET-29(a) + plasmid was established by in silico cloning studies. Immune simulations confirmed the potential of VC-8 to trigger robust B-cell, T-cell, cytokine, and antibody-mediated responses, thereby suggesting a promising candidate for the future of HSV-1 prevention.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-04022-6.

Hotspot Wizard-informed engineering of a hyperthermophilic β-glucosidase for enhanced enzyme activity at low temperatures

Hyperthermophilic enzymes serve as an important source of industrial enzymes due to their high thermostability. Unfortunately, most hyperthermophilic enzymes suffer from reduced activity at low temperatures (e.g., ambient temperature), limiting their applicability. In addition, evolving hyperthermophilic enzymes to increase low temperature activity without compromising other desired properties is generally difficult. In the current study, a variant of β-glucosidase from Pyrococcus furiosus (PfBGL) was engineered to enhance enzyme activity at low temperatures through the construction of a saturation mutagenesis library guided by the HotSpot Wizard analysis, followed by its screening for activity and thermostability. From this library construction and screening, one PfBGL mutant, PfBGL-A4 containing Q214S/A264S/F344I mutations, showed an over twofold increase in β-glucosidase activity at 25 and 50°C compared to the wild type, without compromising high-temperature activity, thermostability and substrate specificity. Our experimental and computational characterizations suggest that the findings with PfBGL-A4 may be due to the elevation of local conformational flexibility around the active site, while slightly compacting the global protein structure. This study showcases the potential of HotSpot Wizard-informed engineering of hyperthermophilic enzymes and underscores the interplays among temperature, enzyme activity, and conformational flexibility in these enzymes.

Computational design of novel chimeric multiepitope vaccine against bacterial and viral disease in tilapia (Oreochromis sp.)

Regarding several infectious diseases in fish, multiple vaccinations are not favorable. The chimeric multiepitope vaccine (CMEV) harboring several antigens for multi-disease prevention would enhance vaccine efficiency in terms of multiple disease prevention. Herein, the immunogens of tilapia's seven pathogens including E. tarda, F. columnare, F. noatunensis, S. iniae, S. agalactiae, A. hydrophila, and TiLV were used for CMEV design. After shuffling and annotating the B-cell epitopes, 5,040 CMEV primary protein structures were obtained. Secondary and tertiary protein structures were predicted by AlphaFold2 creating 25,200 CMEV. Proper amino acid alignment in the secondary structures was achieved by the Ramachandran plot. In silico determination of physiochemical and other properties including allergenicity, antigenicity, glycosylation, and conformational B-cell epitopes were determined. The selected CMEV (OSLM0467, OSLM2629, and OSLM4294) showed a predicted molecular weight (MW) of 70 kDa, with feasible sites of N- and O-glycosylation, and a number of potentially conformational B-cell epitope residues. Molecular docking, codon optimization, and in-silico cloning were tested to evaluate the possibility of protein expression. Those CMEVs will further elucidate in vitro and in vivo to evaluate the efficacy and specific immune response. This research will highlight the new era of vaccines designed based on in silico structural vaccine design.

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.

Catalytic specificity and crystal structure of cystathionine γ-lyase from Pseudomonas aeruginosa

The escalating drug resistance among microorganisms underscores the urgent need for innovative therapeutic strategies and a comprehensive understanding of bacteria's defense mechanisms against oxidative stress and antibiotics. Among the recently discovered barriers, the endogenous production of hydrogen sulfide (H2S) via the reverse transsulfuration pathway, emerges as a noteworthy factor. In this study, we have explored the catalytic capabilities and crystal structure of cystathionine γ-lyase from Pseudomonas aeruginosa (PaCGL), a multidrug-opportunistic pathogen chiefly responsible for nosocomial infections. In addition to a canonical L-cystathionine hydrolysis, PaCGL efficiently catalyzes the production of H2S using L-cysteine and/or L-homocysteine as alternative substrates. Comparative analysis with the human enzyme and counterparts from other pathogens revealed distinct structural features within the primary enzyme cavities. Specifically, a distinctly folded entrance loop could potentially modulate the access of substrates and/or inhibitors to the catalytic site. Our findings offer significant insights into the structural evolution of CGL enzymes across different pathogens and provide novel opportunities for developing specific inhibitors targeting PaCGL.

Overexpression of WT1 in all molecular subtypes of breast cancer and its impact on survival: exploring oncogenic and tumor suppressor roles of distinct WT1 isoforms

BACKGROUND

Breast cancer is a highly heterogeneous solid tumor, posing challenges in developing targeted therapies effective for all mammary carcinoma subtypes. WT1 emerges as a promising target for breast cancer therapy due to its potential oncogenic role in various cancer types. Previous works have yielded inconsistent results. Therefore, further studies are needed to clarify the behavior of this complex gene in breast cancer.

METHODS AND RESULTS

In this study, we examined WT1 expression in both Formalin Fixed Paraffin Embedded breast tumors (n = 41) and healthy adjacent tissues (n = 41) samples from newly diagnosed cases of ductal invasive breast cancer. The fold change in gene expression between the tumor and healthy tissue was determined by calculating 2-∆∆Ct. Disease-free survival analysis was computed using the Kaplan-Meier method. To identify the expression levels of different WT1 isoforms, we explored the ISOexpresso database. Relative quantification of the WT1 gene revealed an overexpression of WT1 in most cases. The percentage of patients surviving free of disease at 8 years of follow-up was lower in the group overexpressing WT1 compared to the group with down-regulated WT1.

CONCLUSIONS

Interestingly, this overexpression was observed in all molecular subtypes of invasive breast cancer, underscoring the significance of WT1 as a potential target in all these subtypes. The observed WT1 down-expression in a few cases of invasive breast cancer, associated with better survival outcomes, may correspond to the down-regulation of a particular WT1-KTS (-) isoform: the WT1 A isoform (EX5-/KTS-). The co-expression of this WT1 oncogenic isoform with a regulated WT1- tumor suppressor isoform, such as the major WT1 F isoform (EX5-/KTS +), could also explain such survival outcomes. Due to its capacity to adopt dual roles, it becomes imperative to conduct individual molecular expression profiling of the WT1 gene. Such an approach holds great promise in the development of personalized treatment strategies for breast cancer.

Click, Compute, Create: A Review of Web-based Tools for Enzyme Engineering

Enzyme engineering, though pivotal across various biotechnological domains, is often plagued by its time-consuming and labor-intensive nature. This review aims to offer an overview of supportive in silico methodologies for this demanding endeavor. Starting from methods to predict protein structures, to classification of their activity and even the discovery of new enzymes we continue with describing tools used to increase thermostability and production yields of selected targets. Subsequently, we discuss computational methods to modulate both, the activity as well as selectivity of enzymes. Last, we present recent approaches based on cutting-edge machine learning methods to redesign enzymes. With exception of the last chapter, there is a strong focus on methods easily accessible via web-interfaces or simple Python-scripts, therefore readily useable for a diverse and broad community.

In Silico Insights into the Arsenic Binding Mechanism Deploying Application of Computational Biology-Based Toolsets

An assortment of environmental matrices includes arsenic (As) in its different oxidation states, which is often linked to concerns that pose a threat to public health worldwide. The current difficulty lies in addressing toxicological concerns and achieving sustained detoxification of As. Multiple conventional degradation methods are accessible; however, they are indeed labor-intensive, expensive, and reliant on prolonged laboratory evaluations. Molecular interaction and atomic level degradation mechanisms for enzyme-As exploration are, however, underexplored in those approaches. A feasible approach in this case for tackling this accompanying concern of As might be to cope with undertaking multivalent computational methodologies and tools. This work aimed to provide molecular-level insight into the enzyme-aided As degradation mechanism. AutoDock Vina, CABS-flex 2.0, and Desmond high-performance molecular dynamics simulation (MDS) were utilized in the current investigation to simulate multivalent molecular processes on two protein sets: arsenate reductase (ArsC) and laccase (LAC) corresponding arsenate (ART) and arsenite (AST), which served as model ligands to comprehend binding, conformational, and energy attributes. The structural configurations of both proteins exhibited variability in flexibility and structure framework within the range of 3.5-4.5 Å. The LAC-ART complex exhibited the lowest calculated binding affinity, measuring -5.82 ± 0.01 kcal/mol. Meanwhile, active site residues ILE-200 and HIS-206 were demonstrated to engage in H-bonding with the ART ligand. In contrast to ArsC, the ligand binding affinity of this bound complex was considerably greater. Additional validation of docked complexes was carried out by deploying Desmond MDS of 100 ns to capture protein and ligand conformation behavior. The system achieved stability during the 100 ns simulation run, as confirmed by the average P-L RMSD, which was ∼1 Å. As a preliminary test of the enzyme's ability to catalyze As species, corresponding computational insights might be advantageous for bridging gaps and regulatory consideration.

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.

Searching for Frataxin Function: Exploring the Analogy with Nqo15, the Frataxin-like Protein of Respiratory Complex I from Thermus thermophilus

Nqo15 is a subunit of respiratory complex I of the bacterium Thermus thermophilus, with strong structural similarity to human frataxin (FXN), a protein involved in the mitochondrial disease Friedreich's ataxia (FRDA). Recently, we showed that the expression of recombinant Nqo15 can ameliorate the respiratory phenotype of FRDA patients' cells, and this prompted us to further characterize both the Nqo15 solution's behavior and its potential functional overlap with FXN, using a combination of in silico and in vitro techniques. We studied the analogy of Nqo15 and FXN by performing extensive database searches based on sequence and structure. Nqo15's folding and flexibility were investigated by combining nuclear magnetic resonance (NMR), circular dichroism, and coarse-grained molecular dynamics simulations. Nqo15's iron-binding properties were studied using NMR, fluorescence, and specific assays and its desulfurase activation by biochemical assays. We found that the recombinant Nqo15 isolated from complex I is monomeric, stable, folded in solution, and highly dynamic. Nqo15 does not share the iron-binding properties of FXN or its desulfurase activation function.

Antibody targeting of conserved sites of vulnerability on the SARS-CoV-2 spike receptor-binding domain

Given the continuous emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VoCs), immunotherapeutics that target conserved epitopes on the spike (S) glycoprotein have therapeutic advantages. Here, we report the crystal structure of the SARS-CoV-2 S receptor-binding domain (RBD) at 1.95 Å and describe flexibility and distinct conformations of the angiotensin-converting enzyme 2 (ACE2)-binding site. We identify a set of SARS-CoV-2-reactive monoclonal antibodies (mAbs) with broad RBD cross-reactivity including SARS-CoV-2 Omicron subvariants, SARS-CoV-1, and other sarbecoviruses and determine the crystal structures of mAb-RBD complexes with Ab246 and CR3022 mAbs targeting the class IV site, WRAIR-2134, which binds the recently designated class V epitope, and WRAIR-2123, the class I ACE2-binding site. The broad reactivity of class IV and V mAbs to conserved regions of SARS-CoV-2 VoCs and other sarbecovirus provides a framework for long-term immunotherapeutic development strategies.

Structure prediction of linear and cyclic peptides using CABS-flex

The structural modeling of peptides can be a useful aid in the discovery of new drugs and a deeper understanding of the molecular mechanisms of life. Here we present a novel multiscale protocol for the structure prediction of linear and cyclic peptides. The protocol combines two main stages: coarse-grained simulations using the CABS-flex standalone package and an all-atom reconstruction-optimization process using the Modeller program. We evaluated the protocol on a set of linear peptides and two sets of cyclic peptides, with cyclization through the backbone and disulfide bonds. A comparison with other state-of-the-art tools (APPTEST, PEP-FOLD, ESMFold and AlphaFold implementation in ColabFold) shows that for most cases, AlphaFold offers the highest resolution. However, CABS-flex is competitive, particularly when it comes to short linear peptides. As demonstrated, the protocol performance can be further improved by combination with the residue-residue contact prediction method or more efficient scoring. The protocol is included in the CABS-flex standalone package along with online documentation to aid users in predicting the structure of peptides and mini-proteins.

A3D Model Organism Database (A3D-MODB): a database for proteome aggregation predictions in model organisms

Protein aggregation has been associated with aging and different pathologies and represents a bottleneck in the industrial production of biotherapeutics. Numerous past studies performed in Escherichia coli and other model organisms have allowed to dissect the biophysical principles underlying this process. This knowledge fuelled the development of computational tools, such as Aggrescan 3D (A3D) to forecast and re-design protein aggregation. Here, we present the A3D Model Organism Database (A3D-MODB) http://biocomp.chem.uw.edu.pl/A3D2/MODB, a comprehensive resource for the study of structural protein aggregation in the proteomes of 12 key model species spanning distant biological clades. In addition to A3D predictions, this resource incorporates information useful for contextualizing protein aggregation, including membrane protein topology and structural model confidence, as an indirect reporter of protein disorder. The database is openly accessible without any need for registration. We foresee A3D-MOBD evolving into a central hub for conducting comprehensive, multi-species analyses of protein aggregation, fostering the development of protein-based solutions for medical, biotechnological, agricultural and industrial applications.

Characterization of exclusively non-commensal Neisseria gonorrhoeae pangenome to prioritize globally conserved and thermodynamically stable vaccine candidates using immune-molecular dynamic simulations

Neisseria gonorrhoeae (Ngo) has emerged as a global threat leading to one of the most common sexually transmitted diseases in the world. It has also become one of the leading antimicrobial resistant organisms, resulting in fewer treatment options and an increased morbidity. Therefore, in recent years, there has been an increased focus on the development of new treatments and preventive strategies to combat its infection. In this study, we have combined the most conserved epitopes from the completely assembled strains of Ngo to develop a universal and a thermodynamically stable vaccine candidate. For our vaccine design, the epitopes were selected for their high immunogenicity, non-allergenicity and non-cytotoxicity, making them the ideal candidates for vaccine development. For the screening process, several reverse vaccinology tools were employed to rigorously extract non-homologous and immunogenic epitopes from the selected proteins. Consequently, a total number of 3 B-cell epitopes and 6 T-cell epitopes were selected and joined by multiple immune-modulating adjuvants and linkers to generate a promiscuous immune response. Additionally, the stability and flexible nature of the vaccine construct was confirmed using various molecular dynamic simulation tools. Overall, the vaccine candidate showed promising binding affinity to various HLA alleles and TLR receptors; however, further studies are needed to assess its efficacy in-vivo. In this way, we have designed a multi-subunit vaccine candidate to potentially combat and control the spread of N. gonorrhoeae.

Oncoinformatic screening of the gene clusters involved in the HER2-positive breast cancer formation along with the in silico pharmacodynamic profiling of selective long-chain omega-3 fatty acids as the metastatic antagonists

The HER2-positive patients occupy ~ 30% of the total breast cancer patients globally where no prevalent drugs are available to mitigate the frequent metastasis clinically except lapatinib and neratinib. This scarcity reinforced researchers' quest for new medications where natural substances are significantly considered. Valuing the aforementioned issues, this research aimed to study the ERBB2-mediated string networks that work behind the HER2-positive breast cancer formation regarding co-expression, gene regulation, GAMA-receptor-signaling pathway, cellular polarization, and signal inhibition. Following the overexpression, promotor methylation, and survivability profiles of ERBB2, the super docking position of HER2 was identified using the quantum tunneling algorithm. Supramolecular docking was conducted to study the target specificity of EPA and DHA fatty acids followed by a comprehensive molecular dynamic simulation (100 ns) to reveal the RMSD, RMSF, Rg, SASA, H-bonds, and MM/GBSA values. Finally, potential drug targets for EPA and DHA in breast cancer were constructed to determine the drug-protein interactions (DPI) at metabolic stages. Considering the values resulting from the combinational models of the oncoinformatic, pharmacodynamic, and metabolic parameters, long-chain omega-3 fatty acids like EPA and DHA can be considered as potential-targeted therapeutics for HER2-positive breast cancer treatment.

Nanobody-peptide-conjugate (NPC) for passive immunotherapy against SARS-CoV-2 variants of concern (VoC): a prospective pan-coronavirus therapeutics

The COVID-19 crisis, incited by the zoonotic SARS-CoV-2 virus, has quickly escalated into a catastrophic public health issue and a grave threat to humankind owing to the advent of mutant viruses. Multiple pharmaceutical therapies or biologics envision stopping the virus from spreading further; however, WHO has voiced concerns about the variants of concern (VoCs) inability to respond. Nanobodies are a new class of antibody mimics with binding affinity and specificity similar to classical mAbs, as well as the privileges of a small molecular weight, ease of entry into solid tissues, and binding cryptic epitopes of the antigen. Herein, we investigated multiple putative anti-SARS-CoV-2 nanobodies targeting the Receptor binding domain of the WHO-listed SARS-CoV-2 variants of concern using a comprehensive computational immunoinformatics methodology. With affinity maturation via alanine scanning mutagenesis, we remodeled an ultrapotent nanobody with substantial breadth and potency, exhibiting pico-molar binding affinities against all the VoCs. An antiviral peptide with specificity for ACE-2 receptors was affixed to make it multispecific and discourage viral entry. Collectively, we constructed a broad-spectrum therapeutic biparatopic nanobody-peptide conjugate (NPC) extending coverage to SARS-CoV-2 VoCs RBDs. We PEGylated the biparatopic construct with 20kD maleimide-terminated PEG (MAL-(PEG)n-OMe) to improve its clinical efficacy limiting rapid renal clearance, and performed in silico cloning to facilitate future experimental studies. Our findings suggest that combining biparatopic nanobody conjugate with standard treatment may be a promising bivariate tool for combating viral entry during COVID-19 illness.

Molecular docking-based virtual screening, drug-likeness, and pharmacokinetic profiling of some anti-Salmonella typhimurium cephalosporin derivatives

Objective

The rising cases of resistance to existing antibiotic therapies in Salmonella typhimurium has made it necessary to search for novel drug candidates. The present study employed the molecular docking technique to screen a set of antibacterial cephalosporin analogues against penicillin-binding protein 1a (PBP1a) of the bacterium. This is the first study to screen cephalosporin analogues against PBP1a, a protein central to peptidoglycan synthesis in S. typhimurium.

Methods

Some cephalosporin analogues were retrieved from a drug repository. The structures of the molecules were optimized using the semi-empirical method of Spartan 14 software and were subsequently docked against the active sites of PBP1a using AutoDock vina software. The most potent ligands were chosen as the most promising leads and subsequently subjected to absorption, distribution, metabolism, excretion, and toxicity (ADMET) profiling using the SwissADME online server and DataWarrior chemoinformatics program. The CABSflex 2.0 server was used to carry out molecular dynamics (MD) simulation on the most stable ligand-protein complex.

Results

Compounds 3, 23, and 28 with binding affinity (ΔG) values of -9.2, -8.7, and -8.9 kcal/mol, respectively, were selected as the most promising leads. The ligands bound to the active sites of PBP1a via hydrophobic bonds, hydrogen bonds, and electrostatic interactions. Furthermore, ADMET analyses of the ligands revealed that they exhibited sound pharmacokinetic and toxicity profiles. In addition, an MD study revealed that the most active ligand bound favorably and dynamically to the target protein.

Conclusion

The findings of this research could provide an excellent platform for the discovery and rational design of novel antibiotics against S. typhimurium. Additional in vitro and in vivo studies should be carried out on the drug candidates to validate the findings of this study.

Proteome based analysis of circulating SARS-CoV-2 variants: approach to a universal vaccine candidate

The discovery of the first infectious variant in Wuhan, China, in December 2019, has posed concerns over global health due to the spread of COVID-19 and subsequent variants. While the majority of patients experience flu-like symptoms such as cold and fever, a small percentage, particularly those with compromised immune systems, progress from mild illness to fatality. COVID-19 is caused by a RNA virus known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Our approach involved utilizing immunoinformatic to identify vaccine candidates with multiple epitopes and ligand-binding regions in reported SARS-CoV-2 variants. Through analysis of the spike glycoprotein, we identified dominant epitopes for T-cells and B-cells, resulting in a vaccine construct containing two helper T-cell epitopes, six cytotoxic T-cell epitopes, and four linear B-cell epitopes. Prior to conjugation with adjuvants and linkers, all epitopes were evaluated for antigenicity, toxicity, and allergenicity. Additionally, we assessed the vaccine Toll-Like Receptors complex (2, 3, and 4). The vaccine construct demonstrated antigenicity, non-toxicity, and non-allergenicity, thereby enabling the host to generate antibodies with favorable physicochemical characteristics. Furthermore, the 3D structure of the B-cell construct exhibited a ProSA-web z-score plot with a value of -1.71, indicating the reliability of the designed structure. The Ramachandran plot analysis revealed that 99.6% of the amino acid residues in the vaccine subunit were located in the high favored observation region, further establishing its strong candidacy as a vaccination option.

Novel COX11 Mutations Associated with Mitochondrial Disorder: Functional Characterization in Patient Fibroblasts and Saccharomyces cerevisiae

Genetic defects in the nuclear encoded subunits and assembly factors of cytochrome c oxidase (mitochondrial complex IV) are very rare and are associated with a wide variety of phenotypes. Biallelic pathogenic variants in the COX11 protein were previously identified in two unrelated children with infantile-onset mitochondrial encephalopathies. Through comprehensive clinical, genetic and functional analyses, here we report on a new patient harboring novel heterozygous variants in COX11, presenting with Leigh-like features, and provide additional experimental evidence for a direct correlation between COX11 protein expression and sensitivity to oxidative stress. To sort out the contribution of the single mutations to the phenotype, we employed a multi-faceted approach using Saccharomyces cerevisiae as a genetically manipulable system, and in silico structure-based analysis of human COX11. Our results reveal differential effects of the two novel COX11 mutations on yeast growth, respiration, and cellular redox status, as well as their potential impact on human protein stability and function. Strikingly, the functional deficits observed in patient fibroblasts are recapitulated in yeast models, validating the conservation of COX11's role in mitochondrial integrity across evolutionarily distant organisms. This study not only expands the mutational landscape of COX11-associated mitochondrial disorders but also underscores the continued translational relevance of yeast models in dissecting complex molecular pathways.

AI-Based Homology Modelling of Fatty Acid Transport Protein 1 Using AlphaFold: Structural Elucidation and Molecular Dynamics Exploration

Fatty acid transport protein 1 (FATP1) is an integral transmembrane protein that is involved in facilitating the translocation of long-chain fatty acids (LCFA) across the plasma membrane, thereby orchestrating the importation of LCFA into the cell. FATP1 also functions as an acyl-CoA ligase, catalyzing the ATP-dependent formation of fatty acyl-CoA using LCFA and VLCFA (very-long-chain fatty acids) as substrates. It is expressed in various types of tissues and is involved in the regulation of crucial signalling pathways, thus playing a vital role in numerous physiological and pathological conditions. Structural insight about FATP1 is, thus, extremely important for understanding the mechanism of action of this protein and developing efficient treatments against its anomalous expression and dysregulation, which are often associated with pathological conditions such as breast cancer. As of now, there has been no prior prediction or evaluation of the 3D configuration of the human FATP1 protein, hindering a comprehensive understanding of the distinct functional roles of its individual domains. In our pursuit to unravel the structure of the most commonly expressed isoforms of FATP1, we employed the cutting-edge ALPHAFOLD 2 model for an initial prediction of the entire protein's structure. This prediction was complemented by molecular dynamics simulations, focusing on the most promising model. We predicted the structure of FATP1 in silico and thoroughly refined and validated it using coarse and molecular dynamics in the absence of the complete crystal structure. Their relative dynamics revealed the different properties of the characteristic FATP1.

First 3-D structural evidence of a native-like intertwined dimer in the acylphosphatase family

Acylphosphatase (AcP, EC 3.6.1.7) is a small model protein conformed by a ferredoxin-like fold, profoundly studied to get insights into protein folding and aggregation processes. Numerous studies focused on the aggregation and/or amyloidogenic properties of AcPs suggest the importance of edge-β-strands in the process. In this work, we present the first crystallographic structure of Escherichia coli AcP (EcoAcP), showing notable differences with the only available NMR structure for this enzyme. EcoAcP is crystalised as an intertwined dimer formed by replacing a single C-terminal β-strand between two protomers, suggesting a flexible character of the C-terminal edge of EcoAcP. Despite numerous works where AcP from different sources have been used as a model system for protein aggregation, our domain-swapped EcoAcP structure is the first 3-D structural evidence of native-like aggregated species for any AcP reported to date, providing clues on molecular determinants unleashing aggregation.

Prediction of Protein Aggregation Propensity via Data-Driven Approaches

Protein aggregation occurs when misfolded or unfolded proteins physically bind together and can promote the development of various amyloid diseases. This study aimed to construct surrogate models for predicting protein aggregation via data-driven methods using two types of databases. First, an aggregation propensity score database was constructed by calculating the scores for protein structures in the Protein Data Bank using Aggrescan3D 2.0. Moreover, feature- and graph-based models for predicting protein aggregation have been developed by using this database. The graph-based model outperformed the feature-based model, resulting in an R2 of 0.95, although it intrinsically required protein structures. Second, for the experimental data, a feature-based model was built using the Curated Protein Aggregation Database 2.0 to predict the aggregated intensity curves. In summary, this study suggests approaches that are more effective in predicting protein aggregation, depending on the type of descriptor and the database.

Multi-epitope vaccine candidates based on mycobacterial membrane protein large (MmpL) proteins against Mycobacterium ulcerans

Buruli ulcer (BU) is a neglected tropical disease. It is caused by the bacterium Mycobacterium ulcerans and is characterized by skin lesions. Several studies were performed testing the Bacillus Calmette-Guérin (BCG) vaccine in human and animal models and M. ulcerans-specific vaccines in animal models. However, there are currently no clinically accepted vaccines to prevent M. ulcerans infection. The aim of this study was to identify T-cell and B-cell epitopes from the mycobacterial membrane protein large (MmpL) proteins of M. ulcerans. These epitopes were analysed for properties including antigenicity, immunogenicity, non-allergenicity, non-toxicity, population coverage and the potential to induce cytokines. The final 8 CD8+, 12 CD4+ T-cell and 5 B-cell epitopes were antigenic, non-allergenic and non-toxic. The estimated global population coverage of the CD8+ and CD4+ epitopes was 97.71%. These epitopes were used to construct five multi-epitope vaccine constructs with different adjuvants and linker combinations. The constructs underwent further structural analyses and refinement. The constructs were then docked with Toll-like receptors. Three of the successfully docked complexes were structurally analysed. Two of the docked complexes successfully underwent molecular dynamics simulations (MDS) and post-MDS analysis. The complexes generated were found to be stable. However, experimental validation of the complexes is required.

In silico screening of herbal phytochemicals to develop a Rasayana for immunity against Nipah virus

BACKGROUND

The first emergence of the Nipah virus (NiV) in 1998 from Malaysia became a major concern when it came to light and resurfaced on different occasions thereafter. NiV is a bat-borne zoonotic and pleomorphic virus that causes severe infection in human and animal hosts. Studies revealed fruit bats are the major reservoirs as natural hosts and pigs as intermediate hosts for the spread of this infection. This became a major concern as the disease was characterized by high pathogenicity varying from 40% to 80% depending on its acuteness. Moreover, the solemnity lies in the fact that the infection transcends from being a mere mild illness to an acute respiratory infection leading to fatal encephalitis with a reportedly high mortality rate. Currently, there is no treatment or vaccine available against the NiV. Many antiviral drugs have been explored and developed but with limited efficacy.

METHODOLOGY

In search of high-affinity ayurvedic alternatives, we conducted a pan-proteome in silico exploration of the NiV proteins for their interaction with the best-suited phytoconstituents. The toxicity prediction of thirty phytochemicals based on their LD50 value identified thirteen potential candidates. Molecular docking studies of those thirteen phytochemicals with five important NiV proteins identified Tanshinone I as the potential compound with a high binding affinity.

RESULTS

The pharmacokinetics and pharmacodynamics studies also aided in determining the absorption, distribution, metabolism, excretion, and toxicity of the selected phytoconstituent. Interestingly, docking studies also revealed Rosmariquinone as a potent alternative to the antiviral drug Remdesivir binding the same pocket of RNA-dependent RNA polymerase of the NiV. A molecular dynamics simulation study of the surface glycoprotein of NiV against Tanshinone I showed a stable complex formation and significant allosteric changes in the protein structure, implying that these phytochemicals could be a natural alternative to synthetic drugs against NiV.

CONCLUSION

This study provides preliminary evidence based on in silico analysis that the herbal molecules showed an effect against NiV. However, it is essential to further evaluate the efficacy of this approach through cell-based experiments, organoid models, and eventually clinical trials.

Anti-COVID-19, Anti-Inflammatory, and Anti-Osteoarthritis Activities of Sesamin from Sesamum indicum L

During the COVID-19 (coronavirus disease 2019) outbreak, many people were infected, and the symptoms may persist for several weeks or months for recovering patients. This is also known as "long COVID" and includes symptoms such as fatigue, joint pain, muscle pain, et cetera. The COVID-19 virus may trigger hyper-inflammation associated with cytokine levels in the body. COVID-19 can trigger inflammation in the joints, which can lead to osteoarthritis (OA), while long-term COVID-19 symptoms may lead to joint damage and other inflammation problems. According to several studies, sesame has potent anti-inflammatory properties due to its major constituent, sesamin. This study examined sesamin's anti-inflammatory, anti-osteoarthritis, and anti-COVID-19 effects. Moreover, in vivo and in vitro assays were used to determine sesamin's anti-inflammatory activity against the RAW264.7 and SW1353 cell lines. Sesamin had a dose-dependent effect (20 mg/kg) in a monoiodoacetic acid (MIA)-induced osteoarthritis rat model. Sesamin reduced paw swelling and joint discomfort. In addition, the findings indicated that sesamin suppressed the expression of iNOS (inducible nitric oxide synthase) and COX-2 (cyclooxygenase-2) in the RAW264.7 cell line within the concentration range of 6.25-50 μM. Furthermore, sesamin also had a suppressive effect on MMP (matrix metalloproteinase) expression in chondrocytes and the SW1353 cell line within the same concentration range of 6.25-50 μM. To examine the anti-viral activity, an in silico analysis was performed to evaluate sesamin's binding affinity with SARS-CoV-2 RdRp (severe acute respiratory syndrome coronavirus 2 RNA-dependent RNA polymerase) and human ACE2 (angiotensin-converting enzyme 2). Compared to the controls, sesamin exhibited strong binding affinities towards SARS-CoV-2 RdRp and human ACE2. Furthermore, sesamin had a higher binding affinity for the ACE2 target protein. This study suggests that sesamin shows potential anti-SARS-CoV-2 activity for drug development.

Immunoinformatics aided designing of a next generation poly-epitope vaccine against uropathogenic Escherichia coli to combat urinary tract infections

Urinary tract infections (UTIs) are the second most prevalent bacterial infections and uropathogenic Escherichia coli (UPEC) stands among the primary causative agents of UTIs. The usage of antibiotics is the routine therapy being used in various countries to treat UTIs but becoming ineffective because of increasing antibiotic resistance among UPEC strains. Thus, there must be the development of some alternative treatment strategies such as vaccine development against UPEC. In the following study, pan-genomics along with reverse vaccinology approaches is used under the framework of bioinformatics for the identification of core putative vaccine candidates, employing 307 UPEC genomes (complete and draft), available publicly. A total of nine T-cell epitopes (derived from B-cells) of both MHC classes (I and II), were prioritized among three potential protein candidates. These epitopes were then docked together by using linkers (GPGPG and AAY) and an adjuvant (Cholera Toxin B) to form a poly-valent vaccine construct. The chimeric vaccine construct was undergone by molecular modelling, further refinement and energy minimization. We predicted positive results of the vaccine construct in immune simulations with significantly high levels of immune cells. The protein-protein docking analysis of vaccine construct with toll-like receptors predicted efficient binding, which was further validated by molecular dynamics simulation of vaccine construct with TLR-2 and TLR-4 at 120 ns, resulting in stable complexes' conformation throughout the simulation run. Overall, the vaccine construct demonstrated positive antigenic response. In future, this chimeric vaccine construct or the identified epitopes could be experimentally validated for the development of UPEC vaccines against UTIs.Communicated by Ramaswamy H. Sarma.

Modulation of Neuron and Astrocyte Dopamine Receptors via Receptor-Receptor Interactions

Dopamine neurotransmission plays critical roles in regulating complex cognitive and behavioral processes including reward, motivation, reinforcement learning, and movement. Dopamine receptors are classified into five subtypes, widely distributed across the brain, including regions responsible for motor functions and specific areas related to cognitive and emotional functions. Dopamine also acts on astrocytes, which express dopamine receptors as well. The discovery of direct receptor-receptor interactions, leading to the formation of multimeric receptor complexes at the cell membrane and providing the cell decoding apparatus with flexible dynamics in terms of recognition and signal transduction, has expanded the knowledge of the G-protein-coupled receptor-mediated signaling processes. The purpose of this review article is to provide an overview of currently identified receptor complexes containing dopamine receptors and of their modulatory action on dopamine-mediated signaling between neurons and between neurons and astrocytes. Pharmacological possibilities offered by targeting receptor complexes in terms of addressing neuropsychiatric disorders associated with altered dopamine signaling will also be briefly discussed.

Identification of Potential Phytochemical Inhibitors From Conium maculatum Targeting the Epidermal Growth Factor Receptor in Metastatic Colorectal Cancer via Molecular Docking Analysis

Background Metastatic colorectal cancer (mCRC) continues to rank as the second deadliest cancer on the global scale. CRC diagnosed at metastatic (stage IV) makes treatment strategies more challenging. Even though there are numerous therapeutic options available, the side effects of these treatments threaten the human health. Therefore, we are in the phase of searching new molecules that are less harmful and cost-effective. The common source of many pharmaceutical medications is plants. This study focuses on virtually screening phytochemicals from Conium maculatum as potential inhibitors of the epidermal growth factor receptor (EGFR), a crucial target in cancer therapy. Methods and materials C. maculatum was selected due to its phytochemicals and prior indications of its anticancer properties. In silico investigations encompass druglikeness screening, pharmacokinetics assessment, molecular docking, toxicity prediction, molecular target screening, and molecular dynamics simulations. A comprehensive analysis led to the identification of promising lead compounds. Results A total of 25 compounds exhibited favorable pharmacokinetic and drug-like characteristics. Among them, 12 compounds displayed a high affinity for EGFR as determined by molecular docking experiments. Further safety assessment using ProTox-II revealed that seven compounds had no anticipated toxicity, affirming their safety profiles.  Conclusion These findings collectively predicted the efficacy of seven phytochemicals from C. maculatum as EGFR inhibitors in mCRC. Further experimental investigations and optimization of the identified leads were needed to validate the efficacy and safety of identified lead compounds and explore their therapeutic potential in CRC.

Disordered regions endow structural flexibility to shell proteins and function towards shell-enzyme interactions in 1,2-propanediol utilization microcompartment

Intrinsically disordered regions in proteins have been functionally linked to the protein-protein interactions and genesis of several membraneless organelles. Depending on their residual makeup, hydrophobicity or charge distribution they may remain in extended form or may assume certain conformations upon biding to a partner protein or peptide. The present work investigates the distribution and potential roles of disordered regions in the integral proteins of 1,2-propanediol utilization microcompartments. We use bioinformatics tools to identify the probable disordered regions in the shell proteins and enzyme of the 1,2-propanediol utilization microcompartment. Using a combination of computational modelling and biochemical techniques we elucidate the role of disordered terminal regions of a major shell protein and enzyme. Our findings throw light on the importance of disordered regions in the self-assembly, providing flexibility to shell protein and mediating its interaction with a native enzyme.Communicated by Ramaswamy H. Sarma.

Investigation of the impact of nonsynonymous mutations on thyroid peroxidase dimer

Congenital hypothyroidism is one of the most common preventable endocrine disorders associated with thyroid dysgenesis or dyshormonogenesis. Thyroid peroxidase (TPO) gene defect is mainly responsible for dyshormonogenesis; a defect in the thyroid hormone biosynthesis pathway. In Bangladesh, there is limited data regarding the genetic etiology of Congenital Hypothyroidism (CH). The present study investigates the impact of the detected mutations (p.Ala373Ser, and p.Thr725Pro) on the TPO dimer protein. We have performed sequential molecular docking of H2O2 and I- ligands with both monomers of TPO dimer to understand the iodination process in thyroid hormone biosynthesis. Understanding homodimer interactions at the atomic level is a critical challenge to elucidate their biological mechanisms of action. The docking results reveal that mutations in the dimer severely disrupt its catalytic interaction with essential ligands. Molecular dynamics simulation has been performed to validate the docking results, thus realizing the consequence of the mutation in the biological system's mimic. The dynamics results expose that mutations destabilize the TPO dimer protein. Finally, principal component analysis exhibits structural and energy profile discrepancies in wild-type and mutant dimers. The findings of this study highlight that the mutations in TPO protein can critically affect the dimer structure and loss of enzymatic activity is persistent. Other factors also might influence the hormone synthesis pathway, which is under investigation.

In silico Analysis of Two Novel Variants in the Pyruvate Carboxylase (PC) Gene Associated with the Severe Form of PC Deficiency

Background

Inborne errors of metabolism are a common cause of neonatal death. This study evaluated the acute early-onset metabolic derangement and death in two unrelated neonates.

Methods

Whole-exome sequencing (WES), Sanger sequencing, homology modeling, and in silico bioinformatics analysis were employed to assess the effects of variants on protein structure and function.

Results

WES revealed a novel homozygous variant, p.G303Afs*40 and p.R156P, in the pyruvate carboxylase (PC) gene of each neonate, which both were confirmed by Sanger sequencing. Based on the American College of Medical Genetics and Genomics guidelines, the p.G303Afs*40 was likely pathogenic, and the p.R156P was a variant of uncertain significance (VUS). Nevertheless, a known variant at position 156, the p.R156Q, was also a VUS. Protein secondary structure prediction showed changes in p.R156P and p.R156Q variants compared to the wild-type protein. However, p.G303Afs*40 depicted significant changes at C-terminal. Furthermore, comparing the interaction of wild-type and variant proteins with the ATP ligand during simulations, revealed a decreased affinity to the ATP in all the variants. Moreover, analysis of Single nucleotide polymorphism impacts on PC protein using Polyphen-2, SNAP2, FATHMM, and SNPs&GO servers predicted both R156P and R156Q as damaging variants. Likewise, free energy calculations demonstrated the destabilizing effect of both variants on PC.

Conclusion

This study confirmed the pathogenicity of both variants and suggested them as a cause of type B Pyruvate carboxylase deficiency. The results of this study would provide the family with prenatal diagnosis and expand the variant spectrum in the PC gene,which is beneficial for geneticists and endocrinologists.

In silico and immunoinformatics based multiepitope subunit vaccine design for protection against visceral leishmaniasis

Visceral leishmaniasis (VL) is a vector-borne neglected tropical protozoan disease with high fatality and no certified vaccine. Conventional vaccine preparation is challenging and tedious. Here in this work, we created a global multiepitope subunit vaccination against VL utilizing innovative immunoinformatics technique based on the extensively conserved epitopic regions of the PrimPol protein of Leishmania donovani consisting of four subunits which were analyzed and studied, out of which DNA primase large subunit and DNA polymerase α subunit B were evaluated as antigens by Vaxijen 2.0. The multiepitope vaccine design includes a single adjuvant β-defensins, eight CTL epitopes, eight HTL epitopes, seven linear BCL epitopes and one discontinuous BCL epitope to induce innate, cellular and humoral immune responses against VL. The Expasy ProtParam tool characterized the physiochemical parameters of the vaccine. At the same time, SOLpro evaluated our vaccine constructs to be soluble upon expression. We also modeled the stable tertiary structure of our vaccine construct through Robetta modeling for molecular docking studies with toll-like receptor proteins through HADDOCK 2.4. Simulations based on molecular dynamics revealed an intact vaccine and TLR8 complex, supporting our vaccine design's immunogenicity. Also, the immune simulation of our vaccine by the C-ImmSim server demonstrated the potency of the multiepitope vaccine construct to induce proper immune response for host defense. Codon optimization and in silico cloning of our vaccine further assured high expression. The outcomes of our study on multiepitope vaccine design significantly produced a potential candidate against VL and can potentially eradicate the disease in the future after clinical investigations.Communicated by Ramaswamy H. Sarma.

Computational Identification of Most Deleterious Missense Mutations in Human PD-1 Gene

Traditional cancer treatment approaches are often hindered by the presence of toxic side effects and the high rate of relapse observed in treated organs. In contrast, novel immunotherapeutic strategies targeting immune checkpoint inhibitors, particularly PD-1, have demonstrated promising results with minimal adverse effects. However, the emergence of immunotherapeutic-resistant tumors, predominantly caused by intrinsic mutations, poses a significant obstacle to successful treatment outcomes. Consequently, the primary objective of this study was to screen for the most detrimental missense mutations in the PD-1 gene associated with immunotherapeutic resistance. To achieve this aim, a comprehensive screening process utilizing 20 web servers, incorporating both sequence- and structure-based methodologies, was undertaken. Through meticulous analysis and mutual disease association sorting, four specific missense mutations were successfully identified. These mutations, namely, R38C, D61V, R94C, and D117V, emerged as the leading contributors to genetic cancer progression and immunotherapeutic resistance against PD-1 blockers. The findings presented in this study are supported by multiple lines of evidence. A thorough examination of protein topology, structural alignment, docking interactions with PD-L1, and protein flexibility collectively confirmed the pathogenic nature of these sorted mutations. By considering these various aspects, we have gained a comprehensive understanding of the underlying mechanisms driving immunotherapeutic resistance. In conclusion, the comprehensive screening process undertaken in this study has successfully identified R38C, D61V, R94C, and D117V as the primary mutations contributing to genetic cancer progression and immunotherapeutic resistance against PD-1 blockers. The integration of protein topology analysis, structural alignment, docking studies with PD-L1, and assessment of protein flexibility have collectively provided robust evidence to support the pathogenic significance of these mutations.

POU6F2 mutation in humans with pubertal failure alters GnRH transcript expression

Idiopathic hypogonadotropic hypogonadism (IHH) is characterized by the absence of pubertal development and subsequent impaired fertility often due to gonadotropin-releasing hormone (GnRH) deficits. Exome sequencing of two independent cohorts of IHH patients identified 12 rare missense variants in POU6F2 in 15 patients. POU6F2 encodes two distinct isoforms. In the adult mouse, expression of both isoform1 and isoform2 was detected in the brain, pituitary, and gonads. However, only isoform1 was detected in mouse primary GnRH cells and three immortalized GnRH cell lines, two mouse and one human. To date, the function of isoform2 has been verified as a transcription factor, while the function of isoform1 has been unknown. In the present report, bioinformatics and cell assays on a human-derived GnRH cell line reveal a novel function for isoform1, demonstrating it can act as a transcriptional regulator, decreasing GNRH1 expression. In addition, the impact of the two most prevalent POU6F2 variants, identified in five IHH patients, that were located at/or close to the DNA-binding domain was examined. Notably, one of these mutations prevented the repression of GnRH transcripts by isoform1. Normally, GnRH transcription increases as GnRH cells mature as they near migrate into the brain. Augmentation earlier during development can disrupt normal GnRH cell migration, consistent with some POU6F2 variants contributing to the IHH pathogenesis.

Understanding the Molecular Regulation of Serotonin Receptor 5-HTR1B-β-Arrestin1 Complex in Stress and Anxiety Disorders

The serotonin receptor subtype 5-HTR1B is widely distributed in the brain with an important role in various behavioral implications including neurological conditions and psychiatric disorders. The neuromodulatory action of 5-HTR1B largely depends upon its arrestin mediated signaling pathway. In this study, we tried to investigate the role of unusually long intracellular loop 3 (ICL3) region of the serotonin receptor 5-HTR1B in interaction with β-arrestin1 (Arr2) to compensate for the absence of the long cytoplasmic tail. Molecular modeling and docking tools were employed to obtain a suitable molecular conformation of the ICL3 region in complex with Arr2 which dictates the specific complex formation of 5-HTR1B with Arr2. This reveals the novel molecular mechanism of phosphorylated ICL3 mediated GPCR-arrestin interaction in the absence of the long cytoplasmic tail. The in-cell disulfide cross-linking experiments and molecular dynamics simulations of the complex further validate the model of 5-HTR1B-ICL3-Arr2 complex. Two serine residues (Ser281 and Ser295) within the 5-HTR1B-ICL3 region were found to be occupying the electropositive pocket of Arr2 in our model and might be crucial for phosphorylation and specific Arr2 binding. The alignment studies of these residues showed them to be conserved only across 5-HTR1B mammalian species. Thus, our studies were able to predict a molecular conformation of 5-HTR1B-Arr2 and identify the role of long ICL3 in the signaling process which might be crucial in designing targeted drugs (biased agonists) that promote GPCR-Arr2 signaling to deter the effects of stress and anxiety-like disorders.

Membrane Cholesterol Inhibits Progesterone-Mediated Sperm Function through the Possible Involvement of ABHD2

Abhydrolase domain containing 2-acylglycerol lipase (ABHD2) was recently claimed as the membrane receptor of progesterone (P4) in sperm cells, mediating cell processes such as sperm chemotaxis and acrosome reaction. Here, we investigated the role of membrane cholesterol (Chol) on ABHD2-mediated human sperm chemotaxis. Human sperm cells were obtained from twelve normozoospemic healthy donors. ABHD2-Chol interaction was modelled by computational molecular-modelling (MM). Sperm membrane Chol content was depleted by incubating cells with cyclodextrin (CD) or augmented by the incubation with the complex between CD and Chol (CD:Chol). Cell Chol levels were quantified by liquid chromatography-mass spectrometry. Sperm migration upon P4 gradient was evaluated through the accumulation assay in a specific migration device. Motility parameters were evaluated by sperm class analyzer, whilst intracellular calcium concentration, acrosome reaction and mitochondrial membrane potential were evaluated with calcium orange, FITC-conjugated anti-CD46 antibody and JC-1 fluorescent probes, respectively. MM analysis showed the possible stable binding Chol to ABHD2, resulting in to major impact on the protein backbone flexibility. The treatment with CD was associated with a dose-dependent increase in sperm migration in a 160 nM P4 gradient, together with increase in sperm motility parameters and levels of acrosome reaction. The treatment with CD:Chol was associated with essentially opposite effects. Chol was, thus, suggested to inhibit P4-mediated sperm function through the possible inhibition of ABHD2.

Attenuation of fluoride-induced hepatorenal oxidative stress by ferulic acid in vivo: An approach with in-silico analysis and interaction informatics of ferulic acid

BACKGROUND

Chronic fluoride toxicity induces oxidative strain and lipid peroxidation and imparts deleterious effects on human metabolic organs.

AIM

The present study aimed to expose the defensive impact of ferulic acid against sodium fluoride (NaF) induced hepatorenal dysfunction at the biochemical and antioxidative systems.

METHODS

In-vivo. Rats were arbitrarily separated into five groups as control, sodium fluoride-treated (200 ppm kg ^-1), vitamin C -as a positive control, and FA co-administered groups with 10 mg kg ^-1 and 20 mg kg ^-1 body weight for 56 days. In the present investigation, we measured antioxidant enzymes, superoxide dismutase, catalase, and lactate dehydrogenase by electrozymographic and spectrophotometric methods. Biochemical assessment of TBARS, conjugated diene, and different serum biomarkers was done for liver and kidney functionality tests. In-silico. An in-silico study was conducted through a molecular docking experiment to evaluate the binding potentiality of FA by employing AutoDock Vina [version 1.5.6] to overcome the abnormality in the activities of catalase, and superoxide dismutase in NaF promoted toxicity of hepatorenal system. In-vitro. An in vitro biochemical experiment was conducted to support the in-silico study.

RESULTS

Superoxide dismutase and catalase were decreased in the intoxicated rat. Ferulic acid (FA) as an antioxidant remarkably defended the NaF-mediated deterioration of the antioxidative status in the hepatorenal system, lowering lipid peroxidation products, malondialdehyde, and conjugated diene. Serum biomarkers, ALT, AST, ALP, urea, and creatinine increased in the intoxicated group than in control. Ferulic acid significantly neutralized the ill effects of NaF on serum lipid profile. In-silico analysis hypothesized the strong interaction of FA with the active side of catalase and superoxide dismutase that prevented the binding of NaF at the active site of these mentioned enzymes and this was further validated by in-vitro assay.

CONCLUSION

However, FA modulates free radical generation and protected these metabolic organs against sodium fluoride-induced injury.

TRAP1 S-nitrosylation as a model of population-shift mechanism to study the effects of nitric oxide on redox-sensitive oncoproteins

S-nitrosylation is a post-translational modification in which nitric oxide (NO) binds to the thiol group of cysteine, generating an S-nitrosothiol (SNO) adduct. S-nitrosylation has different physiological roles, and its alteration has also been linked to a growing list of pathologies, including cancer. SNO can affect the function and stability of different proteins, such as the mitochondrial chaperone TRAP1. Interestingly, the SNO site (C501) of TRAP1 is in the proximity of another cysteine (C527). This feature suggests that the S-nitrosylated C501 could engage in a disulfide bridge with C527 in TRAP1, resembling the well-known ability of S-nitrosylated cysteines to resolve in disulfide bridge with vicinal cysteines. We used enhanced sampling simulations and in-vitro biochemical assays to address the structural mechanisms induced by TRAP1 S-nitrosylation. We showed that the SNO site induces conformational changes in the proximal cysteine and favors conformations suitable for disulfide bridge formation. We explored 4172 known S-nitrosylated proteins using high-throughput structural analyses. Furthermore, we used a coarse-grained model for 44 protein targets to account for protein flexibility. This resulted in the identification of up to 1248 proximal cysteines, which could sense the redox state of the SNO site, opening new perspectives on the biological effects of redox switches. In addition, we devised two bioinformatic workflows ( https://github.com/ELELAB/SNO_investigation_pipelines ) to identify proximal or vicinal cysteines for a SNO site with accompanying structural annotations. Finally, we analyzed mutations in tumor suppressors or oncogenes in connection with the conformational switch induced by S-nitrosylation. We classified the variants as neutral, stabilizing, or destabilizing for the propensity to be S-nitrosylated and undergo the population-shift mechanism. The methods applied here provide a comprehensive toolkit for future high-throughput studies of new protein candidates, variant classification, and a rich data source for the research community in the NO field.

Deciphering Houttuynia cordata extract as electron shuttles with anti-COVID-19 activity and its performance in microbial fuel cells

Background

Traditional herbal medicines usually contain electron shuttle (ES)-like structures compounds which are potential candidates for antiviral compounds selection. Houttuynia cordata is applied as a biomaterial to decipher its potential applications in bioenergy extraction in microbial fuel cells (MFCs) and anti-COVID-19 via molecular docking evaluation.

Methods

H. cordata leaves extracts by water and 60% ethanol solvent were analyzed for total polyphenols, antioxidant activity, cyclic voltammetry (CV), and MFCs. The bioactive compounds of H. cordata leaves extracts were assayed via LC/MS analysis. Identification of the marker substances for potential antiviral activity using a molecular docking model was provided.

Significant findings

60% ethanol extract exhibits the highest total polyphenols and antioxidant activity compared with water extracts. Bioenergy extraction in MFCs showed that 60% ethanol extracts could give 1.76-fold more power generation compared to the blank. Flavonoids and their sugar-to-glycan ratios increased after CV scanning and they are expected to be effective ES substances. Quercitrin, from the H. cordata extract that shares an ES-like structure, was found to exhibit strong binding affinities towards ACE2 and RdRp. This indicated the potential of H. cordata leaves as a promising antiviral herb.

Target-specificity of different amyrin subunits in impeding HCV influx mechanism inside the human cells considering the quantum tunnel profiles and molecular strings of the CD81 receptor: a combined in silico and in vivo study

HCV is a hepatotropic RNA virus recognized for its frequent virulence and fatality worldwide. Despite many vaccine development programs underway, researchers are on a quest for natural bioactive compounds due to their multivalent efficiencies against viral infections, considering which the current research aimed to figure out the target-specificity and therapeutic potentiality of α, β, and δ subunits of amyrin, as novel bioactive components against the HCV influx mechanism. Initially, the novelty of amyrin subunits was conducted from 203 pharmacophores, comparing their in-silico pharmacokinetic and pharmacodynamic profiles. Besides, the best active site of CD81 was determined following the quantum tunneling algorithm. The molecular dynamic simulation was conducted (100 ns) following the molecular docking steps to reveal the parameters- RMSD (Å); Cα; RMSF (Å); MolSA (Ų); Rg (nm); PSA (Å); SASA (Ų), and the MM-GBSA dG binding scores. Besides, molecular strings of CD81, along with the co-expressed genes, were classified, as responsible for encoding CD81-mediated protein clusters during HCV infection, resulting in the potentiality of amyrins as targeted prophylactics in HCV infection. Finally, in vivo profiling of the oxidative stress marker, liver-specific enzymes, and antioxidant markers was conducted in the DMN-induced mice model, where β-amyrin scored the most significant values in all aspects.

MERS virus spike protein HTL-epitopes selection and multi-epitope vaccine design using computational biology

MERS-CoV, a zoonotic virus, poses a serious threat to public health globally. Thus, it is imperative to develop an effective vaccination strategy for protection against MERS-CoV. Immunoinformatics and computational biology tools provide a faster and more cost-effective strategy to design potential vaccine candidates. In this work, the spike proteins from different strains of MERS-CoV were selected to predict HTL-epitopes that show affinity for T-helper MHC-class II HTL allelic determinant (HLA-DRB1:0101). The antigenicity and conservation of these epitopes among the selected spike protein variants in different MERS-CoV strains were analyzed. The analysis identified five epitopes with high antigenicity: QSIFYRLNGVGITQQ, DTIKYYSIIPHSIRS, PEPITSLNTKYVAPQ, INGRLTTLNAFVAQQ and GDMYVYSAGHATGTT. Then, a multi-epitope vaccine candidate was designed using linkers and adjuvant molecules. Finally, the vaccine construct was subjected to molecular docking with TLR5 (Toll-like receptor-5). The proposed vaccine construct had strong binding energy of -32.3 kcal/mol when interacting with TLR5.Molecular dynamics simulation analysis showed that the complex of the vaccine construct and TLR5 is stable. Analysis using in silico immune simulation also showed that the prospective multi-epitope vaccine design had the potential to elicit a response within 70 days, with the immune system producing cytokines and immunoglobulins. Finally, codon adaptation and in silico cloning analysis showed that the candidate vaccine could be expressed in the Escherichia coli K12 strain. Here we also designed support vaccine construct MEV-2 by using B-cell and CD8+ CTL epitopes to generate the complete immunogenic effect. This study opens new avenues for the extension of research on MERS vaccine development.Communicated by Ramaswamy H. Sarma.

Structural differences between hypoxanthine phosphoribosyltransferase family members highlight opportunities for antiparasitic drug design in neglected diseases

Approaches to identify novel druggable targets for treating neglected diseases include computational studies that predict possible interactions of drugs and their molecular targets. Hypoxanthine phosphoribosyltransferase (HPRT) plays a central role in the purine salvage pathway. This enzyme is essential for the survival of the protozoan parasite T. cruzi, the causal agent of Chagas disease, and other parasites related to neglected diseases. Here we found dissimilar functional behaviours between TcHPRT and the human homologue, HsHPRT, in the presence of substrate analogues that can lie in differences in their oligomeric assemblies and structural features. To shed light on this issue, we carried out a comparative structural analysis between both enzymes. Our results show that HsHPRT is considerably more resistant to controlled proteolysis than TcHPRT. Moreover, we observed a variation in the length of two key loops depending on the structural arrangement of each protein (groups D1T1 and D1T1'). Such variations might be involved in inter-subunit communication or influencing the oligomeric state. Besides, to understand the molecular basis that govern D1T1 and D1T1' folding groups, we explored the distribution of charges on the interaction surfaces of TcHPRT and HsHPRT, respectively. To know whether the rigidity degree bears effect on the active site, we studied the flexibility of both proteins. The analysis performed here illuminates the underlying reasons and significance behind each protein's preference for one or the other quaternary arrangement that can be exploited for therapeutic approaches.

Heparan sulfate regulates IL-21 bioavailability and signal strength that control germinal center B cell selection and differentiation

In antibody responses, mutated germinal center B (B_GC) cells are positively selected for reentry or differentiation. As the products from GCs, memory B cells and antibody-secreting cells (ASCs) support high-affinity and long-lasting immunity. Positive selection of B_GC cells is controlled by signals received through the B cell receptor (BCR) and follicular helper T (T_FH) cell-derived signals, in particular costimulation through CD40. Here, we demonstrate that the T_FH cell effector cytokine interleukin-21 (IL-21) joins BCR and CD40 in supporting B_GC selection and reveal that strong IL-21 signaling prioritizes ASC differentiation in vivo. B_GC cells, compared with non-B_GC cells, show significantly reduced IL-21 binding and attenuated signaling, which is mediated by low cellular heparan sulfate (HS) sulfation. Mechanistically, N-deacetylase and N-sulfotransferase 1 (Ndst1)-mediated N-sulfation of HS in B cells promotes IL-21 binding and signal strength. Ndst1 is down-regulated in B_GC cells and up-regulated in ASC precursors, suggesting selective desensitization to IL-21 in B_GC cells. Thus, specialized biochemical regulation of IL-21 bioavailability and signal strength sets a balance between the stringency and efficiency of GC selection.

Facilitating CG Simulations with MAD: The MArtini Database Server

The MArtini Database (MAD - https://mad.ibcp.fr) is a web server designed for the sharing of structures and topologies of molecules parametrized with the Martini coarse-grained (CG) force field. MAD can also convert atomistic structures into CG structures and prepare complex systems (including proteins, lipids, etc.) for molecular dynamics (MD) simulations at the CG level. It is dedicated to the generation of input files for Martini 3, the most recent version of this popular CG force field. Specifically, the MAD server currently includes tools to submit or retrieve CG models of a wide range of molecules (lipids, carbohydrates, nanoparticles, etc.), transform atomistic protein structures into CG structures and topologies, with fine control on the process and assemble biomolecules into large systems, and deliver all files necessary to start simulations in the GROMACS MD engine.

In Silico Study of Coumarins: Wedelolactone as a Potential Inhibitor of the Spike Protein of the SARS-CoV-2 Variants

Despite the rigorous global efforts to control SARS-CoV-2 transmission, it continues to pose a serious threat to humans with the frequent emergence of new variants. Thus, robust therapeutics to combat the virus are a desperate need. The SARS-CoV-2 spike (S) protein is an important target protein as it mediates the entry of the virus inside the host cells, which is initiated by the binding of the receptor-binding domain (RBD) to its cognate receptor, angiotensin-converting enzyme 2 (ACE-2). Herein, the inhibition potential of several naturally occurring coumarins was investigated against the spike proteins of SARS-CoV-2 variants using computational approaches. Molecular docking studies revealed 26 coumarins with better binding energies than the reference ligands, molnupiravir and ceftazidime, against the S-RBD of the omicron variant. The top 10 best-docked coumarins were further analyzed to understand their binding interactions against the spike proteins of other variants (wild-type, Alpha, Beta, Gamma, and Delta), and these studies also demonstrated decent binding energies. Physicochemical, QSAR, and pharmacokinetics analyses of the coumarins revealed wedelolactone as the best inhibitor of the spike protein with ideal Lipinski’s drug-likeness and optimal ADMET properties. Furthermore, coarse-grained molecular dynamics (MD) simulation studies of spike protein-wedelolactone complexes validated the stable binding of wedelolactone in the respective binding pockets. As an outcome, wedelolactone could be utilized to develop a potent drug candidate against COVID-19 by blocking the viral entry into the host cell.