# Role of integrative structural biology in understanding transcriptional initiation

Michael Trnka, Riccardo Pellarin, Philip Robinson

► To cite this version:

Michael Trnka, Riccardo Pellarin, Philip Robinson. Role of integrative structural biology in understanding transcriptional initiation. Methods, 2019, 159-160, pp.4-22. 10.1016/j.ymeth.2019.03.009. pasteur-02413976

HAL Id: pasteur-02413976

https://pasteur.hal.science/pasteur-02413976

Submitted on 16 Dec 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# Role of integrative structural biology in understanding transcriptional initiation

<figure><img src="image_2.png" /></figure>

Michael J. Trnka<sup>a</sup>, Riccardo Pellarin<sup>b</sup>, Philip J. Robinson<sup>c,*</sup>

## ARTICLE INFO

### Keywords:

Integrative structural determination
RNA polymerase II
Transcriptional initiation machinery
Cross-linking mass spectrometry
Multi-subunit complex architecture
Computational molecular modelling
Cryo-EM

## ABSTRACT

Integrative structural biology combines data from multiple experimental techniques to generate complete structural models for the biological system of interest. Most commonly cross-linking data sets are employed alongside electron microscopy maps, crystallographic structures, and other data by computational methods that integrate all known information and produce structural models at a level of resolution that is appropriate to the input data. The precision of these modelled solutions is limited by the sparseness of cross-links observed, the length of the cross-linking reagent, the ambiguity arisen from the presence of multiple copies of the same protein, and structural and compositional heterogeneity. In recent years integrative structural biology approaches have been successfully applied to a range of RNA polymerase II complexes. Here we will provide a general background to integrative structural biology, a description of how it should be practically implemented and how it has furthered our understanding of the biology of large transcriptional assemblies. Finally, in the context of recent breakthroughs in microscope and direct electron detector technology, where increasingly EM is capable of resolving structural features directly without the aid of other structural techniques, we will discuss the future role of integrative structural techniques.

## 1. Introduction

The goal of structural biology is to derive detailed functional and mechanistic information on a biomolecule from the arrangement of its constituent atoms. Since their introduction in the 20th century the two techniques at the forefront of structural biology, X-ray crystallography and electron microscopy (EM), have undergone seismic changes leading to an explosion in the numbers of scientists practicing structural biology and the pursuit of biological targets of ever-increasing complexity. X-ray crystallography has trended from focusing on small soluble proteins and individual protein domains [1–3] to large multi-subunit complexes and challenging membrane-associated targets [4–7]. In contrast electron microscopy has seen an evolution from massive, highly symmetrical macromolecules [8,9], to ever-smaller and asymmetric targets [10]. This evolution reflects an important underlying principle of biology; namely, that biological processes often occur in the context of large and asymmetrical protein, lipid-associated or nucleoprotein assemblies. The function of these assemblies is often dependent on a large set of more transient biomolecule interactions and on significant levels

of domain motion, which act as sources of compositional and conformational heterogeneity, respectively. Such biological heterogeneity has provided significant technological challenges and has frequently exposed the limitations of individual structural biology techniques. For example, the formation of a highly populated and ordered macromolecular crystal lattice, which forms the basis for high-resolution X-ray diffraction data, is dependent on the deposition of a population of macromolecules that are identical in respect to both their composition and conformation. Sources of heterogeneity lead to lattice ‘poisoning’ that attenuate the intensity of useful Bragg diffraction. Likewise, in the context of electron micrographs, which capture images of individual macromolecules with low signal-to-noise, sources of macromolecular heterogeneity lead to increased errors in the assignment of projection angle [11,12] and put extra demands on the computational routines designed to separate particles into distinct structural classes [13].

Due to this limitation of individual structural techniques, studies of large macromolecular assemblies have often been limited to distinct subassemblies, which suffer less from problems due to heterogeneity and therefore can be solved at higher resolution. Integrative structural

<footnote><sup>a</sup>Department of Pharmaceutical Chemistry, University of California San Francisco, San Francisco, CA, USA</footnote>

<footnote><sup>b</sup>Institut Pasteur, Structural Bioinformatics Unit, Department of Structural Biology and Chemistry, CNRS UMR 3528, C3BI USR 3756 CNRS & IP, Paris, France</footnote>

<footnote><sup>c</sup>Department of Biological Sciences, Birkbeck University of London, Institute of Structural and Molecular Biology, London, United Kingdom</footnote>

<footnote>* Corresponding author at: Department of Biological Sciences, Birkbeck University of London, Institute of Structural and Molecular biology, Malet Street, London WC1E 7HX, United Kingdom.</footnote>

<footnote>E-mail address: p.robinson@mail.cryst.bbk.ac.uk (P.J. Robinson).</footnote>

<footnote>https://doi.org/10.1016/j.ymeth.2019.03.009</footnote>

<footnote>Received 15 January 2019; Received in revised form 14 March 2019; Accepted 15 March 2019</footnote>

<footnote>Available online 16 March 2019</footnote>

<footnote>1046-2023/© 2019 The Authors. Published by Elsevier Inc. This is an open access article under the CC BY license</footnote>

<footnote>(http://creativecommons.org/licenses/BY/4.0/).</footnote>
biology approaches were developed to overcome the shortfalls of in-
dividual structural techniques and provide a framework for combining
the data from multiple structural approaches to form a more complete
picture of dynamic biological assemblies. A branch of this integrative
approach, Integrative Structural Determination (ISD), attempts to use
as much of the relevant biochemical and biophysical data about a
macromolecular complex as possible to generate three-dimensional
structures (models). ISD merges data sets that individually cannot lead
to the unambiguous (e.g. atomistic) structural determination of the
whole macromolecular complex. ISD exploits the mutual synergy and
consistency of the datasets in such a way that the resulting model
precision is higher (and therefore more informative) than the precision
of the models generated by each individual dataset. In the ISD ap-
proach, datasets are encoded into a scoring function which is able to
rank models according to how compatible they are to the input data.
The configurational space of the represented complex is then ex-
tensively and thoroughly sampled. Best-scoring models (the ensemble
of solutions) are validated and analyzed to assess their quality and build
testable hypotheses [14,15].

Importantly, two main problems can limit an effective ISD ap-
proach: the lack of data and the incomplete sampling of the models. The
former can be due to the instability of the experimental sample or
technological limitations, while the latter is related to the huge number
of degrees of freedom of large protein complexes. To overcome these
issues, the ISD approach can exploit two important principles that are
intrinsic to the hierarchical architecture of assemblies: the *transferability* of datasets and the *modularity* of the structure (Fig. 4). The transfer-
ability assumes that the bulk of the data is robust, irrespective of the
functional state of the complex. Using this principle, one can collect
data on a functional state where the protein complex is more stable, and
use the data to model the functional state of interest. The modularity
principle assumes that the complex is formed of several stable and in-
dependent sub-complexes, whose architecture varies only modestly
upon changes in the functional state. The modularity allows in-
vestigators to assign the data to distinct components of the system in a
divide-and-conquer fashion, allowing a more efficient model sampling.

Over the last decade the ISD approach has been remarkably suc-
cessful in modeling the architectures of large, fundamentally important
complexes whose structures seemed intractable to single techniques
such as those involved in transcription (pre-initiation complex,
Mediator, TFIIH, etc.), translation (ribosome, eIF3) and transit across
the nucleus (nuclear pore complex) [16–24]. Many of these applications
assumed the modularity and transferability principles.

In this review we will provide a general description of the in-
tegrative modeling pipeline and then review how it has been im-
plemented to provide novel details of a number of complexes involved
in RNA polymerase II (RNAPII) transcription. Finally, we will provide a
discussion of the future of integrative modeling in the context of recent
transformational developments in cryo-EM.

## 2. RNAPII transcription initiation machinery

The initiation of transcription from eukaryotic genes depends on the
assembly at gene promoters of a massive pre-initiation complex (PIC)
involving RNAPII and the general transcription factors (GTFs) TFIIA, -B,
-D, -E, -F, -H and -S. This 31-protein assembly includes proteins that
recognise and bind conserved regulatory promoter DNA elements up-
stream of the transcriptional start site, distort the DNA to direct its path
along the active site cleft of RNAPII and apply DNA helical torsion in an
ATP-dependent manner to melt the DNA duplex and promote single-
stranded template DNA binding within the RNAPII active site. The
Megadalton Mediator complex, traditionally classified as a transcrip-
tional coactivator complex, also plays an essential role in the initiation
pathway [25,26]. Mediator binds transcriptional activator proteins and
supports the stimulation of transcription at inducible gene promoters in
response to activator UAS binding [27,28]. The augmented

transcriptional response to activators depends on the integrity of the
unstructured C-terminal domain (CTD) of RNAPII [29–31], which is
highly conserved across eukaryotes and composed of tandem hepta-
peptide repeats. Mediator binds to the CTD and in doing so acts as a
bridge between cellular regulatory signals and the basal transcriptional
machinery. Mediator-RNAPII interactions are responsible for increasing
RNAPII levels and transcriptional output at gene promoters in vivo
[32,33]. Mediator binds to unmodified CTD repeats [34] and RNAPII
molecules with unmodified CTD sequence are exclusively recruited to
the PIC [35]. Upon full PIC assembly the CTD is brought together with
the CTD kinase TFIIK (a submodule of TFIIH), and the resulting CTD
phosphorylation correlates with loss of the Mediator-RNAPII interac-
tion and the transition to the elongation phase of the transcription cycle
(RNAPII promoter escape) [36,37]. A striking feature of the gene acti-
vation pathway is the modular organisation of the transcriptional ap-
paratus involved. Most of the polypeptides constituting the PIC are
members of multi-subunit sub-assemblies that must be brought together
as ‘ingredients’ for efficient transcriptional initiation. This modularity
provides ample opportunities to regulate this early stage of transcrip-
tion, such as in the example described above, where post-translational
modification of the CTD is used to regulate the stability and residency
of factors at the promoter. A further layer of regulation may exist
through the compartmentalization of the transcription apparatus into
so-called membrane-less compartments within the nucleus. Recent data
have shown that the Mediator complex plays a critical role in bridging
enhancer and promoter elements in the genome [38,39]. An emerging
view in the field is that in Metazoans the colocalization of super-en-
hancer-bound transcription factors, coactivators such as Mediator, and
the promoter-associated general transcription machinery may in part be
driven through a process termed biomolecular condensation, where
multivalent interactions involving intrinsically disordered domains
contribute to a phase-separated state producing regions with a high
local concentration of essential transcription components. Support for
this idea comes from recent observations of nuclear puncta in which
Mediator and RNAPII colocalise with super-enhancer elements in
structures that have properties of phase-separated condensates such as
sensitivity to a non-specific aliphatic alcohol, 1,6-hexanediol [40–42].
A comprehensive biochemical and structural dissection of the interac-
tion between the yeast activator Gcn4 and its Mediator target, Med15,
has described a large heterogenous “fuzzy” complex comprised of
multiple low affinity hydrophobic interactions that interact additively
in a non-specific free-for-all manner to increase the overall binding
affinity [43,44]. The fact that the activation domains of transcription
factors such as Gcn4 are archetypical intrinsically disordered regions
brings into question the contribution of the “fuzzy” hydrophobic in-
terface to the formation of membrane-less transcriptional condensates
in vivo [45].

## 3. Novel biological findings from integrative structural biology-based studies of transcription complexes

### 3.1. Introduction to integrative structural determination

The ISD workflow follows a series of well-defined steps (Fig. 1)
[15,46]. First, all available primary data for the system of interest are
gathered (see Section 4.4 for details). Datasets generally used by ISD
include low-resolution three-dimensional EM reconstructions, small-
angle X-ray scattering (SAXS), atomic-level data for subunits or do-
mains, cross-linking mass spectrometry (CLMS), affinity pulldowns,
labeling and co-localization experiments. The second stage involves
generating a multi-scale representation of all components of the multi-
subunit assembly and translating data into spatial restraints (see
Section 4.5 for details). For certain domains in which atomistic struc-
tural information in the form of either homology models or crystal
structures are available, these regions will be represented at the atomic
or residue level as fixed rigid-body structures. Other domains will have
<figure><img src="image_3.png" /><figcaption>Fig. 1. Representation of the multiple stages of integrative structural determination.</figcaption></figure>

no atomic resolution structural information and will be represented
with chains of coarse-grained beads, with a suitable excluded volume,
and connected by the continuity of the polypeptide chain. The third
stage in the modeling workflow involves performing exhaustive com-
putational sampling of the configurational space in order to identify a
population of model solutions that best satisfy the full set of input re-
straints (see Section 4.6 for details). In practice this is achieved by
computing parallel modeling trajectories or replicas (Replica Exchange)
[47], each starting from a randomised starting configuration. Upon
convergence of the computational sampling the ensemble of top-scoring
solutions is classified by structural similarity to detect discrete struc-
tural subpopulations (see Sections 4.7 and 4.8 for details). At this stage
subunit densities are calculated for each subpopulation cluster and in-
formation relating to how precisely each subunit in the model has been
localised can be generated and graphically represented. Also, in-
formation relating to protein domain proximities in the complex can be
extracted and graphically represented to reveal interesting biology and

form the basis for future biochemical characterisation. To validate the
resulting models a series of tests are performed that are designed to
check restraint data robustness and the completeness of sampling (see
Section 4.9 for details).

Over the last few years, we exploited the modular properties of the
transcription machinery to model the architecture of increasingly
complex transcription factor assemblies using ISD. We proceeded
hierarchically, building on structural studies of isolated subcomplexes,
and eventually characterized the large Mediator-bound PIC (Med-PIC).

3.2. Mediator Head module structure

Biochemical and structural studies have established that the
Mediator complex is composed of three modules, the Head, Middle and
Tail modules [48,49], that have a roughly equal share of the 21 sub-
units that comprise the core yeast Mediator complex. The subunit
composition of each module was defined by a combination of
biochemical and genetic analyses [49–54]. Genetic screens for suppressors of a truncated RNAPII CTD identified a number of so-called SRB genes, five of which encoded core Mediator subunits [27,50,55]. Four of the five core Mediator SRB gene products were found to reside in the Head module, a 220 kDa assembly of seven proteins, most of which were found to be essential for yeast viability. Taken together with biochemical experiments showing that the full Mediator displays a robust CTD binding response [27,50,51], these genetic observations suggested a functional interaction between the CTD and the Head module with a direct interaction likely. Under non-permissive conditions, conditional mutations in the largest Head module subunit, Med17 (Srb4), lead to a dramatic global reduction of mRNA levels brought about through a destabilisation of the Mediator complex [25,26,56]. Therefore, the Mediator complex acts through the Head module to perform an essential role in the pathway to transcriptional activation at most eukaryotic gene promoters. A number of studies set out to elucidate the structural basis of Head module activities in transcription. The structure of a Med18-Med20 complex, encoded by two of the SRB genes, was the first atomic level information available, but provided little framework for understanding Head module function [57]. Structural information was extended to include the full seven-subunit Head module isolated from recombinant preparations of the yeast complex using X-ray crystallography [58]. From the resulting 4.2 Å density maps an architectural polyalanine model, composed mostly of alpha helices, revealed roughly 60% of the primary sequence (PDB: 3RJ1). The structure showed that the Head module can be divided into three independent structured domains, two termed the “Mobile Jaw” and “Fixed Jaw” based upon EM observations [59], which are separated from a third highly helical “Neck” domain by a mainly unstructured “Joint” region (Fig. 2B). In the absence of full chain connectivity and fully resolved side chain density, the very challenging job of assigning and threading subunit sequence into the secondary structure elements of this architectural model was guided by structure predictions and selenomethionine difference density peaks [58].

Shortly after publication of the recombinant yeast Head module we solved the structure of the native yeast Mediator Head module both alone and bound to the RNAPII CTD [34]. A strong set of phases from multiple isomorphous replacement with anomalous scattering (MIRAS) experiments generated maps that, whilst at a similar resolution to that achieved for the recombinant Head module, revealed features not seen previously such as domain connectivity in the central “Joint” region and large sections of β-sheet within both the “Neck” and the “Fixed Jaw” domains. Within these maps elongated CTD density comprising around four heptad repeats was found to follow a path across a highly conserved surface of the “Neck” domain. Modeling the new Head module density would provide a model that accounted for 80% of the primary sequence. However, the following challenges presented themselves during early model-building efforts: (1) sequence assignment for helices in subunits Med11 (“Neck”) and Med17 (“Fixed Jaw”) conflicted with the existing Head module structure (PDB: 3RJ1), (2) sequence ambiguities for inter-domain connectivity within the “Joint” region remained. To address these challenges, we employed an integrative approach where CLMS derived distance restraints were used to help revise sequence assignments and accurately model domain connectivity within the “Joint” region.

Revising the Med11 sequence assignment involved correctly determining the directionality of Med11N-term helix 1 (α1) and helix 2 (α2), which are core elements in a 10-helix bundle within the “Neck” region. Cross-links within the “Neck” were consistent only with a Med11 assignment in which the directionality of Med11 α1 and α2 were reversed from that proposed in 3RJ1 (Fig. 2B). Such a reversal also positioned the C-terminal end of Med11 α2 within continuous map density that connected to the Med11 C-term, and was consistent with a high-resolution X-ray crystal structure of an isolated Med11-Med22 four-helix bundle [60].

Revising C-terminal Med17 sequence involved, amongst other

things, determining the correct sequence connectivity between helices that were previously modelled as unconnected elements and whose sequence had been assigned based mostly on a sparse collection of selenomethionine difference peaks. CLMS analysis supported the re-assignment of Med17 C-terminal 597-611 and 670-685 helices and surrounding sequence (Fig. 2B). For example, based upon the 3RJ1 Med17 sequence assignments, cross-links between Med17 K589, K601 and K608 with five surrounding lysines gave an average Cβ cross-link distance of 40.9 Å, significantly greater than the typical violation cut-off of 30 Å. After reassignment the average Cβ cross-link distance for the same collection of cross-links was 20.6 Å.

A cluster of SRB mutations map to residues within the “Joint” region of the Head module [50,55]. A first step to understanding the structural basis of these mutations was to model this region of the structure, since it remained completely unmodeled in 3RJ1. Despite the presence of clear and continuous protein backbone densities, the combination of low secondary structure content, elevated B-factors and insufficient map resolution to position side-chain densities made “Joint” region modeling particularly challenging. Ambiguities were resolved by incorporating constraints from CLMS analysis into the molecular modeling workflow. A relatively large number of cross-links mapped to the “Joint” region, with a particular concentration in a Med17 hotspot comprised of six independent Med17 cross-links. These strong cross-linking constraints facilitated modeling of a Med17 β-ribbon that encompassed most of the cross-link hotspot (Fig. 2B). Importantly, modeling of this region of Med17 dramatically simplified the process of modeling into the remaining “Joint” density. To our knowledge, this remains a rare example of the use of cross-link constraints to interpret X-ray crystallography maps and highlights the power of integrative structural approaches to overcome the weaknesses inherent in individual structural approaches.

### 3.3. Unique architectural roles of Med17 and Med14 in the Mediator complex

The last few years have seen huge gains in our understanding of the organisation and molecular structure of the Mediator complex, both as an independent complex and also when engaged with other components of the transcriptional apparatus. Aside from the Head module and a few isolated domains from the Middle [61,62] and Tail modules [43,63], until very recently only a basic view of the organisation and subunit architecture of the 21-subunit core Mediator complex was available. Studies probing subunit interaction networks [49,64,65] as well as subunit co-expression [53,54,66,67] and subcomplex isolation studies [68,69] combined to determine the subunit composition of the Mediator modules. This data was augmented by EM studies using different labelling strategies to coarsely map the termini of various Mediator subunits onto two-dimensional projection views of the complex viewed in negative stain [70,71]. However, until very recently no three-dimensional architectural information was available for the Middle and Tail modules, together accounting for three quarters of the mass of the Mediator complex.

The transition from coarse two-dimensional subunit mapping to the first full three-dimensional Mediator model resulted from the application of an ISD strategy to the structure of the RNAPII Holoenzyme, the complex of Mediator and RNAPII. The Holoenzyme complex is not only a functionally relevant form of Mediator, but is also more soluble and stable than free Mediator alone under the conditions used for cross-linking experiments. Using BS3 cross-linking, about 400 cross-links were identified for the Holoenzyme. Despite the availability of a low-resolution Mediator EM map [70], it was impossible to generate a model of the Holoenzyme using the cross-links without a complete Holoenzyme EM map, or a model of the Mediator using the EM map without the cross-links. As a consequence, we envisioned an ISD strategy for the isolated Mediator complex based on the transferability and the modularity principles described above. First, Holoenzyme
<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Novel RNAPII transcription complex biology determined using integrative structural methods. (A) Architecture of the yeast Mediator complex. Mediator subunit localization densities (left side) determined for the free Mediator using ISD methods, docked into Med-PIC cryo-EM map (grey transparent density). Atomic model of the yeast Head and Middle lacking Med1 (right side) from recent cryo-EM studies [77] oriented for direct comparison to the docked Mediator ISD model. Coloured coordinate axes show the orientation of the models relative to the Med-PIC cryo-EM map (panel C). (B) Integrative structural solution of the Mediator Head module. A cross-link dataset (left side; yellow colour) was used to guide model building into X-ray crystallographic map density at 4.2 Å resolution. Alternative view of the Head module (right side) showing regions in which revisions to sequence assignment or novel model build were supported by CLMS distance restraints (red colour). (C) cryo-EM structure of the complete Mediator-bound transcriptional pre-initiation complex [18]. Density for RNAPII, general transcription factors, promoter DNA and individual Mediator modules is coloured to facilitate map interpretation. Coloured coordinate axes provide a reference to orient structural components in other figure panels. (D) Integrative structural determination of TFIIK architecture. TFIIK subunits (Kin28, Tfb3 and Ccl1) were localized to novel Med-PIC cryo-EM map density (grey mesh) by integrating CLMS and EM restraints (left side). Modeling was supported by docking of a homology model for the TFIIH ATPase RAD3 into the cryo-EM density. Cross-links between TFIIK subunits and RAD3 together with EM restraints supported the integrative structural determination of TFIIK subunit localisation maps (solid densities). The position of homology models that contributed to the multi-scale representations of TFIIK subunits are shown positioned within their corresponding localization densities. Schematic representation of the inter-subunit cross-links that provided distance restraints for TFIIK integrative modeling (right side). (E) Dynamics and GCN4 cross-linking of the Mediator Tail module. Cryo-EM map from a particle sub-population in which the Mediator Tail module was in the ‘Tail-Up’ conformation (left side). In the ‘Tail-Up’ conformation the Tail moves closer to RNAPII and the upstream promoter DNA with bound activators. Schematic representation of the cross-links between the activator GCN4 and its main target within Mediator, Med15. The transcriptional activation domain and DNA-binding bZIP domain of GCN4 are indicated. The four Med15 domains known to be targets of GCN4 (KIX domain and activator binding domains 1–3) are indicated. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)</figcaption></figure>

cross-links were transferred to the Mediator-alone structural determination, guided by the low-resolution map. Second, we enforced the modularity principle, by segmenting the Mediator EM density map into regions corresponding to the Head, Middle, and Tail modules of Mediator (Fig. 4).

The Mediator model generated by this ISD approach positioned each subunit within the low-resolution Mediator envelope and provided an unambiguous picture of the organisation of the Mediator Middle and Tail modules (Fig. 2A) [17]. As well as global architectural details, the approach provided a structural framework for interpreting the role played by individual Mediator subunits in the complex. Central to understanding the organisation of the Mediator complex was shedding light on the role of the two large and essential subunits Med14 and Med17. Originally identified in a genetic screen in yeast for genes required for glucose repression [72], Med14 had been classified as a member of both the Middle and Tail modules. Furthermore, truncation of the C-terminal portion of Med14 was shown to destabilise the interaction between these same Mediator modules [73]. Reconstitution of a stable and functionally active human core Mediator complex (Head and Middle subunits) was shown to depend on the presence of Med14, suggesting an important scaffolding role within the complex [74]. By

providing the three-dimensional subunit architecture of the Mediator complex, the ISD approach was the first to shed light on the structural basis of the unique scaffolding role played by Med14 (Fig. 2A). Subunit density for Med14 was mapped over a vast portion of the Mediator complex, spanning a total distance of 220 Å. Cross-linking patterns show that the N-terminus of Med14 is positioned to interact with Med10 and 19 at one extreme of the Mediator Middle module whilst central regions of Med14 encompass the remainder of the Middle. The cross-linking pattern of the C-terminal portion of Med14 showed a distinctive pattern characterised by colocalization with the Tail module proteins Med15 and 16. This observation helped to localise the junction between the Mediator Middle and Tail modules and hence explained the loss of the Tail module in Mediator preparations harbouring a truncation of the C-terminal portion of Med14 [73]. The proximity of Med14 to subunits from all three Mediator modules provided strong support to the view that Med14 plays an essential architectural role in the complex, helping to stabilise the intermodular interactions that hold the complex together. Recent high-resolution structural studies have since confirmed the localisation of Med14 first provided by integrative modeling and have provided further structural evidence for the role of Med14 to bridge the Head, Middle and Tail modules [75–77].
Likewise, Med17, the largest subunit in the Head module had been shown to act as a scaffold by forming extensive interactions with nearly all of the remaining Head module subunits. Despite playing an essential functional role in yeast (P. Robinson unpublished results), the first 181 N-terminal residues of Med17 were unstructured and unresolved in X-ray studies of the Mediator Head module [34,58,78]. The ISD approach provided an explanation for this missing Med17 crystallographic density, which was found to play a critical architectural role in the Mediator complex. The Med17 N-terminal domain (NTD) was found to extend away from the surface of the Head module to interact closely with a central region of the Mediator Middle module (Fig. 2A). The Med17 NTD was mapped to this region through cross-links to the Middle module subunits Med 7, 21, 4, 9, 14, and could be localised with a precision of roughly 10 Å. By extending away from the Head module to take part in an extensive interaction network with these Middle module subunits, the Med17 NTD was predicted to play an important bridging role between these two distinct Mediator modules. As with Med14, later high-resolution studies confirmed the unique bridging interactions played by Med17 within the Mediator complex [75–77] and demonstrated that the ISD approach provided an accurate low-resolution model for the N-terminal region of Med17 within the Mediator Middle module guided by a collection of CLMS distance restraints.

### 3.4. Middle module architecture

The fact that over half of the subunits comprising the 9-subunit, ~300KDa Middle module are essential for yeast survival corroborates the idea that this Mediator region plays a central role in gene activation. As for the conserved Head module, attempts to understand the structural basis of its core biological role were initiated through the purification and X-ray crystal structure solution of small stable domains such as the Med7N-Med31 and Med7C-Med21 heterodimers [61,62]. However, these models accounted for only a small fraction of the primary sequence of the Middle module and in isolation they provided little biological insight. Attempts to build a more complete picture of the architecture of the Middle module came first from cross-linking mass-spectrometry [79] and later from affinity labelling studies where subunits were coarsely mapped onto 2D negative-stain EM images of Mediator using labelled antibodies or based upon difference maps following subunit deletion [70,71]. Although both of these experimental approaches provided valuable new information, neither were sufficient to provide a complete 3D architectural view of the Middle module. These studies had, however, made some interesting predictions as to the arrangement of Middle module subunits. Based upon: (1) cross-link constraints, (2) evolutionary relationships between the Med7-21 and Med4-9 heterodimers, and (3) end-to-end helical packing interactions observed within crystals of Med7C-Med21, Lariviere and colleagues predicted that an extended arrangement of end-to-end stacked helical bundles would form a central scaffold within the Middle module [79]. Later EM observations corroborated the idea of an extended Middle module architecture but were unable to resolve the internal subunit organisation [70,71]. Although tantalizing, these early modeling efforts were manually derived, speculative in nature and remained unchallenged by rigorous computational validation procedures.

A first complete 3D architectural representation of the Middle module was produced using ISD (Fig. 2A) [17]. Crucially, this approach differed from earlier attempts to define Middle module architecture by being able to combine all available structural information, comprising atomic models, comparative models, cross-link restraints and a low-resolution EM envelope in an unbiased computational approach to exhaustively sample configurational space. Strikingly, the computational approach confirmed that the elongated Middle module density was defined by an end-to-end packing of the 4-helix Med7C-21 and Med4N-9 heterodimers. This central scaffold was localised within the EM density at a high level of precision, with an average root means square fluctuation (RMSF) of around 10 Å within the cluster of top-scoring

model solutions. The architectural model provided details of the arrangement of Med10 and Med19 at one extreme of the elongated Middle module and their binding interface with the 4-helix bundle of Med7C-21. At the other extreme the large subunits Med14 and Med1 were shown to interact on opposite surfaces of the equivalent helical bundle formed by Med4N-9. As described above, the model also provided details of an extensive network of interactions positioned close to the central Med7C-21/Med4N-9 junction that involves not only these four subunits but also the N-terminal extension of the Head module subunit Med17 and central portions of Med14. Such details provided by the architectural model of the Middle module produced by integrative modeling have since been validated by high resolution structural studies using X-ray crystallography and cryo-EM (Fig. 2A) [75–77]. Such validation confirms that integrative modeling techniques can produce a very accurate description of the internal 3D organisation of multi-subunit assemblies and form a strong foundation from which to design further structural and functional analyses.

### 3.5. TFIIK location in the pre-initiation complex

TFIIK is a trimeric subassembly within the large multi-subunit general transcription factor TFIIH. Whereas the remainder of TFIIH functions as an ATP-dependent DNA helicase/translocase with roles in promoter DNA opening [80] and DNA damage repair [81], the TFIIK trimer performs a disparate role as a protein kinase. Yeast TFIIK is composed of the cyclin-dependent kinase Kin28, its cyclin Ccl1 and a third protein Tfb3, which together phosphorylate the C-terminal domain of the largest RNAPII subunit Rpb1 [82]. Phosphorylation of Serine 5 residues within the highly conserved tandem CTD heptapeptide repeats correlates with disruption of the Mediator-RNAPII interaction at the gene promoter [37] and transition of RNAPII to the elongation phase of the transcription cycle [38,39]. Therefore, in the context of the promoter-associated pre-initiation complex, TFIIH simultaneously makes two important contributions to transcriptional initiation: (1) ATP-dependent DNA translocase activity to melt the duplex promoter DNA and allow engagement of the template strand in the active cleft of RNAPII and (2) destabilisation of the PIC through CTD phosphorylation leading to RNAPII promoter escape. Understanding the structural basis of these two functions and how they are coordinated in the context of the pre-initiation complex has been a key focus of a number of recent structural studies. The architecture of the 31-subunit PIC, lacking a Mediator, was first determined using *in vitro* reconstitution and cryo-EM analysis of the human [83], and subsequently the yeast [84] promoter assemblies. These studies demonstrated that the PIC has a bilobal structure with one lobe comprised of RNAPII with a number of the general transcription factors including TBP, TFIIA, -B, and -F and the second comprised of TFIIH and regions of TFIIE. In these studies, promoter DNA could be traced through the structure, with large distortions introduced through interactions with TFIIB and TBP, such that prior to DNA melting, the double stranded DNA is directed along the active site cleft of RNAPII before exiting along a path that allows downstream contacts with the Ssl2 translocase subunit of TFIIH. Negative stain EM with affinity labelling [85] as well as ISD analyses [16] had previously determined the course arrangement of subunits within the yeast TFIIH complex. These analyses suggested that the two TFIIH ATPase subunits Ssl2 and Rad3 are associated at opposite ends of a central multi-subunit core (comprising Tfb1/2/4/5 & Ssl1), with the TFIIK trimer at a more peripheral location, adjacent to Rad3. This prior information helped to interpret TFIIH density in the PIC cryo-EM maps and position the core and associated Ssl2 and Rad3 ATPase subunits. Ssl2, the TFIIH translocase with a prominent role in transcriptional initiation, was positioned to interact with the downstream promoter DNA, and at the other end of the TFIIH core a Rad3 homology model could be unambiguously docked into corresponding EM density [84]. However, no density was found corresponding to the ~120KDa TFIIK trimer in the vicinity of Rad3 or otherwise, and
therefore TFIIK was assumed to be highly flexible and averaged out
during image processing.

The location of TFIIK within the PIC was determined using an in-
tegrative structural biology approach in which cryo-EM was combined
with CLMS and ISD to calculate and interpret EM density for an *in vitro*
assembled 52-protein Mediator-bound pre-initiation complex (Fig. 2C)
[18]. Mediator and PIC EM densities could be docked into the Med-PIC
map with little evidence of structural perturbation upon complex for-
mation (Fig. 2A). However, such docking highlighted a region of PIC
density immediately adjacent to Rad3 that had not been observed in the
analysis of PIC complexes lacking Mediator. As the Mediator complex had
earlier been shown to stimulate the CTD kinase activity of TFIIK [27], and
since TFIIK was missing in the earlier PIC reconstructions [83,84], the
TFIIK trimer seemed an obvious candidate for the extra PIC density,
which appeared to become ordered upon Mediator binding. An ISD ap-
proach was used to test whether this extra density, segmented from the
Med-PIC EM map, was consistent with TFIIK cross-link restraints
(Fig. 2D). The analysis was performed by fixing the adjacent Rad3
homology model at its corresponding position in the Med-PIC EM map as
an anchor and then searching configuration space for solutions for TFIIK
subunit localisations that best satisfy both TFIIK cross-links and EM
density restraints. Importantly, all top scoring models from these trial
simulations represented a single self-consistent TFIIK subunit architecture
that localised to the extra EM density whilst fully satisfying TFIIK cross-
links. The architectural solution demonstrated that direct interactions
between Rad3 and Tfb3 likely underpin the connectivity of TFIIK to the
TFIIH core whilst proximity to subunits in the Mediator Middle and Head
modules likely play a role in positioning TFIIK in space. Localising TFIIK
within the PIC provided a clear explanation of the role of Mediator in
stimulating the CTD kinase activity of TFIIK. Unbiased docking of the
Head-CTD crystal structure within the PIC shows that Head interactions
position the CTD in the immediate vicinity of TFIIK. Therefore, the
structural evidence strongly suggests that Mediator acts to simultaneously
position both the RNAPII CTD and TFIIK trimer within the PIC and bring
them into close proximity in order to stimulate CTD kinase activity.

3.6. *Tail module architecture, dynamics and activator interactions*

The Tail module is composed of non-essential protein subunits with
higher levels of sequence divergence than observed within the Head and
Middle modules. Rather than contributing to the basal transcription in-
ition mechanism the Tail module is instead responsible for mediating
the regulatory response to sequence-specific DNA-binding transcription
factors [66]. Indeed, the yeast Tail module Med2-Med3-Med15 triad has
been shown to be a common target of transcriptional activator proteins
including the extensively studied Gal4 and Gcn4 [68,86]. Studies ex-
ploring the nature of the Gcn4-Mediator interaction have described a
"fuzzy" binding mechanism whereby numerous low-specificity and low-
affinity hydrophobic interactions function additively to amplify the
overall binding affinity between activator and Mediator [43]. A recent
study combining cross-linking with NMR demonstrated that interactions
between the Gcn4 activation domains (ADs) and multiple activator-
binding domains (ABDs) within Med15 are very heterogenous in nature
with all possible AD-ABD combinations detected and contributing to a
large "fuzzy" interaction surface [44]. The generality of this binding
mechanism is still an open question but it is interesting to note that
transcriptional coactivators often harbour multiple ABDs, and such a
model could potentially explain how a collection of otherwise unrelated
transcriptional activator proteins can converge on a limited number of
transcriptional coactivator targets. A further open question is how the
transcriptional response is modulated by transcription factor-Mediator
interactions. Do such interactions simply increase the residency of Med-
PIC complexes at the gene promoter, as some have suggested [87], or
does transcription factor binding have an allosteric effect, leading to
conformational changes that modulate transcription levels? Evidence for
the latter has come from EM studies in which activator binding events

have been correlated with structural changes within the Mediator com-
plex [88]. To date, the Tail module remains the least well characterised
of the three Mediator modules. For example, high-resolution structural
studies of the Tail have been limited to a modest collection of small
protein domains [43,89]. Low-resolution mapping of subunit termini
within EM projection images provided some insight into the relative
location of the various Tail subunits, but this information was restricted
to two dimensions [70,71].

The first 3D architectural details of the Tail module were provided alongside those for the Middle module in an ISD study describing the subunit architecture of the full Mediator complex [17]. The model that resulted from this study localised each Tail module subunit within the 3D EM envelope and provided a structural explanation for the EM labelling studies, protein-protein interaction data and sub-complex isolation observations reported previously [49,54,64,65,67–69]. In the case of the Tail module, the architectural solution was strongly influenced by the inclusion of a 540-residue homology model for the N-terminal portion of Med16, identified by strong sequence similarities to the 7-bladed β-propeller of the *S. cerevisiae* vesicle coat protein Sec31. Positioning of this homology model through cross-linking and EM constraints helped to determine the proper location of all other Tail module components. Specifically, Med2, Med3 and Med15 that make up a Tail module triad commonly targeted by activators, were found to colocalize with the C-terminal region of Med14 at the junction between the Tail and Middle modules. The N-terminal portions of Med2 and Med3, both of which are strongly predicted to form a coiled-coil motif, were colocalized with particularly high precision to a single region of the Tail (Fig. 3A), supporting the idea of coiled-coil interactions and explaining the observation of Med2-3-15 sub-complex isolation in earlier biochemical studies [68]. Likewise, the observed Med16 β-propeller-Med5 interface explained earlier reports of Med5-Med16 N-term sub-complex isolation [70].

Information on the architecture and dynamics of the Tail module was expanded in the integrative structural studies of Med-PIC, which revealed some surprising novel features of Tail module biology (Fig. 2E). Firstly, all prior structural studies on the Mediator complex had suggested that only a single inter-module connection existed between the Middle and Tail modules, for which the C-terminal domain of Med14 was shown to be a critical component [70,71]. However, evidence from cryo-EM and cross-linking on Mediator-bound to the PIC showed that in fact a second Middle-Tail connection existed and was brought about by a direct interaction between Med1 and Med5 from the Middle and Tail modules, respectively (Fig. 2E). Such a direct interaction had been picked up in yeast two-hybrid screens but overlooked due to a lack of supporting structural evidence [49]. Second, the Tail module was observed in two distinct structural states. The most populated and best resolved state, termed the “Tail Down” state (Fig. 2C) has the Med2-Med3-Med15 triad positioned a significant distance from the upstream transcriptional activator sequences and out of range to form cross-links to other components of the transcriptional apparatus. However, a second lower populated state, termed the “Tail Up” state, involves a 50 Å movement of the whole ~450KDa Tail module toward the upstream promoter DNA bringing the Tail in close proximity to the back side of RNAPII (Fig. 2E). Cross-linking analysis on Med-PIC, which in contrast to the cryo-EM, was conducted in the presence of the transcriptional activator Gcn4, corroborated EM observations of a “Tail Up” conformation with numerous Tail-RNAPII and Tail-Gcn4 cross-links observed in the CLMS dataset. Such cross-links were only consistent with a Tail module that had transitioned from the “Tail-down” to “Tail-up” conformations, which raised the tantalising possibility that activator interactions may be involved in stabilising such a structural transition. A comparison of the Tail-Gcn4 cross-link pattern observed either within the Med-PIC complex or for isolated Gcn4 and Med15 domains, highlights some interesting points about the mechanism of Gcn4 activation at gene promoters. The Gcn4-Med15 interaction is driven by hydrophobic interactions between residues in the intrinsically disordered N-terminal activation domain of Gcn4 and hydrophobic
surfaces found within four Med15 regions, an N-terminal KIX domain and three ABDs (Fig. 2E) [90]. Biochemical studies have scrutinised these interactions at the primary sequence level [91] and a solution structure of Med15 ABD1 bound by residues from the Gcn4 activation domain demonstrated that Gcn4 residues adopt a helical conformation upon binding to allow aromatic residues to insert into a shallow hydrophobic Med15 cleft [43]. Interactions between residues within the Gcn4 activation domain and Med15 KIX domain and ABDs, have been consistently observed both in the context of the Med-PIC complex as well as more simple systems limited to just the Gcn4 activation domain and truncated Med15 constructs [44]. However, cross-link analysis within the Med-PIC assembly portrays a more complex picture for

Gcn4-Mediator interactions. Specifically, the Gcn4-Mediator cross-link with the highest overall spectral count was not between a Gcn4 residue within the N-terminal activation domain, but instead between a C-terminal Gcn4 residue located at the terminus of the DNA-binding bZIP domain. Importantly, within Mediator this cross-link mapped to Med15 ABD2, which was also found to be heavily cross-linked to Gcn4 residues within the N-terminal activation domain, as expected. Indeed, a striking picture emerged where residues from both the Gcn4 activation and bZIP domains formed an overlapping pattern of cross-links within the KIX and ABDs of Med15 (Fig. 2E). No such pattern has been observed in studies probing the interaction of defined Gcn4 and Med15 domains simply because the Gcn4 bZIP domain was not present [44]. One

A

<figure><img src="image_5.png" /></figure>

B

<figure><img src="image_6.png" /></figure>

C

<figure><img src="image_7.png" /></figure>

D

<figure><img src="image_8.png" /></figure>

(caption on next page)
**Fig. 3.** Cross-linking MS guided integrative structure determination. (A) Cross-linking coverage plays a key role in determining how precisely a model component can be localized. Med2 and Med3 both have numerous inter-protein cross-links constraining their N-terminal regions in the Mediator model (represented on bottom panel by grey lines on blue schematics), yet none were found after the first 130 residues. Coarse grained beads representing 40 residues of Med 2 and 3 are localized in the final model with RMSF precision between 10 and 20 Å from the N-terminal region (blue shades in upper heat map representation). Beyond this, the precision rapidly deteriorates (red regions of heat map). (B) Cross-linked peptides can only be measured by a mass spectrometer if the precursor ion is successfully selected for dissociation and product ion analysis. Tryptic digests of large, purified protein assemblies contain hundreds of peptide ions co-eluting from the nano-LC column at any given moment even after fractionation to enrich for cross-linked peptides. The upper panel shows the total ion current corresponding to a single fraction of a size exclusion chromatography-based cross-link enrichment. The precursor ion scan (MS1) taken at 59 min shows hundreds of ion signals co-eluting. Cross-link identification requires isolation of the quintuple charged precursor ion at 532.31 m/z for gas phase dissociation and measurement of the product ion (MS2) spectrum. Factors influencing whether the cross-link will be identified include: the scan rate of the mass spectrometer, the extent of enrichment/fractionation in the 1st- chromatographic dimension, the peak capacity and gradient length in the 2nd-chromatographic dimension, and the ionization efficiency of the cross-linked peptides and their intensity relative to co-eluting peptides. (C) After measuring a cross-linked product ion spectrum, successful identification requires sufficient product ion formation to identify both members of the cross-linked peptide pair. The upper and lower panel show the best and second-best hits to a cross-linked product ion spectrum after a database search. The correct hit, corresponding to an intra-protein cross-link on Rpb1, differs from the incorrect hit, a potentially more interesting inter-modular cross-link between Rpb1 and Med1, by a single blue y3-ion. In both cases, almost all of the product ion signals are matched to the sequence (green lines), with the longer, red peptide accounting for most of these. (D) Mapping cross-links to existing high-resolution structure provides one method of validating a dataset. Common reagents such as BS3 and DSS leave a suberate bridge between adducted lysine residues. To account for the flexibility of the lysine side-chain, distances are measured between Cα, giving an expected span of about 25 Å. In practice, 30 or 35 Å is typically used as the violation distance to account for protein dynamics and uncertainty in the structural coordinates. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)

intriguing interpretation is that the Gcn4 activation and bZIP domains are colocalised at the Med15 binding interface, which could only result from some unexpected extension of the Gcn4 N-terminal activation domain away from the DNA towards the C-terminal end of the bZIP coiled-coil domain. Further structural studies of Med-PIC with Gcn4 bound upstream will be required to further address such mechanistic questions. Although the “fuzzy” nature of the Gcn4-Med15 interface would make high-resolution studies of Gcn4-Tail interactions extremely challenging, many important questions could be addressed by visualising the interaction even at more modest resolutions. For example, such reconstructions could help address whether bending of the upstream DNA is required to bring bound activators in range of

interactions with the Tail module. Even in the “Tail-Up” state, Med15 was found to be spaced roughly 80 Å from the site of Gcn4 DNA binding meaning that DNA bending or further Tail motion could well be involved in bringing Gcn4 and Med15 together.

# 4. Practical approaches to integrative modeling of the transcription initiation machinery

## 4.1. Cross-link mass spectrometry and the modularity of the proteome

While integrative structure determination (ISD) aims to utilize as much experimental data as possible, in recent practice, CLMS, alongside

<figure><img src="image_9.png" /></figure>

**Fig. 4.** Exploiting structural modularity and data transferability. The assembly is a stable aggregate of individual macromolecules (circles), organized into modules (blue, yellow, green and coral colours). (A) The architectures of two functional states of the same assembly are related by a global conformational change which rearrange the modules without significant perturbation of their internal structure. The two functional states are for instance obtained by varying the composition of the two samples. Functional state A is used to derive a low-resolution 3D density map, by cryo-EM experiments (grey circles), while functional state B is utilized to collect CLMS data, which can be grouped into intra- and inter-module cross-links (coloured lines). (B) If the target of the modeling is the structural determination of the functional state A, the different proteins and components belonging to the modules are attributed to the corresponding regions of the 3D density map, after a suitable segmentation of the density (black lines). (C) To allow the violation of the inter-module distance restraints, two different unknown weights are attributed to the inter- and intra-module cross-links (dashed and continuous lines, respectively), which were collected for functional state B and transferred to functional state A. The weights are determined during the modeling, using the Bayesian scoring function. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)
approximately 5 and 30 residues and bearing some positive charge are ideally detected by reverse-phase chromatography coupled online to positive mode electrospray ionization and fragmented by collisional dissociation. Hence, trypsin has always been the ideal protease for proteomics experiments as the distribution of lysine and arginine residues generates favourable peptides. The use of other enzymes, such as chymotrypsin, Asp-N, or Glu-C may provide some complementary cross-link coverage for difficult to access regions, but in general these enzymes do not provide such favourable peptides for mass spectrometry analysis and comparative studies have not seen much benefit [96]. Namely, these enzymes tend to create peptides that are too short and are hard to identify unambiguously, or else peptides with little positive charge that do not ionize well. There is additionally a nice interplay between trypsin selectivity and lysine-directed cross-linking, in that modification at lysine by the cross-linking molecule blocks trypsin cleavage, ensuring that most cross-linked peptides will only have one likely site of modification. A mass spectrometer can also be tuned for peptides that do not match the criteria listed above. For instance negative ion mode can be used for acidic peptides, or mobile phases can be altered to favour hydrophobic or membrane spanning peptides. In practice, most mass spectrometry facilities are highly oversubscribed and changing the operating procedures of a major instrument places a burden on the available resources.

Third, MS-proteomics are typically conducted in data-dependent acquisition (DDA) mode. There are typically many more peptide ions eluting at a given time from the LC-column than can be isolated, dissociated, and measured (Fig. 3B). Therefore, DDA mode programs the instrument to select ions based on certain criteria, such as ion intensity and charge state. In some cases, DDA mode can select for the presence of a pattern of ions that is characteristic for cross-linked peptides, such as a pair of ions separated by a defined mass difference due to the use of isotopically heavy and light versions of a cross-linking reagent [115]. Of the four aspects required to detect a cross-link, this is the most easily modified or optimized. Cross-linked peptides can be enriched by size-exclusion chromatography [96], cation exchange [95], or by affinity tags in the reagent [116] or the protein complex of interest [24]. Furthermore, the reverse phase chromatography that is directly coupled to the mass spectrometer can be optimized for longer gradients using longer packed or monolithic columns [18]. DDA analysis for cross-linked peptides typically programs mass spectrometers to exclude singly and doubly charged precursors to focus on cross-linked peptides which tend to be quadruply charged or higher. Cross-linking chemistries can be employed that increase the charge state further still [117,118] for improved targeting of cross-linked versus linear peptides. An exciting development is the recent increase in coupling of ion mobility separation (IMS) with mass spectrometry [119,120]. IMS provides a millisecond scale separation of ions in the gas phase based on collisional cross-section and hence charge-state and promises to provide a route to further enrichment and isolation of cross-linked peptides.

Finally, cross-linked peptides, after being isolated for dissociation, must produce enough product ions such that search algorithms can identify both of the cross-linked peptides [100]. Collisional based dissociation methods produce fragment ions according to mobile proton theory [121] and whether product ions from both peptides are observed largely depends on the gas-phase basicities and charge distributions of the peptides. In cross-linking-MS, a common situation is that the product ion distribution identifies one peptide very well, but leaves the second peptide poorly identified. Gas-phase cleavable cross-linking reagents circumvent this problem by containing low-energy bonds that dissociate prior to the peptide backbone during MS2. Subsequent isolation of the individual peptide ions and a second round of dissociation in MS3 provide fragmentation information for each peptide [102,103]. This approach provides both the individual precursor masses for each peptide as well as their fragmentation patterns and can aid in correctly identifying the peptides in the cross-link. In practice, these schemes have inherent trade-offs such as an increased mass spectrometry cycle

time resulting in fewer cross-linked precursors being sequenced, as well as insufficiently specific chemistry in MS2 leading to inefficiencies in targeting the correct precursor ions for MS3. Electron-transfer dissociation (ETD) is an alternative to collisional activation which shows less specificity to the physicochemical properties of the peptide analytes [122]. Efficient electron transfer to begin the dissociation process requires positively charged peptide substrates, so that ETD should be the ideal method for sequencing of positively charged, cross-linked peptides to give an equal distribution of product ions from both peptides. This method was infrequently applied however until recently, because ETD implementation has consistently lacked in speed and sensitivity to collisional methods. However, newer implementations [123] are changing this situation and ETD is likely to become the method of choice for cross-linking analysis.

### 4.3. "Noise" in cross-linking data sets and implications for ISD

In practice, cross-links are never identified completely unambiguously. For any non-trivial-sized search space, a given spectrum will match to multiple theoretical cross-linked peptide pairs that are within the *m/z* tolerance of the precursor ion. These matches need to be discriminated on the basis of how well the experimentally observed product ions are explained by the theoretical cross-linked peptides. Typically, the matched product ion series are incomplete, and the experimental spectra are noisy and contain product ion signals that come from either co-eluting peptides with similar precursor *m/z* or otherwise cannot be explained by the cross-linked peptide sequence. Especially for short peptide sequences that are four or five residues long, multiple proteins in a database might contain either identical or isomeric versions of these peptides that differ only in the arrangement of amino acid residues. In these cases, the correct cross-linked peptide match can sometimes differ from the next best match by only a single additional product ion match (Fig. 3C).

Various metrics describing the match of the theoretical cross-linked sequence to the experimental spectrum are reported by proteomics search engines along with an overall statistically-derived classification score. Hence, like other MS proteomics datasets, cross-linking data are expected to contain misidentifications. Datasets are reported alongside a false discovery rate (FDR), which is commonly based on conducting the search against both the target protein sequences and randomized, decoy versions of the same proteins [124]. For this approach to accurately estimate the chance of a false match, a large enough decoy database must be used to reliably model the distribution of incorrect hits. Hence, the typical proteomics approach of using a decoy database equal in size to the target database is not appropriate for small database searches containing a limited number of protein sequences. In our Mediator studies, we used a decoy database that was 10x larger than the target database, making the total number of protein sequences searched ~700.

Additional validation of cross-linking data comes from mapping the identified cross-links onto the atomic structures (when available) for components of the complex. The distance between Cα atoms of the cross-linked residues should be consistent with the length of the cross-linker bridge group plus the length of the adducted amino acid side chains. An additional 5–10 Å are typically added on top of this to account for imprecision in the atomic coordinates and dynamic motions in the protein complex. Hence BS3 or DSS cross-links are generally considered “violations” when the Cα-distance exceeds a limit of 30 or 35 Å, which includes the suberate bridge (~12 Å), the length of two lysine side chains (~12 Å), plus an addition 5 or 10 Å (Fig. 3D). The violation rate for a cross-linking dataset is often similar in value to the FDR. However, the FDR formally describes the proportion of reported residue pairs that are incorrectly assigned, whereas there are a number of other reasons why violations may occur. These involve deviations from the expected crystal structure either due to a conformational or configurational discrepancy between the sample being cross-linked and
the structural model used for assessment. Therefore, a violation rate that is substantially higher than the FDR can indicate either a discrepancy in the conditions used to obtain the different types of structural information, or a broader problem such as protein aggregation under the cross-linking conditions.

With respect to utilizing cross-linking derived restraints in integrative modeling procedures it is important to realize that cross-linking data are expected to contain incorrect or otherwise inapplicable restraints, and that the probability of an individual restraint being incorrect is reflected in a classification score, which can be passed on to the modeling parameters. In this way, cross-links with more certain assignments can be weighted more strongly than others. One aspect of experimental mass spectrometry data that can be taken advantage of is that there might be multiple lines of evidence supporting the assignment of some cross-links. A CLMS experiment generates between thousands and hundreds of thousands of product ion spectra that are all searched for potential cross-links. Spectra that match above some classification threshold are termed Cross-linked Spectral Matches (CSMs). CSMs match a cross-linked peptide pair, and the cross-linked peptide pair identifies a cross-linked residue pair ("cross-link"). Hence a given cross-link is often identified by multiple cross-linked peptide pairs, which are identified by multiple CSMs. It is important to note that this redundancy derives not only from picking the same precursor multiple times over one or more experiments, but from redundancy in the chemical entity identifying the cross-linked residues. For instance, cross-linked peptides can be identified from differently charged precursor ions, or one or both of the peptides might differ with respect to the length of the peptide (due to missed enzymatic cleavage), with oxidation state at methionine residues, or due to the presence of a "dead-end" modification (in addition to a missed cleavage). These factors will produce different precursor ions at different *m/z* values with different product ion spectra. Hence, cross-links identified from multiple cross-linked peptides have multiple lines of evidence pointing to their existence. Cross-links identified from multiple CSMs however, are not necessarily more confidently identified as the search algorithms are deterministic so that incorrectly identified CSMs will be identified consistently leading to redundant incorrect CSMs.

### 4.4. Treatment of the input data for ISD

The first stage of the ISD approach consists of gathering all the available data needed to build the model and validate it. Spatial information about a given system can include data from a variety of experiments, statistical propensities extracted from known homolog sequences or protein structures, and physical laws, such as atomic interactions obtained from a molecular mechanics force field [46]. In our approach input data include the sequences of the subunits and their stoichiometry, the available atomic structures of domains, subunits or subcomplexes, CLMS data with the identification frequency and identification scores of the cross-links, low-resolution cryo-EM density maps of the complex, EM-labeling or EM-deletion experiments, and protein-protein interactions from affinity purification and yeast two-hybrid assays. All the input information is used to build the representation, the scoring function, and the validation toolbox for the system. First, crystallographic structures of the subunits and homology models are used to build the representation. Second, the cross-linking data and the EM density data, the sequence connectivity and physio-chemical data are used to build the scoring function, third, data which are difficult to encode in a scoring function are left out of the modeling procedure and used for validating the final model.

For integrative modeling of the 21-subunit Mediator complex [17], which will form the basis of discussions in Sections 4.4–4.9, the data included:

(1) 298 Mediator cross-links identified from 1900 spectra, consisting of multiple peptide sequences, charge states, or replicates.

(2) The atomic structures, including X-ray crystallography structures and homology models, covering 80% of the Mediator Head module, and overall 23% of the residues in the Mediator complex. It is important to note that whilst the resolution of the output models is influenced by the coverage of atomic structures, a large percentage of atomic structure is not required to produce a descriptive ISD model. Early ISD models for the nuclear pore complex did not contain any atomic structures [125] and the coverage for the Mediator [17] and SEA complexes (35% coverage) [126] were both very low.

(3) An 18 Å resolution cryo-EM density map of the Mediator complex was used to constrain the overall-shape of the complex.

(4) Validation data, such as yeast two-hybrid, immuno-precipitation, subcomplex isolation assays, and protein localization from labeling and domain deletion EM studies. The validation data was not directly used in the modeling, but a *posteriori*, to assess the quality of the models.

All the modelling protocols were scripted using the Python Modelling Interface (https://github.com/salilab/pmi) [15], a library to model macromolecular complexes based on the open source IMP package (http://salilab.org/imp/) [46]. Files containing the input data, scripts, and output models are available at a free repository, see for instance http://salilab.org/mediator and https://zenodo.org/record/802915.

### 4.5. System representation and scoring function development for ISD

The computational representation of a macromolecular system is the collection of all the degrees of freedom that need to be determined based on input information [14]. The representation assigns the system components (e.g. residues, atoms or domains) to geometric primitives with 3D coordinates (e.g. points, spheres or 3D gaussians) needed to compute the score (Fig. 1). This assignment is decided before any other computations, and is fixed along the computation. Recently, protocols to optimize the molecular representation have been introduced [127], where the optimum is the representation with the highest resolution for which sampling is still exhaustive and the models have a precision that is commensurate with the precision of the representation. The radius of the spheres and the variance of the Gaussians were set to describe the average molecular volume and the molecular electron density of poly-peptide segments, respectively. The representation also uses rigid bodies, where the relative positions of the primitives (e.g. residues in a domain) are constrained based on a crystallographic structure. We use coarse-grained beads to encode the uncertainty related to the degrees of freedom of protein regions which are not directly represented by atmospheric structures in the input data. These coarse-grained beads are free to move as spheres on a string, as opposed to the rigid-bodies, and are connected based on the amino acid sequence. The resolution of the coarse graining is commensurate with the resolution of the input information, in our case the cross-links and the cryo-EM density map. As a consequence, the representation is multi-scale, where different parts of a structure are represented by different coarse-graining levels. This strategy maximizes computational efficiency while avoiding model oversimplification.

In our modeling approach we used three scales to represent the system where the rigid-bodies are defined: the first is where each bead corresponded to individual residues, and was centered at the position of the $C_α$ atom, the second is where each bead represented 10-residue segments and was positioned on the center of mass of all atoms of the corresponding segment, and the third scale is the Gaussian mixture model (GMM) approximation of the electron density of the corresponding structure [128]. We adopted a two-scale representation for the flexible strings: the first is where each bead represented either a 20- or 40-residue segment and was positioned on the center of mass of all atoms of the corresponding segment. The second is a spherical
Gaussian. In this case the bead and the Gaussian centers were enforced to be identical.

The gathered data and the representation are next used to build a scoring function, which ranks alternative models based on how well they quantitatively reproduce the input data [129]. In our approach, the scoring function computes: (1) the excluded volume to avoid clashes between distinct residues/domains, (2) the quality of the covalent connectivity of the polypeptide chains, (3) the agreement with both the cryo-EM density map and with the CLMS data. Also, the scoring function takes into account the uncertainty of the input information, and encodes the presence of modularity and whether some data was transferred from another macromolecular sample. In our approach we prefer a Bayesian formulation of the scoring function [130], which is more objective than traditional scoring functions. This formulation can infer unknown quantities, combine different types of information and account for noise in the data. This is particularly important for the treatment of CLMS data, were some assignments might be erroneous.

The Bayesian scoring function is proportional to the posterior probability $p(M|D,I)$ of a model M given the gathered data D and prior information I. The model $M = (X, \{\alpha_i\})$ includes the structural coordinates $X$ and additional parameters $\{\alpha_i\}$ which describe unknown attributes of the data such as the noise or the relative weight of distinct data pieces. Using Bayes' theorem, the posterior probability is $p(M|D,I) \propto p(D|M,I)p(M,I)$, where the likelihood function $p(D|M,I)$ is the probability of observing data $D$, given $I$ and $M$, and the prior is the probability of model $M$, given $I$. The likelihood function for the CLMS data was built from a “forward model” that predicts the formation of a cross-link given the coordinates of a structural model, and a “noise model” that quantifies how much we can tolerate an error of the model (e.g. a cross-link that was observed and reported in the dataset which is not predicted to be formed in the model) [17,24,101]. Finally, the joint likelihood function $p(D|M,I)$ for a dataset $D = \{d_n\}$ of $N_{XL}$ independently observed cross-links is the product of likelihood functions for each cross-link. Both the functions of the forward and the noise model introduce parameters to describe uncertainty. In particular, the likelihood function uses $\psi$, which is the unknown uncertainty that a cross-link is correctly assigned and regulates the weight of the corresponding restraint.

Importantly, the CLMS data was collected on a sample of Med-RNAPII Holoenzyme and transferred to the Mediator system lacking RNAPII. The addition of RNA pol II is observed to induce motions at the junctions between Mediator modules, while the intra-module topology appears unchanged [48,131]. As a consequence, we expected the inter-module cross-links to be less accurate than the intra-module cross-links in describing the apo-Mediator state. To account for transferability and modularity, we assigned the inter- and intra-module cross-links, respectively 8% and 92% of the total, to different classes that were fit with different values of $\psi$ (Fig. 4C). In this way we lessen the impact of inter-module cross-links while benefitting from intra-module ones.

Finally, the scoring function for the cryo-EM data was computed based on the Gaussian Mixture Model (GMM) representation of each domain, using the cross-correlation coefficient between GMM representations of the EM volume and model components [17]. The weights of each GMM component were normalized to the relative mass of the component vs the mass of the module. We exploited the modularity by dividing the EM density map into segments reflecting the Head, Middle and Tail module on the basis of previous approximate localization of subunits within two distinct regions of the EM map [70,71] (Fig. 4B). During modeling, each Middle or Tail subunit density was restrained to the corresponding Middle or Tail EM density segments. The Head module was kept fixed in its initial docked position, except for the sections that were not resolved in the crystal structure and therefore present as coarse-grained representations. Recently, a Bayesian scoring function based on the GMM approach was developed, allowing an objective weighting of the EM density data with respect to the other information [128].

## 4.6. Exhaustive configurational sampling with replica exchange

One important requirement for the ISD approach is the generation of an exhaustive sampling of the model. To achieve this goal, we have to adopt improved sampling schemes and carry out a large number of independent sampling runs. In each run the positions and orientations of rigid bodies and flexible strings of beads, and the values of free Bayesian parameters, are randomly and iteratively perturbed in an effort to satisfy and optimize the scoring function, based on the Metropolis Monte Carlo algorithm. The improved sampling protocol followed the Replica Exchange scheme [130], using 64 replicas, with temperatures ranging between 1.0 and 2.5. The Replica Exchange scheme enhances the sampling by allowing hopping between alternative minima of the scoring function, thus preventing the sampling from becoming stuck in local minima. To confirm that we had sampled conformational space sufficiently to reach model convergence, we compared two independent halves of the solutions to each other and to the entire set, and assessed whether they display similar structural features. If we couldn't find a satisfactory agreement, we double the number of sampling runs, until convergence. For the Mediator modeling, the exhaustive sampling produced a total of 165,523 models in 20 independent runs. In a recent paper, an automated convergence protocol was proposed with a series of statistical tests [132], which are more objective and stringent.

## 4.7. Model assessment and clustering in ISD

Integrative structural determination eventually results in clusters of individual models that best satisfy the input constraints. To identify these clusters, analysis is carried out on a small fraction of the models representing the very best scoring solutions. First, the 500 best scoring models (the solutions) are selected from the whole sampled ensemble and checked for whether they satisfy the input restraints. Second, the solutions are clustered by structural similarity using the root mean squared distance (RMSD) of the beads as a structural metric. The number of clusters is set to be the most parsimonious number that separates the main structural differences between all the solutions. Finally, the representative cluster(s) are selected according to their score, or alternatively the agreement with data that was not used in the scoring function (the validation data).

In the case of the Mediator, the solutions were grouped into four main clusters, with distinct structural properties. Clusters 1, 2, and 4 had the Middle module oriented in the expected orientation (i.e., Med1 and Med19 were at the bottom and at the top of the module, respectively). The three clusters differ in the arrangement of the Tail subunits. Cluster 3 had the Middle module in a ‘flipped’ orientation (i.e., Med1 and Med19 are at the top and at the bottom of the module, respectively). All the solutions from all clusters satisfied excluded volume, sequence connectivity, and EM restraints. Since the CLMS dataset was transferred from the Holoenzyme sample, not all the cross-links were satisfied. In fact, while 95% of intra-module cross-links were fulfilled, only a small portion (10%) of inter-modular cross-links were satisfied. The Bayesian scoring function automatically down-weighted the inter-modular cross-links to satisfy the EM restraint, pointing to an inconsistency of the inter-modular cross-links in the Holoenzyme dataset with the EM density map of the free Mediator. The satisfaction of intra-modular cross-links suggests that the arrangement of Mediator subunits within each module of the Holoenzyme retains the main characteristics to that in the free Mediator. Since the average scores of the different clusters were only marginally different, it was not possible to rank the clusters and identify the configuration which best fit the input data. Therefore, we used the validation data, and in particular the EM localization experiments, to choose the representative cluster. Remarkably, only Cluster 1, which also happened to be the top-scoring cluster, was fully in agreement with all EM localization experiments.
## 4.8. Calculating precision of the solution, subunit density and subunit-proximities in ISD

The degree to which a component of the model, either a coarse-grained bead or residue, fluctuates across the ensemble of individual models indicates how precisely localized that component is in the modeling solution. We used several analysis metrics to represent the global and relative localization of domains and compute their precision: the localization density, the domain precision calculation, the root mean squared fluctuations (RMSF), and protein-protein contact analysis. First, for a given cluster of solutions, we computed the probability of finding any residue of a given protein at any point in space (i.e., the localization density map). The localization density is stored in a 3D grid, similarly to a 3D density map obtained from cryo-EM data, and it is a very convenient way to display the architecture of the whole complex from an ensemble of solutions. All localization density maps of proteins and domains are displayed in the software Chimera by an isosurface [133], which in our case was obtained using a threshold of 0.15. The precision of a domain (or a protein, or the whole complex) in a given cluster of solutions is calculated as the average RMSD between the cluster center and all other solutions. The precision can be used to compute the distance between two clusters, as the average RMSD on the whole complex, between every pair of solutions taken from the two clusters. The RMSF of a given residue is the standard deviation of the set of distances between the position of the residue in each solution of the cluster and its position in the cluster center.

As might be expected, there is a strong relationship between the number of consistent cross-links restraining a given bead and its localization precision. For instance, numerous cross-links were detected at the N-termini of Mediator Tail module subunits Med2 and Med3, with almost none at the C-termini. This was reflected in a RMSF near 10 Å for the N-terminal regions along with a high level of confidence in placing them at the junction between the Middle and Tail modules of Mediator (Fig. 3A). In contrast, after the first 100 residues or so it becomes difficult to localize these subunits. This intuitive finding points to the need to maximize the coverage of cross-links across the complex in question, particularly in regions with little or no atomistic structure available.

Protein-protein contact analysis of the solutions in a cluster can be carried out at the residue or domain level. In the former, the relative contact frequency of a pair of residues is computed based upon how often the two residues contact each other in the cluster, where a contact occurs when the distance between the residues is less than 10 Å. In the latter, two domains are in contact when the surface of any bead in one domain is within 10 Å of the surface of any bead in the second domain. Long sequences are divided into domains of 200 residues. Finally, the residue and domain interaction frequency can be rendered in a heat map that displays in a compact manner the propensity of contacts for a given cluster of solutions. Some of these contacts result simply from the identification of cross-links between the domains, however other contacts emerge indirectly due to the identification of synergistic sets of cross-links.

## 4.9. Model validation

The final and crucial stage of the whole ISD approach involves validing the models. Without validation, it is difficult to trust the models, as they can arise from artefacts of the modeling procedure itself, such as overfitting. First, the ensemble of solutions of the selected representative cluster is assessed in terms of how well it satisfies the data from which it was computed, including the cross-links, excluded volume, sequence connectivity and the three-dimensional EM restraints. This allows us to check the reliability of the modeling procedure: if models partially satisfy the data, one might suspect undersampling, strong data inconsistency, or mistakes in the data-encoding in terms of the scoring function [14]. As mentioned above, the Mediator model fit the EM map, satisfied most of the intra-modular cross-link

restraints, and it was stereochemically sound since it satisfied the excluded volume and sequence connectivity. Second, we used a cross-validation assessment using the CLMS dataset. In the cross-validation approach, we assessed the ensemble of solutions by comparing it with the ensemble of solutions obtained by jackknifing 10% the CLMS dataset in 34 different random ways. We compared the obtained localization maps, and computed the model precision to estimate the structural differences between clusters obtained using the whole CLMS dataset and clusters obtained using jackknifed datasets. Strikingly, the results were similar: both modeling approaches resulted in the same number of clusters, with the same structural features. That result suggested that the CLMS data was robust and accurate, and that the models were not a mere result of overfitting. Third, we validated the models based on unused data published in the literature. In fact, the Mediator model was remarkably consistent with almost all data from previous subunit interaction and subunit localization studies. The models explained subunit interactions inferred from co-expression [54,67,69], pulldowns/immunoprecipitation [54,68], and yeast two-hybrid assays [49,64,65] with only three discrepancies, which were clearly spurious. Furthermore, the Mediator model was validated by comparison with results from two-dimensional EM studies that used subunit labelling or subunit deletions to map subunit locations. [70,71]. The model was consistent with all these pieces of information, confirming the quality and the predictive power of the structure.

## 5. Conclusions and future perspectives on integrative structural techniques in the age of the Cryo-EM 'resolution revolution'

Despite the success of cross-linking and low-resolution EM-derived integrative models at inferring both global architectures and pseudoresidue level structures of protein machines in recent years, there is a seemingly diminishing role for these workflows in the current era where high-resolution EM reconstructions of large protein complexes are becoming increasingly prevalent. What is the value of an imprecise model if we can seemingly wait another half year before an atomistic structure is revealed? Is there room to improve the performance of cross-linking based models? Since the precision of these models depends upon both the distance defined by the cross-linking reagent as well as the density of cross-links that constrain a given amino acid residue three-dimensionally, it should be possible to achieve more refined integrative models using "zero-length", carbodiimide based cross-linkers [134] as well as by increasing the cross-linkable space with acid-directed [114] or other cross-linking chemistries. One approach to the latter aim is through using carbene generating, photochemical reagents [135]. However, the extremely heterogenous reaction mixtures that result from diazirine cross-linking make site-specific assignments of photo-cross-linked residues by mass spectrometry technically challenging and limit this approach. Recent attempts to utilize photo-cross-link derived restraints to aid protein structural prediction in the CASP11 competition showed no benefit [136]. However, technical developments using zero-length cross-linkers may still provide a future avenue to increasing the precision of cross-linking based structural inference. Such improvements may well be important for the structural analysis of protein complexes that continue to resist high-resolution EM reconstruction and thereby ensure that integrative modeling still has a role to play in the near future.

What factors determine the resistance of a protein complex to high resolution EM reconstruction? Even though cryo-EM has advanced significantly, heterogeneity derived from non-discrete, continual motions of domains still provides a significant challenge to algorithms designed to classify particles in 2- and 3-dimensions. This is especially the case when more than one region of a large multi-subunit complex displays continuous motion and the mobile regions represent a large proportion of the mass of the complex [137]. Many biological systems will be difficult to solve at atomic resolution for this reason and so the combination of low-resolution EM with X-ray crystallography and
cross-linking will still play an important role in these circumstances. Some features of the Mediator complex that were first revealed through our integrative model such as the connection between the N-terminus of Med17 and the Middle module were later shown to be highly accurate by high-resolution EM [75] and crystallography [76]. For other features, our model remains the best description of Mediator available. These include the Tail module which is highly dynamic, contains many disordered regions, and participates in "fuzzy" interactions with activators such as Gcn4 (see Section 3.6). Thus, a high-resolution EM reconstruction of the Tail module might not be possible with existing technology, but ISD approaches can infer molecular features at pseudo-residue level precision. For instance, our Tail module model predicts specific structural features such as the N-terminal portion of Med5 extending across the β-propeller domain of Med16 to contact Med15 [17] and associates a precision (RMSF value) with these predictions. Much of the Tail module remains resolved at poor resolution in our model due to lack of adequate cross-linking restraints. Hence, our structural understanding of the Tail module would likely be improved by the application of cross-linking methodologies that increased coverage in these domains.

In addition to the high frequency fluctuations of the Tail module, our data suggest that it also undergoes a more coordinated rearrangement upon PIC binding that places Med5 in contact with Med1 of the Middle module (Section 3.6). This represents a more stable conformational heterogeneity that is resolvable by EM. CLMS can aid in understanding these states through the use of quantitative MS methods. For instance, the Med-PIC can be cross-linked in the presence and absence of Gcn4 using reagents with different isotopic signatures [138]. The cross-linked samples are combined, prepared for mass spectrometry analysis, where the isotopic ratio of cross-linked peptide precursors reports on states that are enriched or depleted in the presence of Gcn4. Quantitative CLMS (qCLMS) can be extended to multiple states along a mechanistic pathway, such as promoter escape, through using multiplexed proteomics methods such as PRM [139] or TMT labelling [140]. Using such methods, sub-conformations of a protein assembly do not need to be homogenous, as long as one state is enriched sufficiently to be quantified by MS. Hence, these sorts of qCLMS approaches as inputs to ISD might be useful in resolving more dynamic and heterogenous populations of protein assemblies as well, although they do require that the individual mechanistic states can be partially stabilized through modulating the sample conditions, without significant changes to the protein composition.

However, the true value of integrative models going forward lies not in traditional structural biology, which focuses on elucidating stable elements of macromolecular tertiary and quaternary structure, but as tools to study the non-classical, sequence-independent and low-affinity interactions known as "fuzzy-interactions" [141,142] and protein quinary structure [143–145]. These emerging areas of study have been implicated in transcriptional regulation through the binding of gene specific activators to their co-activator targets in complexes without a single ground-state conformation, and by governing the formation and organization of cellular structures that sequester large ensembles of transcriptional machinery in sub-nuclear, membrane-less compartments. These processes are inherently conformationally dynamic and heterogenous and thus not easily represented by single-state atomistic molecular models that can be resolved by high-resolution EM or X-ray derived techniques. Instead, we anticipate that multi-state representations of these systems will be accessible via extension of existing integrative modeling methodology. Experimentally, techniques that are well integrated into integrative methods such as CLMS and mixed resolution EM will continue to play an important role, as will other methods such as NMR [146,147], real-time immunofluorescence microscopy [148], proximity labelling-MS [149], and chemical footprinting-MS methods [150] such as hydrogen-deuterium exchange or oxidative footprinting.

Quinary structure represents a fifth level of protein organization

based on transient and weak-affinity protein interactions [143]. These fleeting interactions are difficult to capture experimentally, yet they play a major role in organizing the intracellular environment at a scale larger than a protein complex yet smaller than an organelle. For instance, various multi-step metabolic pathways rely on the spatial localization of multiple macromolecular complexes for efficient biosynthesis. Membrane-less compartments have been shown to form in the nucleus to sequester molecules with related functions [151]. Examples include: Nuclear speckles, which contain mRNA splicing factors, the Nucleolus, which contains ribosome synthesis machinery, and super-enhancers, which bring together extremely high concentrations of transcriptional machines including activators, co-activators and the PIC, at promoters that seem to have particularly important roles in regulating cellular differentiation and oncogenesis [40]. In lieu of phospholipid membranes, these microenvironments are defined by the oligomerization of multi-valent scaffold proteins that show a high prevalence of intrinsically disordered domains with low sequence complexity [152]. For instance, heterochromatin is dynamically sequestered by the multivalent heterochromatin protein HP1α, which contains two structurally disordered regions. Phosphorylation at one of these domains regulates extension at the other to a conformation that favours inter-HP1α interactions that act as a scaffold to condense and isolate transcriptionally silent regions of the genome [153]. Membrane-less compartments exhibit properties of liquid droplets, such as phase separation and droplet fission and fusion [143]. While the free energy of these interactions is weak (1 kcal/mol, an order of magnitude lower than quaternary protein structure), the densely crowded cellular environment (200–400 g/L protein concentration) amplifies the chemical forces exerted by neighbouring molecules through electrostatic, hydrophobic, polar, and hydrogen bonding interactions. Hence, small chemical modifications can lead to large-scale re-organization of the cellular microenvironment.

In addition to contributing to the characterization of membrane-less compartments, ISD will also likely play an important role in understanding transcription complexes at higher levels of organization. In addition to the Mediator-PIC assembly, numerous protein complexes act at a gene promoter to influence its transcription. These include other co-activators such as TFIID and SAGA as well as chromatin re-modeling complexes such as SWI/SNF and RSC, and inhibitory modules such as the Mediator Kinase module. How these complexes coordinate their actions with the Med-PIC, the nature of their super-assemblies if any, and how promoter identity influences mechanism are open questions. These larger assemblies represent a logical next step in continuing our modular approach to ISD which builds models of increasingly larger assemblies based in part on transferability of data acquired on smaller subassemblies. However, unlike previous iterations of this process, these larger assemblies are governed by quinary interactions. Super enhancers are clusters of spatially proximate enhancers that regulate a single promoter and contain elevated concentrations of Mediator, the general transcriptional machinery and other co-activators relative to regular promoters [40]. Various groups have proposed that super enhancers represent phase separated droplets which sequester high concentrations of transcriptional machinery at gene promoters of particular importance to oncogenesis and determination of cellular fate [41,42,154]. Mediator is enriched in super enhancers and antibodies directed at Middle module subunit Med1, which contains a high portion of low complexity domains, are used to characterize super enhancers in ChIP-Seq or fluorescence imaging experiments. The organization and composition of super-enhancers is very much unknown. For instance, what is the nature of the scaffold that creates a phase-separated compartment and what are the client proteins? The stoichiometry of these compartments is also unknown and likely to vary, and it is likely that organization of the individual protein machines within the compartment is governed by weak interactions with multiple modes of binding. Hence, it is unlikely that EM will ever be able to produce a single, atomistic ground state structure of a super-enhancer. However, ISD
approaches which describe biological systems probabilistically and can address multi-state systems, can potentially be developed to model super enhancers and other conglomerations of transcription complexes in a way that yield clues as to the structural principles that influence gene regulation. Since these complexes are good examples of quinary structures that will be very labile in nature, integrative structural approaches will likely play a key role in advancing the structural biology of these complexes, as it did in the emerging picture of the Mediator complex and related transcription complexes.

The value of integrative models in the future is in addressing problems in structural biology that are not tractable by high-resolution methods. These are problems involving heterogeneous structures, highly dynamic conformational fluctuations, intrinsically disordered domains, “fuzzy”-interactions, and higher-order structures within a cell encompassing ensembles of transcription related protein complexes. These areas comprise many emerging themes in transcriptional regulation, and they challenge existing structural biology techniques and paradigms which are best suited for stable, homogenous assemblies.

## Acknowledgments

MJT is supported by funds from the NIH (R01 GM121962, R01 GM044037) and the Dr Miriam & Sheldon G. Adelson Medical Research Foundation. PJR is a supported with a Medical Research Council Career Development Award (MR/R008795/1).

## References

[1] J.C. Kendrew, R.E. Dickerson, B.E. Strandberg, R.G. Hart, D.R. Davies, D.C. Phillips, V.C. Shore, Structure of myoglobin: A three-dimensional Fourier synthesis at 2 Å resolution, Nature 185 (1960) 422–427.

[2] C.C. Blake, D.F. Koenig, G.A. Mair, A.C. North, D.C. Phillips, V.R. Sarma, Structure of hen egg-white lysozyme. A three-dimensional Fourier synthesis at 2 Ångstrom resolution, Nature 206 (1965) 757–761.

[3] M.F. Perutz, H. Miurhead, J.M. Cox, L.C. Goaman, F.S. Mathews, E.L. McGandy, L.E. Webb, Three-dimensional Fourier synthesis of horse oxyhaemoglobin at 2.8 Å resolution: (1) x-ray analysis, Nature 219 (1968) 29–32.

[4] J. Deisenhofer, O. Epp, K. Miki, R. Huber, H. Michel, Structure of the protein subunits in the photosynthetic reaction centre of Rhodopseudomonas viridis at 3Å resolution, Nature 318 (1986) 618–624.

[5] K. Luger, A.W. Mader, R.K. Richmond, D.F. Sargent, T.J. Richmond, Crystal structure of the nucleosome core particle at 2.8 Å resolution, Nature 389 (1997) 251–260, https://doi.org/10.1038/38444.

[6] P. Cramer, D.A. Bushnell, R.D. Kornberg, Structural basis of transcription: RNA polymerase II at 2.8 Ångstrom resolution, Science 292 (2001) 1863–1876, https://doi.org/10.1126/science.1059493.

[7] M. Selmer, C.M. Dunham, F.V. th Murphy, A. Weixlbaumer, S. Petry, A.C. Kelley, J.R. Weir, V. Ramakrishnan, Structure of the 70S ribosome complexed with mRNA and tRNA, Science 313 (2006) 1935–1942, https://doi.org/10.1126/science.1131127.

[8] T.W. Jeng, R.A. Crowther, G. Stubbs, W. Chiu, Visualization of alpha-helices in tobacco mosaic virus by cryo-electron microscopy, J. Mol. Biol. 205 (1989) 251–257.

[9] R.A. Crowther, L.A. Amos, J.T. Finch, D.J. De Rosier, A. Klug, Three dimensional reconstructions of spherical viruses by Fourier synthesis from electron micrographs, Nature 226 (1970) 421–425.

[10] A. Merk, A. Bartesaghi, S. Banerjee, V. Falconieri, P. Rao, M.I. Davis, R. Pragani, M.B. Boxer, L.A. Earl, J.L.S. Milne, S. Subramanian, Breaking cryo-EM resolution barriers to facilitate drug discovery, Cell 165 (2016) 1698–1707, https://doi.org/10.1016/j.cell.2016.05.040.

[11] A.E. Leschziner, E. Nogales, Visualizing flexibility at molecular resolution: analysis of heterogeneity in single-particle electron microscopy reconstructions, Annu. Rev. Biophys. Biomol. Struct. 36 (2007) 43–62, https://doi.org/10.1146/annurev.biophys.36.040306.132742.

[12] C.M. Spahn, P.A. Penczek, Exploring conformational modes of macromolecular assemblies by multiparticle cryo-EM, Curr. Opin. Struct. Biol. 19 (2009) 623–631, https://doi.org/10.1016/j.sbi.2009.08.001.

[13] H.E. White, A. Ignatiou, D.K. Clare, E.V. Orlova, Structural study of heterogeneous biological samples by cryoelectron microscopy and image processing, Biomed. Res. Int. 2017 (2017) 1032432, https://doi.org/10.1155/2017/1032432.

[14] D. Schneidman-Duhovny, R. Pellarin, A. Sali, Uncertainty in integrative structural modeling, Curr. Opin. Struct. Biol. 28 (2014) 96–104, https://doi.org/10.1016/j.sbi.2014.08.001.

[15] B. Webb, S. Viswanath, M. Bonomi, R. Pellarin, C.H. Greenberg, D. Saltzberg, A. Sali, Integrative structure modeling with the integrative modeling platform, Protein Sci. 27 (2018) 245–258, https://doi.org/10.1002/pro.3311.

[16] J. Luo, P. Cimermancic, S. Viswanath, C.C. Ebmeier, B. Kim, M. Dehecq, V. Raman, C.H. Greenberg, R. Pellarin, A. Sali, D.J. Taatjes, S. Hahn, J. Ranish, Architecture of the human and yeast general transcription and DNA repair factor TFIIH, Mol. Cell. 59 (2015) 794–806, https://doi.org/10.1016/j.molcel.2015.07.016.

[17] P.J. Robinson, M.J. Trnka, R. Pellarin, C.H. Greenberg, D.A. Bushnell, R. Davis, A.L. Burlingame, A. Sali, R.D. Kornberg, Molecular architecture of the yeast Mediator complex, Elife 4 (2015), https://doi.org/10.7554/eLife.08719.

[18] P.J. Robinson, M.J. Trnka, D.A. Bushnell, R.E. Davis, P.J. Mattei, A.L. Burlingame, R.D. Kornberg, Structure of a complete mediator-RNA polymerase II pre-initiation complex, Cell 166 (2016) 1411–1422, https://doi.org/10.1016/j.cell.2016.08.050 e16.

[19] P.C. Whitford, A. Ahmed, Y. Yu, S.P. Hennelly, F. Tama, C.M. Spahn, J.N. Onuchic, K.Y. Sanbonmatsu, Excited states of ribosome translocation revealed through integrative molecular modeling, Proc. Natl. Acad. Sci. U.S.A. 108 (2011) 18943–18948, https://doi.org/10.1073/pnas.1108363108.

[20] J. Fernandez-Martinez, S.J. Kim, Y. Shi, P. Upla, R. Pellarin, M. Gagnon, I.E. Chemmama, J. Wang, I. Nudelman, W. Zhang, R. Williams, W.J. Rice, D.L. Stokes, D. Zenklusen, B.T. Chait, A. Sali, M.P. Rout, Structure and function of the nuclear pore complex cytoplasmic mRNA export platform, Cell 167 (2016) 1215–1228, https://doi.org/10.1016/j.cell.2016.10.028 e25.

[21] S.J. Kim, J. Fernandez-Martinez, I. Nudelman, Y. Shi, W. Zhang, B. Raveh, T. Herricks, B.D. Slaughter, J.A. Hogan, P. Upla, I.E. Chemmama, R. Pellarin, I. Echeverria, M. Shivaraju, A.S. Chaudhury, J. Wang, R. Williams, J.R. Unruh, C.H. Greenberg, E.Y. Jacobs, Z. Yu, M.J. de la Cruz, R. Mironska, D.L. Stokes, J.D. Aitchison, M.F. Jarrold, J.L. Gerton, S.J. Ludtke, C.W. Akey, B.T. Chait, A. Sali, M.P. Rout, Integrative structure and functional anatomy of a nuclear pore complex, Nature 555 (2018) 475–482, https://doi.org/10.1038/nature26003.

[22] E. Obayashi, R.E. Luna, T. Nagata, P. Martin-Marcos, H. Hiraishi, C.R. Singh, J.P. Erzberger, F. Zhang, H. Arthanari, J. Morris, R. Pellarin, C. Moore, I. Harmon, E. Papadopoulos, H. Yoshida, M.L. Nasr, S. Unzai, B. Thompson, E. Aube, S. Hustak, F. Stengel, E. Dagraca, A. Ananbandam, P. Gao, T. Urano, A.G. Hinnebusch, G. Wagner, K. Asano, Molecular landscape of the ribosome pre-initiation complex during mRNA scanning: structural role for eIF3c and its control by eIF5, Cell. Rep. 18 (2017) 2651–2663, https://doi.org/10.1016/j.celrep.2017.02.052.

[23] Y. Shi, J. Fernandez-Martinez, E. Tjioe, R. Pellarin, S.J. Kim, R. Williams, D. Schneidman-Duhovny, A. Sali, M.P. Rout, B.T. Chait, Structural characterization by cross-linking reveals the detailed architecture of a coatomer-related heptameric module from the nuclear pore complex, Mol. Cell. Proteomics MCP 13 (2014) 2927–2943, https://doi.org/10.1074/mcp.M114.041673.

[24] Y. Shi, R. Pellarin, P.C. Fridy, J. Fernandez-Martinez, M.K. Thompson, Y. Li, Q.J. Wang, A. Sali, M.P. Rout, B.T. Chait, A strategy for dissecting the architectures of native macromolecular assemblies, Nat. Methods 12 (2015) 1135–1138, https://doi.org/10.1038/nmeth.3617.

[25] F.C. Holstege, E.G. Jennings, J.J. Wyrick, T.I. Lee, C.J. Hengartner, M.R. Green, T.R. Golub, E.S. Lander, R.A. Young, Dissecting the regulatory circuitry of a eukaryotic genome, Cell 95 (1998) 717–728.

[26] Y. Takagi, R.D. Kornberg, Mediator as a general transcription factor, J. Biol. Chem. 281 (2006) 80–89, https://doi.org/10.1074/jbc.M508253200.

[27] Y.J. Kim, S. Bjorklund, Y. Li, M.H. Sayre, R.D. Kornberg, A multiprotein mediator of transcriptional activation and its interaction with the C-terminal repeat domain of RNA polymerase II, Cell 77 (1994) 599–608.

[28] R.D. Kornberg, Mediator and the mechanism of transcriptional activation, Trends Biochem. Sci. 30 (2005) 235–239, https://doi.org/10.1016/j.tibs.2005.03.011.

[29] L.A. Allison, C.J. Ingles, Mutations in RNA polymerase II enhance or suppress mutations in GAL4, Proc. Natl. Acad. Sci. U.S.A. 86 (1989) 2794–2798.

[30] C. Scafe, D. Chao, J. Lopes, J.P. Hirsch, S. Henry, R.A. Young, RNA polymerase II C-terminal repeat influences response to transcriptional enhancer signals, Nature 347 (1990) 491–494, https://doi.org/10.1038/347491a0.

[31] S.M. Liao, I.C. Taylor, R.E. Kingston, R.A. Young, RNA polymerase II carboxy-terminal domain contributes to the response to multiple acidic activators in vitro, Genes Dev. 5 (1991) 2431–2440.

[32] J. Soutourina, S. Wydau, Y. Ambroise, C. Boschiero, M. Werner, Direct interaction of RNA polymerase II and mediator required for transcription in vivo, Science 331 (2011) 1451–1454, https://doi.org/10.1126/science.1200188.

[33] E. Paul, Z.I. Zhu, D. Landsman, R.H. Morse, Genome-wide association of mediator and RNA polymerase II in wild-type and mediator mutant yeast, Mol. Cell. Biol. 35 (2015) 331–342, https://doi.org/10.1128/MCB.00991-14.

[34] P.J. Robinson, D.A. Bushnell, M.J. Trnka, A.L. Burlingame, R.D. Kornberg, Structure of the mediator head module bound to the carboxy-terminal domain of RNA polymerase II, Proc. Natl. Acad. Sci. U. S. A. 109 (2012) 17931–17935, https://doi.org/10.1073/pnas.1215241109.

[35] P.J. Laybourn, M.E. Dahmus, Phosphorylation of RNA polymerase IIA occurs subsequent to interaction with the promoter and before the initiation of transcription, J. Biol. Chem. 265 (1990) 13165–13173.

[36] J.Q. Svejstrup, Y. Li, J. Fellows, A. Gnatt, S. Bjorklund, R.D. Kornberg, Evidence for a mediator cycle at the initiation of transcription, Proc. Natl. Acad. Sci. U.S.A. 94 (1997) 6075–6078.

[37] T. Max, M. Sogaard, J.Q. Svejstrup, Hyperphosphorylation of the C-terminal repeat domain of RNA polymerase II facilitates dissociation of its complex with mediator, J. Biol. Chem. 282 (2007) 14113–14120, https://doi.org/10.1074/jbc.M701345200.

[38] C. Jeronimo, M.F. Langelier, A.R. Bataille, J.M. Pascal, B.F. Pugh, F. Robert, Tail and kinase modules differently regulate core mediator recruitment and function in vivo, Mol. Cell. 64 (2016) 455–466, https://doi.org/10.1016/j.molcel.2016.09.002.

[39] N. Petrenko, Y. Jin, K.H. Wong, K. Struhl, Mediator undergoes a compositional
change during transcriptional activation, Mol. Cell. 64 (2016) 443–454, https://doi.org/10.1016/j.molcel.2016.09.015.

[40] W.A. Whyte, D.A. Orlando, D. Hnisz, B.J. Abraham, C.Y. Lin, M.H. Kagey, P.B. Rahl, T.I. Lee, R.A. Young, Master transcription factors and mediator establish super-enhancers at key cell identity genes, Cell 153 (2013) 307–319, https://doi.org/10.1016/j.cell.2013.03.035.

[41] D. Hnisz, K. Shrinivas, R.A. Young, A.K. Chakraborty, P.A. Sharp, A phase separation model for transcriptional control, Cell 169 (2017) 13–23, https://doi.org/10.1016/j.cell.2017.02.007.

[42] W.K. Cho, J.H. Spille, M. Hecht, C. Lee, C. Li, V. Grube, I. Cisse, Mediator and RNA polymerase II clusters associate in transcription-dependent condensates, Science 361 (2018) 412–415, https://doi.org/10.1126/science.aar4199.

[43] P.S. Brzovic, C.C. Heikaus, L. Kisseylev, R. Vernon, E. Herbig, D. Pacheco, L. Warfield, P. Littlefield, D. Baker, R.E. Klevit, S. Hahn, The acidic transcription activator Gcn4 binds the mediator subunit Gal11/Med15 using a simple protein interface forming a fuzzy complex, Mol. Cell. 44 (2011) 942–953, https://doi.org/10.1016/j.molcel.2011.11.008.

[44] L.M. Tuttle, D. Pacheco, L. Warfield, J. Luo, J. Ranish, S. Hahn, R.E. Klevit, Gcn4-mediator specificity is mediated by a large and dynamic fuzzy protein-protein complex, Cell. Rep. 22 (2018) 3251–3264, https://doi.org/10.1016/j.celrep.2018.02.097.

[45] A. Boijja, I.A. Klein, B.R. Sabari, A. Dall'Agnese, E.L. Coffey, A.V. Zamudio, C.H. Li, K. Shrinivas, J.C. Manteiga, N.M. Hannett, B.J. Abraham, L.K. Afeyan, Y.E. Guo, J.K. Rimel, C.B. Fant, J. Schuijers, T.I. Lee, D.J. Taatjes, R.A. Young, Transcription factors activate genes through the phase-separation capacity of their activation domains, Cell 175 (2018) 1842–1855, https://doi.org/10.1016/j.cell.2018.10.042 e16.

[46] D. Russel, K. Lasker, J. Webb, J. Velazquez-Muriel, E. Tjoei, D. Schneidman-Duhovny, B. Peterson, A. Sali, Putting the pieces together: integrative modeling platform software for structure determination of macromolecular assemblies, e1001244, PLoS Biol. 10 (2012), https://doi.org/10.1371/journal.pbio.1001244.

[47] M. Habeck, M. Nilges, W. Rieping, Replica-exchange Monte Carlo scheme for bayesian data analysis, 018105, Phys. Rev. Lett. 94 (2005), https://doi.org/10.1103/PhysRevLett.94.018105.

[48] F.J. Asturias, Y.W. Jiang, L.C. Myers, C.M. Gustafsson, R.D. Kornberg, Conserved structures of mediator and RNA polymerase II holoenzyme, Science 283 (1999) 985–987.

[49] B. Guglielmi, N.L. van Berkum, B. Klapholz, T. Bijma, M. Boube, C. Boschiro, H.M. Bourbon, F.C. Holstege, M. Werner, A high resolution protein interaction map of the yeast Mediator complex, Nucleic Acids Res. 32 (2004) 5379–5391, https://doi.org/10.1093/nar/gkh878.

[50] C.M. Thompson, A.J. Koleske, D.M. Chao, R.A. Young, A multisubunit complex associated with the RNA polymerase II CTD and TATA-binding protein in yeast, Cell 73 (1993) 1361–1375.

[51] L.C. Myers, C.M. Gustafsson, D.A. Medhnell, M. Lui, H. Erdument-Bromage, P. Tempst, R.D. Kornberg, The Med proteins of yeast and their function through the RNA polymerase II carboxy-terminal domain, Genes Dev. 12 (1998) 45–54.

[52] Y.C. Lee, Y.J. Kim, Requirement for a functional interaction between mediator components Med6 and Srb4 in RNA polymerase II transcription, Mol. Cell. Biol. 18 (1998) 5364–5370.

[53] T.I. Lee, J.J. Wyrick, S.S. Koh, E.G. Jennings, E.L. Gadbois, R.A. Young, Interplay of positive and negative regulators in transcription initiation by RNA polymerase II holoenzyme, Mol. Cell. Biol. 18 (1998) 4455–4462.

[54] J.S. Kang, S.H. Kim, M.S. Hwang, S.J. Han, Y.C. Lee, Y.J. Kim, The structural and functional organization of the yeast Mediator complex, J. Biol. Chem. 276 (2001) 42003–42010, https://doi.org/10.1074/jbc.M105961200.

[55] M.L. Nonet, R.A. Young, Intragenic and extragenetic suppressors of mutations in the heptapeptide repeat domain of Saccharomyces cerevisiae RNA polymerase II, Genetics 123 (1989) 715–724.

[56] C.M. Thompson, R.A. Young, General requirement for RNA polymerase II holoenzymes in vivo, Proc. Natl. Acad. Sci. U.S.A. 92 (1995) 4587–4590.

[57] L. Lariviere, S. Geiger, S. Hoeppner, S. Rother, K. Strasser, P. Cramer, Structure and TBP binding of the Mediator head subcomplex Med8-Med18-Med20, Nat. Struct. Mol. Biol. 13 (2006) 895–901, https://doi.org/10.1038/nsmb1143.

[58] T. Imasaki, G. Galero, G. Cai, K.L. Tsai, K. Yamada, F. Cardelli, H. Erdument-Bromage, P. Tempst, I. Berger, G.L. Kornberg, F.J. Asturias, R.D. Kornberg, Y. Takagi, Architecture of the Mediator head module, Nature 475 (2011) 240–243, https://doi.org/10.1038/nature10162.

[59] G. Cai, T. Imasaki, K. Yamada, F. Cardelli, Y. Takagi, F.J. Asturias, Mediator head module structure and functional interactions, Nat. Struct. Mol. Biol. 17 (2010) 273–279, https://doi.org/10.1038/nsmb.1757.

[60] M. Seizl, L. Lariviere, T. Pfaffeneder, L. Wenzeeck, P. Cramer, Mediator head subcomplex Med11/22 contains a common helix bundle building block with a specific function in transcription initiation complex stabilization, Nucleic Acids Res. 39 (2011) 6291–6304, https://doi.org/10.1093/nar/gkr229.

[61] S. Baumli, S. Hoeppner, P. Cramer, A conserved mediator hinge revealed in the structure of the MED7-MED21 (Med7.Srb7) heterodimer, J. Biol. Chem. 280 (2005) 18171–18178, https://doi.org/10.1074/jbc.M413466200.

[62] T. Koschubs, M. Seizl, L. Lariviere, F. Kurth, S. Baumli, D.E. Martin, P. Cramer, Identification, structure, and functional requirement of the Mediator subcomplex Med7N/31, EMBO J. 28 (2009) 69–80, https://doi.org/10.1038/emboj.2008.254.

[63] J.K. Thakur, H. Arthanari, F. Yang, S.J. Pan, X. Fan, J. Breger, D.P. Frueh, K. Gulshan, D.K. Li, E. Mylonakis, K. Struhl, W.S. Moye-Rowley, B.P. Cormack, G. Wagner, A.M. Naar, A nuclear receptor-like pathway regulating multidrug resistance in fungi, Nature 452 (2008) 604–609, https://doi.org/10.1038/nature06836.

[64] P. Uetz, L. Giot, G. Cagney, T.A. Mansfield, R.S. Judson, J.R. Knight, D. Lockshon, V. Narayan, M. Srinivasan, P. Pochart, A. Qureshi-Emili, Y. Li, B. Godwin, D. Conover, T. Kalbfleisch, G. Vijayadamodar, M. Yang, M. Johnston, S. Fields, J.M. Rothberg, A comprehensive analysis of protein-protein interactions in Saccharomyces cerevisiae, Nature 403 (2000) 623–627, https://doi.org/10.1038/35001009.

[65] T. Ito, T. Chiba, R. Ozawa, M. Yoshida, M. Hattori, Y. Sakaki, A comprehensive two-hybrid analysis to explore the yeast protein interactome, Proc Natl Acad Sci U. A. 98 (2001) 4569–4574, https://doi.org/10.1073/pnas.061034498.

[66] S.S. Koh, A.Z. Ansari, M. Ptashne, R.A. Young, An activator target in the RNA polymerase II holoenzyme, Mol. Cell. 1 (1998) 895–904.

[67] T. Koschubs, K. Lorenzen, S. Baumli, S. Sandstrom, A.J. Heck, P. Cramer, Preparation and topology of the Mediator middle module, Nucleic Acids Res. 38 (2010) 3186–3195, https://doi.org/10.1093/nar/gkq029.

[68] F. Zhang, L. Sumibcay, A.G. Hinnebusch, M.J. Swanson, A triad of subunits from the Gal11/tail domain of Srb mediator is an in vivo target of transcriptional activator Gcn4p, Mol. Cell. Biol. 24 (2004) 6871–6886, https://doi.org/10.1128/MCB.24.15.6871-6886.2004.

[69] J. Beve, G.Z. Hu, L.C. Myers, D. Balciunas, O. Werngren, K. Hultenby, R. Wibom, H. Ronne, C.M. Gustafsson, The structural and functional role of Med5 in the yeast Mediator tail module, J. Biol. Chem. 280 (2005) 41366–41372, https://doi.org/10.1074/jbc.M511181200.

[70] K.L. Tsai, C. Tomomori-Sato, S. Sato, R.C. Conaway, J.W. Conaway, F.J. Asturias, Subunit architecture and functional modular rearrangements of the transcriptional mediator complex, Cell 158 (2014) 463, https://doi.org/10.1016/j.cell.2014.06.036.

[71] X. Wang, Q. Sun, Z. Ding, J. Ji, J. Wang, X. Kong, J. Yang, G. Cai, Redefining the modular organization of the core Mediator complex, Cell. Rep. 24 (2014) 796–808, https://doi.org/10.1038/cr.2014.64.

[72] A. Sakai, Y. Shimizu, S. Kondou, T. Chibazakura, F. Hishinuma, Structure and molecular analysis of RGR1, a gene required for glucose repression of Saccharomyces cerevisiae, Mol. Cell. Biol. 10 (1990) 4130–4138.

[73] Y. Li, S. Bjorklund, Y.W. Jiang, Y.J. Kim, W.S. Lane, D.J. Stilman, R.D. Kornberg, Yeast global transcriptional regulators Sin4 and Rgr1 are components of mediator complex/RNA polymerase II holoenzyme, Proc. Natl. Acad. Sci. U.S.A. 92 (1995) 10864–10868.

[74] M.A. Cevher, Y. Shi, D. Li, B.T. Chait, S. Malik, R.G. Roeder, Reconstitution of active human core Mediator complex reveals a critical role of the MED14 subunit, Nat. Struct. Mol. Biol. 21 (2014) 1028–1034, https://doi.org/10.1038/nsmb.2914.

[75] K.L. Tsai, X. Yu, S. Gopalan, T.C. Chao, Y. Zhang, L. Florens, M.P. Washburn, K. Murakami, R.C. Conaway, J.W. Conaway, F.J. Asturias, Mediator structure and rearrangements required for holoenzyme formation, Nature 544 (2017) 196–201, https://doi.org/10.1038/nature21393.

[76] K. Nozawa, T.R. Schneider, P. Cramer, Core Mediator structure at 3.4 Å model of transcription initiation complex, Nature 545 (2017) 248–251, https://doi.org/10.1038/nature22328.

[77] S. Schilbach, M. Hantsche, D. Tegunov, C. Dienemann, C. Wigge, H. Urlaub, P. Cramer, Structures of transcription pre-initiation complex with TFIIH and Mediator, Nature 551 (2017) 204–209, https://doi.org/10.1038/nature24282.

[78] L. Lariviere, C. Plaschka, M. Seizl, L. Wenzeeck, F. Kurth, P. Cramer, Structure of the Mediator head module, Nature 492 (2012) 448–451, https://doi.org/10.1038/nature11670.

[79] L. Lariviere, C. Plaschka, M. Seizl, E.V. Petrotchenko, L. Wenzeeck, C.H. Borchers, P. Cramer, Model of the Mediator middle module based on protein cross-linking, Nucleic Acids Res. 41 (2013) 9266–9273, https://doi.org/10.1093/nar/gkt704.

[80] E. Guzman, J.T. Lis, Transcription factor TFIIH is required for promoter melting in vivo, Mol. Cell. Biol. 19 (1999) 5652–5658.

[81] G.S. Winkler, S.J. Araujo, U. Fiedler, W. Vermeulen, F. Coin, J.M. Egly, J.H. Hoeijmakers, R.D. Wood, H.T. Timmers, G. Weeda, TFIIH with inactive XPD helicase functions in transcription initiation but is defective in DNA repair, J. Biol. Chem. 275 (2000) 4258–4266.

[82] W.J. Feaver, O. Gileadi, Y. Li, R.D. Kornberg, CTD kinase associated with yeast RNA polymerase II initiation factor b, Cell 67 (1991) 1223–1230.

[83] Y. He, J. Fang, D.J. Taatjes, E. Nogales, Structural visualization of key steps in human transcription initiation, Nature 495 (2013) 481–486, https://doi.org/10.1038/nature11991.

[84] K. Murakami, K.L. Tsai, N. Kalisman, D.A. Bushnell, F.J. Asturias, R.D. Kornberg, Structure of an RNA polymerase II preinitiation complex, Proc. Natl. Acad. Sci. U.S.A. 112 (2015) 13543–13548, https://doi.org/10.1073/pnas.1518255112.

[85] B.J. Gibbons, E.J. Brignole, M. Azubel, K. Murakami, N.R. Voss, D.A. Bushnell, F.J. Asturias, R.D. Kornberg, Subunit architecture of general transcription factor TFIIH, Proc. Natl. Acad. Sci. U.S.A. 109 (2012) 1949–1954, https://doi.org/10.1073/pnas.1105266109.

[86] H.J. Himmelfarb, J. Pearlberg, D.H. Last, M. Ptashne, GAL11P: a yeast mutation that potentiates the effect of weak GAL4-derived activators, Cell 63 (1990) 1299–1309.

[87] M. Tanaka, Modulation of promoter occupancy by cooperative DNA binding and activation-domain function is a major determinant of transcriptional regulation by activators-in vivo, Proc. Natl. Acad. Sci. U.S.A. 93 (1996) 4311–4315.

[88] C. Bernecky, D.J. Taatjes, Activator-mediator binding stabilizes RNA polymerase II orientation within the human mediator-RNA polymerase II-TFIIH assembly, J. Mol. Biol. 417 (2012) 387–394, https://doi.org/10.1016/j.jmb.2012.02.014.

[89] E. Vojnic, A. Mourao, M. Seizl, B. Simon, L. Wenzeeck, L. Lariviere, S. Baumli, K. Baumgart, M. Meisterernst, M. Sattler, P. Cramer, Structure and VP16 binding of the Mediator Med25 activator interaction domain, Nat. Struct. Mol. Biol. 18
(2011) 404–409, https://doi.org/10.1038/nsmb.1997.

[90] E. Herbig, L. Warfield, L. Fish, J. Fishburn, B.A. Knutson, B. Moorefield, D. Pacheco, S. Hahn, Mechanism of Mediator recruitment by tandem Gcn4 activation domains and three Gal11 activator-binding domains, Mol. Cell. Biol. 30 (2010) 2376–2390, https://doi.org/10.1128/MCB.010469-09.

[91] I. Jedidi, F. Zhang, H. Qiu, S.J. Stahl, I. Palmer, J.D. Friedman, P.S. Nadaud, S. Mukherjee, P.T. Wingfield, C.P. Jaroniec, A.G. Hinnebusch, Activator Gcn4 employs multiple segments of Med15/Gal11, including the KIX domain, to recruit mediator to target genes in vivo, J. Biol. Chem. 285 (2010) 2438–2455, https://doi.org/10.1074/jbc.M109.071589.

[92] A. Leitner, T. Walzthoeni, A. Kahraman, F. Herzog, O. Rinner, M. Beck, R. Aebersold, Probing native protein structures by chemical cross-linking, mass spectrometry, and bioinformatics, Mol. Cell. Proteomics 9 (2010) 1634–1649, https://doi.org/10.1074/mcp.R000001-MCP201.

[93] C. Yu, L. Huang, Cross-linking mass spectrometry: an emerging technology for interactomics and structural biology, Anal. Chem. 90 (2018) 144–165, https://doi.org/10.1021/acs.analchem.7b04431.

[94] F. Chu, D.T. Thornton, H.T. Nguyen, Chemical cross-linking in the structural analysis of protein assemblies, Methods 144 (2018) 53–63, https://doi.org/10.1016/j.ymeth.2018.05.023.

[95] O. Rinner, J. Seebacher, T. Walzthoeni, L.N. Mueller, M. Beck, A. Schmidt, M. Mueller, R. Aebersold, Identification of cross-linked peptides from large sequence databases, Nat. Methods 5 (2008) 315–318, https://doi.org/10.1038/nmeth.1192.

[96] Expanding the chemical cross-linking toolbox by the use of multiple proteases and enrichment by size exclusion chromatography | Molecular & Cellular Proteomics, (n.d.). http://www.mcponline.org/content/early/2012/01/27/mcp.M111.014126short (accessed January 6, 2019).

[97] R.H. Perry, R.G. Cooks, R.J. Noll, Orbitrap mass spectrometry: Instrumentation, ion motion and applications, Mass Spectrom. Rev. 27 (2008) 661–699, https://doi.org/10.1002/mas.20186.

[98] T. Walzthoeni, M. Claassen, A. Leitner, F. Herzog, S. Bohn, F. Förster, M. Beck, R. Aebersold, False discovery rate estimation for cross-linked peptides identified by mass spectrometry, Nat. Methods 9 (2012) 901–903, https://doi.org/10.1038/nmeth.2103.

[99] B. Yang, Y.-J. Wu, M. Zhu, S.-B. Fan, J. Lin, K. Zhang, S. Li, H. Chi, Y.-X. Li, H.-F. Chen, S.-K. Luo, Y.-H. Ding, L.-H. Wang, Z. Hao, L.-Y. Xiu, S. Chen, K. Ye, S.-M. He, M.-Q. Dong, Identification of cross-linked peptides from complex samples, Nat. Methods 9 (2012) 904–906, https://doi.org/10.1038/nmeth.2099.

[100] M.J. Trnka, P.R. Baker, P.J.J. Robinson, A.L. Burlingame, R.J. Chalkley, Matching cross-linked peptide spectra: only as good as the worse identification, Mol. Cell. Proteomics 13 (2014) 420–434, https://doi.org/10.1074/mcp.M113.034009.

[101] J.P. Erzberger, F. Stengel, R. Pellarin, S. Zhang, T. Schaefer, C.H.S. Aylert, P. Cimermančič, D. Boehringer, A. Sali, R. Aebersold, N. Ban, Molecular architecture of the 40S-eIF1-eIF3 translation initiation complex, Cell 158 (2014) 1123–1135, https://doi.org/10.1016/j.cell.2014.07.044.

[102] X. Tang, J.E. Bruce, A new cross-linking strategy: protein interaction reporter (PIR) technology for protein – protein interaction studies, Mol. Biosyst. 6 (2010) 939–947, https://doi.org/10.1039/B920876C.

[103] A. Kao, C. Chiu, D. Vellucci, Y. Yang, V.R. Patel, S. Guan, A. Randall, P. Baldi, S.D. Rychnovsky, L. Huang, Development of a novel cross-linking strategy for fast and accurate identification of cross-linked peptides of protein complexes, Mol. Cell. Proteomics 10 (2011), https://doi.org/10.1074/mcp.M110.002122 M110.002122.

[104] D.K. Schweppe, J.D. Chavez, C.F. Lee, A. Caudal, S.E. Kruse, R. Stuppard, D.J. Marcinek, G.S. Shadel, R. Tian, J.E. Bruce, Mitochondrial protein interactome elucidated by chemical cross-linking mass spectrometry, Proc. Natl. Acad. Sci. (2017), https://doi.org/10.1073/pnas.1617220114 2017220.

[105] F. Liu, P. Lössl, B.M. Rabbitts, R.S. Balaban, A.J.R. Heck, The interactome of intact mitochondria by cross-linking mass spectrometry provides evidence for coexisting respiratory supercomplexes, Mol. Cell. Proteomics 17 (2018) 216–232, https://doi.org/10.1074/mcp.RA117.000470.

[106] D. Fasci, H. van Ingen, R.A. Scheltema, A.J.R. Heck, Histone interaction landscapes visualized by crosslinking mass spectrometry in intact cell nuclei, Mol. Cell. Proteomics 17 (2018) 2018–2033, https://doi.org/10.1074/mcp.RA118.000924.

[107] A. Leitner, T. Walzthoeni, R. Aebersold, Lysine-specific chemical cross-linking of protein complexes and identification of cross-linking sites using LC-MS/MS and the xQuest/xProphet software pipeline, Nat. Protoc. 9 (2014) 120–137, https://doi.org/10.1038/nprot.2013.168.

[108] A.-C. Gavin, P. Aloy, P. Grandi, R. Krause, M. Boesche, M. Marzioch, C. Rau, L.J. Jensen, S. Bastuck, B. Dümpfeldel, A. Edelmann, M.-A. Heurtier, V. Hoffman, C. Hoefert, K. Klein, M. Hudak, A.-M. Michon, M. Schelder, M. Schirle, M. Remor, T. Rudi, S. Hooper, A. Bauer, T. Bouwmeester, G. Casari, G. Drewes, G. Neubauer, J.M. Rick, B. Kuster, P. Bork, R.B. Russell, S. Superti-Furga, Proteome survey reveals modularity of the yeast cell machinery, Nature 440 (2006) 631–636, https://doi.org/10.1038/nature04532.

[109] N.J. Krogan, G. Cagney, H. Yu, G. Zhong, N. Guo, A. Ignatchenko, J. Li, S. Pu, N. Datta, A.P. Tikuisis, T. Punna, J.M. Peregrín-Alvarez, M. Shales, X. Zhang, M. Davey, M.D. Robinson, A. Paccanaro, J.E. Bray, A. Sheung, B. Beattie, D.P. Richards, V. Canadien, A. Lalev, F. Mena, P. Wong, A. Starostine, M.M. Canete, J. Vlasblom, S. Wu, C. Orsi, S.R. Collins, S. Chandran, R. Haw, J.J. Rilstone, K. Gandi, N.J. Thompson, G. Musso, P. St Ongé, S. Ghanny, M.H.Y. Lam, G. Butland, A.M. Altaf-Ul, S. Kanaya, A. Shilatifard, E. O'Shea, J.S. Weissman, C.J. Ingles, T.R. Hughes, J. Parkinson, M. Gerster, S.J. Wodak, A. Emili, J.F. Greenblatt, Global landscape of protein complexes in the yeast Saccharomyces cerevisiae, Nature 440 (2006) 637–643, https://doi.org/10.1038/.

nature04670.

[110] C. Plaschka, L. Larivière, L. Wenzek, M. Seizl, M. Hemann, D. Tegunov, E.V. Petrotchenko, C.H. Borchers, W. Baumeister, F. Herzog, E. Villa, P. Cramer, Architecture of the RNA polymerase II-Mediator core integration complex, Nature 518 (2015) 376–380, https://doi.org/10.1038/nature14229.

[111] K. Murakami, H. Elmlund, N. Kalisman, D.A. Bushnell, C.M. Adams, M. Azubel, D. Elmlund, Y. Levi-Kalisman, X. Liu, B.J. Gibbons, M. Levitt, R.D. Kornberg, Architecture of an RNA polymerase II transcription pre-initiation complex, Science 342 (2013) 1238724, https://doi.org/10.1126/science.1238724.

[112] C. Bich, S. Maedler, K. Chiesa, F. DeGiacomo, N. Bogliotti, R. Zenobi, Reactivity and applications of new amine reactive cross-linkers for mass spectrometric detection of protein – protein complexes, Anal. Chem. 82 (2010) 172–179, https://doi.org/10.1021/ac091651r.

[113] E.V. Petrotchenko, C.H. Borchers, Crosslinking combined with mass spectrometry for structural proteomics, Mass Spectrom. Rev. 29 (2010) 862–876, https://doi.org/10.1002/mas.20293.

[114] A. Leitner, L.A. Joachimiak, P. Unverdorben, T. Walzthoeni, J. Frydman, F. Förster, R. Aebersold, Chemical cross-linking/mass spectrometry targeting acidic residues in proteins and protein complexes, Proc. Natl. Acad. Sci. 111 (2014) 9455–9460, https://doi.org/10.1073/pnas.1320298111.

[115] D.R. Müller, P. Schindler, H. Towbin, U. Wirth, H. Voshol, S. Hoving, M.O. Steinmetz, Isotope-tagged cross-linking reagents. A new tool in mass spectrometric protein interaction analysis, Anal. Chem. 73 (2001) 1927–1934, https://doi.org/10.1021/ac001379a.

[116] F. Chu, S. Mahrus, C.S. Craik, A.L. Burlingame, Isotope-coded and affinity-tagged cross-linking (ICATXL): An efficient strategy to probe protein interaction surfaces, J. Am. Chem. Soc. 128 (2006) 10362–10363, https://doi.org/10.1021/ja0614159.

[117] M.A. Lauber, J.P. Reilly, Novel amidating cross-linker for facilitating analyses of protein structures and interactions, Anal. Chem. 82 (2010) 7736–7743, https://doi.org/10.1021/ac101586z.

[118] M.J. Trnka, A.L. Burlingame, Topographic studies of the GroEL/GroES chaperonin complex by chemical crosslinking using diforimyl ethynyladenine (DEB): the power of high resolution electron transfer dissociation (ETD) for determination of both peptide sequences and their attachment sites, mcp.M110.003764, Mol. Cell. Proteomics (2010), https://doi.org/10.1074/mcp.M110.003764.

[119] S. Pfammatter, E. Bonneil, F.P. McManus, S. Prasad, D.J. Bailey, M. Belford, J.-J. Dunyach, P. Thibault, A novel differential ion mobility device expands the depth of proteome coverage and the sensitivity of multiplex proteomic measurements, Mol. Cell. Proteomics 17 (2018) 2051–2067, https://doi.org/10.1074/mcp.TIR118.000862.

[120] F. Meier, A.-D. Brunner, S. Koch, H. Koch, M. Lubeck, M. Krause, N. Goedecke, J. Decker, T. Kosinski, M.A. Park, N. Bache, O. Hoernig, J. Cox, O. Räther, M. Mann, Online parallel accumulation-serial fragmentation (PASEF) with a novel trapped ion mobility mass spectrometer, Mol. Cell. Proteomics 17 (2018) 2534–2545, https://doi.org/10.1074/mcp.TIR118.000900.

[121] Mobile and localized protons: a framework for understanding peptide dissociation - Wysocki - 2000 - Journal of Mass Spectrometry - Wiley Online Library, (n.d.), https://onlinelibrary.wiley.com/doi/full/10.1002/1096-9888/28200012%3D2935%3A12%3C1399%3A%3AAD-JMS86%3E3.0.C0%3B2-R (accessed January 9, 2019).

[122] D.M. Good, M. Wirtala, G.C. McAlister, J.J. Coon, Performance characteristics of electron transfer dissociation mass spectrometry, Mol. Cell. Proteomics 6 (2007) 1942–1951, https://doi.org/10.1074/mcp.M700073-MCP200.

[123] L. Earley, L.C. Anderson, D.L. Bai, C. Mullen, J.E.P. Syka, A.M. English, J.-J. Dunyach, G.C. Stafford, J. Shabanowitz, D.F. Hunt, P.D. Compton, Front-end electron transfer dissociation: a new ionization source, Anal. Chem. 85 (2013) 8385–8390, https://doi.org/10.1021/ac401783f.

[124] Target-decoy search strategy for increased confidence in large-scale protein identifications by mass spectrometry | Nature Methods, (n.d.), https://www.nature.com/articles/nmeth1019 (accessed January 10, 2019).

[125] F. Alber, S. Dokudovskaya, L.M. Veenhoff, W. Zhang, J. Kipper, D. Devos, A. Suprapto, O. Karni-Schmidt, R. Williams, B.T. Chait, A. Sali, M.P. Rout, The molecular architecture of the nuclear pore complex, Nature 450 (2007) 695–701, https://doi.org/10.1038/nature06045.

[126] R. Algret, J. Fernandez-Martinez, Y. Shi, S.J. Kim, R. Pellarin, P. Cimermancic, E. Cochet, A. Sali, B.T. Chait, M.P. Rout, S. Dokudovskaya, Molecular architecture and function of the SEA complex, a modulator of the TORC1 pathway, Mol. Cell. Proteomics MCP 13 (2014) 2855–2870, https://doi.org/10.1074/mcp.M114.039388.

[127] S. Viswanath, A. Sali, Optimizing model representation for integrative structure determination of macromolecular assemblies, Proc. Natl. Acad. Sci. U. S. A. 116 (2019) 540–545, https://doi.org/10.1073/pnas.1814649116.

[128] M. Bonomi, S. Hanot, C.H. Greenberg, A. Sali, M. Nilges, M. Vendruscolo, R. Pellarin, Bayesian weighing of electron cryo-microscopy data for integrative structural modeling, Struct. Lond. Engl. 1993 (27) (2019) 175–188.e6, https://doi.org/10.1016/j.str.2018.09.011.

[129] F. Alber, F. Förster, D. Korkin, M. Topf, A. Sali, Integrating diverse data for structure determination of macromolecular assemblies, Annu. Rev. Biochem. 77 (2008) 443–477, https://doi.org/10.1146/annurev.biochem.77.060407.135530.

[130] W. Rieping, M. Habeck, M. Nilges, Inferential structure determination, Science 309 (2005) 303–306, https://doi.org/10.1126/science.1110428.

[131] J.A. Davis, Y. Takagi, R.D. Kornberg, F.A. Asturias, Structure of the yeast RNA polymerase II holoenzyme: Mediator conformation and polymerase interaction, Mol. Cell. 10 (2002) 409–415.

[132] S. Viswanath, I.E. Chemmama, P. Cimermancic, A. Sali, Assessing exhaustiveness
of stochastic sampling for integrative modeling of macromolecular structures, Biophys. J. 113 (2017) 2344–2353, https://doi.org/10.1016/j.bpj.2017.10.005.

[133] E.F. Pettersen, T.D. Goddard, C.C. Huang, G.S. Couch, D.M. Greenblatt, E.C. Meng, T.E. Ferrin, UCSF Chimera—a visualization system for exploratory research and analysis, J. Comput. Chem. 25 (2004) 1605–1612, https://doi.org/10.1002/jcc.20084.

[134] Z. Grabarek, J. Gergely, Zero-length crosslinking procedure with the use of active esters, Anal. Biochem. 185 (1990) 131–135, https://doi.org/10.1016/0003-2697(90)90267-D.

[135] M. Nassal, 4-(1-Azi-2,2,2-trifluoroethyl)benzoic acid, a highly photolabile carbene generating label readily fixable to biochemical agents, Liebigs Ann. Chem. 1983 (1983) 1510–1523, https://doi.org/10.1002/jlac.198319830907.

[136] M. Schneider, A. Belsom, J. Rappsilber, O. Brock, Blind testing of cross-linking/mass spectrometry hybrid methods in CASP11, Proteins Struct. Funct. Bioinforma. 84 (2016) 152–163, https://doi.org/10.1002/prot.25028.

[137] Y. Takizawa, E. Binshtein, A.L. Erwin, T.M. Pyburn, K.F. Mittendorf, M.D. Ohi, While the revolution will not be crystallized, biochemistry reigns supreme, Protein Sci. 26 (2017) 69–81, https://doi.org/10.1002/pro.3054.

[138] L. Fischer, Z.A. Chen, J. Rappsilber, Quantitative cross-linking/mass spectrometry using isotope-labelled cross-linkers, J. Proteomics. 88 (2013) 120–128, https://doi.org/10.1016/j.jprot.2013.03.005.

[139] A general method for targeted quantitative cross-linking mass spectrometry, (n.d.). https://journals.plos.org/plosone/article?id = 10.1371/journal.pone.0167547 (accessed January 8, 2019).

[140] Developing a multiplexed quantitative cross-linking mass spectrometry platform for comparative structural analysis of protein complexes - analytical chemistry (ACS Publications), (n.d.). https://pubs.acs.org/doi/full/10.1021/acs.analchem.6b03148 (accessed January 8, 2019).

[141] P. Tompa, M. Fuxreiter, Fuzzy complexes: polymorphism and structural disorder in protein–protein interactions, Trends Biochem. Sci. 33 (2008) 2–8, https://doi.org/10.1016/j.tibs.2007.10.003.

[142] R. Sharma, Z. Raduly, M. Miskei, M. Fuxreiter, Fuzzy complexes: Specific binding without complete folding, FEBS Lett. 589 (2015) 2533–2542, https://doi.org/10.1016/j.febslet.2015.07.022.

[143] A.J. Wirth, M. Gruebele, Quinary protein structure and the consequences of crowding in living cells: Leaving the test-tube behind, BioEssays 35 (2013) 984–993, https://doi.org/10.1002/bies.201300080.

[144] E.H. McConkey, Molecular evolution, intracellular organization, and the quinary structure of proteins, Proc. Natl. Acad. Sci. 79 (1982) 3236–3240, https://doi.org/10.1073/pnas.79.10.3236.

[145] R.D. Cohen, G.J. Pielak, A cell is more than the sum of its (dilute) parts: A brief

history of quinary structure, Protein Sci. 26 (2017) 403–413, https://doi.org/10.1002/pro.3092.

[146] T. Ikeya, D. Ban, D. Lee, Y. Ito, K. Kato, C. Griesinger, Solution NMR views of dynamical ordering of biomacromolecules, Biochim. Biophys. Acta Gen. Subj. 2018 (1862) 287–306, https://doi.org/10.1016/j.bbagen.2017.08.020.

[147] P. Cuniasse, P. Tavares, E.V. Orlova, S. Zinn-Justin, Structures of biomolecular complexes by combination of NMR and cryoEM methods, Curr. Opin. Struct. Biol. 43 (2017) 104–113, https://doi.org/10.1016/j.sbi.2016.12.008.

[148] T. Ando, S.P. Bhamidimarri, N. Breending, H. Colin-York, L. Collinson, N. De Jonge, P.J. de Pablo, E. Debroye, C. Eggeling, C. Franck, M. Fritzsche, H. Gerritsen, B.N.G. Giepmans, K. Grunewald, J. Hofkens, J.P. Hoogenboom, K.P.F. Janssen, R. Kaufman, J. Klumpermann, N. Kurniawan, J. Kusch, N. Liv, V. Parekh, D.B. Peckys, F. Rehfeldt, D.C. Reutens, M.B.J. Roeffaers, T. Salditt, I.A.T. Schaap, U.S. Schwarz, P. Verkade, M.W. Vogel, R. Wagner, M. Winterhalter, H. Yuan, G. Zifarelli, The 2018 correlative microscopy techniques roadmap, 443001, J. Phys. Appl. Phys. 51 (2018), https://doi.org/10.1088/1361-6463/aad055.

[149] A.-C. Gingras, K.T. Abe, B. Raught, Getting to know the neighborhood: using proximity-dependent biotinylation to characterize protein complexes and map organelles, Curr. Opin. Chem. Biol. 48 (2019) 44–54, https://doi.org/10.1016/j.cbpa.2018.10.017.

[150] L. Wang, M.R. Chance, Protein footprinting comes of age: mass spectrometry for biophysical structure assessment, Mol. Cell. Proteomics MCP 16 (2017) 706–716, https://doi.org/10.1074/mcp.O116.064386.

[151] I.A. Sawyer, J. Bartek, M. Dundr, Phase separated microenvironments inside the cell nucleus are linked to disease and regulate epigenetic state, transcription and RNA processing, Semin. Cell Dev. Biol. (2018), https://doi.org/10.1016/j.semcdb.2018.07.001.

[152] Who's In and Who's Out—Compositional Control of Biomolecular Condensates - ScienceDirect, (n.d.). https://www.sciencedirect.com/science/article/pii/S0022283618309112 (accessed January 6, 2019).

[153] A.G. Larson, D. Elnatan, M.M. Keenen, M.J. Trnka, J.B. Johnston, A.L. Burlingame, D.A. Agard, S. Redding, G.J. Narlikar, Liquid droplet formation by HP1α suggests a role for phase separation in heterochromatin, Nature 547 (2017) 236–240, https://doi.org/10.1038/nature22822.

[154] B.R. Sabari, A. Dall'Agnese, A. Boija, I.A. Klein, E.L. Coffey, K. Shrinivas, B.J. Abraham, N.M. Hannett, A.V. Zamudio, J.C. Manteiga, C.H. Li, Y.E. Guo, D.S. Day, J. Schuijers, E. Vasile, S. Malik, D. Hnisz, T.I. Lee, I.I. Cisse, R.G. Roeder, P.A. Sharp, A.K. Chakraborty, R.A. Young, Coactivator condensation at super-enhancers links phase separation and gene control, Science 361 (2018) eaar3958, https://doi.org/10.1126/science.aar3958.