Integrative spatiotemporal map of nucleocytoplasmic transport.

The Nuclear Pore Complex (NPC) facilitates rapid and selective nucleocytoplasmic transport of molecules as large as ribosomal subunits and viral capsids. It is not clear how key emergent properties of this transport arise from the system components and their interactions. To address this question, we constructed an integrative coarse-grained Brownian dynamics model of transport through a single NPC, followed by coupling it with a kinetic model of Ran-dependent transport in an entire cell. The microscopic model parameters were fitted to reflect experimental data and theoretical information regarding the transport, without making any assumptions about its emergent properties. The resulting reductionist model is validated by reproducing several features of transport not used for its construction, such as the morphology of the central transporter, rates of passive and facilitated diffusion as a function of size and valency, in situ radial distributions of pre-ribosomal subunits, and active transport rates for viral capsids. The model suggests that the NPC functions essentially as a virtual gate whose flexible phenylalanine-glycine (FG) repeat proteins raise an entropy barrier to diffusion through the pore. Importantly, this core functionality is greatly enhanced by several key design features, including 'fuzzy' and transient interactions, multivalency, redundancy in the copy number of FG nucleoporins, exponential coupling of transport kinetics and thermodynamics in accordance with the transition state theory, and coupling to the energy-reliant RanGTP concentration gradient. These design features result in the robust and resilient rate and selectivity of transport for a wide array of cargo ranging from a few kilodaltons to megadaltons in size. By dissecting these features, our model provides a quantitative starting point for rationally modulating the transport system and its artificial mimics.

Dynamic molecular mechanism of the nuclear pore complex permeability barrier.

Nuclear pore complexes (NPCs) mediate nucleocytoplasmic transport of specific macromolecules while impeding the exchange of unsolicited material. However, key aspects of this gating mechanism remain controversial. To address this issue, we determined the nanoscopic behavior of the permeability barrier directly within yeast S. cerevisiae NPCs at transport-relevant timescales. We show that the large intrinsically disordered domains of phenylalanine-glycine repeat nucleoporins (FG Nups) exhibit highly dynamic fluctuations to create transient voids in the permeability barrier that continuously shape-shift and reseal, resembling a radial polymer brush. Together with cargo-carrying transport factors the FG domains form a feature called the central plug, which is also highly dynamic. Remarkably, NPC mutants with longer FG domains show interweaving meshwork-like behavior that attenuates nucleocytoplasmic transport in vivo. Importantly, the bona fide nanoscale NPC behaviors and morphologies are not recapitulated by in vitro FG domain hydrogels. NPCs also exclude self-assembling FG domain condensates in vivo, thereby indicating that the permeability barrier is not generated by a self-assembling phase condensate, but rather is largely a polymer brush, organized by the NPC scaffold, whose dynamic gating selectivity is strongly enhanced by the presence of transport factors.

Familial Early-Onset Alzheimer's Caused by Novel Genetic Variant and APP Duplication: A Cross-Sectional Study.

BACKGROUND: The clinical characteristics of symptomatic and asymptomatic carriers of early- onset autosomal dominant Alzheimer's (EOADAD) due to a yet-undescribed chromosomal rearrangement may add to the available body of knowledge about Alzheimer's disease and may enlighten novel and modifier genes. We report the clinical and genetic characteristics of asymptomatic and symptomatic individuals carrying a novel APP duplication rearrangement. METHODS: Individuals belonging to a seven-generation pedigree with familial cognitive decline or intracerebral hemorrhages were recruited. Participants underwent medical, neurological, and neuropsychological evaluations. The genetic analysis included chromosomal microarray, Karyotype, fluorescence in situ hybridization, and whole genome sequencing. RESULTS: Of 68 individuals, six females presented with dementia, and four males presented with intracerebral hemorrhage. Of these, nine were found to carry Chromosome 21 copy number gain (chr21:27,224,097-27,871,284, GRCh37/hg19) including the APP locus (APP-dup). In seven, Chromosome 5 copy number gain (Chr5: 24,786,234-29,446,070, GRCh37/hg19) (Chr5-CNG) cosegregated with the APP-dup. Both duplications co-localized to chromosome 18q21.1 and segregated in 25 pre-symptomatic carriers. Compared to non-carriers, asymptomatic carriers manifested cognitive decline in their mid-thirties. A third of the affected individuals carried a diagnosis of a dis-immune condition. CONCLUSION: APP extra dosage, even in isolation and when located outside chromosome 21, is pathogenic. The clinical presentation of APP duplication varies and may be gender specific, i.e., ICH in males and cognitive-behavioral deterioration in females. The association with immune disorders is presently unclear but may prove relevant. The implication of Chr5-CNG co-segregation and the surrounding chromosome 18 genetic sequence needs further clarification.

Mechanical force application to the nucleus regulates nucleocytoplasmic transport.

Mechanical force controls fundamental cellular processes in health and disease, and increasing evidence shows that the nucleus both experiences and senses applied forces. Such forces can lead to the nuclear translocation of proteins, but whether force controls nucleocytoplasmic transport, and how, remains unknown. Here we show that nuclear forces differentially control passive and facilitated nucleocytoplasmic transport, setting the rules for the mechanosensitivity of shuttling proteins. We demonstrate that nuclear force increases permeability across nuclear pore complexes, with a dependence on molecular weight that is stronger for passive than for facilitated diffusion. Owing to this differential effect, force leads to the translocation of cargoes into or out of the nucleus within a given range of molecular weight and affinity for nuclear transport receptors. Further, we show that the mechanosensitivity of several transcriptional regulators can be both explained by this mechanism and engineered exogenously by introducing appropriate nuclear localization signals. Our work unveils a mechanism of mechanically induced signalling, probably operating in parallel with others, with potential applicability across signalling pathways.

Bayesian metamodeling of complex biological systems across varying representations.

Comprehensive modeling of a whole cell requires an integration of vast amounts of information on various aspects of the cell and its parts. To divide and conquer this task, we introduce Bayesian metamodeling, a general approach to modeling complex systems by integrating a collection of heterogeneous input models. Each input model can in principle be based on any type of data and can describe a different aspect of the modeled system using any mathematical representation, scale, and level of granularity. These input models are 1) converted to a standardized statistical representation relying on probabilistic graphical models, 2) coupled by modeling their mutual relations with the physical world, and 3) finally harmonized with respect to each other. To illustrate Bayesian metamodeling, we provide a proof-of-principle metamodel of glucose-stimulated insulin secretion by human pancreatic ß-cells. The input models include a coarse-grained spatiotemporal simulation of insulin vesicle trafficking, docking, and exocytosis; a molecular network model of glucose-stimulated insulin secretion signaling; a network model of insulin metabolism; a structural model of glucagon-like peptide-1 receptor activation; a linear model of a pancreatic cell population; and ordinary differential equations for systemic postprandial insulin response. Metamodeling benefits from decentralized computing, while often producing a more accurate, precise, and complete model that contextualizes input models as well as resolves conflicting information. We anticipate Bayesian metamodeling will facilitate collaborative science by providing a framework for sharing expertise, resources, data, and models, as exemplified by the Pancreatic ß-Cell Consortium.

Two separate functions of NME3 critical for cell survival underlie a neurodegenerative disorder.

We report a patient who presented with congenital hypotonia, hypoventilation, and cerebellar histopathological alterations. Exome analysis revealed a homozygous mutation in the initiation codon of the NME3 gene, which encodes an NDP kinase. The initiation-codon mutation leads to deficiency in NME3 protein expression. NME3 is a mitochondrial outer-membrane protein capable of interacting with MFN1/2, and its depletion causes dysfunction in mitochondrial dynamics. Consistently, the patient's fibroblasts were characterized by a slow rate of mitochondrial dynamics, which was reversed by expression of wild-type or catalytic-dead NME3. Moreover, glucose starvation caused mitochondrial fragmentation and cell death in the patient's cells. The expression of wild-type and catalytic-dead but not oligomerization-attenuated NME3 restored mitochondrial elongation. However, only wild-type NME3 sustained ATP production and viability. Thus, the separate functions of NME3 in mitochondrial fusion and NDP kinase cooperate in metabolic adaptation for cell survival in response to glucose starvation. Given the critical role of mitochondrial dynamics and energy requirements in neuronal development, the homozygous mutation in NME3 is linked to a fatal mitochondrial neurodegenerative disorder.

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.

Structural basis for chemokine recognition by a G protein-coupled receptor and implications for receptor activation.

Chemokines orchestrate cell migration for development, immune surveillance, and disease by binding to cell surface heterotrimeric guanine nucleotide-binding protein (G protein)-coupled receptors (GPCRs). The array of interactions between the nearly 50 chemokines and their 20 GPCR targets generates an extensive signaling network to which promiscuity and biased agonism add further complexity. The receptor CXCR4 recognizes both monomeric and dimeric forms of the chemokine CXCL12, which is a distinct example of ligand bias in the chemokine family. We demonstrated that a constitutively monomeric CXCL12 variant reproduced the G protein-dependent and ß-arrestin-dependent responses that are associated with normal CXCR4 signaling and lead to cell migration. In addition, monomeric CXCL12 made specific contacts with CXCR4 that are not present in the structure of the receptor in complex with a dimeric form of CXCL12, a biased agonist that stimulates only G protein-dependent signaling. We produced an experimentally validated model of an agonist-bound chemokine receptor that merged a nuclear magnetic resonance-based structure of monomeric CXCL12 bound to the amino terminus of CXCR4 with a crystal structure of the transmembrane domains of CXCR4. The large CXCL12:CXCR4 protein-protein interface revealed by this structure identified previously uncharacterized functional interactions that fall outside of the classical "two-site model" for chemokine-receptor recognition. Our model suggests a mechanistic hypothesis for how interactions on the extracellular face of the receptor may stimulate the conformational changes required for chemokine receptor-mediated signal transduction.

Simple rules for passive diffusion through the nuclear pore complex.

Passive macromolecular diffusion through nuclear pore complexes (NPCs) is thought to decrease dramatically beyond a 30-60-kD size threshold. Using thousands of independent time-resolved fluorescence microscopy measurements in vivo, we show that the NPC lacks such a firm size threshold; instead, it forms a soft barrier to passive diffusion that intensifies gradually with increasing molecular mass in both the wild-type and mutant strains with various subsets of phenylalanine-glycine (FG) domains and different levels of baseline passive permeability. Brownian dynamics simulations replicate these findings and indicate that the soft barrier results from the highly dynamic FG repeat domains and the diffusing macromolecules mutually constraining and competing for available volume in the interior of the NPC, setting up entropic repulsion forces. We found that FG domains with exceptionally high net charge and low hydropathy near the cytoplasmic end of the central channel contribute more strongly to obstruction of passive diffusion than to facilitated transport, revealing a compartmentalized functional arrangement within the NPC.

Slide-and-exchange mechanism for rapid and selective transport through the nuclear pore complex.

Nucleocytoplasmic transport is mediated by the interaction of transport factors (TFs) with disordered phenylalanine-glycine (FG) repeats that fill the central channel of the nuclear pore complex (NPC). However, the mechanism by which TFs rapidly diffuse through multiple FG repeats without compromising NPC selectivity is not yet fully understood. In this study, we build on our recent NMR investigations showing that FG repeats are highly dynamic, flexible, and rapidly exchanging among TF interaction sites. We use unbiased long timescale all-atom simulations on the Anton supercomputer, combined with extensive enhanced sampling simulations and NMR experiments, to characterize the thermodynamic and kinetic properties of FG repeats and their interaction with a model transport factor. Both the simulations and experimental data indicate that FG repeats are highly dynamic random coils, lack intrachain interactions, and exhibit significant entropically driven resistance to spatial confinement. We show that the FG motifs reversibly slide in and out of multiple TF interaction sites, transitioning rapidly between a strongly interacting state and a weakly interacting state, rather than undergoing a much slower transition between strongly interacting and completely noninteracting (unbound) states. In the weakly interacting state, FG motifs can be more easily displaced by other competing FG motifs, providing a simple mechanism for rapid exchange of TF/FG motif contacts during transport. This slide-and-exchange mechanism highlights the direct role of the disorder within FG repeats in nucleocytoplasmic transport, and resolves the apparent conflict between the selectivity and speed of transport.

A backdoor to the nucleus that runs in the family?

Yoshimura and colleagues show that HEAT proteins that are involved in diverse cellular functions may facilitate their own translocation through the nuclear pore complex, owing to their structural similarity to nuclear transport receptors of the karyopherin ß family.

Structural basis of eukaryotic cell-cell fusion.

Cell-cell fusion proteins are essential in development. Here we show that the C. elegans cell-cell fusion protein EFF-1 is structurally homologous to viral class II fusion proteins. The 2.6 Å crystal structure of the EFF-1 trimer displays the same 3D fold and quaternary conformation of postfusion class II viral fusion proteins, although it lacks a nonpolar "fusion loop," indicating that it does not insert into the target membrane. EFF-1 was previously shown to be required in both cells for fusion, and we show that blocking EFF-1 trimerization blocks the fusion reaction. Together, these data suggest that whereas membrane fusion driven by viral proteins entails leveraging of a nonpolar loop, EFF-1-driven fusion of cells entails trans-trimerization such that transmembrane segments anchored in the two opposing membranes are brought into contact at the tip of the EFF-1 trimer to then, analogous to SNARE-mediated vesicle fusion, zip the two membranes into one.

Modeling of proteins and their assemblies with the Integrative Modeling Platform.

To understand the workings of the living cell, we need to characterize protein assemblies that constitute the cell (for example, the ribosome, 26S proteasome, and the nuclear pore complex). A reliable high-resolution structural characterization of these assemblies is frequently beyond the reach of current experimental methods, such as X-ray crystallography, NMR spectroscopy, electron microscopy, footprinting, chemical cross-linking, FRET spectroscopy, small angle X-ray scattering, and proteomics. However, the information garnered from different methods can be combined and used to build models of the assembly structures that are consistent with all of the available datasets, and therefore more accurate, precise, and complete. Here, we describe a protocol for this integration, whereby the information is converted to a set of spatial restraints and a variety of optimization procedures can be used to generate models that satisfy the restraints as well as possible. These generated models can then potentially inform about the precision and accuracy of structure determination, the accuracy of the input datasets, and further data generation. We also demonstrate the Integrative Modeling Platform (IMP) software, which provides the necessary computational framework to implement this protocol, and several applications for specific use cases.

Druggable protein-protein interactions--from hot spots to hot segments.

Protein-Protein Interactions (PPIs) mediate numerous biological functions. As such, the inhibition of specific PPIs has tremendous therapeutic value. The notion that these interactions are 'undruggable' has petered out with the emergence of more and more successful examples of PPI inhibitors, expanding considerably the scope of potential drug targets. The accumulated data on successes in the inhibition of PPIs allow us to analyze the features that are required for such inhibition. Whereas it has been suggested and shown that targeting hot spots at PPI interfaces is a good strategy to achieve inhibition, in this review we focus on the notion that the most amenable interactions for inhibition are those that are mediated by a 'hot segment', a continuous epitope that contributes the majority of the binding energy. This criterion is both useful in guiding future target selection efforts, and in suggesting immediate inhibitory candidates--the dominant peptidic segment that mediates the targeted interaction.

Peptide docking and structure-based characterization of peptide binding: from knowledge to know-how.

Peptide-mediated interactions are gaining increased attention due to their predominant roles in the many regulatory processes that involve dynamic interactions between proteins. The structures of such interactions provide an excellent starting point for their characterization and manipulation, and can provide leads for targeted inhibitor design. The relatively few experimentally determined structures of peptide-protein complexes can be complemented with an outburst of modeling approaches that have been introduced in recent years, with increasing accuracy and applicability to ever more systems. We review different methods to address the considerable challenges in modeling the binding of a short yet highly flexible peptide to its partner. These methods apply an array of sampling strategies and draw from a recent amassing of knowledge about the biophysical nature of peptide-protein interactions. We elaborate on applications of these structure-based approaches and in particular on the characterization of peptide binding specificity to different peptide-binding domains and enzymes. Such applications can identify new biological targets and thus complement our current view of protein-protein interactions in living organisms. Accurate peptide-protein docking is of particular importance in the light of increased appreciation of the crucial functional roles of disordered regions and the many linear binding motifs embedded within.

Modeling peptide-protein interactions.

Peptide-protein interactions are prevalent in the living cell and form a key component of the overall protein-protein interaction network. These interactions are drawing increasing interest due to their part in signaling and regulation, and are thus attractive targets for computational structural modeling. Here we report an overview of current techniques for the high resolution modeling of peptide-protein complexes. We dissect this complicated challenge into several smaller subproblems, namely: modeling the receptor protein, predicting the peptide binding site, sampling an initial peptide backbone conformation and the final refinement of the peptide within the receptor binding site. For each of these conceptual stages, we present available tools, approaches, and their reported performance. We summarize with an illustrative example of this process, highlighting the success and current challenges still facing the automated blind modeling of peptide-protein interactions. We believe that the upcoming years will see considerable progress in our ability to create accurate models of peptide-protein interactions, with applications in binding-specificity prediction, rational design of peptide-mediated interactions and the usage of peptides as therapeutic agents.

Rosetta FlexPepDock web server--high resolution modeling of peptide-protein interactions.

Peptide-protein interactions are among the most prevalent and important interactions in the cell, but a large fraction of those interactions lack detailed structural characterization. The Rosetta FlexPepDock web server (http://flexpepdock.furmanlab.cs.huji.ac.il/) provides an interface to a high-resolution peptide docking (refinement) protocol for the modeling of peptide-protein complexes, implemented within the Rosetta framework. Given a protein receptor structure and an approximate, possibly inaccurate model of the peptide within the receptor binding site, the FlexPepDock server refines the peptide to high resolution, allowing full flexibility to the peptide backbone and to all side chains. This protocol was extensively tested and benchmarked on a wide array of non-redundant peptide-protein complexes, and was proven effective when applied to peptide starting conformations within 5.5 Å backbone root mean square deviation from the native conformation. FlexPepDock has been applied to several systems that are mediated and regulated by peptide-protein interactions. This easy to use and general web server interface allows non-expert users to accurately model their specific peptide-protein interaction of interest.

Rosetta FlexPepDock ab-initio: simultaneous folding, docking and refinement of peptides onto their receptors.

Flexible peptides that fold upon binding to another protein molecule mediate a large number of regulatory interactions in the living cell and may provide highly specific recognition modules. We present Rosetta FlexPepDock ab-initio, a protocol for simultaneous docking and de-novo folding of peptides, starting from an approximate specification of the peptide binding site. Using the Rosetta fragments library and a coarse-grained structural representation of the peptide and the receptor, FlexPepDock ab-initio samples efficiently and simultaneously the space of possible peptide backbone conformations and rigid-body orientations over the receptor surface of a given binding site. The subsequent all-atom refinement of the coarse-grained models includes full side-chain modeling of both the receptor and the peptide, resulting in high-resolution models in which key side-chain interactions are recapitulated. The protocol was applied to a benchmark in which peptides were modeled over receptors in either their bound backbone conformations or in their free, unbound form. Near-native peptide conformations were identified in 18/26 of the bound cases and 7/14 of the unbound cases. The protocol performs well on peptides from various classes of secondary structures, including coiled peptides with unusual turns and kinks. The results presented here significantly extend the scope of state-of-the-art methods for high-resolution peptide modeling, which can now be applied to a wide variety of peptide-protein interactions where no prior information about the peptide backbone conformation is available, enabling detailed structure-based studies and manipulation of those interactions.

Autophosphorylation activates Dictyostelium myosin II heavy chain kinase A by providing a ligand for an allosteric binding site in the alpha-kinase domain.

Dictyostelium discoideum myosin II heavy chain kinase A (MHCK A), a member of the atypical α-kinase family, phosphorylates sites in the myosin II tail that block filament assembly. Here we show that the catalytic activity of A-CAT, the α-kinase domain of MHCK A (residues 552-841), is severely inhibited by the removal of a disordered C-terminal tail sequence (C-tail; residues 806-841). The key residue in the C-tail was identified as Thr(825), which was found to be constitutively autophosphorylated. Dephosphorylation of Thr(825) using shrimp alkaline phosphatase decreased A-CAT activity. The activity of a truncated A-CAT lacking Thr(825) could be rescued by P(i), phosphothreonine, and a phosphorylated peptide, but not by threonine, glutamic acid, aspartic acid, or an unphosphorylated peptide. These results focused attention on a P(i)-binding pocket located in the C-terminal lobe of A-CAT. Mutational analysis demonstrated that the P(i)-pocket was essential for A-CAT activity. Based on these results, it is proposed that autophosphorylation of Thr(825) activates ACAT by providing a covalently tethered ligand for the P(i)-pocket. Ab initio modeling studies using the Rosetta FloppyTail and FlexPepDock protocols showed that it is feasible for the phosphorylated Thr(825) to dock intramolecularly into the P(i)-pocket. Allosteric activation is predicted to involve a conformational change in Arg(734), which bridges the bound P(i) to Asp(762) in a key active site loop. Sequence alignments indicate that a comparable regulatory mechanism is likely to be conserved in Dictyostelium MHCK B-D and metazoan eukaryotic elongation factor-2 kinases.

Allosteric regulation of glycogen synthase kinase 3ß: a theoretical study.

Glycogen synthase kinase 3ß (GSK-3ß) is a serine-threonine kinase belonging to the CMGC family that plays a key role in many biological processes, such as glucose metabolism, cell cycle regulation, and proliferation. Like most protein kinases, GSK-3ß is regulated via multiple pathways and sites. We performed all-atom molecular dynamics simulations on the unphosphorylated and phosphorylated unbound GSK-3ß and the phosphorylated GSK-3ß bound to a peptide substrate, its product, and a derived inhibitor. We found that GSK-3ß autophosphorylation at residue Tyr(216) results in widening of the catalytic groove, thereby facilitating substrate access. In addition, we studied the interactions of the phosphorylated GSK-3ß with a substrate and peptide inhibitor located at the active site and observed higher affinity of the inhibitor to the kinase. Furthermore, we detected a potential remote binding site which was previously identified in other kinases. In agreement with experiments we observed that binding of specific peptides at this remote site leads to stabilization of the activation loop located in the active site. We speculate that this stabilization could enhance the catalytic activity of the kinase. We point to this remote site as being structurally conserved and suggest that the allosteric phenomenon observed here may occur in the protein kinase superfamily.