Stochastic dynamics and electron-conformation interactions in proteins

A review of modern concepts of conformational motions in proteins is presented. Problems relevant to the theory of the dynamic behavior of different structural elements of biomacromolecules and a protein globule as a whole are discussed. A theory of electron-conformation interactions is presented and its applications to electron-transport reactions, the enzymic act, and the diffusion of the ligands in the proteins are discussed.     

Stochastic tunnels in evolutionary dynamics

We study a situation that arises in the somatic evolution of cancer. Consider a finite population of replicating cells and a sequence of mutations: type 0 can mutate to type 1, which can mutate to type 2. There is no back mutation. We start with a homogeneous population of type 0. Mutants of type 1 emerge and either become extinct or reach fixation. In both cases, they can generate type 2, which also can become extinct or reach fixation. If mutation rates are small compared to the inverse of the population size, then the stochastic dynamics can be described by transitions between homogeneous populations. A "stochastic tunnel" arises, when the population moves from all 0 to all 2 without ever being all 1. We calculate the exact rate of stochastic tunneling for the case when type 1 is as fit as type 0 or less fit. Type 2 has the highest fitness. We discuss implications for the elimination of tumor suppressor genes and the activation of genetic instability. Although our theory is developed for cancer genetics, stochastic tunnels are general phenomena that could arise in many circumstances.     

Stochastic dynamics of cancer initiation

Most human cancer types result from the accumulation of multiple genetic and epigenetic alterations in a single cell. Once the first change (or changes) have arisen, tumorigenesis is initiated and the subsequent emergence of additional alterations drives progression to more aggressive and ultimately invasive phenotypes. Elucidation of the dynamics of cancer initiation is of importance for an understanding of tumor evolution and cancer incidence data. In this paper, we develop a novel mathematical framework to study the processes of cancer initiation. Cells at risk of accumulating oncogenic mutations are organized into small compartments of cells and proliferate according to a stochastic process. During each cell division, an (epi)genetic alteration may arise which leads to a random fitness change, drawn from a probability distribution. Cancer is initiated when a cell gains a fitness sufficiently high to escape from the homeostatic mechanisms of the cell compartment. To investigate cancer initiation during a human lifetime, a 'race' between this fitness process and the aging process of the patient is considered; the latter is modeled as a second stochastic Markov process in an aging dimension. This model allows us to investigate the dynamics of cancer initiation and its dependence on the mutational fitness distribution. Our framework also provides a methodology to assess the effects of different life expectancy distributions on lifetime cancer incidence. We apply this methodology to colorectal tumorigenesis while considering life expectancy data of the US population to inform the dynamics of the aging process. We study how the probability of cancer initiation prior to death, the time until cancer initiation, and the mutational profile of the cancer-initiating cell depends on the shape of the mutational fitness distribution and life expectancy of the population.     

Stochastic models of neuronal dynamics

Cortical activity is the product of interactions among neuronal populations. Macroscopic electrophysiological phenomena are generated by these interactions. In principle, the mechanisms of these interactions afford constraints on biologically plausible models of electrophysiological responses. In other words, the macroscopic features of cortical activity can be modelled in terms of the microscopic behaviour of neurons. An evoked response potential (ERP) is the mean electrical potential measured from an electrode on the scalp, in response to some event. The purpose of this paper is to outline a population density approach to modelling ERPs.We propose a biologically plausible model of neuronal activity that enables the estimation of physiologically meaningful parameters from electrophysiological data. The model encompasses four basic characteristics of neuronal activity and organization: (i) neurons are dynamic units, (ii) driven by stochastic forces, (iii) organized into populations with similar biophysical properties and response characteristics and (iv) multiple populations interact to form functional networks. This leads to a formulation of population dynamics in terms of the Fokker-Planck equation. The solution of this equation is the temporal evolution of a probability density over state-space, representing the distribution of an ensemble of trajectories. Each trajectory corresponds to the changing state of a neuron. Measurements can be modelled by taking expectations over this density, e.g. mean membrane potential, firing rate or energy consumption per neuron. The key motivation behind our approach is that ERPs represent an average response over many neurons. This means it is sufficient to model the probability density over neurons, because this implicitly models their average state. Although the dynamics of each neuron can be highly stochastic, the dynamics of the density is not. This means we can use Bayesian inference and estimation tools that have already been established for deterministic systems. The potential importance of modelling density dynamics (as opposed to more conventional neural mass models) is that they include interactions among the moments of neuronal states (e.g. the mean depolarization may depend on the variance of synaptic currents through nonlinear mechanisms).Here, we formulate a population model, based on biologically informed model-neurons with spike-rate adaptation and synaptic dynamics. Neuronal sub-populations are coupled to form an observation model, with the aim of estimating and making inferences about coupling among sub-populations using real data. We approximate the time-dependent solution of the system using a bi-orthogonal set and first-order perturbation expansion. For didactic purposes, the model is developed first in the context of deterministic input, and then extended to include stochastic effects. The approach is demonstrated using synthetic data, where model parameters are identified using a Bayesian estimation scheme we have described previously.     

Stochastic dynamics of metastasis formation

Tumor metastasis accounts for the majority of deaths in cancer patients. The metastatic behavior of cancer cells is promoted by mutations in many genes, including activation of oncogenes such as RAS and MYC. Here, we develop a mathematical framework to analyse the dynamics of mutations enabling cells to metastasize. We consider situations in which one mutation is necessary to confer metastatic ability to the cell. We study different population sizes of the main tumor and different somatic fitness values of metastatic cells. We compare mutations that are positively selected in the main tumor with those that are neutral or negatively selected, but faster at forming metastases. We study whether metastatic potential is the property of all (or the majority of) cells in the main tumor or only the property of a small subset. Our theory shows how to calculate the expected number of metastases that are formed by a tumor.     

Vicinal coupling constants and protein dynamics

The effects of motional averaging on the analysis of vicinal spin-spin coupling constants derived from proton NMR studies of proteins have been examined. Trajectories obtained from molecular dynamics simulations of bovine pancreatic trypsin inhibitor and of hen egg white lysozyme were used in conjunction with an expression for the dependence of the coupling constant on the intervening dihedral angle to calculate the time-dependent behavior of the coupling constants. Despite large fluctuations, the time-average values of the coupling constants are not far from those computed for the average structure in the cases where fluctuations occur about a single potential well. The calculated differences show a high correlation with the variation in the magnitude of the fluctuations of individual dihedral angles. For the cases where fluctuations involve multiple sites, large differences are found between the time-average values and the average structure values for the coupling constants. Comparison of the simulation results with the experimental trends suggests that side chains with more than one position are more common in proteins than is inferred from X-ray results. It is concluded that for the main chain, motional effects do not introduce significant errors where vicinal coupling constants are used in structure determinations; however, for side chains, the motional average can alter deductions about the structure. Accurately measured coupling constants are shown to provide information concerning the magnitude of dihedral angle fluctuations.     

Internal dynamics and protein-matrix coupling in trehalose-coated proteins

We review recent studies on the role played by non-liquid, water-containing matrices on the dynamics and structure of embedded proteins. Two proteins were studied, in water-trehalose matrices: a water-soluble protein (carboxy derivative of horse heart myoglobin) and a membrane protein (reaction centre from Rhodobacter sphaeroides). Several experimental techniques were used: Mossbauer spectroscopy, elastic neutron scattering, FTIR spectroscopy, CO recombination after flash photolysis in carboxy-myoglobin, kinetic optical absorption spectroscopy following pulsed and continuous photoexcitation in Q(B) containing or Q(B) deprived reaction centre from R. sphaeroides. Experimental results, together with the outcome of molecular dynamics simulations, concurred to give a picture of how water-containing matrices control the internal dynamics of the embedded proteins. This occurs, in particular, via the formation of hydrogen bond networks that anchor the protein surface to the surrounding matrix, whose stiffness increases by lowering the sample water content. In the conclusion section, we also briefly speculate on how the protein-matrix interactions observed in our samples may shed light on the protein-solvent coupling also in liquid aqueous solutions.     

Sequence of interactions in receptor-G protein coupling

Guanine nucleotide exchange in heterotrimeric G proteins catalyzed by G protein-coupled receptors (GPCRs) is a key event in many physiological processes. The crystal structures of the GPCR rhodopsin and two G proteins as well as binding sites on both catalytically interacting proteins are known, but the temporal sequence of events leading to nucleotide exchange remains to be elucidated. We employed time-resolved near infrared light scattering to study the order in which the Galpha and Ggamma C-terminal binding sites on the holo-G protein interact with the active state of the GPCR rhodopsin (R*) in native membranes. We investigated these key binding sites within mass-tagged peptides and G proteins and found that their binding to R* is mutually exclusive. The interaction of the holo-G protein with R* requires at least one of the lipid modifications of the G protein (i.e. myristoylation of the Galpha N terminus and/or farnesylation of the Ggamma C terminus). A holo-G protein with a high affinity Galpha C terminus shows a specific change of the reaction rate in the GDP release and GTP uptake steps of catalysis. We interpret the data by a sequential fit model where (i) the initial encounter between R* and the G protein occurs with the Gbetagamma subunit, and (ii) the Galpha C-terminal tail then interacts with R* to release bound GDP, thereby decreasing the affinity of R* for the Gbetagamma subunit. The mechanism limits the time in which both C-terminal binding sites of the G protein interact simultaneously with R* to a short lived transitory state.     

Stochastic evolutionary dynamics on two levels

We consider a population that is subdivided into groups. Individuals reproduce proportional to their fitness. When a group reaches a certain size it has a probability to split into two groups while another group is eliminated. In this stochastic process, the number of groups is constant, while the total population size fluctuates between well-defined bounds. We calculate the fixation probability of newly introduced mutants under constant selection. We show that the described population structure acts as a suppressor of selection compared to an unstructured population of the same size. The maximum suppression of selection is obtained, when the number of groups equals the number of individuals per group. We also study opposing selection on two or more levels by analysing the evolutionary dynamics of hierarchically embedded Moran processes.     

Stochastic evolutionary dynamics of bimatrix games

Evolutionary game dynamics of two-player asymmetric games in finite populations is studied. We consider two roles in the game, roles α and β. α-players and β-players interact and gain payoffs. The game is described by a pair of matrices, which is called bimatrix. One's payoff in the game is interpreted as its fecundity, thus strategies are subject to natural selection. In addition, strategies can randomly mutate to others. We formulate a stochastic evolutionary game dynamics of bimatrix games as a frequency-dependent Moran process with mutation. We analytically derive the stationary distribution of strategies under weak selection. Our result provides a criterion for equilibrium selection in general bimatrix games.     

Stochastic evolutionary dynamics of direct reciprocity

Evolutionary game theory is the study of frequency-dependent selection. The success of an individual depends on the frequencies of strategies that are used in the population. We propose a new model for studying evolutionary dynamics in games with a continuous strategy space. The population size is finite. All members of the population use the same strategy. A mutant strategy is chosen from some distribution over the strategy space. The fixation probability of the mutant strategy in the resident population is calculated. The new mutant takes over the population with this probability. In this case, the mutant becomes the new resident. Otherwise, the existing resident remains. Then, another mutant is generated. These dynamics lead to a stationary distribution over the entire strategy space. Our new approach generalizes classical adaptive dynamics in three ways: (i) the population size is finite; (ii) mutants can be drawn non-locally and (iii) the dynamics are stochastic. We explore reactive strategies in the repeated Prisoner''s Dilemma. We perform ‘knock-out experiments’ to study how various strategies affect the evolution of cooperation. We find that ‘tit-for-tat’ is a weak catalyst for the emergence of cooperation, while ‘always cooperate’ is a strong catalyst for the emergence of defection. Our analysis leads to a new understanding of the optimal level of forgiveness that is needed for the evolution of cooperation under direct reciprocity.     

Stochastic models for phytoplankton dynamics in Mediterranean Sea

In this paper, we review some results obtained from three one-dimensional stochastic models, which were used to analyze picophytoplankton dynamics in two sites of the Mediterranean Sea. Firstly, we present a stochastic advection–reaction–diffusion model to describe the vertical spatial distribution of picoeukaryotes in a site of the Sicily Channel. The second model, which is an extended version of the first one, is used to obtain the vertical stationary profiles of two groups of picophytoplankton, i.e. Pelagophytes and Prochlorococcus, in the same marine site as in the previous case. Here, we include intraspecific competition of picophytoplanktonic groups for limiting factors, i.e. light intensity and nutrient concentration. Finally, we analyze the spatio-temporal behaviour of five picophytoplankton populations in a site of the Tyrrhenian Sea by using a reaction–diffusion–taxis model. The study is performed, taking into account the seasonal changes of environmental variables, obtained starting from experimental findings. The multiplicative noise source, present in all three models, mimics the random fluctuations of temperature and velocity field. The vertical profiles of chlorophyll concentration obtained from the stochastic models show a good agreement with experimental data sampled in the two marine sites considered. The results could be useful to devise a new class of models based on a stochastic approach and able to predict future changes in biomass primary production.     

Evaluation of neuronal coupling dynamics

Temporary correlated activity of neuron assemblies is believed to play a substantial role for the brain's pattern recognition ability. To study the underlying principles of such mechanisms, a method is proposed for the characterization of the interneuronal and stimulus-response coupling changes of two periodically driven and simultaneously recorded units. The coupling measure is derived from the cross correlation function by calculating the actual correlation contributions without performing the subsequent time-average (which would give the cross correlation function). Examples are given for simultaneously recorded spike trains from visual cortical units, but the method can be applied equally well to evoked potentials or intracellular recordings.     

Stochastic nonlinear dynamics pattern formation and growth models

Stochastic evolutionary growth and pattern formation models are treated in a unified way in terms of algorithmic models of nonlinear dynamic systems with feedback built of a standard set of signal processing units. A number of concrete models is described and illustrated by numerous examples of artificially generated patterns that closely imitate wide variety of patterns found in the nature.     

Stochastic models for virus and immune system dynamics

New stochastic models are developed for the dynamics of a viral infection and an immune response during the early stages of infection. The stochastic models are derived based on the dynamics of deterministic models. The simplest deterministic model is a well-known system of ordinary differential equations which consists of three populations: uninfected cells, actively infected cells, and virus particles. This basic model is extended to include some factors of the immune response related to Human Immunodeficiency Virus-1 (HIV-1) infection. For the deterministic models, the basic reproduction number, R₀, is calculated and it is shown that if R₀<1, the disease-free equilibrium is locally asymptotically stable and is globally asymptotically stable in some special cases. The new stochastic models are systems of stochastic differential equations (SDEs) and continuous-time Markov chain (CTMC) models that account for the variability in cellular reproduction and death, the infection process, the immune system activation, and viral reproduction. Two viral release strategies are considered: budding and bursting. The CTMC model is used to estimate the probability of virus extinction during the early stages of infection. Numerical simulations are carried out using parameter values applicable to HIV-1 dynamics. The stochastic models provide new insights, distinct from the basic deterministic models. For the case R₀>1, the deterministic models predict the viral infection persists in the host. But for the stochastic models, there is a positive probability of viral extinction. It is shown that the probability of a successful invasion depends on the initial viral dose, whether the immune system is activated, and whether the release strategy is bursting or budding.     

Stochastic and spatial dynamics of nematode parasites in farmed ruminants

Host-parasite systems provide powerful opportunities for the study of spatial and stochastic effects in ecology; this has been particularly so for directly transmitted microparasites. Here, we construct a fully stochastic model of the population dynamics of a macroparasite system: trichostrongylid gastrointestinal nematode parasites of farmed ruminants. The model subsumes two implicit spatial effects: the host population size (the spatial extent of the interaction between hosts) and spatial heterogeneity ('clumping') in the infection process. This enables us to investigate the roles of several different processes in generating aggregated parasite distributions. The necessity for female worms to find a mate in order to reproduce leads to an Allee effect, which interacts nonlinearly with the stochastic population dynamics and leads to the counter-intuitive result that, when rare, epidemics can be more likely and more severe in small host populations. Clumping in the infection process reduces the strength of this Allee effect, but can hamper the spread of an epidemic by making infection events too rare. Heterogeneity in the hosts' response to infection has to be included in the model to generate aggregation at the level observed empirically.