Integrative genome modeling platform reveals essentiality of rare contact events in 3D genome organizations.

A multitude of sequencing-based and microscopy technologies provide the means to unravel the relationship between the three-dimensional organization of genomes and key regulatory processes of genome function. Here, we develop a multimodal data integration approach to produce populations of single-cell genome structures that are highly predictive for nuclear locations of genes and nuclear bodies, local chromatin compaction and spatial segregation of functionally related chromatin. We demonstrate that multimodal data integration can compensate for systematic errors in some of the data and can greatly increase accuracy and coverage of genome structure models. We also show that alternative combinations of different orthogonal data sources can converge to models with similar predictive power. Moreover, our study reveals the key contributions of low-frequency ('rare') interchromosomal contacts to accurately predicting the global nuclear architecture, including the positioning of genes and chromosomes. Overall, our results highlight the benefits of multimodal data integration for genome structure analysis, available through the Integrative Genome Modeling software package.

KymoKnot: A web server and software package to identify and locate knots in trajectories of linear or circular polymers.

The KymoKnot software package and web server identifies and locates physical knots or proper knots in a series of polymer conformations. It is mainly intended as an analysis tool for trajectories of linear or circular polymers, but it can be used on single instances too, e.g. protein structures in PDB format. A key element of the software package is the so-called minimally interfering chain closure algorithm that is used to detect physical knots in open chains and to locate the knotted region in both open and closed chains. The web server offers a user-friendly graphical interface that identifies the knot type and highlights the knotted region on each frame of the trajectory, which the user can visualize interactively from various viewpoints. The dynamical evolution of the knotted region along the chain contour is presented as a kymograph. All data can be downloaded in text format. The KymoKnot package is licensed under the BSD 3-Clause licence. The server is publicly available at http://kymoknot.sissa.it/kymoknot/interactive.php .

Mechanical and assembly units of viral capsids identified via quasi-rigid domain decomposition.

Key steps in a viral life-cycle, such as self-assembly of a protective protein container or in some cases also subsequent maturation events, are governed by the interplay of physico-chemical mechanisms involving various spatial and temporal scales. These salient aspects of a viral life cycle are hence well described and rationalised from a mesoscopic perspective. Accordingly, various experimental and computational efforts have been directed towards identifying the fundamental building blocks that are instrumental for the mechanical response, or constitute the assembly units, of a few specific viral shells. Motivated by these earlier studies we introduce and apply a general and efficient computational scheme for identifying the stable domains of a given viral capsid. The method is based on elastic network models and quasi-rigid domain decomposition. It is first applied to a heterogeneous set of well-characterized viruses (CCMV, MS2, STNV, STMV) for which the known mechanical or assembly domains are correctly identified. The validated method is next applied to other viral particles such as L-A, Pariacoto and polyoma viruses, whose fundamental functional domains are still unknown or debated and for which we formulate verifiable predictions. The numerical code implementing the domain decomposition strategy is made freely available.

The 3.5 ångström X-ray structure of the human connexin26 gap junction channel is unlikely that of a fully open channel.

The permeability of gap junction channels to metabolites, and not simply to small inorganic ions, is likely to play an important role in development, physiology as well as in etiology of several diseases. Here, we combined dual patch clamp and fluorescence imaging techniques with molecular dynamics (MD) simulations to investigate the permeation of calcein, a relatively large fluorescent tracer (MW 622 Da) through homomeric gap junction channels formed by wild type human connexin26 (hCx26wt) protomers. Our experimental data indicate that the unitary flux of calcein driven by a 125 µM concentration difference is Jpore = 226 molecule/s per channel. In the light of Eyring transition state theory adapted for the liquid phase, this value corresponds to an energy barrier of ~20 kBT (where kB is the Boltzmann constant and T is absolute temperature). The barrier predicted by our MD simulations, based on the 3.5 Å X-ray structural model of the hCx26wt gap junction channel, is ~45 kBT. The main contributions to the energetics of calcein permeation originated from the interaction between the permeating molecule and the charged aminoacids lining the channel pore. Assigning a fake zero total charge to the calcein molecule yielded a value for the barrier height compatible with the experimental data. These results can be accounted for by two different (although not mutually exclusive) hypotheses: (1) the X-ray model of the hCx26wt gap junction channel is not representative of a fully open state; (2) post translational modifications affecting the hCx26wt protein in our expression system differed from the modifications undergone by the proteins in the conditions used to obtain the crystal structure. Hypothesis (1) is compatible with data indicating that, only 10% or less of the channels forming a gap junction plaque are in the open state, and therefore the averaging procedure intrinsic in the generation of the crystal structure data more closely reflects that of a closed channel. Hypothesis (2) is compatible with recent mass spectrometry data and implies that the charge of several amino acid side chains may have been altered, thus modifying substantially the permeation properties of the channels in living cells.

Evaluating the role of the nuclear microenvironment in gene function by population-based modeling.

The nuclear folding of chromosomes relative to nuclear bodies is an integral part of gene function. Here, we demonstrate that population-based modeling-from ensemble Hi-C data-provides a detailed description of the nuclear microenvironment of genes and its role in gene function. We define the microenvironment by the subnuclear positions of genomic regions with respect to nuclear bodies, local chromatin compaction, and preferences in chromatin compartmentalization. These structural descriptors are determined in single-cell models, thereby revealing the structural variability between cells. We demonstrate that the microenvironment of a genomic region is linked to its functional potential in gene transcription, replication, and chromatin compartmentalization. Some chromatin regions feature a strong preference for a single microenvironment, due to association with specific nuclear bodies in most cells. Other chromatin shows high structural variability, which is a strong indicator of functional heterogeneity. Moreover, we identify specialized nuclear microenvironments, which distinguish chromatin in different functional states and reveal a key role of nuclear speckles in chromosome organization. We demonstrate that our method produces highly predictive three-dimensional genome structures, which accurately reproduce data from a variety of orthogonal experiments, thus considerably expanding the range of Hi-C data analysis.

Quantitative imaging of chromatin decompaction in living cells.

Chromatin organization is highly dynamic and regulates transcription. Upon transcriptional activation, chromatin is remodeled and referred to as "open," but quantitative and dynamic data of this decompaction process are lacking. Here, we have developed a quantitative high resolution-microscopy assay in living yeast cells to visualize and quantify chromatin dynamics using the GAL7-10-1 locus as a model system. Upon transcriptional activation of these three clustered genes, we detect an increase of the mean distance across this locus by >100 nm. This decompaction is linked to active transcription but is not sensitive to the histone deacetylase inhibitor trichostatin A or to deletion of the histone acetyl transferase Gcn5. In contrast, the deletion of SNF2 (encoding the ATPase of the SWI/SNF chromatin remodeling complex) or the deactivation of the histone chaperone complex FACT lead to a strongly reduced decompaction without significant effects on transcriptional induction in FACT mutants. Our findings are consistent with nucleosome remodeling and eviction activities being major contributors to chromatin reorganization during transcription but also suggest that transcription can occur in the absence of detectable decompaction.

Integrative modelling of cellular assemblies.

A wide variety of experimental techniques can be used for understanding the precise molecular mechanisms underlying the activities of cellular assemblies. The inherent limitations of a single experimental technique often requires integration of data from complementary approaches to gain sufficient insights into the assembly structure and function. Here, we review popular computational approaches for integrative modelling of cellular assemblies, including protein complexes and genomic assemblies. We provide recent examples of integrative models generated for such assemblies by different experimental techniques, especially including data from 3D electron microscopy (3D-EM) and chromosome conformation capture experiments, respectively. We highlight general concepts in integrative modelling and discuss the need for careful formulation and merging of different types of information.

Optimal Self-Assembly of Linked Constructs and Catenanes via Spatial Confinement.

How to direct the self-assembly of simple templates toward constructs with complex shape and topology is still an open problem. Recent advancements have made it possible to self-assemble various types of knotted constructs, but targeting general multicomponent topologies, that is, links and catenanes, has proved much harder. Here, we study how the yield and complexity of self-assembled links depends on both intrinsic and extrinsic system properties and, particularly, template shape and spatial confinement. We show that slit confinement does not necessarily suppress linking but can rather enhance it significantly thanks to entropic effects. We also found that only a limited set of binary links are recurrent for different template shapes. These privileged topologies include all those experimentally realized so far plus a few additional ones, such as the 772 and 782 links that, hence, ought to be ideal candidates for broadening the current class of constructs with addressable topology.

Self-assembling knots of controlled topology by designing the geometry of patchy templates.

The self-assembly of objects with a set of desired properties is a major goal of material science and physics. A particularly challenging problem is that of self-assembling structures with a target topology. Here we show by computer simulation that one may design the geometry of string-like rigid patchy templates to promote their efficient and reproducible self-assembly into a selected repertoire of non-planar closed folds including several knots. In particular, by controlling the template geometry, we can direct the assembly process so as to strongly favour the formation of constructs tied in trefoil or pentafoil, or even of more exotic torus knots. Polydisperse and racemic mixtures of helical fragments of variable composition add further tunability in the topological self-assembly we discovered. Our results should be relevant to the design of new ways to synthesize molecular knots, which may prove, for instance, to be efficient cargo-carriers due to their mechanical stability.

SPECTRUS: A Dimensionality Reduction Approach for Identifying Dynamical Domains in Protein Complexes from Limited Structural Datasets.

Identifying dynamical, quasi-rigid domains in proteins provides a powerful means for characterizing functionally oriented structural changes via a parsimonious set of degrees of freedom. In fact, the relative displacements of few dynamical domains usually suffice to rationalize the mechanics underpinning biological functionality in proteins and can even be exploited for structure determination or refinement purposes. Here we present SPECTRUS, a general scheme that, by solely using amino acid distance fluctuations, can pinpoint the innate quasi-rigid domains of single proteins or large complexes in a robust way. Consistent domains are usually obtained by using either a pair of representative structures or thousands of conformers. The functional insights offered by the approach are illustrated for biomolecular systems of very different size and complexity such as kinases, ion channels, and viral capsids. The decomposition tool is available as a software package and web server at spectrus.sissa.it.

Permeation pathway of homomeric connexin 26 and connexin 30 channels investigated by molecular dynamics.

Mutations in the genes GJB2 and GJB6 encoding human connnexin26 (hCx26) and connexin30 (hCx30), respectively, are the leading cause of non-syndromic prelingual deafness in several human populations. In this work, we exploited the high degree (77%) of sequence similarity shared by hCx26 and hCx30 to create atomistic models of homomeric hCx26 and hCx30 connexons starting from the X-ray crystallographic structure of an intercellular channel formed by hCx26 protomers at 3.5-å resolution. The equilibrium dynamics of the two protein complexes was followed for 40 ns each by Molecular Dynamics (MD) simulations. Our results indicate that, in hCx26, positively charged Lys41 residues establish a potential barrier within the fully open channel, hindering ion diffusion in the absence of an electrochemical gradient. A similar role is played, in hCx30, by negatively charged Glu49 residues. The different position and charge of these two ion sieves account for the differences in unitary conductance observed experimentally. Our results are discussed in terms of present models of voltage gating in connexin channels.