INTRODUCTION
Living organisms capable of emitting light have been known to mankind since ancient times (
Harvey, 1952;
Lee, 2008). Bioluminescent organisms such as bacteria, fireflies, jellyfish, worms, fungi, and fish, are widely dispersed on the phylogenetic tree, with a vast majority of species being marine inhabitants. Believed to emerge “independently” many times during evolution, bioluminescence serves vital functions ranging from defense to reproduction, yet in many cases its survival value remains a puzzle (
Haddock et al., 2010). The bioluminescence is an enzymatic reaction, where an enzyme, generically referred to as luciferase, catalyzes oxidation of a substrate, generically named luciferin. Both luciferases and luciferins of luminous species belonging to different groups generally display unrelated structures. All known luciferin-luciferase reactions share one common feature, a requirement of oxygen. The chemical mechanism involves oxidation of the luciferin to an electronically excited product, which relaxes to its ground state with radiation of light (
Harvey, 1952;
Shimomura, 2006). Conversion of the reaction energy into photons goes with a remarkable efficiency.
In addition to luciferases many bioluminescence systems contain supplemental proteins which can be divided into two groups on the basis of their function: “accessory” proteins provide luciferase with luciferin, and “antenna” proteins modulate the bioluminescence color by means of energy transfer. One type of accessory protein is a luciferin-binding protein, which binds luciferin noncovalently in a protein cavity where it is protected from autooxidation by solubilized oxygen. Proper delivery of this luciferin to the luciferase is accomplished through conformational changes of the binding protein to render the luciferin available for oxidation in the active site of luciferase. Conformational changes of the binding protein being sensitive to Ca
2+ concentration as in the case of
Renilla, or pH value as in the case of dinoflagellates, allow placing the luminescence reactions under control of the neural system (
Hastings and Morin, 1969;
Anderson et al., 1974;
Fogel and Hastings, 1972;
Morse et al., 1989). Another accessory protein, NADPH:FMN oxidoreductase, supplies FMNH
2 for reaction with oxygen and aldehyde in the active site of bacterial luciferase (
Jeffers et al., 2003).
The function of bioluminescence color modulation belongs to the so-called antenna proteins, such as green-fluorescent proteins (GFP) in the bioluminescent coelenterates (
Morin and Hastings, 1971;
Ward and Cormier, 1979), or lumazine protein and yellow-fluorescent protein in bioluminescent bacteria (
O’Kane et al., 1985;
Daubner et al., 1987;
Macheroux et al., 1987). Fluorescence properties of antenna proteins are due to a bound fluorophore, a molecule known to be involved in intermediary metabolism (bacterial antenna proteins), or a moiety autocatalytically synthesized from amino acid residues of the apo-protein (GFP). The radiative transitions of the antenna protein and the luciferin product of the light reaction undergo weak dipole-dipole coupling of the Förster type, and consequently the antenna protein radiates the energy as its own fluorescence, resulting in a shift of the bioluminescence spectral maximum and change of quantum yield of the bioluminescence reaction. Antenna proteins broaden the range of colors utilized by bioluminescent organisms.
Several lines of evidence indicate that delivery of luciferin from the binding protein to the luciferase as well as modulation of bioluminescence spectra as a result of energy transfer to the antenna proteins, involves protein-protein interactions (
Morise et al., 1974;
Ward and Cormier, 1976,
1978;
Nicolas et al., 1991;
Cutler, 1995;
Schultz et al., 2005;
Titushin et al., 2008,
2010;
Stepanyuk et al., 2009;
Markova et al., 2010). The interactions are of a transient nature, with K
D in the range of 10
−6–10
−3 mol/L, which highly impedes direct detection of the complexes, or their structural studies with X-ray crystallography or NMR methods. For those bioluminescent systems where crystal structures of the individual proteins are known, the topologies of the protein-protein complexes were solved by means of data-driven computer modeling. Among the structures available so far, there are those of the lumazine protein-luciferase complex from the bacteria
Photobacterium (
Sato et al., 2010), the GFP-photoprotein clytin complex from the jellyfish
Clytia (
Titushin et al., 2010), and the coelenterazine-binding protein-luciferase complex from the soft coral
Renilla (
Stepanyuk et al., 2009). The main goal of this review is to cover the mentioned bioluminescence systems, where extensive biochemical studies of the protein-protein interactions have been complemented with structural data, resulting in new insight into the bioluminescence mechanism with regard to the role of the protein-protein complexes. It is also emphasized that computer modeling of the complexes, even being based on comprehensive biophysical data, defines the scope of reliability of the structures, necessitating further structural studies.
BACTERIA
The bioluminescent bacteria are the most abundant and widely distributed of all bioluminescent organisms and are found in marine, freshwater, and terrestrial environments. These bacteria emit continuous light and are usually found symbiotically associated with a host (fish, squid), commonly as a monoculture contained in a light organ (
Hastings and Nealson, 1977;
Wilson and Hastings, 1998;
Haddock et al., 2010). Bacterial luciferase, one of the most extensively studied flavin monooxygenases, catalyzes the reaction of FMNH
2, O
2, and an aliphatic aldehyde to form a high-energy intermediate, which decomposes to FMN and the corresponding carboxylic acid with the emission of blue-green light, a broad spectrum with an emission maximum in the range of 490–509 nm depending on the type of luciferase (Fig. 1A and 1D). Luciferase is provided with reduced FMN by an accessory NADPH:FMN oxidoreductase, with which it forms a 1:1 protein-protein complex (
Jeffers et al., 2003).
The emission spectra of some bioluminescent bacteria are found blue or red shifted from that of the purified luciferase
in vitro reaction, due to accompanying antenna proteins (Fig. 1B and 1C) (
Daubner et al., 1987;
Lee, 1993). These proteins, lumazine protein and yellow-fluorescent protein (YFP), contain a highly fluorescent ligand (6,7-dimethyl-8-ribityl-lumazine, FMN, or riboflavin), and apparently associate with a high-energy intermediate in the luciferase reaction so that the excitation step involves a Förster-type coupling to the antenna transition, and the emission becomes that of the bound lumazine or flavin (
Petushkov and Lee, 1997). It is such protein-protein interactions from the structural perspective that are the subject of the present review.
Structure of bacterial luciferase
Bacterial luciferase (PDB code 3FGC) is a heterodimer of homologous α- (~40 kDa) and β- (~35 kDa) subunits, both of which assume the TIM (β/α)
8 barrel fold, characteristic of the bacterial luciferase family (
Fisher et al., 1995,
1996;
Li et al., 2008;
Campbell et al., 2009). The catalytic site resides exclusively on the α-subunit. Binding of FMN occurs in a large open cavity accessible to solvent via an opening located at the C-terminal ends of the β-strands of the TIM-barrel structure (Fig. 1E). One wall of the cavity is lined predominantly with hydrophobic residues, while the other side is lined with charged polar residues, which places flavin in the right orientation, with the isoalloxazine ring facing the hydrophobic surface and the quininoid portion being proximal to the charged polar surface (Fig. 1F). The isoalloxazine ring of the FMN is held in place almost entirely through backbone contacts, including that with the
cis-A74-A75f peptide bond of the β3-strand important in maintaining the right conformation (
Fisher et al., 1996;
Campbell et al., 2009). The 5′-phosphate binding site comprises the side chains of R107, R125, E175, S176, and T179, and the backbone amide of E175 (Fig. 1F) (
Campbell et al., 2009). Bacterial luciferase is an uncommon flavoprotein for it employs FMNH
2 as a substrate rather than a cofactor. Oxidation of reduced flavin and long chain aldehyde by molecular oxygen provides a reaction energy of 60–80 kcal/mol (
Hastings et al., 1973;
Baldwin and Ziegler, 1992), sufficient to populate the excited state of a highly fluorescent luciferase-intermediate, the origin of the bioluminescence emission having a maximum around 490 nm. The impact of the β-subunit is thought to stabilize the active form of the α-subunit (
Baldwin and Ziegler, 1992;
Sinclair et al., 1993;
Sparks and Baldwin, 2001;
Inlow and Baldwin, 2002), and is achieved through a specific contact between the subunits involving a residue within the mobile loop (F272 on the α-subunit and Y151 on the β-subunit). This mobile loop is the most highly conserved region of the luciferase sequence, and it was suggested that it adopts discrete “open-close” conformations occurring as a slow millisecond event (
Francisco et al., 1998;
Campbell et al., 2009,
2010) to allow for stabilization of the high-energy intermediate. The mobile loop conformational changes might be also associated with binding of the antenna protein.
Structure of lumazine protein
The bioluminescence maximum (470 nm) and spectral distribution of most species of the genus
Photobacterium is the same as the fluorescence of lumazine protein (LumP; 21 kDa) (
O’Kane et al., 1985,
1991;
Lee, 1993). LumP purified from extracts of these bacteria contains a noncovalently bound fluorophore, 6,7-dimethyl-8-ribityl-lumazine, or simply lumazine, the precursor of riboflavin in its biosynthesis process. If included in the
in vitro bacterial luciferase reaction, LumP blue-shifts the bioluminescence from its 490-nm maximum and also enhances the bioluminescence intensity (Fig. 1B) (
O’Kane et al., 1991;
Petushkov et al., 1996b). Lumazine protein (PDB code 3A3G) is a monomer with a spatial structure consisting of two barrels with nearly identical folds (Cα RMSD 0.90 Å, sequence identity 27%). Both β-barrel domains are composed of six antiparallel β-strands and two or three small α-helices, coinciding well with the structure of riboflavin synthase, for which the same lumazine derivative is a substrate (Fig. 1G) (
Liao et al., 2001;
Chatwell et al., 2008). The lumazine molecule can only bind to the N-terminal domain of LumP, in a shallow groove comprised of residues of strands β4 and β5 and of helix α2 (Fig. 1H). Binding of the lumazine in the C-terminal domain, unlike in the synthase, is prevented by steric hindrance from H145 (
Sato et al., 2010). Stabilization of the ring system and the ribityl moiety of the lumazine in the active site of LumP is achieved through a tight hydrogen bond network involving water molecules. The LumP structure is well consistent with previous biophysical studies on LumP. A polar environment of the lumazine binding site was previously suggested from the fluorescence spectral positions, fluorescence quenching, circular dichroism, and NMR data (
Kulinski et al., 1987;
Vervoort et al., 1990;
Lee et al., 1992;
Lee, 1993). Evidence suggested that the site has a high surface exposure (
Lee et al., 1985) and contains no Trp, Tyr, Cys, or His residues which otherwise would quench the high quantum yield of lumazine fluorescence. Fluorescence polarization studies showed that the lumazine derivative attaches itself by burying the ribityl group into the protein, with the lumazine electronic system exposed on the surface (
Lee et al., 1992).
Another antenna protein found in bioluminescent bacteria is yellow-fluorescent protein (YFP; 23 kDa), isolated from a unique strain Y-1, of
Vibrio fischeri. YFP contains riboflavin or flavin as a fluorophore and modulates the blue bioluminescence color of the
in vitro Y-1 luciferase reaction to the yellow (maximum 542 nm) (
Daubner et al., 1987;
Macheroux et al., 1987;
Baldwin et al., 1990;
Petushkov and Lee, 1997). LumP and YFP are homologous (37%) (
O’Kane et al., 1991) and also share homology with riboflavin synthase, which binds two molecules of the lumazine (
O’Kane and Prasher, 1992). A remarkable property of YFP is that the bound FMN exhibits a long fluorescence lifetime (7.6 ns) and high quantum yield (0.6) commensurate with its function, as most flavoproteins are hardly fluorescent at all (
Visser et al., 1997).
Lumazine protein-luciferase interaction
Addition of lumazine protein to the
in vitro reaction of luciferase from
Photobacterium shifts the bioluminescence emission to a shorter wavelength (maximum 470 nm) in a concentration dependent manner, such as even at micromolar amounts of proteins the emission spectrum corresponds to the fluorescence of the added lumazine protein (Fig. 1B). The mechanism by which the excitation energy of the luciferase intermediate can directly populate the fluorescence state of the lumazine molecule in LumP is attributed to Förster resonance energy transfer (FRET) (
Wu and Brand, 1994). For FRET to take place there must be a significant spectral overlap between the fluorescence spectrum of the donor and the absorption spectrum of the acceptor, as well as a donor-acceptor separation less than about 100 Å. At micromolar concentrations of luciferase and LumP (where the average distance between unassociated molecules exceeds 1000 Å), the energy transfer has to involve the formation of a LumP-luciferase complex to bring the lumazine into proximity with the luciferase reaction site (
Lee et al., 1991;
Lee, 1993). The existence of the YFP-luciferase complex is evidenced in an analogous manner (
Petushkov et al., 1996a,
1996b).
The LumP-luciferase association is directly shown from dynamic fluorescence experiments. The fluorescence of LumP (maximum 475 nm upon excitation at 375 nm) is due to the strongly bound lumazine molecule, and has a correlation time of anisotropy decay φ = 20 ns (4°C), consistent with the 21-kDa mass of LumP. Addition of Photobacterium luciferase increases the correlation time to φ = 100 ns indicating a higher mass of the rotator (75 + 21 kDa), which in the present system corresponds to the LumP-luciferase complex (Fig. 1K). Analysis of anisotropy decay titrations suggested a K
D in the 0.2–1 mmol/L range, an interaction much weaker than expected from the only micromolar concentrations required to affect the bioluminescence spectrum (
Visser and Lee, 1982). It was proposed that a tighter interaction might be found with a luciferase reaction product such as the so-called “Fluorescent Transient,” formed following the reaction of the luciferase-bound FMNH
2 with O
2 and aldehyde (
Matheson and Lee, 1983). This fluorescent transient shows an expected φ = 82 ns (4°C) but it is reduced to 4 ns, not increased to 100 ns, on addition of LumP (Fig. 1L). This effect is attributed to energy transfer as a main channel of anisotropy loss concomitant with the bioluminescence spectrum matching that of the LumP fluorescence. Again, analysis of a titration based on this fluorescence anisotropy loss reveals a K
D < 10 µmol/L, consistent with the micromolar concentrations required for the bioluminescence spectral shift (
Lee et al., 1989,
1991;
Lee, 1993). For the YFP-luciferase interaction, the rate of energy transfer was calculated as k
ET = 10 ns
−1 (ten times faster than with LumP k
ET ≈ 1 ns
−1), and K
D = 0.7 µmol/L, evidence that YFP (and LumP in a less strong manner) has a catalytic effect on the luminescence reaction (
Eckstein et al., 1990;
Petushkov et al., 1995,
1996a,
1996b).
The spatial structure of a possible LumP-luciferase complex was derived from computational modeling using structures of bacterial luciferase with FMN (PDB code 3FGC), and LumP (PDB code 3A3G) as starting structures (
Sato et al., 2010). The structure of the complex is shown in Fig. 1M. In this model structure, LumP contacts the reentrant surface of luciferase near the FMN site. The ring system of the lumazine is located close enough to the isoalloxazine ring of the FMN on the luciferase (approximately 10 Å) to allow their electronic transitions to be coupled and hence to account for moderately efficient energy transfer from the fluorescent transient to the bound lumazine. The electrostatic potentials of the molecular surfaces show that the lumazine binding site of LumP is covered by negatively charged residues (D19, D20, D62, and D64) (Fig. 1H), while the FMN site of luciferase forms a concave shape, which is surrounded by positively charged residues (R107, R125, K259, K290, and K291) (Fig. 1F). This countercharge distribution is characteristic of a variety of weak transient protein-protein complexes (
Prudêncio and Ubbink, 2004;
Kiel et al., 2004;
Reichmann et al., 2007) and is believed to provide long-range electrostatic forces for the protein-protein interaction, which impacts on the association rate constant k
A for complex formation. It is highly possible that long-range forces pre-orient LumP and luciferase molecules within an “encounter” complex, and that the successive formation of the fluorescent transient would “glue” the complex by decreasing its dissociation rate constant k
D (
Petushkov et al., 1996b;
Schreiber et al., 2006). The fluorescent transient’s breakdown would lead to dissociation of the complex allowing for the next turnover of luciferase with FMNH
2, O
2 and aldehyde. It also has to be emphasized, that a luciferase complex must be formed prior to the actual bioluminescence excitation step so that the LumP is available to trap the energy before it is emitted at the longer wavelength (energy transfer rate k
ET ≈ 1 ns
−1) (
Visser and Lee, 1982;
Lee et al., 1991).
JELLYFISH CLYTIA
C. gregaria is a pelagic hydrozoan jellyfish commonly found in the eastern Pacific Ocean (Fig. 2C). Its green bioluminescence (maximum 509 nm), thought to have a defensive function, is emitted upon mechanical stimulation from light organs (photophores) distributed evenly along the edge of the umbrella, and containing the photoprotein clytin and its cognate green-fluorescent protein (Clytia GFP).
Clytia utilizes coelenterazine, the most widely found luciferin among marine animals, as a substrate for its bioluminescence reaction. Like the photoprotein aequorin from the jellyfish
Aequorea victoria, clytin is a Ca
2+-regulating photoprotein, and tightly binds and stabilizes a hydroperoxy-substituted coelenterazine within an inner cavity. Binding of Ca
2+ to these photoproteins triggers the bioluminescence reaction to produce the excited state product coelenteramide and CO
2 (Fig. 2A). In the GFP-clytin complex, this excitation step involves again a Förster-type transition coupling with the green fluorescence of the Clytia GFP’s fluorophore (maximum 500 nm). With no Clytia GFP presence, the bioluminescence emission has a maximum around 470 nm (Fig. 2B) (
Levine and Ward, 1982;
Inouye and Tsuji, 1993;
Markova et al., 2010). That these photoproteins are activated or regulated by Ca
2+ implies that the bioluminescence system is under control of the neural system of the animal (
Dunlap et al., 1987).
Structure of the photoprotein clytin
Clytin (22.4 kDa, PDB code 3KPX) shares high structural and sequence similarity with the other Ca
2+-regulated photoproteins, obelin from
Obelia (
Liu et al., 2000) and aequorin from
Aequorea (
Head et al., 2000). Photoproteins can be regarded as a luciferase except with its high-energy reaction intermediate, 2-hydroperoxycoelenterazine, frozen in place (
Hastings and Gibson, 1963). The spatial structure of the protein contains four helix-turn-helix motives known as EF-hands, wrapped to form a two-bowls like globular structure. This structure protects and stabilizes the 2-hydroperoxycoelenterazine, and allows binding of Ca
2+ to trigger the bioluminescence reaction. EF-hands I and II are in the N-terminal, and III and IV are in the C-terminal (Fig. 2E). Only three EF-hands I, III, and IV bind calcium, while EF-hand II lacks the canonical sequence in the loop to permit calcium binding (
Deng et al., 2005;
Vysotski and Lee, 2007). The hydrophobic substrate-binding cavity is comprised of highly conserved residues from almost every α-helix (A–H). Tryptophans W95, W182, W117, and W138 form a “sandwich” structure to orient the 2-hydroperoxycoelenterazine. Residues Y193 and H178 form a hydrogen bond network holding the highly unstable peroxy group of 2-hydroperoxycoelenterazine. The mechanism of the bioluminescence reaction involves conformational changes of clytin upon binding of Ca
2+ in the EF-hand loops. This disturbs the hydrogen bond network around the peroxy group, which in turn triggers decarboxylation of the 2-hydroperoxycoelenterazine, ending up in formation of the product excited coelenteramide and emission of blue light (Fig. 2D) (
Head et al., 2000;
Liu et al., 2006). Clytin completes only one turnover
in vitro because after the bioluminescence reaction the product remains tightly bound.
Structure of Clytia GFP
Clytia GFP from
C. gregaria (PDB code 2HPW) forms the classical β-barrel structure (30 × 40 Å) characteristic of this family of fluorescent proteins. Despite a low sequence identity, the structure of Clytia GFP highly resembles GFPs from
Aequorea (
Ormö et al., 1996) and
Renilla (
Loening et al., 2007). Eleven parallel β-strands (S1–S11) comprise the “walls” of the barrel, and loops cover the barrel from the “top” and “bottom” blocking solvent access to the internal cavity (Fig. 2G). The chromophore of Clytia GFP, 5-(
p-hydroxybenzylidene)imidazolin-4-one, is a product of post-translational autocatalytic cyclization of S68, Y69, and G70 amino acid residues (
Remington, 2006;
Wachter, 2007). Covalently attached to the core α-helix, the fluorophore resides in the center of the protein cylinder in a polar environment, very well protected from the solvent. This apparently minimizes energy dissipation upon fluorescence excitation to favor a high quantum yield of the Clytia GFP fluorescence (0.86; maximum 500 nm upon excitation at 470 nm) (
Markova et al., 2010). From analytical ultracentrifugation, size-exclusion chromatography and fluorescence anisotropy data, Clytia GFP is an obligate dimer (
Titushin et al., 2010;
Malikova et al., 2011). The dimerization interface of the Clytia GFP monomer reveals hydrogen bonding and an accessible surface area value (1370 Å
2) favorable for strong dimerization (
Jones and Thornton, 1996). Within the dimer, the two fluorophores have a coupled electronic state transition (
Malikova et al., 2011).
Clytia GFP-clytin interaction
Isolated clytin converts 2-hydroperoxycoelenterazine into coelenteramide with emission of a broad blue bioluminescence (maximum 470 nm). Upon titration with Clytia GFP, a second longer wavelength band appears and at only micromolar concentrations of both proteins, the bioluminescence (maximum 497 nm) becomes identical to that of the Clytia GFP fluorescence (Fig. 2B). As in the case of bacterial luciferase and lumazine protein, the energy transfer in the Clytia GFP-clytin system is rationalized by the FRET mechanism. The requirement of donor-acceptor separation to be less than 100 Å at micromolar concentrations of clytin and Clytia GFP can only be achieved assuming a tight protein-protein interaction (
Förster, 1960;
Wu and Brand, 1994;
Titushin et al., 2010;
Markova et al., 2010). Although the bioluminescence spectrum shift is complete at micromolar amounts of proteins, implying K
D of Clytia GFP-clytin complexation to be in the micromolar range, no evidence for protein-protein association could be found from analytical ultracentrifugation, size-exclusion chromatography, surface-plasmon resonance, or polarization fluorometry experiments, which are capable of detecting complexation with a K
D below 10
-4 mol/L (
Titushin et al., 2010;
Malikova et al., 2011).
However, a weak complex could be detected with NMR chemical shift perturbation mapping, a technique that is capable of measuring weak protein-protein interactions up to K
D = 10
−2 mol/L, owing to the remarkable sensitivity of
15N-
1H-HSQC chemical shifts to subtle conformational effects (
Zuiderweg, 2002;
Vaynberg and Qin, 2006). When an
15N-labeled clytin sample is titrated with unlabeled Clytia GFP, about 15 peaks on the
15N-
1H HSQC spectrum of clytin experience concentration dependent chemical shift perturbation assignable to interaction with Clytia GFP (Fig. 2I and 2K). Analogously, the
15N-
1H HSQC spectrum of
15N-labeled Clytia GFP reveals approximately 20 peaks with chemical shift altered upon titration with unlabeled clytin (Fig. 2J and 2L). The perturbed shifts are small in amount and also exist in the fast chemical exchange regime on the NMR-time scale; therefore no significant conformational changes of the proteins occur upon the interaction. Subsequent to assignment of the HSQC spectra, the perturbed peaks could be attributed to certain amino acid residues of both clytin and Clytia GFP, and accordingly, located on their 3-D spatial structures as belonging to the interaction surfaces of the proteins.
NMR chemical shift perturbation mapping indicated a continuous path on the clytin surface involving the amino-terminus, α-helix D, and the carboxy-terminus, which are contiguous in the spatial structure (Fig. 2F). Residues of Clytia GFP under “perturbation” clustered on loops S6-S7 and S10-S11, comprising the “top” of the Clytia GFP barrel (Fig. 2H). The crystal structures of clytin and Clytia GFP along with the above designated interaction surfaces were the basis for computing the spatial structure of the Clytia GFP-clytin complex by means of modeling using the “soft-docking” program HADDOCK2.0 (
Dominguez et al., 2003;
de Vries et al., 2007). The resulting lowest energy structure is shown in Fig. 3A. This Clytia GFP-clytin structure bears the characteristics of a weak protein-protein complex, with a relatively small number of intermolecular contacts, an average interaction surface area (
Lo Conte et al., 1999), and a notable charge complementarity of the interfaces. Basic amino acid residues (K11, K13, K100, and K104) introduce positive charge into the interface of clytin, and complementarity to it (also in shape) lies in the acidic loop S10-S11 of Clytia GFP, comprised of residues D211, D213, D214, D215, and E216 (Fig. 3B and 3C). Apart from electrostatic forces, the complexation is governed by a number of hydrophobic contacts and a hydrogen bond network. The fluorophores of clytin and Clytia GFP lie 45 Å apart which, together with the very favorable spectral overlap of clytin bioluminescence and Clytia GFP absorbance (J = 1.3 × 10
−13 L·mol
−1·cm
−3), would account for efficient energy transfer from the excited state coelenteramide to the fluorophore of Clytia GFP. Each monomer of the Clytia GFP homodimer can accommodate one molecule of clytin, which would argue for a tetrameric complex of 2 clytins per Clytia GFP (dimer). The supposition of a heterotetrameric complex was advanced for
in vivo aequorin-Aequorea GFP (
Cutler and Ward, 1997).
Both clytin and Ca
2+-loaded clytin (that is “before” and “after” the bioluminescence reaction) interact only weakly with Clytia GFP as observed from NMR titration, polarization and dynamic fluorescence, and isothermal titration calorimetry experiments (K
D, ~0.9 mmol/L), whereas the bioluminescence effect of Clytia GFP is complete at micromolar concentrations of proteins. To rationalize this contradiction it was hypothesized that on triggering by binding Ca
2+, clytin adopts an intermediate conformation state regarded as a “pre-excited state” clytin, in which its affinity to Clytia GFP is enhanced several orders of magnitude. Within this short lived complex, the chemiexcitation step occurs to populate the coelenteramide-GFP coupled electronic state, resulting in an S
1→S
0 radiative transition corresponding to the fluorescence of the acceptor, Clytia GFP. This proposal bears similarity to the bacterial lumazine protein-luciferase interaction discussed above, and the property of the energized intermediate of the luciferase reaction in enhancing the lumazine protein-luciferase interaction (
Petushkov et al., 1996a). Hence, the structure of the Clytia GFP-clytin complex derived from NMR mapping very probably only represents the structure of the encounter complex of low affinity, in which clytin (of the “before” luminescence reaction state) and Clytia GFP are preoriented, and then this is followed by formation of the pre-excited state complex. Ionic strength sensitivity of the Clytia GFP-clytin interaction confirms the significant role of electrostatic forces known to be an attribute of an encounter complex (
Berg and von Hippel, 1985;
Sheinerman et al., 2000;
Schreiber et al., 2006). Supporting the validity of the Clytia GFP-clytin structure were mutational experiments indicating that substitutions of the residues on the clytin interface decrease both the complex affinity and the energy transfer efficiency, with substitutions disturbing polar contacts having the largest effect.
SOFT CORAL RENILLA
The sea pansy
Renilla found in shallow sandy littoral bottoms emits bright green bioluminescence (maximum 505 nm) probably as a defensive reaction to mechanical disturbance, putatively to deter predators from nibbling at the polyps (Fig. 4C) (
Cormier et al., 1974,
1975). The bioluminescence originates from single polyps, wherein the system is packaged in specialized light emitting cells (photocytes) that contain luciferase, Renilla GFP, and the substrate coelenterazine bound to the Ca
2+-triggered coelenterazine-binding protein (CBP) (Fig. 4A) (
Anderson and Cormier, 1973;
Cormier and Charbonneau, 1977;
Cormier, 1978;
Charbonneau and Cormier, 1979). When the Ca
2+ concentration rises upon nerve stimulation, Ca
2+ binds to CBP to make coelenterazine available for reaction with luciferase and O
2. It is proposed that the bioluminescence reaction occurs within a CBP-luciferase complex so that coelenterazine does not necessarily dissociate (
Charbonneau and Cormier, 1979;
Titushin et al., 2008;
Stepanyuk et al., 2009). Following oxygen addition to coelenterazine in the active site of luciferase, the decarboxylation reaction deposits excitation again into the coupled electronic state transitions of the coelenteramide-GFP fluorophore pair, from which S
1→S
0 radiation as GFP’s green fluorescence takes place (509 nm). In the absence of GFP the light is blue (480 nm) (Fig. 4B) (
DeLuca et al., 1971;
Hart et al., 1979;
Shimomura, 2006). To enable the energy transfer, Renilla luciferase and Renilla GFP are required to form a protein-protein complex (
Ward and Cormier, 1976,
1978;
Ward, 1979), and so luciferase must interact simultaneously with both proteins, CBP and GFP. Despite solid biophysical data on Renilla GFP-Renilla luciferase complexation, no structural information is available for this interaction at the present time.
Structure of Renilla luciferase
Renilla luciferase (37 kDa, PDB code 2PSJ) is a monomeric protein, which catalyzes the decarboxylation of coelenterazine in the presence of molecular oxygen, resulting in the product coelenteramide, carbon dioxide, and blue light (maximum 480 nm) (
Karkhanis and Cormier, 1971;
Matthews et al., 1977a,
1977b;
Lorenz et al., 1991). The crystal structure was solved for the stabilized mutant form of luciferase from
R. reniformis and is shown in Fig. 4F. It has a characteristic α/β-hydrolase fold sequence at its core, that is an α/β-sheet of eight β-sheets connected by α-helices (
Ollis et al., 1992), and shares the conserved catalytic triad of residues D120/E144/H285A employed by the dehalogenases, with D120 being a part of a characteristic “nucleophile elbow.” (
Holmquist, 2000;
Loening et al., 2006) However, in Renilla luciferase this triad is employed in an oxidation reaction, while in dehalogenase it catalyzes a hydrolysis reaction, which makes a high level of primary sequence and tertiary structure similarity between them very intriguing (
Loening et al., 2007). This triad and residues N53/W121/P220 comprise the catalytic core clustered toward the bottom of the bowl-shaped active site. Enzymatic assay, mutagenesis data, and docking computation reveal that catalytic triad D120/E144/H285A is directly involved in the decarboxylation reaction of coelenterazine, apparently coordinating the attack of the coelenterazine molecule by molecular oxygen (
Woo et al., 2008). Around this core is a ring of aromatic residues, which assure specificity of binding and right orientation of the hydrophobic coelenterazine. The bioluminescence reaction of Renilla luciferase is dependent upon the flexible cap domain, which is thought to assist binding of coelenterazine, and may also be responsible for an immediate relocation of the product coelenteramide in the active site after the enzymatic reaction (
Loening et al., 2007).
Structure of coelenterazine-binding protein (CBP)
Coelenterazine-binding protein is a Ca
2+-regulated protein (21 kDa, PDB code 2HPS), containing coelenterazine bound within its inner cavity (
Charbonneau and Cormier, 1979;
Kumar et al., 1990;
Inouye, 2007;
Stepanyuk et al., 2009). Despite low sequence identity, the structure of CBP closely resembles structures of photoproteins obelin, aequorin and clytin. Four sets of helix-turn-helix structural motives, which are characteristic of EF-hand calcium-binding domains, comprise again two “bowls” folded over the coelenterazine molecule (Fig. 4H). Unlike in photoproteins, coelenterazine in CBP is found unsubstituted, without the peroxy group. This is consistent with their difference in function, i.e. CBP does not catalyze oxidation of coelenterazine, but rather undergoes conformational changes triggered by binding of calcium, resulting in release of coelenterazine for oxidation and decarboxylation in the active center of luciferase (
Stepanyuk et al., 2009). In this regard, a photoprotein can be thought of having the combined functions as a substrate binding protein and as a luciferase, with a Ca
2+-sensitive link to a neural network. The CBP holding coelenterazine also protects it from autooxidation and rapid diffusion otherwise allowed by the hydrophobic nature of the molecule (
Shimomura, 1997). Compared to photoproteins, the substrate-binding cavity of CBP is more hydrophilic, contains no tryptophan residues, and coordinates coelenterazine by hydrogen bonding with residues Y36, R22, R19, K139, D183, and five water molecules, which altogether account for weaker affinity of coelenterazine to the binding cavity. Availability of spatial structures of CBP and apo-CBP loaded with Ca
2+ reveals that binding of Ca
2+ to three EF-hand loops of CBP stretches their conformation, so that helices C, D, and E are moved apart to close a small aperture on the protein surface, and a new and bigger one opens on the opposite side (Fig. 4I). Displacement of the cavity comprising residues “pushes out” the coelenterazine to the newly opened aperture, which would look like a perfect gate for coelenterazine to escape, for it is no longer tightly bound in the cavity due to the loss of the original contacts after shift of the hydrophobic side chains (Phe, Met, Val, and Leu), and because the R19 “lid” becomes opened up (
Stepanyuk et al., 2009).
CBP-Renilla luciferase interaction
Indirect evidence for this CBP-luciferase association is derived from kinetic characteristics of the
in vitro luciferase reaction. The bioluminescence reaction of Renilla luciferase is triggered either by injection of coelenterazine into a solution of luciferase, or by injection of calcium into a mixture of luciferase and CBP (
Matthews et al., 1977a;
Titushin et al., 2008). Both reactions have very similar bioluminescence spectra (Fig. 4B), implying that the chemical mechanism of the reaction and structure of the bioluminescence emitter must be identical. However, kinetic parameters of the reactions are different, such as luciferase is able to turn over with coelenterazine enclosed within CBP six times faster than with free coelenterazine, and emits with a two-fold higher quantum yield (Fig. 4D and 4E) (
Titushin et al., 2008). Knowing that apo-CBP does not enhance luciferase activity, a higher bioluminescence quantum yield for Renilla luciferase with CBP as its substrate, would be difficult to explain if coelenterazine had to first dissociate from Ca
2+-loaded CBP into solution. The absorbance spectra of CBP, Ca
2+-loaded CBP, and an aqueous solution of coelenterazine, differ in the position of the visible maximum (444 nm, 434 nm, and 415 nm respectively), revealing that for these three cases, the chemical environment of the coelenterazine molecule must be different. Besides, Ca
2+-loaded CBP is weakly fluorescent, while an aqueous solution of coelenterazine is non-fluorescent (
Shimomura and Teranishi, 2000), which altogether implies that coelenterazine does not necessarily dissociate, but rather stays bound, at least for several minutes, within the cavity of Ca
2+-loaded CBP. Then for the bioluminescence reaction to proceed, coelenterazine should be delivered to the active center of luciferase within a transient protein-protein complex, which dissociation should be rapid enough to allow the luciferase turnover to be 2–5 s
−1 (
Hart et al., 1979;
Loening et al., 2006). The characteristic for such a transient complex is usually an electrostatic complementarity of interaction surfaces of proteins (
Prudêncio and Ubbink, 2004;
Kiel et al., 2004;
Reichmann et al., 2007). This complementarity is observed between the path around the active site of luciferase and the path around the aperture of Ca
2+-loaded CBP (Fig. 4G and 4I). Apart from this complementarity, the 2-(
p-hydroxybenzyl) group of coelenterazine, the region of oxygen attack catalyzed by luciferase, also becomes exposed through the newly opened aperture of the Ca
2+-loaded CBP.
To obtain the spatial structure of the Ca
2+-loaded CBP-luciferase complex, a docking simulation was performed using structures of the luciferase from
R. muelleri modeled in SWISS-MODEL (
Arnold et al., 2006) on the highly homologous structure of the luciferase from
R. reniformis (PDB code 2PHW), and the structure of Ca
2+-loaded CBP of
R. muelleri (PDB code 2HQ8). Both an interaction-restraints driving docking (HADDOCK2.0) and an
ab initio docking software (PatchDoc) (
Schneidman-Duhovny et al., 2005) resulted in lowest energy structures with very similar topology (
Stepanyuk et al., 2009;
Titushin, 2009). On these structures, as well as on the best energy score structure shown in Fig. 4J, the opening on the surface of the Ca
2+-loaded CBP structure interacts with good shape complementarity, with the path flanking the active site of luciferase. On the interface of Ca
2+-loaded CBP, residues R9, K91, K93, K97, and D66, E80, D87 interact complementarily in charge, with corresponding residues D158, E160, E161, and K4, K189, K193, K282 of luciferase (Fig. 4G and 4I). The number of hydrogen bonds and intermolecular contacts imply a weak interaction with a moderate size interaction surface area (~1800 Å
2). There is a path observed in the structure, which might serve as a route for coelenterazine to relocate from the substrate-binding cavity of Ca
2+-loaded CBP towards the active site of luciferase to undergo the catalytic oxidation and decarboxylation (Fig. 4K). The resulting product of the reaction, excited state coelenteramide, may be regarded as partly protected from the solution by the protein-protein complex environment, and is thus subjected to a lesser energy dissipation, which would explain the higher quantum yield of the reaction of luciferase with CBP rather than with free coelenterazine. It should be noted that the spatial structure of the complex needs more verification, and more attempts should be tried to apply NMR-mapping and X-ray crystallography if we are to have a more substantial understanding of the mechanism of Renilla bioluminescence.
CONCLUDING REMARKS
It has been suggested that coelenterazine may have a role in the metabolism of marine creatures other than as a substrate for bioluminescence (
Shimomura, 2006;
Haddock et al., 2010). If this is the case, then it would be another example like 6,7-dimethyl-8-ribityllumazine (lumazine), of a metabolite being seconded in the process of evolution, to perform a novel function. The chemistry of the light reaction limits the natural bioluminescence spectra to maxima in the range 465–495 nm, but the presence of the GFPs extends this to slightly longer wavelengths, apparently to better suit the photic environment closer to surface regions. In the bioluminescent bacteria, the antenna proteins LumP and YFP also alter the color of bioluminescence over about the same range, but additionally, participate in the bioluminescence reaction within the protein-protein intermediate, influencing the kinetics of the reaction. Renilla coelenterazine-binding protein stores coelenterazine protecting it from oxidation and, being a calcium-regulated protein, links the luminescence response to a neural signal. Delivery of coelenterazine to luciferase within the putative CBP-luciferase complex remarkably enhances the bioluminescence quantum yield.
The above mentioned bioluminescence protein-protein complexes are characterized by being weak and transient with affinity constants K
D in the millimolar range and higher. This provides its own set of impediments on how to detect and structurally characterize the interactions. So far none of the bioluminescence complexes have provided crystals with good diffraction, and in the best case, we are left with knowledge of the interaction surfaces derived from biophysical experiments including NMR, which could be mapped onto the crystal structures of individual proteins. Interaction restraints are input to a docking program to derive a number of lowest free energy structures representing the overall topology of the protein-protein complex. Although the highest limit of molecular size in NMR studies is constantly increasing, specifically defining interaction restraints with NOE NMR in the complexes with weak protein-protein interactions is still challenging, and therefore other NMR techniques such as residual dipolar couplings (RDC) and chemical shift perturbations (CSPs)-based docking, are quite important in the studies of these weak transient complexes (
Wang et al., 2011).
Despite that the spatial structures available for several bioluminescence protein-protein complexes need more verification, some of their properties are obviously common. They all share the property of charge complementarity of the interaction surfaces, a relatively low number of hydrogen bonds and intermolecular contacts, and with interaction surface areas having values about average among the known protein-protein complexes. Another interesting similarity is that the interactions with antenna proteins are also governed by some intermediate events. In the bacterial system, a metastable luciferase-bound intermediate is required to increase the affinity to lumazine protein over the unreacted luciferase, and its breakdown promotes the complex dissociation. The photoprotein clytin from Clytia is triggered for conformational changes upon binding of Ca2+, followed by formation of a “high energy intermediate” which, it is proposed, might have a shifted protein conformation now a thousand times more appealing to Clytia GFP than clytin “before” or “after” the bioluminescence.
Bioluminescence protein-protein interactions are not confined to the ones discussed from the structural perspective in this review. Bacterial NADPH:FMN oxidoreductase delivers newly synthesized FMNH
2 to luciferase within a protein-protein complex, complicating the organization of the protein-protein interaction network involved in the bioluminescence system of bacteria. In
Renilla, the energy transfer mediating interaction is known to occur between luciferase and GFP and, even though the complex was reported to be cocrystallized (
Loening et al., 2007), the structural basis for this interaction still remains unknown. It is of particular interest to know how the three proteins of the
Renilla system, luciferase, CBP, and GFP, which are found in photocytes in approximately equal molar concentrations (
Ward and Cormier, 1978), are organized and operate within a multi-protein complex.
Higher Education Press and Springer-Verlag Berlin Heidelberg 2011