Structurally constrained protein evolution: results from a lattice simulation

We simulate the evolution of a protein-like sequence subject to point mutations, imposing conservation of the ground state, thermodynamic stability and fast folding. Our model is aimed at describing neutral evolution of natural proteins. We use a cubic lattice model of the protein structure and test the neutrality conditions by extensive Monte Carlo simulations. We observe that sequence space is traversed by neutral networks, i.e. sets of sequences with the same fold connected by point mutations. Typical pairs of sequences on a neutral network are nearly as different as randomly chosen sequences. The fraction of neutral neighbors has strong sequence to sequence variations, which influence the rate of neutral evolution. In this paper we study the thermodynamic stability of different protein sequences. We relate the high variability of the fraction of neutral mutations to the complex energy landscape within a neutral network, arguing that valleys in this landscape are associated to high values of the neutral mutation rate. We find that when a point mutation produces a sequence with a new ground state, this is likely to have a low stability. Thus we tentatively conjecture that neutral networks of different structures are typically well separated in sequence space. This results indicates that changing significantly a protein structure through a biologically acceptable chain of point mutations is a rare, although possible, event.Comment: added reference, to appear on European Physical Journal

Statistical properties of neutral evolution

Neutral evolution is the simplest model of molecular evolution and thus it is most amenable to a comprehensive theoretical investigation. In this paper, we characterize the statistical properties of neutral evolution of proteins under the requirement that the native state remains thermodynamically stable, and compare them to the ones of Kimura's model of neutral evolution. Our study is based on the Structurally Constrained Neutral (SCN) model which we recently proposed. We show that, in the SCN model, the substitution rate decreases as longer time intervals are considered, and fluctuates strongly from one branch of the evolutionary tree to another, leading to a non-Poissonian statistics for the substitution process. Such strong fluctuations are also due to the fact that neutral substitution rates for individual residues are strongly correlated for most residue pairs. Interestingly, structurally conserved residues, characterized by a much below average substitution rate, are also much less correlated to other residues and evolve in a much more regular way. Our results could improve methods aimed at distinguishing between neutral and adaptive substitutions as well as methods for computing the expected number of substitutions occurred since the divergence of two protein sequences.Comment: 17 pages, 11 figure

Spreading of infections on random graphs: A percolation-type model for COVID-19

We introduce an epidemic spreading model on a network using concepts from percolation theory. The model is motivated by discussing the standard SIR model, with extensions to describe effects of lockdowns within a population. The underlying ideas and behavior of the lattice model, implemented using the same lockdown scheme as for the SIR scheme, are discussed in detail and illustrated with extensive simulations. A comparison between both models is presented for the case of COVID-19 data from the USA. Both fits to the empirical data are very good, but some differences emerge between the two approaches which indicate the usefulness of having an alternative approach to the widespread SIR model

Scaling model for a speed-dependent vehicle noise spectrum

Abstract Considering the well-known features of the noise emitted by moving sources, a number of vehicle characteristics such as speed, unladen mass, engine size, year of registration, power and fuel were recorded in a dedicated monitoring campaign performed in three different places, each characterized by different number of lanes and the presence of nearby reflective surfaces. A full database of 144 vehicles (cars) was used to identify statistically relevant features. In order to compare the vehicle transit noise in different environmental condition, all 1/3-octave band spectra were normalized and analysed. Unsupervised clustering algorithms were employed to group together spectrum levels with similar profiles. Our results corroborate the well-known fact that speed is the most relevant characteristic to discriminate between different vehicle noise spectrum. In keeping with this fact, we present a new approach to predict analytically noise spectra for a given vehicle speed. A set of speed-dependent analytical functions are suggested in order to fit the normalized average spectrum profile at different speeds. This approach can be useful for predicting vehicle speed based purely on its noise spectrum pattern. The present work is complementary to the accurate analysis of noise sources based on the beamforming technique

Nanocrystalline versus microcrystalline Lo2O:B2O 3 composites: Anomalous ionic conductivities and percolation theory

We study ionic transport in nano- and microcrystalline (1−x)Li2O:xB2O3 composites using standard impedance spectroscopy. In the nanocrystalline samples (average grain size of about 20 nm), the ionic conductivity σdc increases with increasing content x of B2O3 up to a maximum at x≈0.5. Above x≈0.92, σdc vanishes. By contrast, in the microcrystalline samples (grain size about 10μm), σdc decreases monotonically with x and vanishes above x≈0.55. We can explain this strikingly different behavior by a percolation model that assumes an enhanced conductivity at the interfaces between insulating and conducting phases in both materials and explicitly takes into account the different grain sizes. © 2000 The American Physical Society

A protein evolution model with independent sites that reproduces site-specific amino acid distributions from the Protein Data Bank.

BACKGROUND: Since thermodynamic stability is a global property of proteins that has to be conserved during evolution, the selective pressure at a given site of a protein sequence depends on the amino acids present at other sites. However, models of molecular evolution that aim at reconstructing the evolutionary history of macromolecules become computationally intractable if such correlations between sites are explicitly taken into account. RESULTS: We introduce an evolutionary model with sites evolving independently under a global constraint on the conservation of structural stability. This model consists of a selection process, which depends on two hydrophobicity parameters that can be computed from protein sequences without any fit, and a mutation process for which we consider various models. It reproduces quantitatively the results of Structurally Constrained Neutral (SCN) simulations of protein evolution in which the stability of the native state is explicitly computed and conserved. We then compare the predicted site-specific amino acid distributions with those sampled from the Protein Data Bank (PDB). The parameters of the mutation model, whose number varies between zero and five, are fitted from the data. The mean correlation coefficient between predicted and observed site-specific amino acid distributions is larger than = 0.70 for a mutation model with no free parameters and no genetic code. In contrast, considering only the mutation process with no selection yields a mean correlation coefficient of = 0.56 with three fitted parameters. The mutation model that best fits the data takes into account increased mutation rate at CpG dinucleotides, yielding = 0.90 with five parameters. CONCLUSION: The effective selection process that we propose reproduces well amino acid distributions as observed in the protein sequences in the PDB. Its simplicity makes it very promising for likelihood calculations in phylogenetic studies. Interestingly, in this approach the mutation process influences the effective selection process, i.e. selection and mutation must be entangled in order to obtain effectively independent sites. This interdependence between mutation and selection reflects the deep influence that mutation has on the evolutionary process: The bias in the mutation influences the thermodynamic properties of the evolving proteins, in agreement with comparative studies of bacterial proteomes, and it also influences the rate of accepted mutations.RIGHTS : This article is licensed under the BioMed Central licence at http://www.biomedcentral.com/about/license which is similar to the 'Creative Commons Attribution Licence'. In brief you may : copy, distribute, and display the work; make derivative works; or make commercial use of the work - under the following conditions: the original author must be given credit; for any reuse or distribution, it must be made clear to others what the license terms of this work are