Folding, Design and Determination of Interaction Potentials Using Off-Lattice Dynamics of Model Heteropolymers

We present the results of a self-consistent, unified molecular dynamics study of simple model heteropolymers in the continuum with emphasis on folding, sequence design and the determination of the interaction parameters of the effective potential between the amino acids from the knowledge of the native states of the designed sequences.Comment: 8 pages, 3 Postscript figures, uses RevTeX. Submitted to Physical Review Letter

Geometry-controlled kinetics

It has long been appreciated that transport properties can control reaction kinetics. This effect can be characterized by the time it takes a diffusing molecule to reach a target -- the first-passage time (FPT). Although essential to quantify the kinetics of reactions on all time scales, determining the FPT distribution was deemed so far intractable. Here, we calculate analytically this FPT distribution and show that transport processes as various as regular diffusion, anomalous diffusion, diffusion in disordered media and in fractals fall into the same universality classes. Beyond this theoretical aspect, this result changes the views on standard reaction kinetics. More precisely, we argue that geometry can become a key parameter so far ignored in this context, and introduce the concept of "geometry-controlled kinetics". These findings could help understand the crucial role of spatial organization of genes in transcription kinetics, and more generally the impact of geometry on diffusion-limited reactions.Comment: Submitted versio

Working Group on Biological Parameters (WGBIOP) 2021

The main objective of the Working Group on Biological Parameters (WGBIOP) is to review the status, issues, developments, and quality assurance of biological parameters used in assessment and management. WGBIOP (1) plans workshops, exchanges, and validation studies on a range of biological varia-bles to review the quality of information supplied for stock assessment and improve quality as-surance and training; (2) investigates data availability and develops documentation and methods to improve communication between data collectors and end-users; (3) delivers new and im-proved functionality for the SmartDots platform. Four otolith exchanges and two workshops were completed in 2020–2021 using SmartDots— eight further exchanges are ongoing. Proposed future exchanges and workshops were reviewed and approved. The development of the SmartDots platform proceeded with the inclusion of the maturity, eggs, atresia, fecundity, and larval identification modules into the software version. A live SmartDots tutorial for event coordinators was conducted. Work to further develop quality assurance guidelines—and review national applications of these—progressed. Age and maturity validation studies were reviewed and a new method for prioritizing future validation work was proposed. Progress with the Stock Identification Database (SID) was reviewed, and the potential for creating a WGBIOP library collection and active involvement of WGBIOP in updating FishBase.org data were evaluated. The importance of identifying and documenting links be-tween all relevant databases and document repositories was identified, and a task to address this was initiated. Work on improving the feedback loop between data collectors and stock assessors on the usage and quality of biological parameters in stock assessment continued. Moving forward, WGBIOP aims to continue collaboration with WGALES and WGSMART on the development of the SmartDots platform, encouraging cross-group sharing of skills and ex-perience to optimize results. WGBIOP aims to improve accessibility to its outputs through up-dates to SID and FishBase.org, and the potential creation of a WGBIOP library collection. WGBIOP hopes to improve two-way communication between data collectors and end-users around the quality and utility of biological parameters used in assessment. WGBIOP also aims to amalgamate all validation activities into one coherent workstream.ICE

The Influence of Transcription Factor Competition on the Relationship between Occupancy and Affinity

Transcription factors (TFs) are proteins that bind to specific sites on the DNA and regulate gene activity. Identifying where TF molecules bind and how much time they spend on their target sites is key to understanding transcriptional regulation. It is usually assumed that the free energy of binding of a TF to the DNA (the affinity of the site) is highly correlated to the amount of time the TF remains bound (the occupancy of the site). However, knowing the binding energy is not sufficient to infer actual binding site occupancy. This mismatch between the occupancy predicted by the affinity and the observed occupancy may be caused by various factors, such as TF abundance, competition between TFs or the arrangement of the sites on the DNA. We investigated the relationship between the affinity of a TF for a set of binding sites and their occupancy. In particular, we considered the case of the transcription factor lac repressor (lacI) in E.coli, and performed stochastic simulations of the TF dynamics on the DNA for various combinations of lacI abundance and competing TFs that contribute to macromolecular crowding. We also investigated the relationship of site occupancy and the information content of position weight matrices (PWMs) used to represent binding sites. Our results showed that for medium and high affinity sites, TF competition does not play a significant role for genomic occupancy except in cases when the abundance of the TF is significantly increased, or when the PWM displays relatively low information content. Nevertheless, for medium and low affinity sites, an increase in TF abundance (for both cognate and non-cognate molecules) leads to an increase in occupancy at several sites. © 2013 Zabet et al

Expression-Dependent Folding of Interphase Chromatin

Multiple studies suggest that chromatin looping might play a crucial role in organizing eukaryotic genomes. To investigate the interplay between the conformation of interphase chromatin and its transcriptional activity, we include information from gene expression profiles into a polymer model for chromatin that incorporates genomic loops. By relating loop formation to transcriptional activity, we are able to generate chromosome conformations whose structural and topological properties are consistent with experimental data. The model particularly allows to reproduce the conformational variations that are known to occur between highly and lowly expressed chromatin regions. As previously observed in experiments, lowly expressed regions of the simulated polymers are much more compact. Due to the changes in loop formation, the distributions of chromatin loops are also expression-dependent and exhibit a steeper decay in highly active regions. As a results of entropic interaction between differently looped parts of the chromosome, we observe topological alterations leading to a preferential positioning of highly transcribed loci closer to the surface of the chromosome territory. Considering the diffusional behavior of the chromatin fibre, the simulations furthermore show that the higher the expression level of specific parts of the chromatin fibre is, the more dynamic they are. The results exhibit that variations of loop formation along the chromatin fibre, and the entropic changes that come along with it, do not only influence the structural parameters on the local scale, but also effect the global chromosome conformation and topology

SMC complexes differentially compact mitotic chromosomes according to genomic context

Structural maintenance of chromosomes (SMC) protein complexes are key determinants of chromosome conformation. Using Hi-C and polymer modelling, we study how cohesin and condensin, two deeply conserved SMC complexes, organize chromosomes in the budding yeast Saccharomyces cerevisiae. The canonical role of cohesin is to co-align sister chromatids, while condensin generally compacts mitotic chromosomes. We find strikingly different roles for the two complexes in budding yeast mitosis. First, cohesin is responsible for compacting mitotic chromosome arms, independently of sister chromatid cohesion. Polymer simulations demonstrate that this role can be fully accounted for through cis-looping of chromatin. Second, condensin is generally dispensable for compaction along chromosome arms. Instead, it plays a targeted role compacting the rDNA proximal regions and promoting resolution of peri-centromeric regions. Our results argue that the conserved mechanism of SMC complexes is to form chromatin loops and that distinct SMC-dependent looping activities are selectively deployed to appropriately compact chromosomes

An Estimate of the Numbers and Density of Low-Energy Structures (or Decoys) in the Conformational Landscape of Proteins

The conformational energy landscape of a protein, as calculated by known potential energy functions, has several minima, and one of these corresponds to its native structure. It is however difficult to comprehensively estimate the actual numbers of low energy structures (or decoys), the relationships between them, and how the numbers scale with the size of the protein.We have developed an algorithm to rapidly and efficiently identify the low energy conformers of oligo peptides by using mutually orthogonal Latin squares to sample the potential energy hyper surface. Using this algorithm, and the ECEPP/3 potential function, we have made an exhaustive enumeration of the low-energy structures of peptides of different lengths, and have extrapolated these results to larger polypeptides.We show that the number of native-like structures for a polypeptide is, in general, an exponential function of its sequence length. The density of these structures in conformational space remains more or less constant and all the increase appears to come from an expansion in the volume of the space. These results are consistent with earlier reports that were based on other models and techniques

Structure and Age Jointly Influence Rates of Protein Evolution

What factors determine a protein's rate of evolution are actively debated. Especially unclear is the relative role of intrinsic factors of present-day proteins versus historical factors such as protein age. Here we study the interplay of structural properties and evolutionary age, as determinants of protein evolutionary rate. We use a large set of one-to-one orthologs between human and mouse proteins, with mapped PDB structures. We report that previously observed structural correlations also hold within each age group – including relationships between solvent accessibility, designabililty, and evolutionary rates. However, age also plays a crucial role: age modulates the relationship between solvent accessibility and rate. Additionally, younger proteins, despite being less designable, tend to evolve faster than older proteins. We show that previously reported relationships between age and rate cannot be explained by structural biases among age groups. Finally, we introduce a knowledge-based potential function to study the stability of proteins through large-scale computation. We find that older proteins are more stable for their native structure, and more robust to mutations, than younger ones. Our results underscore that several determinants, both intrinsic and historical, can interact to determine rates of protein evolution