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.

Effect of ultrasound treatment on porcine myofibrillar protein binding furan flavor compounds at different salt concentrations

The effects of ultrasound (500 W) on the interaction of porcine myofibrillar protein (MP) with furan flavor compounds at different salt concentrations (0.6 %, 1.2 % and 2.4 %) were investigated. With the increase of salt concentration, the particle size of MP decreased, and the surface hydrophobicity and active sulfhydryl content increased due to the unfolding and depolymerization of MP. At the same time, ultrasound promoted the exposure of hydrophobic binding sites and hydrogen bonding sites of MP in different salt concentration systems, thus improving the binding ability of MP with furan compounds by 2 % to 22 %, among which MP had the strongest binding capacity of 2-pentylfuran. In conclusion, ultrasound could effectively promote the unfolding of the secondary structure of MP, which was beneficial to the combination of MP and furan flavor compounds under different salt concentrations.

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.

SM protein Sly1 and a SNARE Habc domain promote membrane fusion through multiple mechanisms

SM proteins including Sly1 are essential cofactors of SNARE-mediated membrane fusion. Using SNARE and Sly1 mutants and chemically defined in vitro assays, we separate and assess proposed mechanisms through which Sly1 augments fusion: (i) opening the closed conformation of the Qa-SNARE Sed5; (ii) close-range tethering of vesicles to target organelles, mediated by the Sly1-specific regulatory loop; and (iii) nucleation of productive trans-SNARE complexes. We show that all three mechanisms are important and operate in parallel, and that close-range tethering promotes trans-complex assembly when cis-SNARE assembly is a competing process. Further, we demonstrate that the autoinhibitory N-terminal Habc domain of Sed5 has at least two positive activities: it is needed for correct Sed5 localization, and it directly promotes Sly1-dependent fusion. "Split Sed5," with Habc presented solely as a soluble fragment, can function both in vitro and in vivo. Habc appears to facilitate events leading to lipid mixing rather than promoting opening or stability of the fusion pore.

Effective expression and characterization of the receptor binding domains in SARS-CoV-2 Spike proteins from original strain and variants of concern using Bombyx mori nucleopolyhedrovirus in silkworm

A new coronavirus, known as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is responsible for the global pandemic of COVID-19 in 2020. Through structural analysis, it was found that several amino acid residues in the human angiotensin-converting enzyme-2 (hACE2) receptor directly interact with those in the receptor binding domain (RBD) of the spike glycoprotein (S-protein). Various cell lines, including HEK293, HeLa cells, and the baculovirus expression vector system (BEVS) with the insect cell line Sf9, have been utilized to produce the RBD. In this study, we investigated the use of Bombyx mori nucleopolyhedrovirus (BmNPV) and BEVS. For efficient production of a highly pure recombinant RBD protein, we designed it with two tags (His tag and STREP tag) at the C-terminus and a solubilizing tag (SUMO) at the N-terminus. After expressing the protein using BmNPV and silkworm and purifying it with a HisTrap excel column, the eluted protein was digested with SUMO protease and further purified using a Strep-Tactin Superflow column. As a result, we obtained the RBD as a monomer with a yield of 2.6 mg/10 mL serum (equivalent to 30 silkworms). The RBD showed an affinity for the hACE2 receptor. Additionally, the RBDs from the Alpha, Beta, Gamma, Delta, and Omicron variants were expressed and purified using the same protocol. It was found that the RBD from the Alpha, Beta, Gamma, and Delta variants could be obtained with yields of 1.4-2.6 mg/10 mL serum and had an affinity to the hACE2 receptor.

Rational synthesis of a heparan sulfate saccharide that promotes the activity of BMP2

Heparan sulfate (HS) is a glycosaminoglycan (GAG) found throughout nature and is involved in a wide range of functions including modulation of cell signalling via sequestration of growth factors. Current consensus is that the specificity of HS motifs for protein binding are individual for each protein. Given the structural complexity of HS the synthesis of libraries of these compounds to probe this is not trivial. Herein we present the synthesis of an HS decamer, the design of which was undertaken rationally from previously published data for HS binding to the growth factor BMP-2. The biological activity of this HS decamer was assessed in vitro, showing that it had the ability to both bind BMP-2 and increase its thermal stability as well as enhancing the bioactivity of BMP-2 in vitro in C2C12 cells. At the same time no undesired anticoagulant effect was observed. This decamer was then analysed in vivo in a rabbit model where higher bone formation, bone mineral density (BMD) and trabecular thickness were observed over an empty defect or collagen implant alone. This indicated that the HS decamer was effective in promoting bone regeneration in vivo.

Agonists or positive allosteric modulators of α7 nicotinic acetylcholine receptor prevent interaction of SARS-Cov-2 receptor-binding domain with astrocytoma cells

SARS-Cov-2, the virus causing COVID-19, penetrates host target cells via the receptor of angiotensin-converting enzyme 2 (ACE2). Disrupting the virus interaction with ACE2 affords a plausible mechanism for prevention of cell penetration and inhibiting dissemination of the virus. Our studies demonstrate that ACE2 interaction with the receptor binding domain of SARS-Cov-2 spike protein (RBD) can be impaired by modulating the α7 nicotinic acetylcholine receptor (α7 nAChR) contiguous with ACE2. U373 cells of human astrocytoma origin were shown to bind both ACE2-specific antibody and recombinant RBD in Cell-ELISA. ACE2 was found to interact with α7 nAChR in U373 cell lysates studied by Sandwich ELISA. Our studies demonstrate that inhibition of RBD binding to ACE2-expressing U373 cells were defined with α7 nAChR agonists choline and PNU282987, but not a competitive antagonist methyllicaconitine (MLA). Additionally, the type 2 positive allosteric modulator (PAM2) PNU120596 and hydroxyurea (HU) also inhibited the binding. Our studies demonstrate that activation of α7 AChRs has efficacy in inhibiting the SARS-Cov-2 interaction with the ACE2 receptor and in such a way can prevent virus target cell penetration. These studies also help to clarify the consistent efficacy and positive outcomes for utilizing HU in treating COVID-19.

Multi-spectroscopic and molecular docking studies for the pH-dependent interaction of β-lactoglobulin with (-)-epicatechin gallate and/or piceatannol: Influence on antioxidant activity and stability

(-)-Epicatechin gallate (ECG) and piceatannol (PIC) are commonly polyphenols with excellent biological activities. β-Lactoglobulin (BLG) is a food-grade globule protein and its morphologies are sensitive to pH. This study used experimental and computational methods to determine the interaction of single or combined ECG and PIC with BLG at different pHs. The static quenching process was determined through fluorescence and ultraviolet-visible spectroscopy. Compared with ECG, PIC could significantly bind to BLG with higher affinity. Their binding affinity for BLG with different morphologies followed the tendency of monomer > dimer > tetramer. The negative contribution of van der Waals forces, electrostatic interactions, and hydrogen bonds to ΔHo exceeded the positive contribution of hydrophobic interactions in the spontaneous and exothermic process. The reduced binding affinity in the ternary systems demonstrated the competitive binding between ECG and PIC on BLG, and the hinder effect of ECG or PIC was enhanced with increasing pH. Molecular docking studies revealed the same binding sites of ECG and PIC on various conformations of BLG and identical driven forces as thermodynamic results. Tryptophan and tyrosine were the main participators in the BLG + ECG and BLG + PIC systems, respectively. The conformational changes in the binary and ternary systems could be ascertained through synchronous fluorescence, circular dichroism, and dynamic light scattering. Furthermore, the effects of pH and BLG encapsulation on the antioxidant capacity and stability of ECG or PIC were also implemented. ECG or PIC was the most stable in the (BLG + PIC) + ECG system at pH 6.0. This study could clarify the interaction mechanism between ECG/PIC and BLG and elucidate the pH effect on their binding information. The results will provide basic support for their usage in food processing and applications.

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.

Raman spectroscopy study of 7,8-dihydrofolate inhibition on the Wuhan strain SARS-CoV-2 binding to human ACE2 receptor

Emerging evidence suggests that elevated levels of folic acid in the bloodstream may confer protection against Wuhan-SARS-CoV-2 infection and mitigate its associated symptoms. Notably, two comprehensive studies of COVID-19 patients in Israel and UK uncovered a remarkable trend, wherein individuals with heightened folic acid levels exhibited only mild symptoms and necessitated no ventilatory support. In parallel, research has underscored the potential connection between decreased folic acid levels and the severity of Covid-19 among hospitalized patients. Yet, the underlying mechanisms governing this intriguing inhibition remain elusive. In a quest to elucidate these mechanisms, we conducted a molecular dynamics simulation approach followed by a Raman spectroscopy study to delve into the intricate interplay between the folic acid metabolite, 7,8-dihydrofolate (DHF), and the angiotensin-converting enzyme ACE2 receptor, coupled with its interaction with the receptor-binding domain (RBD) of the Wuhan strain of SARS-CoV-2. Through a meticulous exploration, we scrutinized the transformation of the ACE2 + RBD complex, allowing these reactants to form bonds. This was juxtaposed with a similar investigation where ACE2 was initially permitted to react with DHF, followed by the exposure of the ACE2 + DHF complex to RBD. We find that DHF, when bonded to ACE2, functions as a physical barrier, effectively inhibiting the binding of the Wuhan strain RBD. This physicochemical process offers a cogent explanation for the observed inhibition of host cell infection in subjects receiving supplementary folic acid doses, as epidemiologically substantiated in multiple studies. This study not only sheds light on a potential avenue for mitigating SARS-CoV-2 infection but also underscores the crucial role of folic acid metabolites in host-virus interactions. This research paves the way for novel therapeutic strategies in the battle against COVID-19 and reinforces the significance of investigating the molecular mechanisms underlying the protective effects of folic acid in the context of viral infections.

Computational methods in the analysis of SARS-CoV-2 in mammals: A systematic review of the literature

SARS-CoV-2 is an enveloped RNA virus that causes severe respiratory illness in humans and animals. It infects cells by binding the Spike protein to the host's angiotensin-converting enzyme 2 (ACE2). The bat is considered the natural host of the virus, and zoonotic transmission is a significant risk and can happen when humans come into close contact with infected animals. Therefore, understanding the interconnection between human, animal, and environmental health is important to prevent and control future coronavirus outbreaks. This work aimed to systematically review the literature to identify characteristics that make mammals suitable virus transmitters and raise the main computational methods used to evaluate SARS-CoV-2 in mammals. Based on this review, it was possible to identify the main factors related to transmissions mentioned in the literature, such as the expression of ACE2 and proximity to humans, in addition to identifying the computational methods used for its study, such as Machine Learning, Molecular Modeling, Computational Simulation, between others. The findings of the work contribute to the prevention and control of future outbreaks, provide information on transmission factors, and highlight the importance of advanced computational methods in the study of infectious diseases that allow a deeper understanding of transmission patterns and can help in the development of more effective control and intervention strategies.

Studying the interaction of glycans with intact virions and virus-like particles by ligand-observed NMR spectroscopy

Virus-glycan interactions play a crucial role in the infection process of many viruses. NMR spectroscopy has emerged as a powerful tool for studying these interactions at the molecular level. In this article, we review several published papers and reports that have highlighted the application of NMR spectroscopy in understanding the complex questions of how viruses engage with and bind to receptor glycans. The use of saturation transfer difference (STD) NMR spectroscopy has demonstrated itself as highly advantageous in investigating the interaction between glycans and intact virions or virus-like particles (VLPs). The broad NMR signal linewidth of virions and VLPs allows efficient saturation without affecting the glycan signals. The advantage of this approach is that the viral capsid environment in protein organization and function is not ignored and therefore provides a more biologically relevant model for exploring the interactions between the virus and the host cell glycans. We will review some examples of using NMR spectroscopy to study influenza cell tropism, rotaviruses, and noroviruses.

Disruption of Ca2+/calmodulin:KSR1 interaction lowers ERK activation

KSR1, a key scaffold protein for the MAPK pathway, facilitates ERK activation upon growth factor stimulation. We recently demonstrated that KSR1 binds the Ca2+-binding protein calmodulin (CaM), thereby providing an intersection between KSR1-mediated and Ca2+ signaling. In this study, we set out to generate a KSR1 point mutant with reduced Ca2+/CaM binding in order to unravel the functional implications of their interaction. To do so, we solved the structural determinants of complex formation. Using purified fragments of KSR1, we showed that Ca2+/CaM binds to the CA3 domain of KSR1. We then used in silico molecular modeling to predict contact residues for binding. This approach identified two possible modes of interaction: (1) binding of extended Ca2+/CaM to a globular conformation of KSR1-CA3 via electrostatic interactions or (2) binding of collapsed Ca2+/CaM to α-helical KSR1-CA3 via hydrophobic interactions. Experimentally, site-directed mutagenesis of the predicted contact residues for the two binding models favored that where collapsed Ca2+/CaM binds to the α-helical conformation of KSR1-CA3. Importantly, replacing KSR1-Phe355 with Asp reduces Ca2+/CaM binding by 76%. The KSR1-F355D mutation also significantly impairs the ability of EGF to activate ERK, which reveals that Ca2+/CaM binding promotes KSR1-mediated MAPK signaling. This work, by uncovering structural insight into the binding of KSR1 to Ca2+/CaM, identifies a KSR1 single-point mutant as a bioreagent to selectively study the crosstalk between Ca2+ and KSR1-mediated signaling.

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.

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 structural origins of complex charge heterogeneity in therapeutic ACE2Fc fusion protein facilitated by free-flow isoelectric focusing

Fc Fusion protein represents a versatile molecular platform with considerable potential as protein therapeutics of which the charge heterogeneity should be well characterized according to regulatory guidelines. Angiotensin-converting enzyme 2 Fc fusion protein (ACE2Fc) has been investigated as a potential neutralizing agent to various coronaviruses, including the lingering SARS-CoV-2, as this coronavirus must bind to ACE2 to allow for its entry into host cells. ACE2Fc, an investigational new drug developed by Henlius (Shanghai China), has passed the Phase I clinical trial, but its huge amount of charge isoforms and complicated charge heterogeneity posed a challenge to charge variant investigation in pharmaceutical development. We employed offline free-flow isoelectric focusing (FF-IEF) fractionation, followed by detailed characterization of enriched ACE2Fc fractions, to unveil the structural origins of charge heterogeneity in ACE2Fc expressed by recombinant CHO cells. We adopted a well-tuned 3-component separation medium for ACE2Fc fractionation, the highest allowable voltage to maximize the FF-IEF separation window and a mild Protein A elution method for preservation of protein structural integrity. Through peptide mapping and other characterizations, we revealed that the intricate profiles of ACE2Fc charge heterogeneity are mainly caused by highly sialylated multi-antenna N-glycosylation. In addition, based on fraction characterization and in silico glycoprotein model analysis, we discovered that the large acidic glycans at N36, N73, and N305 of ACE2Fc were able to decrease the binding activity towards Spike (S) protein of SARS-CoV-2. Our study exemplifies the value of FF-IEF in highly complex fusion protein characterization and revealed a quantitative sialylation-activity relationship in ACE2Fc.

Identification of Factor XIII β-Sandwich Residues Mediating Glutamine Substrate Binding and Activation Peptide Cleavage

BACKGROUND

 Factor XIII (FXIII) forms covalent crosslinks across plasma and cellular substrates and has roles in hemostasis, wound healing, and bone metabolism. FXIII activity is implicated in venous thromboembolism (VTE) and is a target for developing pharmaceuticals, which requires understanding FXIII - substrate interactions. Previous studies proposed the β-sandwich domain of the FXIII A subunit (FXIII-A) exhibits substrate recognition sites.

MATERIAL AND METHODS

 Recombinant FXIII-A proteins (WT, K156E, F157L, R158Q/E, R171Q, and R174E) were generated to identify FXIII-A residues mediating substrate recognition. Proteolytic (FXIII-A*) and non-proteolytic (FXIII-A°) forms were analyzed for activation and crosslinking activities toward physiological substrates using SDS-PAGE and MALDI-TOF MS.

RESULTS

 All FXIII-A* variants displayed reduced crosslinking abilities compared to WT for Fbg αC (233 - 425), fibrin, and actin. FXIII-A* WT activity was greater than A°, suggesting the binding site is more exposed in FXIII-A*. With Fbg αC (233 - 425), FXIII-A* variants R158Q/E, R171Q, and R174E exhibited decreased activities approaching those of FXIII-A°. However, with a peptide substrate, FXIII-A* WT and variants showed similar crosslinking suggesting the recognition site is distant from the catalytic site. Surprisingly, FXIII-A R158E and R171Q displayed slower thrombin activation than WT, potentially due to loss of crucial H-bonding with neighboring activation peptide (AP) residues.

CONCLUSION

 In conclusion, FXIII-A residues K156, F157, R158, R171, and R174 are part of a binding site for physiological substrates [fibrin (α and γ) and actin]. Moreover, R158 and R171 control AP cleavage during thrombin activation. These investigations provide new molecular details on FXIII - substrate interactions that control crosslinking abilities.

Interaction of T-cell-specific adapter protein with Src- and Tec-family kinases

Adapter proteins are flexible and dynamic modulators of cellular signalling that are important for immune cell function. One of these, the T-cell-specific adapter protein (TSAd), interacts with the non-receptor tyrosine kinases Src and Lck of the Src family kinases (SFKs) and Itk of the Tec family kinases (TFKs). Three tyrosine residues in the TSAd C-terminus are phosphorylated by Lck and serve as docking sites for the Src homology 2 (SH2) domains of Src and Lck. The TSAd proline-rich region (PRR) binds to the Src homology 3 (SH3) domains found in Lck, Src and Itk. Despite known interactors, the role TSAd plays in cellular signalling remains largely unknown. TSAd's ability to bind both SFKs and TFKs may point to its function as a general scaffold for both kinase families. Using GST-pulldown as well as peptide array experiments, we found that both the SH2 and SH3 domains of the SFKs Fyn and Hck, as well as the TFKs Tec and Txk, interact with TSAd. This contrasts with Itk, which interacts with TSAd only through its SH3 domain. Although our analysis showed that TSAd is both co-expressed and may interact with Fyn, we were unable to co-precipitate Fyn with TSAd from Jurkat cells, as detected by Western blotting and affinity purification mass spectrometry. This may suggest that TSAd-Fyn interaction in intact cells may be limited by other factors, such as the subcellular localization of the two molecules or the co-expression of competing binding partners.

Deciphering ACE2-RBD binding affinity through peptide scanning: A molecular dynamics simulation approach

Rapid discovery of target information for protein-protein interactions (PPIs) is significant in drug design, diagnostics, vaccine development, antibody therapy, etc. Peptide microarray is an ideal tool for revealing epitope information of PPIs. In this work, the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) spike receptor-binding domain (RBD) and the host cell receptor angiotensin-converting enzyme 2 (ACE2) were introduced as a model to study the epitope information of RBD-specific binding to ACE2 via a combination of theoretical calculations and experimental validation. Through dock and molecular dynamics simulations, it was found that among the 22 peptide fragments that consist of RBD, #14 (YNYLYRLFRKSNLKP) has the highest binding strength. Subsequently, the experiments of peptide microarray constructed based on plasmonic materials chip also confirmed the theoretical calculation data. Compared to other methods, such as phage display technology and surface plasmon resonance (SPR), this method is rapid and cost-effective, providing insights into the investigation of pathogen invasion processes and the timely development of peptide drugs and other fields.

The complexities of ligand/receptor interactions: Exploring the role of molecular vibrations and quantum tunnelling

Molecular vibrations and quantum tunneling may link ligand binding to the function of pharmacological receptors. The well-established lock-and-key model explains a ligand's binding and recognition by a receptor; however, a general mechanism by which receptors translate binding into activation, inactivation, or modulation remains elusive. The Vibration Theory of Olfaction was proposed in the 1930s to explain this subset of receptor-mediated phenomena by correlating odorant molecular vibrations to smell, but a mechanism was lacking. In the 1990s, inelastic electron tunneling was proposed as a plausible mechanism for translating molecular vibration to odorant physiology. More recently, studies of ligands' vibrational spectra and the use of deuterated ligand analogs have provided helpful information to study this admittedly controversial hypothesis in metabotropic receptors other than olfactory receptors. In the present work, based in part on published experiments from our laboratory using planarians as an experimental organism, I will present a rationale and possible experimental approach for extending this idea to ligand-gated ion channels.

Ser252Asn Mutation Introduces a New N-Linked Glycosylation Site and Causes Type IIb Protein C Deficiency

BACKGROUND

 Protein C (PC) is a vitamin K-dependent anticoagulant serine protease zymogen which upon activation by the thrombin-thrombomodulin (TM) complex downregulates the coagulation cascade by degrading cofactors Va and VIIIa by limited proteolysis. We identified a thrombosis patient who carried a heterozygous mutation c.881G > A, p.Ser252Asn (S252N) in PROC. This mutation was originally described in a report of novel mutations in patients presenting with defective PC anticoagulant activity in Paris. The research identified PC-S252N (the "Paris" mutation) in a propositus and her family members and highlighted the critical role of Ser252 in the anticoagulation process of activated PC (APC).

MATERIAL AND METHODS

 We expressed the PC-S252N mutant in mammalian cells and characterized the properties in coagulation assays to decipher the molecular basis of anticoagulant defect of this mutation.

RESULTS

 We demonstrated that PC-S252N had a diminished ability to TM binding, which resulted in its impaired activation by the thrombin-TM complex. However, APC-S252N exhibited a slightly stronger cleavage capacity for the chromogenic substrate. Meanwhile, the catalytic activity of APC-S252N toward FVa was significantly reduced. Sequence analysis revealed that Ser252 to Asn substitution introduced a new potential N-linked glycosylation site (252NTT254) in the catalytic domain of PC, which adversely affected both the activation process of PC and anticoagulant activity of APC.

CONCLUSION

 The new N-glycosylation site (252NTT254) resulting from the mutation of Ser252 to Asn252 in PROC affects the overall structure of the protease, thereby adversely affecting the anticoagulant function of protein C. This modification has a negative impact on both TM-promoted activation of protein C and APC cleavage of FVa, ultimately leading to thrombosis in the patient.

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.

Light-dependent expression and accumulation of miR408-encoded peptide, miPEP408, is regulated by HY5 in Arabidopsis

Recent studies propose that primary transcripts of miRNAs (pri-miRNAs) contain small Open Reading Frames (ORFs) capable of encoding miRNA-encoded peptides (miPEPs). These miPEPs can function as transcriptional regulators for their corresponding pri-miRNAs, ultimately enhancing mature miRNA accumulation. Notably, pri-miR408 encodes the functional peptide miPEP408, regulating expression of miR408 and its target genes, providing plant tolerance to stresses. While miPEPs are crucial regulators, the factors governing them are have not been studied in detail. Here, we explored the light-dependent regulation of miPEP408 in Arabidopsis. Expression analysis during dark-light transitions revealed light-induced transcription and accumulation of the miPEP408. As the promoter of miR408 contains cis-acting elements responsible for binding to the bZIP-type transcription factor ELONGATED HYPOCOTYL5 (HY5), known for light-mediated regulation in plants, we studied its involvement in the regulation of miR408. Analysis of HY5 mutant (hy5-215), complemented line (HY5OX/hy5), and CONSTITUTIVE PHOTOMORPHOGENIC 1 mutant (cop1-4) plants supported HY5's positive regulation of miPEP408. Grafting and GUS assays further suggested the role of HY5 as a shoot-root mobile signal inducing light-dependent miPEP408 expression. This study underscores the regulatory impact of light on small peptides, exemplified by miPEP408, mediated by the key transcription factor HY5.

Identification of the potassium-binding site in serotonin transporter

Clearance of serotonin (5-hydroxytryptamine, 5-HT) from the synaptic cleft after neuronal signaling is mediated by serotonin transporter (SERT), which couples this process to the movement of a Na+ ion down its chemical gradient. After release of 5-HT and Na+ into the cytoplasm, the transporter faces a rate-limiting challenge of resetting its conformation to be primed again for 5-HT and Na+ binding. Early studies of vesicles containing native SERT revealed that K+ gradients can provide an additional driving force, via K+ antiport. Moreover, under appropriate conditions, a H+ ion can replace K+. Intracellular K+ accelerates the resetting step. Structural studies of SERT have identified two binding sites for Na+ ions, but the K+ site remains enigmatic. Here, we show that K+ antiport can drive substrate accumulation into vesicles containing SERT extracted from a heterologous expression system, allowing us to study the residues responsible for K+ binding. To identify candidate binding residues, we examine many cation binding configurations using molecular dynamics simulations, predicting that K+ binds to the so-called Na2 site. Site-directed mutagenesis of residues in this site can eliminate the ability of both K+ and H+ to drive 5-HT accumulation into vesicles and, in patch clamp recordings, prevent the acceleration of turnover rates and the formation of a channel-like state by K+ or H+. In conclusion, the Na2 site plays a pivotal role in orchestrating the sequential binding of Na+ and then K+ (or H+) ions to facilitate 5-HT uptake in SERT.

Exhaustive identification of genome-wide binding events of transcriptional regulators

Genome-wide binding assays aspire to map the complete binding pattern of gene regulators. Common practice relies on replication-duplicates or triplicates-and high stringency statistics to favor false negatives over false positives. Here we show that duplicates and triplicates of CUT&RUN are not sufficient to discover the entire activity of transcriptional regulators. We introduce ICEBERG (Increased Capture of Enrichment By Exhaustive Replicate aGgregation), a pipeline that harnesses large numbers of CUT&RUN replicates to discover the full set of binding events and chart the line between false positives and false negatives. We employed ICEBERG to map the full set of H3K4me3-marked regions, the targets of the co-factor β-catenin, and those of the transcription factor TBX3, in human colorectal cancer cells. The ICEBERG datasets allow benchmarking of individual replicates, comparing the performance of peak calling and replication approaches, and expose the arbitrary nature of strategies to identify reproducible peaks. Instead of a static view of genomic targets, ICEBERG establishes a spectrum of detection probabilities across the genome for a given factor, underlying the intrinsic dynamicity of its mechanism of action, and permitting to distinguish frequent from rare regulation events. Finally, ICEBERG discovered instances, undetectable with other approaches, that underlie novel mechanisms of colorectal cancer progression.

Genomic transcription factor binding site selection is edited by the chromatin remodeling factor CHD4

Biologically precise enhancer licensing by lineage-determining transcription factors enables activation of transcripts appropriate to biological demand and prevents deleterious gene activation. This essential process is challenged by the millions of matches to most transcription factor binding motifs present in many eukaryotic genomes, leading to questions about how transcription factors achieve the exquisite specificity required. The importance of chromatin remodeling factors to enhancer activation is highlighted by their frequent mutation in developmental disorders and in cancer. Here, we determine the roles of CHD4 in enhancer licensing and maintenance in breast cancer cells and during cellular reprogramming. In unchallenged basal breast cancer cells, CHD4 modulates chromatin accessibility. Its depletion leads to redistribution of transcription factors to previously unoccupied sites. During cellular reprogramming induced by the pioneer factor GATA3, CHD4 activity is necessary to prevent inappropriate chromatin opening. Mechanistically, CHD4 promotes nucleosome positioning over GATA3 binding motifs to compete with transcription factor-DNA interaction. We propose that CHD4 acts as a chromatin proof-reading enzyme that prevents unnecessary gene expression by editing chromatin binding activities of transcription factors.

SadNet: a novel multimodal fusion network for protein-ligand binding affinity prediction

Protein-ligand binding affinity prediction plays an important role in the field of drug discovery. Existing deep learning-based approaches have significantly improved the efficiency of protein-ligand binding affinity prediction through their excellent inductive bias capability. However, these methods only focus on fragmented three-dimensional data, which truncates the integrity of pocket data, leading to the neglect of potential long-range interactions. In this paper, we propose a dual-stream framework, with amino acid sequence assisting the atomic data fusion for graph neural network (termed SadNet), to fuse both 3D atomic data and sequence data for more accurate prediction results. In detail, SadNet consists of a pocket module and a sequence module. The sequence module expands the "receptive field" of the pocket module through a mid-term virtual node fusion. To better integrate sequence-level information from the sequence module and 3D structural information from the pocket module, we incorporate structural information for each amino acid within the sequence module. Besides, to better understand the intrinsic relationship between sequences and 3D atomic information, our SadNet utilizes information stacking from both the early stage and later stage. Experimental results on publicly available benchmark datasets demonstrate the superiority of the proposed dual-stream approach over the state-of-the-art alternatives. The code of this work is available online at https://github.com/wardhong/SadNet.

Identification of G-quadruplex-interacting proteins in living cells using an artificial G4-targeting biotin ligase

G-quadruplexes (G4s) are noncanonical nucleic acid structures pivotal to cellular processes and disease pathways. Deciphering G4-interacting proteins is imperative for unraveling G4's biological significance. In this study, we developed a G4-targeting biotin ligase named G4PID, meticulously assessing its binding affinity and specificity both in vitro and in vivo. Capitalizing on G4PID, we devised a tailored approach termed G-quadruplex-interacting proteins specific biotin-ligation procedure (PLGPB) to precisely profile G4-interacting proteins. Implementing this innovative strategy in live cells, we unveiled a cohort of 149 potential G4-interacting proteins, which exhibiting multifaceted functionalities. We then substantiate the directly binding affinity of 7 candidate G4-interacting-proteins (SF3B4, FBL, PP1G, BCL7C, NDUV1, ILF3, GAR1) in vitro. Remarkably, we verified that splicing factor 3B subunit 4 (SF3B4) binds preferentially to the G4-rich 3' splice site and the corresponding splicing sites are modulated by the G4 stabilizer PDS, indicating the regulating role of G4s in mRNA splicing procedure. The PLGPB strategy could biotinylate multiple proteins simultaneously, which providing an opportunity to map G4-interacting proteins network in living cells.

DNA fragility at topologically associated domain boundaries is promoted by alternative DNA secondary structure and topoisomerase II activity

CCCTC-binding factor (CTCF) binding sites are hotspots of genome instability. Although many factors have been associated with CTCF binding site fragility, no study has integrated all fragility-related factors to understand the mechanism(s) of how they work together. Using an unbiased, genome-wide approach, we found that DNA double-strand breaks (DSBs) are enriched at strong, but not weak, CTCF binding sites in five human cell types. Energetically favorable alternative DNA secondary structures underlie strong CTCF binding sites. These structures coincided with the location of topoisomerase II (TOP2) cleavage complex, suggesting that DNA secondary structure acts as a recognition sequence for TOP2 binding and cleavage at CTCF binding sites. Furthermore, CTCF knockdown significantly increased DSBs at strong CTCF binding sites and at CTCF sites that are located at topologically associated domain (TAD) boundaries. TAD boundary-associated CTCF sites that lost CTCF upon knockdown displayed increased DSBs when compared to the gained sites, and those lost sites are overrepresented with G-quadruplexes, suggesting that the structures act as boundary insulators in the absence of CTCF, and contribute to increased DSBs. These results model how alternative DNA secondary structures facilitate recruitment of TOP2 to CTCF binding sites, providing mechanistic insight into DNA fragility at CTCF binding sites.

AMP-activated protein kinase can be allosterically activated by ADP but AMP remains the key activating ligand

The AMP-activated protein kinase (AMPK) is a sensor of cellular energy status. When activated by increases in ADP:ATP and/or AMP:ATP ratios (signalling energy deficit), AMPK acts to restore energy balance. Binding of AMP to one or more of three CBS repeats (CBS1, CBS3, CBS4) on the AMPK-γ subunit activates the kinase complex by three complementary mechanisms: (i) promoting α-subunit Thr172 phosphorylation by the upstream kinase LKB1; (ii) protecting against Thr172 dephosphorylation; (iii) allosteric activation. Surprisingly, binding of ADP has been reported to mimic the first two effects, but not the third. We now show that at physiologically relevant concentrations of Mg.ATP2- (above those used in the standard assay) ADP binding does cause allosteric activation. However, ADP causes only a modest activation because (unlike AMP), at concentrations just above those where activation becomes evident, ADP starts to cause competitive inhibition at the catalytic site. Our results cast doubt on the physiological relevance of the effects of ADP and suggest that AMP is the primary activator in vivo. We have also made mutations to hydrophobic residues involved in binding adenine nucleotides at each of the three γ subunit CBS repeats of the human α2β2γ1 complex and examined their effects on regulation by AMP and ADP. Mutation of the CBS3 site has the largest effects on all three mechanisms of AMP activation, especially at lower ATP concentrations, while mutation of CBS4 reduces the sensitivity to AMP. All three sites appear to be required for allosteric activation by ADP.

Structural and biochemical analysis of family 92 carbohydrate-binding modules uncovers multivalent binding to β-glucans

Carbohydrate-binding modules (CBMs) are non-catalytic proteins found appended to carbohydrate-active enzymes. Soil and marine bacteria secrete such enzymes to scavenge nutrition, and they often use CBMs to improve reaction rates and retention of released sugars. Here we present a structural and functional analysis of the recently established CBM family 92. All proteins analysed bind preferentially to β-1,6-glucans. This contrasts with the diversity of predicted substrates among the enzymes attached to CBM92 domains. We present crystal structures for two proteins, and confirm by mutagenesis that tryptophan residues permit ligand binding at three distinct functional binding sites on each protein. Multivalent CBM families are uncommon, so the establishment and structural characterisation of CBM92 enriches the classification database and will facilitate functional prediction in future projects. We propose that CBM92 proteins may cross-link polysaccharides in nature, and might have use in novel strategies for enzyme immobilisation.

Structure and dynamics of a pentameric KCTD5/CUL3/Gβγ E3 ubiquitin ligase complex

Heterotrimeric G proteins can be regulated by posttranslational modifications, including ubiquitylation. KCTD5, a pentameric substrate receptor protein consisting of an N-terminal BTB domain and a C-terminal domain, engages CUL3 to form the central scaffold of a cullin-RING E3 ligase complex (CRL3KCTD5) that ubiquitylates Gβγ and reduces Gβγ protein levels in cells. The cryo-EM structure of a 5:5:5 KCTD5/CUL3NTD/Gβ1γ2 assembly reveals a highly dynamic complex with rotations of over 60° between the KCTD5BTB/CUL3NTD and KCTD5CTD/Gβγ moieties of the structure. CRL3KCTD5 engages the E3 ligase ARIH1 to ubiquitylate Gβγ in an E3-E3 superassembly, and extension of the structure to include full-length CUL3 with RBX1 and an ARIH1~ubiquitin conjugate reveals that some conformational states position the ARIH1~ubiquitin thioester bond to within 10 Å of lysine-23 of Gβ and likely represent priming complexes. Most previously described CRL/substrate structures have consisted of monovalent complexes and have involved flexible peptide substrates. The structure of the KCTD5/CUL3NTD/Gβγ complex shows that the oligomerization of a substrate receptor can generate a polyvalent E3 ligase complex and that the internal dynamics of the substrate receptor can position a structured target for ubiquitylation in a CRL3 complex.

The Myc-associated zinc finger protein epigenetically controls expression of interferon-γ-stimulated genes by recruiting STAT1 to chromatin

The MYC-Associated Zinc Finger Protein (MAZ) plays important roles in chromatin organization and gene transcription regulation. Dysregulated expression of MAZ causes diseases, such as glioblastoma, breast cancer, prostate cancer, and liposarcoma. Previously, it has been reported that MAZ controls the proinflammatory response in colitis and colon cancer via STAT3 signaling, suggesting that MAZ is involved in regulating immunity-related pathways. However, the molecular mechanism underlying this regulation remains elusive. Here, we investigate the regulatory effect of MAZ on interferon-gamma (IFN-γ)-stimulated genes via STAT1, a protein that plays an essential role in immune responses to viral, fungal, and mycobacterial pathogens. We demonstrate that about 80% of occupied STAT1-binding sites colocalize with occupied MAZ-binding sites in HAP1/K562 cells after IFN-γ stimulation. MAZ depletion significantly reduces STAT1 binding in the genome. By analyzing genome-wide gene expression profiles in the RNA-Seq data, we show that MAZ depletion significantly suppresses a subset of the immune response genes, which include the IFN-stimulated genes IRF8 and Absent in Melanoma 2. Furthermore, we find that MAZ controls expression of the immunity-related genes by changing the epigenetic landscape in chromatin. Our study reveals an important role for MAZ in regulating immune-related gene expression.

LRMP inhibits cAMP potentiation of HCN4 channels by disrupting intramolecular signal transduction

Lymphoid restricted membrane protein (LRMP) is a specific regulator of the hyperpolarization-activated cyclic nucleotide-sensitive isoform 4 (HCN4) channel. LRMP prevents cAMP-dependent potentiation of HCN4, but the interaction domains, mechanisms of action, and basis for isoform-specificity remain unknown. Here, we identify the domains of LRMP essential for this regulation, show that LRMP acts by disrupting the intramolecular signal transduction between cyclic nucleotide binding and gating, and demonstrate that multiple unique regions in HCN4 are required for LRMP isoform-specificity. Using patch clamp electrophysiology and Förster resonance energy transfer (FRET), we identified the initial 227 residues of LRMP and the N-terminus of HCN4 as necessary for LRMP to associate with HCN4. We found that the HCN4 N-terminus and HCN4-specific residues in the C-linker are necessary for regulation of HCN4 by LRMP. Finally, we demonstrated that LRMP-regulation can be conferred to HCN2 by addition of the HCN4 N-terminus along with mutation of five residues in the S5 region and C-linker to the cognate HCN4 residues. Taken together, these results suggest that LRMP inhibits HCN4 through an isoform-specific interaction involving the N-terminals of both proteins that prevents the transduction of cAMP binding into a change in channel gating, most likely via an HCN4-specific orientation of the N-terminus, C-linker, and S4-S5 linker.

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.

ESPDHot: An Effective Machine Learning-Based Approach for Predicting Protein-DNA Interaction Hotspots

Protein-DNA interactions are pivotal to various cellular processes. Precise identification of the hotspot residues for protein-DNA interactions holds great significance for revealing the intricate mechanisms in protein-DNA recognition and for providing essential guidance for protein engineering. Aiming at protein-DNA interaction hotspots, this work introduces an effective prediction method, ESPDHot based on a stacked ensemble machine learning framework. Here, the interface residue whose mutation leads to a binding free energy change (ΔΔG) exceeding 2 kcal/mol is defined as a hotspot. To tackle the imbalanced data set issue, the adaptive synthetic sampling (ADASYN), an oversampling technique, is adopted to synthetically generate new minority samples, thereby rectifying data imbalance. As for molecular characteristics, besides traditional features, we introduce three new characteristic types including residue interface preference proposed by us, residue fluctuation dynamics characteristics, and coevolutionary features. Combining the Boruta method with our previously developed Random Grouping strategy, we obtained an optimal set of features. Finally, a stacking classifier is constructed to output prediction results, which integrates three classical predictors, Support Vector Machine (SVM), XGBoost, and Artificial Neural Network (ANN) as the first layer, and Logistic Regression (LR) algorithm as the second one. Notably, ESPDHot outperforms the current state-of-the-art predictors, achieving superior performance on the independent test data set, with F1, MCC, and AUC reaching 0.571, 0.516, and 0.870, respectively.

Why Does the E1219V Mutation Expand T-Rich PAM Recognition in Cas9 from Streptococcus pyogenes?

Popular RNA-guided DNA endonuclease Cas9 from Streptococcus pyogenes (SpCas9) recognizes the canonical 5'-NGG-3' protospacer adjacent motif (PAM) and triggers double-stranded DNA cleavage activity. Mutations in SpCas9 were demonstrated to expand the PAM readability and hold promise for therapeutic and genome editing applications. However, the energetics of the PAM recognition and its relation to the atomic structure remain unknown. Using the X-ray structure (precatalytic SpCas9:sgRNA:dsDNA) as a template, we calculated the change in the PAM binding affinity in response to SpCas9 mutations using computer simulations. The E1219V mutation in SpCas9 fine-tunes the water accessibility in the PAM binding pocket and promotes new interactions in the SpCas9:noncanonical T-rich PAM, thus weakening the PAM stringency. The nucleotide-specific interaction of two arginine residues (i.e., R1333 and R1335 of SpCas9) ensured stringent 5'-NGG-3' PAM recognition. R1335A substitution (SpCas9R1335A) completely disrupts the direct interaction between SpCas9 and PAM sequences (canonical or noncanonical), accounting for the loss of editing activity. Interestingly, the double mutant (SpCas9R1335A,E1219V) boosts DNA binding affinity by favoring protein:PAM electrostatic contact in a desolvated pocket. The underlying thermodynamics explain the varied DNA cleavage activity of SpCas9 variants. A direct link between the energetics, structures, and activity is highlighted, which can aid in the rational design of improved SpCas9-based genome editing tools.

Cell-free translation system with artificial lipid-monolayer particles as a unique tool for characterizing lipid-monolayer binding proteins

Methods for functional analysis of proteins specifically localizing to lipid monolayers such as rubber particles and lipid droplets are limited. We have succeeded in establishing a system in which artificially prepared lipid monolayer particles are added to a cell-free translation system to confirm the properties of proteins that specifically bind to lipid monolayers in a translation-coupled manner.

Nonbonded Force Field Parameters from MBIS Partitioning of the Molecular Electron Density Improve Binding Affinity Predictions of the T4-Lysozyme Double Mutant

The use of computer simulation for binding affinity prediction is growing in drug discovery. However, its wider use is constrained by the accuracy of the free energy calculations. The key sources of error are the force fields used to depict molecular interactions and insufficient sampling of the configurational space. To improve the quality of the force field, we developed a Python-based computational workflow. The workflow described here uses the minimal basis iterative stockholder (MBIS) method to determine atomic charges and Lennard-Jones parameters from the polarized molecular density. This is done by performing electronic structure calculations on various configurations of the ligand when it is both bound and unbound. In addition, we validated a simulation procedure that accounts for the protein and ligand degrees of freedom to precisely calculate binding free energies. This was achieved by comparing the self-adjusted mixture sampling and nonequilibrium thermodynamic integration methods using various protein and ligand conformations. The accuracy of predicting binding affinity is improved by using MBIS-derived force field parameters and a validated simulation procedure. This improvement surpasses the chemical precision for the eight aromatic ligands, reaching a root-mean-square error of 0.7 kcal/mol.

Database-Driven Identification of Structurally Similar Protein-Protein Interfaces

Analyzing the similarity of protein interfaces in protein-protein interactions gives new insights into protein function and assists in discovering new drugs. Usually, tools that assess the similarity focus on the interactions between two protein interfaces, while sometimes we only have one predicted interface. Herein, we present PiMine, a database-driven protein interface similarity search. It compares interface residues of one or two interacting chains by calculating and searching tetrahedral geometric patterns of α-carbon atoms and calculating physicochemical and shape-based similarity. On a dedicated, tailor-made dataset, we show that PiMine outperforms commonly used comparison tools in terms of early enrichment when considering interfaces of sequentially and structurally unrelated proteins. In an application example, we demonstrate its usability for protein interaction partner prediction by comparing predicted interfaces to known protein-protein interfaces.

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.

ProInterVal: Validation of Protein-Protein Interfaces through Learned Interface Representations

Proteins are vital components of the biological world and serve a multitude of functions. They interact with other molecules through their interfaces and participate in crucial cellular processes. Disruption of these interactions can have negative effects on organisms, highlighting the importance of studying protein-protein interfaces for developing targeted therapies for diseases. Therefore, the development of a reliable method for investigating protein-protein interactions is of paramount importance. In this work, we present an approach for validating protein-protein interfaces using learned interface representations. The approach involves using a graph-based contrastive autoencoder architecture and a transformer to learn representations of protein-protein interaction interfaces from unlabeled data and then validating them through learned representations with a graph neural network. Our method achieves an accuracy of 0.91 for the test set, outperforming existing GNN-based methods. We demonstrate the effectiveness of our approach on a benchmark data set and show that it provides a promising solution for validating protein-protein interfaces.

Molecular Insights into the Variability in Infection and Immune Evasion Capabilities of SARS-CoV-2 Variants: A Sequence and Structural Investigation of the RBD Domain

As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants continuously emerge, an increasing number of mutations are accumulating in the Spike protein receptor-binding domain (RBD) region. Through sequence analysis of various Variants of Concern (VOC), we identified that they predominantly fall within the ο lineage although recent variants introduce any novel mutations in the RBD. Molecular dynamics simulations were employed to compute the binding free energy of these variants with human Angiotensin-converting enzyme 2 (ACE2). Structurally, the binding interface of the ο RBD displays a strong positive charge, complementing the negatively charged binding interface of ACE2, resulting in a significant enhancement in the electrostatic potential energy for the ο variants. Although the increased potential energy is partially offset by the rise in polar solvation free energy, enhanced electrostatic potential contributes to the long-range recognition between the ο variant's RBD and ACE2. We also conducted simulations of glycosylated ACE2-RBD proteins. The newly emerged ο (JN.1) variant has added a glycosylation site at N-354@RBD, which significantly weakened its binding affinity with ACE2. Further, our interaction studies with three monoclonal antibodies across multiple SARS-CoV-2 strains revealed a diminished neutralization efficacy against the ο variants, primarily attributed to the electrostatic repulsion between the antibodies and RBD interface. Considering the characteristics of the ο variant and the trajectory of emerging strains, we propose that newly developed antibodies against SARS-CoV-2 RBD should have surfaces rich in negative potential and, postbinding, exhibit strong van der Waals interactions. These findings provide invaluable guidance for the formulation of future therapeutic strategies.

Discovery of Covalent Lead Compounds Targeting 3CL Protease with a Lateral Interactions Spiking Neural Network

Covalent drugs exhibit advantages in that noncovalent drugs cannot match, and covalent docking is an important method for screening covalent lead compounds. However, it is difficult for covalent docking to screen covalent compounds on a large scale because covalent docking requires determination of the covalent reaction type of the compound. Here, we propose to use deep learning of a lateral interactions spiking neural network to construct a covalent lead compound screening model to quickly screen covalent lead compounds. We used the 3CL protease (3CL Pro) of SARS-CoV-2 as the screen target and constructed two classification models based on LISNN to predict the covalent binding and inhibitory activity of compounds. The two classification models were trained on the covalent complex data set targeting cysteine (Cys) and the compound inhibitory activity data set targeting 3CL Pro, respected, with good prediction accuracy (ACC > 0.9). We then screened the screening compound library with 6 covalent binding screening models and 12 inhibitory activity screening models. We tested the inhibitory activity of the 32 compounds, and the best compound inhibited SARS-CoV-2 3CL Pro with an IC50 value of 369.5 nM. Further assay implied that dithiothreitol can affect the inhibitory activity of the compound to 3CL Pro, indicating that the compound may covalently bind 3CL Pro. The selectivity test showed that the compound had good target selectivity to 3CL Pro over cathepsin L. These correlation assays can prove the rationality of the covalent lead compound screening model. Finally, covalent docking was performed to demonstrate the binding conformation of the compound with 3CL Pro. The source code can be obtained from the GitHub repository (https://github.com/guzh970630/Screen_Covalent_Compound_by_LISNN).

Molecular Insights into the Differential Effects of Acetylation on the Aggregation of Tau Microtubule-Binding Repeats

Aggregation of tau protein into intracellular fibrillary inclusions is characterized as the hallmark of tauopathies, including Alzheimer's disease and chronic traumatic encephalopathy. The microtubule-binding (MTB) domain of tau, containing either three or four repeats with sequence similarities, plays an important role in determining tau's aggregation. Previous studies have reported that abnormal acetylation of lysine residues displays a distinct effect on the formation of pathological tau aggregates. However, the underlying molecular mechanism remains mostly elusive. In this study, we performed extensive replica exchange molecular dynamics (REMD) simulations of 144 μs in total to systematically investigate the dimerization of four tau MTB repeats and explore the impacts of Lys280 (K280) or Lys321 (K321) acetylation on the conformational ensembles of the R2 or R3 dimer. Our results show that R3 is the most prone to aggregation among the four repeats, followed by R2 and R4, while R1 displays the weakest aggregation propensity with a disordered structure. Acetylation of K280 could promote the aggregation of R2 peptides by increasing the formation of β-sheet structures and strengthening the interchain interaction. However, K321 acetylation decreases the β-sheet content of the R3 dimer, reduces the ability of R3 peptides to form long β-strands, and promotes the stable helix structure formation. The salt bridge and Y310-Y310 π-π stacking interactions of the R3 dimer are greatly weakened by K321 acetylation, resulting in the inhibition of dimerization. This study uncovers the structural ensembles of tau MTB repeats and provides mechanistic insights into the influences of acetylation on tau aggregation, which may deepen the understanding of the pathogenesis of tauopathies.

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/ .