Evolutionary dynamics of predator-prey systems: an ecological perspective

 Evolution takes place in an ecological setting that typically involves interactions with other organisms. To describe such evolution, a structure is needed which incorporates the simultaneous evolution of interacting species. Here a formal framework for this purpose is suggested, extending from the microscopic interactions between individuals – the immediate cause of natural selection, through the mesoscopic population dynamics responsible for driving the replacement of one mutant phenotype by another, to the macroscopic process of phenotypic evolution arising from many such substitutions. The process of coevolution that results from this is illustrated in the context of predator–prey systems. With no more than qualitative information about the evolutionary dynamics, some basic properties of predator–prey coevolution become evident. More detailed understanding requires specification of an evolutionary dynamic; two models for this purpose are outlined, one from our own research on a stochastic process of mutation and selection and the other from quantitative genetics. Much of the interest in coevolution has been to characterize the properties of fixed points at which there is no further phenotypic evolution. Stability analysis of the fixed points of evolutionary dynamical systems is reviewed and leads to conclusions about the asymptotic states of evolution rather different from those of game-theoretic methods. These differences become especially important when evolution involves more than one species. Received 10 November 1993; received in revised form 25 July 1994     

Scaffold mediated regulation of MAPK signaling and cytoskeletal dynamics: a perspective

Cell migration is critical for many physiological processes and is often misregulated in developmental disorders and pathological conditions including cancer and neurodegeneration. MAPK signaling and the Rho family of proteins are known regulators of cell migration that exert their influence on cellular cytoskeleton during cell adhesion and migration. Here we review data supporting the view that localized ERK signaling mediated through recently identified scaffold proteins may regulate cell migration.     

A systems biology perspective on protein structural dynamics and signal transduction

The functional dynamics of signal transduction through protein interaction networks are determined both by network topology and by the signal processing properties of component proteins. In order to understand the emergent properties of signal transduction networks in terms of information processing, storage and decision making, we not only need to map the so-called 'interactome' but, perhaps more importantly, we also have to understand how the structural dynamics of constituent proteins shape non-linear responses through cooperativity and allostery. Several in silico methods have been developed to identify networks of cooperative residues in proteins and help infer their mode of action. Applying this type of analysis to important classes of modular signal transduction domains should, in principle, allow the function of these proteins to be abstracted in terms of their information processing characteristics, permitting better comprehension of the systemic properties of biological networks.     

Beyond structure: emerging approaches to study GPCR dynamics

1. Download : Download high-res image (360KB)
2. Download : Download full-size image

     

Antagonist binding and induced conformational dynamics of GPCR A2A adenosine receptor

The A_2A adenosine receptor (A_2AAR) is a unique G‐protein coupled receptor (GPCR), because besides agonist, its antagonist could also lead to therapeutic relevance. Based on A_2AAR‐antagonist crystal structure, we have studied the binding mechanism of two distinct antagonists, ZM241385 and KW6002, and dynamic behaviors of A_2AAR induced by antagonist binding. Key residues interacting with both antagonists and residues specifically binding to one of them are identified. ZM241385 specifically bound to S67^2.65, M177^5.38, and N253^6.55, while KW6002 binds to F62^2.60, A81^3.29, and H264^7.29. Moreover, interactions with L167^5.28 are found for both antagonists, which were not reported in agonist binding. The dynamic behaviors of antagonist bound holo‐A_2AARs were found to be different from the apo‐A_2AAR in three typical functional switches, (i) the “ionic lock” was in equilibrium between formation and breakage in apo‐A_2AAR, but stayed broken in holo‐A_2AARs; (ii) the “rotamer toggle switch,” T88^3.36/F242^6.44/W246^6.48, adopted different rotameric conformations in apo‐A_2AAR and holo‐A_2AARs; (iii) apo‐A_2AAR preferred α‐helical intracellular loop (IC)2 and flexible IC3, while holo‐A_2AARs had a flexible IC2 and α‐helical IC3. Our results indicated that antagonist binding induced different conformational rearrangements of these characteristic functional switches in apo‐A_2AAR and holo‐A_2AARs. Proteins 2013; 81:1399–1410. © 2013 Wiley Periodicals, Inc.     

Cancer attractors: A systems view of tumors from a gene network dynamics and developmental perspective

Cell lineage commitment and differentiation are governed by a complex gene regulatory network. Disruption of these processes by inappropriate regulatory signals and by mutational rewiring of the network can lead to tumorigenesis. Cancer cells often exhibit immature or embryonic traits and dysregulated developmental genes can act as oncogenes. However, the prevailing paradigm of somatic evolution and multi-step tumorigenesis, while useful in many instances, offers no logically coherent reason for why oncogenesis recapitulates ontogenesis. The formal concept of “cancer attractors”, derived from an integrative, complex systems approach to gene regulatory network may provide a natural explanation. Here we present the theory of attractors in gene network dynamics and review the concept of cell types as attractors. We argue that cancer cells are trapped in abnormal attractors and discuss this concept in the light of recent ideas in cancer biology, including cancer genomics and cancer stem cells, as well as the implications for differentiation therapy.     

Structural insights into human GPCR protein OA1: a computational perspective

Human ocular albinism type 1 protein (OA1)—a member of the G-protein coupled receptor (GPCR) superfamily—is an integral membrane glycoprotein expressed exclusively by intracellular organelles known as melanocytes, and is responsible for the proper biogenesis of melanosomes. Mutations in the Oa1 gene are responsible for the disease ocular albinism. Despite its clinical importance, there is a lack of in-depth understanding of its structure and mechanism of activation due to the absence of a crystal structure. In the present study, homology modeling was applied to predicting OA1 structure following thorough sequence analysis and secondary structure predictions. The predicted model had the signature residues and motifs expected of GPCRs, and was used for carrying out molecular docking studies with an endogenous ligand, l-DOPA and an antagonist, dopamine; the results agreed quite well with the available experimental data. Finally, three sets of explicit molecular dynamics simulations were carried out in lipid bilayer, the results of which not only confirmed the stability of the predicted model, but also helped witness some differences in structural features such as rotamer toggle switch, helical tilts and hydrogen bonding pattern that helped distinguish between the agonist- and antagonist-bound receptor forms. In place of the typical “D/ERY”-motif-mediated “ionic lock”, a hydrogen bond mediated by the “DAY” motif was observed that could be used to distinguish the agonist and antagonist bound forms of OA1. In the absence of a crystal structure, this study helped to shed some light on the structural features of OA1, and its behavior in the presence of an agonist and an antagonist, which might be helpful in the future drug discovery process for ocular albinism.     

Sex and sugar in yeast: two distinct GPCR systems

Although eukaryotic G-protein coupled receptor (GPCR) systems are well known for their ability to detect and mediate rapid responses to extracellular signals, the full range of stimuli to which they respond may not yet have been identified. Activation of GPCRs by hormones, pheromones, odorants, neurotransmitters, light and different taste compounds is well established. However, the recent discovery of a glucose-sensing GPCR system in Saccharomyces cerevisiae has unexpectedly added common nutrients to this list of stimuli. This GPCR system mediates glucose activation of adenylate cyclase during the switch from respirative/gluconeogenic metabolism to fermentation. The GPCR system involved in pheromone signalling in S. cerevisiae has already served as an important model and tool for the study of GPCR systems in higher eukaryotic cell types. Here, we highlight the similarities and differences between these two signalling systems. We also indicate how the new glucose-sensing system can serve as a model for GPCR function and as a tool with which to screen for heterologous components of signalling pathways as well as for novel ligands in high-throughput assays.     

A model of propagating calcium-induced calcium release mediated by calcium diffusion

The effect of sudden local fluctuations of the free sarcoplasmic [Ca++]i in cardiac cells on calcium release and calcium uptake by the sarcoplasmic reticulum (SR) was calculated with the aid of a simplified model of SR calcium handling. The model was used to evaluate whether propagation of calcium transients and the range of propagation velocities observed experimentally (0.05-15 mm s(-1)) could be predicted. Calcium fluctuations propagate by virtue of focal calcium release from the SR, diffusion through the cytosol (which is modulated by binding to troponin and calmodulin and sequestration by the SR), and subsequently induce calcium release from adjacent release sites of the SR. The minimal and maximal velocities derived from the simulation were 0.09 and 15 mm s(-1) respectively. The method of solution involved writing the diffusion equation as a difference equation in the spatial coordinates. Thus, coupled ordinary differential equations in time with banded coefficients were generated. The coupled equations were solved using Gear's sixth order predictor-corrector algorithm for stiff equations with reflective boundaries. The most important determinants of the velocity of propagation of the calcium waves were the diastolic [Ca++]i, the rate of rise of the release, and the amount of calcium released from the SR. The results are consistent with the assumptions that calcium loading causes an increase in intracellular calcium and calcium in the SR, and an increase in the amount and rate of calcium released. These two effects combine to increase the propagation velocity at higher levels of calcium loading.     

RNA structure and dynamics: A base pairing perspective

RNA is now known to possess various structural, regulatory and enzymatic functions for survival of cellular organisms. Functional RNA structures are generally created by three-dimensional organization of small structural motifs, formed by base pairing between self-complementary sequences from different parts of the RNA chain. In addition to the canonical Watson–Crick or wobble base pairs, several non-canonical base pairs are found to be crucial to the structural organization of RNA molecules. They appear within different structural motifs and are found to stabilize the molecule through long-range intra-molecular interactions between basic structural motifs like double helices and loops. These base pairs also impart functional variation to the minor groove of A-form RNA helices, thus forming anchoring site for metabolites and ligands. Non-canonical base pairs are formed by edge-to-edge hydrogen bonding interactions between the bases. A large number of theoretical studies have been done to detect and analyze these non-canonical base pairs within crystal or NMR derived structures of different functional RNA. Theoretical studies of these isolated base pairs using ab initio quantum chemical methods as well as molecular dynamics simulations of larger fragments have also established that many of these non-canonical base pairs are as stable as the canonical Watson–Crick base pairs. This review focuses on the various structural aspects of non-canonical base pairs in the organization of RNA molecules and the possible applications of these base pairs in predicting RNA structures with more accuracy.     

Metabolic regulation: A control analytic perspective

A possible basis for a quantitative theory of metabolic regulation is outlined. Regulation is defined here as the alteration of reaction properties to augment or counteract the mass-action trend in a network reactions. In living systems the enzymes that catalyze these reactions are the handles through which such alteration is effected. It is shown how the elasticity coefficients of an enzyme-catalyzed reaction with respect to substrates and products are the sum of a massaction term and a regulatory kinetic term; these coefficients therefore distinguish between massaction effects and regulatory effects and are recognized as the key to quantifying regulation. As elasticity coefficients are also basic ingredients of metabolic control analysis, it is possible to relate regulation to such concepts as control, signalling, stability, and homeostasis. The need for care in the choice of relative or absolute changes when considering questions of metabolic regulation is stressed. Although the concepts are illustrated in terms of a simple coupled reaction system, they apply equally to more complex systems. When such systems are divided into reaction blocks, co-response coefficients can be used to measure the elasticities of these blocks.I dedicate this paper to Henrik Kacser, co-founder of and guiding light in the field of metabolic control analysis. His recent death leaves us bereft of a fount of wisdom and kindness, but his work remains as a monument along the path of our search for an understanding of metabolic behavior.     

Expanding GPCR homology model binding sites via a balloon potential: A molecular dynamics refinement approach

Homology modeling of G protein-coupled receptors is becoming a widely used tool in drug discovery. However, unrefined models built using the bovine rhodopsin crystal structure as the template, often have binding sites that are too small to accommodate known ligands. Here, we present a novel systematic method to refine model active sites based on a pressure-guided molecular dynamics simulation. A distinct advantage of this approach is the ability to introduce systematic perturbations in model backbone atoms in addition to side chain adjustments. The method is validated on two test cases: (1) docking of retinal into an MD-relaxed structure of opsin and (2) docking of known ligands into a homology model of the CCR2 receptor. In both cases, we show that the MD expansion algorithm makes it possible to dock the ligands in poses that agree with the crystal structure or mutagenesis data.     

Regime shifts in marine communities: a complex systems perspective on food web dynamics

Species composition and habitats are changing at unprecedented rates in the world''s oceans, potentially causing entire food webs to shift to structurally and functionally different regimes. Despite the severity of these regime shifts, elucidating the precise nature of their underlying processes has remained difficult. We address this challenge with a new analytic approach to detect and assess the relative strength of different driving processes in food webs. Our study draws on complexity theory, and integrates the network-centric exponential random graph modelling (ERGM) framework developed within the social sciences with community ecology. In contrast to previous research, this approach makes clear assumptions of direction of causality and accommodates a dynamic perspective on the emergence of food webs. We apply our approach to analysing food webs of the Baltic Sea before and after a previously reported regime shift. Our results show that the dominant food web processes have remained largely the same, although we detect changes in their magnitudes. The results indicate that the reported regime shift may not be a system-wide shift, but instead involve a limited number of species. Our study emphasizes the importance of community-wide analysis on marine regime shifts and introduces a novel approach to examine food webs.     

A calcium-induced calcium release mechanism mediated by calsequestrin

Calcium (Ca²⁺)-induced Ca²⁺ release (CICR) is widely accepted as the principal mechanism linking electrical excitation and mechanical contraction in cardiac cells. The CICR mechanism has been understood mainly based on binding of cytosolic Ca²⁺ with ryanodine receptors (RyRs) and inducing Ca²⁺ release from the sarcoplasmic reticulum (SR). However, recent experiments suggest that SR lumenal Ca²⁺ may also participate in regulating RyR gating through calsequestrin (CSQ), the SR lumenal Ca²⁺ buffer. We investigate how SR Ca²⁺ release via RyR is regulated by Ca²⁺ and calsequestrin (CSQ). First, a mathematical model of RyR kinetics is derived based on experimental evidence. We assume that the RyR has three binding sites, two cytosolic sites for Ca²⁺ activation and inactivation, and one SR lumenal site for CSQ binding. The open probability (P_o) of the RyR is found by simulation under controlled cytosolic and SR lumenal Ca²⁺. Both peak and steady-state P_o effectively increase as SR lumenal Ca²⁺ increases. Second, we incorporate the RyR model into a CICR model that has both a diadic space and the junctional SR (jSR). At low jSR Ca²⁺ loads, CSQs are more likely to bind with the RyR and act to inhibit jSR Ca²⁺ release, while at high SR loads CSQs are more likely to detach from the RyR, thereby increasing jSR Ca²⁺ release. Furthermore, this CICR model produces a nonlinear relationship between fractional jSR Ca²⁺ release and jSR load. These findings agree with experimental observations in lipid bilayers and cardiac myocytes.