Computational methods for the prediction of the structure and interactions of coiled-coil peptides

The past several years have seen significant advances in the development of computational methods for the prediction of the structure and interactions of coiled-coil peptides. These methods are generally based on pairwise correlations of amino acids, helical propensity, thermal melts and the energetics of sidechain interactions, as well as statistical patterns based on Hidden Markov Model (HMM) and Support Vector Machine (SVM) techniques. These methods are complemented by a number of public databases that contain sequences, motifs, domains and other details of coiled-coil structures identified by various algorithms. Some of these computational methods have been developed to make predictions of coiled-coil structure on the basis of sequence information; however, structural predictions of the oligomerisation state of these peptides still remains largely an open question due to the dynamic behaviour of these molecules. This review focuses on existing in silico methods for the prediction of coiled-coil peptides of functional importance using sequence and/or three-dimensional structural data.

Polypurine/polypyrimidine sequences as cis-acting transcriptional regulators

Genome sequence information has generated increasing evidence for the claim that repetitive DNA sequences present within and around genes could play a important role in the regulation of gene expression. Polypurine/polypyrimidine sequences [poly(Pu/Py)] have been observed in the vicinity of promoters and within the transcribed regions of many genes. To understand whether such sequences influence the level of gene expression, we constructed several prokaryotic and eukaryotic expression vectors incorporating poly(Pu/Py) repeats both within and upstream of a reporter gene, lacZ (encoding β-galactosidase), and studied its expression in vivo. We find that, in contrast to the situation in Escherichia coli, the presence of poly(Pu/Py) sequences within the gene does not significantly inhibit gene expression in mammalian cells. On the other hand, the presence of such sequences upstream of lacZ leads to a several-fold reduction of gene expression in mammalian cells. Similar down-regulation was observed when a structural cassette containing poly(Pu/Py) sequences upstream of lacZ was integrated into yeast chromosome V. Sequence analysis of the nine totally sequenced yeast chromosomes shows that a large number of such sequences occur upstream of ORFs. On the basis of our experimental results and DNA sequence analysis, we propose that these sequences can function as cis-acting transcriptional regulators.

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5–C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.

Prediction of heparin binding sites in bone morphogenetic proteins (BMPs)

Heparin is a glycosaminoglycan known to bind bone morphogenetic proteins (BMPs) and the growth and differentiation factors (GDFs) and has strong and variable effects on BMP osteogenic activity. In this paper we report our predictions of the likely heparin binding sites for BMP-2 and 14. The N-terminal sequences upstream of TGF-β-type cysteine-knot domains in BMP-2, 7 and 14 contain the basic residues arginine and lysine, which are key components of the heparin/HS-binding sites, with these residues being highly non-conserved. Importantly, evolutionary conserved surfaces on the beta sheets are required for interactions with receptors and antagonists. Furthermore, BMP-2 has electropositive surfaces on two sides compared to BMP-7 and BMP-14. Molecular docking simulations suggest the presence of high and low affinity binding sites in dimeric BMP-2. Histidines were found to play a role in the interactions of BMP-2 with heparin; however, a pKa analysis suggests that histidines are likely not protonated. This is indicative that interactions of BMP-2 with heparin do not require acidic pH. Taken together, non-conserved amino acid residues in the N-terminus and residues protruding from the beta sheet (not overlapping with the receptor binding sites and the dimeric interface) and not C-terminal are found to be important for heparin–BMP interactions.

A comparative structural bioinformatics analysis of the insulin receptor family ectodomain based on phylogenetic information

The insulin receptor (IR), the insulin-like growth factor 1 receptor (IGF1R) and the insulin receptor-related receptor (IRR) are covalently-linked homodimers made up of several structural domains. The molecular mechanism of ligand binding to the ectodomain of these receptors and the resulting activation of their tyrosine kinase domain is still not well understood. We have carried out an amino acid residue conservation analysis in order to reconstruct the phylogeny of the IR Family. We have confirmed the location of ligand binding site 1 of the IGF1R and IR. Importantly, we have also predicted the likely location of the insulin binding site 2 on the surface of the fibronectin type III domains of the IR. An evolutionary conserved surface on the second leucine-rich domain that may interact with the ligand could not be detected. We suggest a possible mechanical trigger of the activation of the IR that involves a slight ‘twist’ rotation of the last two fibronectin type III domains in order to face the likely location of insulin. Finally, a strong selective pressure was found amongst the IRR orthologous sequences, suggesting that this orphan receptor has a yet unknown physiological role which may be conserved from amphibians to mammals.

Plant regeneration from protoplasts of Capsicum annuum L cv. California Wonder

Plant regeneration from mesophyll protoplasts of pepper, Capsicum annuum L. cv. California Wonder has been demonstrated via shoot organogenesis, Protoplasts isolated from fully expanded leaves of 3-week-old axenic shoots when cultured in TM medium supplemented with 1 mgl(-1) NAA, 1 mgl(-1) 2, 4-D, 0.5 mgl(-1) BAP (CM 1) resulted in divisions with a frequency ranging from 20-25%. Antioxidant ascorbic acid and polyvinylpyrrolidone (PVP) in the medium and incubation in the dark helped overcome browning of protoplasts. Microcalli and macrocalli were formed in TM medium containing 2 mgl(-1) NAA and 0.5 mgl(-1) BAP (CM II) and MS gelled medium containing 2 mgl(-1) NAA and 0.5 mgl(-1) BAP (CM III), respectively, Regeneration of plantlets was possible via caulogenesis, Microshoots, 2-5 per callus appeared on MS gelled medium enriched with 0.5 mgl(-1) IAA, 2 mgl(-1) GA and 10 mgl(-1) BAP (CM IVc). Rooting of microshoots was obtained on half strength gelled medium containing 1 mgl(-1) NAA and 0.5 mgl(-1) BAP, Protoplasts isolated from cotyledons failed to divide and degenerated eventually.

X-ray studies on crystalline complexes involving amino acids and peptides. XXIII. Variability in ionization state, conformation and molecular aggregation in the complexes of succinic acid with DL- and L-lysine

Crystalline complexes of succinic acid with DL- and L-lysine have been prepared and analysed by X-ray diffraction. DL-Lysine complex: C6HIsN202 + 1 2- 1 ~C4H404 .~C4H604, Mr -- 264"2, PI, a = 5"506 (4), =8.070(2), c=14.089(2) A,, a=92.02(1), /3= 100"69 (3), y = 95"85 (3) ~>, Z = 2, Dx = 1"44 g cm -3, R = 0.059 for 2546 observed reflections. Form I of the e-lysine complex: C6HIsN20-, ~ .C4H504, Mr = 264.2, P1, a = 5" 125 (2), b = 8"087 (1), c = 8"689 (1) A,, a = 112.06 (1), /3 = 99.08 (2), y = 93"77(2) °, Z--l, D,,,=1"34(3), Dx=l"34gcm 3 R = 0.033 for 1475 observed reflections. Form II of + I 2- the e-lysine complex: C6H15N202 .,iC4H404 .- 1 I ") 4C4H604.4(C4HsO4""H'"CaH404)" , Mr = 264"2, P1, a = 10.143 (4), b = 10.256 (2), c = 12"916 (3) A,, a = 105.00 (2),/3 = 99-09 (3), y = 92"78 (3)::, Z = 4, Dm= 1"37(4), D,.= 1.38gcm 3, R=0.067 for 2809 observed reflections. The succinic acid molecules in the structures exhibit a variety of ionization states. Two of the lysine conformations found in the complexes have been observed for the first time in crystals containing lysine. Form II of the L-lysine complex is highly pseudosymmetric. In all the complexes, unlike molecules aggregate into separate alternating layers. The basic element of aggregation in the lysine layer in the complexes is an S2-type head-to-tail sequence. This element combines in different ways in the three structures. The basic element of aggre gation in the succinic acid layer in the complexes is a hydrogen-bonded ribbon. The ribbons are interconnected indirectly through amino groups in the lysine layer.

Characterization of Somatic Embryogenesis in Sandalwood (Santalum album L.)

In the synchronous embryogenesis system of sandalwood developed in our laboratory, we observed that the early events of differentiation from freshly induced callus (stage 0) are accomplished in three distinct stages viz., preglobular masses (stage 1), globular embryos (stage 2), and bipolar embryos (stage 3). Transition from stage 0 to 1 was accomplished using 2,4-D and involves a stage specific appearance of two polypeptides of 15 and 30 kDa molecular weight. A 24 kDa polypeptide that was detected as a marked band in extracts of primary callus was not detected in stages 1, 2, and 3. Further, the tissue level of a 50 kDa glycoprotein decreased during transition from stage 2 to stage 3. However, the levels of glycoproteins in the medium were markedly higher in stage 0 cultures compared to those in stage 1. The activities of a protein kinase, glycosidase, and xylanase increased markedly with progressing embryogenesis. Our observations suggest that in addition to being controlled at the level of stage-specific gene expression, somatic embryogenesis in sandalwood is also regulated at the level of controls on cell wall flexibility and posttranslational changes in the pool of preexisting proteins.

Analysis of the amino acid sequences of plant Bowman-Birk inhibitors

Plant seeds contain a large number of protease inhibitors of animal, fungal, and bacterial origin. One of the well-studied families of these inhibitors is the Bowman-Birk family(BBI). The BBIs from dicotyledonous seeds are 8K, double-headed proteins. In contrast, the 8K inhibitors from monocotyledonous seeds are single headed. Monocots also have a 16K, double-headed inhibitor. We have determined the primary structure of a Bowman-Birk inhibitor from a dicot, horsegram, by sequential edman analysis of the intact protein and peptides derived from enzymatic and chemical cleavage. The 76-residue-long inhibitor is very similar to that ofMacrotyloma axillare. An analysis of this inhibitor along with 26 other Bowman-Birk inhibitor domains (MW 8K) available in the SWISSPROT databank revealed that the proteins from monocots and dicots belong to related but distinct families. Inhibitors from monocots show larger variation in sequence. Sequence comparison shows that a crucial disulphide which connects the amino and carboxy termini of the active site loop is lost in monocots. The loss of a reactive site in monocots seems to be correlated to this. However, it appears that this disulphide is not absolutely essential for retention of inhibitory function. Our analysis suggests that gene duplication leading to a 16K inhibitor in monocots has occurred, probably after the divergence of monocots and dicots, and also after the loss of second reactive site in monocots.

The ecology of the mountain hare (Lepus timidus L.) in managed boreal forests: habitat associations at different forest scales

Aims of this thesis This study is part of a larger hare project in Finland, which provides answers to basic ecological questions regarding the mountain hare. This study of the ecology of the mountain hare focuses in particular on different levels of managed boreal forest. The feeding habits and intensity of mountain hares in winter are explored, and the connections between mountain hares versus the forest structure are also studied (e.g. habitat use and the importance of different forest layers for hares). The use of the environment by hares at the landscape level was examined (forest patch structures), and the home ranges of mountain hares were studied. Finally, the productivity and survival rate of mountain hare populations were also studied (discussion e.g. predator effects on hare populations). Conclusions Feeding intensity seemed to be highest in the spring-winter, when home ranges were also largest. Favourable food species are covered by snow in winter and the mobility of hares is highest during late winter. A shortage of suitable food species may be problematic for hares, especially during the winter period. In this study mountain hares preferred a dense shrub layer at local level and deciduous and mixed tree forest over coniferous forest at the landscape level. Food and shelter are vital for hares and the preference for particular habitats may also affect the population dynamics of the mountain hare. It would be possible to improve the quality of food and shelter or at least prevent the most negative habitat changes through forest management. At a local level it is also possible to add supplementary food for hares through the winter period. The intensive clearing of young sapling stands and especially the removal of deciduous shrubs and trees reduces the quality of habitats for the mountain hare. Mountain hares primarily live in forest habitat and it is possible that changes in the forest structure play a crucial role in mountain hare habitat preference. Ecological knowledge of the mountain hare is vital to create habitat structure more suitable for the species. More deciduous trees should be saved in managing forests and the mechanical clearing of the shrub layer should be done carefully.

Lignin biosynthesis in Norway spruce: from a model system to the tree

Lignin is a hydrophobic polymer that is synthesised in the secondary cell walls of all vascular plants. It enables water conduction through the stem, supports the upright growth habit and protects against invading pathogens. In addition, lignin hinders the utilisation of the cellulosic cell walls of plants in pulp and paper industry and as forage. Lignin precursors are synthesised in the cytoplasm through the phenylpropanoid pathway, transported into the cell wall and oxidised by peroxidases or laccases to phenoxy radicals that couple to form the lignin polymer. This study was conducted to characterise the lignin biosynthetic pathway in Norway spruce (Picea abies (L.) Karst.). We focused on the less well-known polymerisation stage, to identify the enzymes and the regulatory mechanisms that are involved. Available data for lignin biosynthesis in gymnosperms is scarce and, for example, the latest improvements in precursor biosynthesis have only been verified in herbaceous plants. Therefore, we also wanted to study in detail the roles of individual gene family members during developmental and stress-induced lignification, using EST sequencing and real-time RT-PCR. We used, as a model, a Norway spruce tissue culture line that produces extracellular lignin into the culture medium, and showed that lignin polymerisation in the tissue culture depends on peroxidase activity. We identified in the culture medium a significant NADH oxidase activity that could generate H2O2 for peroxidases. Two basic culture medium peroxidases were shown to have high affinity to coniferyl alcohol. Conservation of the putative substrate-binding amino acids was observed when the spruce peroxidase sequences were compared with other peroxidases with high affinity to coniferyl alcohol. We also used different peroxidase fractions to produce synthetic in vitro lignins from coniferyl alcohol; however, the linkage pattern of the suspension culture lignin could not be reproduced in vitro with the purified peroxidases, nor with the full complement of culture medium proteins. This emphasised the importance of the precursor radical concentration in the reaction zone, which is controlled by the cells through the secretion of both the lignin precursors and the oxidative enzymes to the apoplast. In addition, we identified basic peroxidases that were reversibly bound to the lignin precipitate. They could be involved, for example, in the oxidation of polymeric lignin, which is required for polymer growth. The dibenzodioxocin substructure was used as a marker for polymer oxidation in the in vitro polymerisation studies, as it is a typical substructure in wood lignin and in the suspension culture lignin. Using immunolocalisation, we found the structure mainly in the S2+S3 layers of the secondary cell walls of Norway spruce tracheids. The structure was primarily formed during the late phases of lignification. Contrary to the earlier assumptions, it appears to be a terminal structure in the lignin macromolecule. Most lignin biosynthetic enzymes are encoded for by several genes, all of which may not participate in lignin biosynthesis. In order to identify the gene family members that are responsible for developmental lignification, ESTs were sequenced from the lignin-forming tissue culture and developing xylem of spruce. Expression of the identified lignin biosynthetic genes was studied using real-time RT-PCR. Candidate genes for developmental lignification were identified by a coordinated, high expression of certain genes within the gene families in all lignin-forming tissues. However, such coordinated expression was not found for peroxidase genes. We also studied stress-induced lignification either during compression wood formation by bending the stems or after Heterobasidion annosum infection. Based on gene expression profiles, stress-induced monolignol biosynthesis appeared similar to the developmental process, and only single PAL and C3H genes were specifically up-regulated by stress. On the contrary, the up-regulated peroxidase genes differed between developmental and stress-induced lignification, indicating specific responses.

Purification and properties of protocatechuate-3,4-dioxygenase from Tecoma stans L

Protocatechuate-3,4-dioxygenase from the leaves of Tecoma stans was purified to near homogeneity and some of its properties studied. It was optimally active at pH 5.2 and at 40°C. Its molecular weight of approx. 150 000 was determined by gel filtration on a Sephadex G-150 column. The Km value for protocatechuate was found to be 330 μM and for ferrous sulfate, 40 μM. The enzyme was highly specific for protocatechuate and did not attack any of the substrate analogues. None of the substrate analogues tested inhibited the enzyme activity. Sulfhydryl reagents inhibited the enzyme activity which could be partially reversed by sulfhydryl compounds. The dioxygenase activity was not associated with polyphenol oxidase activity.