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.

Fuzzy recognition by the prokaryotic transcription factor HigA2 from Vibrio cholerae.

Disordered protein sequences can exhibit different binding modes, ranging from well-ordered folding-upon-binding to highly dynamic fuzzy binding. The primary function of the intrinsically disordered region of the antitoxin HigA2 from Vibrio cholerae is to neutralize HigB2 toxin through ultra-high-affinity folding-upon-binding interaction. Here, we show that the same intrinsically disordered region can also mediate fuzzy interactions with its operator DNA and, through interplay with the folded helix-turn-helix domain, regulates transcription from the higBA2 operon. NMR, SAXS, ITC and in vivo experiments converge towards a consistent picture where a specific set of residues in the intrinsically disordered region mediate electrostatic and hydrophobic interactions while "hovering" over the DNA operator. Sensitivity of the intrinsically disordered region to scrambling the sequence, position-specific contacts and absence of redundant, multivalent interactions, point towards a more specific type of fuzzy binding. Our work demonstrates how a bacterial regulator achieves dual functionality by utilizing two distinct interaction modes within the same disordered sequence.

TripHLApan: predicting HLA molecules binding peptides based on triple coding matrix and transfer learning.

Human leukocyte antigen (HLA) recognizes foreign threats and triggers immune responses by presenting peptides to T cells. Computationally modeling the binding patterns between peptide and HLA is very important for the development of tumor vaccines. However, it is still a big challenge to accurately predict HLA molecules binding peptides. In this paper, we develop a new model TripHLApan for predicting HLA molecules binding peptides by integrating triple coding matrix, BiGRU + Attention models, and transfer learning strategy. We have found the main interaction site regions between HLA molecules and peptides, as well as the correlation between HLA encoding and binding motifs. Based on the discovery, we make the preprocessing and coding closer to the natural biological process. Besides, due to the input being based on multiple types of features and the attention module focused on the BiGRU hidden layer, TripHLApan has learned more sequence level binding information. The application of transfer learning strategies ensures the accuracy of prediction results under special lengths (peptides in length 8) and model scalability with the data explosion. Compared with the current optimal models, TripHLApan exhibits strong predictive performance in various prediction environments with different positive and negative sample ratios. In addition, we validate the superiority and scalability of TripHLApan's predictive performance using additional latest data sets, ablation experiments and binding reconstitution ability in the samples of a melanoma patient. The results show that TripHLApan is a powerful tool for predicting the binding of HLA-I and HLA-II molecular peptides for the synthesis of tumor vaccines. TripHLApan is publicly available at https://github.com/CSUBioGroup/TripHLApan.git.

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.

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.

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.

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.

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.

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.

Real-time MBDi-RPA using methyl-CpG binding protein 2: A real-time detection method for simple and rapid estimation of CpG methylation status.

BACKGROUND: Analysis of CpG methylation is informative for cancer diagnosis. Previously, we developed a novel method to discriminate CpG methylation status in target DNA by blocking recombinase polymerase amplification (RPA), an isothermal DNA amplification technique, using methyl-CpG binding domain (MBD) protein 2 (MBD2). The method was named MBD protein interference-RPA (MBDi-RPA). In this study, MBDi-RPA was performed using methyl-CpG binding protein 2 (MeCP2), another MBD family protein, as the blocking agent. RESULTS: MBDi-RPA using MeCP2 detected low levels of CpG methylation, showing that it had higher sensitivity than MBDi-RPA using MBD2. We also developed real-time RPA, which enabled rapid analysis of DNA amplification without the need for laborious agarose gel electrophoresis and used it in combination with MBDi-RPA. We termed this method real-time MBDi-RPA. The method using MeCP2 could determine the abundance ratio of CpG-methylated target DNA simply and rapidly, although highly sensitive detection was challenging. SIGNIFICANCE AND NOVELTY: Real-time MBDi-RPA using MeCP2 could be potentially useful for estimating CpG methylation status in target DNA prior to more detailed analyses.

A new paradigm for applying deep learning to protein-ligand interaction prediction.

Protein-ligand interaction prediction presents a significant challenge in drug design. Numerous machine learning and deep learning (DL) models have been developed to accurately identify docking poses of ligands and active compounds against specific targets. However, current models often suffer from inadequate accuracy or lack practical physical significance in their scoring systems. In this research paper, we introduce IGModel, a novel approach that utilizes the geometric information of protein-ligand complexes as input for predicting the root mean square deviation of docking poses and the binding strength (pKd, the negative value of the logarithm of binding affinity) within the same prediction framework. This ensures that the output scores carry intuitive meaning. We extensively evaluate the performance of IGModel on various docking power test sets, including the CASF-2016 benchmark, PDBbind-CrossDocked-Core and DISCO set, consistently achieving state-of-the-art accuracies. Furthermore, we assess IGModel's generalizability and robustness by evaluating it on unbiased test sets and sets containing target structures generated by AlphaFold2. The exceptional performance of IGModel on these sets demonstrates its efficacy. Additionally, we visualize the latent space of protein-ligand interactions encoded by IGModel and conduct interpretability analysis, providing valuable insights. This study presents a novel framework for DL-based prediction of protein-ligand interactions, contributing to the advancement of this field. The IGModel is available at GitHub repository https://github.com/zchwang/IGModel.

Molecular dynamics simulations to explore the binding mode between the amyloid-ß protein precursor (APP) and adaptor protein Mint2.

Alzheimer's disease (AD) presents a significant challenge in neurodegenerative disease management, with limited therapeutic options available for its prevention and treatment. At the heart of AD pathogenesis is the amyloid-ß (Aß) protein precursor (APP), with the interaction between APP and the adaptor protein Mint2 being crucial. Despite previous explorations into the APP-Mint2 interaction, the dynamic regulatory mechanisms by which Mint2 modulates APP binding remain poorly understood. This study undertakes molecular dynamics simulations across four distinct systems-free Mint2, Mint2 bound to APP, a mutant form of Mint2, and the mutant form bound to APP-over an extensive 400 ns timeframe. Our findings reveal that the mutant Mint2 experiences significant secondary structural transformations, notably the formation of an α-helix in residues S55-K65 upon APP binding, within the 400 ns simulation period. Additionally, we observed a reduction in the active pocket size of the mutant Mint2 compared to its wild-type counterpart, enhancing its APP binding affinity. These insights hold promise for guiding the development of novel inhibitors targeting the Mints family, potentially paving the way for new therapeutic strategies in AD prevention and treatment.

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.