Search

Nucleotide excision repair (NER) clears genomes of DNA adducts formed by UV light, environmental agents, and antitumor drugs. Gene mutations that lead to defects in the core NER reaction cause the skin cancer-prone disease xeroderma pigmentosum. In NER, DNA lesions are excised within an oligonucleotide of 25-30 residues via a complex, multi-step reaction that is regulated by protein-protein interactions. These interactions were first characterized in the 1990s using pull-down, co-IP and yeast two-hybrid assays. More recently, high-resolution structures and detailed functional studies have started to yield detailed pictures of the progression along the NER reaction coordinate. In this review, we highlight how the study of interactions among proteins by structural and/or functional studies have provided insights into the mechanisms by which the NER machinery recognizes and excises DNA lesions. Furthermore, we identify reported, but poorly characterized or unsubstantiated interactions in need of further validation., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Replication Protein A (RPA) is asingle strandedDNA(ssDNA)binding protein that coordinates diverse DNA metabolic processes including DNA replication, repair, and recombination. RPA is a heterotrimeric protein with six functional oligosaccharide/oligonucleotide (OB) domains and flexible linkers. Flexibility enables RPA to adopt multiple configurations andis thought to modulate its function. Here, usingsingle moleculeconfocal fluorescencemicroscopy combinedwith optical tweezers and coarse-grained molecular dynamics simulations, we investigated the diffusional migration of single RPA molecules on ssDNA undertension.The diffusioncoefficientDis the highest (20,000nucleotides 2 /s) at 3pNtension and in 100 mMKCl and markedly decreases whentensionor salt concentrationincreases. We attribute the tension effect to intersegmental transfer which is hindered by DNA stretching and the salt effect to an increase in binding site size and interaction energy of RPA-ssDNA. Our integrative study allowed us to estimate the size and frequency of intersegmental transfer events that occur through transient bridging of distant sites on DNA by multiple binding sites on RPA. Interestingly, deletion of RPA trimeric core still allowed significant ssDNA binding although the reduced contact area made RPA 15-fold more mobile. Finally, we characterized the effect of RPA crowding on RPA migration. These findings reveal how the high affinity RPA-ssDNA interactions are remodeled to yield access, a key step in several DNA metabolic processes., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Chang-Ting Lin was previously employed by LUMICKS B.V and Olivia Yang is currently employed by LUMICS B.V. Y. L. holds The Morton and Gladys Pickman professional chair in Structural Biology., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Cytochrome b 5 (b 5 ) is known to stimulate some catalytic activities of cytochrome P450 (P450, CYP) enzymes, although mechanisms still need to be defined. The reactions most strongly enhanced by b 5 are the 17,20-lyase reactions of P450 17A1 involved in steroid biosynthesis. We had previously used a fluorescently labeled human b 5 variant (Alexa 488-T70C-b 5 ) to characterize human P450 17A1-b 5 interactions, but subsequent proteomic analyses indicated that lysines in b 5 were also modified with Alexa 488 maleimide in addition to Cys-70, due to disulfide dimerization of the T70C mutant. A series of b 5 variants were constructed with Cys replacements for the identified lysine residues and labeled with the dye. Fluorescence attenuation and the function of b 5 in the steroid lyase reaction depended on the modified position. Apo-b 5 (devoid of heme group) studies revealed the lack of involvement of the b 5 heme in the fluorescence attenuation. A structural model of b 5 with P450 17A1 was predicted using AlphaFold-Multimer algorithms/Rosetta docking, based upon the individual structures, which predicted several new contacts not previously reported, that is, interactions of b 5 Glu-48:17A1 Arg-347, b 5 Glu-49:17A1 Arg-449, b 5 Asp-65:17A1 Arg-126, b 5 Asp-65:17A1 Arg-125, and b 5 Glu-61:17A1 Lys-91. Fluorescence polarization assays with two modified b 5 variants yielded K d values (for b 5 -P450 17A1) of 120 to 380 nM, the best estimate of binding affinity. We conclude that both monomeric and dimeric b 5 can bind to P450 17A1 and stimulate activity. Results with the mutants indicate that several Lys residues in b 5 are sensitive to the interaction with P450 17A1, including Lys-88 and Lys-91., Competing Interests: Conflict of interest The authors declare that they have no conflict of interest with the contents of this article., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

RAS is a small GTPase and acts as a binary molecular switch; the transition from its active to inactive state plays a crucial role in various cell signaling processes. Molecular dynamics simulations at the atomistic level suggest that the absence of cofactor Mg 2+ ion generally leads to pronounced structural changes in the Switch-I than Switch-II regions and assists GTP binding. The presence of the Mg 2+ ion also restricts the rotation of ϒ phosphate and enhances the hydrolysis rate of GTP. Further, the simulations reveal that the stability of the protein is almost uncompromised when Mg 2+ is replaced with Zn 2+ and not the Ca 2+ ion. The specificity of H-RAS to GTP was evaluated by substituting with ATP and CTP, which indicates that the binding pocket tolerates purine bases over pyrimidine bases. However, the D119 residue specifically interacts with the guanine base and serves as one of the primary interactions that leads to the selectivity of GTP over ATP. The ring displacement of 32Y serves as gate dynamics in H-RAS which are important for its interaction with GAP for the nucleotide exchange and is restricted in the presence of ATP. Finally, the point mutations 61, 16, and 32 influence the structural changes, specifically in the Switch-II region, which are expected to impact the GTP hydrolysis and thus are termed oncogenic mutations.

Long noncoding RNAs (lncRNAs) are specifically expressed in different diseases and regulate disease progression. To explore the functions of rheumatoid arthritis (RA)-specific lncRNA, we determined the lncRNA expression profile of fibroblast-like synoviocytes (FLS) obtained from patients with RA and osteoarthritis (OA) using a LncRNA microarray and identified up-regulated LncNFYB in RA as a potential therapeutic target. Using gain- and loss-of-function studies, LncNFYB was proven to promote FLS proliferation and cell cycle progress but not affect their invasion, migration, and apoptotic abilities. Further investigation discovered that LncRNA could combine with annexin A2 (ANXA2) and enhance the level of phospho-ANXA2 (Tyr24) in the plasma membrane area, which induced the activation of ERK1/2 to promote proliferation. These findings provide new insights into the biological functions of LncNFYB on modification of FLS, which may be exploited for the therapy of RA., Competing Interests: Conflict of interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Non-coding variants associated with complex traits can alter the motifs of transcription factor (TF)-deoxyribonucleic acid binding. Although many computational models have been developed to predict the effects of non-coding variants on TF binding, their predictive power lacks systematic evaluation. Here we have evaluated 14 different models built on position weight matrices (PWMs), support vector machines, ordinary least squares and deep neural networks (DNNs), using large-scale in vitro (i.e. SNP-SELEX) and in vivo (i.e. allele-specific binding, ASB) TF binding data. Our results show that the accuracy of each model in predicting SNP effects in vitro significantly exceeds that achieved in vivo. For in vitro variant impact prediction, kmer/gkm-based machine learning methods (deltaSVM_HT-SELEX, QBiC-Pred) trained on in vitro datasets exhibit the best performance. For in vivo ASB variant prediction, DNN-based multitask models (DeepSEA, Sei, Enformer) trained on the ChIP-seq dataset exhibit relatively superior performance. Among the PWM-based methods, tRap demonstrates better performance in both in vitro and in vivo evaluations. In addition, we find that TF classes such as basic leucine zipper factors could be predicted more accurately, whereas those such as C2H2 zinc finger factors are predicted less accurately, aligning with the evolutionary conservation of these TF classes. We also underscore the significance of non-sequence factors such as cis-regulatory element type, TF expression, interactions and post-translational modifications in influencing the in vivo predictive performance of TFs. Our research provides valuable insights into selecting prioritization methods for non-coding variants and further optimizing such models., (© The Author(s) 2024. Published by Oxford University Press.)

The functions of many proteins are associated with interconversions among conformational substates. However, these substates can be difficult to measure experimentally, and determining contributions from hydration changes can be especially difficult. Here, we assessed the use of pressure perturbations to sample the substates accessible to the Escherichia coli lactose repressor protein (LacI) in various liganded forms. In the presence of DNA, the regulatory domain of LacI adopts an Open conformation that, in the absence of DNA, changes to a Closed conformation. Increasing the simulation pressure prevented the transition from an Open to a Closed conformation, in a similar manner to the binding of DNA and anti-inducer, ONPF. The results suggest the hydration of specific residues play a significant role in determining the population of different LacI substates and that simulating pressure perturbation could be useful for assessing the role of hydration changes that accompany functionally-relevant amino acid substitutions., Competing Interests: Declaration of Competing Interest The authors declare no competing interests., (Copyright © 2023. Published by Elsevier B.V.)

Mutations in atrial-enriched genes can cause a primary atrial myopathy that can contribute to overall cardiovascular dysfunction. MYBPHL encodes myosin-binding protein H-like (MyBP-HL), an atrial sarcomere protein that shares domain homology with the carboxy-terminus of cardiac myosin-binding protein-C (cMyBP-C). The function of MyBP-HL and the relationship between MyBP-HL and cMyBP-C is unknown. To decipher the roles of MyBP-HL, we used structured illumination microscopy, immuno-electron microscopy, and mass spectrometry to establish the localization and stoichiometry of MyBP-HL. We found levels of cMyBP-C, a major regulator of myosin function, were half as abundant compared to levels in the ventricle. In genetic mouse models, loss of MyBP-HL doubled cMyBP-C abundance in the atria, and loss of cMyBP-C doubled MyBP-HL abundance in the atria. Structured illumination microscopy showed that both proteins colocalize in the C-zone of the A-band, with MyBP-HL enriched closer to the M-line. Immuno-electron microscopy of mouse atria showed MyBP-HL strongly localized 161 nm from the M-line, consistent with localization to the third 43 nm repeat of myosin heads. Both cMyBP-C and MyBP-HL had less-defined sarcomere localization in the atria compared to ventricle, yet areas with the expected 43 nm repeat distance were observed for both proteins. Isometric force measurements taken from control and Mybphl null single atrial myofibrils revealed that loss of Mybphl accelerated the linear phase of relaxation. These findings support a mechanism where MyBP-HL regulates cMyBP-C abundance to alter the kinetics of sarcomere relaxation in atrial sarcomeres., Competing Interests: Competing interests statement:The authors declare no competing interest.

Staphylococcus aureus (S. aureus) is a serious global pathogen that causes a diverse range of invasive diseases. S. aureus utilizes a family of pore-forming toxins, known as bi-component leukocidins, to evade the host immune response and promote infection. Among these is LukAB (leukocidin A/leukocidin B), a toxin that assembles into an octameric β-barrel pore in the target cell membrane, resulting in host cell death. The established cellular receptor for LukAB is CD11b of the Mac-1 complex. Here, we show that hydrogen voltage-gated channel 1 is also required for the cytotoxicity of all major LukAB variants. We demonstrate that while each receptor is sufficient to recruit LukAB to the plasma membrane, both receptors are required for maximal lytic activity. Why LukAB requires two receptors, and how each of these receptors contributes to pore-formation remains unknown. To begin to resolve this, we performed an alanine scanning mutagenesis screen to identify mutations that allow LukAB to maintain cytotoxicity without CD11b. We discovered 30 mutations primarily localized in the stem domains of LukA and LukB that enable LukAB to exhibit full cytotoxicity in the absence of CD11b. Using crosslinking, electron microscopy, and hydroxyl radical protein footprinting, we show these mutations increase the solvent accessibility of the stem domain, priming LukAB for oligomerization. Together, our data support a model in which CD11b binding unlatches the membrane penetrating stem domains of LukAB, and this change in flexibility promotes toxin oligomerization., Competing Interests: Conflict of interest V. J. T. is an inventor on patents and patent applications filed by New York University, which are currently under commercial license to Janssen Biotech Inc. Janssen Biotech provides research funding and other payments associated with the licensing agreement. All other authors declare no conflict of interest., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

The cistrome consists of all cis-acting regulatory elements recognized by transcription factors (TFs). However, only a portion of the cistrome is active for TF binding in a specific tissue. Resolving the active cistrome in plants remains challenging. In this study, we report the assay sequential extraction assisted-active TF identification (sea-ATI), a low-input method that profiles the DNA sequences recognized by TFs in a target tissue. We applied sea-ATI to seven plant tissues to survey their active cistrome and generated 41 motif models, including 15 new models that represent previously unidentified cis-regulatory vocabularies. ATAC-seq and RNA-seq analyses confirmed the functionality of the cis-elements from the new models, in that they are actively bound in vivo, located near the transcription start site, and influence chromatin accessibility and transcription. Furthermore, comparing dimeric WRKY CREs between sea-ATI and DAP-seq libraries revealed that thermodynamics and genetic drifts cooperatively shaped their evolution. Notably, sea-ATI can identify not only positive but also negative regulatory cis-elements, thereby providing unique insights into the functional non-coding genome of plants., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Background: Transcription factors (TFs) exhibit heterogeneous DNA-binding specificities in individual cells and whole organisms under natural conditions, and de novo motif discovery usually provides multiple motifs, even from a single chromatin immunoprecipitation-sequencing (ChIP-seq) sample. Despite the accumulation of ChIP-seq data and ChIP-seq-derived motifs, the diversity of DNA-binding specificities across different TFs and cell types remains largely unexplored., Results: Here, we applied MOCCS2, our k-mer-based motif discovery method, to a collection of human TF ChIP-seq samples across diverse TFs and cell types, and systematically computed profiles of TF-binding specificity scores for all k-mers. After quality control, we compiled a set of TF-binding specificity score profiles for 2,976 high-quality ChIP-seq samples, comprising 473 TFs and 398 cell types. Using these high-quality samples, we confirmed that the k-mer-based TF-binding specificity profiles reflected TF- or TF-family dependent DNA-binding specificities. We then compared the binding specificity scores of ChIP-seq samples with the same TFs but with different cell type classes and found that half of the analyzed TFs exhibited differences in DNA-binding specificities across cell type classes. Additionally, we devised a method to detect differentially bound k-mers between two ChIP-seq samples and detected k-mers exhibiting statistically significant differences in binding specificity scores. Moreover, we demonstrated that differences in the binding specificity scores between k-mers on the reference and alternative alleles could be used to predict the effect of variants on TF binding, as validated by in vitro and in vivo assay datasets. Finally, we demonstrated that binding specificity score differences can be used to interpret disease-associated non-coding single-nucleotide polymorphisms (SNPs) as TF-affecting SNPs and provide candidates responsible for TFs and cell types., Conclusions: Our study provides a basis for investigating the regulation of gene expression in a TF-, TF family-, or cell-type-dependent manner. Furthermore, our differential analysis of binding-specificity scores highlights noncoding disease-associated variants in humans., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Single-stranded DNA-binding proteins (SSBs) play a crucial role in DNA metabolism by binding and stabilizing single-stranded DNA (ssDNA) intermediates. Through their multifaceted roles in DNA replication, recombination, repair, replication restart, and other cellular processes, SSB emerges as a central player in maintaining genomic integrity. These attributes collectively position SSBs as essential guardians of genomic integrity, establishing interactions with an array of distinct proteins. Unlike Escherichia coli , which contains only one type of SSB, some bacteria have two paralogous SSBs, referred to as SsbA and SsbB. In this study, we identified Staphylococcus aureus SsbA (SaSsbA) as a fresh addition to the roster of the anticancer drug 5-fluorouracil (5-FU) binding proteins, thereby expanding the ambit of the 5-FU interactome to encompass this DNA replication protein. To investigate the binding mode, we solved the complexed crystal structure with 5-FU at 2.3 Å (PDB ID 7YM1). The structure of glycerol-bound SaSsbA was also determined at 1.8 Å (PDB ID 8GW5). The interaction between 5-FU and SaSsbA was found to involve R18, P21, V52, F54, Q78, R80, E94, and V96. Based on the collective results from mutational and structural analyses, it became evident that SaSsbA's mode of binding with 5-FU diverges from that of SaSsbB. This complexed structure also holds the potential to furnish valuable comprehension regarding how 5-FU might bind to and impede analogous proteins in humans, particularly within cancer-related signaling pathways. Leveraging the information furnished by the glycerol and 5-FU binding sites, the complexed structures of SaSsbA bring to the forefront the potential viability of several interactive residues as potential targets for therapeutic interventions aimed at curtailing SaSsbA activity. Acknowledging the capacity of microbiota to influence the host's response to 5-FU, there emerges a pressing need for further research to revisit the roles that bacterial and human SSBs play in the realm of anticancer therapy.

Dominant optic atrophy is one of the leading causes of childhood blindness. Around 60-80% of cases 1 are caused by mutations of the gene that encodes optic atrophy protein 1 (OPA1), a protein that has a key role in inner mitochondrial membrane fusion and remodelling of cristae and is crucial for the dynamic organization and regulation of mitochondria 2 . Mutations in OPA1 result in the dysregulation of the GTPase-mediated fusion process of the mitochondrial inner and outer membranes 3 . Here we used cryo-electron microscopy methods to solve helical structures of OPA1 assembled on lipid membrane tubes, in the presence and absence of nucleotide. These helical assemblies organize into densely packed protein rungs with minimal inter-rung connectivity, and exhibit nucleotide-dependent dimerization of the GTPase domains-a hallmark of the dynamin superfamily of proteins 4 . OPA1 also contains several unique secondary structures in the paddle domain that strengthen its membrane association, including membrane-inserting helices. The structural features identified in this study shed light on the effects of pathogenic point mutations on protein folding, inter-protein assembly and membrane interactions. Furthermore, mutations that disrupt the assembly interfaces and membrane binding of OPA1 cause mitochondrial fragmentation in cell-based assays, providing evidence of the biological relevance of these interactions., (© 2023. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Translationally Controlled Tumor Protein (TCTP) serves as a pro-survival factor in tumor cells, inhibiting the mitochondrial apoptosis pathway by enhancing the function of anti-apoptotic Bcl-2 family members Mcl-1 and Bcl-xL. TCTP specifically binds to Bcl-xL, preventing Bax-dependent Bcl-xL-induced cytochrome c release, and it reduces Mcl-1 turnover by inhibiting its ubiquitination, thereby decreasing Mcl-1-mediated apoptosis. TCTP harbors a BH3-like motif that forms a β-strand buried in the globular domain of the protein. In contrast, the crystal structure of the TCTP BH3-like peptide in complex with the Bcl-2 family member Bcl-xL reveals an α-helical conformation for the BH3-like motif, suggesting significant structural changes upon complex formation. Employing biochemical and biophysical methods, including limited proteolysis, circular dichroism, NMR, and SAXS, we describe the TCTP complex with the Bcl-2 homolog Mcl-1. Our findings demonstrate that full-length TCTP binds to the BH3 binding groove of Mcl-1 via its BH3-like motif, experiencing conformational exchange at the interface on a micro- to milli-second timescale. Concurrently, the TCTP globular domain becomes destabilized, transitioning into a molten-globule state. Furthermore, we establish that the non-canonical residue D16 within the TCTP BH3-like motif reduces stability while enhancing the dynamics of the intermolecular interface. In conclusion, we detail the structural plasticity of TCTP and discuss its implications for partner interactions and future anticancer drug design strategies aimed at targeting TCTP complexes., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Unbiased in vivo selections of diverse capsid libraries can yield engineered capsids that overcome gene therapy delivery challenges like traversing the blood-brain barrier (BBB), but little is known about the parameters of capsid-receptor interactions that govern their improved activity. This hampers broader efforts in precision capsid engineering and is a practical impediment to ensuring the translatability of capsid properties between preclinical animal models and human clinical trials. In this work, we utilize the adeno-associated virus (AAV)-PHP.B-Ly6a model system to better understand the targeted delivery and BBB penetration properties of AAV vectors. This model offers a defined capsid-receptor pair that can be used to systematically define relationships between target receptor affinity and in vivo activity of engineered AAV vectors. Here, we report a high-throughput method for quantifying capsid-receptor affinity and demonstrate that direct binding assays can be used to organize a vector library into families with varied affinity for their target receptor. Our data indicate that efficient central nervous system transduction requires high levels of target receptor expression at the BBB, but it is not a requirement for receptor expression to be limited to the target tissue. We observed that enhanced receptor affinity leads to reduced transduction of off-target tissues but can negatively impact on-target cellular transduction and penetration of endothelial barriers. Together, this work provides a set of tools for defining vector-receptor affinities and demonstrates how receptor expression and affinity interact to impact the performance of engineered AAV vectors in targeting the central nervous system. IMPORTANCE Novel methods for measuring adeno-associated virus (AAV)-receptor affinities, especially in relation to vector performance in vivo , would be useful to capsid engineers as they develop AAV vectors for gene therapy applications and characterize their interactions with native or engineered receptors. Here, we use the AAV-PHP.B-Ly6a model system to assess the impact of receptor affinity on the systemic delivery and endothelial penetration properties of AAV-PHP.B vectors. We discuss how receptor affinity analysis can be used to isolate vectors with optimized properties, improve the interpretation of library selections, and ultimately translate vector activities between preclinical animal models and humans., Competing Interests: The authors declare a conflict of interest. This research was supported by Amicus Therapeutics and Passage Bio. JMW is a paid advisor to and holds equity in iECURE, Scout Bio, Passage Bio, and the Center for Breakthrough Medicines (CBM). He also holds equity in the G2 Bio-associated asset companies. He has sponsored research agreements with Amicus Therapeutics, CBM, Elaaj Bio, FA212, G2 Bio, G2 Bio-associated asset companies, iECURE, Passage Bio, and Scout Bio, which are licensees of Penn technology. JMW and JJS are inventors on patents that have been licensed to various biopharmaceutical companies and for which they may receive payments.

The maturation of RNA from its nascent transcription to ultimate utilization (e.g., translation, miR-mediated RNA silencing, etc.) involves an intricately coordinated series of biochemical reactions regulated by RNA-binding proteins (RBPs). Over the past several decades, there has been extensive effort to elucidate the biological factors that control specificity and selectivity of RNA target binding and downstream function. Polypyrimidine tract binding protein 1 (PTBP1) is an RBP that is involved in all steps of RNA maturation and serves as a key regulator of alternative splicing, and therefore, understanding its regulation is of critical biologic importance. While several mechanisms of RBP specificity have been proposed (e.g., cell-specific expression of RBPs and secondary structure of target RNA), recently, protein-protein interactions with individual domains of RBPs have been suggested to be important determinants of downstream function. Here, we demonstrate a novel binding interaction between the first RNA recognition motif 1 (RRM1) of PTBP1 and the prosurvival protein myeloid cell leukemia-1 (MCL1). Using both in silico and in vitro analyses, we demonstrate that MCL1 binds a novel regulatory sequence on RRM1. NMR spectroscopy reveals that this interaction allosterically perturbs key residues in the RNA-binding interface of RRM1 and negatively impacts RRM1 association with target RNA. Furthermore, pulldown of MCL1 by endogenous PTBP1 verifies that these proteins interact in an endogenous cellular environment, establishing the biological relevance of this binding event. Overall, our findings suggest a novel mechanism of regulation of PTBP1 in which a protein-protein interaction with a single RRM can impact RNA association., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Cyclic AMP receptor proteins (CRPs) are important transcription regulators in many species. The prediction of CRP-binding sites was mainly based on position-weighted matrixes (PWMs). Traditional prediction methods only considered known binding motifs, and their ability to discover inflexible binding patterns was limited. Thus, a novel CRP-binding site prediction model called CRPBSFinder was developed in this research, which combined the hidden Markov model, knowledge-based PWMs and structure-based binding affinity matrixes. We trained this model using validated CRP-binding data from Escherichia coli and evaluated it with computational and experimental methods. The result shows that the model not only can provide higher prediction performance than a classic method but also quantitatively indicates the binding affinity of transcription factor binding sites by prediction scores. The prediction result included not only the most knowns regulated genes but also 1089 novel CRP-regulated genes. The major regulatory roles of CRPs were divided into four classes: carbohydrate metabolism, organic acid metabolism, nitrogen compound metabolism and cellular transport. Several novel functions were also discovered, including heterocycle metabolic and response to stimulus. Based on the functional similarity of homologous CRPs, we applied the model to 35 other species. The prediction tool and the prediction results are online and are available at: https://awi.cuhk.edu.cn/∼CRPBSFinder., (© The Author(s) 2023. Published by Oxford University Press.)

Association of single nucleotide polymorphisms in the IL-23 receptor with several auto-inflammatory diseases, led to the heterodimeric receptor and its cytokine-ligand IL-23, becoming important drug targets. Successful antibody-based therapies directed against the cytokine have been licenced and a class of small peptide antagonists of the receptor have entered clinical trials. These peptide antagonists may offer therapeutic advantages over existing anti-IL-23 therapies, but little is known about their molecular pharmacology. In this study, we use a fluorescent version of IL-23 to characterise antagonists of the full-length receptor expressed by living cells using a NanoBRET competition assay. We then develop a cyclic peptide fluorescent probe, specific to the IL23p19:IL23R interface and use this molecule to characterise further receptor antagonists. Finally, we use the assays to study the immunocompromising C115Y IL23R mutation, demonstrating that the mechanism of action is a disruption of the binding epitope for IL23p19., (© 2023. The Author(s).)

Proliferating cell nuclear antigen (PCNA) is a sliding clamp protein that coordinates DNA replication with various DNA maintenance events that are critical for human health. Recently, a hypomorphic homozygous serine to isoleucine (S228I) substitution in PCNA was described to underlie a rare DNA repair disorder known as PCNA-associated DNA repair disorder (PARD). PARD symptoms range from UV sensitivity, neurodegeneration, telangiectasia, and premature aging. We, and others, previously showed that the S228I variant changes the protein-binding pocket of PCNA to a conformation that impairs interactions with specific partners. Here, we report a second PCNA substitution (C148S) that also causes PARD. Unlike PCNA-S228I, PCNA-C148S has WT-like structure and affinity toward partners. In contrast, both disease-associated variants possess a thermostability defect. Furthermore, patient-derived cells homozygous for the C148S allele exhibit low levels of chromatin-bound PCNA and display temperature-dependent phenotypes. The stability defect of both PARD variants indicates that PCNA levels are likely an important driver of PARD disease. These results significantly advance our understanding of PARD and will likely stimulate additional work focused on clinical, diagnostic, and therapeutic aspects of this severe disease., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Effectively and accurately predicting the effects of interactions between proteins after amino acid mutations is a key issue for understanding the mechanism of protein function and drug design. In this study, we present a deep graph convolution (DGC) network-based framework, DGCddG, to predict the changes of protein-protein binding affinity after mutation. DGCddG incorporates multi-layer graph convolution to extract a deep, contextualized representation for each residue of the protein complex structure. The mined channels of the mutation sites by DGC is then fitted to the binding affinity with a multi-layer perceptron. Experiments with results on multiple datasets show that our model can achieve relatively good performance for both single and multi-point mutations. For blind tests on datasets related to angiotensin-converting enzyme 2 binding with the SARS-CoV-2 virus, our method shows better results in predicting ACE2 changes, may help in finding favorable antibodies. Code and data availability: https://github.com/lennylv/DGCddG.

TBR1 is a neuron-specific transcription factor involved in brain development and implicated in a neurodevelopmental disorder (NDD) combining features of autism spectrum disorder (ASD), intellectual disability (ID) and speech delay. TBR1 has been previously shown to interact with a small number of transcription factors and co-factors also involved in NDDs (including CASK, FOXP1/2/4 and BCL11A), suggesting that the wider TBR1 interactome may have a significant bearing on normal and abnormal brain development. Here, we have identified approximately 250 putative TBR1-interaction partners by affinity purification coupled to mass spectrometry. As well as known TBR1-interactors such as CASK, the identified partners include transcription factors and chromatin modifiers, along with ASD- and ID-related proteins. Five interaction candidates were independently validated using bioluminescence resonance energy transfer assays. We went on to test the interaction of these candidates with TBR1 protein variants implicated in cases of NDD. The assays uncovered disturbed interactions for NDD-associated variants and identified two distinct protein-binding domains of TBR1 that have essential roles in protein-protein interaction., (© The Author(s) 2022. Published by Oxford University Press.)

The transcription factor p63 shares a high sequence identity with the tumour suppressor p53 which manifests itself in high structural similarity and preference for DNA sequences. Mutations in the DNA binding domain (DBD) of p53 have been studied in great detail, enabling a general mechanism-based classification. In this study we provide a detailed investigation of all currently known mutations in the p63 DBD, which are associated with developmental syndromes, by measuring their impact on transcriptional activity, DNA binding affinity, zinc binding capacity and thermodynamic stability. Some of the mutations we have further characterized with respect to their ability to convert human dermal fibroblasts into induced keratinocytes. Here we propose a classification of the p63 DBD mutations based on the four different mechanisms of DNA binding impairment which we identified: direct DNA contact, zinc finger region, H2 region, and dimer interface mutations. The data also demonstrate that, in contrast to p53 cancer mutations, no p63 mutation induces global unfolding and subsequent aggregation of the domain. The dimer interface mutations that affect the DNA binding affinity by disturbing the interaction between the individual DBDs retain partial DNA binding capacity which correlates with a milder patient phenotype., (© 2023. The Author(s).)

Replication protein A (RPA) is a eukaryotic single-stranded (ss) DNA-binding (SSB) protein that is essential for all aspects of genome maintenance. RPA binds ssDNA with high affinity but can also diffuse along ssDNA. By itself, RPA is capable of transiently disrupting short regions of duplex DNA by diffusing from a ssDNA that flanks the duplex DNA. Using single-molecule total internal reflection fluorescence and optical trapping combined with fluorescence approaches, we show that S. cerevisiae Pif1 can use its ATP-dependent 5' to 3' translocase activity to chemomechanically push a single human RPA (hRPA) heterotrimer directionally along ssDNA at rates comparable to those of Pif1 translocation alone. We further show that using its translocation activity, Pif1 can push hRPA from a ssDNA loading site into a duplex DNA causing stable disruption of at least 9 bp of duplex DNA. These results highlight the dynamic nature of hRPA enabling it to be readily reorganized even when bound tightly to ssDNA and demonstrate a mechanism by which directional DNA unwinding can be achieved through the combined action of a ssDNA translocase that pushes an SSB protein. These results highlight the two basic requirements for any processive DNA helicase: transient DNA base pair melting (supplied by hRPA) and ATP-dependent directional ssDNA translocation (supplied by Pif1) and that these functions can be unlinked by using two separate proteins.

It has been recently suggested that a significant fraction of homomer protein-protein interfaces evolve neutrally, without contributing to function, due to a hydrophobic bias in missense mutations. However, the fraction of such gratuitous complexes is currently unknown. Here, we quantified the fraction of homodimers where multimerization is unlikely to contribute to their biochemical function. We show that: 1) ligand binding-site structure predicts whether a homomer is functional or not; the vast majority of homodimers with multichain binding-sites (MBS) are likely to be functional, while in homodimers with single-chain binding-sites (SBS) and small to medium interfaces, quaternary structure is unlikely to be functional in a significant fraction-35%, even up to 42%-of complexes; 2) the hydrophobicity of interfaces changes little with the strength of selection, and the amino acid composition of interfaces is shaped by the "hydrophobic ratchet" in both types, but they are not in a strict equilibrium with mutations; particularly cysteines are much more abundant in mutations than in interfaces or surfaces; 3) in MBS homomers, the interfaces are conserved, while in a high fraction of SBS homomers, the interface is not more conserved than the solvent-accessible surface; and 4) MBS homomer interfaces coevolve more strongly with ligand binding sites than the interfaces of SBS homomers, and MBS complexes have higher capacity to transfer information from ligands across the interfaces than SBS homomers, explaining the enrichment of allostery in the former., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

We report the discovery of the androgen receptor missense mutation V770D, that was found in two sisters suffering from complete androgen insensitivity. Experimental validation of AR V770 variants demonstrated that AR V770D was transcriptionally inactive due to the inability to dimerize and a reduced ligand binding affinity. The more conservative AR V770A variant showed a dimerization defect at low levels of DHT with a partial recovery of the transcriptional activity and of the receptor's ability to dimerize when increasing the DHT levels. With V770 located outside of the proposed LBD dimerization interface of the AR LBD homodimer crystal structure, the effects of the V770A mutation on AR dimerization were unexpected. We therefore explored whether the AR LBD dimerization interface would be better described by an alternative dimerization mode based on available human homodimeric LBD crystal structures of other nuclear receptors. Superposition of the monomeric AR LBD in the homodimeric crystal structures of GR, PR, ER, CAR, TRβ, and HNF-4α showed that the GR-like LBD dimer model was energetically the most stable. Moreover, V770 was a key energy residue in the GR-like LBD dimer while it was not involved in the stabilization of the AR LBD homodimer according to the crystal structure. Additionally, the observation that 4 AIS mutations impacted the stability of the AR LBD dimer while 16 mutations affected the GR-like LBD dimer, suggested that the AR LBD dimer crystal is a snapshot of a breathing AR LBD homodimer that can transition into the GR-like LBD dimer model., (© 2023 The Protein Society.)

Tight regulation of microtubule (MT) dynamics is necessary for proper spindle assembly and chromosome segregation. The MT destabilizing Kinesin-8, Kif18B, controls astral MT dynamics and spindle positioning. Kif18B interacts with importin α/β as well as with the plus-tip tracking protein EB1, but how these associations modulate Kif18B is not known. We mapped the key binding sites on Kif18B, made residue-specific mutations, and assessed their impact on Kif18B function. Blocking EB1 interaction disrupted Kif18B MT plus-end accumulation and inhibited its ability to control MT length on monopolar spindles in cells. Blocking importin α/β interaction disrupted Kif18B localization without affecting aster size. In vitro, importin α/β increased Kif18B MT association by increasing the on-rate and decreasing the off-rate from MTs, which stimulated MT destabilization. In contrast, EB1 promoted MT destabilization without increasing lattice binding in vitro, which suggests that EB1 and importin α/β have distinct roles in the regulation of Kif18B-mediated MT destabilization. We propose that importin α/β spatially modulate Kif18B association with MTs to facilitate its MT destabilization activity. Our results suggest that Ran regulation is important not only to control molecular motor function near chromatin but also to provide a spatial control mechanism to modulate MT binding of nuclear localization signal-containing spindle assembly factors.

The glycolipid transfer protein (GLTP) has been linked to many cellular processes aside from its best-known in vitro function as a lipid transport protein. It has been proposed to act as a sensor and regulator of glycosphingolipid homeostasis in cells. Furthermore, through its previously determined interaction with the endoplasmic reticulum membrane protein VAP-A (vesicle-associated membrane protein-associated protein A), GLTP may also be involved in facilitating vesicular transport in cells. In this study, we characterized the phenotype of CRISPR/Cas9-mediated GLTP KO HeLa cells. We showed that motility, three-dimensional growth, and cellular metabolism were all altered by GLTP knockout. Expression of a GLTP mutant incapable of binding VAP disrupted cell spheroid formation, indicating that the GLTP-VAP interaction is linked to cellular adhesion, cohesion, and three-dimensional growth. Most notably, we found evidence that GLTP, through its interaction with VAP-A, affects vesicular trafficking, marking the first cellular process discovered to be directly impacted by a change in GLTP expression., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

The interaction between human Growth Hormone (hGH) and hGH Receptor (hGHR) has basic relevance to cancer and growth disorders, and hGH is the scaffold for Pegvisomant, an anti-acromegaly therapeutic. For the latter reason, hGH has been extensively engineered by early workers to improve binding and other properties. We are particularly interested in E174 which belongs to the hGH zinc-binding triad; the substitution E174A is known to significantly increase binding, but to now no explanation has been offered. We generated this and several computationally-selected single-residue substitutions at the hGHR-binding site of hGH. We find that, while many successfully slow down dissociation of the hGH-hGHR complex once bound, they also slow down the association of hGH to hGHR. The E174A substitution induces a change in the Circular Dichroism spectrum that suggests the appearance of coiled-coiling. Here we show that E174A increases affinity of hGH against hGHR because the off-rate is slowed down more than the on-rate. For E174Y (and certain mutations at other sites) the slowdown in on-rate was greater than that of the off-rate, leading to decreased affinity. The results point to a link between structure, zinc binding, and hGHR-binding affinity in hGH., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Rajkovic et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

The sequence-based predictors of RNA-binding residues (RBRs) are trained on either structure-annotated or disorder-annotated binding regions. A recent study of predictors of protein-binding residues shows that they are plagued by high levels of cross-predictions (protein binding residues are predicted as nucleic acid binding) and that structure-trained predictors perform poorly for the disorder-annotated regions and vice versa. Consequently, we analyze a representative set of the structure and disorder trained predictors of RBRs to comprehensively assess quality of their predictions. Our empirical analysis that relies on a new and low-similarity benchmark dataset reveals that the structure-trained predictors of RBRs perform well for the structure-annotated proteins while the disorder-trained predictors provide accurate results for the disorder-annotated proteins. However, these methods work only modestly well on the opposite types of annotations, motivating the need for new solutions. Using an empirical approach, we design HybridRNAbind meta-model that generates accurate predictions and low amounts of cross-predictions when tested on data that combines structure and disorder-annotated RBRs. We release this meta-model as a convenient webserver which is available at https://www.csuligroup.com/hybridRNAbind/., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Somatic mutations are highly enriched at transcription factor (TF) binding sites, with the strongest trend being observed for ultraviolet light (UV)-induced mutations in melanomas. One of the main mechanisms proposed for this hypermutation pattern is the inefficient repair of UV lesions within TF-binding sites, caused by competition between TFs bound to these lesions and the DNA repair proteins that must recognize the lesions to initiate repair. However, TF binding to UV-irradiated DNA is poorly characterized, and it is unclear whether TFs maintain specificity for their DNA sites after UV exposure. We developed UV-Bind, a high-throughput approach to investigate the impact of UV irradiation on protein-DNA binding specificity. We applied UV-Bind to ten TFs from eight structural families, and found that UV lesions significantly altered the DNA-binding preferences of all the TFs tested. The main effect was a decrease in binding specificity, but the precise effects and their magnitude differ across factors. Importantly, we found that despite the overall reduction in DNA-binding specificity in the presence of UV lesions, TFs can still compete with repair proteins for lesion recognition, in a manner consistent with their specificity for UV-irradiated DNA. In addition, for a subset of TFs, we identified a surprising but reproducible effect at certain nonconsensus DNA sequences, where UV irradiation leads to a high increase in the level of TF binding. These changes in DNA-binding specificity after UV irradiation, at both consensus and nonconsensus sites, have important implications for the regulatory and mutagenic roles of TFs in the cell.

Background: Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has caused millions of infections and deaths worldwide since its discovery in late 2019 in Wuhan, China. The receptor-binding domain (RBD) of the SARS-CoV-2 spike protein binds to the human angiotensin-converting enzyme-2 (ACE2) receptor, a critical component of the renin-angiotensin system (RAS) that initiates the viral transmission. Most of the critical mutations found in SARS-CoV-2 are associated with the RBD of the spike protein. These mutations have the potential to reduce the efficacy of vaccines and neutralizing antibodies., Methods: In this review, the structural details of ACE2, RBD and their interactions are discussed. In addition, some critical mutations of RBD and their impact on ACE2-RBD interactions are also discussed., Conclusion: Preventing the interaction between Spike RBD and ACE2 is considered a viable therapeutic strategy since ACE2 binding by RBD is the first step in virus infection. Because the interactions between the two entities are critical for both viral transmission and therapeutic development, it is essential to understand their interactions in detail., (© 2022. The Author(s), under exclusive licence to Springer Nature B.V.)

Staphylococcal protein A (SpA) domain B (the basis of affibody) has been widely used in affinity chromatography and found therapeutic applications against inflammatory diseases through targeting the Fc part of immunoglobulin G (IgG). We have performed extensive molecular dynamics simulation of 41 SpA mutants and compared their dynamics and conformations to wild type. The simulations revealed the molecular details of structural and dynamics changes that occurred due to introducing point mutations and helped to explain the SPR results. It was observed in some variants a point mutation caused extensive structural changes far from the mutation site, while an effect of some other mutations was limited to the site of the mutated residue. Also, the pattern of hydrogen bond networks and hydrophobic core arrangements were investigated. We figured out mutations that occurred at positions 128, 136, 150 and 153, affected two hydrophobic cores at the interface as well as mutations introduced at positions 129 and 154 interrupted two hydrogen bond networks of the interface, SPR data showed all of these mutations reduced binding affinity significantly. Overall, by scanning the SpA-Fc interface through the large numbers of introduced mutations, the new insights have been gained which would help to design high- affinity ligands of IgG., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Computational prediction of the RBP bound sites using features learned from existing annotation knowledge is an effective method because high-throughput experiments are complex, expensive and time-consuming. Many methods have been proposed to predict RNA-protein binding sites. However, the partial information of RNA sequence is not fully used. In this study, we propose multiple convolutional neural networks (MCNN) method, which predicts RNA-protein binding sites by integrating multiple convolutional neural networks constructed by RNA sequence information extracted from windows with different lengths. First, MCNN trains multiple CNNs base on RNA sequences extracted by different window lengths. Second, MCNN can extract more binding patterns of RBPs by combining these trained multiple CNNs previously. Third, MCNN only uses RNA base sequence information for RNA-protein binding sites prediction, which extracts sequence binding features and predicts the result with same architecture. This avoids the information loss of feature extraction step. Our proposed MCNN demonstrates a competitive performance comparing with other methods on a large-scale dataset derived from CLIP-seq, which is an effective method for RNA-protein binding sites prediction. The source code of our proposed MCNN method can be found in https://github.com/biomg/MCNN.

The introduction of phosphorothioate (PS) linkages to the backbone of therapeutic nucleic acids substantially increases their stability and potency. It also affects their interactions with cellular proteins, but the molecular mechanisms that underlie this effect are poorly understood. Here, we report structural and biochemical studies of interactions between annexin A2, a protein that does not possess any known canonical DNA binding domains, and phosphorothioate-modified antisense oligonucleotides. We show that a unique mode of hydrophobic interactions between a sulfur atom of the phosphorothioate group and lysine and arginine residues account for the enhanced affinity of modified nucleic acid for the protein. Our results demonstrate that this mechanism of interaction is observed not only for nucleic acid-binding proteins but can also account for the association of PS oligonucleotides with other proteins. Using the anomalous diffraction of sulfur, we showed that preference for phosphorothioate stereoisomers is determined by the hydrophobic environment around the PS linkage that comes not only from protein but also from additional structural features within the ASO such as 5-Me groups on cytosine nucleobases., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The high magnitude zoonotic event has caused by Severe Acute Respitarory Syndrome CoronaVirus-2 (SARS-CoV-2) is Coronavirus Disease-2019 (COVID-19) epidemics. This disease has high rate of spreading than mortality in humans. The human receptor, Angiotensin-Converting Enzyme 2 (ACE2), is the leading target site for viral Spike-protein (S-protein) that function as binding ligands and are responsible for their entry in humans. The patients infected with COVID-19 with comorbidities, particularly cancer patients, have a severe effect or high mortality rate because of the suppressed immune system. Nevertheless, there might be a chance wherein cancer patients cannot be infected with SARS-CoV-2 because of mutations in the ACE2, which may be resistant to the spillover between species. This study aimed to determine the mutations in the sequence of the human ACE2 protein and its dissociation with SARS-CoV-2 that might be rejecting viral transmission. The in silico approaches were performed to identify the impact of SARS-CoV-2 S-protein with ACE2 mutations, validated experimentally, occurred in the patient, and reported in cell lines. The identified changes significantly affect SARS-CoV-2 S-protein interaction with ACE2, demonstrating the reduction in the binding affinity compared to SARS-CoV. The data presented in this study suggest ACE2 mutants have a higher and lower affinity with SARS-Cov-2 S-protein to the wild-type human ACE2 receptor. This study would likely be used to report SARS-CoV-2 resistant ACE2 mutations and can be used to design active peptide development to inactivate the viral spread of SARS-CoV-2 in humans., (© 2023. The Author(s).)

ZBTB7A belongs to a small family of transcription factors having three members in humans (7A, 7B, and 7C). They share a BTB/POZ protein interaction domain at the amino end and a zinc-finger DNA-binding domain at the carboxyl end. They control the transcription of a wide range of genes, having varied functions in hematopoiesis, oncogenesis, and metabolism (in particular glycolysis). ZBTB7A-binding profiles at gene promoters contain a consensus G(a/c)CCC motif, followed by a CCCC sequence in some instances. Structural and mutational investigations suggest that DNA-specific contacts with the four-finger tandem array of ZBTB7A are formed sequentially, initiated from ZF1-ZF2 binding to G(a/c)CCC before spreading to ZF3-ZF4, which bind the DNA backbone and the 3' CCCC sequence, respectively. Here, we studied some mutations found in t(8;21)-positive acute myeloid leukemia patients that occur within the ZBTB7A DNA-binding domain. We determined that these mutations generally impair ZBTB7A DNA binding, with the most severe disruptions resulting from mutations in ZF1 and ZF2, and the least from a frameshift mutation in ZF3 that results in partial mislocalization. Information provided here on ZBTB7A-DNA interactions is likely applicable to ZBTB7B/C, which have overlapping functions with ZBTB7A in controlling primary metabolism., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interests with the contents of this article., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

DNA methylation is critical for regulating gene expression, necessitating its accurate placement by enzymes such as the DNA methyltransferase DNMT3A. Dysregulation of this process is known to cause aberrant development and oncogenesis, yet how DNMT3A is regulated holistically by its three domains remains challenging to study. Here, we integrate base editing with a DNA methylation reporter to perform in situ mutational scanning of DNMT3A in cells. We identify mutations throughout the protein that perturb function, including ones at an interdomain interface that block allosteric activation. Unexpectedly, we also find mutations in the PWWP domain, a histone reader, that modulate enzyme activity despite preserving histone recognition and protein stability. These effects arise from altered PWWP domain DNA affinity, which we show is a noncanonical function required for full activity in cells. Our findings highlight mechanisms of interdomain crosstalk and demonstrate a generalizable strategy to probe sequence-activity relationships of nonessential chromatin regulators., (© 2022. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Zinc (Zn) is an essential trace element; it serves as a cofactor for a great number of enzymes, transcription factors, receptors, and other proteins. Zinc is also an important signaling molecule, which can be released from intracellular stores into the cytosol or extracellular space, for example, during synaptic transmission. Amongst cellular effects of zinc is activation of Kv7 (KCNQ, M-type) voltage-gated potassium channels. Here, we investigated relationships between Kv7 channel inhibition by Ca 2+ /calmodulin (CaM) and zinc-mediated potentiation. We show that Zn 2+ ionophore, zinc pyrithione (ZnPy), can prevent or reverse Ca 2+ /CaM-mediated inhibition of Kv7.2. In the presence of both Ca 2+ and Zn 2+ , the Kv7.2 channels lose most of their voltage dependence and lock in an open state. In addition, we demonstrate that mutations that interfere with CaM binding to Kv7.2 and Kv7.3 reduced channel membrane abundance and activity, but these mutants retained zinc sensitivity. Moreover, the relative efficacy of ZnPy to activate these mutants was generally greater, compared with the WT channels. Finally, we show that zinc sensitivity was retained in Kv7.2 channels assembled with mutant CaM with all four EF hands disabled, suggesting that it is unlikely to be mediated by CaM. Taken together, our findings indicate that zinc is a potent Kv7 stabilizer, which may protect these channels from physiological inhibitory effects of neurotransmitters and neuromodulators, protecting neurons from overactivity., Competing Interests: Conflict of interest The authors declare that they have no conflicts of interest with the contents of this article., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Human Replication Protein A (hRPA) is a multidomain protein that interacts with ssDNA intermediates to provide the latter much-needed stability during DNA metabolism and maintain genomic integrity. Although the ssDNA organization with hRPA was studied recently through experimental means, characterizing the underlying mechanism at the atomic level remains challenging because of the dynamic domain architecture of hRPA and poorly understood heterogeneity of ssDNA-protein interactions. Here, we used a computational framework, precisely tailored to capture protein-ssDNA interactions, and investigated the binding of hRPA with a 60 nt ssDNA. Two distinct binding mechanisms are realized based on the hRPA domain flexibility. For a rigid domain architecture of hRPA, ssDNA binds sequentially with hRPA domains, resulting in slow association kinetics. The binding pathway involves the formation of stable and distinct intermediate states. On contrary, for a flexible domain architecture of hRPA, ssDNA binds synergistically to the A and B domains followed by the rest of hRPA. The domain dynamics in hRPA alleviates the free energy cost of domain orientation necessary for specific binding with ssDNA, leading to fast association kinetics along a downhill binding free energy landscape. An ensemble of free energetically degenerate intermediate states is encountered that makes it arduous to characterize them structurally. An excellent match between our results with the available experimental observations provides new insights into the rich dynamics of hRPA binding to ssDNA and in general paves the way to investigate intricate details of ssDNA-protein interactions, crucial for cellular functioning., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Sangeeta, Bhattacherjee. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Motivation: Cell-type-specific gene expression is maintained in large part by transcription factors (TFs) selectively binding to distinct sets of sites in different cell types. Recent research works have provided evidence that such cell-type-specific binding is determined by TF's intrinsic sequence preferences, cooperative interactions with co-factors, cell-type-specific chromatin landscapes and 3D chromatin interactions. However, computational prediction and characterization of cell-type-specific and shared binding sites is rarely studied., Results: In this article, we propose two computational approaches for predicting and characterizing cell-type-specific and shared binding sites by integrating multiple types of features, in which one is based on XGBoost and another is based on convolutional neural network (CNN). To validate the performance of our proposed approaches, ChIP-seq datasets of 10 binding factors were collected from the GM12878 (lymphoblastoid) and K562 (erythroleukemic) human hematopoietic cell lines, each of which was further categorized into cell-type-specific (GM12878- and K562-specific) and shared binding sites. Then, multiple types of features for these binding sites were integrated to train the XGBoost- and CNN-based models. Experimental results show that our proposed approaches significantly outperform other competing methods on three classification tasks. Moreover, we identified independent feature contributions for cell-type-specific and shared sites through SHAP values and explored the ability of the CNN-based model to predict cell-type-specific and shared binding sites by excluding or including DNase signals. Furthermore, we investigated the generalization ability of our proposed approaches to different binding factors in the same cellular environment., Availability and Implementation: The source code is available at: https://github.com/turningpoint1988/CSSBS., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2022. Published by Oxford University Press.)

Background: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causes the coronavirus disease 2019 (COVID-19), which first appeared in December 2019. Angiotensin I converting enzyme 2 ( ACE2 ) receptor, present on the host cells, interacts with the receptor binding domain (RBD) of spike (S) protein of SARS-CoV-2 and facilitates the viral entry into host cells., Methods: Non-synonymous single nucleotide polymorphisms (nsSNPs) in the ACE2 gene may have an impact on the protein's stability and its function. The deleterious or harmful nsSNPs of the ACE2 gene that can change the strength as well as the pattern of interaction with the RBD of S protein were selected for this study., Results: The ACE2 :RBD interactions were analyzed by protein-protein docking study. The missense mutations A242V, R708W, G405E, D292N, Y633C, F308L, and G405E in ACE2 receptor were found to interact with RBD of Omicron subvariants with stronger binding affinity. Among the other selected nsSNPs of human ACE2 (h ACE2 ), R768W, Y654S, F588S, R710C, R710C, A191P, and R710C were found to have lower binding affinity for RBD of Omicron subvariants., Conclusion: The findings of this study suggest that the nsSNPs present in the human ACE2 gene alter the structure and function of the protein and, consequently, the susceptibility to Omicron subvariants., (Copyright© Bentham Science Publishers; For any queries, please email at epub@benthamscience.net.)

Long QT syndrome (LQTS) is a human inherited heart condition that can cause life-threatening arrhythmia including sudden cardiac death. Mutations in the ubiquitous Ca 2+ -sensing protein calmodulin (CaM) are associated with LQTS, but the molecular mechanism by which these mutations lead to irregular heartbeats is not fully understood. Here, we use a multidisciplinary approach including protein biophysics, structural biology, confocal imaging, and patch-clamp electrophysiology to determine the effect of the disease-associated CaM mutation E140G on CaM structure and function. We present novel data showing that mutant-regulated CaMKIIδ kinase activity is impaired with a significant reduction in enzyme autophosphorylation rate. We report the first high-resolution crystal structure of a LQTS-associated CaM variant in complex with the CaMKIIδ peptide, which shows significant structural differences, compared to the WT complex. Furthermore, we demonstrate that the E140G mutation significantly disrupted Ca v 1.2 Ca 2+ /CaM-dependent inactivation, while cardiac ryanodine receptor (RyR2) activity remained unaffected. In addition, we show that the LQTS-associated mutation alters CaM's Ca 2+ -binding characteristics, secondary structure content, and interaction with key partners involved in excitation-contraction coupling (CaMKIIδ, Ca v 1.2, RyR2). In conclusion, LQTS-associated CaM mutation E140G severely impacts the structure-function relationship of CaM and its regulation of CaMKIIδ and Ca v 1.2. This provides a crucial insight into the molecular factors contributing to CaM-mediated arrhythmias with a central role for CaMKIIδ., Competing Interests: Conflict of interest The authors declare no competing or financial interests., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

Background: Rapidly growing genome-wide ChIP-seq data have provided unprecedented opportunities to explore transcription factor (TF) binding under various cellular conditions. Despite the rich resources, development of analytical methods for studying the interaction among TFs in gene regulation still lags behind., Results: In order to address cooperative TF binding and detect TF clusters with coordinative functions, we have developed novel computational methods based on clustering the sample paths of nonhomogeneous Poisson processes. Simulation studies demonstrated the capability of these methods to accurately detect TF clusters and uncover the hierarchy of TF interactions. A further application to the multiple-TF ChIP-seq data in mouse embryonic stem cells (ESCs) showed that our methods identified the cluster of core ESC regulators reported in the literature and provided new insights on functional implications of transcrisptional regulatory modules., Conclusions: Effective analytical tools are essential for studying protein-DNA relations. Information derived from this research will help us better understand the orchestration of transcription factors in gene regulation processes., (© 2022. The Author(s).)

Nanobodies have emerged as powerful protein-binding tools to uncover protein functions. Using functionalized protein binders, proteins of interest can be visualized, degraded, delocalized, or post-translationally modified in vivo. We recently reported the use of two short peptide tags, 10-aa 127D01 and 14-aa VHH05, and their corresponding nanobodies, Nb127D01 and NbVHH05, for both in vitro and in vivo studies in Drosophila. Here, we provide detailed protocols for nanobody production and for visualization of proteins of interest in either fixed or live samples. In addition, we include protocols for endogenous protein tagging using CRISPR-mediated genome engineering. © 2022 Wiley Periodicals LLC. Basic Protocol 1: Nanobody production in S2 cells Basic Protocol 2: Nanobody expression and purification in bacterial cells Basic Protocol 3: Immunostaining with nanobodies Basic Protocol 4: Immunoblotting with nanobodies Basic Protocol 5: Immunoprecipitation with nanobodies prepared from S2 cells Basic Protocol 6: Immunoprecipitation with nanobodies prepared from bacteria Basic Protocol 7: NbVHH05 and Nb127D01 used as chromobodies Basic Protocol 8: NanoTag trap as a method to alter protein localization Support Protocol: CRISPR-mediated tagging of endogenous genes with NanoTags., (© 2022 Wiley Periodicals LLC.)

Since its first appearance in April 2021, B.1.617.2, also termed variant Delta, catalyzed one major worldwide wave dominating the second year of coronavirus disease 2019 (COVID-19) pandemic. Despite its quick disappearance worldwide, the strong virulence caused by a few point mutations remains an unsolved problem largely. Along with the other two sublineages, the Delta variant harbors an accumulation of Spike protein mutations, including the previously identified L452R, E484Q, and the newly emerged T478K on its receptor binding domain (RBD). We used molecular dynamics (MD) simulations, in combination with free energy perturbation (FEP) calculations, to examine the effects of two combinative mutation sets, L452R + E484Q and L452R + T478K. Our dynamic trajectories reveal an enhancement in binding affinity between mutated RBD and the common receptor protein angiotensin converting enzyme 2 (ACE2) through a net increase in the buried molecular surface area of the binary complex. This enhanced binding, mediated through Gln493, sets the same stage for all three sublineages due to the presence of L452R mutation. The other mutation component, E484Q or T478K, was found to impact the RBD-ACE2 binding and help the variant to evade several monoclonal antibodies (mAbs) in a distinct manner. Especially for L452R + T478K, synergies between mutations are mediated through a complex residual and water interaction network and further enhance its binding to ACE2. Taking together, this study demonstrates that new variants of SARS-CoV-2 accomplish both "attack" (infection) and "defense" (antibody neutralization escape) with the same "polished sword" (mutated Spike RBD).

The identification of transcription factor (TF) target genes is central in biology. A popular approach is based on the location by pattern matching of potential cis-regulatory elements (CREs). During the last few years, tools integrating next-generation sequencing data have been developed to improve the performance of pattern matching. However, such tools have not yet been comprehensively evaluated in plants. Hence, we developed a new streamlined method aiming at predicting CREs and target genes of plant TFs in specific organs or conditions. Our approach implements a supervised machine learning strategy, which allows decision rule models to be learnt using TF ChIP-chip/seq experimental data. Different layers of genomic features were integrated in predictive models: the position on the gene, the DNA sequence conservation, the chromatin state and various CRE footprints. Among the tested features, the chromatin features were crucial for improving the accuracy of the method. Furthermore, we evaluated the transferability of predictive models across TFs, organs and species. Finally, we validated our method by correctly inferring the target genes of key TFs controlling metabolite biosynthesis at the organ level in Arabidopsis. We developed a tool-Wimtrap-to reproduce our approach in plant species and conditions/organs for which ChIP-chip/seq data are available. Wimtrap is a user-friendly R package that supports an R Shiny web interface and is provided with pre-built models that can be used to quickly get predictions of CREs and TF gene targets in different organs or conditions in Arabidopsis thaliana, Solanum lycopersicum, Oryza sativa and Zea mays., (© The Author(s) 2022. Published by Oxford University Press on behalf of Japanese Society of Plant Physiologists. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Determining transcriptional factor binding sites (TFBSs) is critical for understanding the molecular mechanisms regulating gene expression in different biological conditions. Biological assays designed to directly mapping TFBSs require large sample size and intensive resources. As an alternative, ATAC-seq assay is simple to conduct and provides genomic cleavage profiles that contain rich information for imputing TFBSs indirectly. Previous footprint-based tools are inheritably limited by the accuracy of their bias correction algorithms and the efficiency of their feature extraction models. Here we introduce TAMC (Transcriptional factor binding prediction from ATAC-seq profile at Motif-predicted binding sites using Convolutional neural networks), a deep-learning approach for predicting motif-centric TF binding activity from paired-end ATAC-seq data. TAMC does not require bias correction during signal processing. By leveraging a one-dimensional convolutional neural network (1D-CNN) model, TAMC make predictions based on both footprint and non-footprint features at binding sites for each TF and outperforms existing footprinting tools in TFBS prediction particularly for ATAC-seq data with limited sequencing depth., Competing Interests: The authors have declared that no competing interests exist.

Genetic alternations can occur at noncoding regions, but how they contribute to cancer pathogenesis is poorly understood. Here, we established a mutational landscape of cis-regulatory regions (CREs) in acute promyelocytic leukemia (APL) based on whole-genome sequencing analysis of paired tumor and germline samples from 24 patients and epigenetic profiling of 16 patients. Mutations occurring in CREs occur preferentially in active enhancers bound by the complex of master transcription factors in APL. Among significantly enriched mutated CREs, we found a recurrently mutated region located within the third intron of WT1, an essential regulator of normal and malignant hematopoiesis. Focusing on noncoding mutations within this WT1 intron, an analysis on 169 APL patients revealed that somatic mutations were clustered into a focal hotspot region, including one site identified as a germline polymorphism contributing to APL risk. Significantly decreased WT1 expression was observed in APL patients bearing somatic and/or germline noncoding WT1 variants. Furthermore, biallelic WT1 inactivation was recurrently found in APL patients with noncoding WT1 variants, which resulted in the complete loss of WT1. The high incidence of biallelic inactivation suggested the tumor suppressor activity of WT1 in APL. Mechanistically, noncoding WT1 variants disrupted MYB binding on chromatin and suppressed the enhancer activity and WT1 expression through destroying the chromatin looping formation. Our study highlights the important role of noncoding variants in the leukemogenesis of APL., (© 2022 by The American Society of Hematology.)

The lysine (K) tRNA synthetase C-terminal (KTSC) domain containing proteins are widely spread in Bacteria, Archaea and Viruses, but the function of this short domain is unclear. The occurrence of the fusion of KTSC domain to a catalytic domain or domains related to DNA or RNA metabolisms suggests its potential role in DNA or RNA binding. Here, we report the characterization of Mvu8s from Methanolobus vulcani, which consists of a single KTSC domain. Mvu8s binds specifically to ssDNA with an affinity approximately 40- and 10-fold higher than those for dsDNA and ssRNA in vitro, respectively. It shows a slight preference to the G-rich DNA sequence but barely binds the A-stretch. Crystal structure of Mvu8s shows that it forms a homo-tetramer, with each monomer composed of a four-strand antiparallel β-sheet and a helix-turn-helix in the order of β1-β2-β3-α1-α2-β4. Four basic residues (R3, R7, K54 and K58) were found to serve important roles in ssDNA-binding. And, the spiral arrangement of the DNA interfaces in Mvu8s homo-tetramer presumably results in ssDNA wrapping. Our results not only offer clues of the functions of the KTSC domain containing proteins but also expand our knowledge on the non-oligonucleotide-binding (OB) fold single-stranded DNA-binding proteins in Archaea., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)