Comprehensive structure and functional adaptations of the yeast nuclear pore complex.

Nuclear pore complexes (NPCs) mediate the nucleocytoplasmic transport of macromolecules. Here we provide a structure of the isolated yeast NPC in which the inner ring is resolved by cryo-EM at sub-nanometer resolution to show how flexible connectors tie together different structural and functional layers. These connectors may be targets for phosphorylation and regulated disassembly in cells with an open mitosis. Moreover, some nucleoporin pairs and transport factors have similar interaction motifs, which suggests an evolutionary and mechanistic link between assembly and transport. We provide evidence for three major NPC variants that may foreshadow functional specializations at the nuclear periphery. Cryo-electron tomography extended these studies, providing a model of the in situ NPC with a radially expanded inner ring. Our comprehensive model reveals features of the nuclear basket and central transporter, suggests a role for the lumenal Pom152 ring in restricting dilation, and highlights structural plasticity that may be required for transport.

Implications of a multiscale structure of the yeast nuclear pore complex.

Nuclear pore complexes (NPCs) direct the nucleocytoplasmic transport of macromolecules. Here, we provide a composite multiscale structure of the yeast NPC, based on improved 3D density maps from cryogenic electron microscopy and AlphaFold2 models. Key features of the inner and outer rings were integrated into a comprehensive model. We resolved flexible connectors that tie together the core scaffold, along with equatorial transmembrane complexes and a lumenal ring that anchor this channel within the pore membrane. The organization of the nuclear double outer ring reveals an architecture that may be shared with ancestral NPCs. Additional connections between the core scaffold and the central transporter suggest that under certain conditions, a degree of local organization is present at the periphery of the transport machinery. These connectors may couple conformational changes in the scaffold to the central transporter to modulate transport. Collectively, this analysis provides insights into assembly, transport, and NPC evolution.

The structural basis for HIV-1 Vif antagonism of human APOBEC3G.

The APOBEC3 (A3) proteins are host antiviral cellular proteins that hypermutate the viral genome of diverse viral families. In retroviruses, this process requires A3 packaging into viral particles1-4. The lentiviruses encode a protein, Vif, that antagonizes A3 family members by targeting them for degradation. Diversification of A3 allows host escape from Vif whereas adaptations in Vif enable cross-species transmission of primate lentiviruses. How this 'molecular arms race' plays out at the structural level is unknown. Here, we report the cryogenic electron microscopy structure of human APOBEC3G (A3G) bound to HIV-1 Vif, and the hijacked cellular proteins that promote ubiquitin-mediated proteolysis. A small surface explains the molecular arms race, including a cross-species transmission event that led to the birth of HIV-1. Unexpectedly, we find that RNA is a molecular glue for the Vif-A3G interaction, enabling Vif to repress A3G by ubiquitin-dependent and -independent mechanisms. Our results suggest a model in which Vif antagonizes A3G by intercepting it in its most dangerous form for the virus-when bound to RNA and on the pathway to packaging-to prevent viral restriction. By engaging essential surfaces required for restriction, Vif exploits a vulnerability in A3G, suggesting a general mechanism by which RNA binding helps to position key residues necessary for viral antagonism of a host antiviral gene.

From systems to structure - using genetic data to model protein structures.

Understanding the effects of genetic variation is a fundamental problem in biology that requires methods to analyse both physical and functional consequences of sequence changes at systems-wide and mechanistic scales. To achieve a systems view, protein interaction networks map which proteins physically interact, while genetic interaction networks inform on the phenotypic consequences of perturbing these protein interactions. Until recently, understanding the molecular mechanisms that underlie these interactions often required biophysical methods to determine the structures of the proteins involved. The past decade has seen the emergence of new approaches based on coevolution, deep mutational scanning and genome-scale genetic or chemical-genetic interaction mapping that enable modelling of the structures of individual proteins or protein complexes. Here, we review the emerging use of large-scale genetic datasets and deep learning approaches to model protein structures and their interactions, and discuss the integration of structural data from different sources.

Integrative structure and functional anatomy of a nuclear pore complex.

Nuclear pore complexes play central roles as gatekeepers of RNA and protein transport between the cytoplasm and nucleoplasm. However, their large size and dynamic nature have impeded a full structural and functional elucidation. Here we determined the structure of the entire 552-protein nuclear pore complex of the yeast Saccharomyces cerevisiae at sub-nanometre precision by satisfying a wide range of data relating to the molecular arrangement of its constituents. The nuclear pore complex incorporates sturdy diagonal columns and connector cables attached to these columns, imbuing the structure with strength and flexibility. These cables also tie together all other elements of the nuclear pore complex, including membrane-interacting regions, outer rings and RNA-processing platforms. Inwardly directed anchors create a high density of transport factor-docking Phe-Gly repeats in the central channel, organized into distinct functional units. This integrative structure enables us to rationalize the architecture, transport mechanism and evolutionary origins of the nuclear pore complex.

Characterization of an A3G-VifHIV-1-CRL5-CBFß Structure Using a Cross-linking Mass Spectrometry Pipeline for Integrative Modeling of Host-Pathogen Complexes.

Structural analysis of host-pathogen protein complexes remains challenging, largely due to their structural heterogeneity. Here, we describe a pipeline for the structural characterization of these complexes using integrative structure modeling based on chemical cross-links and residue-protein contacts inferred from mutagenesis studies. We used this approach on the HIV-1 Vif protein bound to restriction factor APOBEC3G (A3G), the Cullin-5 E3 ring ligase (CRL5), and the cellular transcription factor Core Binding Factor Beta (CBFß) to determine the structure of the (A3G-Vif-CRL5-CBFß) complex. Using the MS-cleavable DSSO cross-linker to obtain a set of 132 cross-links within this reconstituted complex along with the atomic structures of the subunits and mutagenesis data, we computed an integrative structure model of the heptameric A3G-Vif-CRL5-CBFß complex. The structure, which was validated using a series of tests, reveals that A3G is bound to Vif mostly through its N-terminal domain. Moreover, the model ensemble quantifies the dynamic heterogeneity of the A3G C-terminal domain and Cul5 positions. Finally, the model was used to rationalize previous structural, mutagenesis and functional data not used for modeling, including information related to the A3G-bound and unbound structures as well as mapping functional mutations to the A3G-Vif interface. The experimental and computational approach described here is generally applicable to other challenging host-pathogen protein complexes.

Structural dynamics of the human COP9 signalosome revealed by cross-linking mass spectrometry and integrative modeling.

The COP9 signalosome (CSN) is an evolutionarily conserved eight-subunit (CSN1-8) protein complex that controls protein ubiquitination by deneddylating Cullin-RING E3 ligases (CRLs). The activation and function of CSN hinges on its structural dynamics, which has been challenging to decipher by conventional tools. Here, we have developed a multichemistry cross-linking mass spectrometry approach enabled by three mass spectometry-cleavable cross-linkers to generate highly reliable cross-link data. We applied this approach with integrative structure modeling to determine the interaction and structural dynamics of CSN with the recently discovered ninth subunit, CSN9, in solution. Our results determined the localization of CSN9 binding sites and revealed CSN9-dependent structural changes of CSN. Together with biochemical analysis, we propose a structural model in which CSN9 binding triggers CSN to adopt a configuration that facilitates CSN-CRL interactions, thereby augmenting CSN deneddylase activity. Our integrative structure analysis workflow can be generalized to define in-solution architectures of dynamic protein complexes that remain inaccessible to other approaches.

Na+ coordination at the Na2 site of the Na+/I- symporter.

The sodium/iodide symporter (NIS) mediates active I(-) transport in the thyroid-the first step in thyroid hormone biosynthesis-with a 2 Na(+): 1 I(-) stoichiometry. The two Na(+) binding sites (Na1 and Na2) and the I(-) binding site interact allosterically: when Na(+) binds to a Na(+) site, the affinity of NIS for the other Na(+) and for I(-) increases significantly. In all Na(+)-dependent transporters with the same fold as NIS, the side chains of two residues, S353 and T354 (NIS numbering), were identified as the Na(+) ligands at Na2. To understand the cooperativity between the substrates, we investigated the coordination at the Na2 site. We determined that four other residues-S66, D191, Q194, and Q263-are also involved in Na(+) coordination at this site. Experiments in whole cells demonstrated that these four residues participate in transport by NIS: mutations at these positions result in proteins that, although expressed at the plasma membrane, transport little or no I(-) These residues are conserved throughout the entire SLC5 family, to which NIS belongs, suggesting that they serve a similar function in the other transporters. Our findings also suggest that the increase in affinity that each site displays when an ion binds to another site may result from changes in the dynamics of the transporter. These mechanistic insights deepen our understanding not only of NIS but also of other transporters, including many that, like NIS, are of great medical relevance.

DNA exit ramps are revealed in the binding landscapes obtained from simulations in helical coordinates.

DNA molecules are highly charged semi-flexible polymers that are involved in a wide variety of dynamical processes such as transcription and replication. Characterizing the binding landscapes around DNA molecules is essential to understanding the energetics and kinetics of various biological processes. We present a curvilinear coordinate system that fully takes into account the helical symmetry of a DNA segment. The latter naturally allows to characterize the spatial organization and motions of ligands tracking the minor or major grooves, in a motion reminiscent of sliding. Using this approach, we performed umbrella sampling (US) molecular dynamics (MD) simulations to calculate the three-dimensional potentials of mean force (3D-PMFs) for a Na+ cation and for methyl guanidinium, an arginine analog. The computed PMFs show that, even for small ligands, the free energy landscapes are complex. In general, energy barriers of up to ~5 kcal/mol were measured for removing the ligands from the minor groove, and of ~1.5 kcal/mol for sliding along the minor groove. We shed light on the way the minor groove geometry, defined mainly by the DNA sequence, shapes the binding landscape around DNA, providing heterogeneous environments for recognition by various ligands. For example, we identified the presence of dissociation points or "exit ramps" that naturally would terminate sliding. We discuss how our findings have important implications for understanding how proteins and ligands associate and slide along DNA.

The Acetylation Landscape of the H4 Histone Tail: Disentangling the Interplay between the Specific and Cumulative Effects.

Histone tails, the intrinsically disordered terminal regions of histone proteins, are key modulators of the structure and dynamics of chromatin and, consequently, are central to many DNA template-directed processes including replication, repair, and transcription. Acetylation of histone tails is a major post-translational modification (PTM) involved in regulating chromatin, yet it remains unclear how acetylation modifies the disordered state of histone tails and affects their function. We investigated the consequences of increasing acetylation on the isolated H4 histone tail by characterizing the conformational ensembles of unacetylated, mono-, di-, tri-, and tetra-acetylated H4 histone tails using Replica Exchange Molecular Dynamics (REMD) simulations. We found that progressive acetylation has a cumulative effect on the H4 tail, decreasing conformational heterogeneity, increasing helical propensity, and increasing hydrogen bond occupancies. The monoacetylation of lysine 16, however, has unique and specific effects: drastically decreasing the conformational heterogeneity of the H4 tail and leading to highly localized helical secondary structure and elongated conformations. We describe how the cumulative effects of acetylation arise from the charge reduction and increased hydrophobicity associated with adding acetyl groups, while the specific effects are a consequence of steric interactions that are sequence specific. Additionally, we found that increasing the level of acetylation results in the formation of spatially clustered lysines that could serve as recognition patches for binding of chromatin regulating proteins. Hence, we explore the mechanisms by which different acetylation patterns may result in specific recognition of the H4 histone tails by protein or DNA binding partners.

Concerted dihedral rotations give rise to internal friction in unfolded proteins.

Protein chains undergo conformational diffusion during folding and dynamics, experiencing both thermal kicks and viscous drag. Recent experiments have shown that the corresponding friction can be separated into wet friction, which is determined by the solvent viscosity, and dry friction, where frictional effects arise due to the interactions within the protein chain. Despite important advances, the molecular origins underlying dry friction in proteins have remained unclear. To address this problem, we studied the dynamics of the unfolded cold-shock protein at different solvent viscosities and denaturant concentrations. Using extensive all-atom molecular dynamics simulations we estimated the internal friction time scales and found them to agree well with the corresponding experimental measurements (Soranno et al. Proc. Natl. Acad. Sci. U.S.A. 2012, 109, 17800-17806). Analysis of the reconfiguration dynamics of the unfolded chain further revealed that hops in the dihedral space provide the dominant mechanism of internal friction. Furthermore, the increased number of concerted dihedral moves at physiological conditions suggest that, in such conditions, the concerted motions result in higher frictional forces. These findings have important implications for understanding the folding kinetics of proteins as well as the dynamics of intrinsically disordered proteins.

The Molecular Architecture of the Nuclear Basket.

The nuclear pore complex (NPC) is the sole mediator of nucleocytoplasmic transport. Despite great advances in understanding its conserved core architecture, the peripheral regions can exhibit considerable variation within and between species. One such structure is the cage-like nuclear basket. Despite its crucial roles in mRNA surveillance and chromatin organization, an architectural understanding has remained elusive. Using in-cell cryo-electron tomography and subtomogram analysis, we explored the NPC's structural variations and the nuclear basket across fungi (yeast; S. cerevisiae), mammals (mouse; M. musculus), and protozoa (T. gondii). Using integrative structural modeling, we computed a model of the basket in yeast and mammals that revealed how a hub of Nups in the nuclear ring binds to basket-forming Mlp/Tpr proteins: the coiled-coil domains of Mlp/Tpr form the struts of the basket, while their unstructured termini constitute the basket distal densities, which potentially serve as a docking site for mRNA preprocessing before nucleocytoplasmic transport.

Structural insights into maize viviparous14, a key enzyme in the biosynthesis of the phytohormone abscisic acid.

The key regulatory step in the biosynthesis of abscisic acid (ABA), a hormone central to the regulation of several important processes in plants, is the oxidative cleavage of the 11,12 double bond of a 9-cis-epoxycarotenoid. The enzyme viviparous14 (VP14) performs this cleavage in maize (Zea mays), making it a target for the rational design of novel chemical agents and genetic modifications that improve plant behavior through the modulation of ABA levels. The structure of VP14, determined to 3.2-Å resolution, provides both insight into the determinants of regio- and stereospecificity of this enzyme and suggests a possible mechanism for oxidative cleavage. Furthermore, mutagenesis of the distantly related CCD1 of maize shows how the VP14 structure represents a template for all plant carotenoid cleavage dioxygenases (CCDs). In addition, the structure suggests how VP14 associates with the membrane as a way of gaining access to its membrane soluble substrate.

Integrative structure determination of histones H3 and H4 using genetic interactions.

Integrative structure modeling is increasingly used for determining the architectures of biological assemblies, especially those that are structurally heterogeneous. Recently, we reported on how to convert in vivo genetic interaction measurements into spatial restraints for structural modeling: first, phenotypic profiles are generated for each point mutation and thousands of gene deletions or environmental perturbations. Following, the phenotypic profile similarities are converted into distance restraints on the pairs of mutated residues. We illustrate the approach by determining the structure of the histone H3-H4 complex. The method is implemented in our open-source IMP program, expanding the structural biology toolbox by allowing structural characterization based on in vivo data without the need to purify the target system. We compare genetic interaction measurements to other sources of structural information, such as residue coevolution and deep-learning structure prediction of complex subunits. We also suggest that determining genetic interactions could benefit from new technologies, such as CRISPR-Cas9 approaches to gene editing, especially for mammalian cells. Finally, we highlight the opportunity for using genetic interactions to determine recalcitrant biomolecular structures, such as those of disordered proteins, transient protein assemblies, and host-pathogen protein complexes.

Structure and dynamics of the essential endogenous mycobacterial polyketide synthase Pks13.

Mycobacterium tuberculosis is currently the leading cause of death by any bacterial infection1. The mycolic acid layer of the cell wall is essential for viability and virulence, and the enzymes responsible for its synthesis are therefore front line targets for antimycobacterial drug development2,3. Polyketide synthase 13 (Pks13) is a module comprised of a closely symmetric parallel dimer of chains, each encoding several enzymatic and transport functions, that carries out the condensation of two different very long chain fatty acids to produce mycolic acids that are essential components of the mycobacterial cell wall. Consequently individual enzymatic domains of Pks13 are targets for antimycobacterial drug development4. To understand this machinery, we sought to determine the structure and domain trajectories of the dimeric multi-enzyme Pks13, a 2×198,426 Dalton complex, from protein purified endogenously from mycobacteria under normal growth conditions, to capture it with normal substrates bound trapped 'in action'. Structures of the multi-domain assembly revealed by cryogenic electron microscopy (cryoEM) define the ketosynthase (KS), linker, and acyltransferase (AT) domains, each at atomic resolution (1.8Å), with bound substrates defined at 2.4Å and 2.9Å resolution. Image classification reveals two distinct structures with alternate locations of the N-terminal acyl carrier protein (termed ACP1a, ACP1b) seen at 3.6Å and 4.6Å resolution respectively. These two structures suggest plausible intermediate states, related by a ~60Å movement of ACP1, on the pathway for substrate delivery from the fatty acyl-ACP ligase (FadD32) to the ketosynthase domain. The linking sequence between ACP1 and the KS includes an 11 amino acid sequence with 6 negatively charged side chains that lies in different positively charged grooves on the KS in ACP1a versus ACP1b structures. This charge complementarity between the extended chain and the grooves suggests some stabilization of these two distinct orientations. Other domains are visible at lower resolution and indicate flexibility relative to the KS-AT core. The chemical structures of three bound endogenous long chain fatty acid substrates with their proximal regions defined in the structures were determined by electrospray ionization mass spectrometry. The domain proximities were probed by chemical cross-linking and identified by mass spectrometry. These were incorporated into integrative structure modeling to define multiple domain configurations that transport the very long fatty acid chains throughout the multistep Pks13 mediated synthetic pathway.

Structure and dynamics of the essential endogenous mycobacterial polyketide synthase Pks13.

The mycolic acid layer of the Mycobacterium tuberculosis cell wall is essential for viability and virulence, and the enzymes responsible for its synthesis are targets for antimycobacterial drug development. Polyketide synthase 13 (Pks13) is a module encoding several enzymatic and transport functions that carries out the condensation of two different long-chain fatty acids to produce mycolic acids. We determined structures by cryogenic-electron microscopy of dimeric multi-enzyme Pks13 purified from mycobacteria under normal growth conditions, captured with native substrates. Structures define the ketosynthase (KS), linker and acyl transferase (AT) domains at 1.8 Å resolution and two alternative locations of the N-terminal acyl carrier protein. These structures suggest intermediate states on the pathway for substrate delivery to the KS domain. Other domains, visible at lower resolution, are flexible relative to the KS-AT core. The chemical structures of three bound endogenous long-chain fatty acid substrates were determined by electrospray ionization mass spectrometry.

Disaccharide binding to galectin-1: free energy calculations and molecular recognition mechanism.

Galectin-1, a member of the conserved family of carbohydrate-binding proteins with affinity for ß-galactosides, is a key modulator of diverse cell functions such as immune response and regulation. The binding affinity and specificity of galectin-1 for eight different ß-galactosyl terminal disaccharides was studied using molecular-dynamics simulations in which the ligand was pulled away from the binding site using a mechanical force. We present what we believe to be a novel procedure, based on combinations of multistep trajectories, that was used to estimate the binding free energy (ΔG) of each disaccharide. The computed binding free energy differences show excellent correlation with experimental values determined previously. The small differences in affinity among the disaccharides are the result of an exquisite balance between the strengths of the galectin-sugar H-bonds and the H-bonds the protein and the disaccharides make with the solvent. Analysis of the free energies along the reaction coordinate shows that disaccharide unbinding/binding presents no energetic barrier and, therefore, is diffusion-limited. In addition, the calculations revealed that as the ligand is undocked from the binding site, breaking of protein-disaccharide H-bonds takes place in stages with intermediate states in which the interactions are bridged by water molecules.

Using Integrative Modeling Platform to compute, validate, and archive a model of a protein complex structure.

Biology is advanced by producing structural models of biological systems, such as protein complexes. Some systems are recalcitrant to traditional structure determination methods. In such cases, it may still be possible to produce useful models by integrative structure determination that depends on simultaneous use of multiple types of data. An ensemble of models that are sufficiently consistent with the data is produced by a structural sampling method guided by a data-dependent scoring function. The variation in the ensemble of models quantified the uncertainty of the structure, generally resulting from the uncertainty in the input information and actual structural heterogeneity in the samples used to produce the data. Here, we describe how to generate, assess, and interpret ensembles of integrative structural models using our open source Integrative Modeling Platform program (https://integrativemodeling.org).

Reversible phosphorylation of cyclin T1 promotes assembly and stability of P-TEFb.

The positive transcription elongation factor b (P-TEFb) is a critical coactivator for transcription of most cellular and viral genes, including those of HIV. While P-TEFb is regulated by 7SK snRNA in proliferating cells, P-TEFb is absent due to diminished levels of CycT1 in quiescent and terminally differentiated cells, which has remained unexplored. In these cells, we found that CycT1 not bound to CDK9 is rapidly degraded. Moreover, productive CycT1:CDK9 interactions are increased by PKC-mediated phosphorylation of CycT1 in human cells. Conversely, dephosphorylation of CycT1 by PP1 reverses this process. Thus, PKC inhibitors or removal of PKC by chronic activation results in P-TEFb disassembly and CycT1 degradation. This finding not only recapitulates P-TEFb depletion in resting CD4+ T cells but also in anergic T cells. Importantly, our studies reveal mechanisms of P-TEFb inactivation underlying T cell quiescence, anergy, and exhaustion as well as proviral latency and terminally differentiated cells.

Helix propensities calculations for amino acids in alanine based peptides using Jarzynski's equality.

Jarzynski's equality (Jarzynski, Phys Rev E 1997; 56:5018 and Jarzynski, Phys Rev Lett 1997; 78:2690) relates equilibrium free energy differences between two states A and B to the work done when the system is driven repeatedly and irreversibly from an equilibrium state A to equilibrium state B. We present calculations of helix propensities using a novel procedure based on this equality. In particular, a work probability distribution is built based on combinations of multi-step trajectories that give representative work distributions without requiring computing an unreasonable large number of trajectories between states. A small number of trajectories (15) are used to construct a distribution that contains millions of work values. This distribution is used to calculate DeltaG(AB) using Jarsynski's equation. To apply and test this method, we used as a model system a dodeca-alanine helix, analyzing its extension using mechanical force. This helix was used as the basis of a host guest system in which two of the 12 residues are substituted by some other amino acid (as the guests). The differences between the unfolding free energies of the substituted peptides and the all-alanine peptide provided values for DeltaDeltaG that can be interpreted as the helix propensities of each amino acid. Results show good correlation with the experimental measurements of Baldwin and coworkers (Chakrabartty et al., Protein Sci 1994; 3:843-852).