Source author record

Cristian Micheletti

Cristian Micheletti appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

RNA folding landscapes from explicit solvent all-atom simulations

Atomically detailed simulations of RNA folding have proven very challenging in view of the difficulties of developing realistic force fields and the intrinsic computational complexity of sampling rare conformational transitions. To tackle both these issues, we extend to RNA an enhanced path sampling method previously successfully applied to proteins. In this scheme, the information about the RNA's native structure is harnessed by a soft history-dependent biasing force, which is added to the atomistic force field, thus promoting the generation of productive folding trajectories. Here, we report on the results of simulations in explicit solvent of RNA molecules from 20 to 47 nucleotides long and increasing topological complexity. From a statistical analysis of the folding pathways we infer that, differently from proteins, the underlying free energy landscape is significantly frustrated, even for relatively small chains with a simple topology. The folding mechanism and the thermodynamics are in agreement with the available experiments and some of the existing coarse-grained models. This scheme provides a fully microscopic characterization of RNA folding, relating the kinetics and dynamics of the transition to the chemistry of the chain and its solvent. Therefore, it provides a transferable framework that sets the stage for future translational applications.

preprint2021arXiv

essHi-C: Essential component analysis of Hi-C matrices

Motivation: Hi-C matrices are cornerstones for qualitative and quantitative studies of genome folding, from its territorial organization to compartments and topological domains. The high dynamic range of genomic distances probed in Hi-C assays reflects in an inherent stochastic background of the interactions matrices, which inevitably convolve the features of interest with largely aspecific ones. Results: Here we introduce a discuss essHi-C, a method to isolate the specific, or essential component of Hi-C matrices from the aspecific portion of the spectrum that is compatible with random matrices. Systematic comparisons show that essHi-C improves the clarity of the interaction patterns, enhances the robustness against sequencing depth, allows the unsupervised clustering of experiments in different cell lines and recovers the cell-cycle phasing of single-cells based on Hi-C data. Thus, essHi-C provides means for isolating significant biological and physical features from Hi-C matrices.

preprint2020arXiv

Separation of Geometrical and Topological Entanglement in Confined Polymers Driven Out-Of-Equilibrium

We use Brownian dynamics simulations and advanced topological profiling methods to characterize the out-of-equilibrium evolution of self-entanglement in linear polymers confined into nano-channels and under periodic compression. We distinguish two main forms of entanglement, geometrical and topological. The latter is measured by the number of (essential) crossings of the physical knot detected after a suitable bridging of the chain termini. The former is instead measured as the average number of times a linear chain appears to cross itself when viewed under all projections, and is irrespective of the physical knotted state. The key discovery of our work is that these two forms of entanglement are uncoupled and evolve with distinct dynamics. While geometrical entanglement is typically in phase with the compression-elongation cycles and it is primarily sensitive to its force f, the topological measure is mildly sensitive to cyclic modulation but strongly depends on both compression force f and duration k. The findings could assist the interpretation of experiments using fluorescence molecular tracers to track physical knots in polymers. Furthermore, we identify optimal regions in the experimentally-controllable parameter space in which to obtain more/less topological and geometrical entanglement; this may help designing polymers with targeted topology.

preprint2016arXiv

Hi-C-constrained physical models of human chromosomes recover functionally-related properties of genome organization

Combining genome-wide structural models with phenomenological data is at the forefront of efforts to understand the organizational principles regulating the human genome. Here, we use chromosome-chromosome contact data as knowledge-based constraints for large-scale three-dimensional models of the human diploid genome. The resulting models remain minimally entangled and acquire several functional features that are observed in vivo and that were never used as input for the model. We find, for instance, that gene-rich, active regions are drawn towards the nuclear center, while gene poor and lamina-associated domains are pushed to the periphery. These and other properties persist upon adding local contact constraints, suggesting their compatibility with non-local constraints for the genome organization. The results show that suitable combinations of data analysis and physical modelling can expose the unexpectedly rich functionally-related properties implicit in chromosome-chromosome contact data. Specific directions are suggested for further developments based on combining experimental data analysis and genomic structural modelling.

preprint2016arXiv

How to fold intricately: using theory and experiments to unravel the properties of knotted proteins

Over the years, advances in experimental and computational methods have helped us to understand the role of thermodynamic, kinetic and active (chaperone-aided) effects in coordinating the folding steps required to achieving a knotted native state. Here, we review such developments by paying particular attention to the complementarity of experimental and computational studies. Key open issues that could be tackled with either or both approaches are finally pointed out.

preprint2015arXiv

Elastic network models for RNA: a comparative assessment with molecular dynamics and SHAPE experiments

Elastic network models (ENMs) are valuable and efficient tools for characterizing the collective internal dynamics of proteins based on the knowledge of their native structures. The increasing evidence that the biological functionality of RNAs is often linked to their innate internal motions, poses the question of whether ENM approaches can be successfully extended to this class of biomolecules. This issue is tackled here by considering various families of elastic networks of increasing complexity applied to a representative set of RNAs. The fluctuations predicted by the alternative ENMs are stringently validated by comparison against extensive molecular dynamics simulations and SHAPE experiments. We find that simulations and experimental data are systematically best reproduced by either an all-atom or a three-beads-per-nucleotide representation (sugar-base-phosphate), with the latter arguably providing the best balance of accuracy and computational complexity.

preprint2014arXiv

Synonymous mutations reduce genome compactness in icosahedral ssRNA viruses

Recent studies have shown that single-stranded viral RNAs fold into more compact structures than random RNA sequences with similar chemical composition and identical length. Based on this comparison it has been suggested that wild-type viral RNA may have evolved to be atypically compact so as to aid its encapsidation and assist the viral assembly process. In order to further explore the compactness selection hypothesis, we systematically compare the predicted sizes of more than one hundred wild-type viral sequences with those of their mutants, which are evolved in silico and subject to a number of known evolutionary constraints. In particular, we enforce mutation synonynimity, preserve the codon-bias, and leave untranslated regions intact. It is found that progressive accumulation of these restricted mutations still suffices to completely erase the characteristic compactness imprint of the viral RNA genomes, making them in this respect physically indistinguishable from randomly shuffled RNAs. This shows that maintaining the physical compactness of the genome is indeed a primary factor among ssRNA viruses evolutionary constraints, contributing also to the evidence that synonymous mutations in viral ssRNA genomes are not strictly neutral.

preprint2014arXiv

The unknotted strands of life: knots are very rare in RNA structures

The ongoing effort to detect and characterize physical entanglement in biopolymers has so far established that knots are present in many globular proteins and also abound in viral DNA packaged inside bacteriophages. RNA molecules, on the other hand, have not yet been systematically screened for the occurrence of physical knots. We have accordingly undertaken the systematic profiling of the ~6,000 RNA structures present in the protein data bank. The search identified no more than three deeply-knotted RNA molecules. These are ribosomal RNAs solved by cryo-em and consist of about 3,000 nucleotides. Compared to the case of proteins and viral DNA, the observed incidence of RNA knots is therefore practically negligible. This suggests that either evolutionary selection, or thermodynamic and kinetic folding mechanisms act towards minimizing the entanglement of RNA to an extent that is unparalleled by other types of biomolecules. The properties of the three observed RNA knotting patterns provide valuable clues for designing RNA sequences capable of self-tying in a twist-knot fold.

preprint2013arXiv

Colocalization of coregulated genes: a steered molecular dynamics study of human chromosome 19

The connection between chromatin nuclear organization and gene activity is vividly illustrated by the observation that transcriptional coregulation of certain genes appears to be directly influenced by their spatial proximity. This fact poses the more general question of whether it is at all feasible that the numerous genes that are coregulated on a given chromosome, especially those at large genomic distances, might become proximate inside the nucleus. This problem is studied here using steered molecular dynamics simulations in order to enforce the colocalization of thousands of knowledge-based gene sequences on a model for the gene-rich human chromosome 19. Remarkably, it is found that most (~80%) gene pairs can be brought simultaneously into contact. This is made possible by the low degree of intra-chromosome entanglement and the large number of cliques in the gene coregulatory network, that is the many groups of genes that are all mutually coregulated. The constrained conformations for the model chromosome 19 are further shown to be organised in spatial macrodomains that are similar to those inferred from recent HiC measurements. The findings indicate that gene coregulation and colocalization are largely compatible and that this relationship can be exploited to draft the overall spatial organization of the chromosome in vivo. The more general validity and implications of these findings could be investigated by applying to other eukaryotic chromosomes the general and transferable computational strategy introduced here.

preprint2013arXiv

Folding Pathways of a Knotted Protein with a Realistic Atomistic Force Field

We report on atomistic simulation of the folding of a natively-knotted protein, MJ0366, based on a realistic force field. To the best of our knowledge this is the first reported effort where a realistic force field is used to investigate the folding pathways of a protein with complex native topology. By using the dominant-reaction pathway scheme we collected about 30 successful folding trajectories for the 82-amino acid long trefoil-knotted protein. Despite the dissimilarity of their initial unfolded configuration, these trajectories reach the natively-knotted state through a remarkably similar succession of steps. In particular it is found that knotting occurs essentially through a threading mechanism, involving the passage of the C-terminal through an open region created by the formation of the native beta-sheet at an earlier stage. The dominance of the knotting by threading mechanism is not observed in MJ0366 folding simulations using simplified, native-centric models. This points to a previously underappreciated role of concerted amino acid interactions, including non-native ones, in aiding the appropriate order of contact formation to achieve knotting.

preprint2013arXiv

Spontaneous knotting and unknotting of flexible linear polymers: equilibrium and kinetic aspects

We report on a computational study of the statics and dynamics of long flexible linear polymers that spontaneously knot and unknot. Specifically, the equilibrium self-entanglement properties, such as the knotting probability, knot length and position, are investigated with extensive Monte Carlo sampling of chains of up to 15,000 beads. Tens of such equilibrated chains of up to 4, 096 beads are next used as starting points for Langevin dynamics simulations. The complex interplay of chain dynamics and self-knotting is addressed by monitoring the time evolution of various metric and entanglement properties. In particular, the extensive duration of the simulations allows for observing the spontaneous formation and disappearance of prime and composite physical knots in linear chains. Notably, a sizeable fraction of self-knotting and unknotting events is found to involve regions that are far away from the chain termini. To the best of our knowledge this represents the first instance where spontaneous changes in knotting for linear homopolymers are systematically characterized using unbiased dynamics simulations.

preprint2011arXiv

Multiscale entanglement in ring polymers under spherical confinement

The interplay of geometrical and topological entanglement in semiflexible knotted polymer rings confined inside a spherical cavity is investigated using advanced numerical methods. By using stringent and robust algorithms for locating knots, we characterize how the knot length lk depends on the ring contour length, Lc and the radius of the confining sphere, Rc . In the no- and strong- confinement cases we observe weak knot localization and complete knot delocalization, respectively. We show that the complex interplay of lk, Lc and Rc that seamlessly bridges these two limits can be encompassed by a simple scaling argument based on deflection theory. The same argument is used to rationalize the multiscale character of the entanglement that emerges with increasing confinement.

preprint2011arXiv

Polymers with spatial or topological constraints: theoretical and computational results

In this review we provide an organized summary of the theoretical and computational results which are available for polymers subject to spatial or topological constraints. Because of the interdisciplinary character of the topic, we provide an accessible, non-specialist introduction to the main topological concepts, polymer models, and theoretical/computational methods used to investigate dense and entangled polymer systems. The main body of our review deals with: (i) the effect that spatial confinement has on the equilibrium topological entanglement of one or more polymer chains and (ii) the metric and entropic properties of polymer chains with fixed topological state. These problems have important technological applications and implications for the life-sciences. Both aspects, especially the latter, are amply covered. A number of selected open problems are finally highlighted.

preprint2011arXiv

Probing the entanglement and locating knots in ring polymers: a comparative study of different arc closure schemes

The interplay between the topological and geometrical properties of a polymer ring can be clarified by establishing the entanglement trapped in any portion (arc) of the ring. The task requires to close the open arcs into a ring, and the resulting topological state may depend on the specific closure scheme that is followed. To understand the impact of this ambiguity in contexts of practical interest, such as knot localization in a ring with non trivial topology, we apply various closure schemes to model ring polymers. The rings have the same length and topological state (a trefoil knot) but have different degree of compactness. The comparison suggests that a novel method, termed the minimally-interfering closure, can be profitably used to characterize the arc entanglement in a robust and computationally-efficient way. This closure method is finally applied to the knot localization problem which is tackled using two different localization schemes based on top-down or bottom-up searches.

preprint2011arXiv

Structure and dynamics of ring polymers: entanglement effects because of solution density and ring topology

The effects of entanglement in solutions and melts of unknotted ring polymers have been addressed by several theoretical and numerical studies. The system properties have been typically profiled as a function of ring contour length at fixed solution density. Here, we use a different approach to investigate numerically the equilibrium and kinetic properties of solutions of model ring polymers. Specifically, the ring contour length is maintained fixed, while the interplay of inter- and intra-chain entanglement is modulated by varying both solution density (from infinite dilution up to \approx 40 % volume occupancy) and ring topology (by considering unknotted and trefoil-knotted chains). The equilibrium metric properties of rings with either topology are found to be only weakly affected by the increase of solution density. Even at the highest density, the average ring size, shape anisotropy and length of the knotted region differ at most by 40% from those of isolated rings. Conversely, kinetics are strongly affected by the degree of inter-chain entanglement: for both unknots and trefoils the characteristic times of ring size relaxation, reorientation and diffusion change by one order of magnitude across the considered range of concentrations. Yet, significant topology-dependent differences in kinetics are observed only for very dilute solutions (much below the ring overlap threshold). For knotted rings, the slowest kinetic process is found to correspond to the diffusion of the knotted region along the ring backbone.

preprint2006arXiv

Convergent dynamics in the protease enzymatic superfamily

Proteases regulate various aspects of the life cycle in all organisms by cleaving specific peptide bonds. Their action is so central for biochemical processes that at least 2% of any known genome encodes for proteolytic enzymes. Here we show that selected proteases pairs, despite differences in oligomeric state, catalytic residues and fold, share a common structural organization of functionally relevant regions which are further shown to undergo similar concerted movements. The structural and dynamical similarities found pervasively across evolutionarily distant clans point to common mechanisms for peptide hydrolysis.