Multi-level bioinformatics resources support drug target discovery of protein-protein interactions

Drug discovery often begins with a new target. Protein-protein interactions (PPIs) are crucial to multitudinous cellular processes and offer a promising avenue for drug-target discovery. PPIs are characterized by multi-level complexity: at the protein level, interaction networks can be used to identify potential targets, whereas at the residue level, the details of the interactions of individual PPIs can be used to examine a target's druggability. Much great progress has been made in target discovery through multi-level PPI-related computational approaches, but these resources have not been fully discussed. Here, we systematically survey bioinformatics tools for identifying and assessing potential drug targets, examining their characteristics, limitations and applications. This work will aid the integration of the broader protein-to-network context with the analysis of detailed binding mechanisms to support the discovery of drug targets.

Machine learning and classical MD simulation to identify inhibitors against the P37 envelope protein of monkeypox virus

Monkeypox virus (MPXV) outbreak is a serious public health concern that requires international attention. P37 of MPXV plays a pivotal role in DNA replication and acts as one of the promising targets for antiviral drug design. In this study, we intent to screen potential analogs of existing FDA approved drugs of MPXV against P37 using state-of-the-art machine learning and computational biophysical techniques. AlphaFold2 guided all-atoms molecular dynamics simulations optimized P37 structure is used for molecular docking and binding free energy calculations. Similar to members of Phospholipase-D family , the predicted P37 structure also adopts a β-α-β-α-β sandwich fold, harbouring strongly conserved HxKxxxxD motif. The binding pocket comprises of Tyr48, Lys86, His115, Lys117, Ser130, Asn132, Trp280, Asn240, His325, Lys327 and Tyr346 forming strong hydrogen bonds and dense hydrophobic contacts with the screened analogs and is surrounded by positively charged patches. Loops connecting the two domains and C-terminal region exhibit high degree of flexibility. In some structural ensembles, the partial disorderness in the C-terminal region is presumed to be due to its low confidence score, acquired during structure prediction. Transition from loop to β-strands (244-254 aa) in P37-Cidofovir and its analog complexes advocates the need for further investigations. MD simulations support the accuracy of the molecular docking results, indicating the potential of analogs as potent binders of P37. Taken together, our results provide preferable understanding of molecular recognition and dynamics of ligand-bound states of P37, offering opportunities for development of new antivirals against MPXV. However, the need of in vitro and in vivo assays for confirmation of these results still persists.Communicated by Ramaswamy H. Sarma.

Copper interaction with cystatin C: effects on protein structure and oligomerization

Human cystatin C (hCC), a small secretory protein, has gained attention beyond its classical role as a cysteine protease inhibitor owing to its potential involvement in neurodegenerative disorders. This study investigates the interaction between copper(II) ions [Cu(II)] and hCC, specifically targeting histidine residues known to participate in metal binding. Through various analytical techniques, including mutagenesis, circular dichroism, fluorescence assays, gel filtration chromatography, and electron microscopy, we evaluated the impact of Cu(II) ions on the structure and oligomerization of hCC. The results show that Cu(II) does not influence the secondary and tertiary structure of the studied hCC variants but affects their stability. To explore the Cu(II)-binding site, nuclear magnetic resonance (NMR) and X-ray studies were conducted. NMR experiments revealed notable changes in signal intensities and linewidths within the region 86His-Asp-Gln-Pro-His90, suggesting its involvement in Cu(II) coordination. Both histidine residues from this fragment were found to serve as a primary anchor of Cu(II) in solution, depending on the structural context and the presence of other Cu(II)-binding agents. The presence of Cu(II) led to significant destabilization and altered thermal stability of the wild-type and H90A variant, confirming differentiation between His residues in Cu(II) binding. In conclusion, this study provides valuable insights into the interaction between Cu(II) and hCC, elucidating the impact of copper ions on protein stability and identifying potential Cu(II)-binding residues. Understanding these interactions enhances our knowledge of the role of copper in neurodegenerative disorders and may facilitate the development of therapeutic strategies targeting copper-mediated processes in protein aggregation and associated pathologies.

Constitutive activity of ionotropic glutamate receptors via hydrophobic substitutions in the ligand-binding domain

Neurotransmitter ligands electrically excite neurons by activating ionotropic glutamate receptor (iGluR) ion channels. Knowledge of the iGluR amino acid residues that dominate ligand-induced activation would enable the prediction of function from sequence. We therefore explored the molecular determinants of activity in rat N-methyl-D-aspartate (NMDA)-type iGluRs (NMDA receptors), complex heteromeric iGluRs comprising two glycine-binding GluN1 and two glutamate-binding GluN2 subunits, using amino acid sequence analysis, mutagenesis, and electrophysiology. We find that a broadly conserved aspartate residue controls both ligand potency and channel activity, to the extent that certain substitutions at this position bypass the need for ligand binding in GluN1 subunits, generating NMDA receptors activated solely by glutamate. Furthermore, we identify a homomeric iGluR from the placozoan Trichoplax adhaerens that has utilized native mutations of this crucial residue to evolve into a leak channel that is inhibited by neurotransmitter binding, pointing to a dominant role of this residue throughout the iGluR superfamily.

Studying the Effects and Competitive Mechanisms of YOYO-1 on the Binding Characteristics of DOX and DNA Molecules Based on Surface-Enhanced Raman Spectroscopy and Molecular Docking Techniques

Revealing the interaction mechanisms between anticancer drugs and target DNA molecules at the single-molecule level is a hot research topic in the interdisciplinary fields of biophysical chemistry and pharmaceutical engineering. When fluorescence imaging technology is employed to carry out this kind of research, a knotty problem due to fluorescent dye molecules and drug molecules acting on a DNA molecule simultaneously is encountered. In this paper, based on self-made novel solid active substrates NpAA/(ZnO-ZnCl2)/AuNPs, we use a surface-enhanced Raman spectroscopy method, inverted fluorescence microscope technology, and a molecular docking method to investigate the action of the fluorescent dye YOYO-1 and the drug DOX on calf thymus DNA (ctDNA) molecules and the influencing effects and competitive relationships of YOYO-1 on the binding properties of the ctDNA-DOX complex. The interaction sites and modes of action between the YOYO-1 and the ctDNA-DOX complex are systematically examined, and the DOX with the ctDNA-YOYO-1 are compared, and the impact of YOYO-1 on the stability of the ctDNA-DOX complex and the competitive mechanism between DOX and YOYO-1 acting with DNA molecules are elucidated. This study has helpful experimental guidance and a theoretical foundation to expound the mechanism of interaction between drugs and biomolecules at the single-molecule level.

Comparative study of dietary phenols with Tartary buckwheat protein (2S/13S): impact on structure, binding sites and functionality of protein

BACKGROUND

This research was to investigate the interaction mechanism between 2S albumin and 13S globulin (2S and 13S, the most important storage proteins in Tartary buckwheat seeds) and three phenols (rutin, quercetin and myricetin) regarding the structural and antioxidant properties of their complexes.

RESULTS

There are differences in the binding affinity of phenols for 2S and 13S. Rutin had a higher binding affinity for 2S, myricetin had a higher binding affinity for 13S, and 13S exhibited a higher affinity toward phenols than did 2S. Binding with phenols significantly changed the secondary and tertiary structures of 2S and 13S, decreased the surface hydrophobic value and enhanced the antioxidant capacity. Molecular docking and isothermal titration calorimetry showed that the binding processes were spontaneous and that there were hydrogen bonds, hydrophobic bonds and van der Waals force interactions between phenols and proteins.

CONCLUSION

These findings could provide meaningful guidance for the further application of buckwheat protein complex. © 2023 Society of Chemical Industry.

Determination of the Cannabinoid CB1 Receptor's Positive Allosteric Modulator Binding Site through Mutagenesis Studies

Positive allosteric modulators (PAMs) of the cannabinoid CB1 receptor (CB1) offer potential therapeutic advantages in the treatment of neuropathic pain and addiction by avoiding the adverse effects associated with orthosteric CB1 activation. Here, molecular modeling and mutagenesis were used to identify residues central to PAM activity at CB1. Six putative allosteric binding sites were identified in silico, including novel sites previously associated with cholesterol binding, and key residues within each site were mutated to alanine. The recently determined ZCZ011 binding site was found to be essential for allosteric agonism, as GAT228, GAT229 and ZCZ011 all increased wild-type G protein dissociation in the absence of an orthosteric ligand; activity that was abolished in mutants F191A3.27 and I169A2.56. PAM activity was demonstrated for ZCZ011 in the presence of the orthosteric ligand CP55940, which was only abolished in I169A2.56. In contrast, the PAM activity of GAT229 was reduced for mutants R220A3.56, L404A8.50, F191A3.27 and I169A2.56. This indicates that allosteric modulation may represent the net effect of binding at multiple sites, and that allosteric agonism is likely to be mediated via the ZCZ011 site. This study underlines the need for detailed understanding of ligand receptor interactions in the search for pure CB1 allosteric modulators.

Water Exchange from the Buried Binding Sites of Cytochrome P450 Enzymes 1A2, 2D6, and 3A4 Correlates with Conformational Fluctuations

Human cytochrome P450 enzymes (CYPs) are critical for the metabolism of small-molecule pharmaceuticals (drugs). As such, the prediction of drug metabolism by and drug inhibition of CYP activity is an important component of the drug discovery and design process. Relative to the availability of a wide range of experimental atomic-resolution CYP structures, the development of structure-based CYP activity models has been limited. To better characterize the role of CYP conformational fluctuations in CYP activity, we perform multiple microsecond-scale all-atom explicit-solvent molecular dynamics (MD) simulations on three CYP isoforms, 1A2, 2D6, and 3A4, which together account for the majority of CYP-mediated drug metabolism. The MD simulations employ a variety of positional restraints, ranging from keeping all CYP atoms close to their experimentally determined coordinates to allowing full flexibility. We find that, with full flexibility, large fluctuations in the CYP binding sites correlate with efficient water exchange from these buried binding sites. This is especially true for 1A2, which, when restrained to its crystallographic conformation, is unable to exchange water between the binding site and bulk solvent. These findings imply that, in addition to crystal structures, a representative ensemble of conformational states ought to be included when developing structure-based CYP activity models.

Modelling and Molecular Dynamics Predict the Structure and Interactions of the Glycine Receptor Intracellular Domain

Glycine receptors (GlyRs) are glycine-gated inhibitory pentameric ligand-gated ion channels composed of α or α + β subunits. A number of structures of these proteins have been reported, but to date, these have only revealed details of the extracellular and transmembrane domains, with the intracellular domain (ICD) remaining uncharacterised due to its high flexibility. The ICD is a region that can modulate function in addition to being critical for receptor localisation and clustering via proteins such as gephyrin. Here, we use modelling and molecular dynamics (MD) to reveal details of the ICDs of both homomeric and heteromeric GlyR. At their N and C ends, both the α and β subunit ICDs have short helices, which are major sites of stabilising interactions; there is a large flexible loop between them capable of forming transient secondary structures. The α subunit can affect the β subunit ICD structure, which is more flexible in a 4α2:1β than in a 4α1:1β GlyR. We also explore the effects of gephyrin binding by creating GlyR models bound to the gephyrin E domain; MD simulations suggest these are more stable than the unbound forms, and again there are α subunit-dependent differences, despite the fact the gephyrin binds to the β subunit. The bound models also suggest that gephyrin causes compaction of the ICD. Overall, the data expand our knowledge of this important receptor protein and in particular clarify features of the underexplored ICD.

Spatial Layouts of Low-Entropy Hydration Shells Guide Protein Binding

Protein-protein binding enables orderly biological self-organization and is therefore considered a miracle of nature. Protein‒protein binding is driven by electrostatic forces, hydrogen bonding, van der Waals force, and hydrophobic interactions. Among these physical forces, only hydrophobic interactions can be considered long-range intermolecular attractions between proteins due to the electrostatic shielding of surrounding water molecules. Low-entropy hydration shells around proteins drive hydrophobic attraction among them that essentially coordinate protein‒protein binding. Here, an innovative method is developed for identifying low-entropy regions of hydration shells of proteins by screening off pseudohydrophilic groups on protein surfaces and revealing that large low-entropy regions of the hydration shells typically cover the binding sites of individual proteins. According to an analysis of determined protein complex structures, shape matching between a large low-entropy hydration shell region of a protein and that of its partner at the binding sites is revealed as a universal law. Protein‒protein binding is thus found to be mainly guided by hydrophobic collapse between the shape-matched low-entropy hydration shells that is verified by bioinformatics analyses of hundreds of structures of protein complexes, which cover four test systems. A simple algorithm is proposed to accurately predict protein binding sites.

Elucidation of the binding mode of organic polysulfides on the human TRPA1 receptor

Introduction: Previous studies have established that endogenous inorganic polysulfides have significant biological actions activating the Transient Receptor Potential Ankyrin 1 (TRPA1) receptor. Organic polysulfides exert similar effects, but they are much more stable molecules, therefore these compounds are more suitable as drugs. In this study, we aimed to better understand the mechanism of action of organic polysulfides by identification of their binding site on the TRPA1 receptor. Methods: Polysulfides can readily interact with the thiol side chain of the cysteine residues of the protein. To investigate their role in the TRPA1 activation, we replaced several cysteine residues by alanine via site-directed mutagenesis. We searched for TRPA1 mutant variants with decreased or lost activating effect of the polysulfides, but with other functions remaining intact (such as the effects of non-electrophilic agonists and antagonists). The binding properties of the mutant receptors were analyzed by in silico molecular docking. Functional changes were tested by in vitro methods: calcium sensitive fluorescent flow cytometry, whole-cell patch-clamp and radioactive calcium-45 liquid scintillation counting. Results: The cysteines forming the conventional binding site of electrophilic agonists, namely C621, C641 and C665 also bind the organic polysulfides, with the key role of C621. However, only their combined mutation abolished completely the organic polysulfide-induced activation of the receptor. Discussion: Since previous papers provided evidence that organic polysulfides exert analgesic and anti-inflammatory actions in different in vivo animal models, we anticipate that the development of TRPA1-targeted, organic polysulfide-based drugs will be promoted by this identification of the binding site.

Distribution and solvent exposure of Hsp70 chaperone binding sites across the Escherichia coli proteome

Many proteins must interact with molecular chaperones to achieve their native state in the cell. Yet, how chaperone binding-site characteristics affect the folding process is poorly understood. The ubiquitous Hsp70 chaperone system prevents client-protein aggregation by holding unfolded conformations and by unfolding misfolded states. Hsp70 binding sites of client proteins comprise a nonpolar core surrounded by positively charged residues. However, a detailed analysis of Hsp70 binding sites on a proteome-wide scale is still lacking. Further, it is not known whether proteins undergo some degree of folding while chaperone bound. Here, we begin to address the above questions by identifying Hsp70 binding sites in 2258 Escherichia coli (E. coli) proteins. We find that most proteins bear at least one Hsp70 binding site and that the number of Hsp70 binding sites is directly proportional to protein size. Aggregation propensity upon release from the ribosome correlates with number of Hsp70 binding sites only in the case of large proteins. Interestingly, Hsp70 binding sites are more solvent-exposed than other nonpolar sites, in protein native states. Our findings show that the majority of E. coli proteins are systematically enabled to interact with Hsp70 even if this interaction only takes place during a fraction of the protein lifetime. In addition, our data suggest that some conformational sampling may take place within Hsp70-bound states, due to the solvent exposure of some chaperone binding sites in native proteins. In all, we propose that Hsp70-chaperone-binding traits have evolved to favor Hsp70-assisted protein folding devoid of aggregation.

FcγRIIA-specific DARPins as novel tools in blood cell analysis and platelet aggregation

Fc receptors are involved in a variety of physiologically and disease relevant responses. Among them, FcγRIIA (CD32a) is known for its activating functions in pathogen recognition and platelet biology, and, as potential marker of T lymphocytes latently infected with HIV-1. The latter has not been without controversy due to technical challenges complicated by T-B cell conjugates and trogocytosis as well as a lack of antibodies distinguishing between the closely related isoforms of FcγRII. To generate high-affinity binders specific for FcγRIIA, libraries of designed ankyrin repeat proteins (DARPins) were screened for binding to its extracellular domains by ribosomal display. Counter-selection against FcγRIIB eliminated binders cross-reacting with both isoforms. The identified DARPins bound FcγRIIA with no detectable binding for FcγRIIB. Their affinities for FcγRIIA were in the low nanomolar range and could be enhanced by cleavage of the His-tag and dimerization. Interestingly, complex formation between DARPin and FcγRIIA followed a two-state reaction model and discrimination from FcγRIIB was based on a single amino acid residue. In flow cytometry, DARPin F11 detected FcγRIIA⁺ cells even when they made up less than 1% of the cell population. Image stream analysis of primary human blood cells confirmed that F11 caused dim but reliable cell surface staining of a small subpopulation of T lymphocytes. When incubated with platelets, F11 inhibited their aggregation equally efficient as antibodies unable to discriminate between both FcγRII isoforms. The selected DARPins are unique novel tools for platelet aggregation studies as well as the role of FcγRIIA for the latent HIV-1 reservoir.

Benchmarking biomolecular force field-based Zn2+ for mono- and bimetallic ligand binding sites

Zn²⁺ is one of the most versatile biologically available metal ions, but accurate modeling of Zn²⁺ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn²⁺ models are parameterized in bulk solvent, in-depth knowledge about their performance in a protein environment is limited. Thus, we systematically investigate here the behavior of non-polarizable Zn²⁺ models for their ability to reproduce experimentally determined metal coordination and ligand binding in metalloproteins. The benchmarking is performed in challenging environments, including mono- (carbonic anhydrase II) and bimetallic (metallo-β-lactamase VIM-2) ligand binding sites. We identify key differences in the performance between the Zn²⁺ models with regard to the preferred ligating atoms (charged/non-charged), attraction of water molecules, and the preferred coordination geometry. Based on these results, we suggest suitable simulation conditions for varying Zn²⁺ site geometries that could guide the further development of biomolecular Zn²⁺ models.

Druggable Pockets at the RNA Interface Region of Influenza A Virus NS1 Protein Are Conserved across Sequence Variants from Distinct Subtypes

Influenza A viruses still represent a major health issue, for both humans and animals. One of the main viral proteins of interest to target is the NS1 protein, which counters the host immune response and promotes viral replication. NS1 is a homodimer composed of a dimeric RNA-binding domain (RBD), which is structurally stable and conserved in sequence, and two effector domains that are tethered to the RBD by linker regions. This linker flexibility leads to NS1 polymorphism and can therefore exhibit different forms. Previously, we identified a putative drug-binding site, located in the RBD interface in a crystal structure of NS1. This pocket could be targeted to block RNA binding and inhibit NS1 activities. The objective of the present study is to confirm the presence of this druggable site, whatever the sequence variants, in order to develop a universal therapeutic compound that is insensitive to sequence variations and structural flexibility. Using a set of four NS1 full-length structures, we combined different bioinformatics approaches such as pocket tracking along molecular dynamics simulations, druggability prediction and classification. This protocol successfully confirmed a frequent large binding-site that is highly druggable and shared by different NS1 forms, which is promising for developing a robust NS1-targeted therapy.

Comprehensive Survey of Consensus Docking for High-Throughput Virtual Screening

The rapid advances of 3D techniques for the structural determination of proteins and the development of numerous computational methods and strategies have led to identifying highly active compounds in computer drug design. Molecular docking is a method widely used in high-throughput virtual screening campaigns to filter potential ligands targeted to proteins. A great variety of docking programs are currently available, which differ in the algorithms and approaches used to predict the binding mode and the affinity of the ligand. All programs heavily rely on scoring functions to accurately predict ligand binding affinity, and despite differences in performance, none of these docking programs is preferable to the others. To overcome this problem, consensus scoring methods improve the outcome of virtual screening by averaging the rank or score of individual molecules obtained from different docking programs. The successful application of consensus docking in high-throughput virtual screening highlights the need to optimize the predictive power of molecular docking methods.

Substrate Recognition Properties from an Intermediate Structural State of the UreA Transporter

Through a combination of comparative modeling, site-directed and classical random mutagenesis approaches, we previously identified critical residues for binding, recognition, and translocation of urea, and its inhibition by 2-thiourea and acetamide in the Aspergillus nidulans urea transporter, UreA. To deepen the structural characterization of UreA, we employed the artificial intelligence (AI) based AlphaFold2 (AF2) program. In this analysis, the resulting AF2 models lacked inward- and outward-facing cavities, suggesting a structural intermediate state of UreA. Moreover, the orientation of the W82, W84, N279, and T282 side chains showed a large variability, which in the case of W82 and W84, may operate as a gating mechanism in the ligand pathway. To test this hypothesis non-conservative and conservative substitutions of these amino acids were introduced, and binding and transport assessed for urea and its toxic analogue 2-thiourea, as well as binding of the structural analogue acetamide. As a result, residues W82, W84, N279, and T282 were implicated in substrate identification, selection, and translocation. Using molecular docking with Autodock Vina with flexible side chains, we corroborated the AF2 theoretical intermediate model, showing a remarkable correlation between docking scores and experimental affinities determined in wild-type and UreA mutants. The combination of AI-based modeling with classical docking, validated by comprehensive mutational analysis at the binding region, would suggest an unforeseen option to determine structural level details on a challenging family of proteins.

Deep learning for MYC binding site recognition

Motivation: The definition of the genome distribution of the Myc transcription factor is extremely important since it may help predict its transcriptional activity particularly in the context of cancer. Myc is among the most powerful oncogenes involved in the occurrence and development of more than 80% of different types of pediatric and adult cancers. Myc regulates thousands of genes which can be in part different, depending on the type of tissues and tumours. Myc distribution along the genome has been determined experimentally through chromatin immunoprecipitation This approach, although powerful, is very time consuming and cannot be routinely applied to tumours of individual patients. Thus, it becomes of paramount importance to develop in silico tools that can effectively and rapidly predict its distribution on a given cell genome. New advanced computational tools (DeeperBind) can then be successfully employed to determine the function of Myc in a specific tumour, and may help to devise new directions and approaches to experiments first and personalized and more effective therapeutic treatments for a single patient later on. Results: The use of DeeperBind with DeepRAM on Colab platform (Google) can effectively predict the binding sites for the MYC factor with an accuracy above 0.96 AUC, when trained with multiple cell lines. The analysis of the filters in DeeperBind trained models shows, besides the consensus sequence CACGTG classically associated to the MYC factor, also the other consensus sequences G/C box or TGGGA, respectively bound by the SP1 and MIZ-1 transcription factors, which are known to mediate the MYC repressive response. Overall, our findings suggest a stronger synergy between the machine learning tools as DeeperBind and biological experiments, which may reduce the time consuming experiments by providing a direction to guide them.

In Silico Methods for Identification of Potential Active Sites of Therapeutic Targets

Target identification is an important step in drug discovery, and computer-aided drug target identification methods are attracting more attention compared with traditional drug target identification methods, which are time-consuming and costly. Computer-aided drug target identification methods can greatly reduce the searching scope of experimental targets and associated costs by identifying the diseases-related targets and their binding sites and evaluating the druggability of the predicted active sites for clinical trials. In this review, we introduce the principles of computer-based active site identification methods, including the identification of binding sites and assessment of druggability. We provide some guidelines for selecting methods for the identification of binding sites and assessment of druggability. In addition, we list the databases and tools commonly used with these methods, present examples of individual and combined applications, and compare the methods and tools. Finally, we discuss the challenges and limitations of binding site identification and druggability assessment at the current stage and provide some recommendations and future perspectives.

Estimating the Similarity between Protein Pockets

With the exponential increase in publicly available protein structures, the comparison of protein binding sites naturally emerged as a scientific topic to explain observations or generate hypotheses for ligand design, notably to predict ligand selectivity for on- and off-targets, explain polypharmacology, and design target-focused libraries. The current review summarizes the state-of-the-art computational methods applied to pocket detection and comparison as well as structural druggability estimates. The major strengths and weaknesses of current pocket descriptors, alignment methods, and similarity search algorithms are presented. Lastly, an exhaustive survey of both retrospective and prospective applications in diverse medicinal chemistry scenarios illustrates the capability of the existing methods and the hurdle that still needs to be overcome for more accurate predictions.

On the Mechanism of Selective Chemical Exchange of Bacteriopheophytins in the Reaction Centers of Rhodobacter sphaeroides R-26

To elucidate the mechanism of site-selective chemical replacement of chromophores in the reaction centers (RCs) of photosynthetic bacteria by external pigments, we investigated how the efficiency of incorporation of plant pheophytin a (Pheo) into the binding sites for bacteriopheophytin a molecules (BPheo) in the isolated Rhodobacter sphaeroides R-26 RCs depended on the incubation medium temperature, Pheo aggregation state, and the presence of organic solvent (acetone). When Pheo was in a form of monomers in free detergent micelles in a water-detergent incubation medium, the degree of selective replacement of photochemically inactive BPheo H_B molecules upon incubation of the RC/Pheo mixture at 5°C was ~15%. The exchange efficiency increased to 40% upon incubation at 25°C and reached 100% at the same temperature when 10% acetone was added to the incubation medium. At both 5 and 25°C, the degree of pigment exchange increased approximately twice, when a mixture of Pheo monomers and dimers in the presence of 10% acetone was used as the incubation medium. The removal of acetone from this medium with the preservation of pigment forms led to a significant decrease in the efficiency of Pheo incorporation. The effect of acetone on the pigment exchange was also observed at an elevated incubation temperature (43.5°C), when functionally active BPheo H_A molecules were partially replaced. The results are discussed in terms of the mechanism according to which (i) the temperature-dependent internal movements of the RC protein facilitate the release of the BPheo molecule from the binding site with simultaneous insertion of the Pheo molecule into the same site in a coupled process, (ii) the role of temperature largely depends on the steric accessibility of binding pockets in the RC protein, (iii) the incorporation of Pheo occurs from a pool of monomeric molecules included in the RC-detergent micelles, and (iv) the presence of acetone in the incubation medium facilitates the exchange of Pheo monomers between micelles in the solution and the detergent belt of the RC complex.

Profiling of lysophosphatidylethanolamine molecular species in human serum and in silico prediction of the binding site on albumin

Lysophosphatidylethanolamine (LPE) is a major lysophospholipid produced by phospholipids and binds to human serum albumin (HSA). LPEs may play various roles in vivo depending on the differences in their acyl chains. However, only few reports have been published on the biological functions of LPEs. Hence, we determined the exact relative abundance of the major LPEs in the serum of healthy participants (n = 8) using liquid chromatography-tandem mass spectrometry. Consequently, LPE 18:2 (24.1 ± 5.2%) was found to be the most abundant in serum. To understand the distribution of LPEs, the serum separated via gel-filtration high-performance liquid chromatography was subjected to quantitative measurement. LPEs were more observed in the albumin fraction than the lipoprotein fraction. We also performed a fluorescence displacement assay and an in silico molecular docking experiment using AutoDock to confirm the affinity and binding sites of the LPEs on HSA. The binding affinities of the LPEs for drug sites 1 and 2 on HSA were relatively low, with Ki values of approximately 11 and 3.8 μM, respectively. AutoDock analysis revealed the conformation of the LPEs bound to drug sites and the possibility of LPEs binding to other HSA sites. These findings could help to elucidate the biological and pathological functions of LPEs.

5-HT3 Receptor MX Helix Contributes to Receptor Function

5-HT₃ receptors are members of the family of pentameric ligand-gated ion channels. Each subunit has an extracellular, transmembrane, and intracellular domain. Only part of the intracellular domain structure has been solved, revealing it contains two α-helical segments; one, the MA helix, is an extension of M4, while the other, the MX helix, is formed from residues located close to the end of M3. This MX helix is in distinct locations in open and closed receptor structures, suggesting it may play a role in function. Here, we explore this hypothesis using functional responses of Ala-substituted mutant receptors expressed in HEK293 cells. The data show altering many of the MX residues results in a small decrease in EC₅₀ (up to 5-fold), although in one (H232A) this is increased. Radiolabeled ligand binding on selected mutants showed no change in binding affinity, indicating an effect on gating and not binding. In addition, five mutations (P316A, V317A, P318A, D319A, and H323A) initially resulted in nonfunctional receptors, but the function could be rescued by coexpression with a chaperone protein, suggesting a likely role in assembly or folding. Examination of previously obtained MD simulation data shows that the extent of MX encompassed by membrane lipids differs considerably in the open and closed structures, suggesting that lipid–protein interactions in this region could have a major effect on channel opening propensity. We conclude that the MX helix can modulate the function of the receptor and propose that its interactions with membrane lipids play a major role in this.

Adsorption and Diffusion of Oxygen on Pure and Partially Oxidized Metal Surfaces in Ultrahigh Resolution

The interaction of gas molecules with metal and oxide surfaces plays a critical role in corrosion, catalysis, sensing, and heterogeneous materials. However, insights into the dynamics of O₂ from picoseconds to microseconds have remained unavailable to date. We obtained 3D potential energy surfaces for adsorption of O₂ on 11 common pristine and partially oxidized (hkl) surfaces of Ni and Al in picometer resolution and high accuracy of 0.1 kcal/mol, identified binding sites, and surface mobility from 25 to 300 °C. We explain relative oxidation rates and parameters for oxide growth. We employed over 150 000 molecular mechanics and molecular dynamics simulations with the interface force field (IFF) using structural data from X-ray diffraction (XRD) and low-energy electron diffraction (LEED). The methods reach 10 to 50 times higher accuracy than possible before and are suited to analyze gas interactions with metals up to the micrometer scale including defects and irregular nanostructures.

GSH facilitates the binding and inhibitory activity of novel multidrug resistance protein 1 (MRP1) modulators

MRP1 (ABCC1) is a membrane transporter that confers multidrug resistance in cancer cells by exporting chemotherapeutic agents, often in a reduced glutathione (GSH)-dependent manner. This transport activity can be altered by compounds (modulators) that block drug transport while simultaneously stimulating GSH efflux by MRP1. In MRP1-expressing cells, modulator-stimulated GSH efflux can be sufficient to deplete GSH and increase sensitivity to chemotherapy, enhancing cancer cell death. Further development of clinically useful MRP1 modulators requires a better mechanistic understanding of modulator binding and its relationship to GSH binding and transport. Here, we explore the mechanism of two MRP1 small molecule modulators, 5681014 and 7914321, in relation to a bipartite substrate-binding cavity of MRP1. Binding of these modulators to MRP1 was dependent on the presence of GSH but not its reducing capacity. Accordingly, the modulators poorly inhibited organic anion transport by K332L-mutant MRP1, where GSH binding and transport is limited. However, the inhibitory activity of the modulators was also diminished by mutations that limit E₂ 17βG but spare GSH-conjugate binding and transport (W553A, M1093A, W1246A), suggesting overlap between the E₂ 17βG and modulator binding sites. Immunoblots of limited trypsin digests of MRP1 suggest that binding of GSH, but not the modulators, induces a conformation change in MRP1. Together, these findings support the model, in which GSH binding induces a conformation change that facilitates binding of MRP1 modulators, possibly in a proposed hydrophobic binding pocket of MRP1. This study may facilitate the structure-guided design of more potent and selective MRP1 modulators.

The Mechanism of Dynamic Interaction between Doxorubicin and Calf Thymus DNA at the Single-Molecule Level Based on Confocal Raman Spectroscopy

It is of great fundamental significance and practical application to understand the binding sites and dynamic process of the interaction between doxorubicin (DOX) and DNA molecules. Based on the Confocal Raman spectroscopy, the interaction between DOX and calf thymus DNA has been systemically investigated, and some meaningful findings have been found. DOX molecules can not only interact with all four bases of DNA molecules, i.e., adenine, thymine, cytosine, guanine, and phosphate, but also affect the DNA conformation. Meanwhile, the binding site of DOX and its derivatives such as daunorubicin and epirubicin is certain. Furthermore, the interaction between DOX and DNA molecules is a dynamic process since the intensities of each characteristic peaks of the base, e.g., adenine, cytosine, and phosphate, are all regularly changed with the interaction time. Finally, a dynamic mechanism model of the interaction between DOX and DNA molecules is proposed; that is, there are two kinds of interaction between DOX and DNA molecules: DOX-DNA acts to form a complex, and DOX-DOX acts to form a multimer. The two effects are competitive, as the former compresses DNA molecules, and the latter decompresses these DNA molecules. This work is helpful for accurately understanding and developing new drugs and pathways to improve and treat DOX-induced cytotoxicity and cardiotoxicity.

Cryo-electron microscopy reveals how acetogenins inhibit mitochondrial respiratory complex I

Mitochondrial complex I (NADH:ubiquinone oxidoreductase), a crucial enzyme in energy metabolism, captures the redox potential energy from NADH oxidation/ubiquinone reduction to create the proton motive force used to drive ATP synthesis in oxidative phosphorylation. High-resolution single-particle electron cryo-EM analyses have provided detailed structural knowledge of the catalytic machinery of complex I, but not of the molecular principles of its energy transduction mechanism. Although ubiquinone is considered to bind in a long channel at the interface of the membrane-embedded and hydrophilic domains, with channel residues likely involved in coupling substrate reduction to proton translocation, no structures with the channel fully occupied have yet been described. Here, we report the structure (determined by cryo-EM) of mouse complex I with a tight-binding natural product acetogenin inhibitor, which resembles the native substrate, bound along the full length of the expected ubiquinone-binding channel. Our structure reveals the mode of acetogenin binding and the molecular basis for structure-activity relationships within the acetogenin family. It also shows that acetogenins are such potent inhibitors because they are highly hydrophobic molecules that contain two specific hydrophilic moieties spaced to lock into two hydrophilic regions of the otherwise hydrophobic channel. The central hydrophilic section of the channel does not favor binding of the isoprenoid chain when the native substrate is fully bound but stabilizes the ubiquinone/ubiquinol headgroup as it transits to/from the active site. Therefore, the amphipathic nature of the channel supports both tight binding of the amphipathic inhibitor and rapid exchange of the ubiquinone/ubiquinol substrate and product.

Effect of Affinity between Dechlorane Plus and Human Serum Albumin on Its Serum Concentration in a General Population

Dechlorane Plus (DP) has been detected in human serum around the world. Determining how DP binds to human serum albumin (HSA) could improve our understanding of the effects of DP in humans. The results of a fluorescence titration experiment indicated that DP binds to HSA at the interface between domains IIA and IIB, and HSA has more affinity for anti-DP than syn-DP. The DP concentrations in serum from 33 people living in Wuxi City, China, were determined in five consecutive years (2012-2016). The total concentrations of the two DP isomers were 1.73-16.3 ng·g^-1 lipid weight (mean 4.7 ng·g^-1 lipid weight). No clear increasing or decreasing temporal trend was found in the DP concentrations in serum over the study period. The affinity between DP and HSA will increase as the amount of fatty acids combined with HSA increases. This explained the syn-DP and anti-DP concentrations in serum being higher for overweight people than for normal weight people. Ethanol can occupy the DP binding sites on HSA. This explained the syn-DP and anti-DP concentrations being lower for alcoholic people than for nonalcoholic people.

Artificial intelligence in the prediction of protein-ligand interactions: recent advances and future directions

New drug production, from target identification to marketing approval, takes over 12 years and can cost around $2.6 billion. Furthermore, the COVID-19 pandemic has unveiled the urgent need for more powerful computational methods for drug discovery. Here, we review the computational approaches to predicting protein-ligand interactions in the context of drug discovery, focusing on methods using artificial intelligence (AI). We begin with a brief introduction to proteins (targets), ligands (e.g. drugs) and their interactions for nonexperts. Next, we review databases that are commonly used in the domain of protein-ligand interactions. Finally, we survey and analyze the machine learning (ML) approaches implemented to predict protein-ligand binding sites, ligand-binding affinity and binding pose (conformation) including both classical ML algorithms and recent deep learning methods. After exploring the correlation between these three aspects of protein-ligand interaction, it has been proposed that they should be studied in unison. We anticipate that our review will aid exploration and development of more accurate ML-based prediction strategies for studying protein-ligand interactions.

Mass Spectrometry-Based Structural Proteomics for Metal Ion/Protein Binding Studies

Metal ions are critical for the biological and physiological functions of many proteins. Mass spectrometry (MS)-based structural proteomics is an ever-growing field that has been adopted to study protein and metal ion interactions. Native MS offers information on metal binding and its stoichiometry. Footprinting approaches coupled with MS, including hydrogen/deuterium exchange (HDX), "fast photochemical oxidation of proteins" (FPOP) and targeted amino-acid labeling, identify binding sites and regions undergoing conformational changes. MS-based titration methods, including "protein-ligand interactions by mass spectrometry, titration and HD exchange" (PLIMSTEX) and "ligand titration, fast photochemical oxidation of proteins and mass spectrometry" (LITPOMS), afford binding stoichiometry, binding affinity, and binding order. These MS-based structural proteomics approaches, their applications to answer questions regarding metal ion protein interactions, their limitations, and recent and potential improvements are discussed here. This review serves as a demonstration of the capabilities of these tools and as an introduction to wider applications to solve other questions.

Chimeric Investigations into the Diamide Binding Site on the Lepidopteran Ryanodine Receptor

Alterations to amino acid residues G4946 and I4790, associated with resistance to diamide insecticides, suggests a location of diamide interaction within the pVSD voltage sensor-like domain of the insect ryanodine receptor (RyR). To further delineate the interaction site(s), targeted alterations were made within the same pVSD region on the diamondback moth (Plutella xylostella) RyR channel. The editing of five amino acid positions to match those found in the diamide insensitive skeletal RyR1 of humans (hRyR1) in order to generate a human–Plutella chimeric construct showed that these alterations strongly reduce diamide efficacy when introduced in combination but cause only minor reductions when introduced individually. It is concluded that the sites of diamide interaction on insect RyRs lie proximal to the voltage sensor-like domain of the RyR and that the main site of interaction is at residues K4700, Y4701, I4790 and S4919 in the S1 to S4 transmembrane domains.

A {Ni12}-Wheel-Based Metal-Organic Framework for Coordinative Binding of Sulphur Dioxide and Nitrogen Dioxide

Air pollutions by SO 2 and NO 2 have caused significant risks on the environment and human health. Understanding the mechanism of active sites within capture materials is of fundamental importance to the development of new clean-up technologies. Here we report the crystallographic observation of reversible coordinative binding of SO 2 and NO 2 on open Ni(II) sites in a metal-organic framework (NKU-100) incorporating an unprecedented {Ni 12 }-wheel, which exhibits six open Ni(II) sites on desolvation. Immobilised gas molecules are further stabilised by cooperative host-guest interactions comprised of hydrogen bonds, π ··· π interactions and dipole interactions. At 298 K and 1.0 bar, NKU-100 shows adsorption uptakes of 6.21 and 5.80 mmol g -1 for SO 2 and NO 2 , respectively. Dynamic breakthrough experiments have confirmed the selective retention of SO 2 and NO 2 at low concentrations under dry conditions. This work will inspire the future design of efficient sorbents for the capture of SO 2 and NO 2 .

The Candida glabrata glycogen branching enzyme structure reveals unique features of branching enzymes of the Saccharomycetaceae phylum

Branching enzymes (BE) are responsible for the formation of branching points at the 1,6 position in glycogen and starch, by catalyzing the cleavage of α-1,4-linkages and the subsequent transfer by introducing α-1,6-linked glucose branched points. BEs are found in the large GH13 family, eukaryotic BEs being mainly classified in the GH13_8 subfamily, GH13_9 grouping almost exclusively prokaryotic enzymes. With the aim of contributing to the understanding of the mode of recognition and action of the enzymes belonging to GH13_8, and to the understanding of features distinguishing these enzymes from those belonging to subfamily 13_9 we solved the crystal structure of the glycogen branching enzyme (GBE) from the yeast Candida glabrata, CgGBE, in ligand free forms and in complex with a maltotriose. The structures revealed the presence of a domain already observed in Homo sapiens and Oryza sativa BEs and that we named α-helical N-terminal domain, in addition to the three conserved domains found in BE. We confirmed by phylogenetic analysis that this α-helical N-terminal domain is always present in the GH13_8 enzymes suggesting that it could actually present a signature for this subfamily. We identified two binding sites (BS) in the α-helical N-terminal domain and in the carbohydrate binding module 48 (CBM48), respectively, which show a unique structural organization only present in the Saccharomycotina phylum. Our structural and phylogenetic investigation provides new insight into the structural characterization of GH13_8 GBE revealing unique structural features only present in the Saccharomycotina phylum thereby conferring original properties to this group of enzymes.

Effect of Tetraphenylborate on Physicochemical Properties of Bovine Serum Albumin

The binding interactions of bovine serum albumin (BSA) with tetraphenylborate ions ([B(Ph)4]-) have been investigated by a set of experimental methods (isothermal titration calorimetry, steady-state fluorescence spectroscopy, differential scanning calorimetry and circular dichroism spectroscopy) and molecular dynamics-based computational approaches. Two sets of structurally distinctive binding sites in BSA were found under the experimental conditions (10 mM cacodylate buffer, pH 7, 298.15 K). The obtained results, supported by the competitive interactions experiments of SDS with [B(Ph)4]- for BSA, enabled us to find the potential binding sites in BSA. The first site is located in the subdomain I A of the protein and binds two [B(Ph)4]- ions (logK(ITC)1 = 7.09 ± 0.10; ΔG(ITC)1 = -9.67 ± 0.14 kcal mol-1; ΔH(ITC)1 = -3.14 ± 0.12 kcal mol-1; TΔS(ITC)1 = -6.53 kcal mol-1), whereas the second site is localized in the subdomain III A and binds five ions (logK(ITC)2 = 5.39 ± 0.06; ΔG(ITC)2 = -7.35 ± 0.09 kcal mol-1; ΔH(ITC)2 = 4.00 ± 0.14 kcal mol-1; TΔS(ITC)2 = 11.3 kcal mol-1). The formation of the {[B(Ph)4]-}-BSA complex results in an increase in the thermal stability of the alfa-helical content, correlating with the saturation of the particular BSA binding sites, thus hindering its thermal unfolding.

Accurate Sequence-Based Prediction of Deleterious nsSNPs with Multiple Sequence Profiles and Putative Binding Residues

Non-synonymous single nucleotide polymorphisms (nsSNPs) may result in pathogenic changes that are associated with human diseases. Accurate prediction of these deleterious nsSNPs is in high demand. The existing predictors of deleterious nsSNPs secure modest levels of predictive performance, leaving room for improvements. We propose a new sequence-based predictor, DMBS, which addresses the need to improve the predictive quality. The design of DMBS relies on the observation that the deleterious mutations are likely to occur at the highly conserved and functionally important positions in the protein sequence. Correspondingly, we introduce two innovative components. First, we improve the estimates of the conservation computed from the multiple sequence profiles based on two complementary databases and two complementary alignment algorithms. Second, we utilize putative annotations of functional/binding residues produced by two state-of-the-art sequence-based methods. These inputs are processed by a random forests model that provides favorable predictive performance when empirically compared against five other machine-learning algorithms. Empirical results on four benchmark datasets reveal that DMBS achieves AUC > 0.94, outperforming current methods, including protein structure-based approaches. In particular, DMBS secures AUC = 0.97 for the SNPdbe and ExoVar datasets, compared to AUC = 0.70 and 0.88, respectively, that were obtained by the best available methods. Further tests on the independent HumVar dataset shows that our method significantly outperforms the state-of-the-art method SNPdryad. We conclude that DMBS provides accurate predictions that can effectively guide wet-lab experiments in a high-throughput manner.

miR-496 inhibits proliferation via LYN and AKT pathway in gastric cancer

MicroRNAs (miRNAs) operate as tumor suppressor or carcinogen to regulate cell proliferation, metastasis, invasion, differentiation, apoptosis, and metabolic process. In the present research, we investigated the effect and mechanism of miR-496 in human gastric cancer cells. miR-496 was downregulated in two gastric cancer cell lines, AGS and MKN45, compared with normal gastric epithelial cell line GES-1. miR-496 mimics inhibited the proliferation of AGS cells after the transfection for 48 and 72 h. The migration and invasion of AGS cells were also inhibited by the transfection of miR-496 mimics. miR-496 mimics induced the apoptosis through upregulating the levels of Bax and Active Caspase 3 and downregulating the levels of Bcl-2 and Total Caspase 3. Bioinformatics analysis showed that there was a binding site between miR-496 and Lyn kinase (LYN). miR-496 mimics could inhibit the expression of LYN in AGS cells. LYN overexpression blocked the inhibition of tumor cell growth, as well as the inhibition of AKT/mTOR signaling pathway induced by miR-496. In conclusion, miR-496 inhibited the proliferation through the AKT/mTOR signaling pathway via targeting LYN in gastric cancer cells. Our research provides a new potential target for clinical diagnosis and targeted treatment for gastric cancer.

A study of a hierarchical structure of proteins and ligand binding sites of receptors using the triangular spatial relationship-based structure comparison method and development of a size-filtering feature designed for comparing different sizes of protein structures

The presence of receptors and the specific binding of the ligands determine nearly all cellular responses. Binding of a ligand to its receptor causes conformational changes of the receptor that triggers the subsequent signaling cascade. Therefore, systematically studying structures of receptors will provide insight into their functions. We have developed the triangular spatial relationship (TSR)-based method where all possible triangles are constructed with C_α atoms of a protein as vertices. Every triangle is represented by an integer denoted as a "key" computed through the TSR algorithm. A structure is thereby represented by a vector of integers. In this study, we have first defined substructures using different types of keys. Second, using different types of keys represents a new way to interpret structure hierarchical relations and differences between structures and sequences. Third, we demonstrate the effects of sequence similarity as well as sample size on the structure-based classifications. Fourth, we show identification of structure motifs, and the motifs containing multiple triangles connected by either an edge or a vertex are mapped to the ligand binding sites of the receptors. The structure motifs are valuable resources for the researchers in the field of signal transduction. Next, we propose amino-acid scoring matrices that capture "evolutionary closeness" information based on BLOSUM62 matrix, and present the development of a new visualization method where keys are organized according to evolutionary closeness and shown in a 2D image. This new visualization opens a window for developing tools with the aim of identification of specific and common substructures by scanning pixels and neighboring pixels. Finally, we report a new algorithm called as size filtering that is designed to improve structure comparison of large proteins with small proteins. Collectively, we provide an in-depth interpretation of structure relations through the detailed analyses of different types of keys and their associated key occurrence frequencies, geometries, and labels. In summary, we consider this study as a new computational platform where keys are served as a bridge to connect sequence and structure as well as structure and function for a deep understanding of sequence, structure, and function relationships of the protein family.

Recurrent Loss of abaA, a Master Regulator of Asexual Development in Filamentous Fungi, Correlates with Changes in Genomic and Morphological Traits

Gene regulatory networks (GRNs) drive developmental and cellular differentiation, and variation in their architectures gives rise to morphological diversity. Pioneering studies in Aspergillus fungi, coupled with subsequent work in other filamentous fungi, have shown that the GRN governed by the BrlA, AbaA, and WetA proteins controls the development of the asexual fruiting body or conidiophore. A specific aspect of conidiophore development is the production of phialides, conidiophore structures that are under the developmental control of AbaA and function to repetitively generate spores. Fungal genome sequencing has revealed that some filamentous fungi lack abaA, and also produce asexual structures that lack phialides, raising the hypothesis that abaA loss is functionally linked to diversity in asexual fruiting body morphology. To examine this hypothesis, we carried out an extensive search for the abaA gene across 241 genomes of species from the fungal subphylum Pezizomycotina. We found that abaA was independently lost in four lineages of Eurotiomycetes, including from all sequenced species within the order Onygenales, and that all four lineages that have lost abaA also lack the ability to form phialides. Genetic restoration of abaA from Aspergillus nidulans into Histoplasma capsulatum, a pathogenic species from the order Onygenales that lacks an endogenous copy of abaA, did not alter Histoplasma conidiation morphology but resulted in a marked increase in spore viability. We also discovered that species lacking abaA contain fewer AbaA binding motifs in the regulatory regions of orthologs of some AbaA target genes, suggesting that the asexual fruiting body GRN of organisms that have lost abaA has likely been rewired. Our results provide an illustration of how repeated losses of a key regulatory transcription factor have contributed to the diversity of an iconic fungal morphological trait.

Exploring molecular determinants of polysaccharide Lyase family 6-1 enzyme activity

The Polysaccharide Lyase Family 6 (PL6) represents one of the 41 polysaccharide lyase families classified in the CAZy database with the vast majority of its members being alginate lyases grouped into three subfamilies, PL6_1-3. To decipher the mode of recognition and action of the enzymes belonging to subfamily PL6_1, we solved the crystal structures of Pedsa0632, Patl3640, Pedsa3628 and Pedsa3807, which all show different substrate specificities and mode of action (endo-/exo-lyase). Thorough exploration of the structures of Pedsa0632 and Patl3640 in complex with their substrates as well as docking experiments confirm that the conserved residues in subsites -1 to +3 of the catalytic site form a common platform which can accommodate various types of alginate in a very similar manner but with a series of original adaptations bringing them their specificities of action. From comparative studies with existing structures of PL6_1 alginate lyases, we observe that in the right-handed parallel β-helix fold shared by all these enzymes, the substrate binding site harbors the same overall conserved structures and organization. Despite this apparent similarity, it appears that members of the PL6_1 subfamily specifically accommodate and catalyze the degradation of different alginates suggesting that this common platform is actually a highly adaptable and specific tool.

ATAC-Seq Data for Genome-Wide Profiling of Transcription Factor Binding Sites in the Rice False Smut Fungus Ustilaginoidea virens

Identification of transcription factor binding sites is one of the most important steps in understanding the function of transcription factors and regulatory networks in organisms. The assay for transposase accessible chromatin sequencing (ATAC-seq) is a simple protocol for detection of open chromatin that could be a powerful tool to advance studies of protein-DNA interactions. Although ATAC-seq has been used in systematic identification of cis-regulatory regions in animal and plant genomes, this method has been rarely applied in fungi. Here, we describe a valuable ATAC-seq resource in the genome of an economically important phytopathogen, the rice false smut fungus Ustilaginoidea virens. The ATAC-seq data of U. virens mycelia collected from potato sucrose broth (PSB) and PSB supplied with rice spikelet extract were both generated. This is the first genome-wide profiling of open chromatin and transcription factor binding sites in U. virens.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Recognition of Ion Ligand Binding Sites Based on Amino Acid Features with the Fusion of Energy, Physicochemical and Structural Features

[Background: Rational drug molecular design based on virtual screening requires the ligand binding site to be known. Recently, the recognition of ion ligand binding site has become an important research direction in pharmacology.

METHODS

In this work, we selected the binding residues of 4 acid radical ion ligands (NO₂ ^-, CO₃ ^2-, SO₄ ^2- and PO₄ ^3-) and 10 metal ion ligands (Zn²⁺, Cu²⁺, Fe²⁺, Fe³⁺, Ca²⁺, Mg²⁺, Mn²⁺, Na⁺, K⁺ and Co²⁺) as research objects. Based on the protein sequence information, we extracted amino acid features, energy, physicochemical, and structure features. Then, we incorporated the above features and input them into the MultilayerPerceptron (MLP) and support vector machine (SVM) algorithms.

RESULTS

In the independent test, the best accuracy was higher than 92.5%, which was better than the previous results on the same dataset. In addition, we found that energy information is an important factor affecting the prediction results.

CONCLUSION

Finally, we set up a free web server for the prediction of protein-ion ligand binding sites (http://39.104.77.103:8081/lsb/HomePage/HomePage.html). This study is helpful for molecular drug design.

Discovering unknown human and mouse transcription factor binding sites and their characteristics from ChIP-seq data

Transcription factor binding sites (TFBSs) are essential for gene regulation, but the number of known TFBSs remains limited. We aimed to discover and characterize unknown TFBSs by developing a computational pipeline for analyzing ChIP-seq (chromatin immunoprecipitation followed by sequencing) data. Applying it to the latest ENCODE ChIP-seq data for human and mouse, we found that using the irreproducible discovery rate as a quality-control criterion resulted in many experiments being unnecessarily discarded. By contrast, the number of motif occurrences in ChIP-seq peak regions provides a highly effective criterion, which is reliable even if supported by only one experimental replicate. In total, we obtained 2,058 motifs from 1,089 experiments for 354 human TFs and 163 motifs from 101 experiments for 34 mouse TFs. Among these motifs, 487 have not previously been reported. Mapping the canonical motifs to the human genome reveals a high TFBS density ±2 kb around transcription start sites (TSSs) with a peak at -50 bp. On average, a promoter contains 5.7 TFBSs. However, 70% of TFBSs are in introns (41%) and intergenic regions (29%), whereas only 12% are in promoters (-1 kb to +100 bp from TSSs). Notably, some TFs (e.g., CTCF, JUN, JUNB, and NFE2) have motifs enriched in intergenic regions, including enhancers. We inferred 142 cobinding TF pairs and 186 (including 115 completely) tethered binding TF pairs, indicating frequent interactions between TFs and a higher frequency of tethered binding than cobinding. This study provides a large number of previously undocumented motifs and insights into the biological and genomic features of TFBSs.

ELAVL1 primarily couples mRNA stability with the 3' UTRs of interferon-stimulated genes

Upon pathogen detection, the innate immune system triggers signaling events leading to upregulation of pro-inflammatory and anti-microbial mRNA transcripts. RNA-binding proteins (RBPs) interact with these critical mRNAs and regulate their fates at the post-transcriptional level. One such RBP is ELAVL1. Although significant progress has been made in understanding how embryonic lethal vision-like protein 1 (ELAVL1) regulates mRNAs, its target repertoire and binding distribution within an immunological context remain poorly understood. We overlap four high-throughput approaches to define its context-dependent targets and determine its regulatory impact during immune activation. ELAVL1 transitions from binding overwhelmingly intronic sites to 3′ UTR sites upon immune stimulation of cells, binding previously and newly expressed mRNAs. We find that ELAVL1 mediates the RNA stability of genes that regulate pathways essential to pathogen sensing and cytokine production. Our findings reveal the importance of examining RBP regulatory impact under dynamic transcriptomic events to understand their post-transcriptional regulatory roles within specific biological circuitries.

Rothamel et al. show that upon immune activation, the RNA-binding protein ELAVL1 accumulates in the cytoplasm and redistributes from introns to mRNA 3′ UTRs. 3′ UTR binding confers enrichment and transcript stability. Many top-ranking transcripts are interferon-stimulated genes (ISGs), indicating that ELAVL1 is a positive regulator of an innate immune response.

A tri-functional amino acid enables mapping of binding sites for posttranslational-modification-mediated protein-protein interactions

Posttranslational modification (PTM), through the recruitment of effector proteins (i.e., "readers") that signal downstream events, plays key roles in regulating a variety of cellular processes. To understand how a PTM is recognized, it is necessary to find its readers and, importantly, the location of the binding pockets responsible for PTM recognition. Although various methods have been developed to identify PTM readers, it remains a challenge to directly map the PTM-binding regions, especially for intrinsically disordered domains. Here, we demonstrate a photo-crosslinkable, clickable, and cleavable tri-functional amino acid, ADdis-Cys, that when coupled with mass spectrometry (ADdis-Cys-MS) can not only identify PTM readers from complex proteomes but also simultaneously map their PTM-recognition modules. Using ADdis-Cys-MS, we successfully identify the binding sites of several reader-PTM interactions, among which we discover human C1QBP as a histone chaperone. This robust method should find wide applications in examining other histone or non-histone PTM-mediated protein-protein interactions.

Isolation and Functional Characterization of the Promoters of Miltiradiene Synthase Genes, TwTPS27a and TwTPS27b, and Interaction Analysis with the Transcription Factor TwTGA1 from Tripterygium wilfordii

Miltiradiene synthase (MS) genes, TwTPS27a and TwTPS27b, are the key diterpene synthase genes in the biosynthesis of triptolide, which is an important medicinally active diterpenoid in Tripterygium wilfordii. However, the mechanism underlying the regulation of key genes TwTPS27a/b in triptolide biosynthesis remains unclear. In this study, the promoters of TwTPS27a (1496 bp) and TwTPS27b (1862 bp) were isolated and analyzed. Some hormone-/stress-responsive elements and transcription factor (TF) binding sites were predicted in both promoters, which might be responsible for the regulation mechanism of TwTPS27a/b. The β-glucuronidase (GUS) activity analysis in promoter deletion assays under normal and methyl jasmonate (MeJA) conditions showed that the sequence of −921 to −391 bp is the potential core region of the TwTPS27b promoter. And the TGACG-motif, a MeJA-responsive element found in this core region, might be responsible for MeJA-mediated stress induction of GUS activity. Moreover, the TGACG-motif is also known as the TGA TF-binding site. Yeast one-hybrid and GUS transactivation assays confirmed the interaction between the TwTPS27a/b promoters and the TwTGA1 TF (a MeJA-inducible TGA TF upregulating triptolide biosynthesis in T. wilfordii), indicating that TwTPS27a/b are two target genes regulated by TwTGA1. In conclusion, our results provide important information for elucidating the regulatory mechanism of MS genes, TwTPS27a and TwTPS27b, as two target genes of TwTGA1, in jasmonic acid (JA)-inducible triptolide biosynthesis.

Study of Endogen Substrates, Drug Substrates and Inhibitors Binding Conformations on MRP4 and Its Variants by Molecular Docking and Molecular Dynamics

Multidrug resistance protein-4 (MRP4) belongs to the ABC transporter superfamily and promotes the transport of xenobiotics including drugs. A non-synonymous single nucleotide polymorphisms (nsSNPs) in the ABCC4 gene can promote changes in the structure and function of MRP4. In this work, the interaction of certain endogen substrates, drug substrates, and inhibitors with wild type-MRP4 (WT-MRP4) and its variants G187W and Y556C were studied to determine differences in the intermolecular interactions and affinity related to SNPs using protein threading modeling, molecular docking, all-atom, coarse grained, and umbrella sampling molecular dynamics simulations (AA-MDS and CG-MDS, respectively). The results showed that the three MRP4 structures had significantly different conformations at given sites, leading to differences in the docking scores (DS) and binding sites of three different groups of molecules. Folic acid (FA) had the highest variation in DS on G187W concerning WT-MRP4. WT-MRP4, G187W, Y556C, and FA had different conformations through 25 ns AA-MD. Umbrella sampling simulations indicated that the Y556C-FA complex was the most stable one with or without ATP. In Y556C, the cyclic adenosine monophosphate (cAMP) and ceefourin-1 binding sites are located out of the entrance of the inner cavity, which suggests that both cAMP and ceefourin-1 may not be transported. The binding site for cAMP and ceefourin-1 is quite similar and the affinity (binding energy) of ceefourin-1 to WT-MRP4, G187W, and Y556C is greater than the affinity of cAMP, which may suggest that ceefourin-1 works as a competitive inhibitor. In conclusion, the nsSNPs G187W and Y556C lead to changes in protein conformation, which modifies the ligand binding site, DS, and binding energy.

Trichinella spiralis Calreticulin S-Domain Binds to Human Complement C1q to Interfere With C1q-Mediated Immune Functions

Helminths develop strategies to escape host immune responses that facilitate their survival in the hostile host immune environment. Trichinella spiralis, a tissue-dwelling nematode, has developed a sophisticated strategy to escape complement attack. Our previous study demonstrated that T. spiralis secretes calreticulin (TsCRT) to inhibit host classical complement activation through binding to C1q; however, the C1q binding site in TsCRT and the specific mechanism involved with complement-related immune evasion remains unknown. Using molecular docking modeling and fragment expression, we determined that TsCRT-S, a 153-aa domain of TsCRT, is responsible for C1q binding. Recombinant TsCRT-S protein expressed in Escherichia coli had the same capacity to bind and inhibit human C1q-induced complement and neutrophil activation, as full-length TsCRT. TsCRT-S inhibited neutrophil reactive oxygen species and elastase release by binding to C1q and reduced neutrophil killing of newborn T. spiralis larvae. Binding of TsCRT-S to C1q also inhibited formation of neutrophil extracellular traps (NETs), which are involved in autoimmune pathologies and have yet to be therapeutically targeted. These findings provide evidence that the TsCRT-S fragment, rather than the full-length TsCRT, is a potential target for vaccine or therapeutic development for trichinellosis, as well as for complement-related autoimmune disease therapies.

Ultrasensitive, Specific, and Rapid Fluorescence Turn-On Nitrite Sensor Enabled by Precisely Modulated Fluorophore Binding

The precise regulation of fluorophore binding sites in an organic probe is of great significance toward the design of fluorescent sensing materials with specific functions. In this study, a probe with specific fluorescence properties and nitrite detection ability is designed by precisely modulating benzothiazole binding sites. Only the fluorophore bond at the ortho-position of the aniline moiety can specifically recognize nitrite, which ensures that the reaction products displays a robust green emission. The unique 2-(2-amino-4-carboxyphenyl) benzothiazole (ortho-BT) shows superior nitrite detection performance, including a low detection limit (2.2 fg), rapid detection time (<5 s), and excellent specificity even in the presence of >40 types of strong redox active, colored substances, nitro compounds, and metal ions. Moreover, the probe is highly applicable for the rapid on-site and semiquantitative measurement of nitrite. The proposed probe design strategy is expected to start a new frontier for the exploration of probe design methodology.

Two RNA Tunnel Inhibitors Bind in Highly Conserved Sites in Dengue Virus NS5 Polymerase: Structural and Functional Studies

Dengue virus (DENV) NS5 RNA-dependent RNA polymerase (RdRp), an important drug target, synthesizes viral RNA and is essential for viral replication. While a number of allosteric inhibitors have been reported for hepatitis C virus RdRp, few have been described for DENV RdRp. Following a diverse compound screening campaign and a rigorous hit-to-lead flowchart combining biochemical and biophysical approaches, two DENV RdRp nonnucleoside inhibitors were identified and characterized. These inhibitors show low- to high-micromolar inhibition in DENV RNA polymerization and cell-based assays. X-ray crystallography reveals that they bind in the enzyme RNA template tunnel. One compound (NITD-434) induced an allosteric pocket at the junction of the fingers and palm subdomains by displacing residue V603 in motif B. Binding of another compound (NITD-640) ordered the fingers loop preceding the F motif, close to the RNA template entrance. Most of the amino acid residues that interacted with these compounds are highly conserved in flaviviruses. Both sites are important for polymerase de novo initiation and elongation activities and essential for viral replication. This work provides evidence that the RNA tunnel in DENV RdRp offers interesting target sites for inhibition.IMPORTANCE Dengue virus (DENV), an important arthropod-transmitted human pathogen that causes a spectrum of diseases, has spread dramatically worldwide in recent years. Despite extensive efforts, the only commercial vaccine does not provide adequate protection to naive individuals. DENV NS5 polymerase is a promising drug target, as exemplified by the development of successful commercial drugs against hepatitis C virus (HCV) polymerase and HIV-1 reverse transcriptase. High-throughput screening of compound libraries against this enzyme enabled the discovery of inhibitors that induced binding sites in the RNA template channel. Characterizations by biochemical, biophysical, and reverse genetics approaches provide a better understanding of the biological relevance of these allosteric sites and the way forward to design more-potent inhibitors.

How to use the MEROPS database and website to help understand peptidase specificity

The MEROPS website (https://www.ebi.ac.uk/merops) and database was established in 1996 to present the classification and nomenclature of proteolytic enzymes. This was expanded to include a classification of protein inhibitors of proteolytic enzymes in 2004. Each peptidase or inhibitor is assigned to a distinct identifier, based on its biochemical and biological properties, and homologous sequences are assembled into a family. Families in which the proteins share similar tertiary structures are assembled into a clan. The MEROPS classification is thus a hierarchy with at least three levels (protein-species, family, and clan) showing the evolutionary relationship. Several other data collections have been assembled, which are accessed from all levels in the hierarchy. These include, sequence homologs, selective bibliographies, substrate cleavage sites, peptidase-inhibitor interactions, alignments, and phylogenetic trees. The substrate cleavage collection has been assembled from the literature and includes physiological, pathological, and nonphysiological cleavages in proteins, peptides, and synthetic substrates. In this article, we make recommendations about how best to analyze these data and show analyses to indicate peptidase binding site preferences and exclusions. We also identify peptidases where co-operative binding occurs between adjacent binding sites.