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.

Pathogenic mutations of human phosphorylation sites affect protein-protein interactions.

Despite their lack of a defined 3D structure, intrinsically disordered regions (IDRs) of proteins play important biological roles. Many IDRs contain short linear motifs (SLiMs) that mediate protein-protein interactions (PPIs), which can be regulated by post-translational modifications like phosphorylation. 20% of pathogenic missense mutations are found in IDRs, and understanding how such mutations affect PPIs is essential for unraveling disease mechanisms. Here, we employ peptide-based interaction proteomics to investigate 36 disease-associated mutations affecting phosphorylation sites. Our results unveil significant differences in interactomes between phosphorylated and non-phosphorylated peptides, often due to disrupted phosphorylation-dependent SLiMs. We focused on a mutation of a serine phosphorylation site in the transcription factor GATAD1, which causes dilated cardiomyopathy. We find that this phosphorylation site mediates interaction with 14-3-3 family proteins. Follow-up experiments reveal the structural basis of this interaction and suggest that 14-3-3 binding affects GATAD1 nucleocytoplasmic transport by masking a nuclear localisation signal. Our results demonstrate that pathogenic mutations of human phosphorylation sites can significantly impact protein-protein interactions, offering insights into potential molecular mechanisms underlying pathogenesis.

A unified analysis of evolutionary and population constraint in protein domains highlights structural features and pathogenic sites.

Protein evolution is constrained by structure and function, creating patterns in residue conservation that are routinely exploited to predict structure and other features. Similar constraints should affect variation across individuals, but it is only with the growth of human population sequencing that this has been tested at scale. Now, human population constraint has established applications in pathogenicity prediction, but it has not yet been explored for structural inference. Here, we map 2.4 million population variants to 5885 protein families and quantify residue-level constraint with a new Missense Enrichment Score (MES). Analysis of 61,214 structures from the PDB spanning 3661 families shows that missense depleted sites are enriched in buried residues or those involved in small-molecule or protein binding. MES is complementary to evolutionary conservation and a combined analysis allows a new classification of residues according to a conservation plane. This approach finds functional residues that are evolutionarily diverse, which can be related to specificity, as well as family-wide conserved sites that are critical for folding or function. We also find a possible contrast between lethal and non-lethal pathogenic sites, and a surprising clinical variant hot spot at a subset of missense enriched positions.

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.

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.

Perspectives on Protein-Nanoparticle Interactions at the In Vivo Level.

The distinct features of nanoparticles have provided a vast opportunity of developing new diagnosis and therapy strategies for miscellaneous diseases. Although a few nanomedicines are available in the market or in the translation stage, many important issues are still unsolved. When entering the body, nanomaterials will be quickly coated by proteins from their surroundings, forming a corona on their surface, the so-called protein corona. Studies have shown that the protein corona has many important biological implications, particularly at the in vivo level. For example, they can promote the immune system to rapidly clear these outer materials and prevent nanoparticles from playing their designed role in therapy. In this Perspective, the available techniques for characterizing protein-nanoparticle interactions are critically summarized. Effects of nanoparticle properties and environmental factors on protein corona formation, which can further regulate the in vivo fate of nanoparticles, are highlighted and discussed. Moreover, recent progress on the biomedical application of protein corona-engineered nanoparticles is introduced, and future directions for this important yet challenging research area are also briefly discussed.

De novo design of drug-binding proteins with predictable binding energy and specificity.

The de novo design of small molecule-binding proteins has seen exciting recent progress; however, high-affinity binding and tunable specificity typically require laborious screening and optimization after computational design. We developed a computational procedure to design a protein that recognizes a common pharmacophore in a series of poly(ADP-ribose) polymerase-1 inhibitors. One of three designed proteins bound different inhibitors with affinities ranging from <5 nM to low micromolar. X-ray crystal structures confirmed the accuracy of the designed protein-drug interactions. Molecular dynamics simulations informed the role of water in binding. Binding free energy calculations performed directly on the designed models were in excellent agreement with the experimentally measured affinities. We conclude that de novo design of high-affinity small molecule-binding proteins with tuned interaction energies is feasible entirely from computation.

Considerations Around Structure-Based Drug Discovery for KRAS Using DOCK.

Molecular docking is a popular computational tool in drug discovery. Leveraging structural information, docking software predicts binding poses of small molecules to cavities on the surfaces of proteins. Virtual screening for ligand discovery is a useful application of docking software. In this chapter, using the enigmatic KRAS protein as an example system, we endeavor to teach the reader about best practices for performing molecular docking with UCSF DOCK. We discuss methods for virtual screening and docking molecules on KRAS. We present the following six points to optimize our docking setup for prosecuting a virtual screen: protein structure choice, pocket selection, optimization of the scoring function, modification of sampling spheres and sampling procedures, choosing an appropriate portion of chemical space to dock, and the choice of which top scoring molecules to pick for purchase.

1D and 2D NMR for KRAS:Ligand Binding.

Fragment-based screening by ligand-observed 1D NMR and binding interface mapping by protein-observed 2D NMR are popular methods used in drug discovery. These methods allow researchers to detect compound binding over a wide range of affinities and offer a simultaneous assessment of solubility, purity, and chemical formula accuracy of the target compounds and the 15N-labeled protein when examined by 1D and 2D NMR, respectively. These methods can be applied for screening fragment binding to the active (GMPPNP-bound) and inactive (GDP-bound) states of oncogenic KRAS mutants.

Profiling Complex RAS-Effector Interactions Using NMR Spectroscopy.

Knowledge of how effectors interact with RAS GTPases is key to understanding how these switch-like proteins function in cells. Effectors bind specifically to GTP-loaded RAS using RAS association (RA) or RAS binding domains (RBDs) that show wide-ranging affinities and thermodynamic characteristics. Both normal development and RAS-induced tumorigenesis depend on multiple distinct effector proteins that are frequently co-expressed and co-localized, suggesting an antagonistic nature to signaling whereby multiple proteins compete for a limited pool of activated GTPase. NMR spectroscopy offers a powerful approach to multiplex effectors and/or regulatory enzymes and quantifies their interaction with RAS, expanding our biophysical and systems-level understanding of RAS signaling in a more integrated and physiologically relevant setting. Here we describe a method to directly quantitate GTPase binding to competing effectors, using wild-type KRAS complex with ARAF and PLCε1 as a model. Unlabeled RBD/RA domains are added simultaneously to isotopically labeled RAS, and peak intensities at chemical shifts characteristic of individually bound domains provide quantitation. Similar competition-based assays can be run with small molecule interactors, GEF/GAP domains, or regulatory enzymes that drive posttranslational modifications. Such efforts bring in vitro interaction experiments in line with more complex cellular environments.

Viral Receptor-Binding Protein Evolves New Function through Mutations That Cause Trimer Instability and Functional Heterogeneity.

When proteins evolve new activity, a concomitant decrease in stability is often observed because the mutations that confer new activity can destabilize the native fold. In the conventional model of protein evolution, reduced stability is considered a purely deleterious cost of molecular innovation because unstable proteins are prone to aggregation and are sensitive to environmental stressors. However, recent work has revealed that nonnative, often unstable protein conformations play an important role in mediating evolutionary transitions, raising the question of whether instability can itself potentiate the evolution of new activity. We explored this question in a bacteriophage receptor-binding protein during host-range evolution. We studied the properties of the receptor-binding protein of bacteriophage λ before and after host-range evolution and demonstrated that the evolved protein is relatively unstable and may exist in multiple conformations with unique receptor preferences. Through a combination of structural modeling and in vitro oligomeric state analysis, we found that the instability arises from mutations that interfere with trimer formation. This study raises the intriguing possibility that protein instability might play a previously unrecognized role in mediating host-range expansions in viruses.

Engineering an artificial catch bond using mechanical anisotropy.

Catch bonds are a rare class of protein-protein interactions where the bond lifetime increases under an external pulling force. Here, we report how modification of anchor geometry generates catch bonding behavior for the mechanostable Dockerin G:Cohesin E (DocG:CohE) adhesion complex found on human gut bacteria. Using AFM single-molecule force spectroscopy in combination with bioorthogonal click chemistry, we mechanically dissociate the complex using five precisely controlled anchor geometries. When tension is applied between residue #13 on CohE and the N-terminus of DocG, the complex behaves as a two-state catch bond, while in all other tested pulling geometries, including the native configuration, it behaves as a slip bond. We use a kinetic Monte Carlo model with experimentally derived parameters to simulate rupture force and lifetime distributions, achieving strong agreement with experiments. Single-molecule FRET measurements further demonstrate that the complex does not exhibit dual binding mode behavior at equilibrium but unbinds along multiple pathways under force. Together, these results show how mechanical anisotropy and anchor point selection can be used to engineer artificial catch bonds.

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.

Molecular docking and molecular dynamics simulation decoding molecular mechanism of EDCs binding to hERRγ.

CONTEXT: Human estrogen-related receptor γ (hERRγ) is a key protein involved in various endocrines and metabolic signaling. Numerous environmental endocrine-disrupting chemicals (EDCs) can impact related physiological activities through receptor signaling pathways. Focused on hERRγ with 4-isopropylphenol, bisphenol-F (BPF), and BP(2,2)(Un) complexes, we executed molecular docking and multiple molecular dynamics (MD) simulations along with molecular mechanics/Poisson-Boltzmann surface area (MM-PBSA) and solvation interaction energy (SIE) calculation to study the detailed dynamical structural characteristics and interactions between them. Molecular docking showed that hydrogen bonds and hydrophobic interactions were the prime interactions to keep the stability of BPF-hERRγ and hERRγ-BP(2,2)(Un) complexes. Through MD simulations, we observed that all complexes reach equilibrium during the initial 50 ns of simulation, but these three EDCs lead to local structure changes in hERRγ. Energy results further identified key residues L268, V313, L345, and F435 around the binding pockets through CH-π, π-π, and hydrogen bonds interactions play an important stabilizing role in the recognition with EDCs. And most noticeable of all, hydrophobic methoxide groups in BP(2,2)(Un) is useful for decreasing the binding ability between EDCs and hERRγ. These results may contribute to evaluate latent diseases associated with EDCs exposure at the micro level and find potential substitutes. METHOD: Autodock4.2 was used to conduct the molecular docking, sietraj program was performed to calculate the energy, and VMD software was used to visualize the structure. Amber18 was conducted to perform the MD simulation and other analyses.

In Vitro Determination of Temperature-Dependent DNA Binding of the Evening Complex Using Electrophoretic Mobility Shift Assays.

Electrophoretic mobility shift assays (EMSAs) of DNA-binding proteins and labeled DNA allow the qualitative and quantitative characterization of protein-DNA complex formation using native (nondenaturing) polyacrylamide or agarose gel electrophoresis. By varying the incubation temperature of the protein-DNA binding reaction and maintaining this temperature during electrophoresis, temperature-dependent protein-DNA interactions can be investigated. Here, we provide examples of the binding of a transcriptional repressor complex called the Evening Complex, comprising the DNA-binding protein LUX ARRYTHMO (LUX), the scaffold protein EARLY FLOWERING 3 (ELF3), and the adapter protein ELF4, to its cognate DNA and demonstrate direct detection and visualization of thermoresponsive binding in vitro. As negative controls we use the LUX DNA-binding domain and LUX full length protein, which do not exhibit temperature-dependent DNA binding.

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

Bioisoteres for carboxylic acids: From ionized isosteres to novel unionized replacements.

Carboxylic acids are key pharmacophoric elements in many molecules. They can be seen as a problem by some, due to perceived permeability challenges, potential for high plasma protein binding and the risk of forming reactive metabolites due to acyl-glucuronidation. By others they are viewed more favorably as they can decrease lipophilicity by adding an ionizable center which can be beneficial for solubility, and can add enthalpic interactions with the target protein. However, there are many instances where the replacement of a carboxylic acid with a bioisosteric group is required. This has led to the development of a number of ionizable groups which sufficiently mimic the carboxylic acid functionality whilst improving, for example, the metabolic profile of the molecule in question. An alternative strategy involves replacement of the carboxylate by neutral functional groups. This review initially details carefully selected examples whereby tetrazoles, acyl sulfonamides or isoxazolols have been beneficially utilized as carboxylic acid bioisosteres altering physicohemical properties, interactions with the target and metabolism and/or pharmacokinetics, before delving further into the binding mode of carboxylic acid derivatives with their target proteins. This analysis highlights new ways to consider the replacement of carboxylic acids by neutral bioisosteric groups which either rely on hydrogen bonds or cation-π interactions. It should serve as a useful guide for scientists working in drug discovery.

Structure-based discovery of small molecule inhibitors of FKBP51-Hsp90 protein-protein interaction.

The heat shock protein 90 kDa (Hsp90) molecular chaperone machinery is responsible for the folding and activation of hundreds of important clients such as kinases, steroid hormone receptors, transcription factors, etc. This process is dynamically regulated in an ATP-dependent manner by Hsp90 co-chaperones including a group of tetratricopeptide (TPR) motif proteins that bind to the C-terminus of Hsp90. Among these TPR containing co-chaperones, FK506-binding protein 51 kDa (FKBP51) is reported to play an important role in stress-related pathologies, psychiatric disorders, Alzheimer's disease, and cancer, making FKBP51-Hsp90 interaction a potential therapeutic target. In this study, we report identification of potent and selective inhibitors of FKBP51-Hsp90 protein-protein interaction using a structure-based virtual screening approach. Upon in vitro evaluation, the identified hits show a considerable degree of selectivity towards FKBP51 over other TPR proteins, particularly for highly homologous FKBP52. Tyr355 of FKBP51 emerged as an important contributor to inhibitor's specificity. Additionally, we demonstrate the impact of these inhibitors on cellular energy metabolism, and neurite outgrowth, which are subjects of FKBP51 regulation. Overall, the results from this study highlight a novel pharmacological approach towards regulation of FKBP51 function and more generally, Hsp90 function via its interaction with TPR co-chaperones.

Caenorhabditis elegans telomere-binding proteins TEBP-1 and TEBP-2 adapt the Myb module to dimerize and bind telomeric DNA.

Protecting chromosome ends from misrecognition as double-stranded (ds) DNA breaks is fundamental to eukaryotic viability. The protein complex shelterin prevents a DNA damage response at mammalian telomeres. Mammalian shelterin proteins TRF1 and TRF2 and their homologs in yeast and protozoa protect telomeric dsDNA. N-terminal homodimerization and C-terminal Myb-domain-mediated dsDNA binding are two structural hallmarks of end protection by TRF homologs. Yet our understanding of how Caenorhabditis elegans protects its telomeric dsDNA is limited. Recently identified C. elegans proteins TEBP-1 (also called DTN-1) and TEBP-2 (also called DTN-2) are functional homologs of TRF proteins, but how they bind DNA and whether or how they dimerize is not known. TEBP-1 and TEBP-2 harbor three Myb-containing domains (MCDs) and no obvious dimerization domain. We demonstrate biochemically that only the third MCD binds DNA. We solve the X-ray crystal structure of TEBP-2 MCD3 with telomeric dsDNA to reveal the structural mechanism of telomeric dsDNA protection in C. elegans. Mutagenesis of the DNA-binding site of TEBP-1 and TEBP-2 compromises DNA binding in vitro, and increases DNA damage signaling, lengthens telomeres, and decreases brood size in vivo. Via an X-ray crystal structure, biochemical validation of the dimerization interface, and SEC-MALS analysis, we demonstrate that MCD1 and MCD2 form a composite dimerization module that facilitates not only TEBP-1 and TEBP-2 homodimerization but also heterodimerization. These findings provide fundamental insights into C. elegans telomeric dsDNA protection and highlight how different eukaryotes have evolved distinct strategies to solve the chromosome end protection problem.

A hydrophobic groove in secretagogin allows for alternate interactions with SNAP-25 and syntaxin-4 in endocrine tissues.

Vesicular release of neurotransmitters and hormones relies on the dynamic assembly of the exocytosis/trans-SNARE complex through sequential interactions of synaptobrevins, syntaxins, and SNAP-25. Despite SNARE-mediated release being fundamental for intercellular communication in all excitable tissues, the role of auxiliary proteins modulating the import of reserve vesicles to the active zone, and thus, scaling repetitive exocytosis remains less explored. Secretagogin is a Ca2+-sensor protein with SNAP-25 being its only known interacting partner. SNAP-25 anchors readily releasable vesicles within the active zone, thus being instrumental for 1st phase release. However, genetic deletion of secretagogin impedes 2nd phase release instead, calling for the existence of alternative protein-protein interactions. Here, we screened the secretagogin interactome in the brain and pancreas, and found syntaxin-4 grossly overrepresented. Ca2+-loaded secretagogin interacted with syntaxin-4 at nanomolar affinity and 1:1 stoichiometry. Crystal structures of the protein complexes revealed a hydrophobic groove in secretagogin for the binding of syntaxin-4. This groove was also used to bind SNAP-25. In mixtures of equimolar recombinant proteins, SNAP-25 was sequestered by secretagogin in competition with syntaxin-4. Kd differences suggested that secretagogin could shape unidirectional vesicle movement by sequential interactions, a hypothesis supported by in vitro biological data. This mechanism could facilitate the movement of transport vesicles toward release sites, particularly in the endocrine pancreas where secretagogin, SNAP-25, and syntaxin-4 coexist in both α- and ß-cells. Thus, secretagogin could modulate the pace and fidelity of vesicular hormone release by differential protein interactions.