Experimental Rugged Fitness Landscape in Protein Sequence Space

The fitness landscape in sequence space determines the process of biomolecular evolution. To plot the fitness landscape of protein function, we carried out in vitro molecular evolution beginning with a defective fd phage carrying a random polypeptide of 139 amino acids in place of the g3p minor coat protein D2 domain, which is essential for phage infection. After 20 cycles of random substitution at sites 12–130 of the initial random polypeptide and selection for infectivity, the selected phage showed a 1.7×10(4)-fold increase in infectivity, defined as the number of infected cells per ml of phage suspension. Fitness was defined as the logarithm of infectivity, and we analyzed (1) the dependence of stationary fitness on library size, which increased gradually, and (2) the time course of changes in fitness in transitional phases, based on an original theory regarding the evolutionary dynamics in Kauffman's n-k fitness landscape model. In the landscape model, single mutations at single sites among n sites affect the contribution of k other sites to fitness. Based on the results of these analyses, k was estimated to be 18–24. According to the estimated parameters, the landscape was plotted as a smooth surface up to a relative fitness of 0.4 of the global peak, whereas the landscape had a highly rugged surface with many local peaks above this relative fitness value. Based on the landscapes of these two different surfaces, it appears possible for adaptive walks with only random substitutions to climb with relative ease up to the middle region of the fitness landscape from any primordial or random sequence, whereas an enormous range of sequence diversity is required to climb further up the rugged surface above the middle region

Predictability of evolutionary trajectories in fitness landscapes

Experimental studies on enzyme evolution show that only a small fraction of all possible mutation trajectories are accessible to evolution. However, these experiments deal with individual enzymes and explore a tiny part of the fitness landscape. We report an exhaustive analysis of fitness landscapes constructed with an off-lattice model of protein folding where fitness is equated with robustness to misfolding. This model mimics the essential features of the interactions between amino acids, is consistent with the key paradigms of protein folding and reproduces the universal distribution of evolutionary rates among orthologous proteins. We introduce mean path divergence as a quantitative measure of the degree to which the starting and ending points determine the path of evolution in fitness landscapes. Global measures of landscape roughness are good predictors of path divergence in all studied landscapes: the mean path divergence is greater in smooth landscapes than in rough ones. The model-derived and experimental landscapes are significantly smoother than random landscapes and resemble additive landscapes perturbed with moderate amounts of noise; thus, these landscapes are substantially robust to mutation. The model landscapes show a deficit of suboptimal peaks even compared with noisy additive landscapes with similar overall roughness. We suggest that smoothness and the substantial deficit of peaks in the fitness landscapes of protein evolution are fundamental consequences of the physics of protein folding.Comment: 14 pages, 7 figure

Emergence of terpene cyclization in Artemisia annua

The emergence of terpene cyclization was critical to the evolutionary expansion of chemical diversity yet remains unexplored. Here we report the first discovery of an epistatic network of residues that controls the onset of terpene cyclization in Artemisia annua. We begin with amorpha-4,11-diene synthase (ADS) and (E)-b-farnesene synthase (BFS), a pair of terpene synthases that produce cyclic or linear terpenes, respectively. A library of B27,000 enzymes is generated by breeding combinations of natural amino-acid substitutions from the cyclic into the linear producer. We discover one dominant mutation is sufficient to activate cyclization, and together with two additional residues comprise a network of strongly epistatic interactions that activate, suppress or reactivate cyclization. Remarkably, this epistatic network of equivalent residues also controls cyclization in a BFS homologue from Citrus junos. Fitness landscape analysis of mutational trajectories provides quantitative insights into a major epoch in specialized metabolism

Calculation of the relative metastabilities of proteins in subcellular compartments of Saccharomyces cerevisiae

[abridged] Background: The distribution of chemical species in an open system at metastable equilibrium can be expressed as a function of environmental variables which can include temperature, oxidation-reduction potential and others. Calculations of metastable equilibrium for various model systems were used to characterize chemical transformations among proteins and groups of proteins found in different compartments of yeast cells. Results: With increasing oxygen fugacity, the relative metastability fields of model proteins for major subcellular compartments go as mitochondrion, endoplasmic reticulum, cytoplasm, nucleus. In a metastable equilibrium setting at relatively high oxygen fugacity, proteins making up actin are predominant, but those constituting the microtubule occur with a low chemical activity. A reaction sequence involving the microtubule and spindle pole proteins was predicted by combining the known intercompartmental interactions with a hypothetical program of oxygen fugacity changes in the local environment. In further calculations, the most-abundant proteins within compartments generally occur in relative abundances that only weakly correspond to a metastable equilibrium distribution. However, physiological populations of proteins that form complexes often show an overall positive or negative correlation with the relative abundances of proteins in metastable assemblages. Conclusions: This study explored the outlines of a thermodynamic description of chemical transformations among interacting proteins in yeast cells. The results suggest that these methods can be used to measure the degree of departure of a natural biochemical process or population from a local minimum in Gibbs energy.Comment: 32 pages, 7 figures; supporting information is available at http://www.chnosz.net/yeas

Seasonal variations in the nitrogen isotopic composition of settling particles at station K2 in the western subarctic North Pacific

Intensive observations using hydrographical cruises and moored sediment trap deployments during 2010 and 2012 at station K2 in the North Pacific western subarctic gyre (WSG) revealed seasonal changes in δ15N of both suspended and settling particles. Suspended particles (SUS) were collected from depths between the surface and 200 m; settling particles by drifting traps (DST; 100-200 m) and moored traps (MST; 200 and 500 m). All particles showed higher δ15N values in winter and lower in summer, contrary to the expected by isotopic fractionation during phytoplankton nitrate consumption. We suggest that these observed isotopic patterns are due to ammonium consumption via light-controlled nitrification, which could induce variations in δ15N(SUS) of 0.4-3.1 ‰ in the euphotic zone (EZ). The δ15N(SUS) signature was reflected by δ15 N(DST) despite modifications during biogenic transformation from suspended particles in the EZ. δ15 N enrichment (average: 3.6 ‰) and the increase in C:N ratio (by 1.6) in settling particles suggests year-round contributions of metabolites from herbivorous zooplankton as well as TEPs produced by diatoms. Accordingly, seasonal δ15 N(DST) variations of 2.4-7.0 ‰ showed a significant correlation with primary productivity (PP) at K2. By applying the observed δ15 N(DST) vs. PP regression to δ15 N(MST) of 1.9-8.0 ‰, we constructed the first annual time-series of PP changes in the WSG. Moreover, the monthly export ratio at 500 m was calculated using both estimated PP and measured organic carbon fluxes. Results suggest a 1.6 to 1.8 times more efficient transport of photosynthetically-fixed carbon to the intermediate layers occurs in summer/autumn rather than winter/spring