VADAR: a web server for quantitative evaluation of protein structure quality

Computational strategies in Klebsiella pneumoniae vaccine design: navigating the landscape of in silico insights

The emergence of multidrug-resistant Klebsiella pneumoniae poses a grave threat to global public health, necessitating urgent strategies for vaccine development. In this context, computational tools have emerged as indispensable assets, offering unprecedented insights into klebsiellal biology and facilitating the design of effective vaccines. Here, a review of the application of computational methods in the development of K. pneumoniae vaccines is presented, elucidating the transformative impact of in silico approaches. Through a systematic exploration of bioinformatics, structural biology, and immunoinformatics techniques, the complex landscape of K. pneumoniae pathogenesis and antigenicity was unravelled. Key insights into virulence factors, antigen discovery, and immune response mechanisms are discussed, highlighting the pivotal role of computational tools in accelerating vaccine development efforts. Advancements in epitope prediction, antigen selection, and vaccine design optimisation are examined, highlighting the potential of in silico approaches to update vaccine development pipelines. Furthermore, challenges and future directions in leveraging computational tools to combat K. pneumoniae are discussed, emphasizing the importance of multidisciplinary collaboration and data integration. This review provides a comprehensive overview of the current state of computational contributions to K. pneumoniae vaccine development, offering insights into innovative strategies for addressing this urgent global health challenge.

Catalytic function of the laccase enzyme in response to chlorpyrifos and 2,4-dichlorophenoxyacetic acid: behavior in controlled and simulated environments

Enzymes secreted by white-rot fungi, such as laccase, offer a promising solution for treating xenobiotic compounds dangerous to the environment and human health. This study aimed to perform a comprehensive analysis of the tolerance of Pleurotus pulmonarius LBM 105 and its laccase activity toward the pesticides 2,4-D and chlorpyrifos both in vitro and in silico. The fungal strain was able to grow in different concentrations of the pesticides, showing evident morphological alterations. Laccase activity and a 53 kDa electromorph were present in all treatments, showing significant stability with peak activity achieved at a pH of 5.6 and within a temperature range of 50-60 °C. Three laccase genes were mapped, annotated, and characterized from the genome. PplacI obtained better structural validation and affinity energy of - 5.05 and - 7.65 kcal mol-1 with 2,4-D and chlorpyrifos, respectively. The Molecular Mechanics/Poisson-Boltzmann Surface Area analysis at 250 ns confirmed the docking results, revealing the existence of stronger hydrophobic interactions between laccase and chlorpyrifos and highlighting the importance of the Phe341 residue in stabilizing both complexes. Understanding the impact of pesticides on laccase's catalytic function is key to formulating and applying future biotechnological strategies with this enzyme.

Characterization, expressional and evolutionary analysis of five fish-specific CCRs (CCR4La, CCR4Lc, CCR12a1, CCR12a2, and CCR12b) in largemouth bass (Micropterus salmoides)

CC chemokine receptors (CCRs), the numbers of the G protein-coupled receptor (GPCR) superfamily, had crucial roles in treating infection, inflammation, and tissue damage by binding to their ligands. In this study, five fish-specific CCRs, namely CCR4La, CCR4Lc, CCR12a1, CCR12a2, and CCR12b, were identified in largemouth bass (Micropterus salmoides). The correction of nomenclatures of these CCRs were confirmed by phylogenetic analysis, structural analysis and genomic synteny analysis. Following 1 × 106 CFU/mL and 1 × 107 CFU/mL Edwardsiella piscicida infection, these five CCRs were significantly induced in spleen of largemouth bass, indicating their important roles in the immune response against bacterial infection. Selection pressure analysis revealed that CCR4La, CCR4Lc, CCR12a1, and CCR12a2 underwent negative selection pressure, whereas CCR12b experienced positive selection pressure. Robust selection site detection methods identified that positive selected sites of CCR4La, CCR4Lc, CCR12a1, and CCR12a2 mainly distributed in their extracellular regions, which involved in ligand binding and pathogen interaction. Similarly, the positive selected sites of CCR12b were also located in its extracellular regions. The accuracy of the pressure selected sites were also validated by molecular docking analysis. The potential ligands for these five CCRs were identified by molecular docking analysis, with finding that CCL3 and CCL5 might be the ligands of largemouth bass CCR4La/Lc, and CCL5, CCL8, CCL7, CCL13 and CCL26 might be that of largemouth bass CCR12a1/a2/b. Our results provided basis for elucidating the functions of chemokine-receptor complex in largemouth bass.

Chemical inhibition of stomatal differentiation by perturbation of the master-regulatory bHLH heterodimer via an ACT-Like domain

Selective perturbation of protein interactions with chemical compounds enables dissection and control of developmental processes. Differentiation of stomata, cellular valves vital for plant growth and survival, is specified by the basic-helix-loop-helix (bHLH) heterodimers. Harnessing a new amination reaction, we here report a synthesis, derivatization, target identification, and mode of action of an atypical doubly-sulfonylated imidazolone, Stomidazolone, which triggers stomatal stem cell arrest. Our forward chemical genetics followed by biophysical analyses elucidates that Stomidazolone directly binds to the C-terminal ACT-Like (ACTL) domain of MUTE, a master regulator of stomatal differentiation, and perturbs its heterodimerization with a partner bHLH, SCREAM in vitro and in plant cells. On the other hand, Stomidazolone analogs that are biologically inactive do not bind to MUTE or disrupt the SCREAM-MUTE heterodimers. Guided by structural docking modeling, we rationally design MUTE with reduced Stomidazolone binding. These engineered MUTE proteins are fully functional and confer Stomidazolone resistance in vivo. Our study identifies doubly-sulfonylated imidazolone as a direct inhibitor of the stomatal master regulator, further expanding the chemical space for perturbing bHLH-ACTL proteins to manipulate plant development.

Investigating the Inhibitory Potential of Flavonoids against Aldose Reductase: Insights from Molecular Docking, Dynamics Simulations, and gmx_MMPBSA Analysis

Diabetes mellitus (DM) is a complex metabolic disorder characterized by chronic hyperglycemia, with aldose reductase playing a critical role in the pathophysiology of diabetic complications. This study aimed to investigate the efficacy of flavonoid compounds as potential aldose reductase inhibitors using a combination of molecular docking and molecular dynamics (MD) simulations. The three-dimensional structures of representative flavonoid compounds were obtained from PubChem, minimized, and docked against aldose reductase using Discovery Studio's CDocker module. The top 10 compounds Daidzein, Quercetin, Kaempferol, Butin, Genistein, Sterubin, Baicalein, Pulchellidin, Wogonin, and Biochanin_A were selected based on their lowest docking energy values for further analysis. Subsequent MD simulations over 100 ns revealed that Daidzein and Quercetin maintained the highest stability, forming multiple conventional hydrogen bonds and strong hydrophobic interactions, consistent with their favorable interaction energies and stable RMSD values. Comparative analysis of hydrogen bond interactions and RMSD profiles underscored the ligand stability. MMPBSA analysis further confirmed the significant binding affinities of Daidzein and Quercetin, highlighting their potential as aldose reductase inhibitors. This study highlights the potential of flavonoids as aldose reductase inhibitors, offering insights into their binding interactions and stability, which could contribute to developing novel therapeutics for DM complications.

Rapid Coaggregation of Proteins Without Sequence Similarity: Possible Role of Conformational Complementarity

Despite extensive research on the sequence-determined self-assembly of both pathogenic and nonpathogenic proteins, the question of how the sequence identity would influence the coassembly or cross-seeding of diverse proteins without distinct sequence similarity remains largely unanswered. Here, we demonstrate that the rapid coaggregation of proteins with negligible sequence similarity is fundamentally governed by preferred heteromeric interactions between their partially unfolded states via the gain of additional charge complementarity and hydrophobic interactions. The partial loss of intramolecular interactions and concurrent gain of non-native intrinsically disordered regions with sticky groups become crucial for both aggressive heteromeric primary nucleation and secondary nucleation events. The results signify the direct relevance of sequence-independent conformational cross-talk between diverse proteins to the foundational events required for the growth of biological multiprotein amyloid deposits.

Adaptive evolution and functional significance of the PPARGC1A gene across diverse animal species

Codon-based analyses of the PPARGC1A gene across 38 vertebrate species were deployed to elucidate patterns of evolutionary change. Employing maximum likelihood assessments through MEGA, we scrutinized 447 codon positions addressing the entire coding region, excluding positions mired by gaps or missing data. Distinct codons manifested variance in selection pressures, particularly codons 4, 11, 66, and 123, which exhibited positive dN-dS values suggestive of positive selection. Codon 137 displayed the most pronounced dN-dS value, signifying intensified selective advantage. Meanwhile, codons 30 and 90 portrayed near-neutral scores, indicative of purifying selection. Complementary computational methods (IFEL, REL, FUBAR, and SLAC) confirmed positive selection at specific codon sites, with varying degrees of corroboration. The integration of mixed-effect modeling (MEME) identified episodic diversifying selection, pinpointing codons that underwent selection episodes in certain lineages. Refined codon model selection lent insight into substitution rates, revealing nuanced degrees of evolutionary conservation among different codons. Supporting these genetic insights, the phylogenetic analysis highlighted relationships among the PPARGC1A sequences and domain analysis confirmed conserved features across species, while protein-protein interaction networks suggested a complex web of functional interdependencies. These findings advance our understanding of the PPARGC1A gene's evolutionary trajectory and underscore the gene's potential adaptive significance within diverse vertebrate lineages.

Identification of Potential Tryptase Inhibitors from FDA-Approved Drugs Using Machine Learning, Molecular Docking, and Experimental Validation

This study explores the innovative use of machine learning (ML) to identify novel tryptase inhibitors from a library of FDA-approved drugs, with subsequent confirmation via molecular docking and experimental validation. Tryptase, a significant mediator in inflammatory and allergic responses, presents a therapeutic target for various inflammatory diseases. However, the development of effective tryptase inhibitors has been challenging due to the enzyme's complex activation and regulation mechanisms. Utilizing a machine learning model, we screened an extensive FDA-approved drug library to identify potential tryptase inhibitors. The predicted compounds were then subjected to molecular docking to assess their binding affinity and conformation within the tryptase active site. Experimental validation was performed using RBL-2H3 cells, a rat basophilic leukemia cell line, where the efficacy of these compounds was evaluated based on their ability to inhibit tryptase activity and suppress β-hexosaminidase activity and histamine release. Our results demonstrated that several FDA-approved drugs, including landiolol, laninamivir, and cidofovir, significantly inhibited tryptase activity. Their efficacy was comparable to that of the FDA-approved mast cell stabilizer nedocromil and the investigational agent APC-366. These findings not only underscore the potential of ML in accelerating drug repurposing but also highlight the feasibility of this approach in identifying effective tryptase inhibitors. This research contributes to the field of drug discovery, offering a novel pathway to expedite the development of therapeutics for tryptase-related pathologies.

Use of plasma induced modification of biomolecules (PLIMB) to evaluate hemin dissociation from fish and bovine methemoglobins

Hemin dissociation occurs much faster from fish methemoglobin (metHb) compared to mammalian metHb yet the mechanism remains poorly understood. This may involve enhanced solvent access to His(E7) of fish metHbs by a protonation mechanism. Plasma induced modification of biomolecules (PLIMB) produces free radicals that covalently modify solvent accessible residues of proteins, and so can provide insight regarding accessibility of hydronium ions to protonate His(E7). PLIMB-induced modifications to heme crevice sites of trout IV and bovine metHb were determined using tandem mass spectrometry after generating peptides with Trypsin/Lys-C. αHis(CE3) was more modified in trout attributable to the more dynamic nature of bovine αHis(CE3) from available crystal structures. Although His(E7) was not found to be more modified in trout, aspects of trout peptides containing His(E7) hampered modification determinations. An existing computational structure-based approach was also used to estimate protonation tendencies, suggesting His(E7) of metHbs with low hemin affinity are more protonatable.

Molecular simulation reveals that pathogenic mutations in BTB/ANK domains of Arabidopsis thaliana NPR1 circumscribe the EDS1-mediated immune regulation

The NPR1 (nonexpressor of pathogenesis-related genes 1) is a key regulator of the salicylic-acid-mediated immune response caused by pathogens in Arabidopsis thaliana. Mutations C150Y and H334Y in the BTB/ANK domains of NPR1 inhibit the defense response, and transcriptional co-activity with enhanced disease susceptibility 1 (EDS1) has been revealed experimentally. This study examined the conformational changes and reduced NPR1-EDS1 interaction upon mutation using a molecular dynamics simulation. Initially, BTBC150YNPR1 and ANKH334YNPR1 were categorized as pathological mutations rather than others based on sequence conservation. A distant ortholog was used to map the common residues shared among the wild-type because the mutations were highly conserved. Overall, 179 of 373 residues were determining the secondary structures and fold versatility of conformations. In addition, the mutational hotspots Cys150, Asp152, Glu153, Cys155, His157, Cys160, His334, Arg339 and Lys370 were crucial for oligomer-to-monomer exchange. Subsequently, the atomistic simulations with free energy (MM/PB(GB)SA) calculations predicted structural displacements engaging in the N-termini α5133-178α7 linker connecting the central ANK regions (α13260-290α14 and α18320-390α22), where prominent long helices (α516) and short helices (α310) replaced with β-turns and loops disrupting hydrogen bonds and salt bridges in both mutants implicating functional regulation and activation. Furthermore, the mutation repositions the intact stability of multiple regions (L13C149-N356α20BTB/ANK-α17W301-E357α21N-ter/coiled-coil) compromising a dynamic interaction of NPR1-EDS1. By unveiling the transitions between the distinct functions of mutational perception, this study paves the way for future investigation to orchestrate additive host-adapted transcriptional reprogramming that controls defense-related regulatory mechanisms of NPR1s in plants.

Mispacking of the F87 sidechain drives aggregation-promoting conformational fluctuations in the subunit interfaces of the transthyretin tetramer

Aberrant formation and deposition of human transthyretin (TTR) aggregates causes transthyretin amyloidosis. To initialize aggregation, transthyretin tetramers must first dissociate into monomers that partially unfold to promote entry into the aggregation pathway. The native TTR tetramer (T) is stabilized by docking of the F87 sidechain into an interfacial cavity enclosed by several hydrophobic residues including A120. We have previously shown that an alternative tetramer (T*) with mispacked F87 sidechains is more prone to dissociation and aggregation than the native T state. However, the molecular basis for the reduced stability in T* remains unclear. Here we report characterization of the A120L mutant, where steric hindrance is introduced into the F87 binding site. The x-ray structure of A120L shows that the F87 sidechain is displaced from its docking site across the subunit interface. In A120S, a naturally occurring pathogenic mutant that is less aggregation-prone than A120L, the F87 sidechain is correctly docked, as in the native TTR tetramer. Nevertheless, 19F-NMR aggregation assays show an elevated population of a monomeric aggregation intermediate in A120S relative to a control containing the native A120, due to accelerated tetramer dissociation and slowed monomer tetramerization. The mispacking of the F87 sidechain is associated with enhanced exchange dynamics for interfacial residues. At 298 K, the T* populations of various naturally occurring mutants fall between 4% and 7% (ΔG ~ 1.5-1.9 kcal/mol), consistent with the free energy change expected for undocking and solvent exposure of one of the four F87 sidechains in the tetramer (ΔG ~ 1.6 kcal/mol). Our data provide a molecular-level picture of the likely universal F87 sidechain mispacking in tetrameric TTR that promotes interfacial conformational dynamics and increases aggregation propensity.

Leech, potato, and tomato carboxypeptidase inhibitors against Anopheles stephensi carboxypeptidase B1 and B2

Carboxypeptidase B (CPB) in Anopheles spp. breaks down blood and releases free amino acids, which promote Plasmodium sexual development in the mosquito midgut. Our goal was to computationally assess the inhibitory effectiveness of carboxypeptidase inhibitors obtained from tomato, potato (CPiSt), and leech against the Anopheles stephensi CPBAs1 and CPBAs2 enzymes. The tertiary structures of CPB inhibitors were predicted and their interaction mode with CPBAs1 and CPBAs2 were examined using molecular docking. Next, this data was compared with four licensed medications that are known to reduce the Anopheles' CPB activity. Molecular dynamics simulations were used to evaluate the stability of complexes containing CPiSt and its mutant form. Both CPiSt and its mutant form showed promise as possible candidates for further evaluations in the paratransgenesis technique for malaria control, based on the similar bindings of CPiSt and CPiSt-Mut to the active sites of CPBAs1 and CPBAs2, as well as their binding affinity in comparison to the drugs.

Stable monomers in the ancestral sequence reconstruction of the last opisthokont common ancestor of dimeric triosephosphate isomerase

Function and structure are strongly coupled in obligated oligomers such as Triosephosphate isomerase (TIM). In animals and fungi, TIM monomers are inactive and unstable. Previously, we used ancestral sequence reconstruction to study TIM evolution and found that before these lineages diverged, the last opisthokonta common ancestor of TIM (LOCATIM) was an obligated oligomer that resembles those of extant TIMs. Notably, calorimetric evidence indicated that ancestral TIM monomers are more structured than extant ones. To further increase confidence about the function, structure, and stability of the LOCATIM, in this work, we applied two different inference methodologies and the worst plausible case scenario for both of them, to infer four sequences of this ancestor and test the robustness of their physicochemical properties. The extensive biophysical characterization of the four reconstructed sequences of LOCATIM showed very similar hydrodynamic and spectroscopic properties, as well as ligand-binding energetics and catalytic parameters. Their 3D structures were also conserved. Although differences were observed in melting temperature, all LOCATIMs showed reversible urea-induced unfolding transitions, and for those that reached equilibrium, high conformational stability was estimated (ΔGTot = 40.6-46.2 kcal/mol). The stability of the inactive monomeric intermediates was also high (ΔGunf = 12.6-18.4 kcal/mol), resembling some protozoan TIMs rather than the unstable monomer observed in extant opisthokonts. A comparative analysis of the 3D structure of ancestral and extant TIMs shows a correlation between the higher stability of the ancestral monomers with the presence of several hydrogen bonds located in the "bottom" part of the barrel.

Structural characterization of DNA-binding domain of essential mammalian protein TTF 1

Transcription Termination Factor 1 (TTF1) is a multifunctional mammalian protein with vital roles in various cellular processes, including Pol I-mediated transcription initiation and termination, pre-rRNA processing, chromatin remodelling, DNA damage repair, and polar replication fork arrest. It comprises two distinct functional regions; the N-terminal regulatory region (1-445 aa), and the C-terminal catalytic region (445-859 aa). The Myb domain located at the C-terminal region is a conserved DNA binding domain spanning from 550 to 732 aa (183 residues). Despite its critical role in various cellular processes, the physical structure of TTF1 remains unsolved. Attempts to purify the functional TTF1 protein have been unsuccessful till date. Therefore, we focused on characterizing the Myb domain of this essential protein. We started with predicting a 3-D model of the Myb domain using homology modelling, and ab-initio method. We then determined its stability through MD simulation in an explicit solvent. The model predicted is highly stable, which stabilizes at 200ns. To experimentally validate the computational model, we cloned and expressed the codon optimized Myb domain into a bacterial expression vector and purified the protein to homogeneity. Further, characterization of the protein shows that, Myb domain is predominantly helical (65%) and is alone sufficient to bind the Sal Box DNA. This is the first-ever study to report a complete in silico model of the Myb domain, which is physically characterized. The above study will pave the way towards solving the atomic structure of this essential mammalian protein.

Integrative workflows for the characterization of hydrophobin and cerato-platanin in the marine fungus Paradendryphiella salina

Hydrophobins (HFBs) and cerato-platanins (CPs) are surface-active extracellular proteins produced by filamentous fungi. This study identified two HFB genes (pshyd1 and pshyd2) and one CP gene (pscp) in the marine fungus Paradendryphiella salina. The proteins PsCP, PsHYD2, and PsHYD1 had molecular weights of 12.70, 6.62, and 5.98 kDa, respectively, with isoelectric points below 7. PsHYD1 and PsHYD2 showed hydrophobicity (GRAVY score 0.462), while PsCP was hydrophilic (GRAVY score - 0.202). Stability indices indicated in-solution stability. Mass spectrometry identified 2,922 proteins, including CP but not HFB proteins. qPCR revealed differential gene expression influenced by developmental stage and substrate, with pshyd1 consistently expressed. These findings suggest P. salina's adaptation to marine ecosystems with fewer hydrophobin genes than other fungi but capable of producing surface-active proteins from seaweed carbohydrates. These proteins have potential applications in medical biocoatings, food industry foam stabilizers, and environmental bioremediation.

Anti-filarial efficacy of Centratherum anthelminticum: unravelling the underlying mechanisms through biochemical, HRAMS proteomics and MD simulation approaches

Traditionally, Centratherum anthelminticum (CA) has been reported to be a potent anti-filarial, however no reports are available detailing its mechanism of action against filarial parasites. In this study, we have evaluated the anti-filarial activity of CA against lymphatic filarial parasites Setaria cervi using ex vivo biochemical, proteomics and in silico approaches. The motility and viability of the parasites decreased significantly after treatment with CA concentrations of ≥125 μg mL-1. An increase in lipid peroxidation (51.92%), protein carbonylation (48.99%), NADPH oxidase (88.88%) activity and decrease in the glutathione (GSH) (-39.23%), glutathione reductase (GR) (-60.17%), and glutathione S-transferase (GST) (-50.48%) activity was also observed after CA treatment. The proteomics analysis was performed by two-dimensional gel electrophoresis and high-resolution accurate mass spectrometry (HRAMS). In total, 185 proteins were differentially expressed (DEPs) following CA treatment. The major DEPs were mostly involved in tRNA processing, biosynthetic processes, metabolic activities, protein transport, the tricarboxylic acid cycle, protein translation, and stress response. The UPLC-ESI-MS/MS analysis of CA extract revealed the presence of 40 bioactive compounds. Further the docking analysis showed 10 CA bioactive compounds to have high binding affinity towards antioxidant proteins of filarial parasites. Additionally, MD simulation studies showed stable interactions (RMSF ≤ 10 Å) of 3-O-methylquercitin, quinic acid, gentisic acid, and vanillin with filarial antioxidant enzymes/proteins. To our knowledge, this is the first report detailing the molecular mechanism of anti-filarial activity of CA, which can be further evaluated for the development of new anti-filarial formulations.

Comparative in silico analysis of CNS-active molecules targeting the blood-brain barrier choline transporter for Alzheimer's disease therapy

This study investigated the blood‒brain barrier (BBB) permeability of the central nervous system (CNS)-active compounds donepezil (DON), methionine (MET), and memantine (MEM) by employing a comprehensive in silico approach. These compounds are of particular interest for Alzheimer's disease (AD) therapy. Rigid-flexible molecular docking simulations indicated favorable binding affinities of all the compounds with BBB-ChT, with DON exhibiting the highest binding affinity (ΔGbind = -10.26 kcal/mol), predominantly mediated by significant hydrophobic interactions. In silico kinetic profiling suggested the stability of the DON/BBB-ChT complex, with ligand release prompted by conformational changes. 3D molecular alignment corroborated a minor conformational shift for DON in its minimal binding energy pose. Predictions indicated that active transport mechanisms notably enhance the brain distribution of donepezil compared to that of MET and MEM. Additionally, DON and MEM exhibited low mutagenic probabilities, while MET was identified as highly mutagenic. Overall, these findings highlight the potential of donepezil for superior BBB penetration, primarily through active transport mechanisms, underscoring the need for further validation through in vitro and in vivo studies for effective AD treatment.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-024-00245-w.

Glycine-replaced epinecidin-1 variant bestows better stability and stronger antimicrobial activity against a range of nosocomial pathogenic bacteria

Epinecidin-1 (epi-1), an antimicrobial peptide first identified in marine grouper fish, has multifunctional bioactivities. The present study aims to improve its therapeutic potential via structural modifications that could enhance its antimicrobial activity and stability. To achieve it, we replaced glycine and the first histidine in the parent epi-1 with lysine, which resulted in a peptide with a repeating KXXK motif and improved physiochemical properties related to antimicrobial activity. This modified peptide, referred to as glycine-to-lysine replaced-epi-1, also gained stability and a twofold increase in helical propensity. To produce the active peptide, overlap extension PCR was employed to generate the gene of GK-epi-1 via site-directed mutagenesis, which was then cloned into the pET-32a vector and expressed as a recombinant fusion protein in Escherichia coli C43 (DE3) strain. The recombinant protein was purified and digested with enterokinase to release the active peptide fragment, which was then evaluated for antimicrobial activity and stability. The lysine substitution led to an enhancement in broad-spectrum antimicrobial activity against a wide range of nosocomial pathogenic bacteria.

Three-Step Enzymatic Remodeling of Chitin into Bioactive Chitooligomers

Here we describe a complex enzymatic approach to the efficient transformation of abundant waste chitin, a byproduct of the food industry, into valuable chitooligomers with a degree of polymerization (DP) ranging from 6 to 11. This method involves a three-step process: initial hydrolysis of chitin using engineered variants of a novel fungal chitinase from Talaromyces flavus to generate low-DP chitooligomers, followed by an extension to the desired DP using the high-yielding Y445N variant of β-N-acetylhexosaminidase from Aspergillus oryzae, achieving yields of up to 57%. Subsequently, enzymatic deacetylation of chitooligomers with DP 6 and 7 was accomplished using peptidoglycan deacetylase from Bacillus subtilis BsPdaC. The innovative enzymatic procedure demonstrates a sustainable and feasible route for converting waste chitin into unavailable bioactive chitooligomers potentially applicable as natural pesticides in ecological and sustainable agriculture.

Development of a subunit vaccine against the cholangiocarcinoma causing Opisthorchis viverrini: a computational approach

Opisthorchis viverrini is the etiological agent of the disease opisthorchiasis and related cholangiocarcinoma (CCA). It infects fish-eating mammals and more than 10 million people in Southeast Asia suffered from opisthorchiasis with a high fatality rate. The only effective drug against this parasite is Praziquantel, which has significant side effects. Due to the lack of appropriate treatment options and the high death rate, there is a dire need to develop novel therapies against this pathogen. In this study, we designed a multi-epitope chimeric vaccine design against O. viverrini by using immunoinformatics approaches. Non-allergenic and immunogenic MHC-1, MHC-2, and B cell epitopes of three candidate proteins thioredoxin peroxidase (Ov-TPx-1), cathepsin F1 (Ov-CF-1) and calreticulin (Ov-CALR) of O. viverrini, were predicted to construct a potent multiepitope vaccine. The coverage of the HLA-alleles of these selected epitopes was determined globally. Four vaccine constructs made by different adjuvants and linkers were evaluated in the context of their physicochemical properties, antigenicity, and allergenicity. Protein-protein docking and MD simulation found that vaccines 3 was more stable and had a higher binding affinity for TLR2 and TLR4 immune receptors. In-silico restriction cloning of vaccine model led to the formation of plasmid constructs for expression in a suitable host. Finally, the immune simulation showed strong immunological reactions to the engineered vaccine. These findings suggest that the final vaccine construct has the potential to be validated by in vivo and in vitro experiments to confirm its efficacy against the CCA causing O. viverrini.

Temperature effect in the inhibition of PLA2 activity of Bothrops brazili venom by Rosmarinic and Chlorogenic acids, experimental and computational approaches

Myotoxicity caused by snakebite envenoming emerges as one of the main problems of ophidic accidents as it is not well neutralized by the current serum therapy. A promising alternative is to search for efficient small molecule inhibitors that can act against multiple venom components. Phospholipase A2 (PLA2) is frequently found in snake venom and is usually associated with myotoxicity. Thus it represents an excellent target for the search of new treatments. This work reports the effect of temperature in the inhibition of catalytic properties of PLA2 from Bothrops brazili venom by Rosmarinic (RSM) and Chlorogenic (CHL) acids through experimental and computational approaches. Three temperatures were evaluated (25, 37 and 50 °C). In the experimental section, enzymatic assays showed that RSM is a better inhibitor in all three temperatures. At 50 °C, the inhibition efficiency decayed significantly for both acids. Docking studies revealed that both ligands bind to the hydrophobic channel of the protein dimer where the phospholipid binds in the catalytic process, interacting with several functional residues. In this context, RSM presents better interaction energies due to stronger interactions with chain B of the dimer. Molecular dynamics simulations showed that RSM can establish selective interactions with ARG112B of PLA2, which is located next to residues of the putative Membrane Disruption Site in PLA2-like structures. The affinity of RSM and CHL acids towards PLA2 is mainly driven by electrostatic interactions, especially salt bridge interactions established with residues ARG33B (for CHL) and ARG112B (RSM) and hydrogen bonds with residue ASP89A. The inability of CHL to establish a stable interaction with ARG112B was identified as the reason for its lower inhibition efficiency compared to RSM at the three temperatures. Furthermore, extensive structural analysis was performed to explain the lower inhibition efficiency at 50 °C for both ligands. The analysis performed in this work provides important information for the future design of new inhibitors.Communicated by Ramaswamy H. Sarma.

Identification of Prospective Ebola Virus VP35 and VP40 Protein Inhibitors from Myxobacterial Natural Products

The Ebola virus (EBOV) is a lethal pathogen causing hemorrhagic fever syndrome which remains a global health challenge. In the EBOV, two multifunctional proteins, VP35 and VP40, have significant roles in replication, virion assembly, and budding from the cell and have been identified as druggable targets. In this study, we employed in silico methods comprising molecular docking, molecular dynamic simulations, and pharmacological properties to identify prospective drugs for inhibiting VP35 and VP40 proteins from the myxobacterial bioactive natural product repertoire. Cystobactamid 934-2, Cystobactamid 919-1, and Cittilin A bound firmly to VP35. Meanwhile, 2-Hydroxysorangiadenosine, Enhypyrazinone B, and Sorangiadenosine showed strong binding to the matrix protein VP40. Molecular dynamic simulations revealed that, among these compounds, Cystobactamid 919-1 and 2-Hydroxysorangiadenosine had stable interactions with their respective targets. Similarly, molecular mechanics Poisson-Boltzmann surface area (MMPBSA) calculations indicated close-fitting receptor binding with VP35 or VP40. These two compounds also exhibited good pharmacological properties. In conclusion, we identified Cystobactamid 919-1 and 2-Hydroxysorangiadenosine as potential ligands for EBOV that target VP35 and VP40 proteins. These findings signify an essential step in vitro and in vivo to validate their potential for EBOV inhibition.

Drug repurposing for diabetes mellitus: In silico and in vitro investigation of DrugBank database for α-glucosidase inhibitors

The process of developing novel compounds/drugs is arduous, time-intensive, and financially burdensome, characterized by a notably low success rate and relatively high attrition rates. To alleviate these challenges, compound/drug repositioning strategies are employed to predict potential therapeutic effects for DrugBank-approved compounds across various diseases. In this study, we devised a computational and enzyme inhibitory mechanistic approach to identify promising compounds from the pool of DrugBank-approved substances targeting Diabetes Mellitus (DM). Molecular docking analyses were employed to validate the binding interaction patterns and conformations of the screened compounds within the active site of α-glucosidase. Notably, Asp352 and Glu277 participated in interactions within the α-glucosidase-ligand complexes, mediated by conventional hydrogen bonding and van der Waals forces, respectively. The stability of the docked complexes (α-glucosidase-compounds) was scrutinized through Molecular Dynamics (MD) simulations. Subsequent in vitro analyses assessed the therapeutic potential of the repositioned compounds against α-glucosidase. Kinetic studies revealed that "Forodesine" exhibited a lower IC50 (0.24 ± 0.04 mM) compared to the control, and its inhibitory pattern corresponds to that of competitive inhibitors. In-depth in silico secondary structure content analysis detailed the interactions between Forodesine and α-glucosidase, unveiling significant alterations in enzyme conformation upon binding, impacting its catalytic activity. Overall, our findings underscore the potential of Forodesine as a promising candidate for DM treatment through α-glucosidase inhibition. Further validation through in vitro and in vivo studies is imperative to confirm the therapeutic benefits of Forodesine in conformational diseases such as DM.

Mutations of voltage-gated sodium channel contribute to pyrethroid resistance in Panonychus citri

Insecticide resistance in Panonychus citri is a major obstacle to mite control in citrus orchards. Pyrethroid insecticides are continually used to control mites in China, although resistance to pyrethroids has evolved in some populations. Here, the resistance to the pyrethroid fenpropathrin was investigated and 7 out of 8 field-collected populations of P. citri exhibited a high level of resistance, ranging from 171-fold to 15 391-fold higher than the susceptible (SS) comparison strain. Three voltage-gated sodium channel (VGSC) mutations were identified in the tested populations: L1031V, F1747L, and F1751I. Amplicon sequencing was used to evaluate the frequency of these mutations in the 19 field populations. L1031V and F1751I were present in all populations at frequencies of 11.6%-82.1% and 0.5%-31.8%, respectively, whereas the F1747L mutation was only present in 12 populations from Chongqing, Sichuan, Guangxi, and Yunnan provinces. Introduction of these mutations singly or in combination into transgenic flies significantly increased their resistance to fenpropathrin and these flies also exhibited reduced mortality after exposure to the pyrethroids permethrin and β-cypermethrin. Panonychus citri VGSC homology modeling and ligand docking indicate that F1747 and F1751 form direct binding contacts with pyrethroids, which are lost with mutation, whereas L1031 mutation may diminish pyrethroid effects through an allosteric mechanism. Overall, the results provide molecular markers for monitoring pest resistance to pyrethroids and offer new insights into the basis of pyrethroid actions on sodium channels.

Enhancing tuberculosis vaccine development: a deconvolution neural network approach for multi-epitope prediction

Tuberculosis (TB) a disease caused by Mycobacterium tuberculosis (Mtb) poses a significant threat to human life, and current BCG vaccinations only provide sporadic protection, therefore there is a need for developing efficient vaccines. Numerous immunoinformatic methods have been utilized previously, here for the first time a deep learning framework based on Deconvolutional Neural Networks (DCNN) and Bidirectional Long Short-Term Memory (DCNN-BiLSTM) was used to predict Mtb Multiepitope vaccine (MtbMEV) subunits against six Mtb H37Rv proteins. The trained model was used to design MEV within a few minutes against TB better than other machine learning models with 99.5% accuracy. The MEV has good antigenicity, and physiochemical properties, and is thermostable, soluble, and hydrophilic. The vaccine's BLAST search ruled out the possibility of autoimmune reactions. The secondary structure analysis revealed 87% coil, 10% beta, and 2% alpha helix, while the tertiary structure was highly upgraded after refinement. Molecular docking with TLR3 and TLR4 receptors showed good binding, indicating high immune reactions. Immune response simulation confirmed the generation of innate and adaptive responses. In-silico cloning revealed the vaccine is highly expressed in E. coli. The results can be further experimentally verified using various analyses to establish a candidate vaccine for future clinical trials.

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

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

Computational analysis of non-synonymous single nucleotide polymorphism in the bovine PKLR geneComputational analysis of bovine PKLR gene

Pyruvate kinase (PKLR) is a potential candidate gene for milk production traits in cows. The main aim of this work is to investigate the potentially deleterious non-synonymous single nucleotide polymorphisms (nsSNPs) in the PKLR gene by using several computational tools. In silico tools including SIFT, Polyphen-2, SNAP2 and Panther indicated only 18 nsSNPs out of 170 were considered deleterious. The analysis of proteins' stability change due to amino acid substitution performed by the use of the I-mutant, MUpro, CUPSTAT, SDM and Dynamut confirmed that 9 nsSNPs decreased protein stability. ConSurf analysis predicted that all 18 nsSNPs were evolutionary moderately or highly conserved. Two different domains of PKLR protein were revealed by the InterPro tool with 12 nsSNPs positioned in the Pyruvate Kinase barrel domain and 6 nsSNP present in the Pyruvate Kinase C Terminal. The PKLR 3D model was predicted by MODELLER software and validated via Ramachandran plot and Prosa which indicated a good quality model. The analysis of energy minimizations for the native and mutated structures was performed by SWISS PDB viewer with GROMOS 96 program and showed that 3 structural and 4 functional residues had total energy higher than the native model. These findings indicate that these mutant structures (rs441424814, rs449326723, rs476805413, rs472263384, rs474320860, rs475521477, rs441633284) were less stable than the native model. Molecular Dynamics simulations were performed to confirm the impact of nsSNPs on the protein structure and function. The present study provides useful information about functional SNPs that have an impact on PKLR protein in cattle.Communicated by Ramaswamy H. Sarma.

Homology modeling and molecular docking studies to decrease glutamine affinity of Yarrowia lipolytica L-asparaginase

L-Asparaginase is a key component in the treatment of leukemias and lymphomas. However, the glutamine affinity of this therapeutic enzyme is an off-target activity that causes several side effects. The modeling and molecular docking study of Yarrowia lipolytica L-asparaginase (YL-ASNase) to reduce its l-glutamine affinity and increase its stability was the aim of this study. Protein-ligand interactions of wild-type and different mutants of YL-ASNase against L-asparagine compared to l-glutamine were assessed using AutoDock Vina tools because the crystal structure of YL-ASNase does not exist in the protein data banks. The results showed that three mutants, T171S, T171S-N60A, and T171A-T223A, caused a considerable increase in L-asparagine affinity and a decrease in l-glutamine affinity as compared to the wild-type and other mutants. Then, molecular dynamics simulation and MM/GBSA free energy were applied to assess the stability of protein structure and its interaction with ligands. The three mutated proteins, especially T171S-N60A, had higher stability and interactions with L-asparagine than l-glutamine in comparison with the wild-type. The YL-ASNase mutants could be introduced as appropriate therapeutic candidates that might cause lower side effects. However, the functional properties of these mutated enzymes need to be confirmed by genetic manipulation and in vitro and in vivo studies.

Identification of Bile Acid-Derived Chemical Chaperone(s) Targeting E46K-Mutated Alpha-Synuclein Protein to Treat Parkinson's Disease: Molecular Modelling, Docking, ADME, and Simulation Studies

Aggregated α-synuclein (α-syn) present inside small cytoplasmic inclusions in the substantia nigra region marks the major pathological hallmark of Parkinson's disease (PD) and makes it an attractive target for the drug development process. Certain small-molecule chaperones (such as DCA, UDCA, TUDCA) presented the ability to prevent misfolding and aggregation of α-syn as well as to disentangle mature α-syn amyloid fibrils. However, due to toxicity constraints, these small molecules could not be translated into clinical settings. Computational biology methods and bioinformatics approaches allow virtual screening of a large number of molecules, with reduced side effects and better efficacy. In the present study, a library of 10,928 derivatives was generated using DCA, UDCA, and TUDCA bile acid scaffolds and analysed for their binding affinity, pharmacokinetic properties, and drug likeliness profile, to come up with promising compounds with reduced toxicity and better chaperone ability. Molecular docking revealed that with respect to their free binding energy, C1-C25 have the lowest binding energy and bind significantly to recombinantly assembled E46K α-syn fibrils (PDB ID-6UFR). In silico ADME predictions revealed that all these compounds had minimal toxic effects and had good absorption as well as solubility characteristics. Simulation studies further showed that the imidazole ring-based TUDCA derivatives interacted better with the protein in comparison to the others. The proposed study has identified potent chemical chaperones (C2 and C3) as effective therapeutic agents for Parkinson's disease, and further in vitro and in vivo testing will be undertaken to substantiate their potential as novel drugs.

Structure and dynamics of the cyanobacterial regulator SipA

The small, 78-residue long, regulator SipA interacts with the non-bleaching sensor histidine kinase (NblS). We have solved the solution structure of SipA on the basis of 990 nuclear Overhauser effect- (NOE-) derived distance constraints. The average pairwise root-mean-square deviation (RMSD) for the twenty best structures for the backbone residues, obtained by CYANA, was 1.35 ± 0.21 Å, and 1.90 ± 0.16 Å when all heavy atoms were considered (the target function of CYANA was 0.540 ± 0.08). The structure is that of a β-II class protein, basically formed by a five-stranded β-sheet composed of antiparallel strands following the arrangement: Gly6-Leu11 (β-strand 1), which packs against Leu66-Val69 (β-strand 5) on one side, and against Gly36-Thr42 (β-strand 2) on the other side; Trp50-Phe54 (β-strand 3); and Gly57-Leu60 (β-strand 4). The protein is highly mobile, as shown by measurements of R1, R2, NOE and ηxy relaxation parameters, with an average order parameter () of 0.70; this mobility encompasses movements in different time scales. We hypothesize that this high flexibility allows the interaction with other proteins (among them NblS), and it explains the large conformational stability of SipA.

In silico characterization of the novel SDR42E1 as a potential vitamin D modulator

The short-chain dehydrogenase/reductase (SDR) superfamily encompasses enzymes that play essential roles in the metabolism of steroid hormones and lipids. Despite an enigmatic function, recent genetic studies have linked the novel SDR 42 extended-1 (SDR42E1) gene to 25-hydroxyvitamin D levels. This study investigated the potential SDR42E1 functions and interactions with vitamin D using bioinformatics and molecular docking studies. Phylogenetic analysis unveiled that the nucleotide sequences of human SDR42E1 exhibit high evolutionary conservation across nematodes and fruit flies. Molecular docking analysis identified strong binding affinities between SDR42E1 and its orthologs with vitamin D3 and essential precursors, 8-dehydrocholesterol, followed by 7-dehydrocholesterol and 25-hydroxyvitamin D. The hydrophobic interactions observed between the protein residues and vitamin D compounds supported the predicted transmembrane localization of SDR42E1. Our investigation provides valuable insights into the potential role of SDR42E1 in skin vitamin D biosynthesis throughout species. This provides the foundation for future research and development of targeted therapies for vitamin D deficiency and related health conditions.

XRCC1 Polymorphisms p.Arg194Trp, p.Arg280His, and p.Arg399Gln, Polycyclic Aromatic Hydrocarbons, and Infertility: A Case-Control and In Silico Study

XRCC1 is involved in repair of single-strand breaks generated by mutagenic exposure. Polymorphisms within XRCC1 affect its ability to efficiently repair DNA damage. Polycyclic aromatic hydrocarbons or PAHs are genotoxic compounds which form bulky DNA adducts that are linked with infertility. Few reports suggest combined role of XRCC1 polymorphisms and PAHs in infertility. Present study investigates association of three XRCC1 polymorphisms (p.Arg194Trp, p.Arg280His, p.Arg399Gln) with male and female infertility in a North-West Indian population using case-control approach. Additionally, in silico approach has been used to predict whether XRCC1 polymorphisms effect interaction of XRCC1 with different PAHs. For case-control study, XRCC1 polymorphisms were screened in peripheral blood samples of age- and gender-matched 201 infertile cases (♂-100, ♀-101) and 201 fertile controls (♂-100, ♀-101) using PCR-RFLP method. For in silico study, AutoDock v4.2.6 was used for molecular docking of B[a]P, BPDE-I, ( ±)-anti-BPDE, DB[a,l]P, 1-N, 2-N, 1-OHP, 2-OHF with XRCC1 and assess effect of XRCC1 polymorphisms on their interaction. In case-control study, statistical analysis showed association of XRCC1 p.Arg280His GA genotype (p = 0.027), A allele (p = 0.019) with reduced risk of male infertility. XRCC1 p.Arg399Gln AA genotype (p = 0.021), A allele (p = 0.014) were associated with reduced risk for female primary infertility. XRCC1 p.Arg194Trp T allele was associated with increased risk for female infertility (p = 0.035). In silico analysis showed XRCC1-PAH interaction with non-significant effect of XRCC1 polymorphisms on predicted binding. Therefore, present study concludes that XRCC1 polymorphism-modified risk for male and female infertility in North-West Indians without significant effect on predicted XRCC1-PAH interactions. This is the first report on XRCC1 in female infertility.

Characterization of a CrPME indicates its possible role in determining vindoline accumulation in Catharanthus roseus leaves

The leaf-specific Catharanthus roseus alkaloid, vindoline, is the major bottleneck precursor in the production of scarce and costly anticancer bisindoles (vincristine and vinblastine). The final steps of its biosynthesis and storage occur in the laticifers. Earlier, we have shown that vindoline content is directly related to laticifer number. Pectin remodeling enzymes, like pectin methylesterase (PME), are known to be involved in laticifer development. A search in the croFGD yielded a leaf-abundant CrPME isoform that co-expressed with a few vindoline biosynthetic genes. Full-length cloning, tissue-specific expression profiling, and in silico analysis of CrPME were carried out. It was found to possess all the specific characteristics of a typical plant PME. Transient silencing (through VIGS) and overexpression of CrPME in C. roseus indicated a direct relationship between its expression and vindoline content. Comparative analysis of transcript abundance and enzyme activity in three familial C. roseus genotypes differing significantly in their vindoline content and laticifer count (CIM-Sushil > Dhawal > Nirmal) also corroborated the positive relationship of CrPME expression with vindoline content. This study highlights the possible role of CrPME, a cell wall remodeling enzyme, in modulating laticifer-associated secondary metabolism.

Mispacking of the F87 sidechain drives aggregation-promoting conformational fluctuations in the subunit interfaces of the transthyretin tetramer

Aberrant formation and deposition of human transthyretin (TTR) aggregates causes transthyretin amyloidosis. To initialize aggregation, transthyretin tetramers must first dissociate into monomers that partially unfold to promote entry into the aggregation pathway. The native TTR tetramer (T) is stabilized by docking of the F87 sidechain into an interfacial cavity enclosed by several hydrophobic residues including A120. We have previously shown that an alternative tetramer (T*) with mispacked F87 sidechains is more prone to dissociation and aggregation than the native T state. However, the molecular basis for the reduced stability in T* remains unclear. Here we report characterization of the A120L mutant, where steric hindrance is introduced into the F87 binding site. The X-ray structure of A120L shows that the F87 sidechain is displaced from its docking site across the subunit interface. In A120S, a naturally occurring pathogenic mutant that is less aggregation-prone than A120L, the F87 sidechain is correctly docked, as in the native TTR tetramer. Nevertheless, 19F-NMR aggregation assays show an elevated population of a monomeric aggregation intermediate in A120S relative to a control containing the native A120, due to accelerated tetramer dissociation and slowed monomer tetramerization. The mispacking of the F87 sidechain is associated with enhanced exchange dynamics for interfacial residues. At 298 K, the T* populations of various naturally occurring mutants fall between 4-7% (ΔG ~ 1.5-1.9 kcal/mol), consistent with the free energy change expected for undocking and solvent exposure of one of the four F87 sidechains in the tetramer (ΔG ~ 1.6 kcal/mol). Our data provide a molecular-level picture of the likely universal F87 sidechain mispacking in tetrameric TTR that promotes interfacial conformational dynamics and increases aggregation propensity.

Site-specific mutagenesis on Mnemiopsin 2: Calcium coordination and substrate binding properties of new variants

We used site-specific mutagenesis by targeting E179 and F190 on the structure of photoprotein Mnemiopsin 2 (Mn2) from Mnemiopsis leidyi. The tertiary structure of E179S and F190L mutants was made by the MODELLER program. Far-ultraviolet circular dichroism data showed that the overall secondary structural content of photoprotein is not changed upon mutation, however the helicity and stabilizing interactions in helical structure decreases in mutants as compared with the wild-type (WT) photoprotein. Fluorescence spectra data revealed that the tertiary structure of the mutants is more compact than that of WT Mn2. According to the heat-induced denaturation experiments data, the melting temperature (Tm ) for the unfolding of tertiary structure of the F190L variant increases by 3°C compared with that of the WT and E179S mutant. Interestingly, the conformational enthalpy of the F190L mutant (86 kcal mol-1 ) is considerably lower than those in the WT photoprotein (102 kcal mol-1 ) and E179S mutant (106 kcal mol-1 ). The significant difference in the enthalpy of the thermal unfolding process could be explained by considering that the thermally denatured state of the F190L mutant is structurally less expanded than the WT and E179S variants. Bioluminescence activity data showed that the maximum characteristic wavelengths of the mutants undergo blue shift as compared with the WT protein. Initial intensity of the F190L and E179S variants was recorded to be 137.5% and 55.9% of the WT protein, respectively.

Understanding a point mutation signature D54K in the caspase activation recruitment domain of NOD1 capitulating concerted immunity via atomistic simulation

Point mutation D54K in the human N-terminal caspase recruitment domain (CARD) of nucleotide-binding oligomerization domain -1 (NOD1) abrogates an imperative downstream interaction with receptor-interacting protein kinase (RIPK2) that entails combating bacterial infections and inflammatory dysfunction. Here, we addressed the molecular details concerning conformational changes and interaction patterns (monomeric-dimeric states) of D54K by signature-based molecular dynamics simulation. Initially, the sequence analysis prioritized D54K as a pathogenic mutation, among other variants, based on a sequence signature. Since the mutation is highly conserved, we derived the distant ortholog to predict the sequence and structural similarity between native and mutant. This analysis showed the utility of 33 communal core residues associated with structural-functional preservation and variations, concurrently served to infer the cryptic hotspots Cys39, Glu53, Asp54, Glu56, Ile57, Leu74, and Lys78 determining the inter helical fold forming homodimers for putative receptor interaction. Subsequently, the atomistic simulations with free energy (MM/PB(GB)SA) calculations predicted structural alteration that takes place in the N-terminal mutant CARD where coils changed to helices (45 α3- L4-α4-L6- α683) in contrast to native (45T2-L4-α4-L6-T483). Likewise, the C-terminal helices 93T1-α7105 connected to the loops distorted compared to native 93α6-L7105 may result in conformational misfolding that promotes functional regulation and activation. These structural perturbations of D54K possibly destabilize the flexible adaptation of critical homotypic NOD1CARD-CARDRIPK2 interactions (α4Asp42-Arg488α5 and α6Phe86-Lys471α4) is consistent with earlier experimental reports. Altogether, our findings unveil the conformational plasticity of mutation-dependent immunomodulatory response and may aid in functional validation exploring clinical investigation on CARD-regulated immunotherapies to prevent systemic infection and inflammation.Communicated by Ramaswamy H. Sarma.

Elucidation of the structural dynamics of mutations in PHB2 protein associated with growth suppression and cancer progression

Prohibitin is a multifunctional protein that plays an important role in numerous cellular processes. Membrane-associated mitochondrial prohibitin complex is made up of two subunits, PHB1 and PHB2 which are ubiquitously expressed and analogous to each other. High levels of prohibitin expression have consequently been found in esophageal cancer, endometrial adenocarcinoma, gastric cancer, hepatocellular carcinoma, breast cancer and bladder cancer. The aim of this study is to analyse two-point mutation PHB2_MT1(I → A) and PHB2_MT2(I → P), their effect on PHB2 protein and its effect on formation of mitochondrial complex. It is a residual level study, based on current experimental validation. To establish the effects of the two-point mutations, computational approaches such as molecular modelling, molecular docking, normal mode simulation, molecular dynamics simulations and MM/GBSA were used. An analysis of the energy dynamics of both unbound and complex proteins was conducted to elucidate how mutations impact the energy distribution of PHB2. Our study confirmed that the two mutations decreased the overall stability of PHB2. This was evidenced by heightened atomic fluctuations within the mutated region, accompanied by elevated deviations observed in RMSD and Rg values. Furthermore, these mutations were correlated with a decline in the organization of secondary structural elements. The mutations in PHB2_MT1 and PHB2_MT2 resulted in formation a less stable prohibitin complex. Thus, PHB1 and PHB2 may act as molecular target or novel biomarkers for therapeutic intervention in numerous forms of malignancies.

Exploring the inhibitory properties of biflavonoids on α-glucosidase; computational and experimental approaches

Biflavonoids (BFs) are a group of polyphenols that have a unique biochemical structure. One of the key biomedical mechanisms that BFs can have high potential in managing Diabetes mellitus (DM) is α-glucosidase inhibition. Normally, elevated blood glucose levels are caused by high absorption of glucose in the epithelium of the small intestine. Since α-glucosidase helps increase the absorption of glucose in the small intestine in the final stage of glycan catabolism, inhibition of this essential biochemical process in diabetic patients can be considered a suitable approach in the treatment of this disease. The interaction between the BFs and α-glucosidase are still not clear, and need to be deeply investigated. Herein, the aim is to identify BFs with strong α-glucosidase inhibitory activity. Using docking-based virtual screening approach, the potential binding affinity of 18 selected BFs to α-glucosidase was evaluated. The dynamic activity and stability of α-glucosidase-BFs complexes were then measured by molecular dynamics simulation (MDs). "Strychnobiflavone" showed the best score in α-glucosidase inhibition. Arg315 and Phe303 involved in the interactions of α-glucosidase-strychnobiflavone complex through cation-π and π-π stacking, respectively. Based on in vitro kinetic studies, it was determined that the type of inhibition of "strychnobiflavone" corresponds to the pattern of mixed inhibitors. Furthermore, details of the interactions between strychnobiflavone and α-glucosidase were performed by in silico secondary structure content analysis. The findings showed when "strychnobifone" binds to the enzyme, significant alterations occur in the enzyme conformation affecting its catalytic activity. In general, the findings highlighted the potential of "strychnobiflavone" as a promising candidate for the treatment of diabetes mellitus through α-glucosidase inhibition. Further in vitro and in vivo studies have to confirm the therapeutic benefits of "strychnobiflavone" in conformational diseases such as diabetes mellitus.

Molecular characterization of glutor-GLUT interaction and prediction of glutor's drug-likeness: implications for its utility as an antineoplastic agent

Recent experimental evidence from our and other laboratories has strongly indicated that glutor, a piperazine-2-one derivative, which is a pan-GLUT inhibitor, displays a promising antineoplastic action by hampering glucose uptake owing to its ability to inhibit GLUT1 and GLUT3, which are overexpressed in neoplastic cells. However, the molecular mechanism(s) of the inhibiting action of glutor has remained elusive. Thus, for optimal utilization of the antineoplastic potential of glutor, it is essential to decipher the precise mechanism(s) of its interaction with GLUTs. Therefore, the present investigation was carried out to understand the molecular mechanism(s) of the binding of glutor to GLUT1 and GLUT3 in silico. This study suggests that glutor can effectively bind to GLUTs at the reported binding site. Moreover, the docking of glutor to GLUT was stabilised by several contacts between these two partners as shown by the 200 ns long molecular dynamic simulation carried out using Gromacs, indicating the formation of a stable complex. Moreover, glutor was found to possess all characteristics conducive to its drug-likeness. Hence, these observations suggest that glutor has the potential to be used in antineoplastic therapeutic applications.Communicated by Ramaswamy H. Sarma.

Biochemical, structural and dynamical characterizations of the lactate dehydrogenase from Selenomonas ruminantium provide information about an intermediate evolutionary step prior to complete allosteric regulation acquisition in the super family of lactate and malate dehydrogenases

In this work, we investigated the lactate dehydrogenase (LDH) from Selenomonas ruminantium (S. rum), an enzyme that differs at key amino acid positions from canonical allosteric LDHs. The wild type (Wt) of this enzyme recognises pyuvate as all LDHs. However, introducing a single point mutation in the active site loop (I85R) allows S. Rum LDH to recognize the oxaloacetate substrate as a typical malate dehydrogenase (MalDH), whilst maintaining homotropic activation as an LDH. We report the tertiary structure of the Wt and I85RLDH mutant. The Wt S. rum enzyme structure binds NADH and malonate, whilst also resembling the typical compact R-active state of canonical LDHs. The structure of the mutant with I85R was solved in the Apo State (without ligand), and shows no large conformational reorganization such as that observed with canonical allosteric LDHs in Apo state. This is due to a local structural feature typical of S. rum LDH that prevents large-scale conformational reorganization. The S. rum LDH was also studied using Molecular Dynamics simulations, probing specific local deformations of the active site that allow the S. rum LDH to sample the T-inactive state. We propose that, with respect to the LDH/MalDH superfamily, the S. rum enzyme possesses a specificstructural and dynamical way to ensure homotropic activation.

Exploring Flexibility and Folding Patterns Throughout Time in Voltage Sensors

The voltage-sensing domain (VSD) is a module capable of responding to changes in the membrane potential through conformational changes and facilitating electromechanical coupling to open a pore gate, activate proton permeation pathways, or promote enzymatic activity in some membrane-anchored phosphatases. To carry out these functions, this module acts cooperatively through conformational changes. The VSD is formed by four transmembrane segments (S1-S4) but the S4 segment is critical since it carries positively charged residues, mainly Arg or Lys, which require an aqueous environment for its proper function. The discovery of this module in voltage-gated ion channels (VGICs), proton channels (Hv1), and voltage sensor-containing phosphatases (VSPs) has expanded our understanding of the principle of modularity in the voltage-sensing mechanism of these proteins. Here, by sequence comparison and the evaluation of the relationship between sequence composition, intrinsic flexibility, and structural analysis in 14 selected representatives of these three major protein groups, we report five interesting differences in the folding patterns of the VSD both in prokaryotes and eukaryotes. Our main findings indicate that this module is highly conserved throughout the evolutionary scale, however: (1) segments S1 to S3 in eukaryotes are significantly more hydrophobic than those present in prokaryotes; (2) the S4 segment has retained its hydrophilic character; (3) in eukaryotes the extramembranous linkers are significantly larger and more flexible in comparison with those present in prokaryotes; (4) the sensors present in the kHv1 proton channel and the ciVSP phosphatase, both of eukaryotic origin, exhibit relationships of flexibility and folding patterns very close to the typical ones found in prokaryotic voltage sensors; and (5) archaeal channels KvAP and MVP have flexibility profiles which are clearly contrasting in the S3-S4 region, which could explain their divergent activation mechanisms. Finally, to elucidate the obscure origins of this module, we show further evidence for a possible connection between voltage sensors and TolQ proteins.

Structure-based virtual screening, ADMET analysis, and molecular dynamics simulation of Moroccan natural compounds as candidates for the SARS-CoV-2 inhibitors

The lack of treatments and vaccines effective against SARS-CoV-2 has forced us to explore natural compounds that could potentially inhibit this virus. Additionally, Morocco is renowned for its rich plant diversity and traditional medicinal uses, which inspires us to leverage our cultural heritage and the abundance of natural resources in our country for therapeutic purposes. In this study, an extensive investigation was conducted to gather a collection of phytoconstituents extracted from Moroccan plants, aiming to evaluate their ability to inhibit the proliferation of the SARS-CoV-2 virus. Molecular docking of the studied compounds was performed at the active sites of the main protease (6lu7) and spike (6m0j) proteins to assess their binding affinity to these target proteins. Compounds exhibiting high affinity to the proteins underwent further evaluation based on Lipinski's rule and ADME-Tox analysis to gain insights into their oral bioavailability and safety. The results revealed that the two compounds demonstrated strong binding affinity to the target proteins, making them potential candidates for oral antiviral drugs against SARS-CoV-2. The molecular dynamics results from this computational analysis supported the overall stability of the resulting complex.

Structural informatics approach for designing an epitope-based vaccine against the brain-eating Naegleria fowleri

Primary Amoebic Meningoencephalitis (PAM), a severe lethal brain disease, is caused by a parasite, Naegleria fowleri, also known as the "brain-eating amoeba". The chances of a patient's recovery after being affected by this parasite are very low. Only 5% of people are known to survive this life-threatening infection. Despite the fact that N. fowleri causes a severe, fatal infection, there is no proper treatment available to prevent or cure it. In this context, it is necessary to formulate a potential vaccine that could be able to combat N. fowleri infection. The current study aimed at developing a multi-epitope subunit vaccine against N. fowleri by utilizing immunoinformatics techniques and reverse vaccinology approaches. The T- and B-cell epitopes were predicted by various tools. In order to choose epitopes with the ability to trigger both T- and B-cell-mediated immune responses, the epitopes were put through a screening pipeline including toxicity, antigenicity, cytokine-inductivity, and allergenicity analysis. Three vaccine constructs were designed from the generated epitopes linked with linkers and adjuvants. The modeled vaccines were docked with the immune receptors, where vaccine-1 showed the highest binding affinity. Binding affinity and stability of the docked complex were confirmed through normal mode analysis and molecular dynamic simulations. Immune simulations developed the immune profile, and in silico cloning affirmed the expression probability of the vaccine construct in Escherichia coli (E. coli) strain K12. This study demonstrates an innovative preventative strategy for the brain-eating amoeba by developing a potential vaccine through immunoinformatics and reverse vaccinology approaches. This study has great preventive potential for Primary Amoebic Meningoencephalitis, and further research is required to assess the efficacy of the designed vaccine.

Computational prognostic evaluation of Alzheimer's drugs from FDA-approved database through structural conformational dynamics and drug repositioning approaches

Drug designing is high-priced and time taking process with low success rate. To overcome this obligation, computational drug repositioning technique is being promptly used to predict the possible therapeutic effects of FDA approved drugs against multiple diseases. In this computational study, protein modeling, shape-based screening, molecular docking, pharmacogenomics, and molecular dynamic simulation approaches have been utilized to retrieve the FDA approved drugs against AD. The predicted MADD protein structure was designed by homology modeling and characterized through different computational resources. Donepezil and galantamine were implanted as standard drugs and drugs were screened out based on structural similarities. Furthermore, these drugs were evaluated and based on binding energy (Kcal/mol) profiles against MADD through PyRx tool. Moreover, pharmacogenomics analysis showed good possible associations with AD mediated genes and confirmed through detail literature survey. The best 6 drug (darifenacin, astemizole, tubocurarine, elacridar, sertindole and tariquidar) further docked and analyzed their interaction behavior through hydrogen binding. Finally, MD simulation study were carried out on these drugs and evaluated their stability behavior by generating root mean square deviation and fluctuations (RMSD/F), radius of gyration (Rg) and soluble accessible surface area (SASA) graphs. Taken together, darifenacin, astemizole, tubocurarine, elacridar, sertindole and tariquidar displayed good lead like profile as compared with standard and can be used as possible therapeutic agent in the treatment of AD after in-vitro and in-vivo assessment.

Geraniol Suppresses Oxidative Stress, Inflammation, and Interstitial Collagenase to Protect against Inflammatory Arthritis

Geraniol (GER) is a plant-derived acyclic isoprenoid monoterpene that has displayed anti-inflammatory effects in numerous in vivo and in vitro models. This study was therefore designed to evaluate the antiarthritic potential of GER in complete Freund's adjuvant (CFA)-induced inflammatory arthritis (IA) model in rats. IA was induced by intraplantar injection of CFA (0.1 mL), and a week after CFA administration, rats were treated with various doses of methotrexate (MTX; 1 mg/kg) or GER (25, 50, and 100 mg/kg). Treatments were given on every alternate day, and animals were sacrificed on the 35th day. Paw volume, histopathological, hematological, radiographic, and qPCR analyses were performed to analyze the severity of the disease. GER significantly reduced paw edema after 35 days of treatment, and these results were comparable to the MTX-treated group. GER-treated animals displayed a perfect joint structure with minimal inflammation and no signs of cartilage or bone damage. Moreover, GER restored red blood cell and hemoglobin levels, normalized erythrocyte sedimentation rate, platelet, and c-reactive protein values, and also attenuated the levels of rheumatoid factor. RT-qPCR analysis demonstrated that GER decreased mRNA expression of pro-inflammatory cytokines like tumor necrosis factor-alpha (TNF-α) and interleukin-1 beta. GER also down-regulated the transcript levels of cyclooxygenase-2 (COX-2), microsomal prostaglandin E synthase-1, prostaglandin D2 synthase, and interstitial collagenase (MMP-1). Molecular docking of GER with COX-2, TNF-α, and MMP-1 also revealed that the antiarthritic effects of GER could be due to its direct interactions with these mediators. Based on our findings, it is conceivable that the antiarthritic effects of GER could be attributed to downregulation of pro-inflammatory mediators and protease like MMP-1.

New insights into comprehensive analysis of magnesium transporter (MGT) gene family in rice (Oryza sativa L.)

Magnesium transporters (MGTs) regulate magnesium absorption, transport, and redistribution in higher plants. To investigate the role of the Oryza sativa MGTs gene family members under salt stress, this study analyzed the protein properties, gene structure, phylogenetic relationship, synteny patterns, expression, and co-expression networks of 23 non-redundant OsMGT. The evolutionary relationship of the OsMGT gene family was fully consistent with their functional domain, and were divided into three main classes based on the conserved domain: MMgT, CorA-like, and NIPA. The α/β patterns in the protein structures were highly similar in the CorA-like and NIPA members, with the conserved structures in the Mg2+-binding and catalytic regions. The CorA-like clade-related proteins demonstrated the highest numbers of protein channels with Pro, Ser, Lys, Gly, and Tyr, as the critical binding residues. The expression analysis of OsMGT genes in various tissues showed that MGTs' gene family may possess critical functions during rice development. Gene expression analysis of candidate OsMGT using reverse-transcription quantitative real-time PCR (RT-qPCR) found that four OsMGT genes exhibited different expression patterns in salt-sensitive and salt-tolerant rice genotypes. We hypothesize that the OsMGT gene family members may be involved in responses to salt stress. These findings could be useful for further functional investigation of MGTs as well as defining their involvement in abiotic stress studies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03735-4.

Synthesis, Computational Studies, Antioxidant and Anti-Inflammatory Bio-Evaluation of 2,5-Disubstituted-1,3,4-Oxadiazole Derivatives

The 1,3,4-oxadiazole derivatives Ox-6a-f have been synthesized by incorporating flurbiprofen moiety with the aim to explore the potential of target molecules to decrease the oxidative stress. The title compounds Ox-6a-f were prepared by simple reactions in which a flurbiprofen -COOH group was esterified with methanol in an acid-catalyzed medium, which was then reacted with hydrazine to afford the corresponding hydrazide. The acid hydrazide was then cyclized into 1,3,4-oxadiazole-2-thiol by reacting with CS₂ in the presence of KOH. The title compounds Ox-6a-f were synthesized by the reaction of an -SH group with various alkyl/aryl chlorides, which involves an S-alkylation reaction. The structures of the synthesized Ox-6a-f derivatives were ascertained by spectroscopic data. The in silico molecular docking was performed against target proteins cyclooxygenase-2 COX-2 (PDBID 5KIR) and cyclooxygenase-1 COX-1 (PDBID 6Y3C) to determine the binding affinity of the synthesized compounds with these structures. It has been inferred that most of the synthesized compounds bind well with an active binding site of 5KIR compared to 6Y3C, and especially compound Ox-6f showed excellent binding affinity (7.70 kcal/mol) among all synthesized compounds Ox-6a-f. The molecular dynamic (MD) simulation has also been performed to check the stability of docking complexes of ligands with COX-2 by determining their root mean square deviation and root mean square fluctuation. Little fluctuation was observed in case of Ox-6f, which forms the most stable complex with COX-2. The comprehensive antioxidant potential of the synthesized compounds has been evaluated by determining their free radical scavenging activity, including DPPH, OH, nitric oxide (NO), and iron chelation assay. The derivative Ox-6f showed promising results with 80.23% radical scavenging potential at a dose of 100 µg/mL while ascorbic acid exhibited 87.72% inhibition at the same dose. The anti-inflammatory activity of the final products has also been performed, and inflammatory markers were assayed, such as a thiobarbituric acid-reducing substance, nitric oxide, interleukin-6 (IL-6), and COX-2. The derivatives Ox-6d and Ox-6f displayed higher anti-inflammatory activity, exhibiting 70.56% and 74.16% activity, respectively. The results were compared with standard ibuprofen, which showed 84.31% activity at the same dose, 200 µg/mL. The anti-inflammatory potential has been performed by following the carrageen-induced hind paw edema model, and results showed that derivative Ox-6f exhibited 79.83% reduction in edema volume compared to standard ibuprofen, which reduced 84.31% edema volume. As dry lab and wet lab results confirm each other, it has been deduced that derivative Ox-6f may serve as the lead structure to design potent compounds to address oxidative stress.

Molecular evolution steered structural adaptations in the DNA polymerase III α subunit of halophilic bacterium Salinibacter ruber

A significant portion of the earth has a salty environment, and the literature on bacterial survival mechanisms in salty environments is limited. During molecular evolution, halophiles increase acidic amino acid residues on their protein surfaces which leads to a negatively charged surface potential that helps them to maintain the protein integrity and protect them from denaturation by competing with salt ions. Through protein family analysis, we have investigated the molecular-level adaptive features of DNA polymerase III's catalytic subunit (alpha) and its structure-function relationship. This study throws light on the novel understanding of halophilic bacterial replication and the molecular basis of salt adaptation. Comparisons of the amino acid contents and electronegativity of halophilic and mesophilic bacterial proteins revealed adaptations that allow halophilic bacteria to thrive in high salt concentrations. A significantly lower isoelectric point of halophilic bacterial proteins indicates the acidic nature. Also, an abundance of disordered regions in halophiles suggests the requirement of the salt ions that play a crucial role in their stable protein folding. Despite having similar topology, mesophilic and halophilic proteins, a set of very prominent molecular modifications was observed in the alpha subunit of halophiles.

Design of ion channel blocking, toxin-like Kunitz inhibitor peptides from the tapeworm, Echinococcus granulosus, with potential anti-cancer activity

Over-expression of K+ channels has been reported in human cancers and is associated with the poor prognosis of several malignancies. EAG1, a particular potassium ion channel, is widely expressed in the brain but poorly expressed in other normal tissues. Kunitz proteins are dominant in metazoan including the dog tapeworm, Echinococcus granulosus. Using computational analyses on one A-type potassium channel, EAG1, and in vitro cellular methods, including major cancer cell biomarkers expression, immunocytochemistry and whole-cell patch clamp, we demonstrated the anti-tumor activity of three synthetic small peptides derived from E. granulosus Kunitz4 protease inhibitors. Experiments showed induced significant apoptosis and inhibition of proliferation in both cancer cell lines via disruption in cell-cycle transition from the G0/G1 to S phase. Western blotting showed that the levels of cell cycle-related proteins including P27 and P53 were altered upon kunitz4-a and kunitz4-c treatment. Patch clamp analysis demonstrated a significant increase in spontaneous firing frequency in Purkinje neurons, and exposure to kunitz4-c was associated with an increase in the number of rebound action potentials after hyperpolarized current. This noteworthy component in nature could act as an ion channel blocker and is a potential candidate for cancer chemotherapy based on potassium channel blockage.

An in-silico study to understand the effect of lineage diversity on cold shock response: unveiling protein-RNA interactions among paralogous CSPs of E. coli

Cold shock proteins (CSPs) are small, cytoplasmic, ubiquitous and acidic proteins. They have a single nucleic acid-binding domain and pose as "RNA chaperones" by binding to ssRNA in a low sequence specificity and cooperative manner. They are found in a family of nine homologous CSPs in E. coli. CspA, CspB, CspG and CspI are immensely cold inducible, CspE and CspC are consistently released at usual physiological temperatures and CspD is also induced under nutrient stress. The paralogous protein pairs CSPA/CSPB, CSPC/CSPE, CSPG/CSPI and CSPF/CSPH were first identified. The eight proteins were subjected to molecular modelling and simulation to obtain the most stable conformation in correspondence to their equilibrated RMSD and RMSF graph. The results were compared and it was observed that CSPB, CSPE, CSPF and CSPI were more stable than their paralogous partner conforming to their near equilibrated RMSD curve and low fluctuating RMSF graph. The paralogous proteins were docked with ssRNA and simultaneously binding affinity, interaction types, electrostatic surface potential, hydrophobicity, conformational analysis and SASA were calculated to minutely study and understand the molecular mechanism initiated by these proteins. It was found that CSPB, CSPC, CSPH and CSPI displayed higher affinity towards ssRNA than their paralogous partner. The results further corroborated with ΔGmmgbsa and ΔGfold energy. Between the paralogous pairs CSPC, CSPH and CSPI exhibited higher binding free energy than their partner. Further, CSPB, CSPC and CSPI exhibited higher folding free energy than their paralogous pair. CSPH exhibited highest ΔGmmgbsa of - 522.2 kcal/mol and lowest was displayed by CSPG of around - 309.3 kcal/mol. Highest number of mutations were recognised in CSPF/CSPH and CSPG/CSPI pair. Difference in interaction pattern was maximum in CSPF/CSPH owing to their high number of non-synonymous substitutions. Maximum difference in surface electrostatic potential was observed in case of CSPA, CSPG and CSPF. This research work emphasizes on discerning the molecular mechanism initiated by these proteins with a structural, mutational and functional approach.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03656-2.