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.     

Detection of native-like models for amino acid sequences of unknown three-dimensional structure in a data base of known protein conformations.

We present an approach which can be used to identify native-like folds in a data base of protein conformations in the absence of any sequence homology to proteins in the data base. The method is based on a knowledge-based force field derived from a set of known protein conformations. A given sequence is mounted on all conformations in the data base and the associated energies are calculated. Using several conformations and sequences from the globin family we show that the native conformation is identified correctly. In fact the resolution of the force field is high enough to discriminate between a native fold and several closely related conformations. We then apply the procedure to several globins of known sequence but unknown three dimensional structure. The homology of these sequences to globins of known structures in the data base ranges from 49 to 17%. With one exception we find that for all globin sequences one of the known globin folds is identified as the most favorable conformation. These results are obtained using a force field derived from a data base devoid of globins of known structure. We briefly discuss useful applications in protein structural research and future development of our approach.     

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.     

Phosphorus efficiency of plants

Föhse et al. (1988) have shown that P influx per unit root length in seven plant species growing in a low-P soil varied from 0.6×10^-14 to 4.8×10^-14 mol cm^-1s^-1. The objective of this work was to investigate the reasons for these differences. No correlation was found between P influx and root radius, root hairs, cation-anion balance and Ca uptake. However, when root hairs were included in mathematical model calculations, the differences of P influx could be accounted for. These calculations have shown that in soils low in available P, contribution to P uptake by root hairs was up to 90% of total uptake. The large contribution of root hairs to P uptake was partly due to their surface area, which was similar to that of the root cylinder. However, the main reason for the high P uptake efficiency of root hairs was their small radius (approx. 5×10^-4 cm) and their perpendicular growth into the soil from the root axis. Because of the small radius compared to root axes, P concentration at root hair surfaces decreased at a slower pace and therefore P influx remained higher. Under these conditions higher I_max (maximum influx) or smaller K_m values (Michaelis constant) increased P influx. The main reasons for differences found in P influx among species were the size of I_max and the number and length of root hairs. In a soil low in available P, plant species having more root hairs were able to satisfy a higher proportion of their P demand required for maximum growth.     

Theory of chaperonin action: inertial model for enhancement of prokaryotic Rubisco assembly.

We have performed a computational simulation of the aggregation and chaperonin-dependent reconstitution of dimeric prokaryotic ribulose bisphosphate carboxylase/oxygenase (Rubisco), based on the data of P. Goloubinoff et al. (1989, Nature 342, 884-889) and P. V. Viitanen et al. (1990, Biochemistry 29, 5665-5671). The aggregation is simulated by a set of 12 differential equations representing the aggregation of the Rubisco folding intermediate, Rubisco-I, with itself and with aggregates of Rubisco-I, leading up to dodecamers. Four rate constants, applying to forward or reverse steps in the aggregation process, were included. Optimal values for these constants were determined using the ellipsoid algorithm as implemented by one of us (Ecker, J.G. & Kupferschmid, M., 1988, Introduction to Operations Research, Wiley, New York, pp. 315-322). Intensive exploration of simpler aggregation models did not identify an alternative that could simulate the data as well as this one. The activity of the chaperonin in this system was simulated by using this aggregation model, combined with a model similar to that proposed by Goloubinoff et al. (1989). The model assumes that the chaperonin can bind the folding intermediate rapidly, and that the chaperonin complex releases the Rubisco molecule slowly, permitting time for its spontaneous folding while interacting with the chaperonin. This is followed by self-association of the folded Rubisco monomer to yield the active dimeric Rubisco. A modification of the model that simulates temperature effects was also constructed. The most important results we obtained indicate that the chaperonin-dependent reconstitution of Rubisco can be simulated adequately without invoking any catalysis of folding by the chaperonin. In addition, the simulations predict values for the association rate constant of Rubisco-I with the chaperonin, and other variables, that are subject to experimental verification.     

光钳操纵生物活体的动态监测技术的研究

采用摄像、录像和视屏监控系统及显色偏振装置与光钳系统耦合,从空间分辨、色分辨和时间分辨多方面改善系统品质,实现了光钳捕获与操纵生物活体的动态监测、实时记录、资料保存和屏幕再现的功能,并能测量光钳操纵细胞的位移量和由此计算操纵速度,提高了光钳的自我调整和光钳操纵细胞的精细度。本研究为激光光钳技术在细胞工程等方面的应用研究提供了行之有效的技术手段。     

Bacterioplankton in the acidified Lake Gårdsjön

Total bacterial numbers in different strata of lake water and in inlet and outlet streams have been recorded during a yearly cycle. A calculated mean cell volume of 0.342 µm² has then been used to estimate bacterial biomass in the lake. Change of biomass during the year was substantial and the range was from about 0.1 g · m^–3 to about 1.0–1.2 g · m^–3. The seasonal development included a spring-early summer increase followed by a decrease to the minimum in July–August. Correlation between epi- and hypolimnion was high and in both strata two dominant autumn peaks in biomass appeared. With the exception of the last autumn peak the development of bacterial biomass was closely related to development of phytoplankton biomass and production.     

Similarity of EPR Signal IIf rise and P-680+ decay kinetics in Tris-washed chloroplast Photosystem II preparations as a function of pH

The rise time, of Signal II_f and the decay time of P-680⁺ have been measured kinetically as a function of pH by using EPR. The Photosystem II-enriched preparations which were used as samples were derived from spinach chloroplasts, and they evolved oxygen before Tris washing. The onset kinetics of Signal II_f are in agreement, within experimental error, with the fast component of the decay of an EPR signal attributable to P-680⁺. The signal II_f rise kinetics also show good agreement with published values of the pH dependence of the decay of P-680⁺ measured optically (Conjeaud, H. and Mathis, P. (1980) Biochim. Biophys. Acta 590, 353–359). These results are consistent with a model where the species Z (or D₁) responsible for Signal II_f is the immediate electron donor to P-680⁺ in tris-washed Photosystem II fragments.