Charge Interactions Modulate the Encounter Complex Ensemble of Two Differently Charged Disordered Protein Partners of KIX

Molecular dynamics simulations revealed topological frustration in the binding-wrapping process of eIF4G with eIF4E

Interaction between the cap-binding protein eIF4E and the scaffolding protein eIF4G is essential for the cap-dependent translation initiation in eukaryotes. In the Saccharomyces cerevisiae eIF4G/eIF4E complex, the intrinsically disordered eIF4E-binding domain of eIF4G folds into a bracelet-like structure upon binding to eIF4E. Aiming to unveil the molecular mechanism underlying the binding-wrapping process of eIF4G with eIF4E, we performed extensive coarse-grained molecular dynamics simulations and transition path analysis in this work. The major transition pathway revealed from our simulations showed that docking of the eIF4E-binding motif of eIF4G to the folded core of eIF4E initiates the binding process and then the disordered eIF4G wraps around the N-terminal tail of eIF4E. Additionally, we identified a minor transition pathway which indicates the involvement of topological frustration in the binding process. By manipulating the interaction strength of the wrapping contacts and the latching contacts, we further dissected factors affecting the formation of topological frustration and the binding transition kinetics. Our findings provide new clues for experimental studies on the binding mechanism of eIF4G to eIF4E in the future and exemplify the involvement of topological frustration in the binding process of intrinsically disordered proteins.

Myofilament-associated proteins with intrinsic disorder (MAPIDs) and their resolution by computational modeling

The cardiac sarcomere is a cellular structure in the heart that enables muscle cells to contract. Dozens of proteins belong to the cardiac sarcomere, which work in tandem to generate force and adapt to demands on cardiac output. Intriguingly, the majority of these proteins have significant intrinsic disorder that contributes to their functions, yet the biophysics of these intrinsically disordered regions (IDRs) have been characterized in limited detail. In this review, we first enumerate these myofilament-associated proteins with intrinsic disorder (MAPIDs) and recent biophysical studies to characterize their IDRs. We secondly summarize the biophysics governing IDR properties and the state-of-the-art in computational tools toward MAPID identification and characterization of their conformation ensembles. We conclude with an overview of future computational approaches toward broadening the understanding of intrinsic disorder in the cardiac sarcomere.

Regulation of Poly-E Motif Flexibility by pH, Ca2+ and the PPAK Motif

The disordered PEVK region of titin contains two main structural motifs: PPAK and poly-E. The distribution of these motifs in the PEVK region contributes to the elastic properties of this region, but the specific mechanism of how these motifs work together remains unclear. Previous work from our lab has demonstrated that 28-amino acid peptides of the poly-E motif are sensitive to shifts in pH, becoming more flexible as the pH decreases. We extend this work to longer poly-E constructs, including constructs containing PPAK motifs. Our results demonstrate that longer poly-E motifs have a much larger range of pH sensitivity and that the inclusion of the PPAK motif reduces this sensitivity. We also demonstrate that binding calcium can increase the conformational flexibility of the poly-E motif, though the PPAK motif can block this calcium-dependent change. The data presented here suggest a model where PPAK and calcium can alter the stiffness of the poly-E motif by modulating the degree of charge repulsion in the glutamate clusters.

Perspectives on the landscape and flux theory for describing emergent behaviors of the biological systems

We give a review on the landscape theory of the equilibrium biological systems and landscape-flux theory of the nonequilibrium biological systems as the global driving force. The emergences of the behaviors, the associated thermodynamics in terms of the entropy and free energy and dynamics in terms of the rate and paths have been quantitatively demonstrated. The hierarchical organization structures have been discussed. The biological applications ranging from protein folding, biomolecular recognition, specificity, biomolecular evolution and design for equilibrium systems as well as cell cycle, differentiation and development, cancer, neural networks and brain function, and evolution for nonequilibrium systems, cross-scale studies of genome structural dynamics and experimental quantifications/verifications of the landscape and flux are illustrated. Together, this gives an overall global physical and quantitative picture in terms of the landscape and flux for the behaviors, dynamics and functions of biological systems.

Topological frustration leading to backtracking in a coupled folding-binding process

Intrinsically disordered proteins (IDPs) are abundant in all species. Their discovery challenges the traditional "sequence-structure-function" paradigm of protein science because IDPs play important roles in various biological processes without preformed folded structures. Bioinformatic analysis reveals that the intrinsically conformational disorder of IDPs as well as their conformational transition upon binding to their targets is encoded by their amino acid sequences. The rRNase domain of colicin E3 and the immunity protein Im3 are a pair of proteins involved in bacterial survival. While the N-terminal segment and the central segment of E3 make comparable intermolecular contacts with Im3 in the bound state, binding of E3 with Im3 is dominantly triggered by the central segment of E3. In this work, to further investigate the binding mechanism of disordered E3 with Im3, we performed systematic free energy and transition path analysis through coarse-grained molecular dynamics simulations. We observed backtracking of the N-terminal segment of E3 in the binding process, whose occurrence depends on salt concentration. Conformational analysis revealed that initial binding of the N-terminal segment of E3 to Im3 usually leads to misorientation of a central hairpin of E3 on Im3, which generates topological frustration and results in backtracking of the N-terminal segment. Our results not only provide deeper mechanistic insights into the coupled folding-binding process of the E3/Im3 complex, but also suggest that topological frustration could be present in the coupled folding-binding process of IDPs and play an important role in regulating the binding transition pathways.

Physics of biomolecular recognition and conformational dynamics

Biomolecular recognition usually leads to the formation of binding complexes, often accompanied by large-scale conformational changes. This process is fundamental to biological functions at the molecular and cellular levels. Uncovering the physical mechanisms of biomolecular recognition and quantifying the key biomolecular interactions are vital to understand these functions. The recently developed energy landscape theory has been successful in quantifying recognition processes and revealing the underlying mechanisms. Recent studies have shown that in addition to affinity, specificity is also crucial for biomolecular recognition. The proposed physical concept of intrinsic specificity based on the underlying energy landscape theory provides a practical way to quantify the specificity. Optimization of affinity and specificity can be adopted as a principle to guide the evolution and design of molecular recognition. This approach can also be used in practice for drug discovery using multidimensional screening to identify lead compounds. The energy landscape topography of molecular recognition is important for revealing the underlying flexible binding or binding-folding mechanisms. In this review, we first introduce the energy landscape theory for molecular recognition and then address four critical issues related to biomolecular recognition and conformational dynamics: (1) specificity quantification of molecular recognition; (2) evolution and design in molecular recognition; (3) flexible molecular recognition; (4) chromosome structural dynamics. The results described here and the discussions of the insights gained from the energy landscape topography can provide valuable guidance for further computational and experimental investigations of biomolecular recognition and conformational dynamics.

Introducing intrinsic disorder reduces electrostatic steering in protein-protein interactions

Protein-protein interactions underlie many critical biology functions, such as cellular signaling and gene expression, in which electrostatic interactions can play a critical role in mediating the specificity and stability of protein complexes. A substantial portion of proteins are intrinsically disordered, and the influences of structural disorder on binding kinetics and thermodynamics have been widely investigated. However, whether the effect of electrostatic steering depends on structural disorder remains unexplored. In this work, we addressed the consequence of introducing intrinsic disorder in the electrostatic steering of the E3/Im3 complex using molecular dynamics simulation. Our results recapitulated the experimental observations that the responses of stability and kinetics to salt concentration for the ordered E3/Im3 complex were larger than those for the disordered E3/Im3 complex. Mechanistic analysis revealed that the native contact interactions involved in the encounter state and the transition state were essentially identical for both ordered and disordered E3. Therefore, the observed difference in electrostatic steering between ordered E3 and disordered E3 may result from their difference in conformation rather than their difference in binding mechanism. Because charged residues are frequently involved in protein-protein interactions, our results suggest that increasing structural disorder is expected to generally modulate the effect of electrostatic steering.

Assessing the Role of Calmodulin's Linker Flexibility in Target Binding

Calmodulin (CaM) is a highly-expressed Ca2+ binding protein known to bind hundreds of protein targets. Its binding selectivity to many of these targets is partially attributed to the protein’s flexible alpha helical linker that connects its N- and C-domains. It is not well established how its linker mediates CaM’s binding to regulatory targets yet. Insights into this would be invaluable to understanding its regulation of diverse cellular signaling pathways. Therefore, we utilized Martini coarse-grained (CG) molecular dynamics simulations to probe CaM/target assembly for a model system: CaM binding to the calcineurin (CaN) regulatory domain. The simulations were conducted assuming a ‘wild-type’ calmodulin with normal flexibility of its linker, as well as a labile, highly-flexible linker variant to emulate structural changes that could be induced, for instance, by post-translational modifications. For the wild-type model, 98% of the 600 simulations across three ionic strengths adopted a bound complex within 2 μs of simulation time; of these, 1.7% sampled the fully-bound state observed in the experimentally-determined crystallographic structure. By calculating the mean-first-passage-time for these simulations, we estimated the association rate to be ka= 8.7 × 108 M−1 s−1, which is similar to the diffusion-limited, experimentally-determined rate of 2.2 × 108 M−1 s−1. Furthermore, our simulations recapitulated its well-known inverse relationship between the association rate and the solution ionic strength. In contrast, although over 97% of the labile linker simulations formed tightly-bound complexes, only 0.3% achieved the fully-bound configuration. This effect appears to stem from a difference in the ensembles of extended and collapsed states which are controlled by the linker flexibility. Therefore, our simulations suggest that variations in the CaM linker’s propensity for alpha helical secondary structure can modulate the kinetics of target binding.

Thermodynamic and sequential characteristics of phase separation and droplet formation for an intrinsically disordered region/protein ensemble

Liquid–liquid phase separation (LLPS) of some IDPs/IDRs can lead to the formation of the membraneless organelles in vitro and in vivo, which are essential for many biological processes in the cell. Here we select three different IDR segments of chaperon Swc5 and develop a polymeric slab model at the residue-level. By performing the molecular dynamics simulations, LLPS can be observed at low temperatures even without charge interactions and disappear at high temperatures. Both the sequence length and the charge pattern of the Swc5 segments can influence the critical temperature of LLPS. The results suggest that the effects of the electrostatic interactions on the LLPS behaviors can change significantly with the ratios and distributions of the charged residues, especially the sequence charge decoration (SCD) values. In addition, three different forms of swc conformation can be distinguished on the phase diagram, which is different from the conventional behavior of the free IDP/IDR. Both the packed form (the condensed-phase) and the dispersed form (the dilute-phase) of swc chains are found to be coexisted when LLPS occurs. They change to the fully-spread form at high temperatures. These findings will be helpful for the investigation of the IDP/IDR ensemble behaviors as well as the fundamental mechanism of the LLPS process in bio-systems.

The membraneless organelles caused by liquid–liquid phase separation (LLPS) of IDPs/IDRs are involved in a wide range of biological functions such as RNA processing, ribosome biogenesis, and sequestration of mRNA, proteins, and compacted chromatin. In this study, we focus on the histone H2A-H2B binding partner, Swc5 and investigate the effects of the temperature, sequence length and number of charged residues on the LLPS behaviors. Moreover, we proposed three forms of swc conformation in ensemble according to the phase diagram. These characteristics of conformational changes (from LLPS to no LLPS, from condensed-phase to dilute-phase) are observed in the swc chains with different sequence length and charge pattern, which may be the general property for the IDP/IDR ensembles.

Influence of sequence length and charged residues on Swc5 binding with histone H2A-H2B

SWR is a member of chromatin remodeler family and participates the replacement of histone H2A with H2A.Z. One of the SWR subunits, Swc5, has an intrinsically disordered region and binds to H2A-H2B dimer. Though the binding structure of Swc5 and H2A-H2B has been resolved recently, it is still challenging to investigate the binding mechanism as well as the role of the charge interactions between Swc5 and H2A-H2B. Here we developed a coarse-grained structure-based model and performed molecular dynamics simulations to investigate the binding processes of two Swc5 regions with different lengths (swc5-a and swc5-b) to H2A-H2B. The simulation results suggest a different role of electrostatic interactions between swc5-a/swc5-b and H2A-H2B on binding. The electrostatic interactions between swc5-a/swc5-b and H2A-H2B can not only accelerate the initial capture step of binding, but can also trap the swc5-a/swc5-b at the wrong binding site on H2A. Besides, the conserved DEF/Y-2 motif of Swc5 is important for the binding affinity and the recognition with H2A-H2B at the initial step. Both swc5-a and swc5-b undergo a structural shift before reaching the final bound state. This theoretical study provides important details and the underlying physical mechanisms of the binding processes of swc5-a/swc5-b and H2A-H2B.

Binding and folding in transcriptional complexes

Transcription factors are among the classes of proteins with the highest levels of disorder. Investigation of these regulatory proteins is uncovering not just the mechanisms that underlie gene regulation, but relationships that apply to all intrinsically disordered proteins. Recent studies confirm that binding does not necessarily induce folding but that when it does, it tends to follow induced fit mechanisms. Other work emphasises the importance of electrostatics to interactions involving intrinsically disordered proteins, and roles of intrinsic disorder in phase transitions. All these features help direct transcription factors to target sites in the genome to upregulate or downregulate transcription.