UCSF Chimera--a visualization system for exploratory research and analysis

Rescaling protein-protein interactions improves Martini 3 for flexible proteins in solution

Multidomain proteins with flexible linkers and disordered regions play important roles in many cellular processes, but characterizing their conformational ensembles is difficult. We have previously shown that the coarse-grained model, Martini 3, produces too compact ensembles in solution, that may in part be remedied by strengthening protein-water interactions. Here, we show that decreasing the strength of protein-protein interactions leads to improved agreement with experimental data on a wide set of systems. We show that the 'symmetry' between rescaling protein-water and protein-protein interactions breaks down when studying interactions with or within membranes; rescaling protein-protein interactions better preserves the binding specificity of proteins with lipid membranes, whereas rescaling protein-water interactions preserves oligomerization of transmembrane helices. We conclude that decreasing the strength of protein-protein interactions improves the accuracy of Martini 3 for IDPs and multidomain proteins, both in solution and in the presence of a lipid membrane.

Functionalized azirine based scaffolds as endothelin inhibitors for the selective anti-angiogenic activity

Anti-angiogenic therapy has long been used as an adjunct therapy for the resolution of tumor burden. The current findings describe the synthesis of novel marine-based azirine-containing compounds that exhibit anti-angiogenic mediated anti-tumor activity. Azirine-2-carboxylate inhibited HUVEC-mediated tubulogenesis without causing cell death in a dose-dependent manner. Ex-vivo CAM, in-vivo Matrigel implantation, and ear angiogenesis experiments have all shown that azirine-2-carboxylate effectively inhibits angiogenesis. Furthermore, azirine-2-carboxylate inhibits the migration of ECs without disrupting the preformed tubule network. Azirine-2-carboxylate had adequate intramuscular systemic exposure and inhibited tumor growth in a xenograft mouse model. DARTS analysis, competitive binding assay, and gene expression investigations revealed that azirine-2-carboxylate inhibits endothelin-1-mediated angiogenesis. Overall, the discovery of azirine-2-carboxylate demonstrated a potent inhibition of angiogenesis targeting ET1 and a possible application in anti-angiogenic therapy.

Discovery of Triazolopyrimidine Derivatives as Selective P2X3 Receptor Antagonists Binding to an Unprecedented Allosteric Site as Evidenced by Cryo-Electron Microscopy

The P2X3 receptor (P2X3R), an ATP-gated cation channel predominantly expressed in C- and Aδ-primary afferent neurons, has been proposed as a drug target for neurological inflammatory diseases, e.g., neuropathic pain, and chronic cough. Aiming to develop novel, selective P2X3R antagonists, tetrazolopyrimidine-based hit compound 9 was optimized through structure-activity relationship studies by modifying the tetrazole core as well as side chain substituents. The optimized antagonist 26a, featuring a cyclopropane-substituted triazolopyrimidine core, displayed potent P2X3R-antagonistic activity (IC50 = 54.9 nM), 20-fold selectivity versus the heteromeric P2X2/3R, and high selectivity versus other P2XR subtypes. Noncompetitive P2X3R blockade was experimentally confirmed by calcium influx assays. Cryo-electron microscopy revealed that 26a stabilizes the P2X3R in its desensitized state, acting as a molecular barrier to prevent ions from accessing the central pore. In vivo studies in a rat neuropathic pain model (spinal nerve ligation) showed dose-dependent antiallodynic effects of 26a, thus presenting a novel, promising lead structure.

A new framework for novel analogues of pazopanib as potent and selective human carbonic anhydrase inhibitors: Design, repurposing rational, synthesis, crystallographic, in vivo and in vitro biological assessments

Herein, we describe the design and synthesis of novel aryl pyrimidine benzenesulfonamides APBSs 5a-n, 6a-c, 7a-b, and 8 as pazopanib analogues to explore new potent and selective inhibitors for the CA IX. All APBSs were examined in vitro for their promising inhibition activity against a small panel of hCAs (isoforms I, II, IX, and XII). The X-ray crystal structure of CA I in adduct with a representative APBS analogue was solved. APBS-5m, endowed with the best hCA IX inhibitory efficacy and selectivity, was evaluated for antiproliferative activity against a small panel of different cancer cell lines, SK-MEL-173, MDA-MB-231, A549, HCT-116, and HeLa, and it demonstrated one-digit IC50 values range from 2.93 μM (MDA-MB-231) to 5.86 μM (A549). Furthermore, compound APBS-5m was evaluated for its influence on hypoxia-inducible factor (HIF-1α) production, apoptosis induction, and colony formation in MDA-MB-231 cancer cells. The in vivo efficacy of APBS-5m as an antitumor agent was additionally investigated in an animal model of Solid Ehrlich Carcinoma (SEC). In order to offer perceptions into the conveyed hCA IX inhibitory efficacy and selectivity in silico, a molecular docking investigation was also carried out.

Essential oil components interacting with insect odorant-binding proteins: a molecular modelling approach

Essential oils (EOs) are natural products currently used to control arthropods, and their interaction with insect odorant-binding proteins (OBPs) is fundamental for the discovery of new repellents. This in silico study aimed to predict the potential of EO components to interact with odorant proteins. A total of 684 EO components from PubChem were docked against 23 odorant binding proteins from Protein Data Bank using AutoDock Vina. The ligands and proteins were optimized using Gaussian 09 and Sybyl-X 2.0, respectively. The nature of the protein-ligand interactions was characterized using LigandScout 4.0, and visualization of the binding mode in selected complexes was carried out by Pymol. Additionally, complexes with the best binding energy in molecular docking were subjected to 500 ns molecular dynamics simulations using Gromacs. The best binding affinity values were obtained for the 1DQE-ferutidine (-11 kcal/mol) and 2WCH-kaurene (-11.2 kcal/mol) complexes. Both are natural ligands that dock onto those proteins at the same binding site as DEET, a well-known insect repellent. This study identifies kaurene and ferutidine as possible candidates for natural insect repellents, offering a potential alternative to synthetic chemicals like DEET.

Synthesis, in vitro anticancer activity, and pharmacokinetic profiling of the new tetrahydropyrimidines: Part I

Different vanillin-based aldehydes were used to synthesize novel tetrahydropyrimidines (THPMs) via conventional Biginelli reaction. The THPMs were tested against human normal cells (MRC-5) and cancer cell lines (HeLa, K562, and MDA-MB-231). With IC50 values of 10.65, 10.70, and 12.76 µM, compounds 4g, 4h, and 4i exerted the strongest cytotoxic effects against K562 cells. The best activity was achieved for 4g on MDA-MB-231 cells (IC50 = 9.20 ± 0.14 µM). The effects of compounds 4g, 4h, and 4i on the cell-cycle phase distribution of K562 cells were analyzed. Principal component analysis was carried out for the chemometrics analysis to comprehend the relationship between the anticancer activity of the THPMs, pharmacokinetic properties, and partition coefficients, as well as the relationship between the chromatographic behavior and retention parameters. The highest retention rates are found for molecules 4g, 4h, and 4i, which have the longest carbon chains, indicating that the length of the alkyl chain positively affects the molecule's anticancer activity but only if the number of carbon atoms is not higher than seven. Additionally, molecular docking analysis was performed to determine the preferred binding modes of the investigated ligands (4g, 4h, and 4i) with a DNA dodecamer and bovine serum albumin.

Structural basis for linker histone H5-nucleosome binding and chromatin fiber compaction

The hierarchical packaging of chromatin fibers plays a critical role in gene regulation. The 30-nm chromatin fibers, a central-level structure bridging nucleosomal arrays to higher-order organizations, function as the first level of transcriptional dormant chromatin. The dynamics of 30-nm chromatin fiber play a crucial role in biological processes related to DNA. Here, we report a 3.6-angstrom resolution cryogenic electron microscopy structure of H5-bound dodecanucleosome, i.e., the chromatin fiber reconstituted in the presence of linker histone H5, which shows a two-start left-handed double helical structure twisted by tetranucleosomal units. An atomic structural model of the H5-bound chromatin fiber, including an intact chromatosome, is built, which provides structural details of the full-length linker histone H5, including its N-terminal domain and an HMG-motif-like C-terminal domain. The chromatosome structure shows that H5 binds the nucleosome off-dyad through a three-contact mode in the chromatin fiber. More importantly, the H5-chromatin structure provides a fine molecular basis for the intra-tetranucleosomal and inter-tetranucleosomal interactions. In addition, we systematically validated the physiological functions and structural characteristics of the tetranucleosomal unit through a series of genetic and genomic studies in Saccharomyces cerevisiae and in vitro biophysical experiments. Furthermore, our structure reveals that multiple structural asymmetries of histone tails confer a polarity to the chromatin fiber. These findings provide structural and mechanistic insights into how a nucleosomal array folds into a higher-order chromatin fiber with a polarity in vitro and in vivo.

Rational Design, Synthesis, and Computational Investigation of Dihydropyridine [2,3-d] Pyrimidines as Polyphenol Oxidase Inhibitors with Improved Potency

Polyphenol oxidase (PPO) is an industrially important enzyme associated with browning reactions. In the present study, a set of ten new dihydropyridine [2,3-d] pyrimidines (TD-Hid-1-10) were synthesized and was found to be proven characteristically by 1H NMR, 13C NMR, IR, elemental analysis, and assessed as possible PPO inhibitors. PPO was purified from banana using three-phase partitioning, achieving an 18.65-fold purification and 136.47% activity recovery. Enzyme kinetics revealed that the compounds TD-Hid-6 and TD-Hid-7 are to be the most potent inhibitors, exhibiting mixed-type inhibition profile with IC50 values of 1.14 µM, 5.29 µM respectively against purified PPO enzyme. Electronic structure calculations at the B3LYP/PBE0 level of theories using def-2 SVP, def2-TZVP basis sets with various molecular descriptors characterized the electronic behavior of studied derivatives TD-Hid-1-10. Molecular electrostatic potential (MEP) and reduced density gradient analyses of RDG-NCI provided insights into charge distributions and weak intermolecular interactions. Docking study simulations predicted binding poses within crucial amino acid sequence in the 2y9x enzyme's active site, which is typically similar in sequence to the PPO form is not allowed. Ligands were analysed in terms of binding energies, inhibitor concentrations (mM) and various molecular interactions such as H-bonds, H-carbon, π-carbon, π-sigma, π-sigma, π-π T-shaped, π-π stacked, π-alkyl, Van der Waals and Cu interactions. The lowest binding energy (-7.83 kcal/mol) and the highest inhibitory effect (1.83 mM) were shown by the ligand Td-Hid-6, which forms H-bonds with Met280 and Asn260, exhibits π-sigma interactions with His61 and π-alkyl interactions with Val283. Other ligands also showed different interactions with various amino acids; for example, the Td-Hid-1 ligand formed H-bonds with His244 and showed π-sigma interactions with His244 and Val283.

Structural basis of the obligatory exchange mode of human neutral amino acid transporter ASCT2

ASCT2 is an obligate exchanger of neutral amino acids, contributing to cellular amino acid homeostasis. ASCT2 belongs to the same family (SLC1) as Excitatory Amino Acid Transporters (EAATs) that concentrate glutamate in the cytosol. The mechanism that makes ASCT2 an exchanger rather than a concentrator remains enigmatic. Here, we employ cryo-electron microscopy and molecular dynamics simulations to elucidate the structural basis of the exchange mechanism of ASCT2. We establish that ASCT2 binds three Na+ ions per transported substrate and visits a state that likely acts as checkpoint in preventing Na+ ion leakage, both features shared with EAATs. However, in contrast to EAATs, ASCT2 retains one Na+ ion even under Na+-depleted conditions. We demonstrate that ASCT2 cannot undergo the structural transition in TM7 that is essential for the concentrative transport cycle of EAATs. This structural rigidity and the high-affinity Na+ binding site effectively confine ASCT2 to an exchange mode.

Edaravone N-benzyl pyridinium derivatives: BACE-1 inhibition, kinetics and in silico binding pose determination

BACE-1 plays a pivotal role in the production of β-amyloid (Aβ) peptides, implicated in Alzheimer's Disease (AD) pathology. We previously described edaravone N-benzyl pyridinium derivatives (EBPDs) that exhibited multifunctional activity against multiple AD targets. In this study we explored the EBPDs BACE-1 inhibitory activity to potentially enhance the compounds therapeutic profile. The EBPDs exhibited moderate BACE-1 inhibitory activity (IC50 = 44.10 µM - 123.70 µM) and obtained IC50 values between 2.0 and 5.8-fold greater than resveratrol, a known BACE-1 inhibitor (IC50 = 253.20 µM), in this assay. Compound 3 was the most potent inhibitor with an IC50 of 44.10 µM and a Ki of 19.96 µM and a mixed-type mode of inhibition that favored binding in a competitive manner. Molecular docking identified crucial interactions with BACE-1 active site residues, supported by 100 ns MD simulations. The study highlighted the EBPDs therapeutic potential as BACE-1 inhibitors and multifunctional anti-AD therapeutic agents.

Mechanism of PCNA loading by Ctf18-RFC for leading-strand DNA synthesis

The proliferating cell nuclear antigen (PCNA) clamp encircles DNA to hold DNA polymerases (Pols) to DNA for processivity. The Ctf18-RFC PCNA loader, a replication factor C (RFC) variant, is specific to the leading-strand Pol (Polε). We reveal here the underlying mechanism of Ctf18-RFC specificity to Polε using cryo-electron microscopy and biochemical studies. We found that both Ctf18-RFC and Polε contain specific structural features that direct PCNA loading onto DNA. Unlike other clamp loaders, Ctf18-RFC has a disordered ATPase associated with a diverse cellular activities (AAA+) motor that requires Polε to bind and stabilize it for efficient PCNA loading. In addition, Ctf18-RFC can pry prebound Polε off of DNA, then load PCNA onto DNA and transfer the PCNA-DNA back to Polε. These elements in both Ctf18-RFC and Polε provide specificity in loading PCNA onto DNA for Polε.

Capsid structure of bacteriophage ΦKZ provides insights into assembly and stabilization of jumbo phages

Jumbo phages are a group of tailed bacteriophages with large genomes and capsids. As a prototype of jumbo phage, ΦKZ infects Pseudomonas aeruginosa, a multi-drug-resistant (MDR) opportunistic pathogen leading to acute or chronic infection in immunocompromised individuals. It holds potential to be used as an antimicrobial agent and as a model for uncovering basic phage biology. Although previous low-resolution structural studies have indicated that jumbo phages may have more complicated capsid structures than smaller phages such as HK97, the detailed structures and the assembly mechanism of their capsids remain largely unknown. Here, we report a 3.5-Å-resolution cryo-EM structure of the ΦKZ capsid. The structure unveiled ten minor capsid proteins, with some decorating the outer surface of the capsid and the others forming a complex network attached to the capsid's inner surface. This network seems to play roles in driving capsid assembly and capsid stabilization. Similar mechanisms of capsid assembly and stabilization are probably employed by many other jumbo viruses.

Computational studies on the catalytic potential of the double active site for enzyme engineering

Proteins possessing double active sites have the potential to revolutionise enzyme design strategies. This study extensively explored an enzyme that contains both a natural active site (NAS) and an engineered active site (EAS), focusing on understanding its structural and functional properties. Metadynamics simulations were employed to investigate how substrates interacted with their respective active sites. The results revealed that both the NAS and EAS exhibited similar minimum energy states, indicating comparable binding affinities. However, it became apparent that the EAS had a weaker binding site for the substrate due to its smaller pocket and constrained conformation. Interestingly, the EAS also displayed dynamic behaviour, with the substrate observed to move outside the pocket, suggesting the possibility of substrate translocation. To gain further insights, steered molecular dynamics (SMD) simulations were conducted to study the conformational changes of the substrate and its interactions with catalytic residues. Notably, the substrate adopted distinct conformations, including near-attack conformations, in both the EAS and NAS. Nevertheless, the NAS demonstrated superior binding minima for the substrate compared to the EAS, reinforcing the observation that the engineered active site was less favourable for substrate binding due to its limitations. The QM/MM (Quantum mechanics and molecular mechanics) analyses highlight the energy disparity between NAS and EAS. Specifically, EAS exhibited elevated energy levels due to its engineered active site being located on the surface. This positioning exposes the substrate to solvents and water molecules, adding to the energy challenge. Consequently, the engineered enzyme did not provide a significant advantage in substrate binding over the single active site protein. Further, the investigation of internal channels and tunnels within the protein shed light on the pathways facilitating transport between the two active sites. By unravelling the complex dynamics and functional characteristics of this double-active site protein, this study offers valuable insights into novel strategies of enzyme engineering. These findings establish a solid foundation for future research endeavours aimed at harnessing the potential of double-active site proteins in diverse biotechnological applications.

Structural insights into activation mechanisms on NADase of the bacterial DSR2 anti-phage defense system

As a sirtuin (SIR2) family protein, defense-associated sirtuin2 (DSR2) has been demonstrated to participate in bacterial anti-phage resistance via depleting nicotinamide adenine dinucleotide (NAD+) of infected cells, which can be activated by tail tube protein (TTP) and inhibited by DSR anti-defense 1 (DSAD1) of diverse phages. However, the regulating mechanism remains elusive. Here, we determined the cryo-electron microscopy structure of apo DSR2, as well as the respective complex structures with TTP and DSAD1. Structural analyses and biochemical studies reveal that DSR2 forms a tetramer with a SIR2 central core and two distinct conformations. Monomeric TTP preferentially binds to the closed conformation of DSR2, inducing conformational distortions on SIR2 tetramer assembly to activate its NADase activity. DSAD1 combines with the open conformation of DSR2, directly or allosterically inhibiting TTP activation on DSR2 NAD+ hydrolysis. Our findings decipher the detailed molecule mechanisms for DSR2 NADase activity regulation and lay a foundation for in-depth understanding of the DSR2 anti-phage defense system.

Structural basis for VLDLR recognition by eastern equine encephalitis virus

Eastern equine encephalitis virus (EEEV) is the most virulent alphavirus that infects humans, and many survivors develop neurological sequelae, including paralysis and intellectual disability. Alphavirus spike proteins comprise trimers of heterodimers of glycoproteins E2 and E1 that mediate binding to cellular receptors and fusion of virus and host cell membranes during entry. We recently identified very-low density lipoprotein receptor (VLDLR) and apolipoprotein E receptor 2 (ApoER2) as cellular receptors for EEEV and a distantly related alphavirus, Semliki Forest virus (SFV). Here, we use single-particle cryo-electron microscopy (cryo-EM) to determine structures of the EEEV and SFV spike glycoproteins bound to the VLDLR ligand-binding domain and found that EEEV and SFV interact with the same cellular receptor through divergent binding modes. Our studies suggest that the ability of LDLR-related proteins to interact with viral spike proteins through very small footprints with flexible binding modes results in a low evolutionary barrier to the acquisition of LDLR-related proteins as cellular receptors for diverse sets of viruses.

Engineering of a mammalian VMAT2 for cryo-EM analysis results in non-canonical protein folding

Vesicular monoamine transporter 2 (VMAT2) belongs to the major facilitator superfamily (MFS), and mediates cytoplasmic monoamine packaging into presynaptic vesicles. Here, we present two cryo-EM structures of VMAT2, with a frog VMAT2 adopting a canonical MFS fold and an engineered sheep VMAT2 adopting a non-canonical fold. Both VMAT2 proteins mediate uptake of a selective fluorescent VMAT2 substrate into cells. Molecular docking, substrate binding and transport analysis reveal potential substrate binding mechanism in VMAT2. Meanwhile, caution is advised when interpreting engineered membrane protein structures.

Identification and Characterization of CC-AMP1-like and CC-AMP2-like Peptides in Capsicum spp

Antimicrobial peptides (AMPs) are compounds with a variety of bioactive properties. Especially promising are their antibacterial activities, often toward drug-resistant pathogens. Across different AMP sources, AMPs expressed within plants are relatively underexplored with a limited number of plant AMP families identified. Recently, we identified the novel AMPs CC-AMP1 and CC-AMP2 in ghost pepper plants (Capsicum chinense x frutescens), exerting promising antibacterial activity and not classifying into any known plant AMP family. Herein, AMPs related to CC-AMP1 and CC-AMP2 were identified within both Capsicum annuum and Capsicum baccatum. In silico predictions throughout plants were utilized to illustrate that CC-AMP1-like and CC-AMP2-like peptides belong to two broader AMP families, with three-dimensional structural predictions indicating that CC-AMP1-like peptides comprise a novel subfamily of α-hairpinins. The antibacterial activities of several closely related CC-AMP1-like peptides were compared with a truncated version of CC-AMP1 possessing significantly more activity than the full peptide. This truncated peptide was further characterized to possess broad-spectrum antibacterial activity against clinically relevant ESKAPE pathogens. These findings illustrate the value in continued study of plant AMPs toward characterization of novel AMP families, with CC-AMP1-like peptides possessing promising bioactivity.

Pseudophosphorylation of single residues of the J-domain of DNAJA2 regulates the holding/folding balance of the Hsc70 system

The Hsp70 system is essential for maintaining protein homeostasis and comprises a central Hsp70 and two accessory proteins that belong to the J-domain protein (JDP) and nucleotide exchange factor families. Posttranslational modifications offer a means to tune the activity of the system. We explore how phosphorylation of specific residues of the J-domain of DNAJA2, a class A JDP, regulates Hsc70 activity using biochemical and structural approaches. Among these residues, we find that pseudophosphorylation of Y10 and S51 enhances the holding/folding balance of the Hsp70 system, reducing cochaperone collaboration with Hsc70 while maintaining the holding capacity. Truly phosphorylated J domains corroborate phosphomimetic variant effects. Notably, distinct mechanisms underlie functional impacts of these DNAJA2 variants. Pseudophosphorylation of Y10 induces partial disordering of the J domain, whereas the S51E substitution weakens essential DNAJA2-Hsc70 interactions without a large structural reorganization of the protein. S51 phosphorylation might be class-specific, as all cytosolic class A human JDPs harbor a phosphorylatable residue at this position.

The association between genetic variations and morphology-based brain networks changes in Alzheimer's disease

Alzheimer's disease (AD) is a highly heritable disease. The morphological changes of cortical cortex (such as, cortical thickness and surface area) in AD always accompany by the change of the functional connectivity to other brain regions and influence the short- and long-range brain network connections, causing functional deficits of AD. In this study, the first hypothesis is that genetic variations might affect morphology-based brain networks, leading to functional deficits; the second hypothesis is that protein-protein interaction (PPI) between the candidate proteins and known interacting proteins to AD might exist and influence AD. 600 470 variants and structural magnetic resonance imaging scans from 175 AD patients and 214 healthy controls were obtained from the Alzheimer's Disease Neuroimaging Initiative-1 database. A co-sparse reduced-rank regression model was fit to study the relationship between non-synonymous mutations and morphology-based brain networks. After that, PPIs between selected genes and BACE1, an enzyme that was known to be related to AD, are explored by using molecular dynamics (MD) simulation and co-immunoprecipitation (Co-IP) experiments. Eight genes affecting morphology-based brain networks were identified. The results of MD simulation showed that the PPI between TGM4 and BACE1 was the strongest among them and its interaction was verified by Co-IP. Hence, gene variations influence morphology-based brain networks in AD, leading to functional deficits. This finding, validated by MD simulation and Co-IP, suggests that the effect is robust.

Genome-Wide Characterization of Fennel (Anethum foeniculum) MiRNome and Identification of its Potential Targets in Homo sapiens and Arabidopsis thaliana: An Inter and Intra-species Computational Scrutiny

MicroRNAs could be promising biomarkers for various diseases, and small RNA drugs have already been FDA approved for clinical use. This area of research is rapidly expanding and has significant potential for the future. Fennel (Anethum foeniculum) is a highly esteemed spice plant with economic and medicinal benefits, making it an invaluable asset in the pharmaceutical industry. To characterize the fennel miRNAs and their Arabidopsis thaliana and Homo sapience targets with functional enrichment analysis and human disease association. A homology-based computational approach characterized the MiRnome of the Anethum foeniculum genome and assessed its impact on Arabidopsis thaliana and Homo sapience transcriptomes. In addition, functional enrichment analysis was evaluated for both species' targets. Moreover, PPI network analysis, hub gene identification, and MD simulation analysis of the top hub node with fennel miRNA were incorporated. We have identified 100 miRNAs of fennel and their target genes, which include 2536 genes in Homo sapiens and 1314 genes in Arabidopsis thaliana. Functional enrichment analysis reveals 56 Arabidopsis thaliana targets of fennel miRNAs showed involvement in metabolic pathways. Highly enriched human KEGG pathways were associated with several diseases, especially cancer. The protein-protein interaction network of human targets determined the top ten nodes; from them, seven hub nodes, namely MAPK1, PIK3R1, STAT3, EGFR, KRAS, CDC42, and SMAD4, have shown their involvement in the pancreatic cancer pathway. Based on the Blast algorithm, 21 fennel miRNAs are homologs to 16 human miRNAs were predicted; from them, the CSPP1 target was a common target for afo-miR11117a-3p and has-miR-6880-5p homologs miRNAs. Our results are the first to report the 100 fennel miRNAs, and predictions for their endogenous and human target genes provide a basis for further understanding of Anethum foeniculum miRNAs and the biological processes and diseases with which they are associated.

Immunoinformatic-guided novel mRNA vaccine designing to elicit immunogenic responses against the endemic Monkeypox virus

Monkeypox virus (MPXV) is an orthopoxvirus, causing zoonotic infections in humans with smallpox-like symptoms. The WHO reported MPXV cases in May 2022 and the outbreak caused significant morbidity threats to immunocompromised individuals and children. Currently, no clinically validated therapies are available against MPXV infections. The present study is based on immunoinformatics approaches to design mRNA-based novel vaccine models against MPXV. Three proteins were prioritized based on high antigenicity, low allergenicity, and toxicity values to predict T- and B-cell epitopes. Lead T- and B-cell epitopes were used to design vaccine constructs, linked with epitope-specific linkers and adjuvant to enhance immune responses. Additional sequences, including Kozak sequence, MITD sequence, tPA sequence, Goblin 5', 3' UTRs, and a poly(A) tail were added to design stable and highly immunogenic mRNA vaccine construct. High-quality structures were predicted by molecular modeling and 3D-structural validation of the vaccine construct. Population coverage and epitope-conservancy speculated broader protection of designed vaccine model against multiple MPXV infectious strains. MPXV-V4 was eventually prioritized based on its physicochemical and immunological parameters and docking scores. Molecular dynamics and immune simulations analyses predicted significant structural stability and binding affinity of the top-ranked vaccine model with immune receptors to elicit cellular and humoral immunogenic responses against the MPXV. The pursuance of experimental and clinical follow-up of these prioritized constructs may lay the groundwork to develop safe and effective vaccine against MPXV.Communicated by Ramaswamy H. Sarma.

Anticonvulsant activity of Tetrahydrolinalool: behavioral, electrophysiological, and molecular docking approaches

Tetrahydrolinalool (THL) is an acyclic monoterpene alcohol, produced during linalol metabolism and also a constituent of essential oils. As described in the literature, many monoterpenes present anticonvulsant properties, and thus we became interested in evaluating the anticonvulsant activity of Tetrahydrolinalool using in mice model as well as in silico approaches. Our results demonstrated that THL increased latency to seizure onset and also reduced the mortality, in picrotoxin induced seizure tests. The results may be related to GABAergic regulation, which was also suggested in seizure testing induced by 3-mercapto-propionic acid. In the strychnine-induced seizure testing, none of the groups pretreated with THL modulated the parameters indicative of anticonvulsant effect. The electrophysiological results revealed that THL treatment reduces seizures induced by pentylenetetrazole. The in silico molecular docking studies showed that the interaction between THL and a GABAA receptor model formed a stable complex, in comparison to the crystaligraphic structure of diazepam, a structurally related ligand. In conclusion, all the evidences showed that THL presents effective anticonvulsant activity related to the GABAergic pathway, being a candidate for treatment of epileptic syndromes.

Identification of Four Novel Candidate Genes for Non-syndromic Intellectual Disability in Pakistani Families

Intellectual disability, a genetically and clinically varied disorder and is a significant health problem, particularly in less developed countries due to larger family size and high ratio of consanguineous marriages. In the current genetic study, we investigate and find the novel disease causative factors in the four Pakistani families with severe type of non-syndromic intellectual disability. For genetic analysis whole-exome sequencing (WES) and Sanger sequencing was performed. I-TASSER and Cluspro tools were used for Protein modeling and Protein-protein docking. Sanger sequencing confirms the segregation of novel homozygous variants in all the families i.e., c.245 T > C; p.Leu82Pro in SLC50A1 gene in family 1, missense variant c.1037G > A; p.Arg346His in TARS2 gene in family 2, in family 3 and 4, nonsense mutation c.234G > A; p.Trp78Term and missense mutation c.2200G > A; p.Asp734Asn in TBC1D3 and ANAPC2 gene, respectively. In silico functional studies have found the drastic effect of these mutations on protein structure and its interaction properties. Substituted amino acids were highly conserved and present on highly conserved region throughout the species. The discovery of pathogenic variants in SLC50A1, TARS2, TBC1D1 and ANAPC2 shows that the specific pathways connected with these genes may be important in cognitive impairment. The decisive role of pathogenic variants in these genes cannot be determined with certainty due to lack of functional data. However, exome sequencing and segregation analysis of all filtered variants revealed that the currently reported variants were the only variations from the respective families that segregated with the phenotype in the family.

Cumulative asparagine to aspartate deamidation fails to perturb γD-crystallin structure and stability

Deamidation frequently is invoked as an important driver of crystallin aggregation and cataract formation. Here, we characterized the structural and biophysical consequences of cumulative Asn to Asp changes in γD-crystallin. Using NMR spectroscopy, we demonstrate that N- or C-terminal domain-confined or fully Asn to Asp changed γD-crystallin exhibits essentially the same 1H-15N HSQC spectrum as the wild-type protein, implying that the overall structure is retained. Only a very small thermodynamic destabilization for the overall Asn to Asp γD-crystallin variants was noted by chaotropic unfolding, and assessment of the colloidal stability, by measuring diffusion interaction parameters, yielded no substantive differences in association propensities. Furthermore, using molecular dynamics simulations, no significant changes in dynamics for proteins with Asn to Asp or iso-Asp changes were detected. Our combined results demonstrate that substitution of all Asn by Asp residues, reflecting an extreme case of deamidation, did not affect the structure and biophysical properties of γD-crystallin. This suggests that these changes alone cannot be the major determinant in driving cataract formation.

The phosphatase RosC from Streptomyces davaonensis is used for roseoflavin biosynthesis and has evolved to largely prevent dephosphorylation of the important cofactor riboflavin-5'-phosphate

The antibiotic roseoflavin is a riboflavin (vitamin B2) analog. One step of the roseoflavin biosynthetic pathway is catalyzed by the phosphatase RosC, which dephosphorylates 8-demethyl-8-amino-riboflavin-5'-phosphate (AFP) to 8-demethyl-8-amino-riboflavin (AF). RosC also catalyzes the potentially cell-damaging dephosphorylation of the AFP analog riboflavin-5'-phosphate also called "flavin mononucleotide" (FMN), however, with a lower efficiency. We performed X-ray structural analyses and mutagenesis studies on RosC from Streptomyces davaonensis to understand binding of the flavin substrates, the distinction between AFP and FMN and the catalytic mechanism of this enzyme. This work is the first structural analysis of an AFP phosphatase. Each monomer of the RosC dimer consists of an α/β-fold core, which is extended by three specific elongated strand-to-helix sections and a specific N-terminal helix. Altogether these segments envelope the flavin thereby forming a novel flavin-binding site. We propose that distinction between AFP and FMN is provided by substrate-induced rigidification of the four RosC specific supplementary segments mentioned above and by an interaction between the amino group at C8 of AFP and the β -carboxylate of D166. This key amino acid is involved in binding the ring system of AFP and positioning its ribitol phosphate part. Accordingly, site-specific exchanges at D166 disturbed the active site geometry of the enzyme and drastically reduced the catalytic activity. Based on the structure of the catalytic core we constructed a whole series of RosC variants but a disturbing, FMN dephosphorylating "killer enzyme", was not generated.

Molecular mechanism of ligand gating and opening of NMDA receptor

Glutamate transmission and activation of ionotropic glutamate receptors are the fundamental means by which neurons control their excitability and neuroplasticity1. The N-methyl-D-aspartate receptor (NMDAR) is unique among all ligand-gated channels, requiring two ligands-glutamate and glycine-for activation. These receptors function as heterotetrameric ion channels, with the channel opening dependent on the simultaneous binding of glycine and glutamate to the extracellular ligand-binding domains (LBDs) of the GluN1 and GluN2 subunits, respectively2,3. The exact molecular mechanism for channel gating by the two ligands has been unclear, particularly without structures representing the open channel and apo states. Here we show that the channel gate opening requires tension in the linker connecting the LBD and transmembrane domain (TMD) and rotation of the extracellular domain relative to the TMD. Using electron cryomicroscopy, we captured the structure of the GluN1-GluN2B (GluN1-2B) NMDAR in its open state bound to a positive allosteric modulator. This process rotates and bends the pore-forming helices in GluN1 and GluN2B, altering the symmetry of the TMD channel from pseudofourfold to twofold. Structures of GluN1-2B NMDAR in apo and single-liganded states showed that binding of either glycine or glutamate alone leads to distinct GluN1-2B dimer arrangements but insufficient tension in the LBD-TMD linker for channel opening. This mechanistic framework identifies a key determinant for channel gating and a potential pharmacological strategy for modulating NMDAR activity.

BAT2: an Open-Source Tool for Flexible, Automated, and Low Cost Absolute Binding Free Energy Calculations

Absolute binding free energy (ABFE) calculations with all-atom molecular dynamics (MD) have the potential to greatly reduce costs in the first stages of drug discovery. Here, we introduce BAT2, the new version of the Binding Affinity Tool (BAT.py), designed to combine full automation of ABFE calculations with high-performance MD simulations, making it a potential tool for virtual screening. We describe and test several changes and new features that were incorporated into the code, such as relative restraints between the protein and the ligand instead of using fixed dummy atoms, support for the OpenMM simulation engine, a merged approach to the application/release of restraints, support for cobinders and proteins with multiple chains, and many others. We also reduced the simulation times for each ABFE calculation, assessing the effect on the expected robustness and accuracy of the calculations.

An integrated in silico approach for the identification of novel potential drug target and chimeric vaccine against Neisseria meningitides strain 331401 serogroup X by subtractive genomics and reverse vaccinology

Neisseria meningitidis, commonly known as the meningococcus, leads to substantial illness and death among children and young adults globally, revealing as either epidemic or sporadic meningitis and/or septicemia. In this study, we have designed a novel peptide-based chimeric vaccine candidate against the N. meningitidis strain 331,401 serogroup X. Through rigorous analysis of subtractive genomics, two essential cytoplasmic proteins, namely UPI000012E8E0(UDP-3-O-acyl-GlcNAc deacetylase) and UPI0000ECF4A9(UDP-N-acetylglucosamine acyltransferase) emerged as potential drug targets. Additionally, using reverse vaccinology, the outer membrane protein UPI0001F4D537 (Membrane fusion protein MtrC) identified by subcellular localization and recognized for its known indispensable role in bacterial survival was identified as a novel chimeric vaccine target. Following a careful comparison of MHC-I, MHC-II, T-cell, and B-cell epitopes, three epitopes derived from UPI0001F4D537 were linked with three types of linkers-GGGS, EAAAK, and the essential PADRE-for vaccine construction. This resulted in eight distinct vaccine models (V1-V8). Among them V1 model was selected as the final vaccine construct. It exhibits exceptional immunogenicity, safety, and enhanced antigenicity, with 97.7 % of its residues in the Ramachandran plot's most favored region. Subsequently, the vaccine structure was docked with the TLR4/MD2 complex and six different HLA allele receptors using the HADDOCK server. The docking resulted in the lowest HADDOCK score of 39.3 ± 9.0 for TLR/MD2. Immune stimulation showed a strong immune response, including antibodies creation and the activation of B-cells, T Cytotoxic cells, T Helper cells, Natural Killer cells, and interleukins. Furthermore, the vaccine construct was successfully expressed in the Escherichia coli system by reverse transcription, optimization, and ligation in the pET-28a (+) vector for the expression study. The current study proposes V1 construct has the potential to elicit both cellular and humoral responses, crucial for the developing an epitope-based vaccine against N. meningitidis strain 331,401 serogroup X.

Molecular basis for differential recognition of an allosteric inhibitor by receptor tyrosine kinases

Understanding kinase-inhibitor selectivity continues to be a major objective in kinase drug discovery. We probe the molecular basis of selectivity of an allosteric inhibitor (MSC1609119A-1) of the insulin-like growth factor-I receptor kinase (IGF1RK), which has been shown to be ineffective for the homologous insulin receptor kinase (IRK). Specifically, we investigated the structural and energetic basis of the allosteric binding of this inhibitor to each kinase by combining molecular modeling, molecular dynamics (MD) simulations, and thermodynamic calculations. We predict the inhibitor conformation in the binding pocket of IRK and highlight that the charged residues in the histidine-arginine-aspartic acid (HRD) and aspartic acid-phenylalanine-glycine (DFG) motifs and the nonpolar residues in the binding pocket govern inhibitor interactions in the allosteric pocket of each kinase. We suggest that the conformational changes in the IGF1RK residues M1054 and M1079, movement of the ⍺C-helix, and the conformational stabilization of the DFG motif favor the selectivity of the inhibitor toward IGF1RK. Our thermodynamic calculations reveal that the observed selectivity can be rationalized through differences observed in the electrostatic interaction energy of the inhibitor in each inhibitor/kinase complex and the hydrogen bonding interactions of the inhibitor with the residue V1063 in IGF1RK that are not attained with the corresponding residue V1060 in IRK. Overall, our study provides a rationale for the molecular basis of recognition of this allosteric inhibitor by IGF1RK and IRK, which is potentially useful in developing novel inhibitors with improved affinity and selectivity.

A SNF1-related protein kinase regulatory subunit functions as a molecular tuner

Metabolic processes in prokaryotic and eukaryotic organisms are often modulated by kinases which are in turn, dependent on Ca2+ and the cyclic mononucleotides cAMP and cGMP. It has been established that some proteins have both kinase and cyclase activities and that active cyclases can be embedded within the kinase domains. Here, we identified phosphodiesterase (PDE) sites, enzymes that hydrolyse cAMP and cGMP, to AMP and GMP, respectively, in some of these proteins in addition to their kinase/cyclase twin-architecture. As an example, we tested the Arabidopsis thaliana KINγ, a subunit of the SnRK2 kinase, to demonstrate that all three enzymatic centres, adenylate cyclase (AC), guanylate cyclase (GC) and PDE, are catalytically active, capable of generating and hydrolysing cAMP and cGMP. These data imply that the signal output of the KINγ subunit modulates SnRK2, consequently affecting the downstream kinome. Finally, we propose a model where a single protein subunit, KINγ, is capable of regulating cyclic mononucleotide homeostasis, thereby tuning stimulus specific signal output.

In silico identification of D449-0032 compound as a putative SARS-CoV-2 Mpro inhibitor

The SARS-CoV-2 pandemic originated the urgency in developing therapeutic resources for the treatment of COVID-19. Despite the current availability of vaccines and some antivirals, the occurence of severe cases of the disease and the risk of the emergence of new virus variants still motivate research in this field. In this context, this study aimed at the computational prospection of likely inhibitors of the main protease (Mpro) of SARS-CoV-2 since inhibiting this enzyme leads to disruption of the viral replication process. The virtual screening of the antiviral libraries Asinex, ChemDiv, and Enamine targeting SARS-CoV-2 Mpro was performed, indicating the D449-0032 compound as a promising inhibitor. Molecular dynamics simulations showed the stability of the protein-ligand complex and in silico predictions of toxicity and pharmacokinetic parameters indicated the probable drug-like behavior of the compound. In vitro and in vivo studies are essential to confirm the Mpro inhibition by the D449-0032.Communicated by Ramaswamy H. Sarma.

Bioinformatics-based identification of GH12 endoxyloglucanases in citrus-pathogenic Penicillium spp

Millions of tons of citrus peel waste are produced every year as a byproduct of the juice industry. Citrus peel is rich in pectin and xyloglucan, but while the pectin is extracted for use in the food industry, the xyloglucan is currently not valorized. To target hydrolytic degradation of citrus peel xyloglucan into oligosaccharides, we have used bioinformatics to identify three glycoside hydrolase 12 (GH12) endoxyloglucanases (EC 3.2.1.151) from the citrus fruit pathogens Penicillium italicum GL-Gan1 and Penicillium digitatum Pd1 and characterized them on xyloglucan obtained by alkaline extraction from citrus peel. The enzymes displayed pH-temperature optima of pH 4.6-5.3 and 35-37°C. PdGH12 from P. digitatum and PiGH12A from P. italicum share 84% sequence identity and displayed similar kinetics, although kcat was highest for PdGH12. In contrast, PiGH12B from P. italicum, which has the otherwise conserved Trp in subsite -4 replaced with a Tyr, displayed a 3 times higher KM and a 4 times lower kcat/KM than PiGH12A, but was the most thermostable enzyme of the three Penicillium-derived endoxyloglucanases. The benchmark enzyme AnGH12 from Aspergillus nidulans was more thermally stable and had a higher pH-temperature optimum than the enzymes from Penicillum spp. The difference in structure of the xyloglucan oligosaccharides extracted from citrus peel xyloglucan and tamarind xyloglucan by the new endoxyloglucanases was determined by LC-MS. The inclusion of citrus peel xyloglucan demonstrated that the endoxyloglucanases liberated fucosylated xyloglucan oligomers, implying that these enzymes have the potential to upgrade citrus peel residues to produce oligomers useful as intermediates or bioactive compounds.

Structure-based modeling to assess binding and endocrine disrupting potential of polycyclic aromatic hydrocarbons in Daniorerio

Environmental contaminants, such as polycyclic aromatic hydrocarbons (PAHs), have raised concerns regarding their potential endocrine-disrupting effects on aquatic organisms, including fish. In this study, molecular docking and molecular dynamics techniques were employed to evaluate the endocrine-disrupting potential of PAHs in zebrafish, as a model organism. A virtual screening with 72 PAHs revealed a correlation between the number of PAH aromatic rings and their binding affinity to proteins involved in endocrine regulation. Furthermore, PAHs with the highest binding affinities for each protein were identified: cyclopenta[cd]pyrene for AR (-9.7 kcal/mol), benzo(g)chrysene for ERα (-11.5 kcal/mol), dibenzo(a,e)pyrene for SHBG (-8.7 kcal/mol), dibenz(a,h)anthracene for StAR (-11.2 kcal/mol), and 2,3-benzofluorene for TRα (-9.8 kcal/mol). Molecular dynamics simulations confirmed the stability of the protein-ligand complexes formed by the PAHs with the highest binding affinities throughout the simulations. Additionally, the effectiveness of the protocol used in this study was demonstrated by the receiver operating characteristic curve (ROC) analysis, which effectively distinguished decoys from true ligands. Therefore, this research provides valuable insights into the endocrine-disrupting potential of PAHs in fish, highlighting the importance of assessing their impact on aquatic ecosystems.

How does a Saccharomyces cerevisiae extract influence the components of isolated rotavirus particles from stool samples collected in a clinical setting from children?

Human Rotavirus (HRV) is the causative pathogen of severe acute enteric infections that cause mortality among children worldwide. This study focuses on developing a new and effective treatment for rotavirus infection using an extract from Saccharomyces cerevisiae, aiming to make this treatment easily accessible to everyone. 15 antigens and 26 antibodies were detected in serum and stool using ELISA. The titers of HRVq1, HRVq2, HRVC1, and HRVC2 on Vero cells were determined to be 1.2x106, 3.0x106, 4.2x106, and 7.5x105 (Plaque forming unit, PFU/ml) four days after infection, respectively. The HRVq1 isolate induced cytopathic effects, i.e., forming multinucleated, rounded, enlarged, and expanding gigantic cells. RT-PCR identified this isolate, and the accession number 2691714 was assigned to GeneBank. The molecular docking analysis revealed that nonstructural proteins (NSPs) NSP1, NSP2, NSP3, NSP4, NSP5, and NSP6 exhibited significant binding with RNA. NSP2 demonstrated the highest binding affinity and the lowest binding energy (-8.9 kcal/mol). This affinity was maintained via hydrophobic interactions and hydrogen bonds spanning in length from 1.12 Å to 3.11 Å. The ADMET and bioactivity predictions indicated that the yeast extract possessed ideal solubility, was nontoxic, and did not cause cancer. The inhibitory constant values predicted for the S. cerevisiae extract in the presence of HRV vital proteins varied from 5.32 to 7.45 mM, indicating its potential as a viable drug candidate. Saccharomyces cerevisiae extract could be utilized as a dietary supplement to combat HRV as an alternative dietary supplement.

Influence of the interaction between p53 and ZNF568 on mitochondrial oxidative phosphorylation

The tumor suppressor p53 plays important roles in suppressing the development and progression of cancer by responding to various stress signals. In addition, p53 can regulate the metabolic pathways of cancer cells by regulating energy metabolism and oxidative phosphorylation. Here, we present a mechanism for the interaction between p53 and ZNF568. Initially, we used X-ray crystallography to determine the irregular loop structure of the ZNF568 KRAB domain; this loop plays an important role in the interaction between p53 and ZNF568. In addition, Cryo-EM was used to examine how the p53 DBD and ZNF568 KRAB domains bind together. The function of ZNF568 on p53-mediated mitochondrial respiration was confirmed by measuring glucose consumption and lactate production. These findings show that ZNF568 can reduce p53-mediated mitochondrial respiratory activity by binding to p53 and inhibiting the transcription of SCO2. SIGNIFICANCE: ZNF568 can directly bind to the p53 DBD and transcriptionally regulate the SCO2 gene. SCO2 transcriptional regulation by interaction between ZNF568 and p53 may regulate the balance between mitochondrial respiration and glycolysis.

Molecular structure and characterization of the Thermochromatium tepidum light-harvesting 1 photocomplex produced in a foreign host

Purple phototrophic bacteria possess light-harvesting 1 and reaction center (LH1-RC) core complexes that play a key role in converting solar energy to chemical energy. High-resolution structures of LH1-RC and RC complexes have been intensively studied and have yielded critical insight into the architecture and interactions of their proteins, pigments, and cofactors. Nevertheless, a detailed picture of the structure and assembly of LH1-only complexes is lacking due to the intimate association between LH1 and the RC. To study the intrinsic properties and structure of an LH1-only complex, a genetic system was constructed to express the Thermochromatium (Tch.) tepidum LH1 complex heterologously in a modified Rhodospirillum rubrum mutant strain. The heterologously expressed Tch. tepidum LH1 complex was isolated in a pure form free of the RC and exhibited the characteristic absorption properties of Tch. tepidum. Cryo-EM structures of the LH1-only complexes revealed a closed circular ring consisting of either 14 or 15 αβ-subunits, making it the smallest completely closed LH1 complex discovered thus far. Surprisingly, the Tch. tepidum LH1-only complex displayed even higher thermostability than that of the native LH1-RC complex. These results reveal previously unsuspected plasticity of the LH1 complex, provide new insights into the structure and assembly of the LH1-RC complex, and show how molecular genetics can be exploited to study membrane proteins from phototrophic organisms whose genetic manipulation is not yet possible.

Identification of novel chemical scaffolds against kinase domain of cancer causing human epidermal growth factor receptor 2: a systemic chemoinformatic approach

The Human epidermal growth factor receptor 2 (HER2) is expressed in high magnitude in several cancers. Designing new drug molecules that target kinase domain of the HER2 enzyme might provide an appealing platform. Considering this, herein, a multi-phase bioinformatic approach is applied to screen diverse natural and chemical scaffolds to identify compounds that fit best at the kinase domain of HER2. By doing so, three compounds; LAS_51187157, LAC_51217113, LAC_51390233 were pointed with docking score of -11.4 kcal/mol, -11.3 kcal/mol and -11.2 kcal/mol, respectively. In molecular dynamic simulation, the complexes behaved in a stable dynamic with no major local/global structural variations. The intermolecular binding free energies were further estimated that concluded LAC_51390233 complex was the most stable and has less entropy energy. The good docked affinity of LAC_51390233 with HER2 was confirmed by WaterSwap absolute binding free energy. The entropy energy demonstrated that LAC_51390233 has less freedom energy compared to others. Similarly, all three compounds revealed very favorable druglike properties and pharmacokinetics. All the selected three compounds were also non-carcinogenic, non-immunotoxicity, non-mutagenicity, and non-cytotoxic. In a nutshell, the compounds are interesting scaffolds and might be subjected to extensive experimental testing to reveal their real biological potency.Communicated by Ramaswamy H. Sarma.

A novel multi-epitope peptide vaccine candidate targeting hepatitis E virus: An in silico approach

Hepatitis E virus (HEV) is a foodborne virus transmitted through the faecal-oral route that causes viral hepatitis in humans worldwide. Ever since its discovery as a zoonotic agent, HEV was isolated from several species with an expanding range of hosts. HEV possesses several features of other RNA viruses but also has certain HEV-specific traits that make its viral-host interactions inimitable. HEV leads to severe morbidity and mortality in immunocompromised people and pregnant women across the world. The situation in underdeveloped countries is even more alarming. Even after creating a menace across the world, we still lack an effective vaccine against HEV. Till date, there is only one licensed vaccine for HEV available only in China. The development of an anti-HEV vaccine that can reduce HEV-induced morbidity and mortality is required. Live attenuated and killed vaccines against HEV are not accessible due to the lack of a tolerant cell culture system, slow viral replication kinetics and varying growth conditions. Thus, the main focus for anti-HEV vaccine development is now on the molecular approaches. In the current study, we have designed a multi-epitope vaccine against HEV through a reverse vaccinology approach. For the first time, we have used viral ORF3, capsid protein and polyprotein altogether for epitope prediction. These are crucial for viral replication and persistence and are major vaccine targets against HEV. The proposed in silico vaccine construct comprises of highly immunogenic and antigenic T-cell and B-cell epitopes of HEV proteins. The construct is capable of inducing an effective and long-lasting host immune response as evident from the simulation results. In addition, the construct is stable, non-allergic and antigenic for the host. Altogether, our findings suggest that the in silico vaccine construct may be useful as a vaccine candidate for preventing HEV infections.

A Visual Environment for Data Driven Protein Modeling and Validation

In structural biology, validation and verification of new atomic models are crucial and necessary steps which limit the production of reliable molecular models for publications and databases. An atomic model is the result of meticulous modeling and matching and is evaluated using a variety of metrics that provide clues to improve and refine the model so it fits our understanding of molecules and physical constraints. In cryo electron microscopy (cryo-EM) the validation is also part of an iterative modeling process in which there is a need to judge the quality of the model during the creation phase. A shortcoming is that the process and results of the validation are rarely communicated using visual metaphors. This work presents a visual framework for molecular validation. The framework was developed in close collaboration with domain experts in a participatory design process. Its core is a novel visual representation based on 2D heatmaps that shows all available validation metrics in a linear fashion, presenting a global overview of the atomic model and provide domain experts with interactive analysis tools. Additional information stemming from the underlying data, such as a variety of local quality measures, is used to guide the user's attention toward regions of higher relevance. Linked with the heatmap is a three-dimensional molecular visualization providing the spatial context of the structures and chosen metrics. Additional views of statistical properties of the structure are included in the visual framework. We demonstrate the utility of the framework and its visual guidance with examples from cryo-EM.

Circulating rotavirus strains in children with acute gastroenteritis in Iran, 1986 to 2023 and their genetic/antigenic divergence compared to approved vaccines strains (Rotarix, RotaTeq, ROTAVAC, ROTASIIL) before mass vaccination: Clues for vaccination policy makers

In the present study, first, rotaviruses that caused acute gastroenteritis in children under five years of age during the time before the vaccine was introduced in Iran (1986 to 2023) are reviewed. Subsequently, the antigenic epitopes of the VP7 and VP4/VP8 proteins in circulating rotavirus strains in Iran and that of the vaccine strains were compared and their genetic differences in histo-blood group antigens (HBGAs) and the potential impact on rotavirus infection susceptibility and vaccine efficacy were discussed. Overall data indicate that rotavirus was estimated in about 38.1 % of samples tested. The most common genotypes or combinations were G1 and P[8], or G1P[8]. From 2015 to 2023, there was a decline in the prevalence of G1P[8], with intermittent peaks of genotypes G3P[8] and G9P[8]. The analyses suggested that the monovalent Rotarix vaccine or monovalent vaccines containing the G1P[8] component might be proper in areas with a similar rotavirus genotype pattern and genetic background as the Iranian population where the G1P[8] strain is the most predominant and has the ability to bind to HBGA secretors. While the same concept can be applied to RotaTeq and RotasIIL vaccines, their complex vaccine technology, which involves reassortment, makes them less of a priority. The ROTASIIL vaccine, despite not having the VP4 arm (P[5]) as a suitable protection option, has previously shown the ability to neutralize not only G9-lineage I strains but also other G9-lineages at high titers. Thus, vaccination with the ROTASIIL vaccine may be more effective in Iran compared to RotaTeq. However, considering the rotavirus genotypic pattern, ROTAVAC might not be a good choice for Iran. Overall, the findings of this study provide valuable insights into the prevalence of rotavirus strains and the potential effectiveness of different vaccines in the Iranian and similar populations.

Molecular dynamics simulations shed light into the donor substrate specificity of vertebrate poly-alpha-2,8-sialyltransferases ST8Sia IV

BACKGROUND

Sialic acids are essential monosaccharides influencing several biological processes and disease states. The sialyltransferases catalyze the transfer of Sia residues to glycoconjugates playing critical roles in cellular recognition and signaling. Despite their importance, the molecular mechanisms underlying their substrate specificity, especially between different organisms, remain poorly understood. Recently, the human ST8Sia IV, a key enzyme in the synthesis of polysialic acids, was found to accept only CMP-Neu5Ac as a sugar-donor, whereas the whitefish Coregonus maraena enzyme showed a wider donor substrate specificity, accepting CMP-Neu5Ac, CMP-Neu5Gc, and CMP-Kdn. However, what causes these differences in donor substrate specificity is unknown.

METHODS

Computational approaches were used to investigate the structural and biochemical determinants of the donor substrate specificity in ST8Sia IV. Accurate structural models of the human and fish ST8Sia IV catalytic domains and their complexes with three sialic acid donors (CMP-Neu5Ac, CMP-Neu5Gc, and CMP-Kdn) were generated. Subsequently, molecular dynamics simulations were conducted to analyze the stability and interactions within these complexes and identify differences in complex stability and substrate binding sites between the two ST8Sia IV.

RESULTS

Our MD simulations revealed that the human enzyme effectively stabilizes CMP-Neu5Ac, whereas CMP-Neu5Gc and CMP-Kdn are unstable and explore different conformations. In contrast, the fish ST8Sia IV stabilizes all three donor substrates. Based on these data, we identified the key interacting residues for the different Sias parts of the substrate donors.

GENERAL SIGNIFICANCE

This work advances our knowledge of the enzymatic mechanisms governing sialic acid transfer, shedding light on the evolutionary adaptations of sialyltransferases.

Computational biophysical characterization of the effect of gatekeeper mutations on the binding of ponatinib to the FGFR kinase

Fibroblast Growth Factor Receptor (FGFR) is connected to numerous downstream signalling cascades regulating cellular behavior. Any dysregulation leads to a plethora of illnesses, including cancer. Therapeutics are available, but drug resistance driven by gatekeeper mutation impedes the treatment. Ponatinib is an FDA-approved drug against BCR-ABL kinase and has shown effective results against FGFR-mediated carcinogenesis. Herein, we undertake molecular dynamics simulation-based analysis on ponatinib against all the FGFR isoforms having Val to Met gatekeeper mutations. The results suggest that ponatinib is a potent and selective inhibitor for FGFR1, FGFR2, and FGFR4 gatekeeper mutations. The extensive electrostatic and van der Waals interaction network accounts for its high potency. The FGFR3_VM mutation has shown resistance towards ponatinib, which is supported by their lesser binding affinity than wild-type complexes. The disengaged molecular brake and engaged hydrophobic spine were believed to be the driving factors for weak protein-ligand interaction. Taken together, the inhibitory and structural characteristics exhibited by ponatinib may aid in thwarting resistance based on Val-to-Met gatekeeper mutations at an earlier stage of treatment and advance the design and development of other inhibitors targeted at FGFRs harboring gatekeeper mutations.

Comparative CO2 and SiO2 hydratase activity of an enzyme from the siliceous demosponge Suberitesdomuncula

Silicase, an enzyme that catalyzes the hydrolysis of silicon-oxygen bonds, is a crucial player in breaking down silicates into silicic acid, particularly in organisms like aquatic sponges with siliceous skeletons. Despite its significance, our understanding of silicase remains limited. This study comprehensively examines silicase from the demosponge Suberites domuncula, focusing on its kinetics toward CO2 as a substrate, as well as its silicase and esterase activity. It investigates inhibition and activation profiles with a range of inhibitors and activators belonging to various classes. By comparing its esterase activity to human carbonic anhydrase II, we gain insights into its enzymatic properties. Moreover, we investigate silicase's inhibition and activation profiles, providing valuable information for potential applications. We explore the evolutionary relationship of silicase with related enzymes, revealing potential functional roles in biological systems. Additionally, we propose a biochemical mechanism through three-dimensional modeling, shedding light on its catalytic mechanisms and structural features for both silicase activity and CO2 hydration. We highlight nature's utilization of enzymatic expertise in silica metabolism. This study enhances our understanding of silicase and contributes to broader insights into ecosystem functioning and Earth's geochemical cycles, emphasizing the intricate interplay between biology and the environment.

Discovery of novel FGF trap small molecules endowed with anti-myeloma activity

Fibroblast growth factors (FGFs) act as proangiogenic and mitogenic cytokines in several cancers, including multiple myeloma (MM). Indeed, corrupted FGF autocrine and paracrine secretion induces an aberrant activation of the FGF receptor (FGFR) signaling sustaining cancer cell spreading and resistance to pharmacological treatments. Thus, FGF traps may represent a promising anti-cancer strategy to hamper the ligand-dependent activation of the FGF/FGFR system. We previously identified NSC12 as the first orally available small molecule FGF trap able to inhibit the growth and progression of several FGF-dependent tumor models. NSC12 is a pregnenolone derivative carrying a 1,1-bis-trifluoromethyl-1,3-propanediol chain in position 17 of the steroid nucleus. Investigation of structure-activity relationships (SARs) provided more potent and specific NSC12 steroid derivatives and highlighted that the C17-side chain is pivotal for the FGF trap activity. Here, a scaffold hopping approach allowed to obtain two FGF trap compounds (22 and 57) devoid of the steroid nucleus and able to efficiently bind FGF2 and to inhibit FGFR activation in MM cells. Accordingly, these compounds exert a potent anti-tumor activity on MM cell lines both in vitro and in vivo and on MM patient-derived primary cells, strongly affecting the survival of both proteasome-inhibitor sensitive and resistant MM cells. These results propose a new therapeutic option for relapsed/refractory MM patients and set the bases for the development of novel FGF traps prone to chemical diversification to be used in the clinic for the treatment of those tumors in which the FGF/FGFR system plays a pivotal role, including MM.

A contribution to lipid digestion of Odobenidae family: Computational analysis of gastric and pancreatic lipases from walrus (Odobenus rosmarus divergens)

Triacylglycerols (TAGs) are a primary energy source for marine mammals during lipid digestion. Walruses (Odobenus rosmarus divergens) consume prey with a high content of long-chain polyunsaturated fatty acids; however, their digestive physiology and lipid digestion remain poorly studied. The present study aims to model and characterize the gastric (PWGL) and pancreatic (PWPL) lipases of Pacific walruses using an in-silico approach. The confident 3D models of PWGL and PWPL were obtained via homology modeling and protein threading and displayed the structural features of lipases. Molecular docking analysis demonstrated substrate selectivity for long-chain TAG (Trieicosapentaenoin; TC20:5n-3) in PWGL and short-chain TAG (Trioctanoin; TC8:0) in PWPL. Molecular dynamics simulations demonstrate that PWGL bound to tridocosahexaenoin (TC22:6n-3), the protein is considerably stable at all three salinity conditions, but fluctuations are observed in the regions associated with catalytic sites and the lid, indicating the potential hydrolysis of the substrate. This is the first study to report on the digestion of TAGs in walruses, including modeling and lipases characterization and proposing a digestive tract for pinnipeds.

Identification of therapeutic drug target of Shigella Flexneri serotype X through subtractive genomic approach and in-silico screening based on drug repurposing

Shigellosis, induced by Shigella flexneri, constitutes a significant health burden in developing nations, particularly impacting socioeconomically disadvantaged communities. Designated as the second most prevalent cause of diarrheal illness by the World Health Organization (WHO), it precipitates an estimated 212,000 fatalities annually. Within the spectrum of S. flexneri strains, serotype X is notably pervasive and resilient, yet its comprehensive characterization remains deficient. The present investigation endeavors to discern potential pharmacological targets and repurpose existing drug compounds against S. flexneri serotype X. Employing the framework of subtractive genomics, the study interrogates the reference genome of S. flexneri Serotype X (strain 2,002,017; UP000001884) to delineate its proteome into categories of non-homologous, non-paralogous, essential, virulent, and resistant constituents, thereby facilitating the identification of therapeutic targets. Subsequently, a screening of approximately 9000 compounds from the FDA library against the identified drug target aims to delineate efficacious agents for combating S. flexneri serotype X infections. The application of subtractive genomics methodology yields prognostic insights, unveiling non-paralogous proteins (n = 4122), non-homologues (n = 1803), essential (n = 1246), drug-like (n = 389), resistant (n = 167), alongside 42 virulent proteins within the reference proteome. This iterative process culminates in the identification of Serine O-acetyltransferase as a viable drug target. Subsequent virtual screening endeavors to unearth FDA-approved medicinal compounds capable of inhibiting Serine O-acetyltransferase. Noteworthy candidates such as DB12983, DB15085, DB16098, DB16185, and DB16262 emerge, exhibiting potential for mitigating S. flexneri Serotype X. Despite the auspicious findings, diligent scrutiny is imperative to ascertain the efficacy and safety profile of the proposed drug candidates vis-à-vis S. flexneri.

Supercharged coiled-coil protein with N-terminal decahistidine tag boosts siRNA complexation and delivery efficiency of a lipoproteoplex

Short interfering RNA (siRNA) therapeutics have soared in popularity due to their highly selective and potent targeting of faulty genes, providing a non-palliative approach to address diseases. Despite their potential, effective transfection of siRNA into cells requires the assistance of an accompanying vector. Vectors constructed from non-viral materials, while offering safer and non-cytotoxic profiles, often grapple with lackluster loading and delivery efficiencies, necessitating substantial milligram quantities of expensive siRNA to confer the desired downstream effects. We detail the recombinant synthesis of a diverse series of coiled-coil supercharged protein (CSP) biomaterials systematically designed to investigate the impact of two arginine point mutations (Q39R and N61R) and decahistidine tags on liposomal siRNA delivery. The most efficacious variant, N8, exhibits a twofold increase in its affinity to siRNA and achieves a twofold enhancement in transfection activity with minimal cytotoxicity in vitro. Subsequent analysis unveils the destabilizing effect of the Q39R and N61R supercharging mutations and the incorporation of C-terminal decahistidine tags on α-helical secondary structure. Cross-correlational regression analyses reveal that the amount of helical character in these mutants is key in N8's enhanced siRNA complexation and downstream delivery efficiency.

Synthesis and evaluation of di-heterocyclic benzazole compounds as potential antibacterial and anti-biofilm agents against Staphylococcus aureus

Cumulative escalation in antibiotic-resistant pathogens necessitates the quest for novel antimicrobial agents, as current options continue to diminish bacterial resistance. Herein, we report the synthesis of di-heterocyclic benzazole structures (12-19) and their in vitro evaluation for some biological activities. Compounds 16 and 17 demonstrated potent antibacterial activity (MIC = 7.81 μg/mL) against Staphylococcus aureus, along with significant anti-biofilm activity. Noteworthy is the capability of Compound 17 to inhibit biofilm formation by at least 50% at sub-MIC (3.90 μg/mL) concentration. Furthermore, both compounds exhibited the potential to inhibit preformed biofilm by at least 50% at the MIC concentration (7.81 μg/mL). Additionally, Compounds 16 and 17 were examined for cytotoxic effects in HFF-1 cells, using the MTT method, and screened for binding interactions within the active site of S. aureus DNA gyrase using in silico molecular docking technique, employing AutoDock 4.2.6 and Schrödinger Glidse programs. Overall, our findings highlight Compounds 16 and 17 as promising scaffolds warranting further optimization for the development of effective antibacterial and anti-biofilm agents.

Structural insights into the human NuA4/TIP60 acetyltransferase and chromatin remodeling complex

The human NuA4/TIP60 co-activator complex, a fusion of the yeast SWR1 and NuA4 complexes, both incorporates the histone variant H2A.Z into nucleosomes and acetylates histones H4/H2A/H2A.Z to regulate gene expression and maintain genome stability. Our cryo-electron microscopy studies show that, within the NuA4/TIP60 complex, the EP400 subunit serves as a scaffold holding the different functional modules in specific positions, creating a unique arrangement of the ARP module. EP400 interacts with the TRRAP subunit using a footprint that overlaps with that of the SAGA acetyltransferase complex, preventing the formation of a hybrid complex. Loss of the TRRAP subunit leads to mislocalization of NuA4/TIP60, resulting in the redistribution of H2A.Z and its acetylation across the genome, emphasizing the dual functionality of NuA4/TIP60 as a single macromolecular assembly.

Neuronal and non-neuronal functions of the synaptic cell adhesion molecule neurexin in Nematostella vectensis

The evolutionary transition from diffusion-mediated cell-cell communication to faster, targeted synaptic signaling in animal nervous systems is still unclear. Genome sequencing analyses have revealed a widespread distribution of synapse-related genes among early-diverging metazoans, but how synaptic machinery evolved remains largely unknown. Here, we examine the function of neurexins (Nrxns), a family of presynaptic cell adhesion molecules with critical roles in bilaterian chemical synapses, using the cnidarian model, Nematostella vectensis. Delta-Nrxns are expressed mainly in neuronal cell clusters that exhibit both peptidergic and classical neurotransmitter signaling. Knockdown of δ-Nrxn reduces spontaneous peristalsis of N. vectensis polyps. Interestingly, gene knockdown and pharmacological studies suggest that δ-Nrxn is involved in glutamate- and glycine-mediated signaling rather than peptidergic signaling. Knockdown of the epithelial α-Nrxn reveals a major role in cell adhesion between ectodermal and endodermal epithelia. Overall, this study provides molecular, functional, and cellular insights into the pre-neural function of Nrxns, as well as key information for understanding how and why they were recruited to the synaptic machinery.