Computational mechanics analysis of the hidden conformational dynamics within a molecular trajectory

DUE TO COPYRIGHT RESTRICTIONS ONLY AVAILABLE FOR CONSULTATION AT ASTON UNIVERSITY LIBRARY AND INFORMATION SERVICES WITH PRIOR ARRANGEMENT

Description and training of neural network dynamics

Attractor properties of a popular discrete-time neural network model are illustrated through numerical simulations. The most complex dynamics is found to occur within particular ranges of parameters controlling the symmetry and magnitude of the weight matrix. A small network model is observed to produce fixed points, limit cycles, mode-locking, the Ruelle-Takens route to chaos, and the period-doubling route to chaos. Training algorithms for tuning this dynamical behaviour are discussed. Training can be an easy or difficult task, depending whether the problem requires the use of temporal information distributed over long time intervals. Such problems require training algorithms which can handle hidden nodes. The most prominent of these algorithms, back propagation through time, solves the temporal credit assignment problem in a way which can work only if the relevant information is distributed locally in time. The Moving Targets algorithm works for the more general case, but is computationally intensive, and prone to local minima.

The dynamics of a genetic algorithm under stabilizing selection

A formalism recently introduced by Prugel-Bennett and Shapiro uses the methods of statistical mechanics to model the dynamics of genetic algorithms. To be of more general interest than the test cases they consider. In this paper, the technique is applied to the subset sum problem, which is a combinatorial optimization problem with a strongly non-linear energy (fitness) function and many local minima under single spin flip dynamics. It is a problem which exhibits an interesting dynamics, reminiscent of stabilizing selection in population biology. The dynamics are solved under certain simplifying assumptions and are reduced to a set of difference equations for a small number of relevant quantities. The quantities used are the population's cumulants, which describe its shape, and the mean correlation within the population, which measures the microscopic similarity of population members. Including the mean correlation allows a better description of the population than the cumulants alone would provide and represents a new and important extension of the technique. The formalism includes finite population effects and describes problems of realistic size. The theory is shown to agree closely to simulations of a real genetic algorithm and the mean best energy is accurately predicted.

The learning dynamics of a universal approximator

The learning properties of a universal approximator, a normalized committee machine with adjustable biases, are studied for on-line back-propagation learning. Within a statistical mechanics framework, numerical studies show that this model has features which do not exist in previously studied two-layer network models without adjustable biases, e.g., attractive suboptimal symmetric phases even for realizable cases and noiseless data.

Modelling the dynamics of genetic algorithms using statistical mechanics

A formalism for modelling the dynamics of Genetic Algorithms (GAs) using methods from statistical mechanics, originally due to Prugel-Bennett and Shapiro, is reviewed, generalized and improved upon. This formalism can be used to predict the averaged trajectory of macroscopic statistics describing the GA's population. These macroscopics are chosen to average well between runs, so that fluctuations from mean behaviour can often be neglected. Where necessary, non-trivial terms are determined by assuming maximum entropy with constraints on known macroscopics. Problems of realistic size are described in compact form and finite population effects are included, often proving to be of fundamental importance. The macroscopics used here are cumulants of an appropriate quantity within the population and the mean correlation (Hamming distance) within the population. Including the correlation as an explicit macroscopic provides a significant improvement over the original formulation. The formalism is applied to a number of simple optimization problems in order to determine its predictive power and to gain insight into GA dynamics. Problems which are most amenable to analysis come from the class where alleles within the genotype contribute additively to the phenotype. This class can be treated with some generality, including problems with inhomogeneous contributions from each site, non-linear or noisy fitness measures, simple diploid representations and temporally varying fitness. The results can also be applied to a simple learning problem, generalization in a binary perceptron, and a limit is identified for which the optimal training batch size can be determined for this problem. The theory is compared to averaged results from a real GA in each case, showing excellent agreement if the maximum entropy principle holds. Some situations where this approximation brakes down are identified. In order to fully test the formalism, an attempt is made on the strong sc np-hard problem of storing random patterns in a binary perceptron. Here, the relationship between the genotype and phenotype (training error) is strongly non-linear. Mutation is modelled under the assumption that perceptron configurations are typical of perceptrons with a given training error. Unfortunately, this assumption does not provide a good approximation in general. It is conjectured that perceptron configurations would have to be constrained by other statistics in order to accurately model mutation for this problem. Issues arising from this study are discussed in conclusion and some possible areas of further research are outlined.

Modelling Apoptotic cell clearance within the atherosclerotic plaque

The aged population have an increased susceptibility to infection, therefore function of the innate immune system may be impaired as we age. Macrophages, and their precursors monocytes, play an important role in host defence in the form of phagocytosis, and also link the innate and adaptive immune system via antigen presentation. Classically-activated 'M1' macrophages are pro-inflammatory, which can be induced by encountering pathogenic material or pro-inflammatory mediators. Alternatively activated 'M2' macrophages have a largely reparative role, including clearance of apoptotic bodies and debris from tissues. Despite some innate immune receptors being implicated in the clearance of apoptotic cells, the process has been observed to have a dominant anti-inflammatory phenotype with cytokines such as IL-10 and TGF-ß being implicated. The atherosclerotic plaque contains recruited monocytes and macrophages, and is a highly inflammatory environment despite high levels of apoptosis. At these sites, monocytes differentiate into macrophages and gorge on lipoproteins, resulting in formation of 'foam cells' which then undergo apoptosis, recruiting further monocytes. This project seeks to understand why, given high levels of apoptosis, the plaque is a pro-inflammatory environment. This phenomenon may be the result of the aged environment or an inability of foam cells to elicit an anti-inflammatory effect in response to dying cells. Here we demonstrate that lipoprotein treatment of macrophages in culture results in reduced capacity to clear apoptotic cells. The effect of lipoprotein treatment on apoptotic cell-mediated immune modulation of macrophage function is currently under study.

Dynamics of on-line learning in radial basis function networks

On-line learning is examined for the radial basis function network, an important and practical type of neural network. The evolution of generalization error is calculated within a framework which allows the phenomena of the learning process, such as the specialization of the hidden units, to be analyzed. The distinct stages of training are elucidated, and the role of the learning rate described. The three most important stages of training, the symmetric phase, the symmetry-breaking phase, and the convergence phase, are analyzed in detail; the convergence phase analysis allows derivation of maximal and optimal learning rates. As well as finding the evolution of the mean system parameters, the variances of these parameters are derived and shown to be typically small. Finally, the analytic results are strongly confirmed by simulations.

Structure, dynamics, and energetics of siRNA-cationic vector complexation:a molecular dynamics study

The design and synthesis of safe and efficient nonviral vectors for gene delivery has attracted significant attention in recent years. Previous experiments have revealed that the charge density of a polycation (the carrier) plays a crucial role in complexation and the release of the gene from the complex in the cytosol. In this work, we adopt an atomistic molecular dynamics simulation approach to study the complexation of short strand duplex RNA with six cationic carrier systems of varying charge and surface topology. The simulations reveal detailed molecular-level pictures of the structures and dynamics of the RNA-polycation complexes. Estimates for the binding free energy indicate that electrostatic contributions are dominant followed by van der Waals interactions. The binding free energy between the 8(+)polymers and the RNA is found to be larger than that of the 4(+)polymers, in general agreement with previously published data. Because reliable binding free energies provide an effective index of the ability of the polycationic carrier to bind the nucleic acid and also carry implications for the process of gene release within the cytosol, these novel simulations have the potential to provide us with a much better understanding of key mechanistic aspects of gene-polycation complexation and thereby advance the rational design of nonviral gene delivery systems.

The dynamics of HRM systems in Indian BPO firms

This paper highlights the context within which business process outsourcing (BPO) has rapidly grown in India and the critical need to investigate the dynamics of human resource management (HRM) practices and systems in this sector. Using a mixed-method approach involving both in-depth interviews and self-completing questionnaires, we analyze the nature of HRM systems in BPO organizations operating in India. The analysis is based on a sample of 51 BPO companies, a majority of which are located near the capital of New Delhi. The results focus on the nature and structure of work and organization of Indian BPOs, as well as the strategic role played by HRM in such organizations. Furthermore, the findings highlight the way specific HRM practices such as recruitment, performance appraisal, training and development, and compensations are implemented. Our study suggests the existence of formal, structured, and rationalized HRM systems in Indian BPOs. A number of insights related to HRM policies and practices are shared by the HR managers interviewed shedding more light on the inner workings of the Indian BPO companies and their challenges. The analysis provides original and useful information to both academics and practitioners and opens avenues for future research on the nature of HRM systems and practices in the Indian BPO industry.

Using system dynamics in a discrete-event simulation study of a manufacturing plant

Purpose - To provide an example of the use of system dynamics within the context of a discrete-event simulation study. Design/methodology/approach - A discrete-event simulation study of a production-planning facility in a gas cylinder-manufacturing plant is presented. The case study evidence incorporates questionnaire responses from sales managers involved in the order-scheduling process. Findings - As the project progressed it became clear that, although the discrete-event simulation would meet the objectives of the study in a technical sense, the organizational problem of "delivery performance" would not be solved by the discrete-event simulation study alone. The case shows how the qualitative outcomes of the discrete-event simulation study led to an analysis using the system dynamics technique. The system dynamics technique was able to model the decision-makers in the sales and production process and provide a deeper understanding of the performance of the system. Research limitations/implications - The case study describes a traditional discrete-event simulation study which incorporated an unplanned investigation using system dynamics. Further, case studies using a planned approach to showing consideration of organizational issues in discrete-event simulation studies are required. Then the role of both qualitative data in a discrete-event simulation study and the use of supplementary tools which incorporate organizational aspects may help generate a methodology for discrete-event simulation that incorporates human aspects and so improve its relevance for decision making. Practical implications - It is argued that system dynamics can provide a useful addition to the toolkit of the discrete-event simulation practitioner in helping them incorporate a human aspect in their analysis. Originality/value - Helps decision makers gain a broader perspective on the tools available to them by showing the use of system dynamics to supplement the use of discrete-event simulation. © Emerald Group Publishing Limited.

Marketing strategy : making sense of industry dynamics:is there still value in strategic group research?

This paper looks selectively at strategic group theory and seeks to explore the benefits and limitations of modern strategic group analysis within the context of practical strategy making in the pharmaceutical industry. The rise and fall of strategic group research is reviewed and suggestions advanced as to the reasons why strategic group research suffered criticism and frequently produced conflicting results. The paper concludes that strategic group research offers a valuable way to classify firms by their strategy and provides some suggestions as to how industry strategists may benefit from strategic group analysis and avoid the pitfalls exposed by previous research.

The effect of pH on PAMAM dendrimer-siRNA complexation:endosomal considerations as determined by molecular dynamics simulation

Intracellular degradation of genes, most notably within the endo-lysosomal compartment is considered a significant barrier to (non-viral) gene delivery in vivo. Previous reports based on in vitro studies claim that carriers possessing a mixture of primary, secondary and tertiary amines are able to buffer the acidic environment within the endosome, allowing for timely release of their contents, leading to higher transfection rates. In this report, we adopt an atomistic molecular dynamics (MD) simulation approach, comparing the complexation of 21-bp siRNA with low-generation polyamidoamine (PAMAM) dendrimers (G0 and G1) at both neutral and acidic pHs, the latter of which mimics the degradative environment within maturing 'late-endosomes'. Our simulations reveal that the time taken for the dendrimer-gene complex (dendriplex) to reach equilibrium is appreciably longer at low pH and this is accompanied by more compact packaging of the dendriplex, as compared to simulations performed at neutral pH. We also note larger absolute values of calculated binding free energies of the dendriplex at low pH, indicating a higher dendrimer-nucleic acid affinity in comparison with neutral pH. These novel simulations provide a more detailed understanding of low molecular-weight polymer-siRNA behavior, mimicking the endosomal environment and provide input of direct relevance to the "proton sponge theory", thereby advancing the rational design of non-viral gene delivery systems.

The spatial distribution and temporal dynamics of brain regions activated during the perception of object and non-object patterns

Both animal and human studies suggest that the efficiency with which we are able to grasp objects is attributable to a repertoire of motor signals derived directly from vision. This is in general agreement with the long-held belief that the automatic generation of motor signals by the perception of objects is based on the actions they afford. In this study, we used magnetoencephalography (MEG) to determine the spatial distribution and temporal dynamics of brain regions activated during passive viewing of object and non-object targets that varied in the extent to which they afforded a grasping action. Synthetic Aperture Magnetometry (SAM) was used to localize task-related oscillatory power changes within specific frequency bands, and the time course of activity within given regions-of-interest was determined by calculating time-frequency plots using a Morlet wavelet transform. Both single subject and group-averaged data on the spatial distribution of brain activity are presented. We show that: (i) significant reductions in 10-25 Hz activity within extrastriate cortex, occipito-temporal cortex, sensori-motor cortex and cerebellum were evident with passive viewing of both objects and non-objects; and (ii) reductions in oscillatory activity within the posterior part of the superior parietal cortex (area Ba7) were only evident with the perception of objects. Assuming that focal reductions in low-frequency oscillations (< 30 Hz) reflect areas of heightened neural activity, we conclude that: (i) activity within a network of brain areas, including the sensori-motor cortex, is not critically dependent on stimulus type and may reflect general changes in visual attention; and (ii) the posterior part of the superior parietal cortex, area Ba7, is activated preferentially by objects and may play a role in computations related to grasping. © 2006 Elsevier Inc. All rights reserved.

Flow field analysis of some mixing and conveying screw element regions, within a closely intermeshing, co-rotating twin-screw extruder.

Grafting of antioxidants and other modifiers onto polymers by reactive extrusion, has been performed successfully by the Polymer Processing and Performance Group at Aston University. Traditionally the optimum conditions for the grafting process have been established within a Brabender internal mixer. Transfer of this batch process to a continuous processor, such as an extruder, has, typically, been empirical. To have more confidence in the success of direct transfer of the process requires knowledge of, and comparison between, residence times, mixing intensities, shear rates and flow regimes in the internal mixer and in the continuous processor.The continuous processor chosen for the current work in the closely intermeshing, co-rotating twin-screw extruder (CICo-TSE). CICo-TSEs contain screw elements that convey material with a self-wiping action and are widely used for polymer compounding and blending. Of the different mixing modules contained within the CICo-TSE, the trilobal elements, which impose intensive mixing, and the mixing discs, which impose extensive mixing, are of importance when establishing the intensity of mixing. In this thesis, the flow patterns within the various regions of the single-flighted conveying screw elements and within both the trilobal element and mixing disc zones of a Betol BTS40 CICo-TSE, have been modelled using the computational fluid dynamics package Polyflow. A major obstacle encountered when solving the flow problem within all of these sets of elements, arises from both the complex geometry and the time-dependent flow boundaries as the elements rotate about their fixed axes. Simulation of the time dependent boundaries was overcome by selecting a number of sequential 2D and 3D geometries, used to represent partial mixing cycles. The flow fields were simulated using the ideal rheological properties of polypropylene and characterised in terms of velocity vectors, shear stresses generated and a parameter known as the mixing efficiency. The majority of the large 3D simulations were performed on the Cray J90 supercomputer situated at the Rutherford-Appleton laboratories, with pre- and postprocessing operations achieved via a Silicon Graphics Indy workstation. A mechanical model was constructed consisting of various CICo-TSE elements rotating within a transparent outer barrel. A technique has been developed using coloured viscous clays whereby the flow patterns and mixing characteristics within the CICo-TSE may be visualised. In order to test and verify the simulated predictions, the patterns observed within the mechanical model were compared with the flow patterns predicted by the computational model. The flow patterns within the single-flighted conveying screw elements in particular, showed good agreement between the experimental and simulated results.

Investigation of chaotic mixing in laminar fluid systems by the use of computational fluid dynamics

This work presents significant development into chaotic mixing induced through periodic boundaries and twisting flows. Three-dimensional closed and throughput domains are shown to exhibit chaotic motion under both time periodic and time independent boundary motions, A property is developed originating from a signature of chaos, sensitive dependence to initial conditions, which successfully quantifies the degree of disorder withjn the mixing systems presented and enables comparisons of the disorder throughout ranges of operating parameters, This work omits physical experimental results but presents significant computational investigation into chaotic systems using commercial computational fluid dynamics techniques. Physical experiments with chaotic mixing systems are, by their very nature, difficult to extract information beyond the recognition that disorder does, does not of partially occurs. The initial aim of this work is to observe whether it is possible to accurately simulate previously published physical experimental results through using commercial CFD techniques. This is shown to be possible for simple two-dimensional systems with time periodic wall movements. From this, and subsequent macro and microscopic observations of flow regimes, a simple explanation is developed for how boundary operating parameters affect the system disorder. Consider the classic two-dimensional rectangular cavity with time periodic velocity of the upper and lower walls, causing two opposing streamline motions. The degree of disorder within the system is related to the magnitude of displacement of individual particles within these opposing streamlines. The rationale is then employed in this work to develop and investigate more complex three-dimensional mixing systems that exhibit throughputs and time independence and are therefore more realistic and a significant advance towards designing chaotic mixers for process industries. Domains inducing chaotic motion through twisting flows are also briefly considered. This work concludes by offering possible advancements to the property developed to quantify disorder and suggestions of domains and associated boundary conditions that are expected to produce chaotic mixing.