Matrix formulation of the transition from a statistical coil to an intramolecular antiparallel beta sheet

A tractible matrix formulation is developed for the formation of intramolecular antiparallel β sheets in a homopolymer chain molecule. The formulation is applicable to chains with a finite degree of polymerization. It can readily be extended to treat specific-sequence heteropolymers. Individual sheets may contain any number of strands, the number of residues per strand can range upward from two, and there is no artificial constraint linking the numbers of residues in adjacent strands. The weighting scheme utilizes two end-effect parameters, denoted by τ and δ. The first parameter is associated with each residue that does not have a partner in a proceding strand, and the latter is associated with each β bend. A third parameter, t, is associated with every residue in the sheet. Conditions are described which lead to the formation of different types of sheets: (1) “sheets” comprised of isolated extended strands; (2) cross-β fibers in which a sheet contains a large number of very short strands; (3) fibers in which a few very long strands run parallel to the fiber axis; (4) sheets comprised of several strands in which the average strand contains five residues. The fourth type of sheet resembles those found in globular proteins. It is formed when τ and δ are both small, with the ratio, τ/δ, being slightly less than one.     

1H NMR structure of an antifungal γ-thionin protein SIα1: Similarity to scorpion toxins

The three-dimensional structure of the Sorghum bicolor seed protein γ-thionin SIα1 has been determined by 2D ¹H nuclear magnetic resonance (NMR) spectroscopy. The secondary structure of this 47-residue antifungal protein with four disulphide bridges consists of a three-stranded antiparallel sheet and one helix. The helix is tethered to the sheet by two disulphide bridges which link two successive turns of the helix to alternate residues i, i + 2 in one strand. Possible binding sites for antifungal activity are discussed. The same fold has been observed previously in several scorpion toxins. Proteins 32:334–349, 1998. © 1998 Wiley-Liss, Inc.     

Helix–coil stability constants for amino acids in nonaqueous solvents. Studies of random poly(γ-benzyl-L-glutamate-co-L-alanine) and poly(γ-benzyl-L-glutamate-co-L-leucine) in a mixed solvent

The host–guest technique has been applied to the determination of the helix–coil stability constants of two naturally occurring amino acids, L -alanine and L -leucine, in a nonaqueous solvent system. Random copolymers containing L -alanine and L -leucine, respectively, as guest residues and γ-benzyl-L -glutamate as the host residue were synthesized. The polymers were fractionated and characterized for their amino acid content, molecular weight, and helix–coil transition behavior in a dichloroacetic acid (DCA)–1,2-dichloroethane (DCE) mixture. Two types of helix–coil transitions were carried out on the copolymers: solvent-induced transitions in DCA–DCE mixtures at 25°C and thermally induced transitions in a 82:18 (wt %) DCA–DCE mixture. The thermally induced transitions were analyzed by statistical mechanical methods to determine the Zimm-Bragg parameters, σ and s, of the guest residues. The experimental data indicate that, in the nonaqueous solvent, the L -alanine residue stabilizes the α-helical conformation more than the L -leucine residue does. This is in contrast to their behavior in aqueous solution, where the reverse is true. The implications of this finding for the analysis of helical structures in globular proteins are discussed.     

Kinetics of coil to α-helix to β-structure transitions of poly(L-ornithine) in low concentrations of sodium dodecyl sulfate

Conformational changes of poly(L-ornithine) [(Orn)_n] were studied in a sodium dodecyl sulfate (NaDodSO₄) solution by CD. (Orn)_n adopted an unstable and a stable helical structure below and above the NaDodSO₄ concentration range where β-structure was favored, respectively. CD stopped-flow was used to monitor the transitions from coil to the unstable helix, from the helix to β-structure, and from coil to β-structure. Only the rate of the helix to β-structure transition was accelerated by an increase in NaDodSO₄ concentration, whereas the rates of the others were independent of NaDodSO₄ concentration. The fractions of coil, α-helix, and β-structure in each conformation of (Orn)_n caused by NaDodSO₄ were computed by simulating a mixed spectrum of typical CD spectra for these structures to the experimentally obtained spectrum. The contents of the unstable and stable helical structures were less than 50 and 73%, respectively.     

Geometrical principles of homomeric β‐barrels and β‐helices: Application to modeling amyloid protofilaments

Examples of homomeric β‐helices and β‐barrels have recently emerged. Here we generalize the theory for the shear number in β‐barrels to encompass β‐helices and homomeric structures. We introduce the concept of the “β‐strip,” the set of parallel or antiparallel neighboring strands, from which the whole helix can be generated giving it n‐fold rotational symmetry. In this context, the shear number is interpreted as the sum around the helix of the fixed register shift between neighboring identical β‐strips. Using this approach, we have derived relationships between helical width, pitch, angle between strand direction and helical axis, mass per length, register shift, and number of strands. The validity and unifying power of the method is demonstrated with known structures including α‐hemolysin, T4 phage spike, cylindrin, and the HET‐s(218‐289) prion. From reported dimensions measured by X‐ray fiber diffraction on amyloid fibrils, the relationships can be used to predict the register shift and the number of strands within amyloid protofilaments. This was used to construct models of transthyretin and Alzheimer β(40) amyloid protofilaments that comprise a single strip of in‐register β‐strands folded into a “β‐strip helix.” Results suggest both stabilization of an individual β‐strip helix and growth by addition of further β‐strip helices can involve the same pair of sequence segments associating with β‐sheet hydrogen bonding at the same register shift. This process would be aided by a repeat sequence. Hence, understanding how the register shift (as the distance between repeat sequences) relates to helical dimensions will be useful for nanotube design.     

Stabilization of short helices by intramolecular cluster formation

The intramolecular formation of multiple clusters of interacting helices has been characterized in a homopolymer. The configuration partition function permits the formation of clusters in which the number of interacting helices may be as large as the greatest integer in n/2, where n denotes the number of amino acid residues in the chain. The theoretical formulation has its origin in a recent [Mattice, W. L. & Scheraga, H. A. (1984) Biopolymers 23 , 1701–1724], tractable matrix expression for the configuration partition function for intramolecular antiparallel β-sheet formation. Reassignment of the expression for one of the n(n+3)/2 elements in the sparse statistical weight matrix, along with a simple change in notation, converts that treatment into a matrix formulation of the configuration partition function for a chain containing multiple clusters of interacting antiparallel helices. The five statistical weights used are δ, f_l, w, and the Zimm-Bragg σ and s. Each tight bend that connects two interacting helices contributes a factor of δ, f_l is used in the weight for larger loops between interacting helices, and w arises from helix–helix interaction. The influence of the helix–helix interaction is well illustrated by two helix–coil transitions in a chain with n = 156 and σ = 0.001. In the absence of helix–helix interaction, the transition occurs by the nucleation and subsequent elongation of a small number of helices. When helix–helix interaction is attractive, the transition can occur by a different mechanism. Formation of a single pair of interacting helices is followed by addition of new helices to the initial cluster. In the latter process, individual helices experience relatively little growth after they are formed.     

Effect of external electrostatic field on the stability of β sheet structures

To explore the effect of an external electrostatic field (EEF) on the stability of protein conformations, the molecular dynamic modeling approach was applied to evaluate the effect of an EEF along the x or y direction on a water cluster containing a parallel or antiparallel β sheet structure. The β sheet structure contained two strands with a (Gly)₃ sequence separated by a distance d along the x direction. The mean forces between the two strands along the x direction were computed from the trajectories of molecular dynamics simulations. In the absence of the EEF, the forces between the two strands in vacuum were repulsive and attractive in the parallel and antiparallel β sheet structures, respectively. In contrast, the mean forces between the two strands in water were attractive in both the parallel and antiparallel β sheet structures. This is because the electric interactions between the two strands were shielded by water, and the hydrophobic effect dominated the interaction between the two strands. When an EEF >50 MV/cm was applied to the water cluster, the attractive force between the two strands in the parallel and antiparallel β sheet structures decreased and increased, respectively. Further, the binding affinity between the two strands in the parallel and antiparallel β sheet structures also decreased and increased, respectively. This is because the large EEF leads to dielectric saturation, and consequently reduces the effects of the dielectric shielding and hydrophobic interactions. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 861–870, 2014.     

Comparison of the three‐dimensional structures of a humanized and a chimeric Fab of an anti‐γ‐interferon antibody

The objective of this work is to compare the three‐dimensional structures of “humanized” and mouse–human chimeric forms of a murine monoclonal antibody elicited against human γ‐interferon. It is also to provide structural explanations for the small differences in the affinities and biological interactions of the two molecules for this antigen. Antigen‐binding fragments (Fabs) were produced by papain hydrolysis of the antibodies and crystallized with polyethylene glycol (PEG) 8000 by nearly identical microseeding procedures. Their structures were solved by X‐ray analyses at 2.9 Å resolution, using molecular replacement methods and crystallographic refinement. Comparison of these structures revealed marked similarities in the light (L) chains and near identities of the constant (C) domains of the heavy (H) chains. However, the variable (V) domains of the heavy chains exhibited substantial differences in the conformations of all three complementarity‐determining regions (CDRs), and in their first framework segments (FR1). In FR1 of the humanized V_H, the substitution of serine for proline in position 7 allowed the N‐terminal segment (designated strand 4‐1) to be closely juxtaposed to an adjacent strand (4‐2) and form hydrogen bonds typical of an antiparallel β‐pleated sheet. The tightening of the humanized structure was relayed in such a way as to decrease the space available for the last portion of HFR1 and the first part of HCDR1. This compression led to the formation of an α‐helix involving residues 25–32. With fewer steric constraints, the corresponding segment in the chimeric Fab lengthened by at least 1 Å to a random coil which terminated in a single turn of 3₁₀ helix. In the humanized Fab, HCDR1, which is sandwiched between HCDR2 and HCDR3, significantly influenced the structures of both regions. HCDR2 was forced into a bent and twisted orientation different from that in the chimeric Fab, both at the crown of the loop (around proline H52a) and at its base. As in HCDR1, the last few residues of HCDR2 in the humanized Fab were compressed into a space‐saving α‐helix, contrasting with a more extended 3₁₀ helix in the chimeric form. HCDR3 in the humanized Fab was also adjusted in shape and topography. The observed similarities in the functional binding activities of the two molecules can be rationalized by limited induced fit adjustments in their structures on antigen binding. While not perfect replicas, the two structures are testimonials to the progress in making high affinity monoclonal antibodies safe for human use. Copyright © 1999 John Wiley & Sons, Ltd.     

Self‐assembly of regenerated silk fibroin from random coil nanostructures to antiparallel β‐sheet nanostructures

In this work, we studied the effects of incubation concentration and time on the self‐assembly behaviors of regenerated silk fibroin (RSF). Our results showed the assembly ways of RSF were concentration‐dependent and there were four self‐assembly ways of RSF: (i) At relatively low concentration (≤0.015%), RSF molecules assembled into protofilaments (random coil), and then the thickness decreased and the secondary conformation changed to antiparallel β‐sheet; (ii) at the concentration of 0.015%, RSF molecules assembled into protofilaments (random coil), and then assembled into protofibrils (antiparallel β‐sheet). The protofibrils experienced the appearance and disappearance of phase periodic intervals in turn; (iii) at the concentration of 0.03%, RSF molecules assembled into bead‐like oligomers (random coil), and then assembled into protofibrils (antiparallel β‐sheet), and finally the height and phase periodic intervals of RSF protofibrils disappeared in turn; and (iv) at the relatively high concentration (≥0.15%), RSF molecules assembled into protofilaments (random coil), then aggregated into blurry cuboid‐like micelles (random coil), and finally self‐arranged to form smooth and clear cuboid‐like micelles (antiparallel β‐sheet). These results provide useful insights into the process by which the RSF molecules self‐assemble into protofilaments, protofibrils and micelles. Furthermore, our work will be beneficial to basic understanding of the nanoscale structure formations in different silk‐based biomaterials. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 1181–1192, 2014.     

Helix–coil stability constants for the naturally occurring amino acids in water. XXI. Glutamine parameters from random poly(hydroxypropylglutamine-co-L-glutamine) and poly(hydroxybutylglutamine-co-L-glutamine)

Water-soluble, random copolymers containing L -glutamine and either N⁵-(3-hydroxypropyl)-L -glutamine or N⁵-(4-hydroxybutyl)-L -glutamine were synthesized, fractionated, and characterized. The thermally induced helix–coil transitions of these copolymers were studied in water. A short-range interaction theory was used to deduce the Zimm-Bragg parameters σ and s for the helix–coil transition in poly(L -glutamine) in water from an analysis of the melting curves of the copolymers in the manner described in earlier papers. The computed values of s indicate that L -glutamine is helix-indifferent at low temperature and a helix-destabilizing residue at high temperature in water. At all temperatures in the range of 0–70°C, the glutamine residue promotes helix–coil boundaries since the computed value of σ is large.     

Synthesis and conformational study of poly(0,0′-dicarbobenzoxy-L-β-3,4-dihydroxyphenyl-α-alanine)

High-molecular-weight poly(0,0′-dicarbobenzoxy-L -β-3,4-dihydroxyphenyl-α-alanine) was prepared by the N-carboxyanhydride method. From the results obtained by a study of the optical rotation, nuclear magnetic resonance, and solution infrared absorption, the conformation of poly(0,0′-dicarbobenzoxy-L -β-3,4-dihydroxyphenyl-α-alanine) depended greatly on the solvent taking a right-handed helix with [θ]₂₂₅ = ?13,600 ～ ?18,900 in alkyl halides, a left-handed helix with [θ]₂₂₈ = 22,100 ～ 24,800 in cyclic ethers or trimethylphosphate, and a random coil structure in dichloroacetic acid, trifluoroacetic acid, or hexafluoroacetone sesquihydrate. The polypeptide underwent a right-handed helix-coil transition in chloroform/dichloroacetic acid (or trifluoroacetic acid) mixed solvents and a left-handed helix-coil transition in dioxane/dichloroacetic acid (or trifluoroacetic acid) mixed solvents. The results were compared with those of poly(0-carbobenzoxy-L -tyrosine).     

Slow formation of aggregation‐resistant β‐sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding‐related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large β‐helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all‐β‐sheet protein allows detailed analysis of the formation of β‐sheet structure in larger proteins. Using a combination of fluorescence and far‐UV circular dichroism spectroscopy, we show that the pertactin β‐helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and β‐sheet‐rich topology, pertactin refolding is reversible and not complicated by off‐pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate‐limiting step. Furthermore, site‐specific labeling experiments indicate that the β‐helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, β‐sheet‐rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, β‐sheet‐rich refolding intermediates. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Alanine substitutions of noncysteine residues in the cysteine‐stabilized αβ motif

The protein scaffold is a peptide framework with a high tolerance of residue modifications. The cysteine‐stabilized αβ motif (CSαβ) consists of an α‐helix and an antiparallel triple‐stranded β‐sheet connected by two disulfide bridges. Proteins containing this motif share low sequence identity but high structural similarity and has been suggested as a good scaffold for protein engineering. The Vigna radiate defensin 1 (VrD1), a plant defensin, serves here as a model protein to probe the amino acid tolerance of CSαβ motif. A systematic alanine substitution is performed on the VrD1. The key residues governing the inhibitory function and structure stability are monitored. Thirty‐two of 46 residue positions of VrD1 are altered by site‐directed mutagenesis techniques. The circular dichroism spectrum, intrinsic fluorescence spectrum, and chemical denaturation are used to analyze the conformation and structural stability of proteins. The secondary structures were highly tolerant to the amino acid substitutions; however, the protein stabilities were varied for each mutant. Many mutants, although they maintained their conformations, altered their inhibitory function significantly. In this study, we reported the first alanine scan on the plant defensin containing the CSαβ motif. The information is valuable to the scaffold with the CSαβ motif and protein engineering.     

Carbon (δ13C) and nitrogen (δ15N) stable isotope composition in plant and soil in Southern Patagonia's native forests

Stable isotope natural abundance measurements integrate across several biogeochemical processes in ecosystem N and C dynamics. Here, we report trends in natural isotope abundance (δ¹³C and δ¹⁵N in plant and soil) along a climosequence of 33 Nothofagus forest stands located within Patagonia, Southern Argentina. We measured 28 different abiotic variables (both climatic variables and soil properties) to characterize environmental conditions at each of the 33 sites. Foliar δ¹³C values ranged from ?35.4‰ to ?27.7‰, and correlated positively with foliar δ¹⁵N values, ranging from ?3.7‰ to 5.2‰. Soil δ¹³C and δ¹⁵N values reflected the isotopic trends of the foliar tissues and ranged from ?29.8‰ to ?25.3‰, and ?4.8‰ to 6.4‰, respectively, with no significant differences between Nothofagus species (Nothofagus pumilio, Nothofagus antarctica, Nothofagus betuloides). Principal component analysis and multiple regressions suggested that mainly water availability variables (mean annual precipitation), but not soil properties, explained between 42% and 79% of the variations in foliar and soil δ¹³C and δ¹⁵N natural abundance, which declined with increased moisture supply. We conclude that a decline in water use efficiency at wetter sites promotes both the depletion of heavy C and N isotopes in soil and plant biomass. Soil δ¹³C values were higher than those of the plant tissues and this difference increased as annual precipitation increased. No such differences were apparent when δ¹⁵N values in soil and plant were compared, which indicates that climatic differences contributed more to the overall C balance than to the overall N balance in these forest ecosystems.     

Structural versatility of peptides from Cα, α-disubstituted glycines: Crystal-state conformational analysis of homopeptides from Cα-methyl,Cα-benzylglycine [(αMe)Phe]n

The molecular and crystal structures of one derivative and three homopeptides (from the di-to the tetrapeptide level) of the chiral, C^{α, α}-disubstituted glycine C^α-methyl, C^α-benzylglycine [(αMe)Phe], have been determined by x-ray diffraction. The derivative is mClAc-D -(αMe)Phe-OH, and the peptides are pBrBz-[D -(αMe)Phe]₂-NHMe, pBrBz-[D -(αMe)Phe]₃-OH hemihydrate, and pBrBz-[D -(αMe)Phe]₄-OtBu sesquihydrate. All (αMe)Phe residues prefer ?,ψ torsion angles in the helical region of the conformational map. The dipeptide methylamide and the tripeptide carboxylic acid adopt a β-turn conformation with a 1 ← 4 C?O…?H? N intramolecular H bond. The structure of the tripeptide carboxylic acid is further stabilized by a 1 ← 4 C?O…?H? O intramolecular H bond, forming an “oxy-analogue” of a β-turn. The tetrapeptide ester is folded in a regular (incipient) 3₁₀-helix. In general, the relationship between (αMe)Phe chirality and helix screw sense is opposite to that exhibited by protein amino acids. A comparison is made with the conclusions extracted from published work on homopeptides from other C^α-methylated α-amino acids. © 1993 John Wiley & Sons, Inc.     

An amino acid code for β‐sheet packing structure

To understand the relationship between protein sequence and structure, this work extends the knob‐socket model in an investigation of β‐sheet packing. Over a comprehensive set of β‐sheet folds, the contacts between residues were used to identify packing cliques: sets of residues that all contact each other. These packing cliques were then classified based on size and contact order. From this analysis, the two types of four‐residue packing cliques necessary to describe β‐sheet packing were characterized. Both occur between two adjacent hydrogen bonded β‐strands. First, defining the secondary structure packing within β‐sheets, the combined socket or XY:HG pocket consists of four residues i, i+2 on one strand and j, j+2 on the other. Second, characterizing the tertiary packing between β‐sheets, the knob‐socket XY:H+B consists of a three‐residue XY:H socket (i, i+2 on one strand and j on the other) packed against a knob B residue (residue k distant in sequence). Depending on the packing depth of the knob B residue, two types of knob‐sockets are found: side‐chain and main‐chain sockets. The amino acid composition of the pockets and knob‐sockets reveal the sequence specificity of β‐sheet packing. For β‐sheet formation, the XY:HG pocket clearly shows sequence specificity of amino acids. For tertiary packing, the XY:H+B side‐chain and main‐chain sockets exhibit distinct amino acid preferences at each position. These relationships define an amino acid code for β‐sheet structure and provide an intuitive topological mapping of β‐sheet packing. Proteins 2014; 82:2128–2140. © 2014 Wiley Periodicals, Inc.     

Theory of protein molecule self-organization. II. A comparison of calculated thermodynamic parameters of local secondary structures with experiments

Constants of the helix–coil transition for all natural amino acid residues are evaluated on the basis of thermodynamic parameters obtained in paper I of this series. The specific effects at the termini of the helices are also considered as well as the parameters controlling the formation of β-bends in the unfolded protein chain. Evaluated s constants of the helix–coil transition agree with the experimental data on helix–coil transitions of synthetic polypeptides in water. Only a very qualitative correlation exists between s constants (both experimental and theoretical) and the occurrence of corresponding residues in internal turns of α-helices in globular proteins: residues with s > 1 occur in helices as a rule more often than residues with s < 1. At the same time a direct correlation is demonstrated between theoretical parameters of residue incorporation into α-helical termini and β-bends in an unfolded polypeptide chain and the occurrence of residues in corresponding positions of the globular protein secondary structures.     

Optical properties of polypeptides in the β-conformation

The rotational strengths and oscillator strengths of the nπ* band and ππ* exciton bands have been calculated for antiparallel and parallel β-structures of varying length and width. The results are compared with experiment and with previous theoretical treatments of β-structures. The generally good agreement of calculations on the antiparallel β-structure with experimental results on poly-L -lysine and poly-L -serine indicates that these systems are indeed in the antiparallel conformation. It is found that the exciton component strongest in absorption shifts to longer wavelengths as the width of an antiparallel structure increases, and it is suggested that the position of the ππ* absorption band may be a useful criterion of sheet width. The results also reconcile the linear dichroism measurements of Rosenheck and Sommer on poly-L -lysine films with an anti-parallel structure. Calculations on parallel β-structures indicate that the CD spectra of this form will be rather similar to that of the antiparallel form. However, the major absorption band in the antiparallel form is associated with a small positive CD band, while in the parallel form it coincides with a large negative CD band. Finally, it is pointed out that the large positive CD bands predicted for single-stranded parallel and antiparallel β-structures at about 200 mμ render unlikely the suggestion that random-coil polypeptides contain a substantial fraction of extended chain.     

Large scale conformational transitions in β‐structural motif of gramicidin A: kinetic analysis based on CD and FT‐IR data

Gramicidin A (gA) is a polypeptide antibiotic, which forms dimeric channels specific for monovalent cations in artificial and biological membranes. It is a polymorphic molecule that adopts a unique variety of helical conformations, including antiparallel double‐stranded ↑↓β5.6 or ↑↓β7.2 helices (number of residues per turn) and a single‐stranded β6.3 helix (the ‘channel form’). The behavior of gA‐Cs⁺ complex in the micelles of TX‐100 was studied in this work. Transfer of the complex into the micelles activates a cascade of sequential conformational transitions monitored by CD and FT‐IR spectroscopy: At the first step after Cs⁺ removal, the RH ↑↓β5.6 helix is formed, which has been discussed so far only hypothetically. Kinetics of the transitions was measured, and the activation parameters were determined. The activation energies of the ↑↓β5.6 → β‐helical monomer transition in dioxane and dioxane/water solutions were also measured for comparison. The presence of water raises the transition rate constant ~10³ times but does not lead to crucial fall of the activation energy. All activation energies were found in the 20–25 kcal/mol range, i.e. much lower than would be expected for unwinding of the double helix (when 28 H‐bonds are broken simultaneously). These results can be accounted for in the light of local unfolding (or ‘cracking’) model for large scale conformational transitions developed by the P. G.Wolynes team [Miyashita O, Onuchic JN, Wolynes PG. Proc. Natl. Acad. Sci. USA 2003; 100: 12570‐12575.]. Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.     

Reconstructing the δ18O of atmospheric water vapour via the CAM epiphyte Tillandsia usneoides: seasonal controls on δ18O in the field and large‐scale reconstruction of δ18Oa

Using both oxygen isotope ratios of leaf water (δ¹⁸O_L) and cellulose (δ¹⁸O_C) of Tillandsia usneoides in situ, this paper examined how short‐ and long‐term responses to environmental variation and model parameterization affected the reconstruction of the atmospheric water vapour (δ¹⁸O_a). During sample‐intensive field campaigns, predictions of δ¹⁸O_L matched observations well using a non‐steady‐state model, but the model required data‐rich parameterization. Predictions from the more easily parameterized maximum enrichment model (δ¹⁸O_L–M) matched observed δ¹⁸O_L and observed δ¹⁸O_a when leaf water turnover was less than 3.5 d. Using the δ¹⁸O_L–M model and weekly samples of δ¹⁸O_L across two growing seasons in Florida, USA, reconstructed δ¹⁸O_a was ?12.6 ± 0.3‰. This is compared with δ¹⁸O_a of ?12.4 ± 0.2‰ resolved from the growing‐season‐weighted δ¹⁸O_C. Both of these values were similar to δ¹⁸O_a in equilibrium with precipitation, ?12.9‰. δ¹⁸O_a was also reconstructed through a large‐scale transect with δ¹⁸O_L and the growing‐season‐integrated δ¹⁸O_C across the southeastern United States. There was considerable large‐scale variation, but there was regional, weather‐induced coherence in δ¹⁸O_a when using δ¹⁸O_L. The reconstruction of δ¹⁸O_a with δ¹⁸O_C generally supported the assumption of δ¹⁸O_a being in equilibrium with precipitation δ¹⁸O (δ¹⁸O_ppt), but the pool of δ¹⁸O_ppt with which δ¹⁸O_a was in equilibrium – growing season versus annual δ¹⁸O_ppt – changed with latitude.