Homology Modeling Research Articles

Rational In Silico Design of Selective TMPRSS6 Peptidomimetic Inhibitors via Exploitation of the S2 Subpocket.

TMPRSS6 is a potential therapeutic target for the treatment of iron overload due to its role in regulating levels of hepcidin. Although potent TMPRSS6 inhibitors have been previously developed, their lack of specificity requires optimization to avoid potential side effects before pursuing preclinical development with in vivo models. Here, using computer-aided drug design based on a TMPRSS6 homology model, we reveal that the S2 position of TMPRSS6 offers a potential avenue to achieve selectivity against other members of the TTSP family. Accordingly, we synthesized novel peptidomimetic molecules containing lipophilic amino acids at the P2 position to exploit this unexplored pocket. This enabled us to identify TMPRSS6-selective small molecules with low nanomolar affinity. Finally, pharmacokinetic parameters were determined, and a compound was found to be potent in cellulo toward its primary target while retaining TTSP-subtype selectivity and showing no signs of alteration in in vitro TEER experiments.

Computational 3D Models of Fe2+/3+-Aβ1-42 Complexes Associated with Alzheimer's Disease.

The interaction between iron and amyloid-beta (Aβ) peptides has received significant attention in Alzheimer's disease (AD) research due to its potential implications in developing this pathology. However, the coordination preferences of iron and Aβ1-42 have not been thoroughly investigated or remain unknown. This study employs a computational protocol that combines homology modeling techniques with quantum mechanics (DTF-xTB) calculations to build and evaluate several 3D models of Fe2+/3+-Aβ1-42. Our results reveal well-defined complexes for both the metal and peptide moieties, and we discuss the molecular interactions stabilizing these complexes by elucidating the coordinating environments and binding preferences. These proposed models offer valuable insights into the role of iron in Alzheimer's disease (AD) pathology.

Structural impacts of two disease-linked ADAR1 mutants: a molecular dynamics study.

Adenosine deaminases acting on RNA (ADARs) are pivotal RNA-editing enzymes responsible for converting adenosine to inosine within double-stranded RNA (dsRNA). Dysregulation of ADAR1 editing activity, often arising from genetic mutations, has been linked to elevated interferon levels and the onset of autoinflammatory diseases. However, understanding the molecular underpinnings of this dysregulation is impeded by the lack of an experimentally determined structure for the ADAR1 deaminase domain. In this computational study, we utilized homology modeling and the AlphaFold2 to construct structural models of the ADAR1 deaminase domain in wild-type and two pathogenic variants, R892H and Y1112F, to decipher the structural impact on the reduced deaminase activity. Our findings illuminate the critical role of structural complementarity between the ADAR1 deaminase domain and dsRNA in enzyme-substrate recognition. That is, the relative position of E1008 and K1120 must be maintained so that they can insert into the minor and major grooves of the substrate dsRNA, respectively, facilitating the flipping-out of adenosine to be accommodated within a cavity surrounding E912. Both amino acid replacements studied, R892H at the orthosteric site and Y1112F at the allosteric site, alter K1120 position and ultimately hinder substrate RNA binding.

Investigation of the influence of Candida tropicalis on bioethanol production using sugarcane bagasse: stochastic and in silico analysis.

This study investigated the impact of Candida tropicalis NITCSK13 on sugarcane bagasse (SCB) consolidated bioprocessing (CSB) using various parameters, such as pH, steam explosion (STEX) pretreatment, and temperature (at two different temperatures, cellulose hydrolysis and ethanol fermentation). The backpropagation neural network (BPNN) method simulated the optimal CSB conditions, achieving a maximum ethanol yield of 44 ± 0.32g/L (0.443g of ethanol/g of SCB) from STEX pretreated SCB within 48h at 55°C for cellulose hydrolysis and 33°C for ethanol fermentation and pH 3.5. The simulated conditions were experimentally validated and showed an R2 value of 0.998 and absolute average deviation (AAD) of 1.23%. The strain NITCSK13 also exhibited a high ethanol tolerance of 16% (v/v). The interactions between the inhibitors, cellobiose, furfural, and thermocellulase were assessed through molecular docking. The results revealed a maximum inhibitory constant of 3.7mM for furfural against the endoglucanase (EnG) of Humicola insolens (2ENG) at 50°C. Acremonium chrysogenum endoglucanase (5M2D) exhibited a maximum of 88.7µM for cellobiose at 50°C. The SWISS homology model of EnG from Candida viswanathii exhibited inhibitory effects similar to those of EnG from Thermoascus and Thermotoga, indicating that the moderately thermophilic yeast Candida sp. cellulase may be capable of efficiently tolerating inhibitors and could be a promising candidate for consolidated bioprocessing of cellulosic ethanol.

In silico homology modeling of dengue virus non-structural 4B (NS4B) protein and its molecular docking studies using triterpenoids

BackgroundDengue fever has become a significant worldwide health concern, because of its high morbidity rate and the potential for an increase in mortality rates due to lack of adequate treatment. There is an immediate need for the development of effective medication for dengue fever.MethodsHomology modeling of dengue virus (DENV) non-structural 4B (NS4B) protein was performed by SWISS-MODEL to predict the 3D structure of the protein. Structure validation was conducted using PROSA, PROCHECK, Ramachandran plot, and VERIFY-3D. MOE software was used to find out the in-Silico inhibitory potential of the five triterpenoids against the DENV-NS4B protein.ResultsThe SWISS-MODEL was employed to predict the three-dimensional protein structure of the NS4B protein. Through molecular docking, it was found that the chosen triterpenoid NS4B protein had a high binding affinity interaction. It was observed that the NS4B protein binding energy for 15-oxoursolic acid, betulinic acid, ursolic acid, lupeol, and 3-o-acetylursolic acid were − 7.18, − 7.02, − 5.71, − 6.67 and − 8.00 kcal/mol, respectively.ConclusionsNS4B protein could be a promising target which showed good interaction with tested triterpenoids which can be developed as a potential antiviral drug for controlling dengue virus pathogenesis by inhibiting viral replication. However, further investigations are necessary to validate and confirm their efficacy.

In silico evidence of bitopertin's broad interactions within the SLC6 transporter family.

The Glycine Transporter Type 1 (GlyT1) significantly impacts central nervous system functions, influencing glycinergic and glutamatergic neurotransmission. Bitopertin, the first GlyT1 inhibitor in clinical trials, was developed for schizophrenia treatment but showed limited efficacy. Despite this, bitopertin's repositioning could advance treating various pathologies. This study aims to understand bitopertin's mechanism of action using computational methods, exploring off-target effects, and providing a comprehensive pharmacological profile. Similarity Ensemble Approach (SEA) and SwissTargetPrediction initially predicted targets, followed by molecular modeling on SWISS-MODEL and GalaxyWeb servers. Binding sites were identified using PrankWeb, and molecular docking was performed with DockThor and GOLD software. Molecular dynamics analyses were conducted on the Visual Dynamics platform. Reverse screening on SEA and SwissTargetPrediction identified GlyT1 (SLC6A9), GlyT2 (SLC6A5), PROT (SLC6A7), and DAT (SLC6A3) as potential bitopertin targets. Homology modeling on SwissModel generated high-resolution models, optimized further on GalaxyWeb. PrankWeb identified similar binding sites in GlyT1, GlyT2, PROT, and DAT, indicating potential interaction. Docking studies suggested bitopertin's interaction with GlyT1 and proximity to GlyT2 and PROT. Molecular dynamics confirmed docking results, highlighting bitopertin's target stability beyond GlyT1. The study concludes that bitopertin potentially interacts with multiple SLC6 family targets, indicating a broader pharmacological property.

In silico evidence of bitopertin's broad interactions within the SLC6 transporter family.

The Glycine Transporter Type 1 (GlyT1) significantly impacts central nervous system functions, influencing glycinergic and glutamatergic neurotransmission. Bitopertin, the first GlyT1 inhibitor in clinical trials, was developed for schizophrenia treatment but showed limited efficacy. Despite this, bitopertin's repositioning could advance treating various pathologies. This study aims to understand bitopertin's mechanism of action using computational methods, exploring off-target effects, and providing a comprehensive pharmacological profile. Similarity Ensemble Approach (SEA) and SwissTargetPrediction initially predicted targets, followed by molecular modeling on SWISS-MODEL and GalaxyWeb servers. Binding sites were identified using PrankWeb, and molecular docking was performed with DockThor and GOLD software. Molecular dynamics analyses were conducted on the Visual Dynamics platform. Reverse screening on SEA and SwissTargetPrediction identified GlyT1 (SLC6A9), GlyT2 (SLC6A5), PROT (SLC6A7), and DAT (SLC6A3) as potential bitopertin targets. Homology modeling on SwissModel generated high-resolution models, optimized further on GalaxyWeb. PrankWeb identified similar binding sites in GlyT1, GlyT2, PROT, and DAT, indicating potential interaction. Docking studies suggested bitopertin's interaction with GlyT1 and proximity to GlyT2 and PROT. Molecular dynamics confirmed docking results, highlighting bitopertin's target stability beyond GlyT1. The study concludes that bitopertin potentially interacts with multiple SLC6 family targets, indicating a broader pharmacological property.

Biallelic variants in POPDC2 cause a novel autosomal recessive syndrome presenting with cardiac conduction defects and variable hypertrophic cardiomyopathy.

ABSTRACTPOPDC2encodes for the Popeye domain-containing protein 2 which has an important role in cardiac pacemaking and conduction, due in part to its cAMP-dependent binding and regulation of TREK-1 potassium channels. Loss ofPopdc2in mice results in sinus pauses and bradycardia and morpholino knockdown ofpopdc2 inzebrafish results in atrioventricular (AV) block. We identified bi-allelic variants inPOPDC2in 4 families that presented with a phenotypic spectrum consisting of sinus node dysfunction, AV conduction defects and hypertrophic cardiomyopathy. Using homology modelling we show that the identifiedPOPDC2variants are predicted to diminish the ability of POPDC2 to bind cAMP. Inin vitroelectrophysiological studies we demonstrated that, while co-expression of wild-type POPDC2 with TREK-1 increased TREK-1 current density, POPDC2 variants found in the patients failed to increase TREK-1 current density. While patient muscle biopsy did not show clear myopathic disease, it showed significant reduction of the expression of both POPDC1 and POPDC2, suggesting that stability and/or membrane trafficking of the POPDC1–POPDC2 complex is impaired by pathogenic variants in any of the two proteins. Single-cell RNA sequencing from human hearts demonstrated that co-expression of POPDC1 and 2 was most prevalent in AV node, AV node pacemaker and AV bundle cells. Sinoatrial node cells expressed POPDC2 abundantly, but expression of POPDC1 was sparse. Together, these results concur with predisposition to AV node disease in humans with loss-of-function variants inPOPDC1and POPDC2 and presence of sinus node disease inPOPDC2,but not inPOPDC1related disease in human. Using population-level genetic data of more than 1 million individuals we showed that none of the familial variants were associated with clinical outcomes in heterozygous state, suggesting that heterozygous family members are unlikely to develop clinical manifestations and therefore might not necessitate clinical follow-up. Our findings provide evidence forPOPDC2as the cause of a novel Mendelian autosomal recessive cardiac syndrome, consistent with previous work showing that mice and zebrafish deficient in functionalPOPDC2display sinus and AV node dysfunction.GRAPHICAL ABSTRACT

Structural Studies of the Taurine Transporter: A Potential Biological Target from the GABA Transporter Subfamily in Cancer Therapy.

The taurine transporter (TauT, SLC6A6) is a member of the solute carrier 6 (SLC6) family, which plays multiple physiological roles. The SLC6 family is divided into four subfamilies: GABA (γ-aminobutyric acid), monoamine, glycine and neutral amino acid transporters. Proteins from the GABA group, including the taurine transporter, are primarily considered therapeutic targets for treating central nervous system disorders. However, recent studies have suggested that inhibitors of SLC6A6 could also serve as anticancer agents. Overexpression of TauT has been associated with the progression of colon and gastric cancer. The pool of known ligands of this transporter is limited and the exact spatial structure of taurine transporter remains unsolved. Understanding its structure could aid in the development of novel inhibitors. Therefore, we utilized homology modelling techniques to create models of TauT. Docking studies and molecular dynamics simulations were conducted to describe protein-ligand interactions. We compared the obtained information for TauT with literature data on other members of the GABA transporter group. Our in silico analysis allowed us to characterize the transporter structure and point out amino acids crucial for ligand binding: Glu406, Gly62 and Tyr138. The significance of selected residues was confirmed through structural studies of mutants. These results will aid in the development of novel taurine transporter inhibitors, which can be explored as anticancer agents.

Genome-wide identification and expression analysis of calmodulin and calmodulin-like genes in passion fruit (Passiflora edulis) and their involvement in flower and fruit development

BackgroundThe calmodulin (CaM) and calmodulin-like (CML) proteins play regulatory roles in plant growth and development, responses to biotic and abiotic stresses, and other biological processes. As a popular fruit and ornamental crop, it is important to explore the regulatory mechanism of flower and fruit development of passion fruit.ResultsIn this study, 32 PeCaM/PeCML genes were identified from passion fruit genome and were divided into 9 groups based on phylogenetic analysis. The structural analysis, including conserved motifs, gene structure and homologous modeling, illustrates that the PeCaM/PeCML in the same subgroup have relative conserved structural features. Collinearity analysis suggested that the expansion of the CaM/CML gene family likely took place mainly by segmental duplication, and the whole genome replication events were closely related with the rapid expansion of the gene group. PeCaM/PeCMLs were potentially required for different floral tissues development. Significantly, PeCML26 had extremely high expression levels during ovule and fruit development compared with other PeCML genes, suggesting that PeCML26 had potential functions involved in the development of passion fruit flowers and fruits. The co-presence of various cis-elements associated with growth and development, hormone responsiveness, and stress responsiveness in the promoter regions of these PeCaM/PeCMLs might contribute to their diverse regulatory roles. Furthermore, PeCaM/PeCMLs were also induced by various abiotic stresses. This work provides a comprehensive understanding of the CaM/CML gene family and valuable clues for future studies on the function and evolution of CaM/CML genes in passion fruit.ConclusionA total of 32 PeCaM/PeCML genes were divided into 9 groups. The PeCaM/PeCML genes showed differential expression patterns in floral tissues at different development stages. It is worth noting that PeCML26, which is highly homologous to AtCaM2, not only interacts with multiple BBR-BPC TFs, but also has high expression levels during ovule and fruit development, suggesting that PeCML26 had potential functions involved in the development of passion fruit flowers and fruits. This research lays the foundation for future investigations and validation of the potential function of PeCaM/PeCML genes in the growth and development of passion fruit.

Computational bioprospecting of phytoconstituents as potential inhibitors for peptide deformylase from Streptococcus oralis: An opportunistic pathogen

Streptococcus oralis an opportunistic bacterium has been reported to be involved in various blood borne infections like subacute bacterial endocarditis, septicemia, bacterial meningitis and in some cases dental caries too. Among various targets the peptide deformylase, of S.oralis appears to be most potent druggable target as it is involved in protein synthesis is opted for the current study. Due to unavailability of PDB structure of peptide deformylase from S. oralis the study initiates with homology modelling of the protein and 6OW2 of S pneumoniae is considered as the template. Thereafter, Molecular docking, Molecular dynamic simulation, ADME analysis, and MMPBSA analysis was carried out to explore the inhibitory potential of phyto-constituents as potential inhibitors for Peptide deformylase from S.oralis. Actinonin was considered as reference drug. Among 2370 phyto compounds the best observations were recorded for A1-Barrigenol (IMPHY010984) with binding affinity of −8.5 kcal/mol. Calculated RMSD, RMSF, Binding Free Energy for IMPHY010984 averaged at about 0.10 ± 0.03 nm, 0.08 ± 0.05 nm, 131 ± 21 kJ/mol respectively whereas the RMSD, RMSF, Binding Free Energy recorded for reference drug averaged at about 0.19 ± 0.04 nm, 0.11 ± 0.08 nm, −94 ± 18 kJ/mol respectively. Based on in silico observations IMPHY010984 is proved out as superior candidate over reference drug. The study reflects the potential of IMPHY010984 as prophylactic therapeutics for S.oralis.

Phenylalanine Butyramide: A Butyrate Derivative as a Novel Inhibitor of Tyrosinase.

Metabolites resulting from the bacterial fermentation of dietary fibers, such as short-chain fatty acids, especially butyrate, play important roles in maintaining gut health and regulating various biological effects in the skin. However, butyrate is underutilized due to its unpleasant odor. To circumvent this organoleptic unfavorable property, phenylalanine butyramide (PBA), a butyrate precursor, has been synthesized and is currently available on the market. We evaluated the inhibition of mushroom tyrosinase by butyrate and PBA through in vitro assays, finding IC50 values of 34.7 mM and 120.3 mM, respectively. Docking calculations using a homology model of human tyrosinase identified a putative binding mode of PBA into the catalytic site. The anti-aging and anti-spot efficacy of topical PBA was evaluated in a randomized, double-blind, parallel-arm, placebo-controlled clinical trial involving 43 women affected by photo-damage. The results of this study showed that PBA significantly improved skin conditions compared to the placebo and was well tolerated. Specifically, PBA demonstrated strong skin depigmenting activity on both UV and brown spots (UV: -12.7% and -9.9%, Bs: -20.8% and -17.7% after 15 and 30 days, respectively, p < 0.001). Moreover, PBA brightened and lightened the skin (ITA°: +12% and 13% after 15 and 30 days, respectively, p < 0.001). Finally, PBA significantly improved skin elasticity (Ua/Uf: +12.4% and +32.3% after 15 and 30 days, respectively, p < 0.001) and firmness (Uf: -3.2% and -14.9% after 15 and 30 days, respectively, p < 0.01).

Involvement of the cytochrome P450 genes CYP6DW3 and CYP4C64 in afidopyropen resistance in Bemisiatabaci Mediterranean (Q Biotype)

The tobacco whitefly, Bemisia tabaci, is a notorious pest affecting various crops globally, and it exhibits high levels of resistance to various insecticides. Afidopyropen is a recently commercialized pyropene insecticide for B. tabaci control with high selectivity and a novel mode of action. We previously identified a high level of afidopyropen resistance in a field-collected population after selection in the lab, and named it the HD-Afi strain. In the present study, minimal cross-resistance in the HD-Afi strain was found between afidopyropen and other common chemical agents. However, the P450 enzyme activity in HD-Afi was 2.18 times the level in susceptible strain HD-S. Expression analysis revealed that two of 12 candidate P450 genes, namely CYP6DW3 and CYP4C64, were significantly up-regulated in HD-Afi. Silencing CYP6DW3 and CYP4C64 by RNA interference (RNAi) substantially increased the susceptibility of whitefly adults, confirming their involvement in afidopyropen resistance. Homology modeling and molecular docking analyses demonstrated stable binding of afidopyropen to CYP6DW3 and CYP4C64, with binding free energies of -6.87 and -6.11 kcal mol-1, respectively. The findings of this study suggest that the induction of CYP6DW3 and CYP4C64 facilitates afidopyropen detoxification, contributing to the development of resistance in B. tabaci.

Biocompatible ionic liquid self-assembled transdermal delivery system (P-M IL/AKE) for significantly improved skin photoaging and reduced wrinkles

Most peptide drugs have good biological activity and safety, but their stability and permeability are poor, limiting their applications in the field of local drug delivery. Snake venom peptide (AKE) is a tripeptide with good anti-wrinkle effects, but its weak esterophilicity, poor penetration, and low bioavailability limit its full anti-wrinkle potential. To address this issue, a bioactive ionic liquid (P-M IL) was prepared and an efficient self-assembled peptide drug delivery system was constructed through AKE loading, which increased the cumulative absorption of AKE in the skin by more than five-fold, resulting in good biocompatibility. Further, P-M IL/AKE can significantly increase collagen protein secretion by human fibroblasts, improved skin photoaging and reduce the areas and volume of wrinkle in the eye tail lines. In addition, for the first time, quantum chemical calculations, homology modeling, molecular docking and binding free energy calculations were used to analyze the structural stability of the system, construct the 3D model, explore the interaction mode and reveal the theoretical basis of anti-wrinkle effect. Therefore, this supramolecular self-assembled transdermal drug delivery system has potential advantages, providing a novel approach for the highly efficient and non-invasive delivery of peptide drugs.

Ammonia monooxygenase-mediated cometabolic biotransformation of volatile 4-chlorophenol in nitrifying counter-diffused biofilms: A combined molecular dynamics simulation, DFT calculation and experimental study

Ammonia monooxygenase (AMO)-mediated cometabolism of organic pollutants has been widely observed in biological nitrogen removal process. However, its molecular mechanism remains unclear, hindering its practical application. Furthermore, conventional nitrification systems encounter significant challenges such as air pollution and the loss of ammonia-oxidizing bacteria, when dealing with wastewater containing volatile organic pollutants. This study developed a nitrifying membrane-aerated biofilm reactor (MABR) to enhance the biodegradation of volatile 4-chlorophenol (4-CP). Results showed that 4-CP was primarily removed via Nitrosomonas nitrosa-mediated cometabolism in the presence of NH4+-N, supported by the increased nicotinamide adenine dinucleotide (NADH) and adenosine triphosphate (ATP) content, AMO activity and the related genes abundance. Hydroquinone, detected for the first time and produced via oxidative dechlorination, as well as 4-chlorocatechol was primary transformation products of 4-CP. Nitrosomonas nitrosa AMO structural model was constructed for the first time using homology modeling. Molecular dynamics simulation suggested that the ortho-carbon in the benzene ring of 4-CP was more prone to metabolismcompared to the ipso-carbon. Density functional theory calculation revealed that 4-CP was metabolized by AMO via H-abstraction-OH-rebound reaction, with a significantly higher rebound barrier at the ipso-carbon (16.37 kcal·mol-1) as compared to the ortho-carbon (6.7 kcal·mol-1). This study fills the knowledge gap on the molecular mechanism of AMO-mediated cometabolism of organic pollutants, providing practical and theoretical foundations for improving volatile organic pollutants removal through nitrifying MABR.

Simulation screening of umami peptides in the microbiota of soy sauce based on molecular docking and molecular dynamics

Simulation screening and molecular simulation techniques are gradually being applied to the study of bioactive peptides. In this study, we explored and simulated the screening process using metagenomic sequencing data of microorganisms from soy sauce. An LSTM-based screening model was constructed to identify the top 10 scoring peptides from a sample dataset of 19,335 raw peptides, primarily derived from Bacillus and Weissella confusa. Through homology modeling, molecular docking, and molecular dynamics simulations, Ser, Asp, Asn, Glu, and Gln were identified as the main amino acids involved in ligand-receptor binding. Hydrogen bonding, van der Waals forces, and hydrophobic interactions were revealed as the primary forces mediating receptor and ligand binding. The umami thresholds of ISWCFTY and QISPYRRI were verified through wet experiments to be 0.13 mg/mL and 0.11 mg/mL, respectively, exhibiting high umami intensity. Overall, this study presents a high-throughput screening method for predicting umami peptides, offering cost reduction, simplified steps, and saved manpower.

Potential Antidiabetic Activity of β-sitosterol from Zingiber roseum Rosc. via Modulation of Peroxisome Proliferator-activated Receptor Gamma (PPARγ).

To evaluate the antidiabetic potential of β-sitosterol from Zingiber roseum. Diabetes mellitus is a cluster of metabolic disorders, and 90% of diabetic patients are affected with Type II diabetes (DM2). For the treatment of DM2, thiazolidinedione drugs (TZDs) were proposed, but recent studies have shown that TZDs have several detrimental effects, such as weight gain, kidney enlargement (hypertrophy), fluid retention, increased risk of bone fractures, and potential harm to the liver (hepatotoxicity). That is why a new molecule is needed to treat DM2. The current research aimed to assess the efficacy of β-Sitosterol from methanolic extract of Zingiber roseum in managing diabetes via PPARγ modulation. Zingiber roseum was extracted using methanol, and GC-MS was employed to analyze the extract. Through homology modeling, PPARγ structure was predicted. Molecular docking, MD simulation, free binding energies, QSAR, ADMET, and bioactivity and toxicity scores were all used during the in-depth computer-based research. Clinically, agonists of synthetic thiazolidinedione (TZDs) have been used therapeutically to treat DM2, but these TZDs are associated with significant risks. Hence, GC-MS identified phytochemicals to search for a new PPAR-γ agonist. Based on the in-silico investigation, β-sitosterol was found to have a higher binding affinity (-8.9 kcal/mol) than standard drugs. MD simulations and MMGBSA analysis also demonstrated that β-sitosterol bound to the PPAR-γ active site stably. It can be concluded that β-sitosterol from Z. roseum attenuates Type-II diabetes by modulating PPARγ activity.

Optimizing Bioethanol Production from Hassawi Rice Straw with Aspergillus sp. NAS51 cellulosic Enzyme and In Silico homology Modeling

Researchers' attention has been greatly focused on the use of lignocellulosic wastes to produce bioethanol and biogas due to the depletion of non-renewable energy resources. In our effort to find a potent cellulase-producing fungal strain, the fungus NAS51 was isolated among eight isolates from a sponge collected from the Red Sea, Jeddah, and selected as it displayed the highest cellulase activity including (FP-ase, CMC-ase, and β-glucosidase enzymes) at levels of 3.13 U/ml, 2.52 U/ml, and 0.69 U/ml, respectively. The fungus was identified morphologically and genetically by sequencing its 18SrRNA gene as Aspergillus sp. NAS51. The cellulase activity of Aspergillus sp. NAS51 was optimized and maximum enzyme production was obtained at initial pH7, temp 30oC, incubation period 11 days, moisture content 70%, urea as a nitrogen source, and K2HPO4 (2g/L). The crude cellulases from strain NAS51 were characterized, and results showed that the enzyme is stable over a wide pH range (6-10), with peak activity at pH 7.0 and 40°C. The cellulase gene has been sequenced and the protein 3D structure was generated via in silico homology modeling. Determination of binding sites and biological annotations of the constructed protein was carried out via COACH and COFACTOR based on the I-TASSER structure prediction. To reach the maximum enzyme hydrolysis, the rice straw collected from Al-Ahsa, Kingdom of Saudi Arabia was pretreated with NaOH 1.5% to remove lignin and to enhance the saccharification process by Cellulase. The saccharified product was measured using HPLC, fermented by S. cerevisiae and the bioethanol yield produced from the fermentation was 0.454 ml ethanol/g fermentable sugars. This study indicates the potential application of fungal enzymes such as cellulase enzymes in biofuel generation and waste management.

AutoDock and molecular dynamics-based therapeutic potential prediction of flavonoids for primary Sjögren's syndrome

Primary Sjögren's Syndrome (pSS) is a systemic autoimmune disease that leads to reduced saliva production, primarily affecting women due to estrogen deficiency. The estrogen receptor α (ERα) plays a crucial role in mediating the expression of the aquaporin 5 (AQP5) gene through the estrogen response element-dependent signaling pathway, making ERα a key drug target for pSS. Several flavonoids have been reported to have the potential to treat pSS. This study aimed to screen and compare flavonoids binding to ERα using AutoDock, providing a basis for treating pSS with flavonoids. The estrogenic potential of six representative flavonoids was examined in this study. Molecular docking revealed that the binding energy of all six flavonoids to ERα was less than −5.6 kcal/mol. Apigenin, naringenin, and daidzein were the top three flavonoids with even lower binding energies of −7.8, −8.09, and −8.59 kcal/mol, respectively. Similar to the positive control estradiol, apigenin, naringenin, and daidzein showed hydrogen bond interactions with GLU353, GLY521, and HIS524 at the active site. The results of luciferase reporter assays demonstrated that apigenin, naringenin, and daidzein significantly enhanced the transcription of estrogen receptor element (ERE) in the PGL3/AQP5 promoter. Furthermore, molecular dynamics simulations using GROMACS for a time scale of 100 ns revealed relatively stable binding of apigenin-ERα, naringenin-ERα, and daidzein-ERα. Mechanistically, homology modeling indicated that GLU353, GLY521, and HIS524 were the key residues of ERα exerting an estrogenic effect. The therapeutic effect of apigenin on dry mouth in pSS models was further validated. In conclusion, these results indicate the estrogenic and pSS therapeutic potential of apigenin, naringenin, and daidzein.

A Comparison of Antibody-Antigen Complex Sequence-to-Structure Prediction Methods and their Systematic Biases.

The ability to accurately predict antibody-antigen complex structures from their sequences could greatly advance our understanding of the immune system and would aid in the development of novel antibody therapeutics. There have been considerable recent advancements in predicting protein-protein interactions (PPIs) fueled by progress in machine learning (ML). To understand the current state of the field, we compare six representative methods for predicting antibody-antigen complexes from sequence, including two deep learning approaches trained to predict PPIs in general (AlphaFold-Multimer, RoseTTAFold), two composite methods that initially predict antibody and antigen structures separately and dock them (using antibody-mode ClusPro), local refinement in Rosetta (SnugDock) of globally docked poses from ClusPro, and a pipeline combining homology modeling with rigid-body docking informed by ML-based epitope and paratope prediction (AbAdapt). We find that AlphaFold-Multimer outperformed other methods, although the absolute performance leaves considerable room for improvement. AlphaFold-Multimer models of lower-quality display significant structural biases at the level of tertiary motifs (TERMs) towards having fewer structural matches in non-antibody containing structures from the Protein Data Bank (PDB). Specifically, better models exhibit more common PDB-like TERMs at the antibody-antigen interface than worse ones. Importantly, the clear relationship between performance and the commonness of interfacial TERMs suggests that scarcity of interfacial geometry data in the structural database may currently limit application of machine learning to the prediction of antibody-antigen interactions.