Microscopic Self-Organization of Enzymic Reactions in Small Volumes

Pattern formation and localized structures in reaction-diffusion systems with non-Fickian transport

We study the robust dynamical behaviors of reaction-diffusion systems where the transport gives rise to non-Fickian diffusion. A prototype model describing the deposition of molecules in a surface is used to show the generic appearance of Turing structures which can coexist with homogeneous states giving rise to localized structures through the pinning mechanism. The characteristic lengths of these structures are in the nanometer region in agreement with recent experimental observations.

Functional conformational motions in the turnover cycle of cholesterol oxidase

Reexamining experimental data of single-molecule fluorescence correlation spectroscopy for cholesterol oxidase, we find that the existing Michaelis-Menten models with dynamical disorder cannot explain strong correlations between subsequent turnover cycles revealed in the diagonal feature in the joint statistical distribution of adjacent "on" times of this enzyme. We suggest that functional conformational motions representing ordered sequences of transitions between a set of conformational substates are involved, along with equilibrium conformational fluctuations in the turnover cycle of cholesterol oxidase. A two-channel model of single-enzyme dynamics, including a slow functional conformational motion in one of the channels, is proposed that allows us to reproduce such strong correlations.

Interplay between chemical reactions and transport in structured spaces

The main motivation behind this study is to understand the interplay between the reactions and transport in a geometries that are not compact. Typical examples of compact geometries are a box or a sphere. A network made of containers C(1) , C(2),..., C(N) and tubes is an example of the space that is structured and noncompact. In containers, particles react with the rate lambda. Tubes connecting containers allow for the exchange of chemicals with the transport rate D. A situation is considered where a number of reactants is small and kinetics is noise dominated. A method is presented that can be used to calculate the average and higher moments of the reaction time. A number of different chemical reactions are studied and their performance compared in various ways. Two schemes are discussed in general, the reaction on a fixed geometry ensemble (ROGE) and the geometry on a fixed reaction ensemble, examples are given in the ROGE case. The most important findings are as follows. (i) There is a large number of reactions that run faster in a networklike geometry. Such reactions contain antagonistic catalytic influences in the intermediate stages of a reaction scheme that are best dealt with in a networklike structure. (ii) Antagonistic catalytic influences are hard to identify since they are strongly connected to the pattern of injected molecules (inject pattern) and depend on the choice of molecules that have to be synthesized at the end (task pattern). (iii) The reaction time depends strongly on the details of the inject and task patterns.

Nonequilibrium phase transition in a self-activated biological network

We present a lattice model for a two-dimensional network of self-activated biological structures with a diffusive activating agent. The model retains basic and simple properties shared by biological systems at various observation scales, so that the structures can consist of individuals, tissues, cells, or enzymes. Upon activation, a structure emits a new mobile activator and remains in a transient refractory state before it can be activated again. Varying the activation probability, the system undergoes a nonequilibrium second-order phase transition from an active state, where activators are present, to an absorbing, activator-free state, where each structure remains in the deactivated state. We study the phase transition using Monte Carlo simulations and evaluate the critical exponents. As they do not seem to correspond to known values, the results suggest the possibility of a separate universality class.

Conformational-relaxation models of single-enzyme kinetics

Fluorescent spectroscopy experiments with single-enzyme molecules yield a large volume of statistical data that can be analyzed and interpreted using stochastic models of enzyme action. Here, we present two models, each based on the mechanism that an enzyme molecule must pass through a sequence of conformational transformations to complete its catalytic turnover cycle. In the simplest model, only one path leading to the release of product is present. In contrast to this, two different catalytic paths are possible in the second considered model. If a cycle is started from an active state, immediately after the previous product release, it follows a different conformational route and is much shorter. Our numerical investigations show that both models generate non-Markovian molecular statistics. However, their memory landscapes and distributions of cycle times are significantly different. The memory landscape of the double-path model bears strong similarity to the recent experimental data for horseradish peroxidase.

Microfluidic device for combinatorial fusion of liposomes and cells

We describe an electrofusion-based technique for combinatorial synthesis of individual liposomes. A prototype device with containers for liposomes of different compositions and a fusion container was constructed. The sample containers had fluid contact with the fusion container through microchannels. Optical trapping was used to transport individual liposomes and cells through the microchannels into the fusion container. In the fusion container, selected pairs of liposomes were fused together using microelectrodes. A large number of combinatorially synthesized liposomes with complex compositions and reaction systems can be obtained from small sets of precursor liposomes. The order of different reaction steps can be specified and defined by the fusion sequence. This device could also facilitate single cell-cell electrofusions (hybridoma production). This is exemplified by fusion of transported red blood cells.

Oscillatory behavior of a simple kinetic model for proteolysis during cell invasion

Extracellular proteolysis during cell invasion is thought to be tightly organized, both temporally and spatially. This work presents a simple kinetic model that describes the interactions between extracellular matrix (ECM) proteins, proteinases, proteolytic fragments, and integrins. Nonmonotonous behavior arises from enzyme de novo synthesis consecutive to integrin binding to fragments or entire proteins. The model has been simulated using realistic values for kinetic constants and protein concentrations, with fibronectin as the ECM protein. The simulations show damped oscillations of integrin-complex concentrations, indicating alternation of maximal adhesion periods with maximal mobility periods. Comparisons with experimental data from the literature confirm the similarity between this system behavior and cell invasion. The influences on the system of cryptic functions of ECM proteins, proteinase inhibitors, and soluble antiadhesive peptides were examined. The first critical parameter for oscillation is the discrepancy between integrin affinity for intact ECM proteins and the respective proteolytic fragments, thus emphasizing the importance of cryptic functions of ECM proteins in cell invasion. Another critical parameter is the ratio between proteinase and the initial ECM protein concentration. These results suggest new insights into the organization of the ECM degradation during cell invasion.

Modelling Escherichia coli. The concept of competitive coherence

Many of the characteristics of life are plainly evident. Others are not. In principle, computer models might be devised to test whether a hypothetical characteristic is likely to be a characteristic of real life. Here, I propose that the bacterium Escherichia coli can be considered as passing through a series of states in which a distinct set of its constituent molecules or 'elements' are active. The activity of these elements is determined by a competition between two processes. One of these processes depends on the previous cell state whilst the other depends on the internal coherence of the developing state. The simultaneous operation of these two processes I term 'competitive coherence'. To clarify and eventually test these related hypotheses, a possible computer model is outlined.

Self-organizing molecular networks

Strong diffusional mixing and short delivery times typical for micrometer and sub-micrometer reaction volumes lead to a special situation where the turnover times of individual enzyme molecules become the largest characteristic time scale of the chemical kinetics. Under these conditions, populations of cross-regulating allosteric enzymes form molecular networks that exhibit various kinds of self-organized coherent collective dynamics.

Enzymatic reactions in small spatial volumes: comment on a model of Hess and Mikhailov

Recently Hess and Mikhailov pointed out that in small subcellular compartments diffusion is so fast that mixing is instantaneous on the time scale of many enzymatic reactions. This opens the possibility for synchronizing individual reaction events. To illustrate this fact they discuss as example an irreversible enzymatic reaction with allosteric product activation. Under appropriate conditions their model shows coherent spiking in the number of product molecules, caused by the strong correlation between reaction events. In this model only substrate binding is an indeterministic process, all other subsequent transitions between different enzyme states being deterministic, contrary to real processes. The purpose of the present paper was to investigate this interesting phenomenon by means of a more realistic modification of the original model, with only probabilistic transitions. In an attempt to obtain spiking, which was not observed under these conditions, the model was extended to make a clear distinction between allosteric high and low affinity substrate binding, in contrast to the original model using a product dependent mean binding probability. However no periodic signal was detectable in the indeterministic version of the Hess Mikhailov model or the extended version, either by means of direct visualization or on autocorrelation or Fourier analysis. Reasons why spiking is not observed in indeterministic enzyme models are discussed.