Estimation of microscopic averages from metadynamics

With the help of metadynamics it is possible to calculate efficiently the free energy of systems displaying high energy barriers as a function of few selected "collective variables". In doing this, the contribution of all the other degrees of freedom ("microscopic" variables) is averaged out and, thus, lost. In the following, it is shown that it is possible to calculate the thermal average of these microscopic degrees of freedom during the metadynamics, not loosing this piece of information

Microbial Biofilm as a Smart Material

Microbial biofilm colonies will in many cases form a smart material capable of responding to external threats dependent on their size and internal state. The microbial community accordingly switches between passive, protective, or attack modes of action. In order to decide which strategy to employ, it is essential for the biofilm community to be able to sense its own size. The sensor designed to perform this task is termed a quorum sensor, since it only permits collective behaviour once a sufficiently large assembly of microbes have been established. The generic quorum sensor construct involves two genes, one coding for the production of a diffusible signal molecule and one coding for a regulator protein dedicated to sensing the signal molecules. A positive feedback in the signal molecule production sets a well-defined condition for switching into the collective mode. The activation of the regulator involves a slow dimerization, which allows low-pass filtering of the activation of the collective mode. Here, we review and combine the model components that form the basic quorum sensor in a number of Gram-negative bacteria, e.g., Pseudomonas aeruginosa

Compact phases of polymers with hydrogen bonding

We propose an off-lattice model for a self-avoiding homopolymer chain with two different competing attractive interactions, mimicking the hydrophobic effect and the hydrogen bond formation respectively. By means of Monte Carlo simulations, we are able to trace out the complete phase diagram for different values of the relative strength of the two competing interactions. For strong enough hydrogen bonding, the ground state is a helical conformation, whereas with decreasing hydrogen bonding strength, helices get eventually destabilized at low temperature in favor of more compact conformations resembling $\beta$-sheets appearing in native structures of proteins. For weaker hydrogen bonding helices are not thermodynamically relevant anymore.Comment: 5 pages, 3 figures; revised version published in PR

Information theoretical quantification of cooperativity in signalling complexes

<p>Abstract</p> <p>Background</p> <p>Intra-cellular information exchange, propelled by cascades of interacting signalling proteins, is essential for the proper functioning and survival of cells. Now that the interactome of several organisms is being mapped and several structural mechanisms of cooperativity at the molecular level in proteins have been elucidated, the formalization of this fundamental quantity, i.e. information, in these very diverse biological contexts becomes feasible.</p> <p>Results</p> <p>We show here that Shannon's mutual information quantifies information in biological system and more specifically the cooperativity inherent to the assembly of macromolecular complexes. We show how protein complexes can be considered as particular instances of noisy communication channels. Further we show, using a portion of the p27 regulatory pathway, how classical equilibrium thermodynamic quantities such as binding affinities and chemical potentials can be used to quantify information exchange but also to determine engineering properties such as channel noise and channel capacity. As such, this information measure identifies and quantifies those protein concentrations that render the biochemical system most effective in switching between the active and inactive state of the intracellular process.</p> <p>Conclusion</p> <p>The proposed framework provides a new and original approach to analyse the effects of cooperativity in the assembly of macromolecular complexes. It shows the conditions, provided by the protein concentrations, for which a particular system acts most effectively, i.e. exchanges the most information. As such this framework opens the possibility of grasping biological qualities such as system sensitivity, robustness or plasticity directly in terms of their effect on information exchange. Although these parameters might also be derived using classical thermodynamic parameters, a recasting of biological signalling in terms of information exchange offers an alternative framework for visualising network cooperativity that might in some cases be more intuitive.</p

Calculation of accurate small angle X-ray scattering curves from coarse-grained protein models

<p>Abstract</p> <p>Background</p> <p>Genome sequencing projects have expanded the gap between the amount of known protein sequences and structures. The limitations of current high resolution structure determination methods make it unlikely that this gap will disappear in the near future. Small angle X-ray scattering (SAXS) is an established low resolution method for routinely determining the structure of proteins in solution. The purpose of this study is to develop a method for the efficient calculation of accurate SAXS curves from coarse-grained protein models. Such a method can for example be used to construct a likelihood function, which is paramount for structure determination based on statistical inference.</p> <p>Results</p> <p>We present a method for the efficient calculation of accurate SAXS curves based on the Debye formula and a set of scattering form factors for dummy atom representations of amino acids. Such a method avoids the computationally costly iteration over all atoms. We estimated the form factors using generated data from a set of high quality protein structures. No <it>ad hoc </it>scaling or correction factors are applied in the calculation of the curves. Two coarse-grained representations of protein structure were investigated; two scattering bodies per amino acid led to significantly better results than a single scattering body.</p> <p>Conclusion</p> <p>We show that the obtained point estimates allow the calculation of accurate SAXS curves from coarse-grained protein models. The resulting curves are on par with the current state-of-the-art program CRYSOL, which requires full atomic detail. Our method was also comparable to CRYSOL in recognizing native structures among native-like decoys. As a proof-of-concept, we combined the coarse-grained Debye calculation with a previously described probabilistic model of protein structure, TorusDBN. This resulted in a significant improvement in the decoy recognition performance. In conclusion, the presented method shows great promise for use in statistical inference of protein structures from SAXS data.</p

Scale Free Cluster Distributions from Conserving Merging-Fragmentation Processes

We propose a dynamical scheme for the combined processes of fragmentation and merging as a model system for cluster dynamics in nature and society displaying scale invariant properties. The clusters merge and fragment with rates proportional to their sizes, conserving the total mass. The total number of clusters grows continuously but the full time-dependent distribution can be rescaled over at least 15 decades onto a universal curve which we derive analytically. This curve includes a scale free solution with a scaling exponent of -3/2 for the cluster sizes.Comment: 4 pages, 3 figure

Mean first passage time analysis reveals rate-limiting steps, parallel pathways and dead ends in a simple model of protein folding

We have analyzed dynamics on the complex free energy landscape of protein folding in the FOLD-X model, by calculating for each state of the system the mean first passage time to the folded state. The resulting kinetic map of the folding process shows that it proceeds in jumps between well-defined, local free energy minima. Closer analysis of the different local minima allows us to reveal secondary, parallel pathways as well as dead ends.Comment: 7 page

Quantifying information transfer by protein domains: Analysis of the Fyn SH2 domain structure

info:eu-repo/semantics/publishe

SNPeffect: a database mapping molecular phenotypic effects of human non-synonymous coding SNPs

Single nucleotide polymorphisms (SNPs) are an increasingly important tool for genetic and biomedical research. However, the accumulated sequence information on allelic variation is not matched by an understanding of the effect of SNPs on the functional attributes or ‘molecular phenotype’ of a protein. Towards this aim we developed SNPeffect, an online resource of human non-synonymous coding SNPs (nsSNPs) mapping phenotypic effects of allelic variation in human genes. SNPeffect contains 31 659 nsSNPs from 12 480 human proteins. The current release of SNPeffect incorporates data on protein stability, integrity of functional sites, protein phosphorylation and glycosylation, subcellular localization, protein turnover rates, protein aggregation, amyloidosis and chaperone interaction. The SNP entries are accessible through both a search and browse interface and are linked to most major biological databases. The data can be displayed as detailed descriptions of individual SNPs or as an overview of all SNPs for a given protein. SNPeffect will be regularly updated and can be accessed at http://snpeffect.vib.be/