共查询到20条相似文献,搜索用时 0 毫秒
1.
We examined properties of adaptive walks by the fittest on “rough Mt. Fuji-type” fitness landscapes, which are modeled by
superposing small uncorrelated random component on an additive fitness landscape. A single adaptive walk is carried out by
repetition of the evolution cycle composed of (1) mutagenesis process that produces random d-fold point mutants of population size N and (2) selection process that picks out the fittest mutant among them. To comprehend trajectories of the walkers, the fitness
landscape is mapped into a (x, y, z)-space, where x, y and z represent, respectively, normalized Hamming distance from the peak on the additive fitness landscape, scaled additive fitness
and scaled non-additive fitness. Thus a single adaptive walk is expressed as the dynamics of a particle in this space. We
drew the “hill-climbing” vector field, where each vector represents the most probable step for a walker in a single step.
Almost all of the walkers are expected to move along streams of vectors existing on a particular surface that overlies the
(x, y)-plane, toward the neighborhood of a characteristic point at which a mutation-selection-random drift balance is reached.
We could theoretically predict this reachable point in the case of random sampling search strategy.
Received: 1 March 2000 / Published online: 3 August 2000 相似文献
2.
We have theoretically studied the statistical properties of adaptive walks (or hill-climbing) on a Mt. Fuji-type fitness landscape in the multi-dimensional sequence space through mathematical analysis and computer simulation. The adaptive walk is characterized by the "mutation distance" d as the step-width of the walker and the "population size" N as the number of randomly generated d-fold point mutants to be screened. In addition to the fitness W, we introduced the following quantities analogous to thermodynamical concepts: "free fitness" G(W) is identical with W+T x S(W), where T is the "evolutionary temperature" T infinity square root of d/lnN and S(W) is the entropy as a function of W, and the "evolutionary force" X is identical with d(G(W)/T)/dW, that is caused by the mutation and selection pressure. It is known that a single adaptive walker rapidly climbs on the fitness landscape up to the stationary state where a "mutation-selection-random drift balance" is kept. In our interpretation, the walker tends to the maximal free fitness state, driven by the evolutionary force X. Our major findings are as follows: First, near the stationary point W*, the "climbing rate" J as the expected fitness change per generation is described by J approximately L x X with L approximately V/2, where V is the variance of fitness distribution on a local landscape. This simple relationship is analogous to the well-known Einstein relation in Brownian motion. Second, the "biological information gain" (DeltaG/T) through adaptive walk can be described by combining the Shannon's information gain (DeltaS) and the "fitness information gain" (DeltaW/T). 相似文献
3.
Uchiyama H Inaoka T Ohkuma-Soyejima T Togame H Shibanaka Y Yoshimoto T Kokubo T 《Journal of biochemistry》2000,128(3):441-447
Prolyl endopeptidase is the only endopeptidase that specifically cleaves peptides at proline residues. Although this unique specificity is advantageous for application in protein chemistry, the stability of the enzyme is lower than those of commonly used peptidases such as subtilisin and trypsin. Therefore, we attempted to apply a directed evolution system to improve the thermostability of the enzyme. First, an efficient expression system for the enzyme in Escherichia coli was established using the prolyl endopeptidase gene from Flavobacterium meningosepticum. Then, a method for screening thermostable variants was developed by combining heat treatment with active staining on membrane filters. Random mutagenesis by error-prone PCR and screening was repeated three times, and as a result the thermostability of the enzyme was increased step by step as the amino acid substitutions accumulated. The most thermostable mutant obtained after the third cycle, PEP-407, showed a half-life of 42 min at 60 degrees C, which was 60 times longer than that of the wild-type enzyme. The thermostable mutant was also more stable with a high concentration of glycerol, which is a necessary condition for in vitro amidation. 相似文献
4.
Prolyl endopeptidases (PEPs) are believed to be involved in the metabolism of neuropeptide hormones (reviewed in Mentlein [1988]). Genes encoding PEPs have been isolated from various species, but their expression patterns during development have not been determined. In this study, we isolated a gene encoding a predicted PEP from the fruitfly Drosophila melanogaster. The gene encodes a predicted 756-amino acid protein having extensive sequence similarity to human PEP. We demonstrated that the Drosophila gene (DPEP) is expressed in a spatially restricted pattern in imaginal discs and the larval brain. Our results suggest a role for DPEP in the regional specification of larval tissues. They also provide a starting point for a genetic analysis of the function of this enzyme during development. 相似文献
5.
Overexpression and characterization of a prolyl endopeptidase from the hyperthermophilic archaeon Pyrococcus furiosus. 总被引:1,自引:0,他引:1
下载免费PDF全文

V J Harwood J D Denson K A Robinson-Bidle H J Schreier 《Journal of bacteriology》1997,179(11):3613-3618
The maltose-regulated mlr-2 gene from the hyperthermophilic archaeon Pyrococcus furiosus having homology to bacterial and eukaryal prolyl endopeptidase (PEPase) was cloned and overexpressed in Escherichia coli. Extracts from recombinant cells were capable of hydrolyzing the PEPase substrate benzyloxycarbonyl-Gly-Pro-p-nitroanilide (ZGPpNA) with a temperature optimum between 85 and 90 degrees C. Denaturing gel electrophoresis of purified PEPase showed that enzyme activity was associated with a 70-kDa protein, which is consistent with that predicted from the mlr-2 sequence. However, an apparent molecular mass of 59 kDa was obtained from gel permeation studies. In addition to ZGPpNA (K(Mapp) of 53 microM), PEPase was capable of hydrolyzing azocasein, although at a low rate. No activity was detected when ZGPpNA was replaced by substrates for carboxypeptidase A and B, chymotrypsin, subtilisin, and neutral endopeptidase. N-[N-(L-3-trans-Carboxirane-2-carbonyl)-L-Leu]-agmatine (E-64) and tosyl-L-Lys chloromethyl ketone did not inhibit PEPase activity. Both phenylmethylsulfonyl fluoride and diprotin A inhibited ZGPpNA cleavage, the latter doing so competitively (K(lapp) of 343 microM). At 100 degrees C, the enzyme displayed some tolerance to sodium dodecyl sulfate treatment. Stability of PEPase over time was dependent on protein concentration; at temperatures above 65 degrees C, dilute samples retained most of their activity after 24 h while the activity of concentrated preparations diminished significantly. This decrease was found to be due, in part, to autoproteolysis. Partially purified PEPase from P. furiosus exhibited the same temperature optimum, molecular weight, and kinetic characteristics as the enzyme overexpressed in E. coli. Extracts from P. furiosus cultures grown in the presence of maltose were approximately sevenfold greater in PEPase activity than those grown without maltose. Activity could not be detected in clarified medium obtained from maltose-grown cultures. We conclude that mlr-2, now called prpA, encodes PEPase; the physiological role of this protease is presently unknown. 相似文献
6.
Huang CH Suen CS Lin CT Chien CH Lee HY Chung KM Tsai TY Jiaang WT Hwang MJ Chen X 《Journal of biochemistry》2011,149(6):685-692
Fibroblast activation protein (FAP) is a prolyl-cleaving endopeptidase proposed as an anti-cancer drug target. It is necessary to define its cleavage-site specificity to facilitate the identification of its in vivo substrates and to understand its biological functions. We found that the previously identified substrate of FAP, α(2)-anti-plasmin, is not a robust substrate in vitro. Instead, an intracellular protein, SPRY2, is cleavable by FAP and more suitable for investigation of its substrate specificity in the context of the full-length globular protein. FAP prefers uncharged residues, including small or bulky hydrophobic amino acids, but not charged amino acids, especially acidic residue at P1', P3 and P4 sites. Molecular modelling analysis shows that the substrate-binding site of FAP is surrounded by multiple tyrosine residues and some negatively charged residues, which may exert least preference for substrates with acidic residues. This provides an explanation why FAP cannot cleave interleukins, which have a glutamate at either P4 or P2', despite their P3-P2-P1 sites being identical to SPRY2 or α-AP. Our study provided new information on FAP cleavage-site specificity, which differs from the data obtained by profiling with a peptide library or with the denatured protein, gelatin, as the substrate. Furthermore, our study suggests that negatively charged residues should be avoided when designing FAP inhibitors. 相似文献
7.
8.
A comparison has been made of the specificity of the mammalian neutral metalloendopeptidase, endopeptidase 24.11, with that of the bacterial neutral metalloendopeptidase thermolysin. A series of synthetic oligopeptides which have previously been studied as substrates for thermolysin and used in computer modeling were examined as substrates for the mammalian enzyme. It was found that P1, P2, and P'3 subsite interactions in the mammalian enzyme, although similar to those found in thermolysin, are less restrictive spatially and are considerably less dependent on hydrophobic interactions. This difference was maximally expressed with the synthetic substrate dansyl-D-alanylglycylnitrophenylalanylglycine which is a substrate for the mammalian enzyme, but not for the bacterial enzyme. A comparison of substrates in the free acid form with their corresponding amides showed that binding to the mammalian enzyme is dependent in part on an ionic interaction between the substrate carboxylate group and the enzyme. Such an ionic interaction was not observed with the bacterial enzyme. 相似文献
9.
S Chevallier P Goeltz P Thibault D Banville J Gagnon 《The Journal of biological chemistry》1992,267(12):8192-8199
A prolyl endopeptidase was purified from Flavobacterium meningosepticum. It was digested with trypsin. Two oligonucleotides, based on tryptic peptide sequences and used in PCR experiments, amplified a 300-base pair (bp) fragment. A 2.4-kilobase EcoRI fragment that hybridized to the 300-bp probe was cloned in lambda ZAP and sequenced from both strands. It contains a reading frame of 2115 bp, encoding the complete protein sequence of 705 amino acids. Ion-spray mass spectrometry experiments demonstrated the presence of an NH2-terminal signal peptide: the periplasmic mature protease is 685 residues in length for a molecular mass of 76784 Da. The prolyl endopeptidase showed no general sequence homology with known protein sequences except with that of porcine brain prolyl endopeptidase. In order to identify the active-site serine, the prolyl endopeptidase was labeled with [3H]diisopropyl fluorophosphate. One labeled peptide was purified and sequenced. The active-site serine was located in position 536 within the sequence GRSNGG. This sequence is different from the active-site sequence of the trypsin (GDSGGP) and subtilisin (GTSMAS) families. 相似文献
10.
11.
Tryptic peptides from hemoglobin (Hb) beta-chains were used as model substrates for limited proteolysis by prolyl endopeptidase (EC 3.4.21.26) from porcine muscle. From the physicochemical and enzymatic properties of prolyl endopeptidase the conditions for routine digestion were established as follows: the molar ratio of enzyme to substrate was 1 to 100, and the reaction was carried out in sodium phosphate buffer (pH 6.4) at 37 degrees C for 4 h. Under these conditions the peptide bonds on the carboxyl terminal sides of proline and alanine residues in the tryptic peptides from Hb beta-chains (with Mr values of less than 2100) were hydrolyzed by the enzyme with the exception of the amino terminal alanyl bond and aminoacyl alanyl bond. In addition, one of five seryl bonds was cleaved by the enzyme. However, the Hb beta-chain itself, Mr 16,600, and its two CNBr-peptides with Mr 10,200 and Mr 6400, respectively, were not hydrolyzed. Under the same conditions a prolyl bond in oxidized B-chains of insulin, Mr 3400, was partially digested, and an alanyl bond was not hydrolyzed. The data indicate that the prolyl endopeptidase is useful for the limited proteolysis of peptides with relative masses of less than 3000 at both prolyl and alanyl bonds. 相似文献
12.
J Rosén B Tomkinson G Pettersson O Zetterqvist 《The Journal of biological chemistry》1991,266(6):3824-3834
The present work describes the detection, purification, and characterization of a serine endopeptidase with preference for a phosphoserine in the P1' position of the substrate. During probing for the enzyme in crude extracts, as well as during its 64,000-fold purification, 32P-labeled guanidovaleryl-Arg-Ala-Ser(P)-isobutyl amide (I) was used to measure the cleavage of the Ala-Ser(P) bond. With this substrate, kcat was 1.7 s-1 and Km was 30 microM at the pH optimum, 7.5. The enzyme was classified as a serine peptidase from its reaction with a set of inhibitors, among which diisopropyl fluorophosphate was effective at low (20 microM) concentration. The endopeptidase showed an Mr of 74,000 under native as well as denaturing and reducing conditions, indicating that the native enzyme consists of only one major polypeptide chain. The molecular size and inhibition profile suggested identity of this enzyme with prolyl endopeptidase (EC 3.4.21.26). This was supported by its activity against specific substrates, such as succinyl-Gly-Pro-Leu-Pro-7-amido-4-methylcoumarin (kcat = 7.2 s-1 and Km = 290 microM), and by the inhibition of the latter activity by I. Compared with the cleavage of 100 microM I, Gly-Val-Leu-Arg-Arg-Ala-Ser-Val-Ala-Gln-Leu, after phosphorylation by cAMP-dependent protein kinase, was cleaved at the Ala-Ser(P) bond at a relative rate of 0.43, while cleavage of the Ala-Ser bond of the unphosphorylated undecapeptide was undetectable, i.e. less than 0.03. The pentapeptide Arg-Arg-Pro-Ser-Val was rapidly cleaved at the Pro-Ser bond (relative rate, 2.2). Still, the cleavage of the Pro-Ser(P) bond of the corresponding phosphorylated pentapeptide was even higher (relative rate, 4.0). These data suggest that phosphorylation of a serine residue in the P1' position of at least a few substrates of prolyl endopeptidase will increase the rate of their cleavage. 相似文献
13.
Animal movement, whether for foraging, mate-seeking, predator avoidance, dispersal, or migration, is a fundamental aspect of ecology that shapes spatial abundance distributions, genetic compositions, and dynamics of populations. A variety of movement models have been used for predicting the effects of natural or human-caused landscape changes, invading species, or other disturbances on local ecology. Here we introduce the flow network—a general modeling framework for population dynamics and movement in a metapopulation representing a network of habitat sites (nodes). Based on the principles of physical transport phenomena such as fluid flow through pipes (Pouiselle’s Law) and analogously, the flow of electric current across a circuit (Ohm’s Law), the flow network provides a novel way of modeling movement, where flow rates are functions of relative node pressures and the resistance to movement between them. Flow networks offer the flexibility of incorporating abiotic and biotic conditions that affect either pressures, resistance, or both. To illustrate an application of the flow network, we present a theoretical invasion scenario. We consider the effects of spatial structure on the speed of invasion by varying the spatial regularity of node arrangement. In the context of invasion, we model management actions targeting nodes or edges, and consider the effects on speed of invasion, node occupation, and total abundance. The flow network approach offers the flexibility to incorporate spatial heterogeneity in both rates of flow and site pressures and offers an intuitive approach to connecting population dynamics and landscape features to model movement. 相似文献
14.
T Benchetrit M C Fournié-Zaluski B P Roques 《Biochemical and biophysical research communications》1987,147(3):1034-1040
The inhibitory potency of separate enantiomers of thiorphan and retrothiorphan has shown that several particularities of the active site of thermolysin are also present in the neutral endopeptidase 24.11, "enkephalinase", such as its ability: i) to recognize a retroamide bond as well as a standard amide bond, ii) to interact similarly with residues in P1' position of either R or S configuration in the thiorphan series but contrastingly to discriminate between the R and S isomers in the retrothiorphan series. These four inhibitors were modellized in the thermolysin active site and their spatial arrangement compared with that of a thiol inhibitor co-crystallized with thermolysin. In all cases, the essential interactions involved in the stabilization of the bound inhibitor were conserved. However, the bound (R) retrothiorphan displayed unfavorable intramolecular contacts, accounting for its lower inhibitory potency for the two metallopeptidases. 相似文献
15.
In relating genotypes to fitness, models of adaptation need to both be computationally tractable and qualitatively match observed data. One reason that tractability is not a trivial problem comes from a combinatoric problem whereby no matter in what order a set of mutations occurs, it must yield the same fitness. We refer to this as the bookkeeping problem. Because of their commutative property, the simple additive and multiplicative models naturally solve the bookkeeping problem. However, the fitness trajectories and epistatic patterns they predict are inconsistent with the patterns commonly observed in experimental evolution. This motivates us to propose a new and equally simple model that we call stickbreaking. Under the stickbreaking model, the intrinsic fitness effects of mutations scale by the distance of the current background to a hypothesized boundary. We use simulations and theoretical analyses to explore the basic properties of the stickbreaking model such as fitness trajectories, the distribution of fitness achieved, and epistasis. Stickbreaking is compared to the additive and multiplicative models. We conclude that the stickbreaking model is qualitatively consistent with several commonly observed patterns of adaptive evolution. 相似文献
16.
In vitro molecular evolution is regarded as a hill-climbing on a fitness landscape in sequence space, where the 'fitness' is a quantitative measure of a certain physicochemical property of a biopolymer. We analyzed a 'cross-section' of the enzymatic activity landscape of dihydrofolate reductase (DHFR) by using a method of analysis of a fitness landscape. We limited the sequence space of interest to the five-dimensional sequence space, where the coordinate corresponds to the 1st, 16th, 20th, 42nd and 92nd site in the DHFR sequence. Thirty six mutants mapped into the limited sequence space were taken in the analysis. As a result, the cross-section is of the rough Mt Fuji type based on the mutational additivity. The ratio of the mean slope to the roughness is 2.8 and the Z-score of the original ratio against a distribution of random references is 7.0, which indicates a large statistical significance. The existence of such a cross-section was discussed in terms of the occurrence probability of sets of five sites distantly separated from each other on the DHFR 3D structure. Our results support the effectiveness of the evolution strategy which exploits the accumulation of advantageous single point mutations in such a cross-section. 相似文献
17.
Nevozhay D Adams RM Van Itallie E Bennett MR Balázsi G 《PLoS computational biology》2012,8(4):e1002480
Gene expression actualizes the organismal phenotypes encoded within the genome in an environment-dependent manner. Among all encoded phenotypes, cell population growth rate (fitness) is perhaps the most important, since it determines how well-adapted a genotype is in various environments. Traditional biological measurement techniques have revealed the connection between the environment and fitness based on the gene expression mean. Yet, recently it became clear that cells with identical genomes exposed to the same environment can differ dramatically from the population average in their gene expression and division rate (individual fitness). For cell populations with bimodal gene expression, this difference is particularly pronounced, and may involve stochastic transitions between two cellular states that form distinct sub-populations. Currently it remains unclear how a cell population's growth rate and its subpopulation fractions emerge from the molecular-level kinetics of gene networks and the division rates of single cells. To address this question we developed and quantitatively characterized an inducible, bistable synthetic gene circuit controlling the expression of a bifunctional antibiotic resistance gene in Saccharomyces cerevisiae. Following fitness and fluorescence measurements in two distinct environments (inducer alone and antibiotic alone), we applied a computational approach to predict cell population fitness and subpopulation fractions in the combination of these environments based on stochastic cellular movement in gene expression space and fitness space. We found that knowing the fitness and nongenetic (cellular) memory associated with specific gene expression states were necessary for predicting the overall fitness of cell populations in combined environments. We validated these predictions experimentally and identified environmental conditions that defined a "sweet spot" of drug resistance. These findings may provide a roadmap for connecting the molecular-level kinetics of gene networks to cell population fitness in well-defined environments, and may have important implications for phenotypic variability of drug resistance in natural settings. 相似文献
18.
Meng-Xiang Wang Chan Zhong Qiu-Feng Cai Guang-Ming Liu Ling Zhang Kenji Hara Wen-Jin Su Min-Jie Cao 《Process Biochemistry》2012,47(12):2211-2218
A prolyl endopeptidase (PEP) was purified to homogeneity from the skeletal muscle of common carp using a procedure involving ammonium sulfate fractionation and column chromatography involving DEAE-Sephacel, Phenyl-Sepharose, DEAE-Sepharose Fast Flow, and hydroxyapatite. The molecular weight of the PEP was 82 kDa as determined by SDS-PAGE. Using Suc-Gly-Pro-MCA as a substrate, the optimal pH and temperature of the purified enzyme were pH 6.0 and 35 °C, respectively, and the Km and kcat were 8.33 μM and 1.71 S?1, respectively. The activity of the PEP was inhibited by SUAM-14746, a specific inhibitor of prolyl endopeptidases, and was partially inhibited by the serine proteinase inhibitors PMSF and Pefabloc SC. According to peptide mass fingerprinting, 12 peptide fragments with a total of 134 amino acid residues were obtained, which were highly identical to prolyl endopeptidases from zebrafish (Danio rerio) and sponge (Amphimedon queenslandica), confirming the purified enzyme was a prolyl endopeptidase. Our present study for the first time reported the existence of a prolyl endopeptidase in fish muscle. 相似文献
19.
H. -A. Eggers-Schumacher G. Forkmann K. Wöhrmann 《TAG. Theoretical and applied genetics. Theoretische und angewandte Genetik》1977,49(4):187-195
Summary A new fitness function is proposed. It is a balance function which results from a combination of a turnover function defining the expression of favourable genetic factors and a cost function describing the costs of metabolism. A favourable genetic factor is defined as the heterozygote state at a locus. Five loci with two alleles each are considered. In computer simulations various kinds of this function are investigated. Their influence on the genetic composition of a population in equilibrium, the number and types of equilibria and the mean fitness are reported. Free recombination and tight linkage are considered. 相似文献
20.
The NK model of rugged fitness landscapes and its application to maturation of the immune response 总被引:2,自引:0,他引:2
Adaptive evolution is, to a large extent, a complex combinatorial optimization process. Such processes can be characterized as "uphill walks on rugged fitness landscapes". Concrete examples of fitness landscapes include the distribution of any specific functional property such as the capacity to catalyze a specific reaction, or bind a specific ligand, in "protein space". In particular, the property might be the affinity of all possible antibody molecules for a specific antigenic determinant. That affinity landscape presumably plays a critical role in maturation of the immune response. In this process, hypermutation and clonal selection act to select antibody V region mutant variants with successively higher affinity for the immunizing antigen. The actual statistical structure of affinity landscapes, although knowable, is currently unknown. Here, we analyze a class of mathematical models we call NK models. We show that these models capture significant features of the maturation of the immune response, which is currently thought to share features with general protein evolution. The NK models have the important property that, as the parameter K increases, the "ruggedness" of the NK landscape varies from a single peaked "Fujiyama" landscape to a multi-peaked "badlands" landscape. Walks to local optima on such landscapes become shorter as K increases. This fact allows us to choose a value of K that corresponds to the experimentally observed number of mutational "steps", 6-8, taken as an antibody sequence matures. If the mature antibody is taken to correspond to a local optimum in the model, tuning the model requires that K be about 40, implying that the functional contribution of each amino acid in the V region is affected by about 40 others. Given this value of K, the model then predicts several features of "antibody space" that are in qualitative agreement with experiment: (1) The fraction of fitter variants of an initial "roughed in" germ line antibody amplified by clonal selection is about 1-2%. (2) Mutations at some sites of the mature antibody hardly affect antibody function at all, but mutations at other sites dramatically decrease function. (3) The same "roughed in" antibody sequence can "walk" to many mature antibody sequences. (4) Many adaptive walks can end on the same local optimum. (5) Comparison of different mature sequences derived from the same initial V region shows evolutionary hot spots and parallel mutations. All these predictions are open to detailed testing by obtaining monoclonal antibodies early in the immune response and carrying out in vitro mutagenesis and adaptive hill climbing with respect to affinity for the immunizing antigen. 相似文献