Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulating topological domains in human chromosomes with a fitting-free model
68
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We discuss a polymer model for the 3D organization of human chromosomes. A chromosome is represented by a string of beads, with each bead being colored according to 1D bioinformatic data (e.g., chromatin state, histone modification, GC content). Individual spheres (representing bi- and multi-valent transcription factors) can bind reversibly and selectively to beads with the appropriate color. During molecular dynamics simulations, the factors bind, and the string spontaneously folds into loops, rosettes, and topologically-associating domains (TADs). This organization occurs in the absence of any specified interactions between distant DNA segments, or between transcription factors. A comparison with Hi-C data shows that simulations predict the location of most boundaries between TADs correctly. The model is fitting-free in the sense that it does not use Hi-C data as an input; consequently, one of its strengths is that it can -- in principle -- be used to predict the 3D organization of any region of interest, or whole chromosome, in a given organism, or cell line, in the absence of existing Hi-C data. We discuss how this simple model might be refined to include more transcription factors and binding sites, and to correctly predict contacts between convergent CTCF binding sites.
rate research
Read More
The probability of two loci, separated by a certain genome length, being in contact can be inferred using the Chromosome Conformation Capture (3C) method and related Hi-C experiments. How to go from the contact map, a matrix listing the mean contact probabilities between a large number of pairs of loci, to an ensemble of three-dimensional structures is an open problem. A solution to this problem, without assuming an assumed energy function, would be the first step in understanding the way nature has solved the packaging of chromosomes in tight cellular spaces. We created a theory, based on polymer physics characteristics of chromosomes and the maximum entropy principles, referred to as HIPPS (Hi-C-Polymer-Physics-Structures) method, that allows us to calculate the 3D structures solely from Hi-C contact maps. We created an ensemble of 3D structures for the 23 chromosomes from lymphoblastoid cells using the measured contact maps as inputs. The HIPPS method shows that conformations of chromosomes are heterogeneous even in a single cell type. The differences in the conformational heterogeneity of the same chromosome in different cell types (normal as well as cancerous cells) can also be quantitatively discerned using our theory. We validate the method by showing that the calculated volumes of the 23 chromosomes from the predicted 3D structures are in good agreement with experimental estimates. Because the method is general, the 3D structures for any species may be calculated directly from the contact map without the need to assume a specific polymer model, as is customarily done.
We propose a stochastic model for gene transcription coupled to DNA supercoiling, where we incorporate the experimental observation that polymerases create supercoiling as they unwind the DNA helix, and that these enzymes bind more favourably to regions where the genome is unwound. Within this model, we show that when the transcriptionally induced flux of supercoiling increases, there is a sharp crossover from a regime where torsional stresses relax quickly and gene transcription is random, to one where gene expression is highly correlated and tightly regulated by supercoiling. In the latter regime, the model displays transcriptional bursts, waves of supercoiling, and up-regulation of divergent or bidirectional genes. It also predicts that topological enzymes which relax twist and writhe should provide a pathway to down-regulate transcription. This article has been published in Physical Review Letters, May 2016.
In multi-resolution simulations, different system components are simultaneously modelled at different levels of resolution, these being smoothly coupled together. In the case of enzyme systems, computationally expensive atomistic detail is needed in the active site to capture the chemistry of substrate binding. Global properties of the rest of the protein also play an essential role, determining the structure and fluctuations of the binding site; however, these can be modelled on a coarser level. Similarly, in the most computationally efficient scheme only the solvent hydrating the active site requires atomistic detail. We present a methodology to couple atomistic and coarse-grained protein models, while solvating the atomistic part of the protein in atomistic water. This allows a free choice of which protein and solvent degrees of freedom to include atomistically, without loss of accuracy in the atomistic description. This multi-resolution methodology can successfully model stable ligand binding, and we further confirm its validity via an exploration of system properties relevant to enzymatic function. In addition to a computational speedup, such an approach can allow the identification of the essential degrees of freedom playing a role in a given process, potentially yielding new insights into biomolecular function.
Genomes contain rare guanine-rich sequences capable of assembling into four-stranded helical structures, termed G-quadruplexes, with potential roles in gene regulation and chromosome stability. Their mechanical unfolding has only been reported to date by all-atom simulations, which cannot dissect the major physical interactions responsible for their cohesion. Here, we propose a mesoscopic model to describe both the mechanical and thermal stability of DNA G-quadruplexes, where each nucleotide of the structure, as well as each central cation located at the inner channel, is mapped onto a single bead. In this framework we are able to simulate loading rates similar to the experimental ones, which are not reachable in simulations with atomistic resolution. In this regard, we present single-molecule force-induced unfolding experiments by a high-resolution optical tweezers on a DNA telomeric sequence capable of forming a G-quadruplex conformation. Fitting the parameters of the model to the experiments we find a correct prediction of the rupture-force kinetics and a good agreement with previous near equilibrium measurements. Since G-quadruplex unfolding dynamics is halfway in complexity between secondary nucleic acids and tertiary protein structures, our model entails a nanoscale paradigm for non-equilibrium processes in the cell.
The elastic network (EN) is a prime model that describes the long-time dynamics of biomolecules. However, the use of harmonic potentials renders this model insufficient for studying large conformational changes of proteins (e.g. stretching of proteins, folding and thermal unfolding). Here, we extend the capabilities of the EN model by using a harmonic approximation described by Lennard-Jones (LJ) interactions for far contacts and native contacts obtained from the standard overlap criterion as in the case of Go-like models. While our model is validated against the EN model by reproducing the equilibrium properties for a number of proteins, we also show that the model is suitable for the study of large conformation changes by providing various examples. In particular, this is illustrated on the basis of pulling simulations that predict with high accuracy the experimental data on the rupture force of the studied proteins. Furthermore, in the case of DDFLN4 protein, our pulling simulations highlight the advantages of our model with respect to Go-like approaches, where the latter fail to reproduce previous results obtained by all-atom simulations that predict an additional characteristic peak for this protein. In addition, folding simulations of small peptides yield different folding times for alpha-helix and beta-hairpin, in agreement with experiment, in this way providing further opportunities for the application of our model in studying large conformational changes of proteins. In contrast to the EN model, our model is suitable for both normal mode analysis and molecular dynamics simulation. We anticipate that the proposed model will find applications in a broad range of problems in biology, including, among others, protein folding and thermal unfolding.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
