Multidomain Protein Research Articles

Domain-wise dissection of thermal stability enhancement in multidomain proteins

Stability is critical for the proper functioning of all proteins. Optimization of protein thermostability is a key step in the development of industrial enzymes and biologics. Herein, we demonstrate that multidomain proteins can be stabilized significantly using domain-based engineering followed by the recombination of the optimized domains. Domain-level analysis of designed protein variants with similar structures but different thermal profiles showed that the independent enhancement of the thermostability of a constituent domain improves the overall stability of the whole multidomain protein. The crystal structure and AlphaFold-predicted model of the designed proteins via domain-recombination provided a molecular explanation for domain-based stepwise stabilization. Our study suggests that domain-based modular engineering can minimize the sequence space for calculations in computational design and experimental errors, thereby offering useful guidance for multidomain protein engineering.

Conserved linker length in double dsRBD proteins from plants restricts interdomain motion

Double stranded RNA binding domains (dsRBDs) are ubiquitous in all kingdoms of life. They can participate both in RNA and protein recognition and are usually present in multiple copies in multidomain proteins. We analyzed the linkers between dsRBDs in different proteins and found that sequences corresponding to plant proteins have a highly conserved linker length. In order to assess the importance of linker length in the conformational freedom of double dsRBD plant proteins, we introduced lanthanide binding tags (LBTs) in different positions of the dsRBD containing protein HYL1 from Arabidopsis thaliana. These constructs were used to obtain conformational restraints from Double electron–electron resonance (DEER) measurements on doubly labeled proteins and from paramagnetic relaxation enhancement (PRE) in single labeled samples. Fitting the experimental datasets to a computational model of the ensemble created by allowing freedom to the linker region we found that the domains tend to explore a particular region of the allowed conformational space. The high conservation in linker length suggests that this restricted conformational sampling is functional, possibly hindering HYL1-dsRBD2 from contacting the substrate dsRNA and allowing it to participate in protein-protein interactions.

Inter-domain distance prediction based on deep learning for domain assembly.

AlphaFold2 achieved a breakthrough in protein structure prediction through the end-to-end deep learning method, which can predict nearly all single-domain proteins at experimental resolution. However, the prediction accuracy of full-chain proteins is generally lower than that of single-domain proteins because of the incorrect interactions between domains. In this work, we develop an inter-domain distance prediction method, named DeepIDDP. In DeepIDDP, we design a neural network with attention mechanisms, where two new inter-domain features are used to enhance the ability to capture the interactions between domains. Furthermore, we propose a data enhancement strategy termed DPMSA, which is employed to deal with the absence of co-evolutionary information on targets. We integrate DeepIDDP into our previously developed domain assembly method SADA, termed SADA-DeepIDDP. Tested on a given multi-domain benchmark dataset, the accuracy of SADA-DeepIDDP inter-domain distance prediction is 11.3% and 21.6% higher than trRosettaX and trRosetta, respectively. The accuracy of the domain assembly model is 2.5% higher than that of SADA. Meanwhile, we reassemble 68 human multi-domain protein models with TM-score ≤ 0.80 from the AlphaFold protein structure database, where the average TM-score is improved by 11.8% after the reassembly by our method. The online server is at http://zhanglab-bioinf.com/DeepIDDP/.

Unfolding Individual Domains of BmrA, a Bacterial ABC Transporter Involved in Multidrug Resistance.

The folding and stability of proteins are often studied via unfolding (and refolding) a protein with urea. Yet, in the case of membrane integral protein domains, which are shielded by a membrane or a membrane mimetic, urea generally does not induce unfolding. However, the unfolding of α-helical membrane proteins may be induced by the addition of sodium dodecyl sulfate (SDS). When protein unfolding is followed via monitoring changes in Trp fluorescence characteristics, the contributions of individual Trp residues often cannot be disentangled, and, consequently, the folding and stability of the individual domains of a multi-domain membrane protein cannot be studied. In this study, the unfolding of the homodimeric bacterial ATP-binding cassette (ABC) transporter Bacillus multidrug resistance ATP (BmrA), which comprises a transmembrane domain and a cytosolic nucleotide-binding domain, was investigated. To study the stability of individual BmrA domains in the context of the full-length protein, the individual domains were silenced by mutating the existent Trps. The SDS-induced unfolding of the corresponding constructs was compared to the (un)folding characteristics of the wild-type (wt) protein and isolated domains. The full-length variants BmrAW413Y and BmrAW104YW164A were able to mirror the changes observed with the isolated domains; thus, these variants allowed for the study of the unfolding and thermodynamic stability of mutated domains in the context of full-length BmrA.

DomBpred: Protein Domain Boundary Prediction Based on Domain-Residue Clustering Using Inter-Residue Distance.

Domain boundary prediction is one of the most important problems in the study of protein structure and function, especially for large proteins. At present, most domain boundary prediction methods have low accuracy and limitations in dealing with multi-domain proteins. In this study, we develop a sequence-based protein domain boundary prediction, named DomBpred. In DomBpred, the input sequence is first classified as either a single-domain protein or a multi-domain protein through a designed effective sequence metric based on a constructed single-domain sequence library. For the multi-domain protein, a domain-residue clustering algorithm inspired by Ising model is proposed to cluster the spatially close residues according inter-residue distance. The unclassified residues and the residues at the edge of the cluster are then tuned by the secondary structure to form potential cut points. Finally, a domain boundary scoring function is proposed to recursively evaluate the potential cut points to generate the domain boundary. DomBpred is tested on a large-scale test set of FUpred comprising 2549 proteins. Experimental results show that DomBpred better performs than the state-of-the-art methods in classifying whether protein sequences are composed by single or multiple domains, and the Matthew's correlation coefficient is 0.882. Moreover, on 849 multi-domain proteins, the domain boundary distance and normalised domain overlap scores of DomBpred are 0.523 and 0.824, respectively, which are 5.0% and 4.2% higher than those of the best comparison method, respectively. Comparison with other methods on the given test set shows that DomBpred outperforms most state-of-the-art sequence-based methods and even achieves better results than the top-level template-based method. The executable program is freely available at https://github.com/iobio-zjut/DomBpred and the online server at http://zhanglab-bioinf.com/DomBpred/.

Site-specific labeling of SBDS to monitor interactions with the 60S ribosomal subunit.

The Shwachman-Diamond syndrome (SDS) is a rare inherited ribosomopathy that is predominantly caused by mutations in the Shwachman-Bodian-Diamond Syndrome gene (SBDS). SBDS is a ribosomal maturation factor that is essential for the release of eukaryotic translation initiation factor 6 (eIF6) from 60S ribosomal subunits during the late stages of 60S maturation. Release of eIF6 is critical to permit inter-subunit interactions between the 60S and 40S subunits and to form translationally competent 80S monosomes. SBDS has three key domains that are highly flexible and adopt varied conformations in solution. To better understand the domain dynamics of SBDS upon binding to 60S and to assess the effects of SDS-disease specific mutations, we aimed to site-specifically label individual domains of SBDS. Here we detail the generation of a fluorescently labeled SBDS to monitor the dynamics of select domains upon binding to 60S. We describe the incorporation of 4-azido-l-phenylalanine (4AZP), a noncanonical amino acid in human SBDS. Site-specific labeling of SBDS using fluorophore and assessment of 60S binding activity are also described. Such labeling approaches to capture the interactions of individual domains of SBDS with 60S are also applicable to study the dynamics of other multi-domain proteins that interact with the ribosomal subunits.

Crystal structure of the CoV-Y domain of SARS-CoV-2 nonstructural protein 3

Replication of the coronavirus genome starts with the formation of viral RNA-containing double-membrane vesicles (DMV) following viral entry into the host cell. The multi-domain nonstructural protein 3 (nsp3) is the largest protein encoded by the known coronavirus genome and serves as a central component of the viral replication and transcription machinery. Previous studies demonstrated that the highly-conserved C-terminal region of nsp3 is essential for subcellular membrane rearrangement, yet the underlying mechanisms remain elusive. Here we report the crystal structure of the CoV-Y domain, the most C-terminal domain of the SARS-CoV-2 nsp3, at 2.4 Å-resolution. CoV-Y adopts a previously uncharacterized V-shaped fold featuring three distinct subdomains. Sequence alignment and structure prediction suggest that this fold is likely shared by the CoV-Y domains from closely related nsp3 homologs. NMR-based fragment screening combined with molecular docking identifies surface cavities in CoV-Y for interaction with potential ligands and other nsps. These studies provide the first structural view on a complete nsp3 CoV-Y domain, and the molecular framework for understanding the architecture, assembly and function of the nsp3 C-terminal domains in coronavirus replication. Our work illuminates nsp3 as a potential target for therapeutic interventions to aid in the on-going battle against the COVID-19 pandemic and diseases caused by other coronaviruses.

Structural organization of RDGB (retinal degeneration B), a multi-domain lipid transfer protein: a molecular modelling and simulation based approach

Lipid transfer proteins (LTPs) that shuttle lipids at membrane contact sites (MCS) play an important role in maintaining cellular homeostasis. One such important LTP is the Retinal Degeneration B (RDGB) protein. RDGB is localized at the MCS formed between the endoplasmic reticulum (ER) and the apical plasma membrane (PM) in Drosophila photoreceptors where it transfers phosphatidylinositol (PI) during G-protein coupled phospholipase C signalling. Previously, the C-terminal domains of RDGB have been shown to be essential for its function and accurate localization. In this study, using in-silico integrative modelling we predict the structure of entire RDGB protein in complex with the ER membrane protein VAP. The structure of RDGB has then been used to decipher the structural features of the protein important for its orientation at the contact site. Using this structure, we identify two lysine residues in the C-terminal helix of the LNS2 domain important for interaction with the PM. Using molecular docking, we also identify an unstructured region USR1, immediately c-terminal to the PITP domain that is important for the interaction of RDGB with VAP. Overall the 10.06 nm length of the predicted RDGB-VAP complex spans the distance between the PM and ER and is consistent with the cytoplasmic gap between the ER and PM measured by transmission electron microscopy in photoreceptors. Overall our model explains the topology of the RDGB-VAP complex at this ER-PM contact site and paves the way for analysis of lipid transfer function in this setting. Communicated by Ramaswamy H. Sarma

Protein kinase C showcases allosteric control: activation of LRRK1.

Allosteric regulation of multi-domain protein kinases provides a common mechanism to acutely control kinase activity. Protein kinase C serves as a paradigm for multi-domain proteins whose activity is exquisitely tuned by interdomain conformational changes that keep the enzyme off in the absence of appropriate stimuli, but unleash activity in response to second messenger binding. Allosteric regulation of protein kinase C signaling has been optimized not just for itself: Alessi and colleagues discover that protein kinase C phosphorylates LRRK1, a kinase with even more domains, at sites on its CORB GTPase domain to allosterically activate LRRK1.

Novel LRR-ROC Motif That Links the N- and C-terminal Domains in LRRK2 Undergoes an Order–Disorder Transition Upon Activation

Mutations in LRRK2, a large multi-domain protein kinase, create risk factors for Parkinson’s Disease (PD). LRRK2 has seven well-folded domains that include three N-terminal scaffold domains (NtDs) and four C-terminal domains (CtDs). In full-length inactive LRRK2 there is an additional well-folded motif, the LRR-ROC Linker, that lies between the NtDs and the CtDs. This motif, which is stabilized by hydrophobic residues in the LRR and ROC/COR-A domains, is anchored to the C-Lobe of the kinase domain. The LRR-ROC Linker becomes disordered when the NtDs are unleashed from the CtDs following activation by Rab29 or by various PD mutations. A key residue within the LRR-ROC Linker, W1295, sterically blocks access of substrate proteins. The W1295A mutant blocks cis-autophosphorylation of S1292 and reduces phosphorylation of heterologous Rab substrates. GaMD simulations show that the LRR-Linker motif, P + 1 loop and the inhibitory helix in the DYGψ motif are very stable. Finally, in full-length inactive LRRK2 ATP is bound to the kinase domain and GDP:Mg to the GTPase/ROC domain. The fundamentally different mechanisms for binding nucleotide (G-Loop vs P-Loop) are captured by these GaMD simulations. In this model, where ATP binds with low affinity (μM range) to N-Lobe capping residues, the known auto-phosphorylation sites are located in the space that is sampled by the flexible phosphates thus providing a potential mechanism for cis-autophosphorylation.

Functional characterization of Vip3Aa from Bacillus thuringiensis reveals the contributions of specific domains to its insecticidal activity

Microbially derived, protein-based biopesticides offer a more sustainable pest management alternative to synthetic pesticides. Vegetative insecticidal proteins (Vip3), multidomain proteins secreted by Bacillus thuringiensis, represent a second-generation insecticidal toxin that has been preliminarily used in transgenic crops. However, the molecular mechanism underlying Vip3's toxicity is poorly understood. Here, we determine the distinct functions and contributions of the domains of the Vip3Aa protein to its toxicity against Spodoptera frugiperda larvae. We demonstrate that Vip3Aa domains II and III (DII-DIII) bind the midgut epithelium, while DI is essential for Vip3Aa's stability and toxicity inside the protease enriched host insect midgut. DI-DIII can be activated by midgut proteases and exhibits cytotoxicity similar to full-length Vip3Aa. In addition, we determine that DV can bind the peritrophic matrix via its glycan-binding activity, which contributes to Vip3Aa insecticidal activity. In summary, this study provides multiple insights into Vip3Aa's mode-of-action which should significantly facilitate the clarification of its insecticidal mechanism and its further rational development.

Transcriptomic analysis of polyketide synthesis in dinoflagellate, Prorocentrum lima

The benthic dinoflagellate Prorocentrum lima is among the most common toxic morphospecies with a cosmopolitan distribution. P. lima can produce polyketide compounds, such as okadaic acid (OA), dinophysistoxin (DTX) and their analogues, which are responsible for diarrhetic shellfish poisoning (DSP). Studying the molecular mechanism of DSP toxin biosynthesis is crucial for understanding the environmental driver influencing toxin biosynthesis as well as for better monitoring of marine ecosystems. Commonly, polyketides are produced by polyketide synthases (PKS). However, no gene has been confirmatively assigned to DSP toxin production. Here, we assembled a transcriptome from 94,730,858 Illumina RNAseq reads using Trinity, resulting in 147,527 unigenes with average sequence length of 1035 nt. Using bioinformatics analysis methods, we found 210 unigenes encoding single-domain PKS with sequence similarity to type I PKSs, as reported in other dinoflagellates. In addition, 15 transcripts encoding multi-domain PKS (forming typical type I PKSs modules) and 5 transcripts encoding hybrid nonribosomal peptide synthetase (NRPS)/PKS were found. Using comparative transcriptome and differential expression analysis, a total of 16 PKS genes were identified to be up-regulated in phosphorus-limited cultures, which was related to the up regulation of toxin expression. In concert with other recent transcriptome analyses, this study contributes to the building consensus that dinoflagellates may utilize a combination of Type I multi-domain and single-domain PKS proteins, in an as yet undefined manner, to synthesize polyketides. Our study provides valuable genomic resource for future research in order to understand the complex mechanism of toxin production in this dinoflagellate.

Disordered proteins mitigate the temperature dependence of site-specific binding free energies

Biophysical characterization of protein-protein interactions involving disordered proteins is challenging. A common simplification is to measure the thermodynamics and kinetics of disordered site binding using peptides containing only the minimum residues necessary. We should not assume, however, that these few residues tell the whole story. Son of sevenless, a multidomain signaling protein from Drosophila melanogaster, is critical to the mitogen-activated protein kinase pathway, passing an external signal to Ras, which leads to cellular responses. The disordered 55kDaC-terminal domain of Son of sevenless is an autoinhibitor that blocks guanidine exchange factor activity. Activation requires another protein, Downstream of receptor kinase (Drk), which contains two Src homology 3 domains. Here, we utilized NMR spectroscopy and isothermal titration calorimetry to quantify the thermodynamics and kinetics of the N-terminal Src homology 3 domain binding to the strongest sites incorporated into the flanking disordered sequences. Comparing these results to those for isolated peptides provides information about how the larger domain affects binding. The affinities of sites on the disordered domain are like those of the peptides at low temperatures but less sensitive to temperature. Our results, combined with observations showing that intrinsically disordered proteins become more compact with increasing temperature, suggest a mechanism for this effect.

Understanding the molecular basis of folding cooperativity through a comparative analysis of a multidomain protein and its isolated domains

Although cooperativity is a well-established and general property of folding, our current understanding of this feature in multidomain folding is still relatively limited. In fact, there are contrasting results indicating that the constituent domains of a multidomain protein may either fold independently on each other or exhibit interdependent supradomain phenomena. To address this issue, here we present the comparative analysis of the folding of a tandem repeat protein, comprising two contiguous PDZ domains, in comparison to that of its isolated constituent domains. By analyzing in detail the equilibrium and kinetics of folding at different experimental conditions, we demonstrate that despite each of the PDZ domains in isolation being capable of independent folding, at variance with previously characterized PDZ tandem repeats, the full-length construct folds and unfolds as a single cooperative unit. By exploiting quantitatively, the comparison of the folding of the tandem repeat to those observed for its constituent domains, as well as by characterizing a truncated variant lacking a short autoinhibitory segment, we successfully rationalize the molecular basis of the observed cooperativity and attempt to infer some general conclusions for multidomain systems.

A unified approach to protein domain parsing with inter-residue distance matrix.

It is fundamental to cut multi-domain proteins into individual domains, for precise domain-based structural and functional studies. In the past, sequence-based and structure-based domain parsing was carried out independently with different methodologies. The recent progress in deep learning-based protein structure prediction provides the opportunity to unify sequence-based and structure-based domain parsing. Based on the inter-residue distance matrix, which can be either derived from the input structure or predicted by trRosettaX, we can decode the domain boundaries under a unified framework. We name the proposed method UniDoc. The principle of UniDoc is based on the well-accepted physical concept of maximizing intra-domain interaction while minimizing inter-domain interaction. Comprehensive tests on five benchmark datasets indicate that UniDoc outperforms other state-of-the-art methods in terms of both accuracy and speed, for both sequence-based and structure-based domain parsing. The major contribution of UniDoc is providing a unified framework for structure-based and sequence-based domain parsing. We hope that UniDoc would be a convenient tool for protein domain analysis. https://yanglab.nankai.edu.cn/UniDoc/. Supplementary data are available at Bioinformatics online.

From elucidating mechanisms of allosteric signaling to designing allosteric drugs and biologics.

Allostery is a pervasive phenomenon underlies a multitude of cellular processes. We have previously developed a Structure-Based Statistical Mechanical Model of Allostery (SBSMMA) to quantitatively describe the causality and energetics of allosteric signaling due to perturbations such as ligand binding, mutations, disorder-order transition and post-translational modifications. The computational model was implemented in the AlloSigMA server and AlloMAPS database, which allow one to explore allosteric effects of perturbations in form of the exhaustive Allosteric Signaling/Probing Map (ASM/APM), to identify potential allosteric sites, and to facilitate effector design. Here, we demonstrate some recent work on applying our broadly applicable framework to different aspects of allostery. First, analysis of communication between different sites/regions of the Spike protein from SARS-CoV-2 revealed mutations harbored by variants of concern, as well as glycosylation, that originate strong signaling to the distant receptor-binding domain despite the separation. Using reverse perturbation of allosteric communication and ASM/APM, we identified potential druggable sites that can induce tunable allosteric responses in different functional regions of Spike. Second, while the strengths of allosteric medicines are widely recognized, rational design of effectors remains a challenging problem. We demonstrated with GPCRs and kinases that effector design should begin with the search/design for latent/non-natural binding sites, followed by fragment-based design of ligands, and consequently mutual adjustments of site-ligand pairs to obtain and tune the desired binding affinity and allosteric signals. Third, rapid progress in harnessing allostery calls for an understanding of how communication is established by different structural elements and how it can be modified. We obtained a picture of conserved allosteric communication characteristic in major folds, alterations of structure-driven signaling via sequence-determined divergence to distinct functions, and its diversification in multi-domain proteins and oligomers, which will facilitate engineering/design of allosteric communication.

Molecular architecture of LRRK2 in association with cellular organelles

Mutations in leucine-rich repeat kinase 2 (LRRK2) lead to increased risk of Parkinson's disease. LRRK2 is a large multidomain protein that regulates vesicular trafficking, cytoskeleton dynamics, and lysosomal function. However, the exact mechanism of how LRRK2 carries out these various functionalities remains elusive. Elucidating the localization of LRRK2 and mapping its association with microtubules and cellular vesicles will shed light on the mechanism of action of the protein in normal and disease states. Combining cryo-correlative light and electron microscopy, cryo-FIB milling, cryo-electron tomography, subtomogram analysis and integrative modeling we aim to characterize various functional states of LRRK2 on cellular organelles. We first applied this workflow to determine the molecular architecture of LRRK2 as it associates with microtubules. In this study we aim to characterize association of LRRK2 with microtubules and cellular vesicles in human brain microglial cells. We aim to characterize WT as well as PD causing mutations and study the effects of kinase inhibitors on its association with cellular organelles. This study will provide important insight into how LRRK2 operates in its native and disease state.

Understanding the role of the LRRK2 tail in the crosstalk between kinase and the N-terminal domains.

LRRK2 is a multi-domain protein encompassing biology's two most critical molecular switches; a kinase and a GTPase. Mutations in LRRK2 are one of the key players in the pathogenesis of Parkinson's disease (PD). The availability of multiple structures (full-length and truncated) has opened doors to exploring the physiological role of LRRK2. A helix extending from the WD40 domain and stably docking onto the kinase domain is a common feature in all available structures. This C-terminal (Ct) helix is a hub of phosphorylation and organelle-localization motifs and can serve as a multi-functional protein: protein interaction module. To examine its intra-domain interactions, we have recombinantly expressed a stable Ct motif (residues 2480-2527) and are using peptide arrays to identify specific binding sites. We have identified a potential interaction site between the Ct helix and a loop in the CORB region with a combination of Gaussian MD simulations and peptide arrays. This Ct-Motif contains multiple phosphorylation sites including two auto-phosphorylation sites (Thr2083 and Thr2524), a 14-3-3 binding site (Thr2524) and a putative PKA phosphorylation site. These two sites together form a part of a dynamic “CAP” that regulates the N-lobe of the kinase domain. We also explore how the cellular localization of LRRK2RCKW and its docking onto microtubules is altered by mutations in the Ct helix, if at all. We hypothesize that in inactive, full-length LRRK2, the Ct-helix will mediate interactions with the N-terminal armadillo, ankyrin, and LRR domains (NTDs) and that binding of Rab substrates, PD mutations, or kinase inhibitors will unleash the NTDs.

Structural features of EWSR1 affecting Ewing's sarcoma. Computational domain assembly and conformational dynamics studies.

Proteins are key players in the prognosis and therapeutics of various sarcomas by influencing downstream signaling cascades. This work provides insight on the structural and conformational interactive behavior of EWS RNA-Binding Protein 1 (EWSR1) with RNA associated with Ewing's sarcoma. The comparative modeling approach was employed to predict the structural models of EWSR1 domains separately, and the Domain Enhanced MOdeling (DEMO) server was used for multidomain protein structure assembly to predict the structure composed of all three domains. Furthermore, RNA motifs binding to EWSR1 were predicted and the 3D model of RNA bound to EWSR1 was built by using MC-Fold. We used HDOCK docking server to check the conformational interactions between EWSR1 and RNA. Moreover, Molecular Dynamics simulations were performed to check the stability of EWSR1-RNA complex by computing the root mean square deviation and fluctuation (RMSD/F), radius of gyration (Rg) and solvent accessible surface area (SASA). Our results explored the structural insights on EWSR1domains and their interactive binding with RNA. Based on computational assessments, it has been concluded that certain structural and dynamical features of EWSR1-RNA complex could be used as a target in future development of drugs for Ewing's sarcoma.

Hydrophobic residues in disordered and helical regions mediate the oligomerization and phase separation of TDP-43 C-terminal domain.

TAR DNA binding protein 43 (TDP-43) is a nuclear, multidomain protein implicated in RNA metabolism. The cytoplasmic accumulation of TDP-43 can result in the formation of inclusion bodies which are a hallmark of neurodegenerative diseases such as amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). The prion-like, disordered C-terminal domain (CTD) of TDP-43 (spanning residues 267-414) is aggregation-prone and harbors the majority (∼90%) of all ALS-related mutations. Several studies report that CTD can undergo liquid-liquid phase separation (LLPS). The conserved, hydrophobic region (CR, spanning residues 319-341) in CTD populates an α-helical structure which is essential for phase separation. The requirement of CR for phase separation is attributed to its tendency to undergo helix-helix self-association and promote oligomerization. In addition to CR, aromatic residues in the flanking disordered regions (IDR1/2) also play a critical role in CTD self-assembly. However, the relative contributions of critical residue types in CR and flanking regions on the overall phase separation propensity of CTD are unclear. Here, we utilized multiscale simulations coupled with biophysical experiments to determine the role of hydrophobic residues in the CR and IDR1/2 regions on CTD phase separation. Single chain and condensed phase simulations collectively identified the contribution of aromatic residues while also suggesting significant contributions to LLPS by other residue types, which were previously not recognized as essential for phase separation within and outside the CR region. The computational observations were successfully corroborated by in vitro microscopy and saturation concentration measurements. Finally, NMR experiments determined that CTD oligomerization was critically dependent on these residues within CR and IDR1/2 regions. Overall, our work uncovers the collective importance of hydrophobic residues across both helical and disordered regions in modulating the oligomerization and phase separation of CTD.