Using methylation data to improve transcription factor binding prediction

Modelling the regulatory mechanisms that determine cell fate, response to external perturbation, and disease state depends on measuring many factors, a task made more difficult by the plasticity of the epigenome. Scanning the genome for the sequence patterns defined by Position Weight Matrices (PWM) can be used to estimate transcription factor (TF) binding locations. However, this approach does not incorporate information regarding the epigenetic context necessary for TF binding. CpG methylation is an epigenetic mark influenced by environmental factors that is commonly assayed in human cohort studies. We developed a framework to score inferred TF binding locations using methylation data. We intersected motif locations identified using PWMs with methylation information captured in both whole-genome bisulfite sequencing and Illumina EPIC array data for six cell lines, scored motif locations based on these data, and compared with experimental data characterizing TF binding (ChIP-seq). We found that for most TFs, binding prediction improves using methylation-based scoring compared to standard PWM-scores. We also illustrate that our approach can be generalized to infer TF binding when methylation information is only proximally available, i.e. measured for nearby CpGs that do not directly overlap with a motif location. Overall, our approach provides a framework for inferring context-specific TF binding using methylation data. Importantly, the availability of DNA methylation data in existing patient populations provides an opportunity to use our approach to understand the impact of methylation on gene regulatory processes in the context of human disease.

Design, synthesis, and bioevaluation of 1h-pyrrolo[3,2-c]pyridine derivatives as colchicine-binding site inhibitors with potent anticancer activities

A new series of 1H-pyrrolo[3,2-c]pyridine derivatives were designed and synthesised as colchicine-binding site inhibitors. Preliminary biological evaluations showed that most of the target compounds displayed moderate to excellent antitumor activities against three cancer cell lines (HeLa, SGC-7901, and MCF-7) in vitro. Among them, 10t exhibited the most potent activities against three cancer cell lines with IC50 values ranging from 0.12 to 0.21 μM. Tubulin polymerisation experiments indicated that 10t potently inhibited tubulin polymerisation at concentrations of 3 μM and 5 μM, and immunostaining assays revealed that 10t remarkably disrupted tubulin microtubule dynamics at a concentration of 0.12 μM. Furthermore, cell cycle studies and cell apoptosis analyses demonstrated that 10t at concentrations of 0.12 μM, 0.24 μM, and 0.36 μM significantly caused G2/M phase cell cycle arrest and apoptosis. The results of molecular modelling studies suggested that 10t interacts with tubulin by forming hydrogen bonds with colchicine sites Thrα179 and Asnβ349. In addition, the prediction of physicochemical properties disclosed that 10t conformed well to the Lipinski's rule of five.

Synergistic effects of putative Ca2+-binding sites of calmodulin in fungal development, temperature stress and virulence of Aspergillus fumigatus

In pathogenic fungi, calcium-calmodulin-dependent serine-threonine-specific phosphatase calcineurin is involved in morphogenesis and virulence. Therefore, calcineurin and its tightly related protein complexes are attractive antifungal drug targets. However, there is limited knowledge available on the relationship between in vivo Ca2+-binding sites of calmodulin (CaM) and its functions in regulating stress responses, morphogenesis, and pathogenesis. In the current study, we demonstrated that calmodulin is required for hyphal growth, conidiation, and virulence in the human fungal pathogen, Aspergillus fumigatus. Site-directed mutations of calmodulin revealed that a single Ca2+-binding site mutation had no significant effect on A. fumigatus hyphal development, but multiple Ca2+-binding site mutations exhibited synergistic effects, especially when cultured at 42 °C, indicating that calmodulin function in response to temperature stress depends on its Ca2+-binding sites. Western blotting implied that mutations in Ca2+-binding sites caused highly degraded calmodulin fragments, suggesting that the loss of Ca2+-binding sites results in reduced protein stability. Moreover, normal intracellular calcium homeostasis and the nuclear translocation of the transcriptional factor CrzA are dependent on Ca2+-binding sites of AfCaM, demonstrating that Ca2+-binding sites of calmodulin are required for calcium signalling and its major transcription factor CrzA. Importantly, in situ mutations for four Ca2+-binding sites of calmodulin resulted in an almost complete loss of virulence in the Galleria mellonella wax moth model. This study shed more light on the functional characterization of putative calcium-binding sites of calmodulin in the morphogenesis and virulence of A. fumigatus, which enhances our understanding of calmodulin biological functions in cells of opportunistic fungal pathogens.

Recent advances on the structure and the function relationships of the TRPV4 ion channel

The members of the superfamily of Transient Receptor Potential (TRP) ion channels are physiologically important molecules that have been studied for many years and are still being intensively researched. Among the vanilloid TRP subfamily, the TRPV4 ion channel is an interesting protein due to its involvement in several essential physiological processes and in the development of various diseases. As in other proteins, changes in its function that lead to the development of pathological states, have been closely associated with modification of its regulation by different molecules, but also by the appearance of mutations which affect the structure and gating of the channel. In the last few years, some structures for the TRPV4 channel have been solved. Due to the importance of this protein in physiology, here we discuss the recent progress in determining the structure of the TRPV4 channel, which has been achieved in three species of animals (Xenopus tropicalis, Mus musculus, and Homo sapiens), highlighting conserved features as well as key differences among them and emphasizing the binding sites for some ligands that play crucial roles in its regulation.

In vitro investigation of the binding characteristics of dacomitinib to human α 1-acid glycoprotein: Multispectral and computational modeling

Dacomitinib is a highly selective second-generation tyrosine kinase inhibitor that can irreversibly bind to tyrosine kinase and is mainly used in the treatment of lung cancer. The binding characteristics of dacomitinib with human α 1-acid glycoprotein (HAG) were analyzed by multispectral and computational simulation techniques. The fluorescence spectra showed that dacomitinib can quench the fluorescence of HAG by forming the HAG-dacomitinib complex with a molar ratio of 1:1 (static quenching). At the temperature similar to that of the human body, the affinity of dacomitinib to HAG (8.95 × 106 M-1) was much greater than that to BSA (3.39 × 104 M-1), indicating that dacomitinib will give priority to binding onto HAG. Thermodynamics parameters analysis and driving force competition experiments showed that hydrogen bonding and hydrophobic forces were the major sources for keeping the complex of HAG-dacomitinib stable. The experimental outcomes also showed that the binding of dacomitinib can lead to the loosening of the skeleton structure of HAG, which led to a slight change in the secondary structure, and also reduces the hydrophobicity of the microenvironment of Trp and Tyr residues. The binding sites of dacomitinib on HAG and the contribution of key amino acid residues to the binding reaction were determined by molecular docking and molecular dynamics (MD) simulation. In addition, it was found that there was a synergistic effect between dacomitinib and Mg2+ and Co2+ ions. Mg2+ and Co2+ could increase the Kb of dacomitinib to HAG and prolong the half-life of dacomitinib.

Structural and functional analysis of Cyanovirin-N homologs: Carbohydrate binding affinities and antiviral potential of cyanobacterial peptides

Cyanobacteria, a group of photosynthetic prokaryotes, can sinthesize several substances due to their secondary metabolism, with notable properties, such as Cyanovirin-N(CVN), a carbohydrate-binding lectin, that exhibits antiviral activity against several pathogens, due to its ability to bind viral surface carbohydrates such as mannose, thus interfering with the viral entry on the cell. CVN has been described in several cyanobacterial strains and shows biotechnological potential for the development of drugs of pharmaceutical interest. This study focuses on the genomic exploration and characterization of Cyanovirin-N homologs to assess the conservation of carbohydrate-binding affinity within the group. The analysis of their antiviral properties was carried out using bioinformatics tools to study protein models through an in silico pipeline, following the steps of genomic prospection on public databases, homology modeling, docking, molecular dynamics and energetic analysis. Mannose served as the reference ligand, and the lectins' binding affinity with mannose was assessed across Cyanovirin-N homologs. Genomic mining identified 33 cyanobacterial lectin sequences, which underwent structural and functional characterization. The results obtained from this work indicate strong carbohydrate affinity on several homologs, pointing to the conservation of antiviral properties alongside the group. However, this affinity was not uniformly distributed among sequences, exhibiting significant heterogeneity in binding site residues, suggesting potential multi-ligand binding capabilities on the Cyanovirin-N homologs group. Studies focused on the properties involved in these molecules and the investigation of the genetic diversity of Cyanovirin-N homologs could provide valuable insights into the discovery of new drug candidates, harvesting the potential of bioinformatics for large-scale functional and structural analysis.

Uncovering PPAR-γ agonists: An integrated computational approach driven by machine learning

Peroxisome proliferator-activated receptor gamma (PPAR-γ) serves as a nuclear receptor with a pivotal function in governing diverse facets of metabolic processes. In diabetes, the prime physiological role of PPAR-γ is to enhance insulin sensitivity and regulate glucose metabolism. Although PPAR-γ agonists such as Thiazolidinediones are effective in addressing diabetes complications, it is vital to be mindful that they are associated with substantial side effects that could potentially give rise to health challenges. The recent surge in the discovery of selective modulators of PPAR-γ inspired us to formulate an integrated computational strategy by leveraging the promising capabilities of both machine learning and in silico drug design approaches. In pursuit of our objectives, the initial stage of our work involved constructing an advanced machine learning classification model, which was trained utilizing chemical information and physicochemical descriptors obtained from known PPAR-γ modulators. The subsequent application of machine learning-based virtual screening, using a library of 31,750 compounds, allowed us to identify 68 compounds having suitable characteristics for further investigation. A total of four compounds were identified and the most favorable configurations were complemented with docking scores ranging from -8.0 to -9.1 kcal/mol. Additionally, the compounds engaged in hydrogen bond interactions with essential conserved residues including His323, Leu330, Phe363, His449 and Tyr473 that describe the ligand binding site. The stability indices investigated herein for instance root-mean-square fluctuations in the backbone atoms indicated higher mobility in the region of orthosteric site in the presence of agonist with the deviation peaks in the range of 0.07-0.69 nm, signifying moderate conformational changes. The deviations at global level revealed that the average values lie in the range of 0.25-0.32 nm. In conclusion, our identified hits particularly, CHEMBL-3185642 and CHEMBL-3554847 presented outstanding results and highlighted the stable conformation within the orthosteric site of PPAR-γ to positively modulate the activity.

Comparison of binding sites and affinity of flavonol-Cu(II) complexes with the same parent nucleus: Synthesis, DFT prediction, and coordination pattern

This study explores the coordination dynamics between dietary polyphenols, specifically kaempferol, quercetin, and myricetin, and Cu ions in aqueous environments. A novel synthesis method for flavonol-Cu(II) coordination compounds is introduced, effectively reducing interference from free metal ions. Our results reveal consistent binding patterns of Cu ions with flavonols (2:1 ratio of flavonol to Cu(II)), predominantly at the 4,5 sites. Various analytical techniques are used to validate these coordination ratios and sites. The binding affinity of the flavonols for Cu ions follows a descending sequence: myricetin > quercetin > kaempferol. Notably, coordination with Cu ions enhances the free-radical scavenging activities of these flavonols. These findings hold substantial importance for food chemistry, biology, and medicine, providing crucial insights into the way dietary flavonols form stable structures in environments similar to human body fluids and their interactions with metal ions, opening new possibilities for their application and understanding in diverse scientific domains.

Characterization of ultrasound-assisted covalent binding interaction between β-lactoglobulin and dicaffeoylquinic acid: Great potential for the curcumin delivery

The low bioavailability and poor gastrointestinal instability of curcumin hampers its application in pharmaceutical and food industries. Thus, it is essential to explore efficient carrier (e.g. a combination of polyphenols and proteins) for food systems. In this study, covalent β-lactoglobulin (LG)-dicaffeoylquinic acids (DCQAs) complexes were prepared by combining ultrasound and free radical induction methods. Covalent interactions between LG and DCQAs were confirmed by analyzing reactive groups. Variations in secondary or tertiary structure and potential binding sites of covalent complexes were explored using Fourier transform infrared spectroscopy and circular dichroism. Results showed that the β-sheet content decreased and the unordered content increased significantly (P < 0.05). The embedding rate of curcumin in prepared LG-DCQAs complexes using ultrasound could reach 49 % - 62 %, proving that complexes could embed curcumin effectively. This study highlights the benefit of ultrasound application in fabrication of protein-polyphenol complexes for delivering curcumin.

Targeting ATP-binding site of WRN Helicase: Identification of novel inhibitors through pocket analysis and Molecular Dynamics-Enhanced virtual screening

WRN helicase is a critical protein involved in maintaining genomic stability, utilizing ATP hydrolysis to dissolve DNA secondary structures. It has been identified as a promising synthetic lethal target for microsatellite instable (MSI) cancers. However, few WRN helicase inhibitors have been discovered, and their potential binding sites remain unexplored. In this study, we analyzed potential binding sites for WRN inhibitors and focused on the ATP-binding site for screening new inhibitors. Through molecular dynamics-enhanced virtual screening, we identified two compounds, h6 and h15, which effectively inhibited WRN's helicase and ATPase activity in vitro. Importantly, these compounds selectively targeted WRN's ATPase activity, setting them apart from other non-homologous proteins with ATPase activity. In comparison to the homologous protein BLM, h6 exhibits some degree of selectivity towards WRN. We also investigated the binding mode of these compounds to WRN's ATP-binding sites. These findings offer a promising strategy for discovering new WRN inhibitors and present two novel scaffolds, which might be potential for the development of MSI cancer treatment.

Different functional groups of carbon dots influence the formation of protein crowns and pepsin characteristic in vitro digestion

Application of nanomaterials (NMs) in agriculture poses an ingestion risk to humans and may affect the digestive process. Different fates of NMs with differential charges in the gastrointestinal tract should be considered. In this study, the interaction between three carbon dots (CDs) carried with different functional groups (-NH2, -OH, and -COOH) and pepsin was analyzed through an in vitro digestion model. The results showed that CDs significantly reduced pepsin activity. Among them, CDs-NH2 had the greatest effect, following by CDs-OH, and CDs-COOH. Besides, molecular docking demonstrated the specific binding site of CDs to pepsin, while the most stable binding energy (-8.10 kcal/mol) was formed between CDs-NH2 and pepsin. Further, CDs formed a nanomaterial-protein crown structure with pepsin. The present study enriches the functional group properties of CDs in the digestion and provides new ideas for the potential human health of NMs.

Investigating the binding interaction of quinoline yellow with bovine serum albumin and anti-amyloidogenic behavior of ferulic acid on QY-induced BSA fibrils

Protein aggregation induces profound changes in the structure along with the conformation of the protein, and is responsible for the pathogenesis of a number of neurodegenerative conditions such as Huntington's, Creutzfeldt-Jacob, Type II diabetes mellitus, Alzheimer's, etc. Numerous multi-spectroscopic approaches and in-silico experiments were utilized to investigate BSA's biomolecular interaction and aggregation in the presence of quinoline yellow. The present research investigation evaluated the interaction of BSA with the food colorant (QY) at two different pH (7.4 and 2.0). The development of the BSA-QY complex was established with UV visible and fluorescence spectroscopy. The quenching of fluorescence upon the interaction of BSA with QY revealed the static nature of quenching mechanism. The Kb value obtained from our result is 4. 54 × 10-4 M-1. The results from the competitive site marker study infer that quinoline yellow is binding with the sub-domain IB of bovine serum albumin, specifically on site III. Three-dimensional fluorescence and synchronous fluorescence spectroscopy were applied for monitoring the alterations in the microenvironment of BSA upon the addition of quinoline yellow. The results from turbidity and RLS studies showed that higher concentrations of QY (80-400 µM) triggered bovine serum albumin (BSA) aggregation at pH 2.0. At pH 7.4, QY couldn't manage to trigger bovine serum albumin aggregation, perhaps because of the repulsion between negatively charged dye (QY) and anionic bovine serum albumin. The results from far-UV CD, Congo Red, and scanning electron microscopy implicate that the QY-induced aggregates exhibit amyloid fibril-like structures. Molecular docking results revealed that hydrophobic interactions, hydrogen bonding, and Pi-Sulfur interactions contribute to QY-induced aggregation of BSA. Further, the amyloid inhibitory potential of ferulic acid (FA), a phenolic acid on QY-induced aggregation of BSA, has also been assessed. The QY-induced amyloid fibrils are FA-soluble, as confirmed by turbidity, RLS, and far-UV CD studies. Far-UV CD results showed that FA retains α helix and inhibits cross β sheet formation when the BSA samples were pre-incubated with increasing concentrations of FA (0-500 µM). Our findings conclude that QY dye successfully stimulates BSA aggregation, but ferulic acid inhibits QY-induced aggregation of BSA. Thus, FA can serve as a therapeutic agent and can help in the treatment of various amyloid-related conditions.

Thermodynamic and functional changes of alpha-chymotrypsin after interaction with gallic acid

In the present study, the interaction mechanism between gallic acid (GA) and α-Chymotrypsin (α-CT) was investigated by employing a series ofspectroscopic methods, computational docking and molecular dynamic (MD) simulation. Fluorescence spectra analysis indicated the formation of a stable complex between GA and α-CT, where the quenching of the fluorescence emission was predominantly characterized by a static mechanism. TheCA obtained binding constants for the α-CT-GA complex were in the order of 103 M-1, indicating the moderate binding affinity of GA for α-CT. The corresponding CD findings showed that the interaction between GA and α-CT resulted in an alteration of the protein's secondary structure. The findings of the enzyme activity investigation clearly showed that the presence of GA led to a notable decline in the enzymatic activity of α-CT, highlighting GA's function as an effective inhibitor for α-CT. The molecular docking simulations revealed the optimal binding site for the GA molecule within the α-CT structure and MD simulations confirmed the stability of the α-CT-GA complex. This research expands our comprehension regarding the behavior of enzymes in the presence of small-molecule ligands and opens avenues for food safety.

Modulation of protein-ligand interactions in the presence of ZIF-8: Spectroscopy and molecular dynamics simulation

In this paper, we investigated the protein-ligand interactions in the presence of ZIF-8 using multi-spectroscopic approaches and molecular dynamics simulation. Fluorescence experiments and molecular docking results showed that ZIF-8 did not change the type of quenching and interaction force between ciprofloxacin (CIP) and human serum albumin (HSA), but made the binding constant of HSA-CIP to be smaller, suggesting that ZIF-8 maybe accelerate the dissociation of CIP from HSA-CIP complex. Moreover, the effect of ZIF-8 on the physiological function of HSA was explored. Multi-spectroscopic methods revealed that ZIF-8 did not significantly alter the microenvironment of amino acid groups, but cause a slight decrease in the content of α-helical conformation, and a sparse and flexible structure of the protein backbone. These peculiarities might lead to the diminution of HSA's ability to control drugs. In short, ZIF-8 might enhance drug effect due to affecting the binding of drugs to proteins. However, the present study is only a preliminary investigation of the suitability of ZIF-8 as a drug carrier in vitro, and subsequent in vivo experimental studies will be required to further confirm the idea.

GPER binding site detection and description: A flavonoid-based docking and molecular dynamics simulations study

Flavonoids, a phenolic compounds class widely distributed in the plant kingdom, have attracted much interest for their implications on several health and disease processes. Usually, the consumption of this type of compounds is approximately 1 g/d, primarily obtained from cereals, chocolate, and dry legumes ensuring its beneficial role in maintaining the homeostasis of the human body. In this context, in cancer disease prominent data points to the role of flavonoids as adjuvant treatment aimed at the regression of the disease. GPER, an estrogen receptor on the cell surface, has been postulated as a probable orchestrator of the beneficial effects of several flavonoids through modulation/inhibition of various mechanisms that lead to cancer progression. Therefore, applying pocket and cavity protein detection and docking and molecular dynamics simulations (MD), we generate, from a cluster composed of 39 flavonoids, crucial insights into the potential role as GPER ligands, of Puerarin, Isoquercetin, Kaempferol 3-O-glucoside and Petunidin 3-O-glucoside, aglycones whose sugar moiety delimits a new described sugar-acceptor sub-cavity into the cavity binding site on the receptor, as well as of the probable activation mechanism of the receptor and the pivotal residues involved in it. Altogether, our results shed light on the potential use of the aforementioned flavonoids as GPER ligands and for further evaluations in in vitro and in vivo assays to elucidate their probable anti-cancer activity.

The green cupredoxin CopI is a multicopper protein able to oxidize Cu(I)

Anthropogenic activities in agriculture and health use the antimicrobial properties of copper. This has led to copper accumulation in the environment and contributed to the emergence of copper resistant microorganisms. Understanding bacterial copper homeostasis diversity is therefore highly relevant since it could provide valuable targets for novel antimicrobial treatments. The periplasmic CopI protein is a monodomain cupredoxin comprising several copper binding sites and is directly involved in copper resistance in bacteria. However, its structure and mechanism of action are yet to be determined. To study the different binding sites for cupric and cuprous ions and to understand their possible interactions, we have used mutants of the putative copper binding modules of CopI and spectroscopic methods to characterize their properties. We show that CopI is able to bind a cuprous ion in its central histidine/methionine-rich region and oxidize it thanks to its cupredoxin center. The resulting cupric ion can bind to a third site at the N-terminus of the protein. Nuclear magnetic resonance spectroscopy revealed that the central histidine/methionine-rich region exhibits a dynamic behavior and interacts with the cupredoxin binding region. CopI is therefore likely to participate in copper resistance by detoxifying the cuprous ions from the periplasm.

A conserved asparagine residue stabilizes iron binding in Manduca sexta transferrin-1

Transferrin 1 (Tsf1) is an insect-specific iron-binding protein that is abundant in hemolymph and other extracellular fluids. It binds iron tightly at neutral pH and releases iron under acidic conditions. Tsf1 influences the distribution of iron in the body and protects against infection. Elucidating the mechanisms by which Tsf1 achieves these functions will require an understanding of how Tsf1 binds and releases iron. Previously, crystallized Tsf1 from Manduca sexta was shown to have a novel type of iron coordination that involves four iron-binding ligands: two tyrosine residues (Tyr90 and Tyr204), a buried carbonate anion, and a solvent-exposed carbonate anion. The solvent-exposed carbonate anion was bound by a single amino acid residue, a highly conserved asparagine at position 121 (Asn121); thus, we predicted that Asn121 would be essential for high-affinity iron binding. To test this hypothesis, we analyzed the iron-binding and -release properties of five forms of recombinant Tsf1: wild-type, a Y90F/Y204F double mutant (negative control), and three Asn121 mutants (N121A, N121D and N121S). Each of the Asn121 mutants exhibited altered spectral properties, confirming that Asn121 contributes to iron coordination. The N121D and N121S mutations resulted in slightly lower affinity for iron, especially at acidic pH, while iron binding and release by the N121A mutant was indistinguishable from that of the wild-type protein. The surprisingly minor consequences of mutating Asn121, despite its high degree of conservation in diverse insect species, suggest that Asn121 may play a role that is essential in vivo but non-essential for high affinity iron binding in vitro.

Some pyrroles as inhibitors of the pentose phosphate pathways enzymes: An in vitro and molecular docking study

Glucose-6-phosphate dehydrogenase (G6PD) and 6-phosphogluconate dehydrogenase (6PGD) are pentose phosphate pathway enzymes. Compounds with a heterocyclic pyrrole ring system containing this atom can be derivatized with various functional groups into highly effective bioactive agents. In this study, pyrrole derivatives on these enzyme's activity were investigated. The IC50 values of different concentrations of pyrrole derivatives for G6PD were found in the range of 0.022-0.221 mM Ki values 0.021 ± 0.003-0.177 ± 0.021 and for 6PGD IC50 values 0.020-0.147, mM Ki values 0.013 ± 0.002-0.113 ± 0.030 mM. The 2-acetyl-1-methylpyrrole (1g) showed the best inhibition value for G6PD and 6PGD enzymes. In addition, in silico molecular docking experiments were performed to elucidate how these pyrrole derivatives (1a-g) interact with the binding sites of the target enzymes. The study's findings on pyrrole derivatives could be used to create innovative therapeutics that could be a treatment for many diseases, especially cancer manifestations.

Biophysical and structural characterization of tetramethrin serum protein complex and its toxicological implications

Tetramethrin (TMT) is a commonly used insecticide and has a carcinogenic and neurodegenerative effect on humans. The binding mechanism and toxicological implications of TMT to human serum albumin (HSA) were examined in this study employing a combination of biophysical and computational methods indicating moderate binding affinity and potential hepato and renal toxicity. Fluorescence quenching experiments showed that TMT binds to HSA with a moderate affinity, and the binding process was spontaneous and predominantly enthalpy-driven. Circular dichroism spectroscopy revealed that TMT binding did not induce any significant conformational changes in HSA, resulting in no changes in its alpha-helix content. The binding site and modalities of TMT interactions with HSA as computed by molecular docking and molecular dynamics simulations revealed that it binds to Sudlow site II of HSA via hydrophobic interactions through its dimethylcyclopropane carboxylate methyl propanyl group. The structural dynamics of TMT induce proper fit into the binding site creating increased and stabilizing interactions. Additionally, molecular mechanics-Poisson Boltzmann surface area calculations also indicated that non-polar and van der Waals were found to be the major contributors to the high binding free energy of the complex. Quantum mechanics (QM) revealed the conformational energies of the binding confirmation and the degree of deviation from the global minimum energy conformation of TMT. The results of this study provide a comprehensive understanding of the binding mechanism of TMT with HSA, which is important for evaluating the toxicity of this insecticide in humans.

Design, synthesis and biological evaluation of new 2,6-difluorinated phenyl 4-(2-oxoimidazolidin-1-yl)benzenesulfonates as new antimicrotubule agents

We previously discovered a novel family of antimicrotubule agents designated as phenyl 4-(2-oxoimidazolidin-1-yl)benzenesulfonates (PIB-SOs). In this study, we evaluated the effect of the difluorination of the aromatic ring bearing the imidazolidin-2-one moiety (ring A) at positions 3, 5 and 2, 6 on their antiproliferative activity on four cancer cell lines, their ability to disrupt the microtubules and their toxicity toward chick embryos. We thus synthesized, characterized and biologically evaluated 24 new difluorinated PIB-SO derivatives designated as phenyl 3,5-difluoro-4-(2-oxoimidazolidin-1-yl)benzenesulfonates (3,5-PFB-SOs, 4-15) and phenyl 2,6-difluoro-4-(2-oxoimidazolidin-1-yl)benzenesulfonates (2,6-PFB-SOs, 16-27). The concentration of the drug required to inhibit cell growth by 50% (IC50) of 3,5-PFB-SOs is over 1000 nM while most of 2,6-PFB-SOs exhibit IC50 in the nanomolar range (23-900 nM). Furthermore, the most potent 2,6-PFB-SOs 19, 26 and 27 arrest the cell cycle progression in G2/M phase, induce cytoskeleton disruption and impair microtubule polymerization. Docking studies also show that the most potent 2,6-PFB-SOs 19, 21, 24, 26 and 27 have binding affinity toward the colchicine-binding site (C-BS). Moreover, their antiproliferative activity is not affected by antimicrotubule- and multidrug-resistant cell lines. Besides, they exhibit improved in vitro hepatic stability in the mouse, rat and human microsomes compared to their non-fluorinated counterparts. They also showed theoretical pharmacokinetic, physicochemical and drug-like properties suited for further in vivo assays. In addition, they exhibit low to no systemic toxicity toward chick embryos. Finally, our study evidences that PIB-SOs must be fluorinated in specific positions on ring A to maintain both their antiproliferative activity and their biological activity toward microtubules.

Human Serum Albumin Immobilized On Magnetizable Beads: A Rapid Method for Compound HSA Binding Study

Immobilized human serum albumin (HSA) was developed by coupling His-tagged HSA onto Ni2+-coupled magnetizable beads (HSA-beads), allowing the HSA to be easily removed from incubation components. The HSA-beads system provides a rapid and convenient method to study HSA compound binding. In this study, the HSA-beads system was characterized and evaluated as a tool for assessing compound HSA binding properties. The free fraction (fu) values of test compounds measured using HSA-beads were comparable to those determined by equilibrium dialysis (ED), which is commonly used to evaluate albumin binding in vitro. The equilibrium dissociation constant (Kd) values determined for a series of compounds using the HSA-beads method demonstrated good correlation with literature data. This good correlation also suggests that the binding of His-HSA to the beads does not impact the conformations of the two compound binding sites of HSA, as the range of compounds tested encompassed binding to both sites. Furthermore, the Kd values of representative compounds itraconazole and BIRT2584 that were difficult to assess using ED, due to significant cellulose membrane adsorption, were successfully determined. The HSA-beads provide several advantages over ED, such as simple preparation, short assay incubation duration, and the ability to quantify both free and HSA-bound species of the test compound, facilitated by the simple separation of HSA-beads from the solution phase using a magnetic field. These properties render the HSA-beads method suitable for high-throughput studies on compound HSA binding.

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

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

Identification of two novel chemical classes of Autotaxin (ATX) inhibitors using Enalos Asclepios KNIME nodes

Autotaxin is a secreted lysophospholipase D which is a member of the ectonucleotide pyrophosphatase/phosphodiesterase family converting extracellular lysophosphatidylcholine and other non-choline lysophospholipids, such as lysophosphatidylethanolamine and lysophosphatidylserine, to the lipid mediator lysophosphatidic acid. Autotaxin is implicated in various fibroproliferative diseases including interstitial lung diseases, such as idiopathic pulmonary fibrosis and hepatic fibrosis, as well as in cancer. In this study, we present an effort of identifying ATX inhibitors that bind to allosteric ATX binding sites using the Enalos Asclepios KNIME Node. All the available PDB crystal structures of ATX were collected, prepared, and aligned. Visual examination of these structures led to the identification of four crystal structures of human ATX co-crystallized with four known inhibitors. These inhibitors bind to five binding sites with five different binding modes. These five binding sites were thereafter used to virtually screen a compound library of 14,000 compounds to identify molecules that bind to allosteric sites. Based on the binding mode and interactions, the docking score, and the frequency that a compound comes up as a top-ranked among the five binding sites, 24 compounds were selected for in vitro testing. Finally, two compounds emerged with inhibitory activity against ATX in the low micromolar range, while their mode of inhibition and binding pattern were also studied. The two derivatives identified herein can serve as "hits" towards developing novel classes of ATX allosteric inhibitors.

The noncanonical nucleotide binding site 1 of the bile salt export pump is optimized for proper function of the transporter

The bile salt export pump (ABCB11/BSEP) is a hepatocyte plasma membrane-resident protein translocating bile salts into bile canaliculi. The sequence alignment of the four full-length transporters of the ABCB subfamily (ABCB1, ABCB4, ABCB5 and ABCB11) indicates that the NBD-NBD contact interface of ABCB11 differs from that of other members in only four residues. Notably, these are all located in the noncanonical nucleotide binding site 1 (NBS1). Substitution of all four deviant residues with canonical ones (quadruple mutant) significantly decreased the transport activity of the protein. In this study, we mutated two deviant residues in the signature sequence to generate a double mutant (R1221G/E1223Q). Furthermore, a triple mutant (E502S/R1221G/E1223Q) was generated, in which the deviant residues of the signature sequence and Q-loop were mutated concurrently to canonical residues. The double and triple mutants showed 80% and 60%, respectively, of the activity of wild-type BSEP. As expected, an increasing number of mutations gradually impair transport as an intricate network of interactions within the ABC proteins ensures proper functioning.

Computational insights into diverse binding modes of the allosteric modulator and their regulation on dopamine D1 receptor

Allosteric drugs hold the promise of addressing many challenges in the current drug development of GPCRs. However, the molecular mechanism underlying their allosteric modulations remain largely elusive. The dopamine D1 receptor (DRD1), a member of Class A GPCRs, is critical for treating psychiatric disorders, and LY3154207 serves as its promising positive allosteric modulator (PAM). In the work, we utilized extensive Gaussian-accelerated molecular dynamics simulations (a total of 41μs) for the first time probe the diverse binding modes of the allosteric modulator and their regulation effects, based on the DRD1 and LY3154207 as representative. Our simulations identify four binding modes of LY3154207 (one boat mode, two metastable vertical modes and a novel cleft-anchored mode), in which the boat mode is the most stable while there three modes are similar in the stability. However, it is interesting to observed that the most stable boat mode inversely exhibits the weakest positive allosteric effect on influencing the orthosteric ligand binding and maintaining the activity of the transducer binding site. It should result from its induced weaker correlation between the allosteric site and the orthosteric site, and between the orthosteric site and the transducer binding site than the other three binding modes, as well as its weakened interaction between a crucial activation-related residue (S2025.46) and the orthosteric ligand (dopamine). Overall, the work offers atomic-level information to advance our understanding of the complex allosteric regulation on GPCRs, which is beneficial to the allosteric modulator design and development.

Classification of binding property of amyloid β to lipid membranes: Membranomic research using quartz crystal microbalance combined with the immobilization of lipid planar membranes

A biomembrane-related fibrillogenesis of Amyloid β from Alzheimer' disease (Aβ) is closely related to its accumulation behavior. A binding property of Aβ peptides from Alzheimer' disease to lipid membranes was then classified by a quartz crystal microbalance (QCM) method combined with an immobilization technique using thiol self-assembled membrane. The accumulated amounts of Aβ, Δfmax, was determined from the measurement of the maximal frequency reduction using QCM. The plots of Δfmax to Aβ concentration gave the slope and saturated value of Δfmax, (Δfmax)sat that are the parameters for binding property of Aβ to lipid membranes. Therefore, the Aβ-binding property on lipid membranes was classified by the slope and (Δfmax)sat. The plural lipid system was described as X + Y where X = L1, L1/L2, and L1/L2/L3. The slope and (Δfmax)sat values plotted as a function of mixing ratio of Y to X was classified on a basis of the lever principle (LP). The LP violation observed in both parameters resulted from the formation of the crevice or pothole, as Aβ-specific binding site, generated at the boundary between ld and lo phases. The LP violation observed only in the slope resulted from glycolipid-rich domain acting as Aβ-specific binding site. Furthermore, lipid planar membranes indicating strong LP violation favored strong fibrillogenesis. Especially, lipid planar membranes indicating the LP violation only in the slope induced lateral aggregated and spherulitic fibrillar aggregates. Thus, the classification of Aβ binding property on lipid membranes appeared to be related to the fibrillogenesis with a certain morphology.

Modular preparation of biphenyl triazoles via click chemistry as non-competitive hyaluronidase inhibitors

Hyaluronidase is a promising target in drug discovery, given its overexpression in a range of physiological and pathological processes, including tumor migration, skin aging, sagging, and wrinkling, as well as inflammation and bacterial infections. In this study, to identify novel hyaluronidase inhibitors, we applied click chemistry for the modular synthesis of 370 triazoles in 96-well plates, starting with biphenyl azide. Utilizing an optimized turbidimetric screening assay in microplates, we identified Fmoc-containing triazoles 5 and 6, as well as quinoline-containing triazoles 15 and 16, as highly effective hyaluronidase inhibitors. Subsequent research indicated that these triazoles potentially interact with a novel binding site of hyaluronidase. Notably, these inhibitors displayed minimal cytotoxicity and showed promising anti-inflammatory effects in LPS-stimulated macrophages. Remarkably, compound 6 significantly reduced NO release by 74 % at a concentration of 20 μM.

Probing the toxic effect of chlorpyrifos as an environmental pollutant on the structure and biological activity of lysozyme under physiological conditions

The pervasive use of pesticides like chlorpyrifos (CPY) has been associated with deleterious effects on biomolecules, posing significant risks to environmental integrity, public health, and overall ecosystem equilibrium. Accordingly, in this study, we investigated the potential binding interaction between the well-conserved enzyme, lysozyme (LSZ), and CPY through various spectroscopic techniques and molecular modeling. The UV-vis absorption and fluorescence experiments confirmed the complex formation and static quenching of the intrinsic fluorescence intensity. LSZ revealed a singular binding site for CPY, with binding constants around 105 M-1 across different temperature ranges. Analysis of thermodynamic parameters showed the spontaneous nature of the complexation process, while also revealing the pivotal role of hydrophobic interactions in stabilizing the LSZ-CPY system. According to circular dichroism and Fourier transform infrared studies, CPY binding changed the secondary structure of LSZ by boosting α-helix presence and reducing the levels of β-sheet and β-turn content. Further, CPY decreased the stability and activity of LSZ. Computational docking delineated the specific and highly preferred binding site of CPY within the structure of LSZ. Molecular dynamic simulation indicated the enduring stability of the LSZ/CPY complex and revealed structural modifications in the LSZ after binding with CPY. This research provides a detailed understanding of the intermolecular dynamics between CPY and LSZ, concurrently elucidating the molecular-level implications for the potential hazards of pesticides in the natural environment.

Unraveling the Ligand-Binding Sites of CYP3A4 by Molecular Dynamics Simulations with Solvent Probes

Cytochrome P450 3A4 (CYP3A4) is one of the most important drug-metabolizing enzymes in the human body and is well known for its complicated, atypical kinetic characteristics. The existence of multiple ligand-binding sites in CYP3A4 has been widely recognized as being capable of interfering with the active pocket through allosteric effects. The identification of ligand-binding sites other than the canonical active site above the heme is especially important for understanding the atypical kinetic characteristics of CYP3A4 and the intriguing association between the ligand and the receptor. In this study, we first employed mixed-solvent molecular dynamics (MixMD) simulations coupled with the online computational predictive tools to explore potential ligand-binding sites in CYP3A4. The MixMD approach demonstrates better performance in dealing with the receptor flexibility compared with other computational tools. From the sites identified by MixMD, we then picked out multiple sites for further exploration using ensemble docking and conventional molecular dynamics (cMD) simulations. Our results indicate that three extra sites are suitable for ligand binding in CYP3A4, including one experimentally confirmed site and two novel sites.

Structure-Based Protein Assembly Simulations Including Various Binding Sites and Conformations

Many biological functions are mediated by large complexes formed by multiple proteins and other cellular macromolecules. Recent progress in experimental structure determination, as well as in integrative modeling and protein structure prediction using deep learning approaches, has resulted in a rapid increase in the number of solved multiprotein assemblies. However, the assembly process of large complexes from their components is much less well-studied. We introduce a rapid computational structure-based (SB) model, GoCa, that allows to follow the assembly process of large multiprotein complexes based on a known native structure. Beyond existing SB Go̅-type models, it distinguishes between intra- and intersubunit interactions, allowing us to include coupled folding and binding. It accounts automatically for the permutation of identical subunits in a complex and allows the definition of multiple minima (native) structures in the case of proteins that undergo global transitions during assembly. The model is successfully tested on several multiprotein complexes. The source code of the GoCa program including a tutorial is publicly available on Github: https://github.com/ZachariasLab/GoCa. We also provide a web source that allows users to quickly generate the necessary input files for a GoCa simulation: https://goca.t38webservices.nat.tum.de.

Binding Sites of Bicarbonate in Phosphoenolpyruvate Carboxylase

Phosphoenolpyruvate carboxylase (PEPC) is used in plant metabolism for fruit maturation or seed development as well as in the C4 and crassulacean acid metabolism (CAM) mechanisms in photosynthesis, where it is used for the capture of hydrated CO2 (bicarbonate). To find the yet unknown binding site of bicarbonate in this enzyme, we have first identified putative binding sites with nonequilibrium molecular dynamics simulations and then ranked these sites with alchemical free energy calculations with corrections of computational artifacts. Fourteen pockets where bicarbonate could bind were identified, with three having realistic binding free energies with differences with the experimental value below 1 kcal/mol. One of these pockets is found far from the active site at 14 Å and predicted to be an allosteric binding site. In the two other binding sites, bicarbonate is in direct interaction with the magnesium ion; neither sequence alignment nor the study of mutant K606N allowed to discriminate between these two pockets, and both are good candidates as the binding site of bicarbonate in phosphoenolpyruvate carboxylase.

Deciphering S100B Allosteric Signaling: The Role of a Peptide Target, TRTK-12, as an Ensemble Modulator

Allostery is an essential biological phenomenon in which perturbation at one site in a biomolecule elicits a functional response at a distal location(s). It is integral to biological processes, such as cellular signaling, metabolism, and transcription regulation. Understanding allostery is also crucial for rational drug discovery. In this work, we focus on an allosteric S100B protein that belongs to the S100 class of EF-hand Ca2+-binding proteins. The Ca2+-binding affinity of S100B is modulated allosterically by TRTK-12 peptide binding 25 Å away from the Ca2+-binding site. We investigated S100B allostery by carrying out nuclear magnetic resonance (NMR) measurements along with microsecond-long molecular dynamics (MD) simulations on S100B/Ca2+ with/without TRTK-12 at different NaCl salt concentrations. NMR HSQC results show that TRTK-12 reorganizes how S100B/Ca2+ responds to different salt concentrations at both orthosteric and allosteric sites. The MD data suggest that TRTK-12 breaks the dynamic aromatic and hydrogen-bond interactions (not observed in X-ray crystallographic structures) between the hinge/helix and Ca2+-binding EF-hand loop of the two subunits in the homodimeric protein. This triggers rearrangement in the protein network architectures and leads to allosteric communication. Finally, computational studies of S100B at distinct ionic strengths suggest that ligand-bound species are more robust to the changing environment relative to the S100B/Ca2+ complex.

Calcium-Driven In Silico Inactivation of a Human Olfactory Receptor

Conformational changes as well as molecular determinants related to the activation and inactivation of olfactory receptors are still poorly understood due to the intrinsic difficulties in the structural determination of this GPCR family. Here, we perform, for the first time, the in silico inactivation of human olfactory receptor OR51E2, highlighting the possible role of calcium in this receptor state transition. Using molecular dynamics simulations, we show that a divalent ion in the ion binding site, coordinated by two acidic residues at positions 2.50 and 3.39 conserved across most ORs, stabilizes the receptor in its inactive state. In contrast, protonation of the same two acidic residues is not sufficient to drive inactivation within the microsecond timescale of our simulations. Our findings suggest a novel molecular mechanism for OR inactivation, potentially guiding experimental validation and offering insights into the possible broader role of divalent ions in GPCR signaling.

Drug-Online: an online platform for drug-target interaction, affinity, and binding sites identification using deep learning

BACKGROUND

Accurately identifying drug-target interaction (DTI), affinity (DTA), and binding sites (DTS) is crucial for drug screening, repositioning, and design, as well as for understanding the functions of target. Although there are a few online platforms based on deep learning for drug-target interaction, affinity, and binding sites identification, there is currently no integrated online platforms for all three aspects.

RESULTS

Our solution, the novel integrated online platform Drug-Online, has been developed to facilitate drug screening, target identification, and understanding the functions of target in a progressive manner of "interaction-affinity-binding sites". Drug-Online platform consists of three parts: the first part uses the drug-target interaction identification method MGraphDTA, based on graph neural networks (GNN) and convolutional neural networks (CNN), to identify whether there is a drug-target interaction. If an interaction is identified, the second part employs the drug-target affinity identification method MMDTA, also based on GNN and CNN, to calculate the strength of drug-target interaction, i.e., affinity. Finally, the third part identifies drug-target binding sites, i.e., pockets. The method pt-lm-gnn used in this part is also based on GNN.

CONCLUSIONS

Drug-Online is a reliable online platform that integrates drug-target interaction, affinity, and binding sites identification. It is freely available via the Internet at http://39.106.7.26:8000/Drug-Online/ .

Protein characteristics substantially influence the propensity of activity cliffs among kinase inhibitors

Activity cliffs (ACs) are pairs of structurally similar molecules with significantly different affinities for a biotarget, posing a challenge in computer-assisted drug discovery. This study focuses on protein kinases, significant therapeutic targets, with some exhibiting ACs while others do not despite numerous inhibitors. The hypothesis that the presence of ACs is dependent on the target protein and its complete structural context is explored. Machine learning models were developed to link protein properties to ACs, revealing specific tripeptide sequences and overall protein properties as critical factors in ACs occurrence. The study highlights the importance of considering the entire protein matrix rather than just the binding site in understanding ACs. This research provides valuable insights for drug discovery and design, paving the way for addressing ACs-related challenges in modern computational approaches.

A comprehensive dataset of protein-protein interactions and ligand binding pockets for advancing drug discovery

This dataset represents a collection of pocket-centric structural data related to protein-protein interactions (PPIs) and PPI-related ligand binding sites. The dataset includes high-quality structural information on more than 23,000 pockets, 3,700 proteins on more than 500 organisms, and nearly 3500 ligands that can aid researchers in the fields of bioinformatics, structural biology, and drug discovery. It encompasses a diverse set of PPI complexes with more than 1,700 unique protein families including some with associated ligands, enabling detailed investigations into molecular interactions at the atomic level. This article introduces an indispensable resource designed to unlock the full potential of PPIs while pioneering a novel metric for pocket similarity for hypothesizing protein partners repurposing.

MMGAT: a graph attention network framework for ATAC-seq motifs finding

BACKGROUND

Motif finding in Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq) data is essential to reveal the intricacies of transcription factor binding sites (TFBSs) and their pivotal roles in gene regulation. Deep learning technologies including convolutional neural networks (CNNs) and graph neural networks (GNNs), have achieved success in finding ATAC-seq motifs. However, CNN-based methods are limited by the fixed width of the convolutional kernel, which makes it difficult to find multiple transcription factor binding sites with different lengths. GNN-based methods has the limitation of using the edge weight information directly, makes it difficult to aggregate the neighboring nodes' information more efficiently when representing node embedding.

RESULTS

To address this challenge, we developed a novel graph attention network framework named MMGAT, which employs an attention mechanism to adjust the attention coefficients among different nodes. And then MMGAT finds multiple ATAC-seq motifs based on the attention coefficients of sequence nodes and k-mer nodes as well as the coexisting probability of k-mers. Our approach achieved better performance on the human ATAC-seq datasets compared to existing tools, as evidenced the highest scores on the precision, recall, F1_score, ACC, AUC, and PRC metrics, as well as finding 389 higher quality motifs. To validate the performance of MMGAT in predicting TFBSs and finding motifs on more datasets, we enlarged the number of the human ATAC-seq datasets to 180 and newly integrated 80 mouse ATAC-seq datasets for multi-species experimental validation. Specifically on the mouse ATAC-seq dataset, MMGAT also achieved the highest scores on six metrics and found 356 higher-quality motifs. To facilitate researchers in utilizing MMGAT, we have also developed a user-friendly web server named MMGAT-S that hosts the MMGAT method and ATAC-seq motif finding results.

CONCLUSIONS

The advanced methodology MMGAT provides a robust tool for finding ATAC-seq motifs, and the comprehensive server MMGAT-S makes a significant contribution to genomics research. The open-source code of MMGAT can be found at https://github.com/xiaotianr/MMGAT , and MMGAT-S is freely available at https://www.mmgraphws.com/MMGAT-S/ .

Analysis of REST binding sites with canonical and non-canonical motifs in human cell lines

BACKGROUND

Repressor element 1 (RE1) silencing transcription factor (REST) is a transcriptional repressor abundantly expressed in aging human brains. It is known to regulate genes associated with oxidative stress, inflammation, and neurological disorders by binding to a canonical form of sequence motif and its non-canonical variations. Although analysis of genomic sequence motifs is crucial to understand transcriptional regulation by transcription factors (TFs), a comprehensive characterization of various forms of RE1 motifs in human cell lines has not been performed.

RESULTS

Here, we analyzed 23 ENCODE REST ChIP-seq datasets from diverse human cell lines and identified a non-redundant set of 68,975 loci with ChIP-seq peaks. Our systematic characterization of these binding sites revealed that the canonical form of REST binding motif was found primarily in ChIP-seq peaks shared across multiple cell lines, while non-canonical forms of motifs were identified in both cell-line-specific binding sites and those shared across cell lines. Remarkably, we observed a notable prevalence of non-canonical motifs that corresponded to half segments of the canonical motif. Furthermore, our analysis unveiled the presence of cell-line-specific REST binding patterns, as evidenced by the clustering of ChIP-seq experiments according to their respective cell lines. This observation underscores the cell-line specificity of REST binding at certain genomic loci, implying intricate cell-line-specific regulatory mechanisms.

CONCLUSIONS

Overall, our study provides a comprehensive characterization of REST binding motifs in human cell lines and genome-wide RE1 motif profiles. These findings contribute to a deeper understanding of REST-mediated transcriptional regulation and highlight the importance of considering cell-line-specific effects in future investigations.

Genome-scale annotation of protein binding sites via language model and geometric deep learning

Revealing protein binding sites with other molecules, such as nucleic acids, peptides, or small ligands, sheds light on disease mechanism elucidation and novel drug design. With the explosive growth of proteins in sequence databases, how to accurately and efficiently identify these binding sites from sequences becomes essential. However, current methods mostly rely on expensive multiple sequence alignments or experimental protein structures, limiting their genome-scale applications. Besides, these methods haven't fully explored the geometry of the protein structures. Here, we propose GPSite, a multi-task network for simultaneously predicting binding residues of DNA, RNA, peptide, protein, ATP, HEM, and metal ions on proteins. GPSite was trained on informative sequence embeddings and predicted structures from protein language models, while comprehensively extracting residual and relational geometric contexts in an end-to-end manner. Experiments demonstrate that GPSite substantially surpasses state-of-the-art sequence-based and structure-based approaches on various benchmark datasets, even when the structures are not well-predicted. The low computational cost of GPSite enables rapid genome-scale binding residue annotations for over 568,000 sequences, providing opportunities to unveil unexplored associations of binding sites with molecular functions, biological processes, and genetic variants. The GPSite webserver and annotation database can be freely accessed at https://bio-web1.nscc-gz.cn/app/GPSite.

Tissue-specific atlas of trans-models for gene regulation elucidates complex regulation patterns

BACKGROUND

Deciphering gene regulation is essential for understanding the underlying mechanisms of healthy and disease states. While the regulatory networks formed by transcription factors (TFs) and their target genes has been mostly studied with relation to cis effects such as in TF binding sites, we focused on trans effects of TFs on the expression of their transcribed genes and their potential mechanisms.

RESULTS

We provide a comprehensive tissue-specific atlas, spanning 49 tissues of TF variations affecting gene expression through computational models considering two potential mechanisms, including combinatorial regulation by the expression of the TFs, and by genetic variants within the TF. We demonstrate that similarity between tissues based on our discovered genes corresponds to other types of tissue similarity. The genes affected by complex TF regulation, and their modelled TFs, were highly enriched for pharmacogenomic functions, while the TFs themselves were also enriched in several cancer and metabolic pathways. Additionally, genes that appear in multiple clusters are enriched for regulation of immune system while tissue clusters include cluster-specific genes that are enriched for biological functions and diseases previously associated with the tissues forming the cluster. Finally, our atlas exposes multilevel regulation across multiple tissues, where TFs regulate other TFs through the two tested mechanisms.

CONCLUSIONS

Our tissue-specific atlas provides hierarchical tissue-specific trans genetic regulations that can be further studied for association with human phenotypes.

Structural, Spectroscopic, and Computational Insights from Canavanine-Bound and Two Catalytically Compromised Variants of the Ethylene-Forming Enzyme

The ethylene-forming enzyme (EFE) is an Fe(II), 2-oxoglutarate (2OG), and l-arginine (l-Arg)-dependent oxygenase that either forms ethylene and three CO2/bicarbonate from 2OG or couples the decarboxylation of 2OG to C5 hydroxylation of l-Arg. l-Arg binds with C5 toward the metal center, causing 2OG to change from monodentate to chelate metal interaction and OD1 to OD2 switch of D191 metal coordination. We applied anaerobic UV-visible spectroscopy, X-ray crystallography, and computational approaches to three EFE systems with high-resolution structures. The ineffective l-Arg analogue l-canavanine binds to the EFE with O5 pointing away from the metal center while promoting chelate formation by 2OG but fails to switch the D191 metal coordination from OD1 to OD2. Substituting alanine for R171 that interacts with 2OG and l-Arg inactivates the protein, prevents metal chelation by 2OG, and weakens l-Arg binding. The R171A EFE had electron density at the 2OG binding site that was identified by mass spectrometry as benzoic acid. The substitution by alanine of Y306 in the EFE, a residue 12 Å away from the catalytic metal center, generates an interior cavity that leads to multiple local and distal structural changes that reduce l-Arg binding and significantly reduce the enzyme activity. Flexibility analyses revealed correlated and anticorrelated motions in each system, with important distinctions from the wild-type enzyme. In combination, the results are congruent with the currently proposed enzyme mechanism, reinforce the importance of metal coordination by OD2 of D191, and highlight the importance of the second coordination sphere and longer range interactions in promoting EFE activity.

Beyond copper: examining the significance of His-mutations in mycobacterial GroEL1 HRCT for Ni(II) complex stability and formation

Recently, we have studied the coordination chemistry of the Cu(II)-histidine-rich C-terminal tail (HRCT) complex of the mycobacterial GroEL1 protein. The structure of this domain differs significantly compared to the well-known methionine-glycine-rich GroEL chaperonin - it was predicted that mycobacterial GroEL1 could play a significant role in the metal homeostasis of Mycobacteria, especially copper. However, we found that this particular domain's pattern also repeats in a number of Ni(II)-binding proteins. Here, we present the studies concerning the properties of GroEL1 HRCT as a ligand for Ni(II) ions. For this purpose, we chose eight model peptides: L1 - Ac-DHDHHHGHAH, L2 - Ac-DKPAKAEDHDHHHGHAH, and 6 mutants of the latter in the pH range of 2-11. We examined the stoichiometry, stability, and spectroscopic features of copper complexes. We noticed that similar to the Cu(II)-complex, the presence of a Lys5 residue significantly increases the stability of the system. The impact of His mutations was also examined and carefully studied using NMR spectroscopy. His9 and His13 are the crucial residues for Ni(II) binding, whereas His12 has minimal relevance in complex formation.

Bi(III) Binding Stoichiometry and Domain-Specificity Differences Between Apo and Zn(II)-bound Human Metallothionein 1a

Bismuth is a xenobiotic metal with a high affinity to sulfur that is used in a variety of therapeutic applications. Bi(III) induces the cysteine-rich metallothionein (MT), a protein known to form two-domain cluster structures with certain metals such as Zn(II), Cd(II), or Cu(I). The binding of Bi(III) to MTs has been previously studied, but there are conflicting reports on the stoichiometry and binding pathway, which appear to be highly dependent on pH and initial metal-loading status of the MT. Additionally, domain specificity has not been thoroughly investigated. In this paper, ESI-MS was used to determine the binding constants of [Bi(EDTA)]- binding to apo-MT1a and its individual αMT fragment. The results were compared to previous experiments using βMT1a and βαMT3. Domain specificity was investigated using proteolysis methods and the initial cooperatively formed Bi2MT was found to bind to cysteines that spanned across the traditional metal binding domain regions. Titrations of [Bi(EDTA)]- into Zn7MT were performed and were found to result in a maximum stoichiometry of Bi7MT, contrasting the Bi6MT formed when [Bi(EDTA)]- was added to apo-MT. These results show that the initial structure of the apo-MT determines the stoichiometry of new incoming metals and explains the previously observed differences in stoichiometry.

Insight into the molecular interaction between the anticancer drug, enzalutamide and human alpha-2-macroglobulin: Biochemical and biophysical approach

The drug pharmacokinetics is affected upon binding with proteins, thus making drug-protein interactions crucial. This study investigated the interaction between enzalutamide and human major antiproteinase alpha-2-macroglobulin (α2M) by using multi spectroscopic and calorimetric techniques. The spectroscopic techniques such as circular dichroism (CD), intrinsic fluorescence, and UV-visible absorption were used to determine the mechanism of enzalutamide-α2M interaction. Studies on the quenching of fluorescence at three different temperatures showed that the enzalutamide-α2M complex is formed through static quenching mechanism. The change in microenvironment around tyrosine residues in protein was detected through synchronised fluorescence. The secondary structure of α2M was slightly altered by enzalutamide according to far UV-CD spectral analysis. Changes in position of amide I band in FTIR spectra further confirm the secondary structural alteration in α2M. According to thermodynamic characteristics such as fluorescence quenching and isothermal titration calorimetry (ITC), hydrogen bonds and hydrophobic interactions were involved in the interaction machanism. The ITC reiterated the exothermic and spontaneous nature of the interaction. The lower proteinase inhibitory activity of the α2M-enzalutamide conjugate as reflects the disruption of the native α2M structure upon interaction with enzalutamide.

Sm-like protein Rof inhibits transcription termination factor ρ by binding site obstruction and conformational insulation

Transcription termination factor ρ is a hexameric, RNA-dependent NTPase that can adopt active closed-ring and inactive open-ring conformations. The Sm-like protein Rof, a homolog of the RNA chaperone Hfq, inhibits ρ-dependent termination in vivo but recapitulation of this activity in vitro has proven difficult and the precise mode of Rof action is presently unknown. Here, our cryo-EM structures of ρ-Rof and ρ-RNA complexes show that Rof undergoes pronounced conformational changes to bind ρ at the protomer interfaces, undercutting ρ conformational dynamics associated with ring closure and occluding extended primary RNA-binding sites that are also part of interfaces between ρ and RNA polymerase. Consistently, Rof impedes ρ ring closure, ρ-RNA interactions and ρ association with transcription elongation complexes. Structure-guided mutagenesis coupled with functional assays confirms that the observed ρ-Rof interface is required for Rof-mediated inhibition of cell growth and ρ-termination in vitro. Bioinformatic analyses reveal that Rof is restricted to Pseudomonadota and that the ρ-Rof interface is conserved. Genomic contexts of rof differ between Enterobacteriaceae and Vibrionaceae, suggesting distinct modes of Rof regulation. We hypothesize that Rof and other cellular anti-terminators silence ρ under diverse, but yet to be identified, stress conditions when unrestrained transcription termination by ρ may be detrimental.

Unravelling carbohydrate binding module 21 (CBM21) dynamics of interaction with amylose

The carbohydrate binding module 21 (CBM21) from Rhizopus oryzae is a dual-site CBM proposed to disrupt polysaccharide structures. Additionally, it serves as a purification tag in industry. CBM21 crystal structure features a Glc residue in an unusual 1S3 conformation, whose relevance for the CBM mechanism of action is unclear. In this context, we seek to contribute for the understanding of CBM21 mechanism of action by: i) investigating the role of the 1S3 conformation on carbohydrate recognition, and ii) characterize the protein-carbohydrate binding dynamics using molecular dynamics and metadynamics simulations at MM and QM/MM levels. Results indicate the 1S3 Glc conformation is unlikely to occur under biological conditions, being originated from the crystallographic environment. CBM21 binding to small ligands appears transient and unstable, while protein dimerization and polysaccharide chain size influence complex stability. In interactions with amylose, CBM21 exhibits a repeated unbinding followed by re-binding, while simultaneously alternating between binding sites I and II. These results suggest that CBM21 acts through transient interactions, directing carbohydrates to the catalytic center rather than forming strong and long-lasting bonds with carbohydrates. Accordingly, we expect such atomistic depiction of CBM21 mechanism could aid in CBM design targeting biotechnological applications.

Reciprocating RNA Polymerase batters through roadblocks

RNA polymerases must transit through protein roadblocks to produce full-length transcripts. Here we report real-time measurements of Escherichia coli RNA polymerase passing through different barriers. As intuitively expected, assisting forces facilitated, and opposing forces hindered, RNA polymerase passage through lac repressor protein bound to natural binding sites. Force-dependent differences were significant at magnitudes as low as 0.2 pN and were abolished in the presence of the transcript cleavage factor GreA, which rescues backtracked RNA polymerase. In stark contrast, opposing forces promoted passage when the rate of RNA polymerase backtracking was comparable to, or faster than the rate of dissociation of the roadblock, particularly in the presence of GreA. Our experiments and simulations indicate that RNA polymerase may transit after roadblocks dissociate, or undergo cycles of backtracking, recovery, and ramming into roadblocks to pass through. We propose that such reciprocating motion also enables RNA polymerase to break protein-DNA contacts that hold RNA polymerase back during promoter escape and RNA chain elongation. This may facilitate productive transcription in vivo.

NAC and Zuotin/Hsp70 chaperone systems coexist at the ribosome tunnel exit in vivo

The area surrounding the tunnel exit of the 60S ribosomal subunit is a hub for proteins involved in maturation and folding of emerging nascent polypeptide chains. How different factors vie for positioning at the tunnel exit in the complex cellular environment is not well understood. We used in vivo site-specific cross-linking to approach this question, focusing on two abundant factors-the nascent chain-associated complex (NAC) and the Hsp70 chaperone system that includes the J-domain protein co-chaperone Zuotin. We found that NAC and Zuotin can cross-link to each other at the ribosome, even when translation initiation is inhibited. Positions yielding NAC-Zuotin cross-links indicate that when both are present the central globular domain of NAC is modestly shifted from the mutually exclusive position observed in cryogenic electron microscopy analysis. Cross-linking results also suggest that, even in NAC's presence, Hsp70 can situate in a manner conducive for productive nascent chain interaction-with the peptide binding site at the tunnel exit and the J-domain of Zuotin appropriately positioned to drive stabilization of nascent chain binding. Overall, our results are consistent with the idea that, in vivo, the NAC and Hsp70 systems can productively position on the ribosome simultaneously.

5-Aminothiazoles Reveal a New Ligand-Binding Site on Prolyl Oligopeptidase Which is Important for Modulation of Its Protein-Protein Interaction-Derived Functions

A series of novel 5-aminothiazole-based ligands for prolyl oligopeptidase (PREP) comprise selective, potent modulators of the protein-protein interaction (PPI)-mediated functions of PREP, although they are only weak inhibitors of the proteolytic activity of PREP. The disconnected structure-activity relationships are significantly more pronounced for the 5-aminothiazole-based ligands than for the earlier published 5-aminooxazole-based ligands. Furthermore, the stability of the 5-aminothiazole scaffold allowed exploration of wider substitution patterns than that was possible with the 5-aminooxazole scaffold. The intriguing structure-activity relationships for the modulation of the proteolytic activity and PPI-derived functions of PREP were elaborated by presenting a new binding site for PPI modulating PREP ligands, which was initially discovered using molecular modeling and later confirmed through point mutation studies. Our results suggest that this new binding site on PREP is clearly more important than the active site of PREP for the modulation of its PPI-mediated functions.