Identification and characterization of key substructures involved in the early folding events of a (beta/alpha)8-barrel protein as studied by experimental and computational methods

A number of studies have examined the structural properties of late folding intermediates of (beta/alpha)8-barrel proteins involved in tryptophan biosynthesis, whereas there is little information available about the early folding events of these proteins. To identify the contiguous polypeptide segments important to the folding of the (beta/alpha)8-barrel protein Escherichia coli N-(5'-phosphoribosyl)anthranilate isomerase, we structurally characterized fragments and circularly permuted forms of the protein. We also simulated thermal unfolding of the protein using molecular dynamics. Our fragmentation experiments demonstrate that the isolated (beta/alpha)(1-4)beta5 fragment is almost as stable as the full-length protein. The far and near-UV CD spectra of this fragment are indicative of native-like secondary and tertiary structures. Structural analysis of the circularly permutated proteins shows that if the protein is cleaved within the two N-terminal betaalpha modules, the amount of secondary structure is unaffected, whereas, when cleaved within the central (beta/alpha)(3-4)beta5 segment, the protein simply cannot fold. An ensemble of the denatured structures produced by thermal unfolding simulations contains a persistent local structure comprised of beta3, beta4 and beta5. The presence of this three-stranded beta-barrel suggests that it may be an important early-stage folding intermediate. Interactions found in (beta/alpha)(3-4)beta5 may be essential for the early events of ePRAI folding if they provide a nucleation site that directs folding.     

Accelerating all-atom protein folding simulations through reduced dihedral barriers

Protein folding requires extensive changes of backbone and sidechain dihedral angles, whose energy barriers constitute obstacles for the folding kinetics. Folding of small proteins is furthermore thought to be path-independent. Here, we propose that time-consuming all-atom protein folding simulations may be accelerated through a reduction of the dihedral barriers of the force field. In order to investigate this hypothesis, we performed various folding simulations of two small proteins. We report an acceleration towards smaller root-mean-square deviations from the native protein structure using our proposed method.     

Mathematical simulation of the interactions among cyanobacteria, purple sulfur bacteria and chemotrophic sulfur bacteria in microbial mat communities

Abstract: A deterministic one-dimensional reaction diffusion model was constructed to simulate benthic stratification patterns and population dynamics of cyanobacteria, purple and colorless sulfur bacteria as found in marine microbial mats. The model involves the major biogeochemical processes of the sulfur cycle and includes growth metabolism and their kinetic parameters as described from laboratory experimentation. Hence, the metabolic production and consumption processes are coupled to population growth. The model is used to calculate benthic oxygen, sulfide and light profiles and to infer spatial relationships and interactions among the different populations. Furthermore, the model is used to explore the effect of different abiotic and biotic environmental parameters on the community structure. A strikingly clear pattern emerged of the interaction between purple and colorless sulfur bacteria: either colorless sulfur bacteria dominate or a coexistence is found of colorless and purple sulfur bacteria. The model predicts that purple sulfur bacteria only proliferate when the studied environmental parameters surpass well-defined threshold levels. However, once the appropriate conditions do occur, the purple sulfur bacteria are extremely successful as their biomass outweighs that of colorless sulfur bacteria by a factor of up to 17. The typical stratification pattern predicted closely resembles the often described bilayer communities which comprise a layer of purple sulfur bacteria below a cyanobacterial top-layer; colorless sulfur bacteria are predicted to sandwich in between both layers. The profiles of oxygen and sulfide shift on a diel basis similarly as observed in real systems.     

Redox regulation of chloroplast enzymes in Galdieria sulphuraria in view of eukaryotic evolution

Redox modulation is a general mechanism for enzyme regulation, particularly for the post-translational regulation of the Calvin cycle in chloroplasts of green plants. Although red algae and photosynthetic protists that harbor plastids of red algal origin contribute greatly to global carbon fixation, relatively little is known about post-translational regulation of chloroplast enzymes in this important group of photosynthetic eukaryotes. To address this question, we used biochemistry, phylogenetics and analysis of recently completed genome sequences. We studied the functionality of the chloroplast enzymes phosphoribulokinase (PRK, EC 2.7.1.19), NADP-dependent glyceraldehyde 3-phosphate dehydrogenase (NADP-GAPDH, GapA, EC 1.2.1.13), fructose 1,6-bisphosphatase (FBPase, EC 3.1.3.11) and glucose 6-phosphate dehydrogenase (G6PDH, EC 1.1.1.49), as well as NADP-malate dehydrogenase (NADP-MDH, EC 1.1.1.37) in the unicellular red alga Galdieria sulphuraria (Galdieri) Merola. Despite high sequence similarity of G. sulphuraria proteins to those of other photosynthetic organisms, we found a number of distinct differences. Both PRK and GAPDH co-eluted with CP12 in a high molecular weight complex in the presence of oxidized glutathione, although Galdieria CP12 lacks the two cysteines essential for the formation of the N-terminal peptide loop present in higher plants. However, PRK inactivation upon complex formation turned out to be incomplete. G6PDH was redox modulated, but remained in its tetrameric form; FBPase was poorly redox regulated, despite conservation of the two redox-active cysteines. No indication for the presence of plastidic NADP-MDH (and other components of the malate valve) was found.     

Quantifying agonist activity at G protein-coupled receptors

When an agonist activates a population of G protein-coupled receptors (GPCRs), it elicits a signaling pathway that culminates in the response of the cell or tissue. This process can be analyzed at the level of a single receptor, a population of receptors, or a downstream response. Here we describe how to analyze the downstream response to obtain an estimate of the agonist affinity constant for the active state of single receptors. Receptors behave as quantal switches that alternate between active and inactive states (Figure 1). The active state interacts with specific G proteins or other signaling partners. In the absence of ligands, the inactive state predominates. The binding of agonist increases the probability that the receptor will switch into the active state because its affinity constant for the active state (K(b)) is much greater than that for the inactive state (K(a)). The summation of the random outputs of all of the receptors in the population yields a constant level of receptor activation in time. The reciprocal of the concentration of agonist eliciting half-maximal receptor activation is equivalent to the observed affinity constant (K(obs)), and the fraction of agonist-receptor complexes in the active state is defined as efficacy (ε) (Figure 2). Methods for analyzing the downstream responses of GPCRs have been developed that enable the estimation of the K(obs) and relative efficacy of an agonist. In this report, we show how to modify this analysis to estimate the agonist K(b) value relative to that of another agonist. For assays that exhibit constitutive activity, we show how to estimate K(b) in absolute units of M(-1). Our method of analyzing agonist concentration-response curves consists of global nonlinear regression using the operational model. We describe a procedure using the software application, Prism (GraphPad Software, Inc., San Diego, CA). The analysis yields an estimate of the product of K(obs) and a parameter proportional to efficacy (τ). The estimate of τK(obs) of one agonist, divided by that of another, is a relative measure of K(b) (RA(i)). For any receptor exhibiting constitutive activity, it is possible to estimate a parameter proportional to the efficacy of the free receptor complex (τ(sys)). In this case, the K(b) value of an agonist is equivalent to τK(obs)/τ(sys). Our method is useful for determining the selectivity of an agonist for receptor subtypes and for quantifying agonist-receptor signaling through different G proteins.     

Genome organization of Magnaporthe grisea: integration of genetic maps, clustering of transposable elements and identification of genome duplications and rearrangements

 A high-density genetic map of the rice blast fungus Magnaporthe grisea (Guy11×2539) was constructed by adding 87 cosmid-derived RFLP markers to previously generated maps. The new map consists of 203 markers representing 132 independently segregating loci and spans approximately 900 cM with an average resolution of 4.5 cM. Mapping of 33 cosmid probes from the genetic map generated by Sweigard et al. has allowed the integration of two M. grisea maps. The integrated map showed that the linear order of markers along all seven chromosomes in both maps is in good agreement. Thirty of eighty seven markers were derived from cosmid clones that contained the retrotransposon MAGGY (M. grisea gypsy element). Mapping of single-copy DNA sequences associated with the MAGGY cosmids indicated that MAGGY elements are scattered throughout the fungal genome. In eight cases, the probes associated with MAGGY elements showed abnormal segregation patterns. This suggests that MAGGY may be involved in genomic rearrangements. Two RFLP probes linked to MAGGY elements, and another flanking other repetitive DNA elements, identified sequences that were duplicated in the Guy11 genome. Most of the MAGGY cosmids also contained other classes of repetitive DNA suggesting that repetitive DNA sequences tend to cluster in the M. grisea genome. Received: 17 February 1997 / Accepted: 21 February 1997     

Limitations of Metazoan 18S rRNA Sequence Data: Implications for Reconstructing a Phylogeny of the Animal Kingdom and Inferring the Reality of the Cambrian Explosion

We document the phylogenetic behavior of the 18S rRNA molecule in 67 taxa from 28 metazoan phyla and assess the effects of among-site rate variation on reconstructing phylogenies of the animal kingdom. This empirical assessment was undertaken to clarify further the limits of resolution of the 18S rRNA gene as a phylogenetic marker and to address the question of whether 18S rRNA phylogenies can be used as a source of evidence to infer the reality of a Cambrian explosion. A notable degree of among-site rate variation exists between different regions of the 18S rRNA molecule, as well as within all classes of secondary structure. There is a significant negative correlation between inferred number of nucleotide substitutions and phylogenetic information, as well as with the degree of substitutional saturation within the molecule. Base compositional differences both within and between taxa exist and, in certain lineages, may be associated with long branches and phylogenetic position. Importantly, excluding sites with different degrees of nucleotide substitution significantly influences the topology and degree of resolution of maximum-parsimony phylogenies as well as neighbor-joining phylogenies (corrected and uncorrected for among-site rate variation) reconstructed at the metazoan scale. Together, these data indicate that the 18S rRNA molecule is an unsuitable candidate for reconstructing the evolutionary history of all metazoan phyla, and that the polytomies, i.e., unresolved nodes within 18S rRNA phylogenies, cannot be used as a single or reliable source of evidence to support the hypothesis of a Cambrian explosion. Received: 9 December 1997 / Accepted: 23 March 1998     

Purification, characterisation, and inhibition by monoterpenes of acetylcholinesterase from the waxmoth, Galleria mellonella (L.)

Acetylcholinesterase (AChE) was purified from the brain of the waxmoth, Galleria mellonella (L.) by affinity chromatography followed by anion exchange chromatography. It resolved as a single band by polyacrylamide gel electrophoresis (PAGE) both non-denaturing and SDS (silver stained), and as a single peak by high pressure liquid chromatography (HPLC), in an overall yield of 32% representing 283-fold purification. This was a true acetylcholinesterase, with no activity as a non-specific cholinesterase (butyrylcholinesterase). The molecular weight determined by PAGE in the absence and presence of sodium dodecyl sulphate (SDS) was ca. 240,000 Da and 60,000 Da respectively, indicating an arrangement of tetrameric subunits. 2-Heptanone, the honeybee alarm pheromone, reversibly and competitively inhibited the purified Galleria AChE with a K_i value of 1.34×10⁻³ M. Furthermore, five monoterpenes associated with plant defence and representing a range of functional groups, also were reversible competitive inhibitors of the purified AChE from Galleria, which is consistent with previous data for electric eel AChE (Ryan, M.F., Byrne, O., 1988. Plant-insect coevolution and inhibition of acetylcholinesterase. Journal of Chemical Ecology 14, 1965-1975).