Investigating the Role of Large-Scale Domain Dynamics in Protein-Protein Interactions

A coarse-grained model for disordered and multi-domain proteins

Many proteins contain more than one folded domain, and such modular multi-domain proteins help expand the functional repertoire of proteins. Because of their larger size and often substantial dynamics, it may be difficult to characterize the conformational ensembles of multi-domain proteins by simulations. Here, we present a coarse-grained model for multi-domain proteins that is both fast and provides an accurate description of the global conformational properties in solution. We show that the accuracy of a one-bead-per-residue coarse-grained model depends on how the interaction sites in the folded domains are represented. Specifically, we find excessive domain-domain interactions if the interaction sites are located at the position of the Cα atoms. We also show that if the interaction sites are located at the center of mass of the residue, we obtain good agreement between simulations and experiments across a wide range of proteins. We then optimize our previously described CALVADOS model using this center-of-mass representation, and validate the resulting model using independent data. Finally, we use our revised model to simulate phase separation of both disordered and multi-domain proteins, and to examine how the stability of folded domains may differ between the dilute and dense phases. Our results provide a starting point for understanding interactions between folded and disordered regions in proteins, and how these regions affect the propensity of proteins to self-associate and undergo phase separation.

Molecular switching in transcription through splicing and proline-isomerization regulates stress responses in plants

The Arabidopsis thaliana DREB2A transcription factor interacts with the negative regulator RCD1 and the ACID domain of subunit 25 of the transcriptional co-regulator mediator (Med25) to integrate stress signals for gene expression, with elusive molecular interplay. Using biophysical and structural analyses together with high-throughput screening, we reveal a bivalent binding switch in DREB2A containing an ACID-binding motif (ABS) and the known RCD1-binding motif (RIM). The RIM is lacking in a stress-induced DREB2A splice variant with retained transcriptional activity. ABS and RIM bind to separate sites on Med25-ACID, and NMR analyses show a structurally heterogeneous complex deriving from a DREB2A-ABS proline residue populating cis- and trans-isomers with remote impact on the RIM. The cis-isomer stabilizes an α-helix, while the trans-isomer may introduce energetic frustration facilitating rapid exchange between activators and repressors. Thus, DREB2A uses a post-transcriptionally and post-translationally modulated switch for transcriptional regulation.

Specifying conformational heterogeneity of multi-domain proteins at atomic resolution

The conformational landscape of multi-domain proteins is inherently linked to their specific functions. This also holds for polyubiquitin chains that are assembled by two or more ubiquitin domains connected by a flexible linker thus showing a large interdomain mobility. However, molecular recognition and signal transduction are associated with particular conformational substates that are populated in solution. Here, we apply high-resolution NMR spectroscopy in combination with dual-scale MD simulations to explore the conformational space of K6-, K29-, and K33-linked diubiquitin molecules. The conformational ensembles are evaluated utilizing a paramagnetic cosolute reporting on solvent exposure plus a set of complementary NMR parameters. This approach unravels a conformational heterogeneity of diubiquitins and explains the diversity of structural models that have been determined for K6-, K29-, and K33-linked diubiquitins in free and ligand-bound states so far. We propose a general application of the approach developed here to demystify multi-domain proteins occurring in nature.

Decoding Protein Dynamics in Cells Using Chemical Cross-Linking and Hierarchical Analysis

Protein dynamics play a crucial role in their diverse functions. The intracellular environment significantly influences protein dynamics, particularly for intrinsically disordered proteins (IDPs). To comprehensively capture structural information from various proteins within cells and characterize protein dynamics, chemical cross-linking mass spectrometry was employed. In this study, we introduce a hierarchical decoding strategy that enables the investigation of protein dynamics in vivo. Computational analysis based on distance restraints derived from cross-links is used to infer protein dynamics in cells. To facilitate this analysis, we leverage the prior structure obtained from AlphaFold2. By employing this strategy, we can characterize the full-length structure of multi-domain proteins taking into account their distinct dynamic features. Furthermore, by combining restraint sampling with an unbiased sampling and evaluation approach, we can provide a comprehensive description of the intrinsic motion of IDPs. Consequently, the hierarchical strategy we propose holds significant potential in advancing our understanding of the molecular mechanisms that undelie protein functions in cells.

Conformation and Affinity Modulations by Multiple Phosphorylation Occurring in the BIN1 SH3 Domain Binding Site of the Tau Protein Proline-Rich Region

An increase in phosphorylation of the Tau protein is associated with Alzheimer's disease (AD) progression through unclear molecular mechanisms. In general, phosphorylation modifies the interaction of intrinsically disordered proteins, such as Tau, with other proteins; however, elucidating the structural basis of this regulation mechanism remains challenging. The bridging integrator-1 gene is an AD genetic determinant whose gene product, BIN1, directly interacts with Tau. The proline-rich motif recognized within a Tau(210-240) peptide by the SH3 domain of BIN1 (BIN1 SH3) is defined as ₂₁₆PTPP₂₁₉, and this interaction is modulated by phosphorylation. Phosphorylation of T217 within the Tau(210-240) peptide led to a 6-fold reduction in the affinity, while single phosphorylation at either T212, T231, or S235 had no effect on the interaction. Nonetheless, combined phosphorylation of T231 and S235 led to a 3-fold reduction in the affinity, although these phosphorylations are not within the BIN1 SH3-bound region of the Tau peptide. Using nuclear magnetic resonance (NMR) spectroscopy, these phosphorylations were shown to affect the local secondary structure and dynamics of the Tau(210-240) peptide. Models of the (un)phosphorylated peptides were obtained from molecular dynamics (MD) simulation validated by experimental data and showed compaction of the phosphorylated peptide due to increased salt bridge formation. This dynamic folding might indirectly impact the BIN1 SH3 binding by a decreased accessibility of the binding site. Regulation of the binding might thus not only be due to local electrostatic or steric effects from phosphorylation but also to the modification of the conformational properties of Tau.

NMR Provides Unique Insight into the Functional Dynamics and Interactions of Intrinsically Disordered Proteins

Intrinsically disordered proteins are ubiquitous throughout all known proteomes, playing essential roles in all aspects of cellular and extracellular biochemistry. To understand their function, it is necessary to determine their structural and dynamic behavior and to describe the physical chemistry of their interaction trajectories. Nuclear magnetic resonance is perfectly adapted to this task, providing ensemble averaged structural and dynamic parameters that report on each assigned resonance in the molecule, unveiling otherwise inaccessible insight into the reaction kinetics and thermodynamics that are essential for function. In this review, we describe recent applications of NMR-based approaches to understanding the conformational energy landscape, the nature and time scales of local and long-range dynamics and how they depend on the environment, even in the cell. Finally, we illustrate the ability of NMR to uncover the mechanistic basis of functional disordered molecular assemblies that are important for human health.

Improving Martini 3 for Disordered and Multidomain Proteins

Coarse-grained molecular dynamics simulations are a useful tool to determine conformational ensembles of proteins. Here, we show that the coarse-grained force field Martini 3 underestimates the global dimensions of intrinsically disordered proteins (IDPs) and multidomain proteins when compared with small-angle X-ray scattering (SAXS) data and that increasing the strength of protein-water interactions favors more expanded conformations. We find that increasing the strength of interactions between protein and water by ca. 10% results in improved agreement with the SAXS data for IDPs and multidomain proteins. We also show that this correction results in a more accurate description of self-association of IDPs and folded proteins and better agreement with paramagnetic relaxation enhancement data for most IDPs. While simulations with this revised force field still show deviations to experiments for some systems, our results suggest that it is overall a substantial improvement for coarse-grained simulations of soluble proteins.

Conformational ensembles of intrinsically disordered proteins and flexible multidomain proteins

Intrinsically disordered proteins (IDPs) and multidomain proteins with flexible linkers show a high level of structural heterogeneity and are best described by ensembles consisting of multiple conformations with associated thermodynamic weights. Determining conformational ensembles usually involves the integration of biophysical experiments and computational models. In this review, we discuss current approaches to determine conformational ensembles of IDPs and multidomain proteins, including the choice of biophysical experiments, computational models used to sample protein conformations, models to calculate experimental observables from protein structure, and methods to refine ensembles against experimental data. We also provide examples of recent applications of integrative conformational ensemble determination to study IDPs and multidomain proteins and suggest future directions for research in the field.

Single-molecule conformational dynamics of a transcription factor reveals a continuum of binding modes controlling association and dissociation

Binding and unbinding of transcription factors to DNA are kinetically controlled to regulate the transcriptional outcome. Control of the release of the transcription factor NF-κB from DNA is achieved through accelerated dissociation by the inhibitor protein IκBα. Using single-molecule FRET, we observed a continuum of conformations of NF-κB in free and DNA-bound states interconverting on the subseconds to minutes timescale, comparable to in vivo binding on the seconds timescale, suggesting that structural dynamics directly control binding kinetics. Much of the DNA-bound NF-κB is partially bound, allowing IκBα invasion to facilitate DNA dissociation. IκBα induces a locked conformation where the DNA-binding domains of NF-κB are too far apart to bind DNA, whereas a loss-of-function IκBα mutant retains the NF-κB conformational ensemble. Overall, our results suggest a novel mechanism with a continuum of binding modes for controlling association and dissociation of transcription factors.

Intrinsically Disordered Proteins as Regulators of Transient Biological Processes and as Untapped Drug Targets

Intrinsically disordered proteins (IDPs) are critical players in the dynamic control of diverse cellular processes, and provide potential new drug targets because their dysregulation is closely related to many diseases. This review focuses on several medicinal studies that have identified low-molecular-weight inhibitors of IDPs. In addition, clinically relevant liquid-liquid phase separations-which critically involve both intermolecular interactions between IDPs and their posttranslational modification-are analyzed to understand the potential of IDPs as new drug targets.

Molecular recognition of ubiquitin and Lys63-linked diubiquitin by STAM2 UIM-SH3 dual domain: the effect of its linker length and flexibility

Multidomain proteins represent a broad spectrum of the protein landscape and are involved in various interactions. They could be considered as modular building blocks assembled in distinct fashion and connected by linkers of varying lengths and sequences. Due to their intrinsic flexibility, these linkers provide proteins a subtle way to modulate interactions and explore a wide range of conformational space. In the present study, we are seeking to understand the effect of the flexibility and dynamics of the linker involved in the STAM2 UIM-SH3 dual domain protein with respect to molecular recognition. We have engineered several constructs of UIM-SH3 with different length linkers or domain deletion. By means of SAXS and NMR experiments, we have shown that the modification of the linker modifies the flexibility and the dynamics of UIM-SH3. Indeed, the global tumbling of both the UIM and SH3 domain is different but not independent from each other while the length of the linker has an impact on the ps-ns time scale dynamics of the respective domains. Finally, the modification of the flexibility and dynamics of the linker has a drastic effect on the interaction of UIM-SH3 with Lys63-linked diubiquitin with a roughly eight-time weaker dissociation constant.

HuR biological function involves RRM3-mediated dimerization and RNA binding by all three RRMs

HuR/ELAVL1 is an RNA-binding protein involved in differentiation and stress response that acts primarily by stabilizing messenger RNA (mRNA) targets. HuR comprises three RNA recognition motifs (RRMs) where the structure and RNA binding of RRM3 and of full-length HuR remain poorly understood. Here, we report crystal structures of RRM3 free and bound to cognate RNAs. Our structural, NMR and biochemical data show that RRM3 mediates canonical RNA interactions and reveal molecular details of a dimerization interface localized on the α-helical face of RRM3. NMR and SAXS analyses indicate that the three RRMs in full-length HuR are flexibly connected in the absence of RNA, while they adopt a more compact arrangement when bound to RNA. Based on these data and crystal structures of tandem RRM1,2-RNA and our RRM3-RNA complexes, we present a structural model of RNA recognition involving all three RRM domains of full-length HuR. Mutational analysis demonstrates that RRM3 dimerization and RNA binding is required for functional activity of full-length HuR in vitro and to regulate target mRNAs levels in human cells, thus providing a fine-tuning for HuR activity in vivo.

Quantum Molecular Dynamics, Topological, Group Theoretical and Graph Theoretical Studies of Protein-Protein Interactions

BACKGROUND

Protein-protein interactions (PPIs) are becoming increasingly important as PPIs form the basis of multiple aggregation-related diseases such as cancer, Creutzfeldt-Jakob, and Alzheimer's diseases. This mini-review presents hybrid quantum molecular dynamics, quantum chemical, topological, group theoretical, graph theoretical, and docking studies of PPIs. We also show how these theoretical studies facilitate the discovery of some PPI inhibitors of therapeutic importance.

OBJECTIVE

The objective of this review is to present hybrid quantum molecular dynamics, quantum chemical, topological, group theoretical, graph theoretical, and docking studies of PPIs. We also show how these theoretical studies enable the discovery of some PPI inhibitors of therapeutic importance.

METHODS

This article presents a detailed survey of hybrid quantum dynamics that combines classical and quantum MD for PPIs. The article also surveys various developments pertinent to topological, graph theoretical, group theoretical and docking studies of PPIs and highlight how the methods facilitate the discovery of some PPI inhibitors of therapeutic importance.

RESULTS

It is shown that it is important to include higher-level quantum chemical computations for accurate computations of free energies and electrostatics of PPIs and Drugs with PPIs, and thus techniques that combine classical MD tools with quantum MD are preferred choices. Topological, graph theoretical and group theoretical techniques are shown to be important in studying large network of PPIs comprised of over 100,000 proteins where quantum chemical and other techniques are not feasible. Hence, multiple techniques are needed for PPIs.

CONCLUSION

Drug discovery and our understanding of complex PPIs require multifaceted techniques that involve several disciplines such as quantum chemistry, topology, graph theory, knot theory and group theory, thus demonstrating a compelling need for a multi-disciplinary approach to the problem.

Recent advances in automated protein design and its future challenges

INTRODUCTION

Protein function is determined by protein structure which is in turn determined by the corresponding protein sequence. If the rules that cause a protein to adopt a particular structure are understood, it should be possible to refine or even redefine the function of a protein by working backwards from the desired structure to the sequence. Automated protein design attempts to calculate the effects of mutations computationally with the goal of more radical or complex transformations than are accessible by experimental techniques. Areas covered: The authors give a brief overview of the recent methodological advances in computer-aided protein design, showing how methodological choices affect final design and how automated protein design can be used to address problems considered beyond traditional protein engineering, including the creation of novel protein scaffolds for drug development. Also, the authors address specifically the future challenges in the development of automated protein design. Expert opinion: Automated protein design holds potential as a protein engineering technique, particularly in cases where screening by combinatorial mutagenesis is problematic. Considering solubility and immunogenicity issues, automated protein design is initially more likely to make an impact as a research tool for exploring basic biology in drug discovery than in the design of protein biologics.

A mechanism for lipid binding to apoE and the role of intrinsically disordered regions coupled to domain-domain interactions

Relative to the apolipoprotein E (apoE) E3 allele of the APOE gene, apoE4 strongly increases the risk for the development of late-onset Alzheimer's disease. However, apoE4 differs from apoE3 by only a single amino acid at position 112, which is arginine in apoE4 and cysteine in apoE3. It remains unclear why apoE3 and apoE4 are functionally different. Described here is a proposal for understanding the functional differences between these two isoforms with respect to lipid binding. A mechanism is proposed that is based on the full-length monomeric structure of the protein, on hydrogen-deuterium exchange mass spectrometry data, and on the role of intrinsically disordered regions to control protein motions. It is proposed that lipid binds between the N-terminal and C-terminal domains and that separation of the two domains, along with the presence of intrinsically disordered regions, controls this process. The mechanism explains why apoE3 differs from apoE4 with respect to different lipid-binding specificities, why lipid increases the binding of apoE to its receptor, and why specific residues are conserved.

Human EGF-derived direct and reverse short linear motifs: conformational dynamics insight into the receptor-binding residues

Short linear motifs (SLiMs) have been recognized to perform diverse functions in a variety of regulatory proteins through the involvement in protein-protein interactions, signal transduction, cell cycle regulation, protein secretion, etc. However, detailed molecular mechanisms underlying their functions including roles of definite amino acid residues remain obscure. In our previous studies, we demonstrated that conformational dynamics of amino acid residues in oligopeptides derived from regulatory proteins such as alpha-fetoprotein (AFP), carcino-embryonic antigen (CEA), and pregnancy specific β1-glycoproteins (PSGs) contributes greatly to their biological activities. In the present work, we revealed the 22-member linear modules composed of direct and reverse AFP_14-20-like heptapeptide motifs linked by CxxGY/FxGx consensus motif within epidermal growth factor (EGF), growth factors of EGF family and numerous regulatory proteins containing EGF-like modules. We showed, first, the existence of similarity in amino acid signatures of both direct and reverse motifs in terms of their physicochemical properties. Second, molecular dynamics (MD) simulation study demonstrated that key receptor-binding residues in human EGF in the aligned positions of the direct and reverse motifs may have similar distribution of conformational probability densities and dynamic behavior despite their distinct physicochemical properties. Third, we found that the length of a polypeptide chain (from 7 to 53 residues) has no effect, while disulfide bridging and backbone direction significantly influence the conformational distribution and dynamics of the residues. Our data may contribute to the atomic level structure-function analysis and protein structure decoding; additionally, they may provide a basis for novel protein/peptide engineering and peptide-mimetic drug design.