Metapopulation dynamics with quasi-local competition

Stepping-stone models for the ecological dynamics of metapopulations are often used to address general questions about the effects of spatial structure on the nature and complexity of population fluctuations. Such models describe an ensemble of local and spatially isolated habitat patches that are connected through dispersal. Reproduction and hence the dynamics in a given local population depend on the density of that local population, and a fraction of every local population disperses to neighboring patches. In such models, interesting dynamic phenomena, e.g. the persistence of locally unstable predator-prey interactions, are only observed if the local dynamics in an isolated patch exhibit non-equilibrium behavior. Therefore, the scope of these models is limited. Here we extend these models by making the biologically plausible assumption that reproductive success in a given local habitat not only depends on the density of the local population living in that habitat, but also on the densities of neighboring local populations. This would occur if competition for resources occurs between neighboring populations, e.g. due to foraging in neighboring habitats. With this assumption of quasi-local competition the dynamics of the model change completely. The main difference is that even if the dynamics of the local populations have a stable equilibrium in isolation, the spatially uniform equilibrium in which all local populations are at their carrying capacity becomes unstable if the strength of quasi-local competition reaches a critical level, which can be calculated analytically. In this case the metapopulation reaches a new stable state, which is, however, not spatially uniform anymore and instead results in an irregular spatial pattern of local population abundance. For large metapopulations, a huge number of different, spatially non-uniform equilibrium states coexist as attractors of the metapopulation dynamics, so that the final state of the system depends critically on the initial conditions. The existence of a large number of attractors has important consequences when environmental noise is introduced into the model. Then the metapopulation performs a random walk in the space of all attractors. This leads to large and complicated population fluctuations whose power spectrum obeys a red-shifted power law. Our theory reiterates the potential importance of spatial structure for ecological processes and proposes new mechanisms for the emergence of non-uniform spatial patterns of abundance and for the persistence of complicated temporal population fluctuations.     

ATP synthases: insights into their motor functions from sequence and structural analyses

ATP synthases are motor complexes comprised of F₀ and F₁ parts that couple the proton gradient across the membrane to the synthesis of ATP by rotary catalysis. Although a great deal of information has been accumulated regarding the structure and function of ATP synthases, their motor functions are not fully understood. For this reason, we performed the alignments and analyses of the protein sequences comprising the core of the ATP synthase motor complex, and examined carefully the locations of the conserved residues in the subunit structures of ATP synthases. A summary of the findings from this bioinformatic study is as follows. First, we found that four conserved regions in the sequence of subunit are clustered into three patches in its structure. The interactions of these conserved patches with the and subunits are likely to be critical for energy coupling and catalytic activity of the ATP synthase. Second, we located a four-residue cluster at the N-terminal domain of mitochondrial OSCP or bacterial (or chloroplast) subunit which may be critical for the binding of these subunits to F₁. Third, from the localizations of conserved residues in the subunits comprising the rotors of ATP synthases, we suggest that the conserved interaction site at the interface of subunit c and (mitochondria) or (bacteria and chloroplasts) may be important for connecting the rotor of F₁ to the rotor of F₀. Finally, we found the sequence of mitochondrial subunit b to be highly conserved, significantly longer than bacterial subunit b, and to contain a shorter dimerization domain than that of the bacterial protein. It is suggested that the different properties of mitochondrial subunit b may be necessary for interaction with other proteins, e.g., the supernumerary subunits.     

Protein family alignment annotation

For bioscientists studying protein structure and function, the Protein Family Alignment Annotation Tool (Pfaat) is a useful and simple program for annotating collections of proteins. This open-source software includes methods for viewing and aligning protein families, and for annotating sequence structure and residues with known functions. It offers new options to aid the study of proteins, and an extensible annotation tool for bioinformatics developers.     

Observation of a distinct transition in the mode of interconversion of ring pucker conformers in non-crystalline d-ribose-2'-d from 2H NMR spin-alignment

Internal motions of d-ribose selectively ²H-labeled at the 2 position were measured using solid state ²H NMR experiments. A sample of d-ribose-2 -d was prepared in a hydrated, non-crystalline state to eliminate effects of crystal-packing. Between temperatures of –74 and –60°C the C2–H2 bond was observed to undergo two kinds of motions which were similar to those of C2–H2/H2 found previously in crystalline deoxythymidine (Hiyama et al. (1989) J. Am. Chem. Soc., 111, 8609–8613): (1) Nanosecond motion of small angular displacement with an apparent activation energy of 3.6 ± 0.7 kcal mol^–1, and (2) millisecond to microsecond motion of large amplitude with an apparent activation energy 4 kcal mol^–1. At –74°C, the slow, large-amplitude motion was best characterized as a two-site jump with a correlation time on the millisecond time scale, whereas at –60°C it was diffusive on the microsecond time scale. The slow, large-amplitude motions of the C2–H2 bond are most likely from interconversions between C2-endo and C3-endo by way of the O4-endo conformation, whereas the fast, small-amplitude motions are probably librations of the C2–H2 bond within the C2-endo and C3-endo potential energy minima.     

Equilibrium structure and stability in a frequency-dependent,two-population diploid model

We investigate the equilibrium structure for an evolutionary genetic model in discrete time involving two monoecious populations subject to intraspecific and interspecific random pairwise interactions. A characterization for local stability of an equilibrium is found, related to the proximity of this equilibrium with evolutionarily stable strategies (ESS). This extends to a multi-population framework a principle initially proposed for single populations, which states that the mean population strategy at a locally stable equilibrium is as close as possible to an ESS.     

The role of pattern databases in sequence analysis

In the wake of the numerous now-fruitful genome projects, we are entering an era rich in biological data. The field of bioinformatics is poised to exploit this information in increasingly powerful ways, but the abundance and growing complexity both of the data and of the tools and resources required to analyse them are threatening to overwhelm us. Databases and their search tools are now an essential part of the research environment. However, the rate of sequence generation and the haphazard proliferation of databases have made it difficult to keep pace with developments. In an age of information overload, researchers want rapid, easy-to-use, reliable tools for functional characterisation of newly determined sequences. But what are those tools? How do we access them? Which should we use? This review focuses on a particular type of database that is increasingly used in the task of routine sequence analysis--the so-called pattern database. The paper aims to provide an overview of the current status of pattern databases in common use, outlining the methods behind them and giving pointers on their diagnostic strengths and weaknesses.     

The design of a hyperstable mutant of the Abp1p SH3 domain by sequence alignment analysis

We have characterized the thermodynamic stability of the SH3 domain from the Saccharomyces cerevisiae Abp1p protein and found it to be relatively low compared to most other SH3 domains, with a Tm of 60 degrees C and a deltaGu of 3.08 kcal/mol. Analysis of a large alignment of SH3 domains led to the identification of atypical residues at eight positions in the wild-type Abp1p SH3 domain sequence that were subsequently replaced by the residue seen most frequently at that position in the alignment. Three of the eight mutants constructed in this way displayed increases in Tm ranging from 8 to 15 degrees C with concomitant increases in deltaGu of up to 1.4 kcal/mol. The effects of these substitutions on folding thermodynamics and kinetics were entirely additive, and a mutant containing all three was dramatically stabilized with a Tm greater than 90 degrees C and a deltaGu more than double that of the wild-type domain. The folding rate of this hyperstable mutant was 10-fold faster than wild-type, while its unfolding rate was fivefold slower. All of the stabilized mutants were still able to bind a target peptide with wild-type affinity. We have analyzed the stabilizing amino acid substitutions isolated in this study and several other similar sequence alignment based studies. In approximately 25% of cases, increased stability can be explained by enhanced propensity of the substituted residue for the local backbone conformation at the mutagenized site.     

Prediction of amino acid sequence from structure

We have developed a method for the prediction of an amino acid sequence that is compatible with a three-dimensional backbone structure. Using only a backbone structure of a protein as input, the algorithm is capable of designing sequences that closely resemble natural members of the protein family to which the template structure belongs. In general, the predicted sequences are shown to have multiple sequence profile scores that are dramatically higher than those of random sequences, and sometimes better than some of the natural sequences that make up the superfamily. As anticipated, highly conserved but poorly predicted residues are often those that contribute to the functional rather than structural properties of the protein. Overall, our analysis suggests that statistical profile scores of designed sequences are a novel and valuable figure of merit for assessing and improving protein design algorithms.     

Evaluation of PSI-BLAST alignment accuracy in comparison to structural alignments

The PSI-BLAST algorithm has been acknowledged as one of the most powerful tools for detecting remote evolutionary relationships by sequence considerations only. This has been demonstrated by its ability to recognize remote structural homologues and by the greatest coverage it enables in annotation of a complete genome. Although recognizing the correct fold of a sequence is of major importance, the accuracy of the alignment is crucial for the success of modeling one sequence by the structure of its remote homologue. Here we assess the accuracy of PSI-BLAST alignments on a stringent database of 123 structurally similar, sequence-dissimilar pairs of proteins, by comparing them to the alignments defined on a structural basis. Each protein sequence is compared to a nonredundant database of the protein sequences by PSI-BLAST. Whenever a pair member detects its pair-mate, the positions that are aligned both in the sequential and structural alignments are determined, and the alignment sensitivity is expressed as the percentage of these positions out of the structural alignment. Fifty-two sequences detected their pair-mates (for 16 pairs the success was bi-directional when either pair member was used as a query). The average percentage of correctly aligned residues per structural alignment was 43.5+/-2.2%. Other properties of the alignments were also examined, such as the sensitivity vs. specificity and the change in these parameters over consecutive iterations. Notably, there is an improvement in alignment sensitivity over consecutive iterations, reaching an average of 50.9+/-2.5% within the five iterations tested in the current study.     

The identification of conserved interactions within the SH3 domain by alignment of sequences and structures

The SH3 domain, comprised of approximately 60 residues, is found within a wide variety of proteins, and is a mediator of protein-protein interactions. Due to the large number of SH3 domain sequences and structures in the databases, this domain provides one of the best available systems for the examination of sequence and structural conservation within a protein family. In this study, a large and diverse alignment of SH3 domain sequences was constructed, and the pattern of conservation within this alignment was compared to conserved structural features, as deduced from analysis of eighteen different SH3 domain structures. Seventeen SH3 domain structures solved in the presence of bound peptide were also examined to identify positions that are consistently most important in mediating the peptide-binding function of this domain. Although residues at the two most conserved positions in the alignment are directly involved in peptide binding, residues at most other conserved positions play structural roles, such as stabilizing turns or comprising the hydrophobic core. Surprisingly, several highly conserved side-chain to main-chain hydrogen bonds were observed in the functionally crucial RT-Src loop between residues with little direct involvement in peptide binding. These hydrogen bonds may be important for maintaining this region in the precise conformation necessary for specific peptide recognition. In addition, a previously unrecognized yet highly conserved beta-bulge was identified in the second beta-strand of the domain, which appears to provide a necessary kink in this strand, allowing it to hydrogen bond to both sheets comprising the fold.