Multiple protein sequence alignment from tertiary structure comparison: assignment of global and residue confidence levels.

An algorithm is presented for the accurate and rapid generation of multiple protein sequence alignments from tertiary structure comparisons. A preliminary multiple sequence alignment is performed using sequence information, which then determines an initial superposition of the structures. A structure comparison algorithm is applied to all pairs of proteins in the superimposed set and a similarity tree calculated. Multiple sequence alignments are then generated by following the tree from the branches to the root. At each branchpoint of the tree, a structure-based sequence alignment and coordinate transformations are output, with the multiple alignment of all structures output at the root. The algorithm encoded in STAMP (STructural Alignment of Multiple Proteins) is shown to give alignments in good agreement with published structural accounts within the dehydrogenase fold domains, globins, and serine proteinases. In order to reduce the need for visual verification, two similarity indices are introduced to determine the quality of each generated structural alignment. Sc quantifies the global structural similarity between pairs or groups of proteins, whereas Pij' provides a normalized measure of the confidence in the alignment of each residue. STAMP alignments have the quality of each alignment characterized by Sc and Pij' values and thus provide a reproducible resource for studies of residue conservation within structural motifs.     

Recurrent alpha beta loop structures in TIM barrel motifs show a distinct pattern of conserved structural features.

A systematic survey of seven parallel alpha/beta barrel protein domains, based on exhaustive structural comparisons, reveals that a sizable proportion of the alpha beta loops in these proteins--20 out of a total of 49--belong to either one of two loop types previously described by Thornton and co-workers. Six loops are of the alpha beta 1 type, with one residue between the alpha-helix and beta-strand, and 13 are of the alpha beta 3 type, with three residues between the helix and the strand. Protein fragments embedding the identified loops, and termed alpha beta connections since they contain parts of the flanking helix and strand, have been analyzed in detail revealing that each type of connection has a distinct set of conserved structural features. The orientation of the beta-strand relative to the helix and loop portions is different owing to a very localized difference in backbone conformation. In alpha beta 1 connections, the chain enters the beta-strand via a residue adopting an extended conformation, while in alpha beta 3 it does so via a residue in a near alpha-helical conformation. Other conserved structural features include distinct patterns of side chain orientation relative to the beta-sheet surface and of main chain H-bonds in the loop and the beta-strand moieties. Significant differences also occur in packing interactions of conserved hydrophobic residues situated in the last turn of the helix. Yet the alpha-helix surface of both types of connections adopts similar orientations relative to the barrel sheet surface. Our results suggest furthermore that conserved hydrophobic residues along the sequence of the connections, may be correlated more with specific patterns of interactions made with neighboring helices and sheet strands than with helix/strand packing within the connection itself. A number of intriguing observations are also made on the distribution of the identified alpha beta 1 and alpha beta 3 loops within the alpha/beta-barrel motifs. They often occur adjacent to each other; alpha beta 3 loops invariably involve even numbered beta-strands, while alpha beta 1 loops involve preferentially odd beta-strands; all the analyzed proteins contain at least one alpha beta 3 loop in the first half of the eightfold alpha/beta barrel. Possible origins of all these observations, and their relevance to the stability and folding of parallel alpha/beta barrel motifs are discussed.     

Three-dimensional structure of the quinoprotein methylamine dehydrogenase from Paracoccus denitrificans determined by molecular replacement at 2.8 A resolution.

The three-dimensional structure of the quinoprotein methylamine dehydrogenase from Paracoccus dentrificans (PD-MADH) has been determined at 2.8 A resolution by the molecular replacement method combined with map averaging procedures, using data collected from an area detector. The structure of methylamine dehydrogenase from Thio-bacillus versutus, which contains an "X-ray" sequence, was used as the starting search model. MADH consists of 2 heavy (H) and 2 light (L) subunits related by a molecular 2-fold axis. The H subunit is folded into seven four-stranded beta segments, forming a disk-shaped structure, arranged with pseudo-7-fold symmetry. A 31-residue elongated tail exists at the N-terminus of the H subunit in MADH from T. versutus but is partially digested in this crystal form of MADH from P. denitrificans, leaving the H subunit about 18 residues shorter. Each L subunit contains 127 residues arranged into 10 beta-strands connected by turns. The active site of the enzyme is located in the L subunit and is accessible via a hydrophobic channel between the H and L subunits. The redox cofactor of MADH, tryptophan tryptophylquinone is highly unusual. It is formed from two covalently linked tryptophan side chains at positions 57 and 107 of the L subunit, one of which contains an orthoquinone.     

Application of a directed conformational search for generating 3-D coordinates for protein structures from alpha-carbon coordinates.

A directed conformational search algorithm using the program CONGEN (ref. 3), which samples backbone conformers, is described. The search technique uses information from the partially built structures to direct the search process and is tested on the problem of generating a full set of backbone Cartesian coordinates given only alpha-carbon coordinates. The method has been tested on six proteins of known structure, varying in size and classification, and was able to generate the original backbone coordinates with RMSs ranging from 0.30-0.87A for the alpha-carbons and 0.5-0.99A RMSs for the backbone atoms. Cis peptide linkages were also correctly identified. The procedure was also applied to two proteins available with only alpha-carbon coordinates in the Brookhaven Protein Data Bank; thioredoxin (SRX) and triacylglycerol acylhydrolase (TGL). All-atom models are proposed for the backbone of both these proteins. In addition, the technique was applied to randomized coordinates of flavodoxin to assess the effects of irregularities in the data on the final RMS. This study represents the first time a deterministic conformational search was used on such a large scale.     

Developmental changes in barrier web structure under different levels of predation risk inNephila clavipes (Araneae: Tetragnathidae)

Individuals of the orb-weaving spider Nephila clavipesbuild complex webs with a region used for prey capture, the orb, and tangle webs opposite either face, the barrier webs. Barrier webs have been hypothesized to serve a variety of functions, including predator defense, and the primary function of the barrier web should be reflected in the relative size of the barrier to the orb under varying conditions of foraging success and predation risk. To investigate the effects of predation pressure and foraging success on barrier web structure, I conducted a comparative study in three disjunct populations that differed in predation risk and foraging success. Although both the orb web and the barrier webs are silk, there was no indication of a foraging-defense trade-off. Barrier web structure did not change during seasonal shifts in orb web size related to changes in preycapture rate, and barrier web silk density and orb radius were positively correlated. The hypothesis that the construction of barrier webs is in part a response to predation pressure was supported. Barrier webs do deflect attacks by some predators, and barrier webs built by small spiders, suffering frequent predation attempts, had a higher silk density than barrier webs built by larger individuals. Additionally, barrier web complexity decreased at a later age in areas with higher predation risk.     

Size dependency of gametophytes decay inAthyrium brevifrons Nakai during spring desiccation

Gametophyte populations inAthyrium brevifrons were analysed with respect to population size and surviving area (%) of individual thalli in a transplant garden at Sapporo during 5–26 April 1983, to study the safe-microsite for gametophyte establishment in nature. Spores dispersed in August 1982 germinated and grew into thalli of various widths (<10 mm); 10.3% of the thalli matured by early October 1982. Maturation was attained by gametophytes of width 4–7 mm. The number of gametophytes gradually decreased with increasing width. By April 1983, 20.5% of total gametophytes were mature with a mode of 5–6 mm in width. The relative number of gametophytes with surviving area of 2–20% increased and that of 85–100% decreased in accordance with collection days delayed until after snow-melt. Surviving area (%) on gametophyte of all widths decreased with decreasing soil moisture contents. In particular, immature gametophytes of 2–4 mm width showed a significant correlation (P<0.01) between soil moisture content and relative number of gametophytes with 0–20% surviving area and mean surviving area (%) of every width of thalli. The spring desiccation might be a factor that reduces or limits gametophyte populations in nature.     

The gap theory in forest dynamics

Since the late 1970s, ecologists interested in forest dynamics have focused their attention to the responses of individuals, populations and communities to “gaps” which are openings created in the forest canopy. This review intended to introduce some collective knowledge on major subjects of the gap theory in forest dynamics, in relation to gap-disturbance regimes, tree regeneration behaviour and community structure.     

Fish size and habitat depth relationships in headwater streams

Summary Surveys of 262 pools in 3 small streams in eastern Tennessee demonstrated a strong positive relationship between pool depth and the size of the largest fish within a pool (P<0.001). Similarly, the largest colonizers of newly-created deep pools were larger than the colonizers of shallow pools. We explored the role of predation risk in contributing to the bigger fish — deeper habitat pattern, which has been noted by others, by conducting five manipulative field experiments in two streams. Three experiments used stoneroller minnows (Campostoma anomalum); one used creek chubs (Semotilus atromaculatus); and one used striped shiners (Notropis chrysocephalus). The stoneroller experiments showed that survival of fish approximately 100 mm in total length (TL) was much lower in shallow pools (10 cm deep) than in deep (40 cm maximum) pools (19% versus 80% survival over 12 d in one experiment) and added cover markedly increased stoneroller survival in shallow pools (from 49% to 96% in an 11-d experiment). The creek chub experiment showed that, as for stonerollers, pool depth markedly influenced survival: the chubs survived an average of 4.9 d in shallow pools and >10.8 d in deep pools. In the striped shiner experiment in shallow artificial streamside troughs, no individuals 75–100 mm TL survived as long as 13 d, where-as smaller (20–25 mm) fish had 100% survival over 13 d. The results of the experiments show that predation risk from wading/diving animals (e.g., herons and raccoons) is much higher for larger fishes in shallow water than for these fishes in deeper water or for smaller fish in shallow water. We discuss the role of predation risk from two sources (piscivorous fish, which are more effective in deeper habitats, and diving/wading predators, which are more effective in shallow habitats) in contributing to the bigger fish — deeper habitat pattern in streams.     

Current perspectives on plasmodesmata: structure and function

Recent studies on plasmodesmata have shown that these important intercellular passages for communication and transport are much more sophisticated in both structure and regulatory abilities than previously imagined. A complex, but not well understood, substructure has been revealed by a variety of increasingly reliable ultrastructural techniques. Proteinaceous particles are seen within the cytoplasmic sleeve surrounding the desmotubule. Dye-coupling studies have provided experimental evidence for the physical pathway of solute movement, supporting conclusions about substructural dimensions within plasmodesmata drawn from the ultrastructural studies. Calcium has been identified as a major factor in the regulation of intercellular communication via plasmodesmata. Evidence from studies on virus movement through plasmodesmata suggests a direct interaction between virallycoded movement proteins and plasmodesmata in the systemic spread of many viruses. There is increasing evidence, albeit indirect, that in some plant species phloem loading may involve transport of photoassimilate entirely within the symplast from mesophyll cells to the sieve element-companion cell complexes of minor veins.     

Towards an explanation of carboxylation/oxygenation bifunctionality in Rubisco. Transition structure for the carboxylation reaction of 2,3,4-pentanetriol

We construct a theoretical model of the transition structure for the carboxylation reaction of ribulose-1,5-biphosphate catalyzed by Rubisco. This is a first-order saddle point on the energy hypersurface for the nucleophilic attack of carbon dioxide on CH₃-(CHOH)₃-CH₃ at the C2 center.Ab initio analytical gradients methods at a 4-31G basis set level are used.The carbon framework and oxygens of the stationary structure superpose with the corresponding atoms of 2-carboxyarabinitol-1,5-biphosphate, which is a transition state analog that has recently been highly refined with X-ray methods. The hydroxyl group in C3 iscis to the C2 oxygen. The C3 center is somewhat pyramidized, the dienol O2-C2-C3-O3 is not planar.The geometry of the transition state allows for simple explanations of both the enolization of Rubisco's substrate ribulose-1,5-biphosphate, O₃PO-CH2-CO-(CHOH)₂-CH₂-OPO₃ and oxygenation reaction. The former is due to the pyramidal deformation at C3 and out of plane of O2-C2-C3-O3 framework: the enoliation is intramolecular and is probably enhanced by proton tunnelling. The latter is related with the fact that a rotation around an ethylene-like bond brings the triplet state down in energy. The reactive skeleton has a stationary geometry in the triplet state not very different from the one obtained in the global transition structure. There, the triplet is only 9 kcal/mol above the singlet. The spin densities at C2 and C3 centers clearly indicate the place where oxygenation will take place.