Dynamics, Conformational Entropy, and Frustration in Protein-Protein Interactions Involving an Intrinsically Disordered Protein Domain

Frustraevo: a web server to localize and quantify the conservation of local energetic frustration in protein families

According to the Principle of Minimal Frustration, folded proteins can only have a minimal number of strong energetic conflicts in their native states. However, not all interactions are energetically optimized for folding but some remain in energetic conflict, i.e. they are highly frustrated. This remaining local energetic frustration has been shown to be statistically correlated with distinct functional aspects such as protein-protein interaction sites, allosterism and catalysis. Fuelled by the recent breakthroughs in efficient protein structure prediction that have made available good quality models for most proteins, we have developed a strategy to calculate local energetic frustration within large protein families and quantify its conservation over evolutionary time. Based on this evolutionary information we can identify how stability and functional constraints have appeared at the common ancestor of the family and have been maintained over the course of evolution. Here, we present FrustraEvo, a web server tool to calculate and quantify the conservation of local energetic frustration in protein families.

Impact of mutations on the stability of SARS-CoV-2 nucleocapsid protein structure

The nucleocapsid (N) protein of SARS-CoV-2 is known to participate in various host cellular processes, including interferon inhibition, RNA interference, apoptosis, and regulation of virus life cycles. Additionally, it has potential as a diagnostic antigen and/or immunogen. Our research focuses on examining structural changes caused by mutations in the N protein. We have modeled the complete tertiary structure of native and mutated forms of the N protein using Alphafold2. Notably, the N protein contains 3 disordered regions. The focus was on investigating the impact of mutations on the stability of the protein's dimeric structure based on binding free energy calculations (MM-PB/GB-SA) and RMSD fluctuations after MD simulations. The results demonstrated that 28 mutations out of 37 selected mutations analyzed, compared with wild-type N protein, resulted in a stable dimeric structure, while 9 mutations led to destabilization. Our results are important to understand the tertiary structure of the N protein dimer of SARS-CoV-2 and the effect of mutations on it, their behavior in the host cell, as well as for the research of other viruses belonging to the same genus additionally, to anticipate potential strategies for addressing this viral illness․.

Local energetic frustration conservation in protein families and superfamilies

Energetic local frustration offers a biophysical perspective to interpret the effects of sequence variability on protein families. Here we present a methodology to analyze local frustration patterns within protein families and superfamilies that allows us to uncover constraints related to stability and function, and identify differential frustration patterns in families with a common ancestry. We analyze these signals in very well studied protein families such as PDZ, SH3, ɑ and β globins and RAS families. Recent advances in protein structure prediction make it possible to analyze a vast majority of the protein space. An automatic and unsupervised proteome-wide analysis on the SARS-CoV-2 virus demonstrates the potential of our approach to enhance our understanding of the natural phenotypic diversity of protein families beyond single protein instances. We apply our method to modify biophysical properties of natural proteins based on their family properties, as well as perform unsupervised analysis of large datasets to shed light on the physicochemical signatures of poorly characterized proteins such as the ones belonging to emergent pathogens.

Conformational entropy in molecular recognition of intrinsically disordered proteins

Broad conformational ensembles make intrinsically disordered proteins or regions entropically intriguing. Although methodologically challenging and understudied, emerging studies into their changes in conformational entropy (ΔS°conf) upon complex formation have provided both quantitative and qualitative insight. Recent work based on thermodynamics from isothermal titration calorimetry and NMR spectroscopy uncovers an expanded repertoire of regulatory mechanisms, where ΔS°conf plays roles in partner selection, state behavior, functional buffering, allosteric regulation, and drug design. We highlight these mechanisms to display the large entropic reservoir of IDPs for the regulation of molecular communication. We call upon the field to make efforts to contribute to this insight as more studies are needed for forwarding mechanistic decoding of intrinsically disordered proteins and their complexes.

Residue-Frustration-Based Prediction of Protein-Protein Interactions Using Machine Learning

The study of protein-protein interactions (PPIs) is important in understanding the function of proteins. However, it is still a challenge to investigate the transient protein-protein interaction by experiments. Hence, the computational prediction for protein-protein interactions draws growing attention. Statistics-based features have been widely used in the studies of protein structure prediction and protein folding. Due to the scarcity of experimental data of PPI, it is difficult to construct a conventional statistical feature for PPI prediction, and the application of statistics-based features is very limited in this field. In this paper, we explored the application of frustration, a statistical potential, in PPI prediction. By comparing the energetic contribution of the extra stabilization energy from a given residue pair in the native protein with the statistics of the energies, we obtained the residue pair's frustration index. By calculating the number of residue pairs with a high frustration index, the highly frustrated density, a residue-frustration-based feature, was then obtained to describe the tendency of residues to be involved in PPI. Highly frustrated density, as well as structure-based features, were then used to describe protein residues and combined with the long short-term memory (LSTM) neural network to predict PPI residue pairs. Our model correctly predicted 75% dimers when only the top 2‰ residue pairs were selected in each dimer. Our model, which considers the statistics-based features, is significantly different from the models based on the chemical features of residues. We found that frustration can effectively describe the tendency of residue to be involved in PPI. Frustration-based features can replace chemical features to combine with machine learning and realize the better performance of PPI prediction. It reveals the great potential of statistical potential such as frustration in PPI prediction.

FrustratometeR: an R-package to compute local frustration in protein structures, point mutants and MD simulations

SUMMARY

Once folded, natural protein molecules have few energetic conflicts within their polypeptide chains. Many protein structures do however contain regions where energetic conflicts remain after folding, i.e. they are highly frustrated. These regions, kept in place over evolutionary and physiological timescales, are related to several functional aspects of natural proteins such as protein-protein interactions, small ligand recognition, catalytic sites and allostery. Here, we present FrustratometeR, an R package that easily computes local energetic frustration on a personal computer or a cluster. This package facilitates large scale analysis of local frustration, point mutants and molecular dynamics (MD) trajectories, allowing straightforward integration of local frustration analysis into pipelines for protein structural analysis.

AVAILABILITY AND IMPLEMENTATION

https://github.com/proteinphysiologylab/frustratometeR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Exploring the Conformational Landscape of the Neh4 and Neh5 Domains of Nrf2 Using Two Different Force Fields and Circular Dichroism

The nuclear factor erythroid 2-related factor 2 (Nrf2)-ARE transcriptional response pathway plays a critical role in protecting the cell from oxidative stresses via the upregulation of cytoprotective genes. Aberrant activation of Nrf2 in cancer cells can confer this cytoprotectivity, thereby reducing the efficacy of both chemotherapeutics and radiotherapies. Key to this antioxidant pathway is the interaction between Nrf2 and CREB binding protein (CBP), mediated by the Neh4 and Neh5 domains of Nrf2. Disruption of this interaction via small-molecule therapeutics could negate the effects of aberrant Nrf2 upregulation. Due to the disordered nature of these domains, there remains no three-dimensional structure of Neh4 or Neh5, making structure-based drug design a challenge. Here, we performed 48 μs of unbiased molecular dynamics (MD) simulations with the Amber99SB*-ILDNP and CHARMM36m force fields and circular dichroism (CD) spectropolarimetry experiments to elucidate the free-state structures of these domains; no previous data regarding their conformational landscapes exists. There are two main findings: First, we find Neh5 to be markedly more disordered than Neh4, which has nine residues in the middle of the domain showing α-helical propensity, thus pointing to Neh4 and Neh5 having different binding mechanisms. Second, the two force fields show strong differences for the glutamic acid-rich Neh5 peptide but are in reasonable agreement for Neh4, which has no glutamic acid. The CHARMM36m force field agrees more closely with the CD results.

Chemical control of peptide material phase transitions

Progressive solute-rich polymer phase transitions provide pathways for achieving ordered supramolecular assemblies. Intrinsically disordered protein domains specifically regulate information in biological networks via conformational ordering. Here we consider a molecular tagging strategy to control ordering transitions in polymeric materials and provide a proof-of-principle minimal peptide phase network captured with a dynamic chemical network.

Substrate initiated assembly of a dynamic chemical network.

Structural insights into the formation of oligomeric state by a type I Hsp40 chaperone

Molecular chaperones can prevent and repair protein misfolding and aggregation to maintain protein homeostasis in cells. Hsp40 chaperones interact with unfolded client proteins via the dynamic multivalent interaction (DMI) mechanism with their multiple client-binding sites. Here we report that a type I Hsp40 chaperone from Streptococcus pneumonia (spHsp40) forms a concentration-independent polydispersity oligomer state in solution. The crystal structure of spHsp40 determined at 2.75 Å revealed that each monomer has a type I Hsp40 structural fold containing a zinc finger domain and C-terminal domains I and II (CTD I and CTD II). Subsequent quaternary structure analysis using a PISA server generated two dimeric models. The interface mutational analysis suggests the conserved C-terminal dimeric motif as a basis for dimer formation and that the novel dimeric interaction between a client-binding site in CTD I and the zinc finger domain promotes the formation of the spHsp40 oligomeric state. In vitro functional analysis demonstrated that spHsp40 oligomer is fully active and possess the optimal activity in stimulating the ATPase activity of spHsp70. The oligomer state of type I Hsp40 and its formation might be important in understanding Hsp40 function and its interaction with client proteins.

Protein tolerance to random circular permutation correlates with thermostability and local energetics of residue-residue contacts

Adenylate kinase (AK) orthologs with a range of thermostabilities were subjected to random circular permutation, and deep mutational scanning was used to evaluate where new protein termini were nondisruptive to activity. The fraction of circularly permuted variants that retained function in each library correlated with AK thermostability. In addition, analysis of the positional tolerance to new termini, which increase local conformational flexibility, showed that bonds were either functionally sensitive to cleavage across all homologs, differentially sensitive, or uniformly tolerant. The mobile AMP-binding domain, which displays the highest calculated contact energies, presented the greatest tolerance to new termini across all AKs. In contrast, retention of function in the lid and core domains was more dependent upon AK melting temperature. These results show that family permutation profiling identifies primary structure that has been selected by evolution for dynamics that are critical to activity within an enzyme family. These findings also illustrate how deep mutational scanning can be applied to protein homologs in parallel to differentiate how topology, stability, and local energetics govern mutational tolerance.

Characterization of the dynamics and the conformational entropy in the binding between TAZ1 and CTAD-HIF-1α

The interaction between the C-terminal transactivation domain of HIF-1α (CTAD-HIF-1α) and the transcriptional adapter zinc binding 1 (TAZ1) domain of CREB binding protein participate in the initiation of gene transcription during hypoxia. Unbound CTAD-HIF-1α is disordered but undergoes a disorder-to-order transition upon binding to TAZ1. We have here performed NMR side chain and backbone relaxation studies on TAZ1 and side chain relaxation measurements on CTAD-HIF-1α in order to investigate the role of picosecond to nanosecond dynamics. We find that the internal motions are significantly affected upon binding, both on the side chain and the backbone level. The dynamic response corresponds to a conformational entropy change that contributes substantially to the binding thermodynamics for both binding partners. Furthermore, the conformational entropy change for the well-folded TAZ1 varies upon binding to different IDP targets. We further identify a cluster consisting of side chains in bound TAZ1 and CTAD-HIF-1α that experience extensive dynamics and are part of the binding region that involves the N-terminal end of the LPQL motif in CTAD-HIF-1α; a feature that might have an important role in the termination of the hypoxic response.

Features of molecular recognition of intrinsically disordered proteins via coupled folding and binding

The sequence-structure-function paradigm of proteins has been revolutionized by the discovery of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). In contrast to traditional ordered proteins, IDPs/IDRs are unstructured under physiological conditions. The absence of well-defined three-dimensional structures in the free state of IDPs/IDRs is fundamental to their function. Folding upon binding is an important mode of molecular recognition for IDPs/IDRs. While great efforts have been devoted to investigating the complex structures and binding kinetics and affinities, our knowledge on the binding mechanisms of IDPs/IDRs remains very limited. Here, we review recent advances on the binding mechanisms of IDPs/IDRs. The structures and kinetic parameters of IDPs/IDRs can vary greatly, and the binding mechanisms can be highly dependent on the structural properties of IDPs/IDRs. IDPs/IDRs can employ various combinations of conformational selection and induced fit in a binding process, which can be templated by the target and/or encoded by the IDP/IDR. Further studies should provide deeper insights into the molecular recognition of IDPs/IDRs and enable the rational design of IDP/IDR binding mechanisms in the future.

Intrinsically Disordered Transactivation Domains Bind to TAZ1 Domain of CBP via Diverse Mechanisms

CREB-binding protein is a multidomain transcriptional coactivator whose transcriptional adaptor zinc-binding 1 (TAZ1) domain mediates interactions with a number of intrinsically disordered transactivation domains (TADs), including the CREB-binding protein/p300-interacting transactivator with ED-rich tail, the hypoxia inducible factor 1α, p53, the signal transducer and activator of transcription 2, and the NF-κB p65 subunit. These five disordered TADs undergo partial disorder-to-order transitions upon binding TAZ1, forming fuzzy complexes with helical segments. Interestingly, they wrap around TAZ1 with different orientations and occupy the binding sites with various orders. To elucidate the microscopic molecular details of the binding processes of TADs with TAZ1, in this work, we carried out extensive molecular dynamics simulations using a coarse-grained topology-based model. After careful calibration of the models to reproduce the residual helical contents and binding affinities, our simulations were able to recapitulate the experimentally observed flexibility profiles. Although great differences exist in the complex structures, we found similarities between hypoxia inducible factor 1α and signal transducer and activator of transcription 2 as well as between CREB-binding protein/p300-interacting transactivator with ED-rich tail and NF-κB p65 subunit in the binding kinetics and binding thermodynamics. Although the origins of similarities and differences in the binding mechanisms remain unclear, our results provide some clues that indicate that binding of TADs to TAZ1 could be templated by the target as well as encoded by the TADs.

Dynamics of the intrinsically disordered protein NUPR1 in isolation and in its fuzzy complexes with DNA and prothymosin α

Intrinsically disordered proteins (IDPs) explore diverse conformations in their free states and, a few of them, also in their molecular complexes. This functional plasticity is essential for the function of IDPs, although their dynamics in both free and bound states is poorly understood. NUPR1 is a protumoral multifunctional IDP, activated during the acute phases of pancreatitis. It interacts with DNA and other IDPs, such as prothymosin α (ProTα), with dissociation constants of ~0.5 μM, and a 1:1 stoichiometry. We studied the structure and picosecond-to-nanosecond (ps-ns) dynamics by using both NMR and SAXS in: (i) isolated NUPR1; (ii) the NUPR1/ProTα complex; and (iii) the NUPR1/double stranded (ds) GGGCGCGCCC complex. Our SAXS findings show that NUPR1 remained disordered when bound to either partner, adopting a worm-like conformation; the fuzziness of bound NUPR1 was also pinpointed by NMR. Residues with the largest values of the relaxation rates (R₁, R_1ρ, R₂ and η_xy), in the free and bound species, were mainly clustered around the 30s region of the sequence, which agree with one of the protein hot-spots already identified by site-directed mutagenesis. Not only residues in this region had larger relaxation rates, but they also moved slower than the rest of the molecule, as indicated by the reduced spectral density approach (RSDA). Upon binding, the energy landscape of NUPR1 was not funneled down to a specific, well-folded conformation, but rather its backbone flexibility was kept, with distinct motions occurring at the hot-spot region.

Entropy and Information within Intrinsically Disordered Protein Regions

Bioinformatics and biophysical studies of intrinsically disordered proteins and regions (IDRs) note the high entropy at individual sequence positions and in conformations sampled in solution. This prevents application of the canonical sequence-structure-function paradigm to IDRs and motivates the development of new methods to extract information from IDR sequences. We argue that the information in IDR sequences cannot be fully revealed through positional conservation, which largely measures stable structural contacts and interaction motifs. Instead, considerations of evolutionary conservation of molecular features can reveal the full extent of information in IDRs. Experimental quantification of the large conformational entropy of IDRs is challenging but can be approximated through the extent of conformational sampling measured by a combination of NMR spectroscopy and lower-resolution structural biology techniques, which can be further interpreted with simulations. Conformational entropy and other biophysical features can be modulated by post-translational modifications that provide functional advantages to IDRs by tuning their energy landscapes and enabling a variety of functional interactions and modes of regulation. The diverse mosaic of functional states of IDRs and their conformational features within complexes demands novel metrics of information, which will reflect the complicated sequence-conformational ensemble-function relationship of IDRs.

Role of Conformational Entropy in Molecular Recognition by TAZ1 of CBP

The globular transcriptional adapter zinc binding 1 (TAZ1) domain of CREB binding protein participates in protein-protein interactions that are involved in transcriptional regulation. TAZ1 binds numerous targets, of which many are intrinsically disordered proteins that undergo a disorder-to-order transition to various degrees. One such target is the disordered transactivation domain of transcription factor RelA (TAD-RelA), which with its interaction with TAZ1 is involved in transcriptional regulation of genes in NF-κB signaling. We have here performed nuclear magnetic resonance backbone and side-chain relaxation studies to investigate the influence of RelA-TA2 (residues 425-490 in TAD-RelA) binding on the subnanosecond internal motions of TAZ1. We find a considerable dynamic response on both the backbone and side-chain levels, which corresponds to a conformational entropy change that contributes significantly to the binding energetics. We further show that the microscopic origins of the dynamic response of TAZ1 vary depending on the target. This study demonstrates that folded protein domains that are able to interact with various targets are not dynamically passive but can have a significant role in the motional response upon target association.

Role of Backbone Dynamics in Modulating the Interactions of Disordered Ligands with the TAZ1 Domain of the CREB-Binding Protein

The intrinsically disordered transactivation domains of HIF-1α and CITED2 compete for binding of the TAZ1 domain of the CREB-binding protein by a unidirectional allosteric mechanism involving direct competition for shared binding sites, ternary complex formation, and TAZ1 conformational changes. To gain insight into the mechanism by which CITED2 displaces HIF-1α from TAZ1, we used nuclear magnetic resonance spin relaxation methods to obtain an atomic-level description of the picosecond to nanosecond backbone dynamics that contribute to TAZ1 binding and competition. We show that HIF-1α and CITED2 adopt different dynamics in their complexes with TAZ1, with flexibility observed for HIF-1α in regions that would maintain accessibility for CITED2 to bind to TAZ1 and facilitate subsequent HIF-1α dissociation. In contrast, critical regions of CITED2 adopt a rigid structure in its complex with TAZ1, minimizing the ability of HIF-1α to compete for binding. We also find that TAZ1, previously thought to be a rigid scaffold for binding of disordered protein ligands, displays altered backbone dynamics in its various bound states. TAZ1 is more rigid in its CITED2-bound state than in its free state or in complex with HIF-1α, with increased rigidity observed not only in the CITED2 binding site but also in regions of TAZ1 that undergo conformational changes between the HIF-1α- and CITED2-bound structures. Taken together, these data suggest that backbone dynamics in TAZ1, as well as in the HIF-1α and CITED2 ligands, play a role in modulating the occupancy of TAZ1 and highlight the importance of characterizing both binding partners in molecular interactions.

Elucidating binding mechanisms and dynamics of intrinsically disordered protein complexes using NMR spectroscopy

Advances in characterizing complexes of intrinsically disordered proteins (IDPs) have led to the discovery of a remarkably diverse interaction landscape that includes folding-upon-binding, highly dynamic complexes, multivalent interactions as well as regulatory switches controlled by post-translational modifications. Nuclear magnetic resonance (NMR) spectroscopy has in recent years made significant contributions to this field by describing the binding mechanisms and mapping conformational dynamics on multiple time scales. Importantly, this progress has been associated with specific methodological developments in NMR, for example in exchange techniques, allowing challenging biological systems to be studied at atomic resolution. In general, the level of dynamics observed in IDP complexes does not correlate with binding affinities, demonstrating the intricate relationship between conformational dynamics and IDP regulatory function.

Structure and dynamics conspire in the evolution of affinity between intrinsically disordered proteins

In every established species, protein-protein interactions have evolved such that they are fit for purpose. However, the molecular details of the evolution of new protein-protein interactions are poorly understood. We have used nuclear magnetic resonance spectroscopy to investigate the changes in structure and dynamics during the evolution of a protein-protein interaction involving the intrinsically disordered CREBBP (CREB-binding protein) interaction domain (CID) and nuclear coactivator binding domain (NCBD) from the transcriptional coregulators NCOA (nuclear receptor coactivator) and CREBBP/p300, respectively. The most ancient low-affinity "Cambrian-like" [540 to 600 million years (Ma) ago] CID/NCBD complex contained less secondary structure and was more dynamic than the complexes from an evolutionarily younger "Ordovician-Silurian" fish ancestor (ca. 440 Ma ago) and extant human. The most ancient Cambrian-like CID/NCBD complex lacked one helix and several interdomain interactions, resulting in a larger solvent-accessible surface area. Furthermore, the most ancient complex had a high degree of millisecond-to-microsecond dynamics distributed along the entire sequences of both CID and NCBD. These motions were reduced in the Ordovician-Silurian CID/NCBD complex and further redistributed in the extant human CID/NCBD complex. Isothermal calorimetry experiments show that complex formation is enthalpically favorable and that affinity is modulated by a largely unfavorable entropic contribution to binding. Our data demonstrate how changes in structure and motion conspire to shape affinity during the evolution of a protein-protein complex and provide direct evidence for the role of structural, dynamic, and frustrational plasticity in the evolution of interactions between intrinsically disordered proteins.