Quantitative and qualitative analysis of type III antifreeze protein structure and function

Some cold water marine fishes avoid cellular damage because of freezing by expressing antifreeze proteins (AFPs) that bind to ice and inhibit its growth; one such protein is the globular type III AFP from eel pout. Despite several studies, the mechanism of ice binding remains unclear because of the difficulty in modeling the AFP-ice interaction. To further explore the mechanism, we have determined the x-ray crystallographic structure of 10 type III AFP mutants and combined that information with 7 previously determined structures to mainly analyze specific AFP-ice interactions such as hydrogen bonds. Quantitative assessment of binding was performed using a neural network with properties of the structure as input and predicted antifreeze activity as output. Using the cross-validation method, a correlation coefficient of 0.60 was obtained between measured and predicted activity, indicating successful learning and good predictive power. A large loss in the predictive power of the neural network occurred after properties related to the hydrophobic surface were left out, suggesting that van der Waal's interactions make a significant contribution to ice binding. By combining the analysis of the neural network with antifreeze activity and x-ray crystallographic structures of the mutants, we extend the existing ice-binding model to a two-step process: 1) probing of the surface for the correct ice-binding plane by hydrogen-bonding side chains and 2) attractive van der Waal's interactions between the other residues of the ice-binding surface and the ice, which increases the strength of the protein-ice interaction.

Alternative roles for putative ice-binding residues in type I antifreeze protein

Two sets of variants of type I antifreeze protein have been synthesized to investigate the role of Leu and Asn in the activity of this 37-residue alpha-helix. Leu and Asn flank the central two of four regularly spaced ice-binding Thr in the i-1 and i + 3 positions, respectively. All three residues project from the same side of the helix to form the protein's putative ice-adsorption site and are considered in some models to act together as an "ice-binding motif". Replacement of Asn by residues with shorter side chains resulted in either a small loss (Ala) or gain (Thr) of antifreeze activity. However, substitution of Asn by its slightly larger homologue (Gln) abolished thermal hysteresis activity. The Gln-containing peptide was very soluble, largely monomeric, and fully helical. Of the three variants in which Leu was replaced by Ala, two of the three were more active than their Leu-containing counterparts, but all three variants began to precipitate as the peptide concentration increased. None of the seven variants tested showed dramatic differences in ice crystal morphology from that established by the wild type. These results are consistent with a primary role for Leu in preventing peptide aggregation at the antifreeze protein concentrations (10 mg/mL) normally present in fish serum. Similarly the role for Asn may have more to do with enhancing the solubility of these rather hydrophobic peptides than of making a stereospecific hydrogen-bonding match to the ice lattice as traditionally thought. Nevertheless, the dramatic loss of activity in the Asn-to-Gln replacement demonstrates the steric restriction on residues in or near the ice-binding site of the peptide.

Type II fish antifreeze protein accumulation in transgenic tobacco does not confer frost resistance

Type II fish antifreeze protein (AFP) is active in both freezing point depression and the inhibition of ice recrystallization. This extensively disulfide-bonded 14 kDa protein was targeted for accumulation in its pro- and mature forms in the cytosol and apoplast of transgenic tobacco plants. Type II AFP gene constructs under control of a duplicate cauliflower mosaic virus 35S promoter, both with and without a native plant transit peptide sequence, were introduced into tobacco by Agrobacterium tumefaciens-mediated transformation. AFP did not accumulate in the cytosol of transgenic plants, but active AFP was present as 2% the total protein present in the apoplast. Plant-produced AFP was the same size as mature Type II AFP isolated from fish, and was comparable to wild-type AFP in thermal hysteresis activity and its effect on ice crystal morphology. Field trials conducted in late summer on R1 generation transgenic plants showed similar AFP accumulation in plants under field conditions at levels suitable for large-scale production: but no difference in frost resistance was observed between transgenic and wild-type plants during the onset of early fall frosts.

Recent and rapid amplification of the sperm basic nuclear protein genes in winter flounder

The high molecular weight basic nuclear proteins (HMrBNPs), which are tightly bound to sperm chromatin in winter flounder, are made up of imperfect reiterations of simple peptide sequences that contain phosphorylatable DNA-binding motifs. Genomic Southern blots hybridized with probes to the coding and non-coding regions of HMrBNP mRNA showed that HMrBNP sequences form a complex multi-gene family. Previously, one gene (2B) was used to establish an evolutionary link between histone H1 and the HMrBNPs. Further examination of this complex, multi-gene family has now revealed that the majority of the HMrBNP genes are linked as 4.5 kb direct tandem repeats that each contain a 2.8 kb coding region and a 1.7 kb intergenic region (IR). These findings, combined with the cloning of the IR, established that the tandemly repeated genes lack introns and code for the abundant 3 kb HMrBNP mRNAs that produce the prominent 110 kDa HMrBNP. Southern blotting of DNAs from other righteye flounder species showed that HMrBNP multi-gene families were present in closely related species, though with substantial differences in restriction patterns and band intensities, but were not detected in more distantly related flounders. These observations are consistent with recent and rapid elaboration of the HMrBNP gene family.

Increased gene dosage augments antifreeze protein levels in transgenic Drosophila melanogaster

One of the principal environmental adaptations of certain fishes inhabiting polar and northern coastal waters is the synthesis of antifreeze proteins (AFPs). AFPs bind to and prevent the growth of nascent ice crystals, thus depressing the serum freezing point. The transgenic expression of AFP holds great promise for conferring freeze resistance to commercially important plant and animal species. Since fish at the greatest risk of freezing have multiple AFP gene copies in order to synthesize higher levels of this protein, we have evaluated this evolutionary strategy as a way to maximize AFP expression in a model transgenic host, the fruit fly Drosophila melanogaster. A construct in which AFP genes of the Atlantic wolffish are fused to the Drosophila yolk protein 1,2 promoter/enhancer region was transferred to flies through P-element mediated transformation. Several independent transgenic fly lines were used in genetic crosses to obtain multi-insert lines. Haemolymph freezing point depression (thermal hysteresis) was greater in homozygotes relative to heterozygotes for a given insert. Similarly, multi-insert lines consistently displayed greater haemolymph AFP activity than the single insert lines from which they were derived. The thermal hysteresis value obtained with a fly line harboring 8 AFP gene copies, 0.43 degree C, represents the highest such value to date recorded in a transgenic host, and is even higher than the levels found in some AFP-producing fish.

NMR structural studies on antifreeze proteins

Antifreeze proteins (AFPs) are a structurally diverse class of proteins that bind to ice and inhibit its growth in a noncolligative manner. This adsorption-inhibition mechanism operating at the ice surface results in a lowering of the (nonequilibrium) freezing point below the melting point. A lowering of approximately 1 degree C, which is sufficient to prevent fish from freezing in ice-laden seawater, requires millimolar AFP levels in the blood. The solubility of AFPs at these millimolar concentrations and the small size of the AFPs (typically 3-15 kDa) make them ideal subjects for NMR analysis. Although fish AFPs are naturally abundant, seasonal expression, restricted access to polar fishes, and difficulties in separating numerous similar isoforms have made protein expression the method of choice for producing AFPs for structural studies. Expression of recombinant AFPs has also facilitated NMR analysis by permitting isotopic labeling with 15N and 13C and has permitted mutations to be made to help with the interpretation of NMR data. NMR analysis has recently solved two AFP structures and provided valuable information about the disposition of ice-binding side chains in a third. The potential exists to solve other AFP structures, including the newly described insect AFPs, and to use solid-state NMR techniques to address fundamental questions about the nature of the interaction between AFPs and ice.

The ice-binding site of sea raven antifreeze protein is distinct from the carbohydrate-binding site of the homologous C-type lectin

Antifreeze proteins lower the freezing point of their solution by binding to ice and inhibiting its growth. One of several structurally different antifreeze proteins in fishes (type II) is homologous to the carbohydrate-recognition domain of Ca2+-dependent lectins and adopts the same three-dimensional fold. Type II antifreeze proteins from herring and smelt require Ca2+ for binding to ice, whereas this same antifreeze protein in sea raven binds to ice in the absence of Ca2+ and has only two of the five Ca2+-liganding amino acids that are present in the lectin. To locate the ice-binding site, site-directed mutants of the 15 kDa, globular, disulfide-bonded sea raven antifreeze protein were produced by secretion from Pichia pastoris. Pairs of amino acid replacements, insertions, and a peptide loop swap were made in the region equivalent to the sugar-binding site of the lectin that encompasses loops 3 and 4 and beta-sheets 7 and 8. Even the most extensive mutation caused only a 25% decrease in antifreeze activity and demonstrated that the residues corresponding to the Ca2+-binding site are only peripherally involved in ice binding. When adjacent surface residues were mutated, the replacement of one residue, Ser120 by His, caused a 35% decrease in activity by itself and an 80% loss in conjunction with the peptide loop swap mutation. This pivotal sea raven antifreeze protein amino acid does not coincide with the herring ice-binding epicenter, but is located within the region corresponding to the proposed CaCO3-binding surface of a third homologue, the pancreatic stone protein. Intron and exon structure of the sea raven AFP gene also suggests that it might be more closely related to the stone protein gene than to the lectin gene. These results support the notion that this family of proteins has evolved more than one binding surface from the same protein scaffold.

Disulfide bond mapping and structural characterization of spruce budworm antifreeze protein

The 9-kDa, Thr-, Ser-, and Cys-rich thermal hysteresis protein from spruce budworm (sbwTHP) is 10-30 times more effective than fish antifreeze proteins (AFPs) at depressing solution freezing points via ice-crystal growth inhibition. Since this insect protein is only available in microgram quantities from its natural source, recombinant sbwTHP was produced from inclusion bodies in Escherichia coli by a refolding protocol. Incompletely folded forms were removed during ion-exchange and reverse-phase chromatography, resulting in fully active sbwTHP that was indistinguishable in its properties from native sbwTHP. The antifreeze was completely inactivated by reduction, showed no reaction with sulfhydryl reagents, and was not inhibited by EDTA. All eight cysteine residues appear to be involved in disulfide bond formation. Tryptic cleavage and peptide analysis is consistent with linkages between the first and second cysteine residues, the third and fourth, fifth and eighth, and the sixth and seventh. NMR analysis confirmed that the fully folded form of sbwTHP was well structured and had a single conformation. Both NMR and CD spectra indicate the presence of extensive beta structure (70-80%) with little or no alpha helix. The protein maintains antifreeze activity over a broad range of pH values, and its conformation is independent of both temperature (over the range 0 degrees C to 20 degrees C), and the presence of 50% trifluoroethanol.

Binding of an oligopeptide to a specific plane of ice

The alpha-helical antifreeze protein (AFP) from winter flounder inhibits ice growth by binding to a specific set of pyramidal surface planes that are not otherwise macroscopically expressed. The 37-residue AFP contains three 11-amino acid repeats that make a stereo-specific fit to the ice lattice along the <01-12> direction of the (20-21) and equivalent binding planes. When the AFP was shortened to delete two of the three 11-amino acid ice-binding repeats, the resulting 15-residue peptide and its variants were less helical and showed no antifreeze activity. However, when the helicity of the peptide was reinforced by an internal lactam bridge between Glu-7 and Lys-11, the minimized AFP was able to stably express the pyramidal plane (20-21) on the surface of growing ice crystals. This dynamic shaping of the ice surface by a single ice-binding repeat provides evidence that AFP adsorption to the ice lattice is not an "all-or-nothing" interaction. Instead, a partial interaction can help develop the binding site on ice to which the remainder of the AFP (or other AFP molecules) can orient and bind.

The solution structure of type II antifreeze protein reveals a new member of the lectin family

A recombinant form of the sea raven type II antifreeze protein (SRAFP) has been produced using the Pichia pastoris expression system. The antifreeze activity of recombinant SRAFP is indistinguishable from that of the wild-type protein. The global fold of SRAFP has been determined by two-dimensional 1H homonuclear and three-dimensional 1H-¿15N¿ heteronuclear NMR spectroscopy using 785 NOE distance restraints and 47 angular restraints. The molecule folds into one globular domain that consists of two helices and nine beta-strands in two beta-sheets. The structure confirms the proposed existence of five disulfide bonds. The global fold of SRAFP is homologous to C-type lectins and pancreatic stone proteins, even though the sequence identity is only approximately 20%.

The high molecular weight chromatin proteins of winter flounder sperm are related to an extreme histone H1 variant

Unlike mammals, birds, and most other fishes, winter flounder completes spermatogenesis without replacing its germ cell histones with protamines. Instead, during spermiogenesis, these fish produce a family of high molecular weight (80,000-200,000) basic nuclear proteins (HMrBNPs) that bind to sperm chromatin containing the normal complement of histones. These large, basic proteins are built up of tandem iterations of oligopeptide repeats that contain phosphorylatable DNA-binding motifs. Although the HMrBNPs have no obvious homology to histones, protamines, or other sperm-specific chromatin proteins, we report here the isolation of a clone (2B) from a winter flounder genomic DNA library that establishes a link between the HMrBNPs and histone H1. The 2B sequence contains an open reading frame, which, when conceptually translated, encodes a 265-residue protein. At its N terminus the translation product contains numerous simple repeats that match the oligopeptides contained within the HMrBNPs. Unexpectedly, the C terminus of the putative protein shows 66% identity and 76% conservation to the histone H1 globular domain. This connection suggests that the HMrBNPs may have originated from the extended N-terminal tail region of a testis-specific, H1-like linker histone.

Antifreeze proteins bind independently to ice

It has been suggested that cooperative interactions between antifreeze proteins (AFPs) on the ice surfaces are required for complete inhibition of ice crystal growth. To test this hypothesis, a 7-kDa type III AFP was linked through its N-terminus to thioredoxin (12 kDa) or maltose-binding protein (42 kDa). The resultant 20-kDa and 50-kDa fusion proteins were larger in diameter than free AFP and thus precluded any extensive AFP-AFP contacts on the ice surface. Both fusion proteins were at least as active as free AFP at virtually all concentrations tested. By these criteria, AFPs function independently of each other and do not require specific intermolecular interactions to bind tightly to ice.

The effects of steric mutations on the structure of type III antifreeze protein and its interaction with ice

The interaction of proteins with ice is poorly understood and difficult to study, partly because ice is transitory and can present many binding surfaces, and partly because structures have been determined for only two ice-binding proteins. This paper focuses on one of these, a 66-residue antifreeze protein (AFP) from eel pout. The high resolution X-ray structure of this fish AFP demonstrated that the proposed ice-binding surface is remarkably flat for such a small protein. The residues on the planar surface thought to be involved in ice binding are restrained by hydrogen bonds or by tight packing of their side-chains. To probe the requirement for a flat binding surface, a conserved alanine in the center of the AFP planar surface was substituted with larger residues. Six alanine replacement mutants (Ala16 > Cys, Thr, Met, Arg, His and Tyr), designed to disrupt the planarity of the surface and sterically block binding to ice, were characterized by X-ray crystallography and compared with the wild-type AFP. In each case, the detail provided by these crystal structures has helped explain the effects of the mutation on antifreeze activity. The substitutions, Ala16 > His and Ala16 > Tyr, were large enough to shield Gln44, one of the putative ice-binding residues, contributing to their very low thermal hysteresis activity. In addition to sterically hindering the putative ice-binding site, the bulkier residues also caused shifts in the putative ice-binding residues owing to the tight packing of side-chains on the planar surface. This unexpected consequence of the mutations helps account for the severely reduced antifreeze activity. One explanation for residual antifreeze activity in some of the mutants lies in the possibility that AFPs have a role in shaping the site on the ice to which they bind. Thus, side-chain dislocations might be partially accommodated by ice that can freeze around them. It is evident that the disruption of the planarity, by introducing larger residues at the center of the proposed ice-binding site, is not the only factor responsible for the loss of antifreeze activity. There are multiple causes including positional change and steric blockage of some putative ice-binding residues.

A diminished role for hydrogen bonds in antifreeze protein binding to ice

The most abundant isoform (HPLC-6) of type I antifreeze protein (AFP1) in winter flounder is a 37-amino-acid-long, alanine-rich, alpha-helical peptide, containing four Thr spaced 11 amino acids apart. It is generally assumed that HPLC-6 binds ice through a hydrogen-bonding match between the Thr and neighboring Asx residues to oxygens atoms on the {2021} plane of the ice lattice. The result is a lowering of the nonequilibrium freezing point below the melting point (thermal hysteresis). HPLC-6, and two variants in which the central two Thr were replaced with either Ser or Val, were synthesized. The Ser variant was virtually inactive, while only a minor loss of activity was observed in the Val variant. CD, ultracentrifugation, and NMR studies indicated no significant structural changes or aggregation of the variants compared to HPLC-6. These results call into question the role of hydrogen bonds and suggest a much more significant role for entropic effects and van der Waals interactions in binding AFP to ice.

Antifreeze proteins

Antifreeze proteins comprise a structurally diverse class of proteins that inhibit the growth of ice. Recently, new AFP types have been discovered; more active AFPs have been isolated; antecedents have been recognized supporting the notion of recent, multiple origins; and detailed structures have emerged leading to models for their adsorption to ice.

Backbone structure and dynamics of a hemolymph protein from the mealworm beetle Tenebrio molitor

Pheromones play a vital role in the survival of insects and are used for chemical communication between members of the same species by their olfactory system. The selection and transportation of these lipophilic messengers by carrier proteins through the hydrophilic sensillum lymph in the antennae toward their membrane receptors remains the initial step for the signal transduction pathway. A moderately abundant 12.4 kDa hydrophilic protein present in hemolymph from the mealworm beetle Tenebrio molitor is approximately 38% identical to a family of insect pheromone-binding proteins. The backbone structure and dynamics of the 108-residue protein have been characterized using three-dimensional 1H-15N NMR spectroscopy, combined with 15N relaxation and 1H/D exchange measurements. The secondary structure, derived from characteristic patterns of dipolar connectivities between backbone protons, secondary chemical shifts, and homonuclear three-bond JHNH alpha coupling constants, consists of a predominantly disordered N-terminus from residues 1 to 10 and six alpha-helices connected by four 4-7 residue loops and one beta-hairpin structure. The up-and-down arrangement of alpha-helices is stabilized by two disulfide bonds and hydrophobic interactions between amphipathic helices. The backbone dynamics were characterized by the overall correlation time, order parameters, and effective correlation times for internal motions. Overall, a good correlation between secondary structure and backbone dynamics was found. The 15N relaxation parameters T1 and T2 and steady-state NOE values of the six alpha-helices could satisfactorily fit the Lipari-Szabo model. In agreement with their generalized order parameters (> 0.88), residues in helical regions exhibited restricted motions on a picosecond time scale. The stability of this highly helical protein was confirmed by thermal denaturation studies.

Biosynthetic production of type II fish antifreeze protein: fermentation by Pichia pastoris

Sea raven type II antifreeze protein (SRAFP) is one of three different fish antifreeze proteins isolated to date. These proteins are known to bind to the surface of ice and inhibit its growth. To solve the three-dimensional structure of SRAFP, study its ice-binding mechanism, and as a basis for engineering these molecules, an efficient system for its biosynthetic production was developed. Several different expression systems have been tested including baculovirus, Escherichia coli and yeast. The latter, using the methylotrophic organism Pichia pastoris as the host, was the most productive. In shake-flask cultures the levels of SRAFP secreted from Pichia were up to 5 mg/l. The recombinant protein has an identical activity to SRAFP from sea raven serum. In order to increase yields further, four different strategies were tested in 10-l fermentation vessels, including: (1) optimization of pH and dissolved oxygen, (2) mixed feeding of methanol and glycerol with Mut(s) clones, (3) supplementation of amino acid building blocks, and (4) methanol feeding with Mut+ clones. The mixed-feeding/Mut(s) strategy proved to be the most efficient with SRAFP yields reaching 30 mg/l.

The antifreeze potential of the spruce budworm thermal hysteresis protein

Antifreeze proteins (AFP) inhibit ice growth by surface adsorption that results in a depression of the freezing point below the melting point. The maximum level of this thermal hysteresis shown by the four structurally unrelated fish AFP is approximately 1.5 degrees C. In contrast, hemolymph and crude extracts from insects can have 5 degrees to 10 degrees C of thermal hysteresis. Based on the isolation, cloning, and expression of a thermal hysteresis protein (THP) from spruce budworm (Choristoneura fumiferana), the vastly greater activity is attributable to a 9 kDa protein. This novel, threonine- and cysteine-rich THP has striking effects on ice crystal morphology, both before and during freezing. It is also 10 to 30 times more active than any known fish AFP, offering the prospect of superior antifreeze properties in cryoprotective applications.

The effects of truncations of the small subunit on m-calpain activity and heterodimer formation

In order to study subunit interactions in calpain, the effects of small subunit truncations on m-calpain activity and heterodimer formation have been measured. It has been shown previously that active calpain is formed by co-expression of the large subunit (80 kDa) of rat m-calpain with a delta 86 form (21 kDa) of the small subunit. cDNA for the full-length 270 amino acid (28.5 kDa) rat calpain small subunit has now been cloned, both with and without an N-terminal histidine tag (NHis10). The full-length small subunit constructs yielded active calpains on co-expression with the large subunit, and the small subunit was autolysed to 20 kDa on exposure of these calpains to Ca2+. A series of deletion mutants of the small subunit, NHis10-delta 86, -delta 99, -delta 107, and -delta 116, gave active heterodimeric calpains with unchanged specific activities, although in decreasing yield, and with a progressive decrease in stability. NHis10-delta 125 formed a heterodimer which was inactive and unstable. Removal of 25 C-terminal residues from delta 86, leaving residues 87-245, abolished both activity and heterodimer formation. The results show that: (a) generation of active m-calpain in Escherichia coli requires heterodimer formation; (b) small subunit residues between 94 and 116 contribute to the stability of the active heterodimer but do not directly affect the catalytic mechanism; (c) residues in the region 245-270 are essential for subunit binding. Finally, it was shown that an inactive mutant Cys103-->Ser-80k/delta 86 calpain, used in order to preclude autolysis, did not dissociate in the presence of Ca2+, a result which does not support the proposal that Ca(2+)-induced dissociation is involved in calpain activation.

Structure of a calpain Ca(2+)-binding domain reveals a novel EF-hand and Ca(2+)-induced conformational changes

The crystal structure of a Ca(2+)-binding domain (dVI) of rat m-calpain has been determined at 2.3 A resolution, both with and without bound Ca2+. The structures reveal a unique fold incorporating five EF-hand motifs per monomer, three of which bind calcium at physiological calcium concentrations, with one showing a novel EF-hand coordination pattern. This investigation gives us a first view of the calcium-induced conformational changes, and consequently an insight into the mechanism of calcium induced activation in calpain. The crystal structures reveal a dVI homodimer which provides a preliminary model for the subunit dimerization in calpain.

Introns boost transgene expression in Drosophila melanogaster

Using Drosophila as a host, we have examined the effect that the presence of an intron has on the accumulation of a processed transgene mRNA. To provide a model system that was free from major position effects, two fish antifreeze protein (AFP) transgenes were arranged in divergent transcriptional orientation from a central Drosophila yolk protein promoter/enhancer region and introduced into the flies by P-element transformation and/or mobilization. In this way the organization of both the structural genes and the promoter elements mimicked their natural arrangements. When one member of the fish AFP transgene pair had its single 180 bp intron deleted, there was a 2- to 11-fold (average 5-fold) decrease in its mRNA level compared to that generated from the control gene containing the fish AFP intervening sequence. When the deleted intron was replaced by a 70 bp intervening sequence originating from the yolk protein 1 gene, mRNA accumulation was restored to its original level. Even for the streamlined genome of Drosophila, where the intron number and size are generally reduced compared to mammals, the presence of an intervening sequence appears to facilitate mRNA accumulation.

NMR characterization of side chain flexibility and backbone structure in the type I antifreeze protein at near freezing temperatures

The flexibility of the polar side chains in the alpha-helical Type I antifreeze protein (AFP) near the solution freezing temperature was investigated by two-dimensional nuclear magnetic resonance spectroscopy. These experiments were conducted to define the rotameric conformations of the proposed ice-binding groups, threonines and asparagines, in order to probe the molecular mechanism for ice binding. On the basis of the 3J alpha beta 2 NMR coupling constant values of 7.1, 8.5, 8.5, and 6.8 Hz for residues T2, T13, T24, and T35, respectively, it can be calculated that the regularly spaced ice-binding threonines sample many possible rotameric states prior to ice binding. The lack of a dominant side chain rotamer is further corroborated by nuclear Overhauser distance measurements for T13 and T24. N16 and N27, both with 3J alpha beta 2 and 3J alpha beta 3 coupling constants of 8.4 and 4.5 Hz, respectively, show a slight preference for the side chain conformation with a chi 1 of -60 degrees. These data suggest that prior to ice binding the threonine and asparagine side chains are free to rotate and that a unique preformed ice-binding structure in solution is not apparent. These observations do not support the rigid side chain model proposed recently by an X-ray study [Sicheri, F., & Yang, D. S. C. (1995) Nature 375, 427-431].