Retinal-Protein Interactions in Halorhodopsin from Natronomonas pharaonis: Binding and Retinal Thermal Isomerization Catalysis

Halorhodopsin from Natronomonas pharaonis (NpHR) is a member of the retinal protein group and serves as a light-driven chloride pump in which chloride ions are transported through the membrane following light absorption by the retinal chromophore. In this study, we examined two main issues: (1) factors controlling the binding of the retinal chromophore to the NpHR opsin and (2) the ability of the NpHR opsin to catalyze the thermal isomerization of retinal isomers. We have revealed that the reconstitution process of pharaonis HR (NpHR) pigment from its apoprotein and all-trans retinal depends on the pH, and the process has a pK_a of 5.8 ± 0.1. It was proposed that this pK_a is associated with the pK_a of the lysine residue that binds the retinal chromophore (Lys256). The pigment formation is regulated by the concentration of sodium chloride, and the maximum yield was observed at 3.7 M NaCl. The low yield of pigment in a lower concentration of NaCl (< 3 M) may be due to an altered conformation adopted by the apomembrane, which is not capable of forming the pigment. Unexpectedly and unlike the apomembrane of bacteriorhodopsin, NpHR opsin produces pigments with 11-cis retinal and 9-cis retinal owing to the thermal isomerization of these retinal isomers to all-trans retinal. The isomerization rate depends on the pH, and it is faster at a higher pH. The pK_a value of the isomerization process is similar to the pK_a of the binding process of these retinals, which suggests that Lys256 is also involved in the isomerization process. The isomerization is independent of the sodium chloride concentration. However, in the absence of sodium chloride, the apoprotein adopts such a conformation, which does not prevent the isomerization of retinal, but it prevents a covalent bond formation with the lysine residue. The rate and the thermodynamic parameter analysis of the retinal isomerization by NpHR apoprotein led to the conclusion that the apomembrane catalyzes the isomerization via a triplet mechanism.     

Energy cost for the proper ionization of active site residues in 6-phosphogluconate dehydrogenase from T. brucei

The catalytic mechanism of 6-phosphogluconate dehydrogenase requires the inversion of a Lys/Glu couple from its natural ionization state. The pK_a of these residues in free and substrate bound enzymes has been determined measuring by ITC the proton release/uptake induced by substrate binding at different pH values. Wt 6-phosphogluconate dehydrogenase from Trypanosoma brucei and two active site enzyme mutants, K185H and E192Q were investigated. Substrate binding was accompanied by proton release and was dependent on the ionization of a group with pK_a 7.07 which was absent in the E192Q mutant. Kinetic data highlighted two pK_a, 7.17 and 9.64, in the enzyme–substrate complex, the latter being absent in the E192Q mutant, suggesting that the substrate binding shifts Glu192 pK_a from 7.07 to 9.64. A comparison of wt and E192Q mutant appears to show that the substrate binding shifts Lys185 pK_a from 9.9 to 7.17. By comparing differences in proton release and the binding enthalpy of wt and mutant enzymes, the enthalpic cost of the change in the protonation state of Lys185 and Glu192 was estimated at ≈ 6.1 kcal/mol. The change in protonation state of Lys185 and Glu192 has little effect on Gibbs free energy, 240–325 cal/mol. However proton balance evidences the dissociation of other group(s) that can be collectively described by a single pK_a shift from 9.1 to 7.54. This further change in ionization state of the enzyme causes an increase of free energy with a total cost of 1.2–2.3 kcal/mol to set the enzyme into a catalytically competent form.     

The pKa Values of Acidic and Basic Residues Buried at the Same Internal Location in a Protein Are Governed by Different Factors

The pK_a values of internal ionizable groups are usually very different from the normal pK_a values of ionizable groups in water. To examine the molecular determinants of pK_a values of internal groups, we compared the properties of Lys, Asp, and Glu at internal position 38 in staphylococcal nuclease. Lys38 titrates with a normal or elevated pK_a, whereas Asp38 and Glu38 titrate with elevated pK_a values of 7.0 and 7.2, respectively. In the structure of the L38K variant, the buried amino group of the Lys38 side chain makes an ion pair with Glu122, whereas in the structure of the L38E variant, the buried carboxyl group of Glu38 interacts with two backbone amides and has several nearby carboxyl oxygen atoms. Previously, we showed that the pK_a of Lys38 is normal owing to structural reorganization and water penetration concomitant with ionization of the Lys side chain. In contrast, the pK_a values of Asp38 and Glu38 are perturbed significantly owing to an imbalance between favorable polar interactions and unfavorable contributions from dehydration and from Coulomb interactions with surface carboxylic groups. Their ionization is also coupled to subtle structural reorganization. These results illustrate the complex interplay between local polarity, Coulomb interactions, and structural reorganization as determinants of pK_a values of internal groups in proteins. This study suggests that improvements to computational methods for pK_a calculations will require explicit treatment of the conformational reorganization that can occur when internal groups ionize.     

Structural origins of high apparent dielectric constants experienced by ionizable groups in the hydrophobic core of a protein

The side chains of Lys66, Asp66, and Glu66 in staphylococcal nuclease are fully buried and surrounded mainly by hydrophobic matter, except for internal water molecules associated with carboxylic oxygen atoms. These ionizable side chains titrate with pK_a values of 5.7, 8.8, and 8.9, respectively. To reproduce these pK_a values with continuum electrostatics calculations, we treated the protein with high dielectric constants. We have examined the structural origins of these high apparent dielectric constants by using NMR spectroscopy to characterize the structural response to the ionization of these internal side chains. Substitution of Val66 with Lys66 and Asp66 led to increased conformational fluctuations of the microenvironments surrounding these groups, even under pH conditions where Lys66 and Asp66 are neutral. When Lys66, Asp66, and Glu66 are charged, the proteins remain almost fully folded, but resonances for a few backbone amides adjacent to the internal ionizable residues are broadened. This suggests that the ionization of the internal groups promotes a local increase in dynamics on the intermediate timescale, consistent with either partial unfolding or increased backbone fluctuations of helix 1 near residue 66, or, less likely, with increased fluctuations of the charged side chains at position 66. These experiments confirm that the high apparent dielectric constants reported by internal Lys66, Asp66, and Glu66 reflect localized changes in conformational fluctuations without incurring detectable global structural reorganization. To improve structure-based pK_a calculations in proteins, we will need to learn how to treat this coupling between ionization of internal groups and local changes in conformational fluctuations explicitly.     

The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK a

Several small ribozymes employ general acid–base catalysis as a mechanism to enhance site-specific RNA cleavage, even though the functional groups on the ribonucleoside building blocks of RNA have pK_a values far removed from physiological pH. The rate of the cleavage reaction is strongly affected by the identity of the metal cation present in the reaction solution; however, the mechanism(s) by which different cations contribute to rate enhancement has not been determined. Using the Neurospora VS ribozyme, we provide evidence that different cations confer particular shifts in the apparent pK_a values of the catalytic nucleobases, which in turn determines the fraction of RNA in the protonation state competent for general acid–base catalysis at a given pH, which determines the observed rate of the cleavage reaction. Despite large differences in observed rates of cleavage in different cations, mathematical models of general acid–base catalysis indicate that k₁, the intrinsic rate of the bond-breaking step, is essentially constant irrespective of the identity of the cation(s) in the reaction solution. Thus, in contrast to models that invoke unique roles for metal ions in ribozyme chemical mechanisms, we find that most, and possibly all, of the ion-specific rate enhancement in the VS ribozyme can be explained solely by the effect of the ions on nucleobase pK_a. The inference that k₁ is essentially constant suggests a resolution of the problem of kinetic ambiguity in favor of a model in which the lower pK_a is that of the general acid and the higher pK_a is that of the general base.     

Cleavage of peptide bonds bearing ionizable amino acids at P1 by serine proteases with hydrophobic S1 pocket

Enzymatic hydrolysis of the synthetic substrate succinyl-Ala-Ala-Pro-Xxx-pNA (where Xxx = Leu, Asp or Lys) catalyzed by bovine chymotrypsin (CHYM) or Streptomyces griseus protease B (SGPB) has been studied at different pH values in the pH range 3-11. The pH optima for substrates having Leu, Asp, and Lys have been found to be 7.5-8.0, 5.5-6.0, and ∼10, respectively. At the normally reported pH optimum (pH 7-8) of CHYM and SGPB, the substrate with Leu at the reactive site is more than 25,000-fold more reactive than that with Asp. However, when fully protonated, Asp is nearly as good a substrate as Leu. The pK values of the side chains of Asp and Lys in the hydrophobic S₁ pocket of CHYM and SGPB have been calculated from pH-dependent hydrolysis data and have been found to be about 9 for Asp and 7.4 and 9.7 for Lys for CHYM and SGPB, respectively. The results presented in this communication suggest a possible application of CHYM like enzymes in cleaving peptide bonds contributed by acidic amino acids between pH 5 and 6.     

Mutation of a pH-modulating residue in a GH51 α-l-arabinofuranosidase leads to a severe reduction of the secondary hydrolysis of transfuranosylation products

Background

The development of enzyme-mediated glycosynthesis using glycoside hydrolases is still an inexact science, because the underlying molecular determinants of transglycosylation are not well understood. In the framework of this challenge, this study focused on the family GH51 α-l-arabinofuranosidase from Thermobacillus xylanilyticus, with the aim to understand why the mutation of position 344 provokes a significant modification of the transglycosylation/hydrolysis partition.

Methods

Detailed kinetic analysis (k_cat, K_M, pK_a determination and time-course NMR kinetics) and saturation transfer difference nuclear magnetic resonance spectroscopy was employed to determine the synthetic and hydrolytic ability modification induced by the redundant N344 mutation disclosed in libraries from directed evolution.

Results

The mutants N344P and N344Y displayed crippled hydrolytic abilities, and thus procured improved transglycosylation yields. This behavior was correlated with an increased pK_a of the catalytic nucleophile (E298), the pK_a of the acid/base catalyst remaining unaffected. Finally, mutations at position 344 provoked a pH-dependent product inhibition phenomenon, which is likely to be the result of a significant modification of the proton sharing network in the mutants.

Conclusions and general significance

Using a combination of biochemical and biophysical methods, we have studied TxAbf-N344 mutants, thus revealing some fundamental details concerning pH modulation. Although these results concern a GH51 α-l-arabinofuranosidase, it is likely that the general lessons that can be drawn from them will be applicable to other glycoside hydrolases. Moreover, the effects of mutations at position 344 on the transglycosylation/hydrolysis partition provide clues as to how TxAbf can be further engineered to obtain an efficient transfuranosidase.     

Combined DFT and electrostatics study of the proton pumping mechanism in cytochrome c oxidase

Cytochrome c oxidase is a redox-driven proton pump which converts atmospheric oxygen to water and couples the oxygen reduction reaction to the creation of a membrane proton gradient. The structure of the enzyme has been solved; however, the mechanism of proton pumping is still poorly understood. Recent calculations from this group indicate that one of the histidine ligands of enzyme's Cu_B center, His291, may play the role of the pumping element. In this paper, we report on the results of calculations that combined first principles DFT and continuum electrostatics to evaluate the energetics of the key energy generating step of the model—the transfer of the chemical proton to the binuclear center of the enzyme, where the hydroxyl group is converted to water, and the concerted expulsion of the proton from δ-nitrogen of His291 ligand of Cu_B center. We show that the energy generated in this step is sufficient to push a proton against an electrochemical membrane gradient of about 200 mV. We have also re-calculated the pK_a of His291 for an extended model in which the whole Fe_a3-Cu_B center with their ligands is treated by DFT. Two different DFT functionals (B3LYP and PBE0), and various dielectric models of the protein have been used in an attempt to estimate potential errors of the calculations. Although current methods of calculations do not allow unambiguous predictions of energetics in proteins within few pK_a units, as required in this case, the present calculation provides further support for the proposed His291 model of CcO pump and makes a specific prediction that could be targeted in the experimental test.     

An evolutionally conserved Lys122 is essential for function in Rhodospirillum rubrum bona fide RuBisCO and Bacillus subtilis RuBisCO-like protein

Ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO) and RuBisCO-like protein (RLP) catalyze similar enolase-type reactions. Both enzymes have a conserved non-catalytic Lys122 or Arg122 on the β-strand E lying in the interface between the N- and C-terminal domains. We used site-directed mutagenesis to analyze the function of Lys122 in the form II Rhodospirillum rubrum RuBisCO (RrRuBisCO) and Bacillus subtilis RLP (BsRLP). The K122R mutant of RrRuBisCO had a 40% decrease in k_cat for carboxylase activity, a 2-fold increase in K_m for CO₂, and a 1.9-fold increase in K_m for ribulose-1,5-bisphosphate. K122M and K122E mutants of RrRuBisCO were almost inactive. None of the substitutions affected the thermal stability of RrRuBisCO. The K122R mutant of BsRLP had a 32% decrease in k_cat and lower thermal stability than the wild-type enzyme. The K122M and K122E mutants of BsRLP failed to form a catalytic dimer. Our results suggest that the lysine residue is essential for function in both enzymes, although in each case, its role is likely distinct.     

Weakened coupling of conserved arginine to the proteorhodopsin chromophore and its counterion implies structural differences from bacteriorhodopsin

In wild-type proteorhodopsin (pR), titration of the chromophore's counterion Asp⁹⁷ occurs with a pK_a of 8.2±0.1. R94C mutation reduces this slightly to 7.0±0.2, irrespective of treatment with ethylguanidinium. This contrasts with the homologous archaeal protein bacteriorhodopsin (bR), where R82C mutation was previously shown to elevate the pK_a of Asp⁸⁵ by ∼5 units, while reconstitution with ethylguanidinium restores it nearly to the wild-type value of 2.5. We conclude there is much weaker electrostatic coupling between Arg⁹⁴ and Asp⁹⁷ in the unphotolyzed state of pR, in comparison to Arg⁸² and Asp⁸⁵ in bR. Therefore, while fast light-driven H⁺ release may depend on these two residues in pR as in bR, no tightly conserved pre-photolysis configuration of them is required.     

Catalytic role of the calcium ion in GH97 inverting glycoside hydrolase

The role of calcium ion in the active site of the inverting glycoside hydrolase family 97 enzyme, BtGH97a, was investigated through structural and kinetic studies. The calcium ion was likely directly involved in the catalytic reaction. The pH dependence of k_cat/K_m values in the presence or absence of calcium ion indicated that the calcium ion lowered the pK_a of the base catalyst. The significant decreases in k_cat/K_m for hydrolysis of substrates with basic leaving groups in the absence of calcium ion confirmed that the calcium ion facilitated the leaving group departure.     

Catalytic mechanism of inulinase from Arthrobacter sp. S37

Detailed catalytic roles of the conserved Glu323, Asp460, and Glu519 of Arthrobacter sp. S37 inulinase (EnIA), a member of the glycoside hydrolase family 32, were investigated by site-directed mutagenesis and pH-dependence studies of the enzyme efficiency and homology modeling were carried out for EnIA and for D460E mutant. The enzyme efficiency (k_cat/K_m) of the E323A and E519A mutants was significantly lower than that of the wild-type due to a substantial decrease in k_cat, but not due to variations in K_m, consistent with their putative roles as nucleophile and acid/base catalyst, respectively. The D460A mutant was totally inactive, whereas the D460E and D460N mutants were active to some extent, revealing Asp460 as a catalytic residue and demonstrating that the presence of a carboxylate group in this position is a prerequisite for catalysis. The pH-dependence studies indicated that the pK_a of the acid/base catalyst decreased from 9.2 for the wild-type enzyme to 7.0 for the D460E mutant, implicating Asp460 as the residue that interacts with the acid/base catalyst Glu519 and elevates its pK_a. Homology modeling and molecular dynamics simulation of the wild-type enzyme and the D460E mutant shed light on the structural roles of Glu323, Asp460, and Glu519 in the catalytic activity of the enzyme.     

Mutation of non-conserved amino acids surrounding catalytic site to shift pH optimum of Bacillus circulans xylanase

Improvement of enzyme function by engineering pH dependence of enzymatic activity is of importance for industrial application of Bacillus circulans xylanases. Target mutation sites were selected by structural alignment between B. circulans xylanase and other xylanases having different pH optima. We selected non-conserved mutant sites within 8 Å from the catalytic residues, to see whether these residues have some role in modulating pK_as of the catalytic residues. We hypothesized that the non-conserved residues which may not have any role in enzyme catalysis might perturb pK_as of the catalytic residues. Change in pK_a of a titratable group due to change in electrostatic potential of a mutation was calculated and the change in pH optimum was predicted from the change in pK_a of the catalytic residues. Our strategy is proved to be useful in selection of promising mutants to shift the pH optimum of the xylanases towards desired side.     

Electrostatic contributions to the stability of the GCN4 leucine zipper structure

Ion pairs are ubiquitous in X-ray structures of coiled coils, and mutagenesis of charged residues can result in large stability losses. By contrast, pK_a values determined by NMR in solution often predict only small contributions to stability from charge interactions. To help reconcile these results we used triple-resonance NMR to determine pK_a values for all groups that ionize between pH 1 and 13 in the 33 residue leucine zipper fragment, GCN4p. In addition to the native state we also determined comprehensive pK_a values for two models of the GCN4p denatured state: the protein in 6 M urea, and unfolded peptide fragments of the protein in water. Only residues that form ion pairs in multiple X-ray structures of GCN4p gave large pK_a differences between the native and denatured states. Moreover, electrostatic contributions to stability were not equivalent for oppositely charged partners in ion pairs, suggesting that the interactions between a charge and its environment are as important as those within the ion pair. The pH dependence of protein stability calculated from NMR-derived pK_a values agreed with the stability profile measured from equilibrium urea-unfolding experiments as a function of pH. The stability profile was also reproduced with structure-based continuum electrostatic calculations, although contributions to stability were overestimated at the extremes of pH. We consider potential sources of errors in the calculations, and how pK_a predictions could be improved. Our results show that although hydrophobic packing and hydrogen bonding have dominant roles, electrostatic interactions also make significant contributions to the stability of the coiled coil.     

Electrostatic analysis of the hepatitis C virus NS3 helicase reveals both active and allosteric site locations

Multi-conformation continuum electrostatics (MCCE) was used to analyze various structures of the NS3 RNA helicase from the hepatitis C virus in order to determine the ionization state of amino acid side chains and their pK_as. In MCCE analyses of HCV helicase structures that lacked ligands, several active site residues were identified to have perturbed pK_as in both the nucleic acid binding site and in the distant ATP-binding site, which regulates helicase movement. In all HCV helicase structures, Glu493 was unusually basic and His369 was abnormally acidic. Both these residues are part of the HCV helicase nucleic acid binding site, and their roles were analyzed by examining the pH profiles of site-directed mutants. Data support the accuracy of MCCE predicted pK_a values, and reveal that Glu493 is critical for low pH enzyme activation. Several key residues, which were previously shown to be involved in helicase-catalyzed ATP hydrolysis, were also identified to have perturbed pK_as including Lys210 in the Walker-A motif and the DExD/H-box motif residues Asp290 and His293. When DNA was present in the structure, the calculated pK_as shifted for both Lys210 and Asp290, demonstrating how DNA binding might lead to electrostatic changes that stimulate ATP hydrolysis.     

Properties of polyphenol oxidase from tubers of the yam Dioscorea bulbifera

The subunit MW of Dioscorea bulbifera polyphenol oxidase (MW 115 000 ± 2000) determined by SDS-PAGE is ca. 31 000 indicating that the enzyme is an oligomeric protein with four subunits. K_i values of various inhibitors and their modes of inhibition have been determined with catechol and pyrogallol as substrates. p-Nitrophenol, p-cresol, quinoline and resorcinol are competitive inhibitors of catechol binding while only orcinol and p-nitrophenol behave in the same way towards pyrogallol as substrate. From the effect of pH on V_max, groups with pK values ca. 4.7 and 6.8 have been identified to be involved in catalytic activity. The Arrhenius activation energy (E_a) at pH 4.0 is 8.9 kcal/mol between 40–65°. At pH 7.0, the value is 22.1 kcal/mol between 40 and 60°. The enthalpies (ΔH) at pH 4.0 and pH 7.0 are 2.3 kcal/mol and 32.4 kcal/mol respectively. The results are discussed considering the conformational changes of the enzyme during substrate binding.     

Oxidizing substrate specificity of Mycobacterium tuberculosis alkyl hydroperoxide reductase E: kinetics and mechanisms of oxidation and overoxidation

Alkyl hydroperoxide reductase E (AhpE), a novel subgroup of the peroxiredoxin family, comprises Mycobacterium tuberculosis AhpE (MtAhpE) and AhpE-like proteins present in many bacteria and archaea, for which functional characterization is scarce. We previously reported that MtAhpE reacted ~ 10³ times faster with peroxynitrite than with hydrogen peroxide, but the molecular reasons for that remained unknown. Herein, we investigated the oxidizing substrate specificity and the oxidative inactivation of the enzyme. In most cases, both peroxidatic thiol oxidation and sulfenic acid overoxidation followed a trend in which those peroxides with the lower leaving-group pK_a reacted faster than others. These data are in agreement with the accepted mechanisms of thiol oxidation and support that overoxidation occurs through sulfenate anion reaction with the protonated peroxide. However, MtAhpE oxidation and overoxidation by fatty acid-derived hydroperoxides (~ 10⁸ and 10⁵ M^− 1 s^− 1, respectively, at pH 7.4 and 25 °C) were much faster than expected according to the Brønsted relationship with leaving-group pK_a. A stoichiometric reduction of the arachidonic acid hydroperoxide 15-HpETE to its corresponding alcohol was confirmed. Interactions of fatty acid hydroperoxides with a hydrophobic groove present on the reduced MtAhpE surface could be the basis of their surprisingly fast reactivity.     

High-light induced singlet oxygen formation in cytochrome b6f complex from Bryopsis corticulans as detected by EPR spectroscopy

Electron paramagnetic resonance (EPR) spectroscopy was used to detect the light-induced formation of singlet oxygen (¹O₂*) in the intact and the Rieske-depleted cytochrome b₆f complexes (Cyt b₆f) from Bryopsis corticulans, as well as in the isolated Rieske Fe–S protein. It is shown that, under white-light illumination and aerobic conditions, chlorophyll a (Chl a) bound in the intact Cyt b₆f can be bleached by light-induced ¹O₂*, and that the ¹O₂* production can be promoted by D₂O or scavenged by extraneous antioxidants such as l-histidine, ascorbate, β-carotene and glutathione. Under similar experimental conditions, ¹O₂* was also detected in the Rieske-depleted Cyt b₆f complex, but not in the isolated Rieske Fe–S protein. The results prove that Chl a cofactor, rather than Rieske Fe–S protein, is the specific site of ¹O₂* formation, a conclusion which draws further support from the generation of ¹O₂* with selective excitation of Chl a using monocolor red light.     

Solubilization and purification of cytochrome a1 from Nitrosomonas

Cytochrome a₁ was solubilized with Triton X-100 from a membrane-envelope preparation of Nitrosomonas and partially purified by repeated fractionation with (NH₄)₂SO₄. The purified fraction of cytochrome a₁ was enriched over the crude extract by a factor of 16 and 300 with respect to protein and c-type cytochrome, respectively. The cytochrome was characterized as cytochrome a₁ on the basis of (a) reduced absorption maxima at 444 nm and 595 nm, (b) acid acetone extractibility and ether solubility of the heme and (c) absorption maximum of 587 nm of the ferro-hemochrome in alkaline pyridine. The α absorption band shifted from 600 nm to 595 nm upon solubilization of the cytochrome with Triton X-100. Spectral shifts were observed in the presence of cyanide and azide and the cytochrome changed with aging to a form with a reduced absorption band at 422 nm. Cytochrome a₁ was reduced anaerobically in the presence of reduced mammalian cytochrome c and was rapidly reoxidized in the presence of O₂. CO caused a shift in the soret peak of the reduced form but did not prevent reoxidation of cytochrome a₁ in the presence of CO-O₂ (95:5, v/v).     

The influence of cone age on the relative longevity of Banksia seeds

Background and Aims

Serotiny is common in the genus Banksia, so any seed collection is likely to be comprised of seeds that were produced in many different years. This study aimed to determine the impact of cone age and degree of serotiny on longevity in ex situ storage.

Methods

Cones of identifiable age classes were collected from three species of Banksia. Seeds were extracted from cones and the degree of serotiny calculated. An estimate of initial viability (K_i), the time for viability to fall by one probit (σ) and the relative longevity of seeds (p₅₀) for each species and cone age class was determined using a comparative longevity test (50 °C, 63 % relative humidity).

Key Results

The degree of serotiny ranged from moderate (7·9) for Banksia attenuata to strong (40·4) for B. hookeriana. Survival curves for all seed age classes within each species could be described by regressions with a common slope (1/σ), but with different values for K_i. The time taken for viability to fall by one probit (σ) could be described by a common value (29·1 d) for all three species.

Conclusions

Differences in seed longevity between cone age classes and species was related to variation in initial viability (K_i) rather than to differences in σ. While targeting the youngest mature seed cohort on a plant will maximize the viability of seeds collected, a wide range of age classes should be collected (but stored as separate cohorts if possible) for quality conservation/restoration seed collections where genetic diversity is important.