Abstract:
The present invention relates generally to hydrogen production for use in fuel cells, foodstuffs and chemical production, and more particularly, to biologically and photosynthetically produced hydrogen. Specifically, disclosed is a method for producing bacteria and green alga that can produce hydrogen in quantities that exceed four hundred percent of the hydrogen produced by green alga in nature; thus, producing organisms which can serve as hydrogen generators for fuel cells, chemical production and numerous other applications.

Description:
CROSS REFERENCE TO RELATED APPLICATIONS 
     This application is a continuation-in-part of U.S. Pat. No. 8,124,347 issued on Feb. 28, 2012 and filed as U.S. patent application Ser. No. 12/428,471 on Apr. 23, 2009, which claims the benefit of U.S. Provisional Patent Application No. 61/046,989 filed on Apr. 22, 2008, entitled “Photosynthetic Hydrogen Production from the Green Alga  Chlamydomonas Reinhardii ”, the entire disclosures of which are hereby incorporated by reference for all purposes. 
    
    
     INCORPORATION OF SEQUENCE LISTING 
     The entire contents of a paper copy of the “Sequence Listing” and a computer readable form of the sequence listing on diskette, containing the file named Seq_Listing_ST25.txt, which is 489 kilobytes in size and was created on Apr. 23, 2009 are herein incorporated by reference. 
     TECHNICAL FIELD 
     The present invention relates generally to hydrogen production for use in fuel cells, foodstuffs and chemical production, and more particularly to biologically and photosynthetically produced hydrogen. 
     BACKGROUND OF THE INVENTION 
     Evidence has shown that the combustion of fossil fuels is causing a change in the composition of our atmosphere. The resulting increase in average global temperature requires an immediate and global response. A recent British climate change report suggests that we would have to decrease emissions of carbon dioxide and other greenhouse gases by 25% by the year 2050 to avoid as much as a 20% decrease in global Gross Domestic Product (GDP) caused by catastrophic drought, flooding, and disease. Ultimately, an 80% decrease in emissions would be necessary. So, if not fossil fuels, what should we use as a source of energy? If we switched to a hydrogen economy and utilized molecular hydrogen and fuel cells in all of our cars, trucks, trains, etc., a 50% reduction in the emission of carbon monoxide and nitrous oxides is likely. Of course, a decrease in emissions and a concomitant improvement in climate change is dependant on how the hydrogen is produced. Climate change would not occur if we continue to produce hydrogen by the steam reformation of natural gas and coal as this process results in localized emissions, but emissions nonetheless. However, if the hydrogen were produced biologically from water, perhaps by a photosynthetic organism, there would be little or no release of carbon dioxide, nitrous oxides, or methane. 
     Hydrogen is currently produced by steam reforming the hydrogen atoms from coal or natural gas. The reactions are: CH 4 +H 2 O→CO+3H 2  (natural gas) or C+H 2 O→CO+H 2  (coal) and CO+H 2 O→CO 2 +H 2 . Either fuel could be the basis of a national hydrogen economy; however both fuels generate carbon dioxide, which would add greenhouse gases to our atmosphere. If future coal driven hydrogen power plants utilized carbon sequestration, pumping the carbon dioxide into a deep underground location, this problem could be minimized. Alternatively, a carbon neutral hydrogen economy could be realized if hydrogen could be produced from the electrolysis of water where the electricity, the impetus for the reaction, is generated from a nuclear reactor, wind energy, or solar power or through photosynthetic hydrogen generation. 
     The study of biological hydrogen production in green algae began as a curiosity and after 75 years of research, its evolutionary origin still remains an enigma. General progress in the field has been ongoing since Hans Gaffron early 1940s discovery that the green alga  Scenedesmus obliquus  produced hydrogen; however, the last decade is marked by dramatic advances. Specifically, the hydrogenase genes for several species of green algae have been sequenced and the crystal structure determined, for two homologous bacterial hydrogenases,  C. pasteurinum  and  D. desulfuricans . In addition, the mechanism by which a hydrogenase creates molecular hydrogen has been elucidated from extensive research on the structure, assembly, and biological properties of all hydrogenases. 
     Hydrogenases are iron-sulfur proteins, which have played an important role in the energy metabolism of bacteria since the earliest life on Earth. In fact, homologous non-hydrogen producing iron-sulfur proteins are common in most living cells, including humans and pathogenic bacteria. The hydrogenases, however, are different from their evolutionary cousins in that their iron sulfur clusters contain unique cyanide and carbon monoxide ligands ( FIG. 1 ). 
     Hydrogen is produced by enzymatically combining protons with electrons from the photosynthetic electron transport chain. The protons and the electrons are generated from the first step in the photosynthetic cycle, the splitting of water into oxygen and protons. The electrons are immediately energized by a photon (λ=680 nm) in Photosystem II and passed from one compound to another, all of which compose the electron transport chain ( FIG. 2 ). Most of the electron carriers are quinones (Q), plastiquinones (PQ), or cytochromes (Cyt). A second input of light energy (λ=700 nm) occurs during Photosystem I and the energized electrons are passed to the terminal electron carrier, ferredoxin. At this point, the electrons can participate in CO 2  fixation, i.e. cell growth, or be transferred to the hydrogenase to produce hydrogen. 
     SUMMARY OF THE INVENTION 
     In one aspect, the present disclosure provides a method for selecting a chimeric hydrogenase which produces hydrogen at an increased rate relative to the hydrogen production rate of a wild type hydrogenase following transformation of a cell with the chimeric hydrogenase and expression of the chimeric hydrogenase in the transformed cell, the method comprising: calculating a positive to negative electrostatic potential surface area (EPSA) ratio for the chimeric hydrogenase, wherein a positive to negative EPSA ratio of about 1 to about 115 is indicative of an increased hydrogen production rate in the cell by the chimeric hydrogenase relative to the hydrogen production rate of the wild type hydrogenase. In the method, the cell may be an algal cell and the chimeric hydrogenase may be a chimeric algal hydrogenase. The cell may be for example a  Chlamydomonas  cell. Alternatively, the cell may be a bacterial cell and the chimeric hydrogenase may be a chimeric bacterial hydrogenase. In another aspect, the present disclosure also provides a chimeric hydrogenase selected according to the preceding method. 
     In another aspect, the present disclosure provides a method of producing an algae capable of enhanced hydrogen production, the method comprising: making and testing one or more mutated algal hydrogenases for amount of hydrogen production; identifying one or more of the mutated algal hydrogenases as capable of causing enhanced hydrogen production in algae, wherein an amount of hydrogen production higher than that of a wild type hydrogenase is indicative of the ability to cause enhanced hydrogen production in algae; expressing in the algae cells a DNA sequence coding for at least one mutated algal hydrogenase identified as capable of increased hydrogen production; selecting an algae expressing the mutated algal hydrogenase identified as capable of increased hydrogen production, wherein the mutated algal hydrogenase has a positive to negative EPSA ratio in the range of about 1 to about 115, which is indicative that the algae is capable of enhanced hydrogen production. In the method, making one or more mutated algal hydrogenases may comprise: identifying two or more hydrogenase parent proteins; obtaining cDNAs coding for said parent hydrogenases; producing multiple DNA fragments corresponding to segments of each of said cDNA coding for said parent hydrogenases; and reconstructing full-length chimeric hydrogenase cDNAs by putting segments from cDNAs coding for different parent hydrogenases together in the same order as they occur in the parent hydrogenase coding sequences. Full-length chimeric hydrogenase cDNAs encoding mutated algal hydrogenases can be tested for example in a bacterial system or an algal system to identify mutations capable of enhanced hydrogen production. Full-length chimeric hydrogenase cDNAs encoding mutated algal hydrogenases can be analyzed to identify specific mutations that lead to enhanced hydrogen production. The method may further comprise the step of expressing one or more mutated algal hydrogenases by transforming algae using a plasmid comprising the full-length chimeric hydrogenase cDNAs. The plasmid may be for example pSMP. In another aspect, the present disclosure also provides a recombinat algal cell produced according to the preceding method. The recombinant algal cell may be a  Chlamydomonas  cell, such as but not limited to a  C. reinhardtii  cell. 
     In any of the methods, calculating the positive to negative EPSA ratio may comprise computing the positive EPSA, the negative EPSA, and the ratio of the positive EPSA to the negative EPSA. The positive to negative EPSA ratio can be for example from about 2 to about 50, or from about 5 to about 20. 
    
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         FIG. 1  is a diagram illustrating the four-iron four-sulfur (4Fe-4S) cluster and the two-iron two-sulfur (2Fe-2S) active site cluster that is present in Fe-only hydrogenases. 
         FIG. 2  illustrates the photosynthetic electron transport chain, known as the “Z-scheme.” 
         FIG. 3  shows the crystal structure of the  C. pasteurinum  bacterial hydrogenase. 
         FIG. 4  is an example of sequence space; the set of all possible amino acid sequences. According to the disclosure, the protein of interest is just 2 amino acids in length. 
         FIG. 5  is a diagram of a “Family Shuffle.” 
         FIG. 6  is the RACHITT technique for generating a combinatorial library. 
         FIG. 7  is the sequence comparison of the two  C. reinhardtii  hydrogenase genes. 
         FIG. 8  is a plasmid map of the algal expression vectors: A) pSMP1, B) pSMP1c, and C) pSMP2. 
         FIG. 9  is a picture of an agarose gel (1%) showing the results of PCR on the genomic DNA extracts of algal transformants. 
         FIG. 10  shows the RNA expression levels by real time PCR of cDNA from algal transformants. 
         FIG. 11  shows the western blot of one of the transformants and a non-transformed control. 
         FIG. 12  shows the primers (SEQ ID NOS: 3-12) utilized to create restriction sites in the pSMP1, pSMP1c, and pSMP2 plasmids. 
         FIG. 13  illustrates the degenerative oligonucleotide gene shuffling (DOGS), the technique used to create the chimeric library. 
         FIG. 14  illustrates the design of the primers for the original (DOGS) technique. 
         FIG. 15  illustrates the variation of the degenerative oligonucleotide gene shuffling (DOGS) technique that utilizes the exonuclease SapI. 
         FIG. 16  is shows the pET DLS plasmid (8.2 kb) containing the  C. acetobutylicum  HydE accessory protein and the  C. thermocellum  hydrogenase (pET DLS C. thermo plasmid). 
         FIG. 17  illustrates the chimeric hydrogenases that were tested for hydrogen production. 
         FIG. 18  shows the hydrogen production of the specified chimeric hydrogenases as a percentage of the positive control ( C. acetobutylicum ). 
         FIG. 19  shows the primers (SEQ ID NOS: 17-22) utilized to clone the Clostridial hydrogenases into the pET DLS expression vector. 
         FIG. 20  shows the pET DLS Plasmid (8.2 kb) which contains the  C. thermocellum  hydrogenase and the HydE accessory protein. 
         FIG. 21  shows the pCDF Plasmid (6.3 kb) which contains the  C. thermocellum  HydF and HydG accessory proteins. 
         FIG. 22  shows the list of primers (SEQ ID NOS: 210-223) used for isolating the  C. thermocellum  accessory proteins from a genomic DNA preparation and for cloning the proteins into the pET DLS and pCDF plasmids. 
         FIG. 23  illustrates the reaction mechanism for the reaction of molecular hydrogen with the redox dye, resazurin. 
         FIG. 24  are before and after pictures of four cuvettes containing resazurin. 
         FIG. 25  shows the spectra of the various forms of resazurin and resorufin. 
         FIG. 26  illustrates the effects of light on two algal strains. 
         FIG. 27  is a diagram of the metronidazole method of selection for algae with an enhanced hydrogenase. 
         FIG. 28  illustrates the photosynthetic electron transport chain, known as the “Z-scheme”. 
         FIG. 29  shows the proton exchange membrane fuel cell (PEM FC). 
         FIG. 30  shows the reactions that occur in a PEM fuel cell. 
         FIG. 31  is a diagram of the four-iron four-sulfur (4Fe-4S) cluster and the two-iron two sulfur (2Fe-2S) active site cluster that is present in Fe-only hydrogenases. 
         FIG. 32  shows the ratio of positive to negative electrostatic potential surface area (EPSA) plotted against hydrogen production for bacterial mutant hydrogenases. A molecular probe radius of 1.4 Å was used to calculate the EPSA ratios. 
         FIG. 33  is a graph of positive to negative EPSA ratios plotted against hydrogen production for bacterial mutant hydrogenases. A molecular probe radius of 1.0 Å was used to calculate the EPSA ratios. 
         FIG. 34  is a graph of positive to negative EPSA ratios plotted against hydrogen production for algal mutant hydrogenases. A molecular probe radius of 1.0 Å was used to calculate the EPSA ratios. 
         FIG. 35  is a graph of hydrogen production of bacterial hydrogenase mutants assessed by batch reactor headspace hydrogen count following after several hours of production. 
         FIG. 36  a graph of hydrogen production by algal hydrogenase mutants assessed by batch reactor headspace hydrogen count following several hours of production. 
     
    
    
     DETAILED DESCRIPTION 
     The present disclosure is based in part on the surprising demonstration of methods for producing photosynthetic alga capable of enhanced hydrogen production relative to wild type algae. The methods also involve production of genetically-modified bacteria which produce hydrogen, and full-length chimeric libraries of mutant hydrogenases. To produce the alga and bacteria, a technique known as directed evolution is used, whereby mutations are introduced into the DNA of hydrogenases native to the organism, and these tested for an increased rate of hydrogen production. 
     A. Definitions 
     Section headings as used in this section and the entire disclosure herein are not intended to be limiting. 
     As used herein, the singular forms “a,” “an” and “the” include plural referents unless the context clearly dictates otherwise. For the recitation of numeric ranges herein, each intervening number there between with the same degree of precision is explicitly contemplated. For example, for the range 6-9, the numbers 7 and 8 are contemplated in addition to 6 and 9, and for the range 6.0-7.0, the numbers 6.0, 6.1, 6.2, 6.3, 6.4, 6.5, 6.6, 6.7, 6.8, 6.9 and 7.0 are explicitly contemplated. 
     As used herein, the term “about” refers to approximately a +/−10% variation from the stated value. It is to be understood that such a variation is always included in any given value provided herein, whether or not it is specifically referred to. 
     Unless defined otherwise herein, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. 
     The following terms and abbreviations will have the following meanings throughout this disclosure: 
     aa: amino acid. 
     bp: base pair. 
     nt: nucleotide. 
     kD: kilodalton. 
     PCR: polymerase chain reaction. 
     RT-PCR: reverse transcriptase polymerase chain reaction. 
     ssDNA: single stranded deoxyribonucleic acid. 
     DNase I: enzyme that degrades DNA. 
     Restriction enzymes: enzymes that cleave DNA at a specific sequence. XhoI, PacI, NdeI, NheI, EcoRI, PspOMI, AsiSI are the names of restriction enzymes that cleave unique sequences which are listed in the New England Biolabs catalog. 
     cDNA: copy DNA, the same DNA as the original gene of interest, except all the introns, or non-coding material, have been removed. 
     LB (or NZY) and TAP: nutrient broths that are used to grow bacteria and algae, respectively Vector or Plasmid: circular DNA that can be transformed (inserted) into cells to express a gene of interest from its promoter. 
     HydA1 and HydA2: hydrogenase A1 and A2 in  Chlamydomonas reinhardtii.    
     UTR: untranslated region (of the DNA). 
     Intron: an untranslated region of a gene within a translated region. 
     GOI: gene of interest. 
     Hydrons: hydrogen atoms with two electrons (H). 
     pGenD: the name of a plasmid (p) that expresses some cDNA in algae. 
     pSMP1: the name of a plasmid derived from pGenD that expresses hydrogenase cDNA and has a PspOMI restriction site immediately after the HydA1 leader sequence and before the coding sequence for the HydA1. 
     pSMP1c: the same as pSMP1, but the PspOMI restriction site was added 20 nts downstream from the one in pSMP1. 
     pSMP2: the same as pSMP1, but an AsiSI restriction site instead of a PspOMI site was added 45 nts downstream of the HydA1 leader sequence. 
     IPTG: isopropyl-beta-D-thiogalactopyranoside, a chemical that artificially induces expression SDS-PAGE: sodium docecyl sulfate-polyacrylamide gel electrophoresis, a technique used to separate a mixture of multiple proteins. 
     Western Blot: in combination with SDS-PAGE, a technique used to identify one protein from a mixture of multiple proteins. 
     StEP: staggered extension process. 
     ITCHY: iterative truncation for the creation of hybrid enzymes. 
     RATCHITT: random chimeragenesis on transient templates. 
     DOGS: degenerative oligonucleotide gene shuffling. 
     B. Genetically Modified Alga and Bacteria with Enhanced Hydrogen Production 
     Hydrogenases are enzymes that generate hydrogen by combining protons with electrons from the photosynthetic electron transport chain. Molecular hydrogen is then released into the environment. It would be very beneficial to have a biological source of molecular hydrogen, in particular a photosynthetic organism capable of an increased rate of hydrogen production, to use as an energy source that reduces emissions of carbon dioxide, nitrous oxides, and/or methane. For example, fuel cells using an organism capable of sufficient hydrogen generation would efficiently generate power for numerous uses. Molecular hydrogen is the ideal fuel for use in fuel cells, if it could be produced at a cost that is competitive with current sources of energy. 
     Biohydrogen production from photosynthetic algae has the potential to be a viable alternative to hydrogen production from fossil fuels. It would not produce greenhouse gases (H 2 O+sunlight→O 2 +H 2 ); in fact, algae, like most plants, utilizes carbon dioxide for cellular growth, so it would serve as a carbon sink. In addition, a bioreactor would not produce toxic waste, but only algae and wastewater, similar to a fish tank. Also, a bioreactor would likely be about the size of an air conditioner and survive on low amounts of sunlight, so it would occupy a small amount of space and it could be located anywhere. 
     However, present commercial photosynthetic hydrogen production is not viable because of two major problems that prevent the hydrogenase from producing useful amounts of molecular hydrogen. First, the hydrogenase has a short half-life that prevents it from producing hydrogen for longer than a minute. Second, it is necessary that the hydrogenase be tolerant of oxygen. Since all known hydrogenases have a short half-life even in the presence of very low concentrations of oxygen, only a modified hydrogenase with increased hydrogen production and/or decreased oxygen sensitivity will allow for the commercial production of photosynthetically generated hydrogen. 
     Despite the evident challenges, the green alga  Chlamydomonas reinhardtii  ( C. reinhardtii ) has substantial potential as a candidate hydrogen producer. Each  C. reinhardtii  hydrogenase is capable of generating 6000-9000 molecules of molecular hydrogen per second. Once sustainable, a mole of hydrogenases, producing hydrogen at this rate, would generate enough hydrogen to fill the Graf Zeppelin in 10 minutes, or the main tank of the space shuttle in just 2 hours. In addition,  C. reinhardtii , is a common lab research organism whose genome has been sequenced. 
     Additionally, the crystal structures of homologous bacterial Fe-only hydrogenases have recently been described, which has provided insight into how oxygen irreversibly inhibits the enzyme. The active site for the production of hydrogen (2H+2e → H 2 ) is protected by its location deep within the center of the mature hydrogenase ( FIG. 3 ). Each of the reactants and the resulting hydrogen has a channel to the surface. The hydrogenase structures all have different ferredoxin binding motifs, but they all exhibit a chain of four iron, four sulfur clusters (4Fe-4S) that pass the electrons from the surface of the hydrogenase to the active site ( FIG. 1  and  FIG. 3 ). The modular 4Fe-4S clusters are separated by 1-1.5 nm intervals that allow for efficient electron transport to the active site. The protons pass through a putative second channel to reach the active site. This channel is lined with amino acids and protein-bound water molecules that are capable of binding the protons as they pass down the channel. Two putative channels exist for the release of molecular hydrogen resulting from the enzymatic reaction of the protons with the electrons at the active site ( FIG. 3 ). Unfortunately, these channels also allow for the passage of the larger oxygen molecule, which irreversibly inhibits the hydrogenase, probably by oxidizing an iron (II) to an iron (III) in the active site cluster. In brief then, photosynthesis generates oxygen as well as protons and energized electrons, and hydrogen production is dependent upon the photosynthetic process. Hence, an evolutionary enigma lies in the extreme sensitivity of hydrogenase to oxygen. Enzymatic hydrogen production is a short-term shunt to rid the cell of excess electrons. Therefore, the enzyme has probably been exposed too little if any selective pressure, which infers that its evolutionary potential is likely untapped. 
     The present disclosure provides the solution to the two major barriers preventing commercialization of photosynthetic hydrogen production to date: 1) the hydrogenase enzyme that produces the molecular hydrogen is sensitive to oxygen and, 2) the amount of hydrogen production needs to be increased. Since the  C. reinhardtii  hydrogenase genes have been cloned, mutation of the original parental genes is disclosed herein, as well as a method to search among a library of such mutants and select a mutant algal hydrogenase with an improved phenotype with respect to hydrogen production rate. After the mutant genes are created, this disclosure further describes insertion and expression of the mutants in  C. reinhardtii . A novel method of selecting or screening the mutants for enhanced traits in algae is also disclosed. The present disclosure further provides for a chimeric hydrogenase with a mutation(s) that decreases the diameter of the proton channels described above, thereby restricting the larger oxygen molecules by size exclusion, while still allowing for passage of the smaller bio-hydrogen molecules. 
     As will be further apparent from the Examples provided herein below, the present disclosure encompasses methods for mutating bacterial hydrogenases to produce increased amounts of hydrogen after just one round of a combinatorial shuffle. Additionally, exogenous DNA of hydrogenases can be successfully re-introduced into the genome of  C. reinhardtii . Once transformed, the plasmid containing the hydrogenase gene can be successfully expressed at both the RNA and protein levels, thus illustrating that mutated  C. reinhardtii  hydrogenase genes can also be inserted and expressed. Therefore, the present disclosure provides for expression of chimeric hydrogenases in  C. reinhardtii , and expression of chimeric algal hydrogenases containing mutations with improved function. Further disclosed then is a method for using these chimeric algal hydrogenases, when expressed for example in  C. reinhardtii , for photosynthetic hydrogen production. 
     Additionally, disclosed is a model system by which different mutations can readily be created and tested. Further analysis of the hydrogen production of the remainder of the existing mutant library in conjunction with molecular modeling determines what characteristics are shared between the mutants with improved hydrogen production. Considerable improvement frequently results from the comparison and combination of mutations isolated in the first round of shuffling. In addition, future rounds of directed evolution might benefit by the inclusion of the other four Clostridial parent sequences. In order to link the photosynthetic transport chain with an improved hydrogenase, the mutations that result in the most improved bacterial mutant should be re-created in the algal hydrogenase. Since the algal hydrogenase has a disparate codon bias and is merely homologous, not identical, to the bacterial hydrogenase, mutations in the algal enzyme need to be created and tested to determine if the same improvement is realized. Ultimately, an improved chimeric algal hydrogenase, transformed into  C. reinhardtii , is the basis for an economically viable method of hydrogen production. 
     Clostridal hydrogenases were selected because they share homology with their eukaryotic algal  C. reinhardtii  cousins, so that mutations resulting in an enhanced clostridial phenotype had an increased likelihood of producing the same enhancement when the analogous mutation was created in  C. reinhardtii . Therefore, the present disclosure also describes the creation of libraries of mutated clostridial hydrogenase genes, and their subsequent testing for enhanced hydrogen production. 
     To create the chimeric mutant hydrogenase libraries using homologous bacterial hydrogenases, a process known as directed evolution was applied. The process mimics natural evolution in that multiple mutations are created. Albeit on a shorter timescale, proteins of interest are submitted to repeat cycles of evolutionary pressure to produce a variety of mutations. Resulting mutants are then tested for desirable traits, such as in the present case, increased rate of hydrogen production relative to wild type hydrogenases. Mutants demonstrating desirable traits are then transformed into an algae such as  C. reinhardtii.    
     Bacterial hydrogenases with enhanced hydrogen production can be generated using degenerative oligonucleotide gene shuffling (DOGS). (For a description of the DOGS technique, see, for example, M. D. Gibbs et al., Gene 271(1):13-20 (2001), and WO/2002/018629, the entire disclosures of which are incorporated by reference in their entireties). Briefly, DOGS involves assigning at least one segment of at least one gene, typically at least two homologous genes, based on regions of encoded amino acid sequence; amplifying the one or more assigned segments of the gene(s) using primers specific for each segment; and combining, or causing recombination of the one or more amplified segments to form a mutant or chimeric gene. An oligonucleotide primer suitable for use in gene shuffling can be used, wherein the primer has a non-degenerate core based on a segment or template of a gene to be amplified, and the core is flanked by both 5′ and 3′ degenerate ends. As described herein, the method may involve forming a mutant or chimeric gene from two or more genes, in which one or more segments of each gene are assigned based on regions of encoded amino acid sequence; amplifying the one or more segments of the genes; and combining at least some of the amplified segments so as to form a mutant or chimeric gene. The two or more genes may belong to the same gene family encoding the same functional protein. According to the present disclosure, the two or more genes are homologous genes that each encode a hydrogenase. 
     Suitable bacterial parental hydrogenases include but are not limited to those obtained from Clostridial bacteria, such as for example  Clostridium saccharobutylium, Clostridium acetobutylicum, Clostridium pasteurinum, Clostridium perfringens, Clostridium tetani , and  Clostridium thermocellum.    
     The DOGS technique may be applied by dividing each parental hydrogenase gene from the bacteria into the multiple segments of roughly equal length. As described herein, the homologous parental genes are divided into multiple, e.g. eight (8), segments of roughly equal length. The segment boundaries are specifically chosen with positions within regions of high homology amongst the two parental genes. Except for the first and last segments (which in examples described herein are segments #1 and #8) each segment has two unique sticky ends created by a restriction endonuclease such as, in non-limiting example, SapI. The uniqueness of each sticky end allows it to overlap only with its neighboring segments. For example, as described herein, the downstream end of segment #2 could only overlap with the upstream end of segment #3 and the downstream end of segment #3 would overlap with only the upstream end of segment #4, thereby generating a full-length chimeric gene library with, segments from both parental genes. 
     A variation to the DOGS technique may also be used as in which only three base pairs of homology are required as the overlap is created by a restriction enzyme thus allowing elimination of the overlap PCR steps in favor of a simple annealing of the sticky ends generated by the enzyme. 
     Either version of the DOGS technique can be used to produce a full-length chimeric gene library wherein the overlap positions occur in regions of high homology. The gene shuffling technique effectively swaps segments of the parental genes by choosing the segment intersections in advance. It has the advantage of capitalizing on the evolutionary differences amongst the parent genes including the regions of high homology while maintaining the overall length of the gene. 
     As is known in the art, a number of different programs can be used to identify whether a nucleic acid or polypeptide has sequence identity or similarity to a known sequence, after aligning the sequences and introducing gaps, if necessary, to achieve the maximum percent sequence similarity. The sequence alignment for purposes of determining percent amino acid sequence similarity may be achieved in any of a number of established various ways that are known to those of routine skill in the art. Readily available computer software can be used to determine sequence similarity. For example, BLAST, gapped BLAST, BLAST-2, ALIGN or Megalign (DNASTAR) software is generally widely known and used in the art (see, e.g., Altschul et al., J. Mol. Biol. 215, 403-410, (1990); Altschul et al., Methods in Enzymology, 266, 460-480 (1996); Altschul et al., Nucleic Acids Res. 25, 3389-3402 (1997)). A hydrogenase mutant polypeptide will have at least about 80% amino acid sequence similarity, alternatively at least about 60%, 70%, 75%, 80%, 81%, 82%, 83%, 84%, 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, or 99% amino acid sequence similarity, to a full-length wild type hydrogenase (control) sequence as disclosed herein. 
     Expression of mutant hydrogenases can be achieved in a suitable host organism such as but not limited to a bacterial host such as  E. coli , and a green algae such as  C. reinhartdtii , according to techniques well known and described in the art. A wide variety of molecular and biochemical methods and tools including suitable vectors are available for generating mutant hydrogenases, transforming host organisms and expressing mutant hydrogenases as described herein. See, e.g., M OLECULAR  C LONING , A L ABORATORY  M ANUAL  (Sambrook et al., Cold Spring Harbor Laboratory); C URRENT  P ROTOCOLS IN  M OLECULAR  B IOLOGY  (Eds. Ausubel et al., Greene Publ. Assoc., Wiley-Interscience, NY). Alga used according to the present disclosure may be any alga capable of producing hydrogen. Preferably a green or blue-green alga, such as but not limited to  C. reinhardtii , is selected. 
     Hydrogen production by mutant dehydrogenases can be measured using any established method, such as but not limited to the methyl viologen assay, the resazurin assay and the metronidazole assay as described in further detail in the Examples herein below. Alternatively, a batch reactor can be employed, and the headspace above the reaction mixture sampled periodically, for example every hour to obtain a hydrogen content in the headspace using gas chromatography. Efficiency of each mutant can be determined for example by comparing the time required for a mutant to produce the maximum amount of hydrogen obtained from the wild type hydrogenase under standard conditions. 
     The present disclosure describes calculation of a positive to negative electrostatic potential surface area (EPSA) ratio for each mutant hydrogenase. As described herein, the positive to negative EPSA is indicative of hydrogen production capability, wherein bacterial and algal mutants with a positive to negative EPSA ratio from about 1 to about 115, preferably from about 2 to about 50, more preferably about 5 to about 20, yield hydrogen at a faster rate than that observed with a wild type control organism. A positive to negative EPSA ratio for any hydrogenase mutant may be determined according to the method described in detail in Examples 5-7 herein below. Briefly, amino acid sequences are used to construct alpha-helix structures for each mutant. These structures are then energy minimized using an established energy optimization or geometry optimization algorithm, such as but not limited to the simple gradient algorithm, the nonlinear conjugate gradient algorithm, or preferably the physics-based all atom force field for proteins (OPLS-AA). (See, e.g., W. Jorgensen et al., J. Am. Chem. Soc. 118(45): 11225-36 (1996)). Next, the positive and negative electrostatic potential surface areas (EPSA) for each structure are calculated using a selected probe radius of 0.5 to 1.5 Å. The probe radius must be selected to yield the most accurate correlation between the EPSA ratio and hydrogen production as described below. Actual values in the calculations described herein were either 1.0 and 1.4 Å. Overlapping positive and negative regions may or may not be excluded in the calculations. Computations may be performed with computational chemistry software, such as but not limited to the HyperChem 7.5 computational chemistry program. In the calculations herein, a 32 bit computer was used. Each positive to negative EPSA ratio thus obtained may be plotted versus experimentally measured hydrogen production for any mutant, and compared to that observed with a wild type control. Such a plot may be best fit with a Log Normal Gaussian function, which can be used as the prediction model. The selected prediction model suggests that the ratio of positive to negative electrostatic potential surface areas is a measure of how a given protein mutant will fold into its tertiary state and how the final tertiary state affects hydrogen production. Initial results showed that an EPSA ratio in the range of about 15 to about 115, obtained with a probe radius of 1.4 Å, is required to achieve hydrogen production from bacterial hydrogenases greater than that obtained with a wild type bacterial control, and that maximum hydrogen production is achieved with a positive to negative EPSA ratio of about 42 (see  FIG. 32 ). Later more accurate results, using a probe radius of 1.0 Å, showed that an EPSA ratio in the preferable range of about 2 to about 50, and more preferably from about 5 to about 20, to be indicative of the highest rates of hydrogen production observed from bacterial hydrogenases, relative to a wild type control. (See  FIG. 33 ). For example, the model predicts that an EPSA ratio of about 16 will produce hydrogen from a bacterial hydrogenase at a rate of about 530 times that of a wild type hydrogenase. 
     C. Adaptations of the Methods of the Present Disclosure 
     By way of example, not of limitation, examples of the present invention shall now be given. 
     Example 1 
     Generation of Mutant and Chimeric Hydrogenase Libraries 
     Gene shuffling is the process of creating mutant DNA sequences which codes for a protein with an improved or a novel function. Most frequently, the DNA sequence represents a gene and the goal is to modify it by mutating it into many thousands of new genes, each mutated in a different manner, to create a new phenotype with the desired characteristics. A sizeable percentage of these genes will encode for non-functional proteins or for proteins with no improvement over the original. However, several mutant or chimeric proteins will exhibit a significantly improved functionality. The proteins with improved function can be re-shuffled, thereby amplifying the improved function. Thousands or millions of mutants can be created and all of “sequence space” can, in principle, be successfully surveyed ( FIG. 4 ). 
     The sequence space (20 n) is the set of all 20 possible (common) amino acids in each position for a protein of a particular length. (n). A plot of the sequence space versus the desired trait shows local maxima and minima in the desired trait. After several rounds of shuffling, it is common to find that the protein cannot be improved any further. Protein function improvements of 500×−32,000× are known to result from this method in other organisms, thereby confirming that this method can produce significant positive changes in a protein. Consequently, we can use gene shuffling to create a library of chimeric proteins. By selecting for an improved trait from that library, we can direct the evolution of the original parental protein. This technique is referred to herein as “Directed Evolution”. 
     Gene shuffling using Willem Stemmer&#39;s family shuffling technique can be used to effect direct evolution. In Stemmer&#39;s protocol, multiple parent genes were digested. The fragments were recombined using multiple cycles of PCR to form chimeric progeny ( FIG. 5 ). Other techniques soon followed, including: family shuffling with single stranded DNA (ssDNA), staggered extension process (StEP), iterative truncation for the creation of hybrid enzymes (ITCHY), random chimeragenesis on transient templates (RACHITT), and degenerative oligonucleotide gene shuffling (DOGS). Each technique generates a large number of diverse gene sequences, referred to as a library. A selection is then performed on the chimeric proteins that are generated from the library of genes; a selection for an enhanced trait such as enzyme selectivity, stability, or activity. For algal hydrogenase research, organisms with enhanced tolerance for oxygen and/or an increased capacity for the production of hydrogen are favored in the selection process. As described herein, this directed evolution may also be achieved with gene shuffling via a technique known as error-prone PCR. This method of creating a mutant library merely alters the conditions of standard PCR so that they were no longer optimal, thereby forcing the polymerase to make mistakes. This method in practice proved to be not sufficiently random and resulting mutant offspring were not significantly different from the parent genes. 
     In order to find the optimal method to generate a library, several of the above shuffling methods can be used, alone or in combination. The RACHITT method generates a library containing a high percentage of diverse chimeras with little or no representation of the original parental genes ( FIG. 6 ). It does so by choosing one of the parental genes to be a template strand, while the others are designated as donor strands. All the parental DNA is made single-stranded (ssDNA). In addition, the donor strands are digested into pieces of about 200-800 nt. The digested donor DNA oligonucleotides bind to the template strand. Once the gaps are filled in, the template strand is removed by completely digesting it down to the level of individual bases. Since all of the parental DNA was digested, a library of full-length chimeras is created. Next, the standard family shuffle ( FIG. 5 ) was used. Error-prone PCR and the DOGS technique, as described below, are particularly useful in creating chimeric libraries. 
     In order to find the optimal method to generate a library of chimeric hydrogenases, several shuffling methods can be tried, alone or in combination. Here, to create a library of chimeras, several mutational techniques were used. The family shuffle and random chimeragenesis on transient templates (RACHITT) were used initially, but were found to be less effective than other methods. However, two other techniques, random mutagenesis by error-prone PCR and a variation of degenerative oligonucleotide gene shuffling (DOGS) were highly successful and led to the creation of the mutational libraries disclosed herein. 
     Error-prone PCR is similar to standard PCR in that it uses the same primers to replicate the gene of interest (GOI). However, the polymerase has been altered so that errors in the replication process occur more often. This technique generated a mutant library with an average of 10 point mutations per gene. Although libraries of mutated genes are readily created, error-prone PCR is not the most efficient technique for the creation of a protein with improved function, as it is not known which or where a point mutation would be beneficial. 
     The DOGS technique divides each parental GOI into roughly equal segments (see  FIG. 13 ). In this study, the segment boundaries were specifically chosen within regions of high homology amongst the parent genes. Each segment of each gene was generated by an initial PCR. Forward and reverse degenerate primers were employed on both sides of the segment boundary in order to facilitate an overlap of the two segments (see  FIG. 14 ). Each segment boundary is a unique sequence. A second PCR combined the segments in the same segment order as the original parent genes, thereby generating a full-length chimeric library. 
     A variation to the DOGS technique was also used; it required only three base pairs of homology as the overlap is created by a restriction enzyme (see  FIG. 15 ). The overlap PCR steps are eliminated in favor of a simple annealing of the sticky ends generated by the enzyme. 
     Either version of the DOGS technique lead to a full-length chimeric gene library wherein the overlap positions occur in regions of high homology. Therefore, this shuffling technique effectively swaps segments of the parental genes by choosing the segment intersections in advance. It has the advantage of capitalizing on the evolutionary differences amongst the parent genes including the regions of high homology while maintaining the overall length of the gene. 
     Two mutagenesis libraries were created. The first library was generated using a random mutagenesis kit (Stratagene GeneMorph II: #200550) according to the manufacturer&#39;s protocol. The PCR protocol that generated the highest number of mutations was followed: an initial denaturation of 2′ at 95° C. followed by 30 cycles of [0.5° at 95° C., 0.5° at 59° C. 2′ at 72° C.]. 
     A second mutagenesis library was created using a variation of the original DOGS method (see  FIGS. 13 through 15 ). The total length of the gene was first divided into eight segments of equal length. The segment lengths varied from 160 bp to 290 bp, as it was necessary to choose the segment boundaries within regions of homology (see Table 1, SEQ ID NOs. 23-38 (middle column); SEQ ID NOs: 39-54 (right hand column)). Each segment was individually generated by PCR using 30 ng of the plasmid containing the parent hydrogenase. The PCR protocol utilized Phusion polymerase (New England Biolabs (NEB): #F530L) according to the manufacturer&#39;s protocol: an initial denaturation of 30 sec. at 98° C. followed by 35 cycles of [10 sec. at 98° C., 20 sec. at 57° C. Δ7° C., 20 sec. at 72° C.], and a final extension of 10′ at 72° C. The segments were gel purified (Qiagen Qiaex II kit: #20021) and digested with SapI (NEB: #R0569L). 
     The digested segments were column purified (Qiagen Qia-Quick PCR kit: #28104) and 50-100 ng of each segment was mixed together and ligated (Lucigen T4 ligase: #30025-2). To generate the full-length hydrogenases, the ligated segments were amplified by PCR using 1-2 μl of the ligation mixture and primers that complemented the homologous regions just upstream and downstream of segment numbers one and eight, respectively. 
     The SapI variation of the DOGS method could generate ligations of up to four segments, so the library was ultimately created using a combination of the original DOGS method and the SapI variation. The SapI variation on the DOGS method was utilized, as above, to combine the eight individual segments, from both parental genes, into four ¼-length chimeras, each with a combination of two segments, i.e. all the possible combinations of segment 1 and segment 2 were generated as well as segment 3 with segment 4, etc. After ligating the SapI generated segments, the ¼-length chimeras were amplified by PCR using primers appropriate for the particular ligation. 
     The appropriate primers (see Table 1, SEQ ID NOs. 23-38 (middle column); SEQ ID NOs: 39-54 (right hand column)) had extra bases added (see Table 2, SEQ ID NOs. 55-66) that included 10-15 bp of overlap with the ¼-length chimeras on either side. The quarter genes were purified (Qiagen Qiaex II kit: #20021) and then re-amplified into ½-length or full-length chimeras using a two-step program that first allowed the individual ¼-length chimeras to overlap before the outside primers generated the full-length gene: an initial denaturation of 30 sec. at 98° C. followed by two cycles of [10 sec. at 98° C., 30 sec. at 60° C., 30 sec. at 72° C.], a pause at 93° C., 35 cycles of [10 sec. at 98° C., 20 sec. at 59° C., 40 sec. at 72° C.], and a final extension of 10′ at 72° C. 
     The first two cycles allowed the segments to self-prime and extend from the overlapping regions. The outside primers were added during the pause, which lasted approximately three minutes before the full-length chimeras were generated. The resulting full-length chimeras were gel purified, digested with NcoI &amp; AscI (NEB: #R0193S and #R0558S), column purified, and cloned into the pET DLS vector. The above procedure, of digestion with SapI followed by ligation and PCR amplification, can be repeated several times, as necessary, to obtain a full-length chimeric hydrogenase. 
     
       
         
               
             
               
               
               
               
               
               
             
           
               
                 TABLE 1 
               
             
             
               
                   
               
               
                 [SEQ ID NOs. 23-38 (middle column); SEQ ID NOs: 39-54 (right hand column)] 
               
             
          
           
               
                   
                   
                 SapI 
                   
                 Upstream &amp; Downstream Primers 
                 Upstream &amp; Downstream Primers 
               
               
                   
                 Segment 
                 Generated 
                 Overhang 
                 for 
                 for 
               
               
                 Segment 
                 Length 1   
                 Overhang 2   
                 Location 3   
                 
                   C. saccharobutylicum 
                   4 
                 
                 
                   C. acetobutylicum 
                   4 
                 
               
               
                   
               
               
                 1 
                 211 
                 ggt 
                 212-214 
                 5′ cccacgccgaaacaag 3′ 
                 5′ cccacgccgaaacaag 3′ 
               
               
                   
                 (462) 
                   
                   
                 5′ accataccttcttcaacttttg 3′ 
                 5′ accattccatcttcaactttggc 3′ 
               
               
                   
               
               
                 2 
                 229 
                 tgt 
                 441-443 
                 5′ ggtagtaaaaacaaactcag 3′ 
                 5′ ggtaataaacacagaatccgatg 3′ 
               
               
                   
                   
                   
                   
                 5′ acacttagttctgtctattacaat 3′ 
                 5′ acattttgatctgtcaattacaat 3′ 
               
               
                   
               
               
                 3 
                 231 
                 cat 
                 672-674 
                 5′ tgtgtgctatgcggaaga 3′ 
                 5′ tgtgtactatgcggtagatg 3′ 
               
               
                   
                   
                   
                   
                 5′ atgcttattaggatcttctaatg 3′ 
                 5′ atgttttttagggtcattaagag 3′ 
               
               
                   
               
               
                 4 
                 192 
                 aaa 
                 864-866 
                 5′ catgtaatagttgctatggc 3′ 
                 5′ catgtcattgttgcaatggc 3′ 
               
               
                   
                   
                   
                   
                 5′ tttaactctttcaataaactctg 3′ 
                 5′ tttaactctgcctaaaagttcagt 3′ 
               
               
                   
               
               
                 5 
                 284 
                 aga 
                 1148-1150 
                 5′ aaaaataatggaccattccca 3′ 
                 5′ aaaaataatggcccattccctatg 3′ 
               
               
                   
                   
                   
                   
                 5′ tctcttgtagtaagtactgcatc 3′ 
                 5′ tctcttgtagttaaggatgcatcaa 3′ 
               
               
                   
               
               
                 6 
                 214 
                 ata 
                 1362-1364 
                 5′ agaattagcaaaaatgattaaagatg 
                 5′ agagcttgcaaaaatgattaaagatg 
               
               
                   
                   
                   
                   
                 caa 3′ 
                 c 3′ 
               
               
                   
                   
                   
                   
                 5′ tattccttgtaatcctcttatttgtg 
                 5′ tatgcctttaaagcctcttacttcag 3′ 
               
               
                   
                   
                   
                   
                 tat 3′ 
                   
               
               
                   
               
               
                 7 
                 159 
                 gga 
                 1521-1523 
                 5′ ataaaagaggctacagtagaaattgg 
                 5′ ataaaagaagcggaagttgaaattgc 3′ 
               
               
                   
                   
                   
                   
                 tgg 3′ 
                   
               
               
                   
                   
                   
                   
                 5′ tccaccatttacacatccgcc 3′ 
                 5′ tccaccatttatacatccaccag 3′ 
               
               
                   
               
               
                 8 
                 205 
                 N/A 
                 N/A 
                 5′ ggaggacaaccacacgta 3′ 
                 5′ ggaggtcaacctcacgtaaatg 3′ 
               
               
                   
                 (455) 
                   
                   
                 5′ tacgattactttctgttcgactta 3′ 
                 5′ tacgattactttctgttcgactta 3′ 
               
               
                   
               
               
                   1 The first and last segments were longer to allow for nested primers. The total length of the segment is listed in parenthesis. 
               
               
                   2 The overhang was generated at the end of the listed segment. 
               
               
                   3 The location in  C. Saccharobutylicum . 
               
               
                   4 Each of these primers had a SapI “gactgactGCTCTTCt” tail. 
               
             
          
         
       
     
     
       
         
               
             
               
               
               
             
           
               
                 TABLE 2 
               
             
             
               
                   
               
               
                 [SEQ ID NOs. 55-66] 
               
             
          
           
               
                 Overlapping 
                 Overlapping 
                 Forward/Reverse 
               
               
                 Segments 1   
                 Parental Segments 2   
                 Degenerate Overlap Primers 3   
               
               
                   
               
               
                 1&amp;2 with 3&amp;4 
                 Saccharo &amp; Saccharo 
                 5′ gtaatagacagaactaagtgtgtgctatgcggaagatgt 3′ 
               
               
                   
                 Saccharo &amp; Aceto 
                 5′ gtaattgacagatcaaaatgtgtgctatgcggaagatgt 3′ 
               
               
                   
                 Aceto &amp; Aceto 
                 5′ gtaatagacagaactaagtgtgtactatgcggtagatgcg 3′ 
               
               
                   
                 Aceto &amp; Saccharo 
                 5′ gtaattgacagatcaaaatgtgtactatgcggtagatgcg 3′ 
               
               
                   
               
               
                 3&amp;4 with 5&amp;6 
                 Saccharo &amp; Saccharo 
                 5′ cagagtttattgaaagagttaaaaataatggaccattcccaatg 3′ 
               
               
                   
                 Saccharo &amp; Aceto 
                 5′ cagagtttattgaaagagttaaaaataatggcccattccctatg 3′ 
               
               
                   
                 Aceto &amp; Aceto 
                 5′ ctgaacttttaggcagagttaaaaataatggcccattccctatg 3′ 
               
               
                   
                 Aceto &amp; Saccharo 
                 5′ ctgaacttttaggcagagttaaaaataatggaccattcccaatg 3′ 
               
               
                   
               
               
                 5&amp;6 with 7&amp;8 
                 Saccharo &amp; Saccharo 
                 5′ aaataagaggattacaaggaataaaagaggctacagtagaaattg 3′ 
               
               
                   
                 Saccharo &amp; Aceto 
                 5′ aaataagaggattacaaggaataaaagaagcggaagttgaaattgc 3′ 
               
               
                   
                 Aceto &amp; Aceto 
                 5′ gtaagaggctttaaaggcataaaagaagcggaagttgaaattgc 3′ 
               
               
                   
                 Aceto &amp; Saccharo 
                 5′ gtaagaggctttaaaggcataaaagaggctacagtagaaattg 3′ 
               
               
                   
               
               
                   1 The original parental genes were divided into eight segments. 
               
               
                   2 The parents are  C. saccharobutylicum  (Saccharo) and  C. acetobutylicum  (Aceto). 
               
               
                   3 The reverse primers are simply the inverse complement of the listed forward primer. 
               
             
          
         
       
     
     Example 2 
     Hydrogenase Library 
     The error-prone PCR technique created a library of progeny genes that contained an average of 10 point mutations per hydrogenase gene. The DOGS technique created a library of 254 full-length hydrogenase genes, but each of the 254 chimeras was unique and contained sizeable amounts of shuffled genetic material. 
       C. reinhardtii  is known to contain two hydrogenases (HydA1 and HydA2), both of which produce molecular hydrogen. However, both the transcription of the RNA and the mature hydrogenase protein are sensitive to oxygen, a by-product of photosynthesis, which results in only a brief evolution of molecular hydrogen. However, mutations engineered into the native algal hydrogenase will generate chimeric hydrogenases with an increased rate of hydrogen production and/or an increased tolerance to the presence of oxygen. Each type of chimeric hydrogenase would result in an increased amount of hydrogen production. More importantly, subsequent rounds of mutation have the potential for further improvement. In order to successfully create a mutant version of  C. reinhardtii , the goal was to identify, clone, sequence, and express the hydrogenase genes that were responsible for the generation of hydrogen in algae. 
       C. reinhardtii  contains an enzyme (hydrogenase) that produces molecular hydrogen from electrons donated by ferredoxin, an enzyme in the photosynthetic electron transport pathway, and in fact, according to the present disclosure, it was determined that  C. reinhardtii  actually contains two hydrogenases, HydA1 and HydA2. 
     Example 3 
     Vector Creation 
     HydA1 and HydA2 were cloned and a vector that could express the hydrogenases in algae designed and tested. The transformation of heterogeneous DNA, via a vector, into  C. reinhardtii  was difficult due the distinct codon bias of the alga. Subsequent extraction of the transformed DNA was also difficult due to the alga&#39;s predilection for digesting the transformed DNA and incorporating it at random locations within the genome. 
     Probing the cDNA library successfully identified two hydrogenases, designated as HydA1 and HydA2, which were cloned and sequenced ( FIG. 7 ). HydA1 and HydA2 are similar in that they both have the ability to produce molecular hydrogen and are sensitive to molecular oxygen. HydA2 is 68% identical and 74% similar (i.e. a related amino acid) to HydA1. Although functionally similar and approximately the same length, HydA1 (495 aa) and HydA2 (505 aa) are distinct proteins from separate genes with different promoter regions. In addition, the leader sequences in HydA1 and HydA2 are quite distinct, which indicates that they may function in different locations of the chloroplast. First, the un-translated regions (UTR) immediately upstream (5′ UTR) and downstream (3′ UTR) of the gene were utilized to create a vector (pGenD) that is capable of driving the expression of cDNA in algae, including cDNA foreign to the alga. 
     Then apGenD vector was used to create a pGenD+Ble by inserting a Bleomycin antibiotic resistance gene (BLE) into XhoI/PacI downstream of the 3′ un-translated region. Antibiotic resistance confers selectivity for the alga containing the plasmid, allowing it to grow in the presence of Bleomycin antibiotic, while repressing the growth of bacteria or algae without the plasmid. 
     Then, the pGenD+Ble vector was altered to the new vector, named “pSMP”, by inserting the cDNA of one of the  C. reinhardtii  hydrogenases (HydA1) between the NdeI and NheI restriction sites thereby replacing the original coding sequence ( FIG. 8 ). A strepII tag was added at the 3′ terminus of the HydA1 gene between the NheI and EcoRI restrictions sites ( FIG. 8 ). The tag allows ready detection and differentiation of the HydA1 DNA, RNA, and protein from the indigenous hydrogenases. Lastly, a restriction site was created, so as to separate the HydA1 leader sequence from the DNA that codes for the mature protein. 
     The leader sequence contains cellular instructions for the placement of the mature protein in the correct location within the algal chloroplast. It was important to insert mutated coding sequences without changing the leader sequence. Since it was unknown whether or not an internal restriction site would impede the successful expression of the vector, three pSMP vectors (pSMP1, pSMP1c, and pSMP2) were created, each with a slightly different restriction site. In the case of pSMP1, a PspOMI restriction site was inserted, de novo, immediately downstream of the leader sequence ( FIG. 8A ). 
     As it was possible that even a small (6 bp) de novo insertion could interfere with expression, two additional pSMP vectors were created. The vector pSMP1c had a PspOMI site created from a conservative mutation about 20 nucleotides (nt) downstream of the end of the leader sequence ( FIG. 8B ). The vector pSMP2 had an AsiSI site created from a conservative mutation about 45 nt downstream from the end of the leader sequence ( FIG. 8C ). These vectors were tested for expression of HydA1 cDNA in  C. reinhardtii.    
     The vectors each containing the hydrogenase, HydA1, were transformed into  C. reinhardtii . Genomic DNA was extracted from 20 different transformations of the pSMP vectors. The hydrogenase DNA was amplified by polymerase chain reaction (PCR) with primers specific to the 3′ end of the hydrogenase and the transformed vector. The 800 base pair (bp) band indicated that each of the 20 genomic extracts received the plasmid ( FIG. 9 ). There was no difference in the transformation efficiency of the three vectors. 
     Twelve of the above transformants were tested for RNA expression using real-time reverse transcriptase PCR(RT-PCR) and four were positive for expression ( FIG. 10 ). Five of the twelve had been transformed with pSMP1 and two (A &amp; B in  FIG. 10 ) of the five were positive for expression (40%). In addition, three of the twelve were transformed with pSMP1c and two (C &amp; D in  FIG. 10 ) of the three were positive for expression (67%). The remaining four transformants were transformed with pSMP2 and all were negative. The pSMP1 and pSMP1c vectors are therefore useful for expressing the hydrogenases in  C. reinhardtii . Plasmid pSMP1c appears to be the best transformation vector given that two-thirds of the transformants were positive for expression and they had a higher level of expression than the pSMP1 transformants. 
     A transformant containing pSMP1 was tested for protein expression by Western blot ( FIG. 11 ). Two proteins, at approximately 49 kiloDaltons (kD) and 97 kD, were present in the transformed sample and not present in the non-transformed control. The 49 kD band is the correct molecular weight of the HydA1 protein and the faint band at 97 kD (data not shown) is likely a dimer. 
     The results show that pSMP1 and pSMP1c vectors are useful for expressing the hydrogenases in  C. reinhardtii . Plasmid pSMP1c appears to be the best transformation vector given that two-thirds of the transformants were positive for expression and they both had a higher level of expression than the pSMP1 transformants. Because the RT-PCR and Western blot were positive for expression insertion and expression of the wild-type hydrogenase (cDNA) in  C. reinhardtii  was accomplished. 
     In order to isolate the HydA1 and HydA2, a cDNA library was obtained (John Davies, Exelixis Inc.: San Francisco, Calif.), which had been created under anaerobic conditions and believed to contain the hydrogenases. The library consisted of viral particles, known as phage, which contained the cDNA of  C. reinhardtii . In order to detect cDNA that contained a hydrogenase, the phage library was allowed to infect XL-1 Blue bacteria (Stratagene: #211204 and #211203). The infected cells were then mixed with top agarose and plated on LB-agar plates at a density of 5000 plaques per plate. The phage-infected bacteria formed a lawn punctuated with plaques, clear spots where the cDNA containing phage have infected and lysed all the bacteria surrounding the initial infection. The pattern of clear spots on the plates were copied by overlaying the plate with a membrane (ISC BioExpress: #F-3193-82) for three minutes. Afterwards, the plate was sealed with Para film and saved at 4° C. until it could be determined whether any of the plaques contained a hydrogenase. The membrane was treated and washed three times (Promega Technical Bulletin #72) before the DNA was cross-linked to the membrane via exposure to UV light (.lamda.=365 nm). 
     The membrane, which then contained a copy of the plaques on the LB plate, was pre-hybridized (Sigma Perfect-Hyb: #H-3032) for 30 minutes at 68° C. to block non-specific binding sites, before it was hybridized to a chemiluminescent probe for 3 hours at 68° C. The chemiluminescent probe was created by binding a reporter molecule, digoxigenin, to the oligonucleotide of putative hydrogenase DNA mentioned above. First though, the oligonucleotide was purified and concentrated before it was labeled with digoxigenin-dUTP (Roche: #1-573-152) using a random labeling kit (MBI Fermentas: #K0621). The suspected hydrogenase portion of the probe binds homologous cDNA on the membrane, presumably hydrogenase containing cDNA. After the probe bound to the cDNA on the membrane, it was washed, blocked, and treated with an anti-digoxigenin antibody, which is bound to alkaline phosphatase (Roche: #1-093-274). Detection of the bound antibody was accomplished by providing two substrates to the alkaline phosphatase, which results in a purple color wherever the probe is bound to the membrane. (Roche: #1-175-041). 
     The purple spots on the treated membrane corresponded to plaques on the LB plate that contained cDNA, which bound the probe. The bound cDNA were possible hydrogenases, so the designated plaques were removed from the agar with a wide bore pipette tip. Phage were purified from the agar by incubating the plaques in 1 ml of phage buffer while shaking for three hours at room temperature. The purified phage were diluted and allowed to re-infect fresh XL-1 Blue bacteria. Afterwards, they were re-plated at an approximate density of 100 plaques per plate. The purification process was repeated, as above, by copying the new plates with membranes, treating the membranes with the chemiluminescent probe, and isolating the positive plaques. The initial positive plaques were considered “plaque purified” after performing three rounds of plating and purifying. 
     The cDNA were extracted from the purified phage by infecting XPORT bacteria (Stratagene: #211204) and plating with NZY top agarose on NZY plates. The plates were incubated overnight at 37° C. The phage was collected by overlaying the plate with 3 ml of SM buffer (recipe in Appendix D) for 30 minutes, centrifuging, and resuspending in 100 μl of SM Buffer. The concentrated cells were used to infect XLOLR cells (Stratagene: #211204) before plating on LB plates. As a result, cDNA was now contained as a plasmid within the XLOLR bacteria. 
     The  C. reinhardtii  hydrogenase HydA1 gene was amplified by PCR from the cDNA described above. Restriction sites (NdeI and NheI) were added to the primers to facilitate cloning into the vector. The pSMP1, pSMP1c, and pSMP2 plasmids were created by replacing the original gene (PsaD) with the HydA1 DNA. Both the pGenD+Ble plasmid and the HydA1 gene were digested with NdeI and NheI (New England Biolabs (NEB): NdeI: #R0111S and NheI: #R0131S), purified (Qiagen Qiaex II kit: #20021), and ligated (Lucigen T4 ligase: #30025-2). The strepII tag (IBA GmbH: Gottingen, Germany) was inserted, by PCR mutagenesis, at the 3′ terminus of the HydA1 gene between the NheI and EcoRI (NEB: #R0101S) restriction sites. The mutagenesis (Stratagene QuikChange Site-Directed Mutagenesis kit: #200519) was performed according to the manufacturer&#39;s protocol, except the extension time was increased to 2-3 minutes/kb in order to accommodate the large insertion (30 bp) ( FIG. 12 ). The PspOMI (NEB: #V0215S) restriction sites in pSMP1 and pSMP1c and the AsiSI (NEB: #R0630S) restriction site in pSMP2 were created by site directed mutation (Stratagene Quik-change Multi site directed mutagenesis kit: #200515) after an existing PspOMI site was removed from a different location. The PCR protocol and the mutational primers were designed according to the manufacturer&#39;s protocol ( FIG. 12 ). 
     All of the algal strains, including the pSMP1 transformants, were grown under oxygenic conditions in 50 ml of sterile TAP media utilizing 250 ml Erlenmeyer flasks. The flasks were shaken at 100-1.50 rpm and exposed to 100 μE (1 μE=1 μEinstein=1 μmoles photons/m 2 sec=500 foot candles) of light from standard fluorescent fixtures. 
     Each of the three pSMP vectors were transformed into the dw15-1 (courtesy of Barb Sears: Michigan State Univ.), a fast growing variant of the common  C. reinhardtii  cc425 strain, using a BTX ECM 630 electroporator with the following settings: 750 Volts, 25 and 1575 ohms. Immediately after electroporation, the algae were mixed with 3 ml of cornstarch suspended in TAP+60 μM sucrose and plated on Tris-Acetate-Phosphate (TAP) plates containing 10 μg/ml of the antibiotic Zeocin (140). The resulting transformants were grown, as above, for one week before they were transferred to standard TAP plates containing no antibiotics. The transformants were transferred to plates without antibiotics as Zeocin is a presumed mutagen and antibiotics are not necessary to maintain the transformants, since the plasmid DNA is incorporated into the alga&#39;s genome. 
     Algal transformants containing the pSMP1, pSMP1c, or pSMP2 vectors, were inoculated from plates into liquid culture. The cultures were harvested at chlorophyll (Chl) concentrations between 15-18 μg Chl/ml (exponential growth phase) and all 50 ml was centrifuged at 4000 g for 5 minutes. The genomic DNA was isolated (Qiagen DNAeasy Plant kit: #69104) and PCR was performed on the genomic preparations from each of the three different pSMP transformations. The primers were specific to the 3′ end of the hydrogenase cDNA. The PCR conditions were: 1′ initial denaturation at 95° C., 35 cycles of [30 sec. at 95° C., 30 sec at 60° C., and 2′ at 72° C.], and 10′ final extension 72° C. The upstream primer (5′ CCAGCTGCTGCCAGAATTC 3′, SEQ ID NO. 16) and the downstream primer (5′ CCAGCTGCTGCCAGAATTC 3′, SEQ ID NO. 13) amplify the last 800 bases at the 3′ end of the HydA1 hydrogenase. 
     Cells were harvested at chlorophyll (Chl) concentrations between 18-25 μg Chl/ml by centrifuging 10 ml of algal culture at 2000 g for 1 minute. RNA was obtained (Qiagen RNAeasy kit: #74104) and treated with a DNaseI (Ambion Turbo DNA-free kit: #1907) to remove the residual plasmid DNA. The cDNA was generated from the purified RNA (Qiagen Quantitect Reverse Transcription kit: #205311). Real-Time RT-PCR was performed on all of the samples using an upstream primer (5′ GACGAGAAGAAGGCTAGCGC 3′, SEQ ID NO. 15) specific to the 3′ end of the hydrogenase and a downstream primer (5′ CCAGCTGCTGCCAGAATTC 3′, SEQ ID NO. 14) that was specific to the strepII tag (IBA GmbH: Germany). DyNAmo SYBR Green (Applied Biosystems: #4367660) on an Applied Biosystems 7300 Real Time PCR system was used. The real time RTPCR program was: an initial dissociation of 10 minutes at 95° C. followed by 40 cycles of 95° C. for 15 seconds and 60° C. for one minute, and lastly, a dissociation stage of 95° C. for 15 seconds, 60° C. for 30 seconds, and 95° C. for 15 seconds. 
     Cells were harvested at chlorophyll (Chl) concentrations between 16-20 μg Chl/ml by centrifuging 10 ml of algal culture at 2000 g for one minute. The pellet was resuspended in lysis buffer (6% SDS in 1×PBS) plus 10 μl/ml each of Phosphatase Inhibitor Cocktails 1 &amp; 2 (Sigma: #P5726-5 ml and #P2850-5 ml). The samples were vortexed, mixed with loading buffer, and separated by SDS-PAGE under reducing conditions (1% β-mercaptoethanol) using an 8-16% Tris-Glycine gel (Invitrogen: #EC6048BOX). 
     Since SDS interferes with the Bradford total protein assay, the chlorophyll concentration of each of the samples was determined by the method of Harris (76) and equal amounts of chlorophyll (and thereby of protein) were loaded onto the gel. The gel was subjected to 125 Volts for approximately 2 hours, or until the dye front reached the bottom of the gel. The proteins were transferred to a PVDF Hybond-P membrane (Amersham: #RPN303F) for 1.5 hours at 25 Volts. Proof of equal loading of the gel lanes was obtained by staining the gel (for total protein) overnight in Coomasie (Pierce: #24590) (not shown). The stained gel was dried and sealed in cellophane (Invitrogen: #NC2380). 
     The membrane was blocked in 1×PBS containing 0.5% Tween-20, 3% BSA, and 1/1000 Biotin Blocking Buffer (IBA GmbH: #2-0501-002) following the IBA protocol for all washes. The membrane was probed with an anti strepII tag mouse antibody (IBA GmbH: #2-1507-001) at 4° C. overnight at a concentration of 1000 ng/ml (1/200 dilution of the stock concentration). The secondary antibody (donkey anti-mouse conjugated to horse radish peroxidase (HRP)) (Jackson Labs: #715-035-151) was allowed to incubate for one hour at room temperature at a concentration of 67 ng/ml (1:3000 dilution of the stock). Protein was detected using ECL-plus developing solutions (Amersham: #RPN2132) and visualized using a Kodak DS Image Station 440CF using 1D Image Analysis Software. 
     Bacterial hydrogenases were cloned and the best two were used as the parent genes in a combinatorial shuffle that resulted in a library of chimeric hydrogenases. Several chimeric hydrogenases showed hydrogen production, including three chimeras that produced hydrogen at levels three to four times higher than the parent hydrogenases. 
     Since several chimeras produced hydrogen, despite significant differences in their amino acid sequences, several mutational pathways may result in mutants with improved hydrogen production. The remainder of the mutant library can be tested for hydrogen production. The chimeric proteins with an improved hydrogen production can be sequenced and computationally modeled, both of which are capable of identifying the mutations that result in the improved chimeric proteins. Such modeling has been initiated and is described below. When such a comparison is made of the diversity of successful chimeras, it is possible that a combination of the individual mutations will result in a chimera with a dramatically increased level of hydrogen production. Shuffles with additional parent genes or successive shuffles of the improved chimeric hydrogenases from the first shuffle may also result in a highly evolved hydrogenase. 
     Example 4 
     Expression of Bacterial Hydrogenases and Measurement of Hydrogen Production 
     Five clostridial hydrogenases were cloned:  Clostridium saccharobutylicum  ( C. saccharobutylicum ),  Clostridium pasteurinum  ( C. pasteurinum ),  Clostridium perfringens  ( C. perfringens ),  Clostridium tetani  ( C. tetani ), and  Clostridium thermocellum  ( C. thermocellum ). All hydrogenases were cloned via restriction sites (NcoI/AscI) into the pET DLS expression vector as described further below (see also  FIG. 16 ) and subsequently tested for hydrogen production using the Methyl Viologen assay. The Methyl Viologen assay was used because bacteria need an artificial electron donor in place of the photosynthetic machinery that provides electrons to algal hydrogenases. 
     The original Clostridial hydrogenase,  C. acetobutylicum , along with its accessory proteins, were cloned by King et al. into two dual expression vectors, pET Duet and pCDF (Novagen: #TB337 and #TB390). These two vectors were used as the positive control. The Strep II tag “SAWSHPQFEK” (IBA GmbH: Göttingen, Germany) and a stop codon were cloned between the AscI and NotI restriction sites of the pET Duet vector (see  FIG. 16 ). Thus, the Strep II tag is in frame and immediately downstream of the AscI restriction enzyme site which is just downstream of the hydrogenase gene. This plasmid, renamed “pET DLS”, was designed in this manner so that both the hydrogenase and the tag were easy to insert and readily removable from the plasmid backbone. The pCDF plasmid was used without alterations. 
     Genomic DNA for  C. perfringens  was obtained from the American Type Culture Collection (ATCC: #13124D-5). The hydrogenase in  C. perfringens  (CpeI) was isolated by amplifying the DNA by PCR from aliquots of the genomic DNA. 
     The  C. thermocellum  (Cthl) hydrogenase, encoded in a plasmid, was obtained from the Lee Lynd Laboratory (Dartmouth College, NH), but was found to be missing a portion of sequence on the 3′ end. The missing sequence corresponded to an additional two amino acids (lysine and a stop codon), which were identified by comparison to the genomic information at the Joint Genome Institute (JGI: genome.jgi-psf.org/draft_microbes/cloth/cloth.home.html). The missing DNA sequence was added to the end of the gene using PCR and a new downstream primer with the additional bases included (see  FIG. 19 , SEQ ID NOS: 17-22). 
     The upstream and downstream PCR primers (see  FIG. 19 ) were designed to be homologous to the start and stop codons of the published sequences (NCBI:  C. perfringens  #AB035092 and BAA95936;  C. thermocellum #AAD33071 and AF148212). Both primers had extra bases (tails) that corresponded to either a NcoI restriction site (upstream primer (5′ CCAGCTGCTGCCAGAATTC 3′, SEQ ID NO. 13)) or an AscI restriction site (downstream primer (5′ CCAGCTGCTGCCAGAATTC 3′, SEQ ID NO. 14)). The  C. perfringens  PCR product was digested with NcoI and AscI (New England Biolabs (NEB): #R0193S and #R0558S) and cloned into the pET DLS vector (see  FIG. 16 ). Prior to the cloning of the  C. thermocellum  PCR product into pET DLS, it was blunt cloned into a shuttle vector (Lucigen&#39;s Clonesmart HC Kan: #40728-1) and mutated (GCC (Ala)→GCG (Ala) at 1410 bp) (Stratagene Quikchange II site directed mutagenesis kit: #200523) in order to remove an internal NcoI site (see  FIG. 19 ). Post mutation, the  C. thermocellum  PCR product was digested, as above, and cloned into the pET DLS vector. 
     BL21  E. coli  cells (Lucigen: #60300-1) were co-transformed with 30 ng of pCDF and 30 ng of pET DLS ( FIG. 16 ) using an electroporator (Bio-Rad Micropulser: #165-2100) and 1 mm cuvettes (Bio-Rad: #165-2089). The pET DLS vector contains the hydrogenase to be tested and also codes for the  C. acetobutylicum  HydE accessory protein. The pCDF vector codes for the HydF and HydG accessory proteins, also from  C. acetobutylicum . The standard settings for electroporating  E. coli  were: 1800 Volts, 600 Ohms, and 10 μF. After electroporation, the cultures were allowed to grow for one hour and plated (200 μL) on standard LB plates containing Ampicillin (50 μg/ml) and Streptomycin (50 μg/ml). 
     Expression of the hydrogenases was performed using the method of Posewitz et al. (M. C. Posewitz, M. C. et al., Identification of genes required for hydrogenase activity in  chlamydomonas reinhardtii , Biochem. Soc. Trans. 33: 102□4 (2005)). Overnight cultures (10 ml) were started from a single colony and used to inoculate sterile flasks containing 125 ml of LB media supplemented with Ampicillin and Streptomycin. The flasks were shaken (300 rpm) at 37° C. until they reached an optical density of A600=0.6-0.8 (approx. 5-6 hours). The optical density was measured on a spectrophotometer (Nanodrop Technologies: #ND-1000, Software Version 3.1.0). 
     In order to induce expression, IPTG was added to a final concentration of 100 mM and the flasks were transferred to a room temperature shaker (150 rpm) for one hour. The solutions were then transferred to 125 ml serum bottles (Wheaton: #223748), rubber septa (Chemglass: #CG-3022-24) were inserted, and the cultures were thoroughly degassed with Argon 5.0 (General Air: #UN 1006)) for 4-6 hours. The cultures were allowed to sit overnight at room temperature before being analyzed further. 
     The genome of  C. thermocellum  has been sequenced, so its accessory proteins (HydE, HydF, and HydG) were found by electronically searching the published genome (maple.1sd.oml.gov/microbial/cthe/) for comparable sequences to the known accessory protein sequences of  C. acetobutylicum ).  C. thermocellum  HydE, HydF, and HydG were designated as gene 2382 on contig 245, gene 3099 on contig 253, and gene 1213 on contig 221, respectively. 
     All three  C. thermocellum  accessory proteins were cloned into the same locations within the same vectors as their corresponding genes from  C. acetobutylicum . The new pET DLS ( FIG. 20 ) and pCDF plasmids ( FIG. 21 ), with the substituted  C. thermocellum  accessory proteins, were tested for hydrogen production using the Methyl Viologen assay described in Chapter Three. A second new pCDF plasmid was also created wherein the  C. thermocellum  HydG replaced the  C. acetobutylicum  version, but the  C. acetobutylicum  HydF was still intact. This hybrid vector, along with the first pCDF plasmid was constructed to determine if all three of the  C. thermocellum  accessory proteins are necessary for expression or if some combination of  C. thermocellum  and  C. acetobutylicum  proteins is needed. 
     Along with a positive and negative control, the above mentioned plasmids were tested for hydrogen expression. Five different combinations (see Table 3) of plasmids were tested and all were negative. These plasmids were also given to Matthew Posewitz (NREL: Golden, Colo.) and his lab obtained negative results as well. Hence, it seems that the hydrogenase from  C. thermocellum  requires a higher temperature or something yet to be determined, in order to produce molecular hydrogen. 
     
       
         
               
               
               
               
             
           
               
                   
                 TABLE 3 
               
               
                   
                   
               
               
                   
                 pDLS plasmid 
                 pCDF plasmid 
                 Result 
               
               
                   
                   
               
             
             
               
                   
                 Positive Control 
                 Standard pCDF 
                 Positive H2 production 
               
               
                   
                 Standard pDLS 
                 Standard pCDF 
                 No H 2  production 
               
               
                   
                 Plasmid #1 
                 Standard pCDF 
                 No H 2  production 
               
               
                   
                 Plasmid #1 
                 Plasmid #2 
                 No H 2  production 
               
               
                   
                 Plasmid #1 
                 Plasmid #3 
                 No H 2  production 
               
               
                   
                 Standard pDLS 
                 Plasmid #2 
                 No H 2  production 
               
               
                   
                 Standard pDLS 
                 Plasmid #3 
                 No H 2  production 
               
               
                   
                   
               
             
          
         
       
     
     Using primers that were designed for these putative sequences ( FIG. 22 ) and genomic DNA from  C. thermocellum  (American Type Culture Collection: #27405D), the genes for all three accessory proteins were obtained by PCR. KOD polymerase (VWR: #80511-384) was utilized with the following PCR conditions: 2′ initial denaturation at 94° C., 35 cycles of [15 sec. at 94° C., 30 sec at 55-60° C. gradient, and 1.5° at 72° C.], and a 10′ final extension at 72.degree. C 
     The  C. thermocellum  HydE gene was cloned into the  C. thermocellum  pET DLS vector, replacing the  C. acetobutylicum  HydE gene ( FIG. 20 ). However, the first step was to silently mutate the BglII site found within the  C. thermocellum  hydrogenase, so that the  C. acetobutylicum  HydE protein could be ligated only into the NdeI/BglII site ( FIG. 22 ). The Quikchange Multi site directed mutagenesis kit (Stratagene: #200515-5) was employed with the following PCR conditions: 1′ initial denaturation at 95° C., 30 cycles of [1′ at 95° C., 1′ at 55° C., and 16.5° at 65° C.]. Since a silent mutation was made in the  C. thermocellum  hydrogenase, it needed to be re-inserted in the NcoI/AscI site, once the  C. thermocellum  HydE gene was successfully ligated. The resulting plasmid was sequenced and compared to the expected sequence in the published genome. 
     The  C. thermocellum  HydF and HydG genes were cloned into pCDF, also in place of the corresponding  C. acetobutylicum  genes ( FIG. 21 ). As was done for the pET DLS plasmid, a silent mutation was introduced into the  C. acetobutylicum  HydF to remove an internal NdeI site ( FIG. 22 ). Again, the Quikchange Multi site directed mutagenesis kit (Stratagene: #200515-5) was utilized under similar PCR conditions: 1′ initial denaturation at 95° C., 30 cycles of [1′ at 95° C., 1′ at 55° C., and 13′ at 65° C.]. This mutation allowed us to replace the  C. acetobutylicum  HydG gene with the  C. thermocellum  version into the NdeI/BglII site without disrupting the HydF gene. The new  C. thermocellum  HydF was then cloned into the NcoI/BamHI sites, thereby replacing the mutated  C. acetobutylicum  HydF gene. The pCDF plasmid with the  C. acetobutylicum  HydF gene and the  C. thermocellum  HydG gene was kept, so it could be determined if all or just one of the  C. thermocellum  accessory proteins were necessary for successful expression. 
     Two sets of primers were used to create each of the silent mutations in the  C. thermocellum  pET DLS and pCDF vectors. The second set of complementary primers was non-mutational and hybridized, approximately half way around the plasmid from the engineered point mutation. The second set of primers was designed to eliminate any problems due to the processivity of the polymerase, so that it could successfully copy these large plasmids ( FIG. 22 ). 
     Methyl Viologen Assay: 
     Using an argon purged gas-tight syringe (Hamilton: #81430), 1 ml of 2× Methyl Viologen solution and 100 μL of a 100 mM dithionite solution were added through a septum into an anaerobic 10 ml serum bottle (Wheaton: #223739). Subsequently, 1 ml of cells was transferred anaerobically to the serum bottle. The serum bottles were allowed to sit for at least 60 minutes, mixing occasionally by inversion. 
     The headspace of the serum bottles was tested for molecular hydrogen by gas chromatography. A 400 μL sample of the headspace gas was injected into a GC(SRI Instruments: #310C) containing a molecular sieve column (SRI Instruments: #5A 2 feet long) at a constant temperature of 40° C. The results were recorded using the PeakSimple Software (SRI Instruments version 3.29). 
     Three hydrogenases, those isolated from  C. saccharobutylicum, C. pasteurinum , and  C. perfringens , showed hydrogen production that was comparable to the values published for  C. acetobutylicum . However, the hydrogenase isolated from  C. tetani  showed attenuated hydrogen production and the  C. thermocellum  hydrogenase showed no hydrogen production whatsoever. Since  C. thermocellum  is a thermophilic organism (55° C.) and it has the lowest homology of the five Clostridial hydrogenases that were cloned, it is possible that this hydrogenase needs its own accessory proteins in order to fold correctly. 
     Example 5 
     Cloning and Expression of Chimeric Hydrogenases and Measurement of Hydrogen Production 
     Two chimeric libraries, by two separate methods, were created in an attempt to realize an improvement in the hydrogen production and/or oxygen tolerance over the parental Clostridial hydrogenases. First, a random mutagenesis library was created from C. saccharobutylicum. Utilizing the two wild-type hydrogenases that produced the most hydrogen, the hydrogenases from C. saccharobutylicum and  C. acetobutylicum , a second library was created utilizing the variation on the DOGS technique as described herein above (see also  FIGS. 13 through 15 ). 
     Cloning and expression of the resulting chimeras was achieved in  E. coli  substantially as described above in Examples 3 and 4, except that the random mutagenesis library was cloned via restriction sites (NcoI/AscI) into the pET DLS vector in place of the parent hydrogenase. The sequenced clones had an average of six mutations per gene (1.8 kb). A predominance of Adenine (A) &amp; Thymine (T) transitions over Guanine (G) &amp; Cytosine (C) transitions was observed, but was in accordance with the relative percentage of As and Ts to Gs and Cs. None of the mutations were identical. This library was not tested for hydrogen production and would benefit from a high through-put screen. 
     The library created by the DOGS method, using the SapI variation, generated 254 distinct chimeras. Sixteen chimeras, which were representative of the possible mutations in the second half of the protein, were randomly chosen and tested for hydrogen production (see  FIG. 17 ). Several showed normalized hydrogen production at or greater than the positive control (see  FIG. 18 ). The positive control was the non-mutated hydrogenase from  C. acetobutylicum  used in an earlier study (King et al., Functional Studies of [FeFe] Hydrogenase Maturation in an  Escherichia Coli  Biosynthetic System. J. Bacteriol. 188:2163-72 (2006)). Specifically, chimera #2 (protein sequence—SEQ ID NO: 1 and DNA sequence—SEQ ID NO: 2) shows the most hydrogen production (4× the hydrogen production of the positive control). In addition to producing more hydrogen, chimera #2 also produced molecular hydrogen at a faster rate than the other mutants; the rate for chimera #2 dramatically decreased after one hour of incubation. The other chimeras produced 3.0-50% more hydrogen in the second hour of incubation, which likely indicates that they have a slower turnover of the hydrogenase protein. The remaining hydrogenase chimeras, #18A (protein sequence—SEQ ID NO: 85 and DNA sequence—SEQ ID NO: 153), 18S (protein sequence—SEQ ID NO: 86 and DNA sequence—SEQ ID NO: 154), and 28A (protein sequence—SEQ ID NO: 89 and DNA sequence—SEQ ID NO: 157), produced 3×, 2.8×, and 1× the amount of hydrogen of the positive control, respectively. Of the mutants tested, chimera #28S (protein sequence—SEQ ID NO: 158 and DNA sequence—SEQ ID NO: 90) was the only hydrogenase that exhibited attenuated hydrogen production, approximately 10% of the positive control. 
     Thus the approach described herein produced mutant hydrogenases capable of producing hydrogen at a faster rate than the wild-type control, the hydrogenase from  C. acetobutylicum. Chimeras such as # 2, 18A, and 18S showed more hydrogen production than the positive control. In fact, after one hour of incubation in the Methyl Viologen assay, these chimeras show 1-4× the amount of hydrogen production relative to the positive control. Since nearly all of the chimeras tested produced hydrogen above the level of the positive control (Chimeras #2, 18A, 18S, and 28A), or else produced no hydrogen at all (Chimeras #1, (protein sequence—SEQ ID NO: 83 and DNA sequence—SEQ ID NO: 151), 21A (protein sequence—SEQ ID NO: 87 and DNA sequence—SEQ ID NO: 155), 21S (protein sequence—SEQ ID NO: 88 and DNA sequence—SEQ ID NO: 156), 30S (protein sequence—SEQ ID NO: 91 and DNA sequence—SEQ ID NO: 159), and 44 (protein sequence—SEQ ID NO: 92 and DNA sequence—SEQ ID NO: 160)), which indicates that protein folding is critical and impacts hydrogen production in a rather Boolean fashion. 
     Several structure/function relationships can be noted of the chimeras tested thus far. For example, chimera #1 and #18A are identical except for segment #7, and chimeras #18A and #18S are identical to chimeras #21A and #21S, except for segment #8 (see  FIG. 17 ). Since chimeras #18A and #18S are positive for hydrogen production, whereas chimeras #1, #21A and #21S show no hydrogen production, it appears that beneficial results occur when segment #7 is derived from  C. acetobutylicum  and segment #8 be derived from C. saccharobutylicum. However, if both segments #5 and #8 are derived from  C. acetobutylicum  (Chimeras #28A and #28S) instead of just segment #8 (Chimeras #21A and #21S), hydrogen production is restored, although not to the levels of chimeras #18A and #18S. 
     Chimera #2 is a simple amalgam whose first half is derived from C. saccharobutylicum while its second half is derived from  C. acetobutylicum . Again, the simple change of one segment (segment #8) is enough to decrease hydrogen production to zero (Chimera #30S), but changing segments #5 and #6 results in the hydrogen production being restored (Chimera #18S). Chimera #2 produced the most hydrogen of all the mutants tested, while its mirror image (Chimera #1) showed no discernible hydrogen production. 
     Example 6 
     Expression of Hydrogenases in  C. reinhardtii  and Measurement of Hydrogen Production 
     Chimeric and wild-type hydrogenases were cloned into plasmid pSMP1c for expression in  C. reinhardtii . The pSMP vectors were transformed into dw15-1, a fast growing variant of the common  C. reinhardtii  cc425 strain, using a BTX ECM 630 electroporator with the following settings: 750 Volts, 25 μF, and 1575 ohms. Immediately after electroporation, the algae were mixed with 3 ml of cornstarch suspended in TAP+60 μM sucrose and plated on Tris-Acetate-Phosphate (TAP) plates containing 10 μg/ml of the antibiotic Zeocin. The resulting transformants were grown under oxygenic conditions in 50 ml of sterile TAP media utilizing 250 ml Erlenmeyer flasks. The flasks were shaken at 100-150 rpm and exposed to 100 μE (Einstein) of light from standard fluorescent fixtures. After one week, they were transferred to standard TAP plates containing no antibiotics. The transformants were transferred to plates without antibiotics as Zeocin is a presumed mutagen and antibiotics are not necessary to maintain the transformants, since the plasmid DNA is incorporated into the alga&#39;s genome. 
     For hydrogen production, algal transformants were inoculated from plates into liquid 10, culture (TAP, pH=7.2, or TAP plus arginine) and grown as detailed above. Cells were harvested at chlorophyll (Chl) concentrations between 10-18 μg Chl/ml (exponential growth phase), centrifuged at 4000 g for 5 minutes, and resuspended at 16 μg Chl/ml in anaerobic induction buffer (AIB) plus 0.001% (wt) PdC12 and 8 μM Resazurin. AIB consists of 40 mM K2HPO4, 10 mM KH2PO4, and 3 mM MgCl2. A stock suspension of 0.05% (wt) PdCl2 was made by mixing water and powdered PdCl2 (Pressure Chemical: #1735) in a serum bottle, which was capped, flushed with nitrogen gas (General Air: purity 4.8) for 15 min, and then autoclaved. A working solution of 0.01% (wt) Resazurin (Research Chemicals Ltd: #21187) in AIB was created by diluting a 0.1% (wt) stock solution. 
     Resazurin Assay: 
     The resazurin assay was developed to measure the evolution of molecular hydrogen from  C. reinhardtii  cells. Resazurin is a redox dye that reacts with molecular hydrogen ( FIG. 23 ). Since resazurin has two oxidized forms rather than just one, its reaction kinetics are quite complex. For example, molecular hydrogen reacts irreversibly with resazurin (blue, A max  of 600 nm) to form reduced resorufin (colorless). The oxidized form of resorufin (pink, A max  of 580 nm) is in a fast equilibrium with the reduced form. In theory, the reaction between molecular hydrogen and resazurin could be measured by the decrease in absorbance of resazurin or the increase in absorption of the oxidized form of resorufin. 
     Using sterile technique, 2.5 ml of the algal suspensions were added to a sterile anaerobic cuvette, made to order by Allen Scientific Glass (Boulder, Colo.) by fusing the top 2 cm of a serum bottle to about 4 cm of square glass tubing. Each cuvette contained a stir bar (VWR #58949-030) and was capped with a butyl rubber stopper (VWR #80062-438). The cuvettes were wrapped in aluminum foil to inhibit photosynthetic oxygen generation and purged with Argon (General Air: purity 5.0) for 15 minutes to remove the oxygen using 25G 7/8″ needles (VWR #BD305124). The purged cuvettes were incubated at room temperature, in the dark, for four hours in order to induce production of oxygen sensitive hydrogenases. 
     After induction, cuvettes were placed on a stir plate and exposed to blue light (Dolan-Jenner #BG2820) filtered by a solution of 1% CUSO4. The cuvettes were exposed for 6 minutes at 5001 μE. Following illumination, the cuvettes were immediately placed in 50 ml Falcon tubes and centrifuged for 5 min at 5000 g. Spectra (350-800 nm) of the cell-free supernatants were obtained using a Varian CARY 5E UV/Vis/NIR spectrophotometer. 
     Hydrogen gas in the headspace of each cuvette was quantified using a HP 5890 series II gas chromatograph equipped with a molecular sieve column (Supelco 5A 60/80) and a thermal conductivity (TCD) detector. A simple constant temperature (60° C.) program is sufficient to separate H 2  from O 2  and N 2 . 
     With no algae present, resazurin successfully reacted with 10% molecular hydrogen (Balance Nitrogen, General Air) in anaerobic cuvettes (Allen Scientific Glass: Boulder, Colo.). ( FIG. 24 ). Hydrogen gas (100 μl of 10% H2) was injected into cuvette #2 (the second cuvette from the left), which was enough hydrogen to turn the solution slightly pink. Cuvettes #3 and #4 both had an excess of hydrogen injected, which forced the equilibrium into the reduced and colorless form of resorufin. Afterwards, cuvette #4 was opened to the air and thereby completely converted into the oxidized form of resorufin. The resulting UV/Vis spectra were overlapping but distinguishable from each other ( FIG. 25 ). 
     Alternatively, if the alga transformants were exposed to oxygen before they were tested for molecular hydrogen production this assay could also be used to select for alga that contained mutant hydrogenases that are less sensitive to oxygen. If the alga were exposed to a concentration of oxygen at or slightly above their tolerance level, the hydrogenase chimeras that survived and produced hydrogen would have a higher oxygen tolerance. Finally, the enhanced hydrogenase chimeras from both of these selection techniques are sequenced and their individual mutations combined to create a mutant that exhibits both enhancements. 
     As indicated by the reduced size of the 600 nm peak in  FIG. 26 , a yet to be determined compound(s) in the algae is responsible for reducing the resazurin before it could be reduced by hydrogen generated from the algae. This reduced spectra occurred, regardless of when the resazurin/PdCl 2  mixture was added to the cells. Furthermore, this technique is not recommended because the spectrum of hydrogen producing algae (dw15-1,  FIG. 26A ) was not significantly different from the spectrum of non-hydrogen producing algae (sta6 mutant,  FIG. 26B ). A difference of approximately 0.07 absorbance units, at λ=610 nm, between these two strains was common. Duplicate samples produced spectra that were similar but whose absorbencies, at λ=610 nm, deviated slightly from each other. In addition to sta6, several other mutants with retarded hydrogen production were tested. Thus, the present disclosure provides that while spectra similar to sta6 were evident, there is no evident relationship between hydrogen production and the absorbance differences (A 610  of dw15-1 minus the A 610  of a mutant). However, molecular hydrogen definitely reacts with resazurin and the difference is significant. Yet, when algae generated the molecular hydrogen the resulting spectra showed a reaction with molecular hydrogen or a reaction with some unknown compound(s) within the algae. 
     Alternatively, selection of algal mutants can be accomplished using the metronidazole assay. The electrons created by the hydrolysis of water are ultimately transported to ferredoxin ( FIG. 27 ). Ferredoxin is a common compound that is involved in many electron transport pathways. In  C. reinhardtii , as in most plants, ferredoxin normally transfers electrons to the Calvin cycle, the pathway that is responsible for cell growth by converting carbon dioxide from the atmosphere into glucose. In  C. reinhardtii , ferredoxin can also transfer electrons to a hydrogenase when the algae experiences dark and anaerobic conditions. 
     The photo-hydrogen group at the National Renewable Energy Lab (NREL: Golden, Colo.) developed a chemical method to select for oxygen tolerant mutants. The  C. reinhardtii  cells were induced to produce hydrogenases and then treated with metronidazole (MNZ), a toxic substance that non-specifically kills cells by creating oxygen free radicals known as superoxide radicals ( FIG. 27 ). Superoxide radicals spawn organic radicals and hydrogen peroxide, all of which are quite toxic to  C. reinhardtii  cells. When the algal cells are induced to produce hydrogenases, cell growth ceases and a brief period of hydrogen production ensues when the cells are exposed to light. 
     If, however, sufficient concentrations of metronidazole and oxygen are added before the cells are exposed to the light, the hydrogenase is irreversibly inhibited by the oxygen, ferredoxin transfers all of its electrons to MNZ, and the algal cells die. It was theorized that, if little or no oxygen were added, a proportion of the electrons would be transferred to the surviving hydrogenases; hydrogenases that likely have a higher tolerance to oxygen. Given a constant number of available electrons, the algal cells containing hydrogenases with a higher tolerance to oxygen would then see more electrons transferred to the hydrogenase and fewer electrons transferred to MNZ; a lower concentration of toxic radicals and enhanced cell viability would result. The  C. reinhardtii  cells containing enhanced hydrogenases could then be isolated and analyzed. 
     Wild type strains of  Chlamydomonas reinhardtii , dw15-1 (Barb Sears, MSU) and cc425 were grown in tris-acetate-phosphate (TAP), pH=7.2, or TAP plus 50-200 μg/ml arginine (cc425) in sterile 250 ml Erlenmeyer flasks on a shaker rotating at 150 rpm in a 25° C. constant temperature room and constantly illuminated with 100 μE of white fluorescent light (1 μE=1 μEinstein=1 μmoles photons/m 2 sec=500 foot candles). Cells were harvested when they reached a concentration of 2.4×10 6  cells/ml, centrifuged at 2-3000 g for 5 minutes at 22° C., and resuspended to a final concentration of 1×10 7  cells/ml in sterile induction buffer (AIB) plus 20 μl of sterile 0.5M sodium acetate (10 mM final concentration) per ml of cells. AIB consists of 40 mM K 2 HPO 4 , 10 mM KH 2 PO 4 , and 3 mM MgCl 2  and was prepared according to Ghirardi et al. 
     Using sterile technique, 1 ml of the resuspended algae was transferred to a 12 ml serum vial with a flea sized stir bar and sealed with a sterile septum (VWR #80062-438) using sterilized forceps. The serum vials were wrapped in aluminum foil to inhibit photosynthetic oxygen generation and gently purged with Argon (General Air: purity 5.0) for 2 hours using 25G7/8 needles (VWR #BD305124). A manifold was used to flush multiple samples simultaneously. The purged serum vials were incubated at room temperature, in the dark, an additional 2 hours in order to induce production of oxygen sensitive hydrogenases. Metronidazole treatment was started immediately or the samples were refrigerated overnight to minimize the loss of hydrogenase activity. If refrigerated overnight, the vials were allowed to warm to room temperature, approximately 15-30 minutes, and flushed with argon for 15 minutes before the metronidazole selection. 
     During this procedure, the septa were frequently wiped with ethanol to insure sterility. A maximum of 45 ml of fresh Metronidazole stock solution (50 mM Metronidazole &amp; 1.25 mM NaN 3 ) was added to 70 ml serum vials wrapped in aluminum foil. The metronidazole containing vials were sealed with a septum and gently flushed with argon for 30 minutes, as described above. Add 4 ml of the anaerobic MNZ stock, using a 5 ml syringe and a 20 gauge needle. The MNZ-cell mixture was placed on a stir plate until the solution was well mixed. After mixing, an initial 100 μl sample of cells was removed using a 1 ml syringe and a 23 gauge needle in a dark sterile hood in a very dark lab. The cells were dispensed into an eppendorf tube and stored at room temperature in the dark. An aliquot of 100% oxygen (2.2 ml of 100% oxygen results in an approximate final headspace oxygen concentration of 25% given 9 ml of headspace in a 15 ml serum vial) was added to the headspace of the MNZ-cell mixture. The vial was shaken vigorously, by hand, for 5-10 seconds and then allowed to mix vigorously on a stir plate for 4 minutes. Afterwards, the aluminum foil was removed and the serum vial was exposed to 320 μE (approx. 2200 ft. candles) of light using a source (Dolan-Jenner #BG2820) filtered by a solution of 1% CuSO 4 . The cells were exposed to the light for 6 minutes and then a second 100 μl sample of cells was removed to a second Eppendorf tube and store, as above. The serum vials were re-wrapped in aluminum foil and also stored in the dark at room temperature. 
     In a dark laminar flow hood, the time point samples were washed and resuspended in TAP before making the following dilutions: 1/100, 1/1000 and 1/10,000. For each of the dilutions, 100 μl was plated using 1-2 ml of cornstarch solution (25% cornstarch in TAP+60 mM sucrose). In order to obtain all of the surviving algae, the serum vials were washed several times in TAP and inoculated in 50 ml of TAP and grown, as above. 
     MNZ is quite toxic, so gloves were always used when it was handled. The excess MNZ solutions, as well as the MNZ/cell mixtures, were disposed into a waste bottle and submitted to the Environmental Health and Safety department. 
     However, the metronidazole failed to kill more than 90% of the wild type  C. reinhardtii  cells, regardless of the amount of oxygen that was injected into the serum vials. In fact, the cell death rate did not increase even when the cells were bubbled with 100% oxygen. Therefore, the technique was abandoned. Several brands of MNZ were tried with identical outcomes, so possibilities for why this technique didn&#39;t work include that the MNZ was unable to enter the algal cell or that the MNZ entered the cytoplasm, only to be destroyed. Since very little remained of the original NREL MNZ technique and the originally published algal selected mutants contained hydrogenases that were only slightly enhanced over the wild-type, it remains unknown as whether this technique was responsible for the published mutants. 
     More recently, hydrogen has been used in the “hydro-treating” or purification of fuel oil, the hydrogenation of oils to form margarine, in the production of ammonia-based fertilizers, in rocket fuel, and in fuel cells. The chemical reactions forming the basis of the fuel cell were first realized by William Nicholson, Anthony Carlisle, Christian Schoenbein, and William Grove; Swiss and English scientists in the early 1800s who were the first to combine oxygen and hydrogen gases in order to produce water. The term “fuel cell” was coined in 1889 when the very first fuel cell was created by Ludwig Mond and Charles Langer. Francis Bacon invented the first practical fuel cell in 1959. Since his invention, fuel cells have been used to power a diverse range of vehicles from a simple farm tractor to the Apollo mission and the space shuttles. It wasn&#39;t until the year 2000 that practical fuel cells for cars were unveiled by Ballard Power Systems. Perhaps someday soon, John Bockris&#39; 1970 prediction of a “hydrogen economy” will become a reality and a national network of hydrogen energy will become the norm. In order to have a national network, reliable method of producing molecular hydrogen is needed. Hydrogen is currently produced by steam reforming the hydrogen atoms from coal or natural gas. However, the present disclosure providing for biohydrogen production from photosynthetic algae has the potential to be a viable alternative to hydrogen production from fossil fuels. It would not produce greenhouse gases (H 2 O+sunlight→O 2 +H 2 ); in fact, algae like most plants, utilizes carbon dioxide for cellular growth, so it would serve as a carbon sink. In addition, a bioreactor would not produce toxic waste, just algae and wastewater; similar to a fish tank. Also, a bioreactor would likely be about the size of an air conditioner and survive on low amounts of sunlight, so it would occupy a small amount of space and it could be located anywhere in the country. 
     Hydrogen could be generated locally by on-site electrolyzers or more likely by a centralized power plant. If molecular hydrogen is generated in a central location, it could eventually be piped to on site compressors and storage tanks or in the near-term, it could be simply delivered in trucks, just as gasoline is today. Often hydrogen is safer than gasoline in that it is non-toxic and non-poisonous. In addition and unlike gasoline, hydrogen will not contribute to groundwater pollution when it leaks from an underground storage tank. 
     Molecular hydrogen could be used to power cars directly as well as truck and trains, personal residences, and the workplace. In fact, the first area of significant usage of hydrogen power is via stationary fuel cells in industrial applications. Stationary and modular fuel cells already provide the benefit of highly reliable power with a consistent voltage, which is ideal for modern industries, which depend on computers. Fuel cells are also available for portable uses such as remote construction/military sites, laptop computers, and cell phones. To date, fuel cells have also been installed and demonstrated in 50 different types of mobile vehicles from “buses to bicycles.” Lastly, home usage of fuel cells is possible as modular fuel cells, about the size of a common refrigerator, can provide enough power for a residence. 
     There are a myriad of different types of fuel cells, however the proton exchange membrane (PEM) fuel cell is the most common ( FIG. 29 ) and hydrogen is the ultimate fuel for a fuel cell. In principle, fuel cells operate by chemically combining hydrogen with oxygen to form water, electricity, and heat. In fact, a fuel cell operates similarly to a battery that doesn&#39;t require recharging. Specifically, hydrogen gas flows into the anode side of a PEM fuel cell where a platinum catalyst removes the protons of hydrogen from the electrons. The electrons form a current, i.e. electricity, while the protons pass through the membrane that divides the anode from the cathode. Meanwhile, oxygen or air enters the cathode where it is combined with the protons and the electrons in order to form water ( FIG. 30 ). Hydrogen fuel cells perform the reverse electrolysis reaction. 
     Example 7 
     Development of a Computational Predictive Model 
     Six chimeric hydrogenases and a wild type hydrogenase from  C. acetobutylicum  were used to develop a prediction model for other hydrogenase mutants that were not tested for hydrogen production. The hydrogenase mutants yielded hydrogen productions from zero to 4 times that produced by the wild type control. 
     In the first step, amino acid sequences were used to construct alpha-helix structures. These structures were then energy minimized using OPLS molecular mechanics technology. Next, the positive and negative electrostatic potential surface areas (EPSA) for each structure were calculated using a probe radius of 1.4 Å. The computations were performed with the HyperChem 7.5 computational chemistry program using a 32 bit PC computer. 
     In the next step, positive to negative EPSA ratios were calculated (overlapping positive and negative regions were not included in the calculations) for each of the mutant and wild-type hydrogenases. The positive to negative EPSA ratios were plotted versus experimentally measured hydrogen production for the six hydrogenase mutants and the wild type control. The plot was regressed with several mathematical functions using the Microcal Origin 4.1 program. The function that fit the plot with the least error was the Log Normal Gaussian. Hence, this function was chosen as the prediction model. The Log Normal Gaussian function has been successfully used in several other studies relating molecular performance with molecular structure. The plot is shown in  FIG. 32 . 
     The selected prediction model suggests that the ratio of positive to negative EPSA is a measure of how a given protein mutant will fold into its tertiary state and how the final tertiary state affects hydrogen production. These results show that an EPSA ratio in the range of about 15 to about 115 is indicative of the capability to produce hydrogen at a greater rate than that observed with the wild type control (EPSA ratio of 116). The maximum hydrogen production was correlated with a positive to negative EPSA ratio of 42 (see  FIG. 32 ). 
     In the final step of testing the prediction model, fifteen hydrogenase mutants, for which hydrogen productions were not experimentally obtained, were analyzed using the above described procedure. Five mutants exhibited positive to negative EPSA ratios that were within the above-mentioned range. The mutants were named C26S (protein sequence—SEQ ID NO: 207 and DNA sequence—SEQ ID NO: 144), C4S (protein sequence—SEQ ID NO: 187 and DNA sequence—SEQ ID NO: 122), C22S (protein sequence—SEQ ID NO: 203 and DNA sequence—SEQ ID NO: 140), C24S (protein sequence—SEQ ID NO: 205 and DNA sequence—SEQ ID NO: 142), and C25S (protein sequence—SEQ ID NO: 206 and DNA sequence—SEQ ID NO: 143), and each of these are predicted to produce hydrogen at 39.7 times, 32.6 times, 31.8 times, 27.0 times, and 26.8 times, respectively, of that produced by the wild type protein. In contrast, the remaining ten mutants, C6S, (protein sequence—SEQ ID NO: 192 and DNA sequence—SEQ ID NO: 127), C11S (protein sequence—SEQ ID NO: 194 and DNA sequence—SEQ ID NO: 129), C13S (protein sequence—SEQ ID NO: 196 and DNA sequence—SEQ ID NO: 131), C14S (protein sequence—SEQ ID NO: 197 and DNA sequence—SEQ ID NO: 132), C15S (protein sequence—SEQ ID NO: 198 and DNA sequence—SEQ ID NO: 133), C19S (protein sequence—SEQ ID NO: 201 and DNA sequence—SEQ ID NO: 137), C23S (protein sequence—SEQ ID NO: 204 and DNA sequence—SEQ ID NO: 141) and C29S (protein sequence—SEQ ID NO: 209 and DNA sequence—SEQ ID NO: 147) are predicted to yield no significant hydrogen productions since their EPSA ratios were outside the EPSA range of about 15 to about 115. These predictions have proven correct in subsequent tests. 
     Example 8 
     Hydrogen Production by Chimeric Hydrogenases 
     Eleven chimeric bacterial hydrogenases were cloned and expressed in BL21  E. coli  cells essentially as detailed above. Hydrogen production was determined using the Methyl Viologen assay essentially as described above. A wild-type hydrogenase from  Clostridium acetobutylicum  was used as a positive control. The levels of hydrogen produced by the chimeric hydrogenases ranged from about 0 to about 399 times that produced by the wild-type hydrogenase. ( FIG. 33 ) 
     Twenty-one chimeric algal hydrogenases were cloned and expressed in BL21  E. coli  cells essentially as detailed above. Hydrogen production was determined using the Methyl Viologen assay essentially as described above. The amount of hydrogen produced by the chimeric algal hydrogenases ranges from about 0 to about 2.5 times that produced by the wild-type hydrogenase. ( FIG. 34 ) 
     Example 9 
     Computational Modeling and Calculation of EPSA Ratios 
     EPSA ratios of all the bacterial and algal chimeric hydrogenases tested above as well as wild-type hydrogenases were calculated essentially as described above in Example 7 except a probe radius of 1.0 Å was used and the overlapping positive and negative EPSA regions were included in calculating the EPSA ratio. Positive to negative EPSA ratios were plotted against hydrogen production for the bacterial and algal mutant hydrogenases. The plots were regressed with several mathematical functions using the Microcal Origin 4.1 program. The Log Normal function was found to yield the least error for both bacteria ( FIG. 33 ) and algal mutants ( FIG. 34 ). Thus, the Log Normal function was chosen to define the range of mutants that yielded the optimum hydrogen production and for use in finding other optimum mutants. 
     Bacterial mutants with a positive to negative EPSA ratio from about 2 to about 50 yield hydrogen productions greater than that obtained with the wild type control. That is, hydrogen productions of 2.8, 4.0, 3.6 and 1.1 times that for the wild-type control were obtained at EPSA ratios of 5.2, 9.1, 25.4 and 38.3, respectively. The best mutants had EPSA ratios of 15.2 and 18.5, which yielded hydrogen productions of 399 and 13.5 times that of the wild-type control ( FIG. 33 ). The model predicts that an EPSA ratio of about 16 will produce hydrogen of about 530 times that of the wild type. The wild type control has an EPSA ratio of 54.5. 
     The Log Normal model for the algal mutants ( FIG. 34 ) shows that a positive to negative EPSA ratio from about 1 to about 4 yielded hydrogen productions equal to or greater than that produced by the algal wild-type control. The EPSA ratio of the wild-type algal hydrogenase was about 1.6. The maximum algal hydrogen production of 3.4 times that of the wild type occurred with an EPSA ratio of about 2.2. Other results suggest that algal mutants with positive to negative EPSA ratios in range of about 12 to about 22 will yield hydrogen productions similar to those via bacteria mutants, i.e., about 400 to about 530 times that of the wild type. Hence, there appear to be two EPSA ratio ranges in algal hydrogenases where hydrogen production is optimized. 
     Example 10 
     Calculated Positive/Negative EPSA ratios and Corresponding Hydrogen Production 
     The positive to negative EPSA ratios for several additional bacterial and algal hydrogenase sequences were calculated as described herein above. Those which had a positive to negative EPSA ratio in the preferred range were also found to demonstrate high levels of hydrogen production relative to wild type controls. The bacterial sequences (both DNA and protein), corresponding EPSA ratios and experimentally determined rates of hydrogen production for the bacterial hydrogenases are listed below. In addition, three algal sequences are provided below. These sequences were also experimentally analyzed for EPSA ratios and hydrogen production rate, as described herein, but were analyzed without their leader sequences, although the sequence for the leader is provided as well as the actual sequences tested. 
     
       
         
               
               
             
           
               
                 Bacterial Hydrogenase Mutant C16A: 
                   
               
               
                 Positive/Negative EPSA Ratio: 15.2 
               
               
                 Hydrogen Production: 399 times Wild Type 
               
               
                 Amino Acid Sequence (SEQ ID NO: 224): 
               
               
                 MAKTIILNGNEVHTDKDITILELARENNVDIPTLCFLKDCGNFGKCGVCMVEVEG 
               
               
                   
               
               
                 KGFRAACVAKVEDGMVINTESDEVKERIKKRVSMLLDKHEFKCGQCSRRENCEFLKLVI 
               
               
                   
               
               
                 KTKAKASKPFLPEDKDALVDNRSKAIVIDRSKCVLCGRCVAACKQHTSTCSIQFIKKDG 
               
               
                   
               
               
                 QRAVGTVDDVCLDDSTCLLCGQCVIACPVAALKEKSHIEKVQEALNDPKKHVIVAMAP 
               
               
                   
               
               
                 SIRTSMGELFKLGYGVDVTGKLYASMRALGFDKVFDINFGADMTIMEEATEFIERVKNN 
               
               
                   
               
               
                 GPFPMFTSCCPAWVRLAQNYHPELLDNLSSAKSPQQIFGTASKTYYPSISGIAPEDVYTV 
               
               
                   
               
               
                 TIMPCNDKKYEADIPFMETNSLRDIDASLTTRELAKMIKDAKIKFADLEDGEVDPAMGT 
               
               
                   
               
               
                 YSGAGAIFGATGGVMEAAIRSAKDFAENKELENVDYTEVRGFKGIKEAEVEIAGNKLNV 
               
               
                   
               
               
                 AVINGASNFFEFMKSGKMNEKQYHFIEVMACPGGCINGGGQPHVNALDRENVDYRKLR 
               
               
                   
               
               
                 ASVLYNQDKNVLSKRKSHDNPAIIKMYDSYFGKPGEGLAHKLLHVKYTKDKNVSKHE 
               
               
                   
               
               
                 DNA Sequence (SEQ ID NO: 225): 
               
               
                 atggccaaaacaataatcttaaatggcaatgaagtgcatacagataaagatattactatccttgagctagcaagagaaaataatgta 
               
               
                   
               
               
                 gatatcccaacactctgctttttaaaggattgtggcaattttggaaaatgcggagtctgtatggtagaggtagaaggcaagggctttagagctg 
               
               
                   
               
               
                 cttgtgttgccaaagttgaagatggaatggtaataaacacagaatccgatgaagtaaaagaacgaatcaaaaaaagagtttcaatgcttcttg 
               
               
                   
               
               
                 ataagcatgaatttaaatgtggacaatgttctagaagagaaaattgtgaattccttaaacttgtaataaagacaaaagcaaaagcttcaaaacc 
               
               
                   
               
               
                 atttttaccagaagataaggatgctctagttgataatagaagtaaggctattgtaattgacagatcaaaatgtgtactatgcggtagatgcgtag 
               
               
                   
               
               
                 ctgcatgtaaacagcacacaagcacttgctcaattcaatttattaaaaaagatggacaaagggctgttggaactgttgatgatgtttgtcttgatg 
               
               
                   
               
               
                 actcaacatgcttattatgcggtcagtgtgtaatcgcttgtcctgttgctgctttaaaagaaaaatcccatatagaaaaagttcaagaagctctta 
               
               
                   
               
               
                 atgaccctaaaaaacatgtaatagttgctatggcaccatcaatcagaacttctatgggagagttatttaaattaggctatggggttgatgtaactg 
               
               
                   
               
               
                 gaaaattatatgcttcaatgagagcattaggatttgataaggtatttgatattaactttggggctgatatgacaataatggaagaagcaacagag 
               
               
                   
               
               
                 tttattgaaagagttaaaaataatggcccattccctatgtttacatcttgctgtcctgcatgggtaagattagctcaaaattatcatcctgaatta 
               
               
                   
               
               
                 ttagataatctttcatcagcaaaatcaccacaacaaatatttggtactgcatcaaaaacttactatccttcaatttcaggaatagctccagaagat 
               
               
                   
               
               
                 gtttatacagttactatcatgccttgtaatgataaaaaatatgaagcagatattcctttcatggaaactaacagcttaagagatattgatgcatcc 
               
               
                   
               
               
                 ttaactacaagagagcttgcaaaaatgattaaagatgcaaaaattaaatttgcagatcttgaagatggtgaagttgatcctgctatgggtacttac 
               
               
                   
               
               
                 agtggtgctggagctatctttggtgcaaccggtggcgttatggaagctgcaataagatcagctaaagactttgctgaaaataaagaacttgaaaat 
               
               
                   
               
               
                 gttgattacactgaagtaagaggctttaaaggcataaaagaagcggaagttgaaattgctggaaataaactaaacgttgctgttataaatggtgct 
               
               
                   
               
               
                 tctaacttcttcgagtttatgaaatctggaaaaatgaacgaaaaacaatatcactttatagaagtaatggcttgccctggtggatgtataaatggt 
               
               
                   
               
               
                 ggaggtcaacctcacgtaaatgctcttgatagagaaaatgttgattacagaaaactaagagcatcagtattatacaaccaagataaaaatgttc 
               
               
                   
               
               
                 tttcaaagagaaagtcacatgataatccagctattattaaaatgtatgatagctactttggaaaaccaggtgaaggacttgctcacaaattacta 
               
               
                   
               
               
                 cacgtaaaatacacaaaagataaaaatgtttcaaaacatgaa 
               
               
                   
               
               
                 Bacterial Hydrogenase Mutant C15S: 
               
               
                 Positive/Negative EPSA Ratio: 18.5 
               
               
                 Hydrogen Production: 13.5 times Wild Type 
               
               
                 Amino Acid Sequence (SEQ ID NO: 225): 
               
               
                 MAKTIILNGNEVHTDKDITILELARENNVDIPTLCFLKDCGNFGKCGVCMVEVEG 
               
               
                   
               
               
                 KGFRAACVAKVEDGMVINTESDEVKERIKKRVSMLLDKHEFKCGQCSRRENCEFLKLVI 
               
               
                   
               
               
                 KTKAKASKPFLPEDKDALVDNRSKAIVIDRSKCVLCGRCEAACKTKTGTGAISICKSESG 
               
               
                   
               
               
                 RIVQATGGKCFDDINCLLCGQCVAACPVGALTEKTHVDRVKEALEDPNKHVIVAMAPS 
               
               
                   
               
               
                 VRTAMGELFKMGYGKDVTGKLYTALRMLGFDKVFDINFGADMTIMEEATELLGRVKN 
               
               
                   
               
               
                 NGPFPMFTSCCPAWVRQVENYYPEFLENLSSAKSPQQIFGAASKTYYPQISGISAKDVFT 
               
               
                   
               
               
                 VTIMPCTAKKFEADREEMYNEGIKNIDAVLTTRELAKMIKDAKINFANLEDEQADPAMG 
               
               
                   
               
               
                 EYTGAGVIFGATGGVMEAALRTAKDFVEDKDLTDIEYTQIRGLQGIKEATVEIGGENYN 
               
               
                   
               
               
                 VAVINGAANLAEFMNSGKILEKNYHFIEVMACPGGCVNGGGQPHVSAKEREKVDVRTV 
               
               
                   
               
               
                 RASVLYNQDKNLEKRKSHKNTALLNMYYDYMGAPGQGKAHELLHLKYNK 
               
               
                   
               
               
                 DNA Sequence (SEQ ID NO: 226): 
               
               
                 atggccaaaacaataatcttaaatggcaatgaagtgcatacagataaagatattactatccttgagctagcaagagaaaataatgta 
               
               
                   
               
               
                 gatatcccaacactctgctttttaaaggattgtggcaattttggaaaatgcggagtctgtatggtagaggtagaaggcaagggctttagagctg 
               
               
                   
               
               
                 cttgtgttgccaaagttgaagatggaatggtaataaacacagaatccgatgaagtaaaagaacgaatcaaaaaaagagtttcaatgcttcttg 
               
               
                   
               
               
                 ataagcatgaatttaaatgtggacaatgttctagaagagaaaattgtgaattccttaaacttgtaataaagacaaaagcaaaagcttcaaaacc 
               
               
                   
               
               
                 atttttaccagaagataaggatgctctagttgataatagaagtaaggctattgtaattgacagatcaaaatgtgtgctatgcggaagatgtgaag 
               
               
                   
               
               
                 cagcatgtaaaacaaagacaggtacaggagctatttcaatttgtaagagtgaatcaggaagaatagtgcaagcaacaggcggaaagtgctt 
               
               
                   
               
               
                 tgatgatacaaattgtttattatgtggacaatgcgttgcagcatgtccagtaggagctttaactgaaaaaacacacgttgatagagttaaagaag 
               
               
                   
               
               
                 cattagaagatcctaataagcatgtcattgttgcaatggctccatcagtaagaactgctatgggcgaattattcaaaatgggatatggaaaaga 
               
               
                   
               
               
                 tgtaacaggaaaactatatactgcacttagaatgttaggctttgataaagtatttgatataaactttggtgcagatatgactataatggaagaagc 
               
               
                   
               
               
                 tactgaacttttaggcagagttaaaaataatggaccattcccaatgtttacttcatgttgtccggcatgggttagacaagtggaaaattattaccc 
               
               
                   
               
               
                 agaatttttagaaaacttatcatcagctaaatcaccacaacaaatatttggtgcagcaagcaaaacatactatcctcaaatatcaggtataagtg 
               
               
                   
               
               
                 ctaaagatgtatttactgttacaataatgccttgtacagcaaagaaatttgaggctgatagagaagaaatgtataatgagggaattaaaaatata 
               
               
                   
               
               
                 gatgcagtacttactacaagagaattagcaaaaatgattaaagatgcaaagattaattttgctaatttagaagacgaacaagctgatccagcaa 
               
               
                   
               
               
                 tgggagaatacactggggctggagttatattcggagctacaggtggagttatggaagcagcacttagaactgctaaggatttcgttgaagata 
               
               
                   
               
               
                 aagatttaactgatatagaatatacacaaataagaggattacaaggaataaaagaggctacagtagaaattggtggagaaaattataacgtag 
               
               
                   
               
               
                 ctgtaattaatggtgcagcaaacttagctgaattcatgaatagcggtaaaatccttgaaaagaactatcattttattgaagtaatggcttgcccag 
               
               
                   
               
               
                 gcggatgtgtaaatggtggaggacaaccacacgtaagtgcaaaggaaagagaaaaagtagatgttagaactgtaagagcatctgttttatat 
               
               
                   
               
               
                 aaccaagataaaaatttagagaagagaaaatcacataaaaatacagcattattaaatatgtactatgattatatgggagctccaggacaagga 
               
               
                   
               
               
                 aaagctcatgaattattacacttaaaatacaataaa 
               
               
                   
               
               
                   Chlamydomonas reinhardtii  (CrA1) leader sequence: 
               
               
                 Amino Acid Sequence (SEQ ID NO: 227): 
               
               
                 MSALVLKPCAAVSIRGSSCRARQVAPRAPLAASTVRVALATLEAPARRLGNVACAA 
               
               
                   
               
               
                 DNA Sequence (SEQ ID NO: 228): 
               
               
                 atgtcggcgctcgtgctgaagccctgcgcggccgtgtctattcgcggcagctcctgcagggcgcggcaggtcgccccccgcg 
               
               
                   
               
               
                 ctccgctcgcagccagcaccgtgcgtgtagcccttgcaacacttgaggcgcccgcacgccgcctaggcaacgtcgcttgcgcggct 
               
             
          
         
       
     
     
       
         
               
               
             
           
               
                 Algal Hydrogenase Mutant C7: 
                   
               
               
                 Positive/Negative EPSA Ratio: 1.81 
               
               
                 Hydrogen Production: 1.99 times Wild Type (tested w/o leader sequence) 
               
               
                 Amino Acid Sequence (as tested; SEQ ID NO: 229): 
               
               
                 MAAPAAEAPLSHVQQALAELAKPKDDPTRKHVCVQVAPAVRVAIAETLGLAPG 
               
               
                   
               
               
                 ATTPKQLAEGLRRLGFDEVFDTLFGADLTIMEEGSELLHRLTEHLEAHPHSDEPLPMFTS 
               
               
                   
               
               
                 CCPGWIAMLEKSYPDLIPYVSSCKSPQMMLAAMVKSYLAEKKGIAPKDMVMVSIMPCT 
               
               
                   
               
               
                 RKQSEADRDWFCVDADPTLRQLDHVITTVELGNIFKERGINLAELPEGEWDNPMGVGS 
               
               
                   
               
               
                 GAGVLFGTTGGVMEAALRTAYELFTGTPLPRLSLSEVRGMDGIKETNITMVPAPGSKFE 
               
               
                   
               
               
                 ELLKHRAAARAEAAAHGTPGPLAWDGGAGFTSEDGRGGITLRVAVANGLGNAKKLITK 
               
               
                   
               
               
                 MQAGEAKYDFVEIMACPAGCVGGGGQPRSTDKAITQKRQAALYDLDERNTLRRSHEN 
               
               
                   
               
               
                 EAVNQLYKEFLGEPLSHRAHELLHTHYVPGGAEADA 
               
               
                   
               
               
                 DNA Sequence (as tested; SEQ ID NO: 230): 
               
               
                 atggccgcacccgctgcggaggcgcctttgagtcatgtccagcaggcgctcgccgagcttgccaagcccaaggacgacccca 
               
               
                   
               
               
                 cgcgcaagcacgtctgcgtgcaggtggctccggccgttcgtgtcgctattgccgagaccctgggcctggcgccgggcgccaccaccccc 
               
               
                   
               
               
                 aagcagctggccgagggcctccgccgcctcggctttgacgaggtgtttgacacgctgtttggcgccgacctgaccatcatggaggagggc 
               
               
                   
               
               
                 agcgagctgctgcaccgcctcaccgagcacctggaggcccacccgcactccgacgagccgctgcccatgttcaccagctgctgccccg 
               
               
                   
               
               
                 gctggatcgctatgctggagaaatcttacccggacctgatcccctacgtgagcagctgcaagagcccccagatgatgctggcggcaatggt 
               
               
                   
               
               
                 caagtcctacctagcggaaaagaagggcatcgcgccaaaggacatggtcatggtgtccatcatgccctgcacgcgcaagcagtcggagg 
               
               
                   
               
               
                 ctgaccgcgactggttctgtgtggacgccgaccccaccctgcgccagctggaccacgtcatcaccaccgtggagctgggcaacatcttca 
               
               
                   
               
               
                 aggagcgcggcatcaacctggccgagctgcccgagggcgagtgggacaatccaatgggcgtgggctcgggcgccggcgtgctgttcg 
               
               
                   
               
               
                 gcaccaccggcggtgtcatggaggcggcgctgcgcacggcctatgagctgttcacgggcacgccgctgccgcgcctgagcctgagcga 
               
               
                   
               
               
                 ggtgcgcggcatggacggcatcaaggagaccaacatcacaatggtgcccgcgcccgggtccaagtttgaggagctgctgaagcaccgc 
               
               
                   
               
               
                 gccgccgcgcgcgccgaggccgccgcgcacggcacccccgggccgctggcctgggacggcggcgcgggcttcaccagcgaggac 
               
               
                   
               
               
                 ggcaggggcggcatcacactgcgcgtggccgtggccaacgggctgggcaacgccaagaagctgatcaccaagatgcaggccggcga 
               
               
                   
               
               
                 ggccaagtacgactttgtggagatcatggcctgccccgcgggctgtgtgggcggcggcggccagccccgctccaccgacaaggccatc 
               
               
                   
               
               
                 acgcagaagcggcaggcggctctgtacgacctggacgagcgcaacacgctgcgccgcagccacgaaaacgaggcggtcaaccagct 
               
               
                   
               
               
                 gtacaaggagttcctgggcgagcccctgtcccaccgcgcccacgagctgctgcacacccactacgtgccaggcggcgccgaggccgat 
               
               
                   
               
               
                 gcttag 
               
               
                   
               
               
                 Algal Hydrogenase Mutant C80: 
               
               
                 Positive/Negative EPSA Ratio: 1.79 
               
               
                 Hydrogen Production: 2.53 times Wild Type (tested w/o leader sequence) 
               
               
                 Amino Acid Sequence (as tested; SEQ ID NO: 231): 
               
               
                 MAAPAAEAPLSHVQQALAELAKPKDDPTRKHVCVQVAPAVRVAIAETLGLAPG 
               
               
                   
               
               
                 ATTPKQLAEGLRRLGFDEVFDTLFGADLTIMEEGSELLHRLTEHLEAHPHSDEPLPMFTS 
               
               
                   
               
               
                 CCPGWIAMLEKSYPDLIPYVSSCKSPQMMLAAMVKSYLAEKKGIAPKDMVMVSIMPCT 
               
               
                   
               
               
                 RKQSEADRDWFCVDADPTLRQLDHVITTVELGNIFKERGINLAELPEGEWDNPMGVGS 
               
               
                   
               
               
                 GAGVLFGTTGGVMEAALRTAYELFTGTPLPRLSLSEVRGMDGIKETNITMVPAPGSKFE 
               
               
                   
               
               
                 ELLKHRAAARAEAAAHGTPGPLAWDGGAGFTSEDGKGGLKLRVAVANGLGNAKKLIG 
               
               
                   
               
               
                 KMVSGEAKYDFVEIMACPAGCVGGGGQPRSTDKQITQKRQAALYNLDEKSTLRRSHEN 
               
               
                   
               
               
                 PSIRELYDTYLGEPLGHKAHELLHTHYVAGGVEEKDEKK 
               
               
                   
               
               
                 DNA Sequence (as tested; SEQ ID NO: 232): 
               
               
                 atggccgcacccgctgcggaggcgcctttgagtcatgtccagcaggcgctcgccgagcttgccaagcccaaggacgacccca 
               
               
                   
               
               
                 cgcgcaagcacgtctgcgtgcaggtggctccggccgttcgtgtcgctattgccgagaccctgggcctggcgccgggcgccaccaccccc 
               
               
                   
               
               
                 aagcagctggccgagggcctccgccgcctcggctttgacgaggtgtttgacacgctgtttggcgccgacctgaccatcatggaggagggc 
               
               
                   
               
               
                 agcgagctgctgcaccgcctcaccgagcacctggaggcccacccgcactccgacgagccgctgcccatgttcaccagctgctgccccg 
               
               
                   
               
               
                 gctggatcgctatgctggagaaatcttacccggacctgatcccctacgtgagcagctgcaagagcccccagatgatgctggcggcaatggt 
               
               
                   
               
               
                 caagtcctacctagcggaaaagaagggcatcgcgccaaaggacatggtcatggtgtccatcatgccctgcacgcgcaagcagtcggagg 
               
               
                   
               
               
                 ctgaccgcgactggttctgtgtggacgccgaccccaccctgcgccagctggaccacgtcatcaccaccgtggagctgggcaacatcttca 
               
               
                   
               
               
                 aggagcgcggcatcaacctggccgagctgcccgagggcgagtgggacaatccaatgggcgtgggctcgggcgccggcgtgctgttcg 
               
               
                   
               
               
                 gcaccaccggcggtgtcatggaggcggcgctgcgcacggcctatgagagttcacgggcacgccgctgccgcgcctgagcctgagcga 
               
               
                   
               
               
                 ggtgcgcggcatggacggcatcaaggagaccaacatcacaatggtgcccgcgcccgggtccaagtttgaggagctgctgaagcaccgc 
               
               
                   
               
               
                 gccgccgcgcgcgccgaggccgccgcgcacggcaccccegggccgctggcctgggacggcggcgcgggcttcaccagcgaggac 
               
               
                   
               
               
                 ggcaagggcggcctgaagctgcgggtggcggtggcgaacggcctgggcaacgccaagaagctgatcggcaagatggtatctggcgag 
               
               
                   
               
               
                 gccaagtacgacttcgtagaaatcatggcctgccctgccggctgcgtgggcggcggcggccagccccgctccaccgacaagcagatca 
               
               
                   
               
               
                 cccagaagcggcaggcggcgctgtacaacctggacgagaagtccacgctgcgccgcagccacgagaacccgtccatccgcgagctgt 
               
               
                   
               
               
                 acgacacgtacctcggagagccgctgggccacaaggcgcacgagagctgcacacccactacgtggccggcggcgtggaggagaag 
               
               
                   
               
               
                 gacgagaagaagtga 
               
               
                   
               
               
                 Algal Hydrogenase Mutant C81: 
               
               
                 Positive/Negative EPSA Ratio: 2.73 
               
               
                 Hydrogen Production: 1.78 times Wild Type (tested w/o leader sequence) 
               
               
                 Amino Acid Sequence (as tested; SEQ ID NO: 233): 
               
               
                 MAAPAAEAPLSHVQQALAELAKPKDDPTRKHVCVQVAPAVRVAIAETLGLAPG 
               
               
                   
               
               
                 ATTPKQLAEGLRRLGFDEVFDTLFGADLTIMEEGSELLHRLTEHLEAHPHSDEPLPMFTS 
               
               
                   
               
               
                 CCPGWIAMLEKSYPDLIPYVSSCKSPQMMLAAMVKSYLAEKKGIAPKDMVMVSIMPCT 
               
               
                   
               
               
                 RKQSEADRDWFCVDADPTLRQLDHVITTVELGNIFKERGINLAELPEGEWDNPMGVGS 
               
               
                   
               
               
                 GAGVLFGTTGGVMEAALRTAYELFTGTPLPRLSLSEVRGLDGIKEASVTLVPAPGSKFA 
               
               
                   
               
               
                 ELVAERLAHKVEEAAAAEAAAAVEGAVKPPIAYDGGQGFSTDDGRGGITLRVAVANGL 
               
               
                   
               
               
                 GNAKKLITKMQAGEAKYDFVEIMACPAGCVGGGGQPRSTDKAITQKRQAALYNLDEKS 
               
               
                   
               
               
                 TLRRSHENPSIRELYDTYLGEPLGHKAHELLHTHYVAGGVEEKDEKK 
               
               
                   
               
               
                 DNA Sequence (as tested; SEQ ID NO: 234): 
               
               
                 atggccgcacccgctgcggaggcgcctttgagtcatgtccagcaggcgctcgccgagcttgccaagcccaaggacgacccca 
               
               
                   
               
               
                 cgcgcaagcacgtctgcgtgcaggtggctccggccgttcgtgtcgctattgccgagaccctgggcctggcgccgggcgccaccaccccc 
               
               
                   
               
               
                 aagcagctggccgagggcctccgccgcctcggctttgacgaggtgtttgacacgctgtttggcgccgacctgaccatcatggaggagggc 
               
               
                   
               
               
                 agcgagctgctgcaccgcctcaccgagcacctggaggcccacccgcactccgacgagccgctgcccatgttcaccagctgctgccccg 
               
               
                   
               
               
                 gctggatcgctatgctggagaaatcttacccggacctgatccectacgtgagcagctgcaagagcccccagatgatgctggcggcaatggt 
               
               
                   
               
               
                 caagtcctacctagcggaaaagaagggcatcgcgccaaaggacatggtcatggtgtccatcatgccctgcacgcgcaagcagtcggagg 
               
               
                   
               
               
                 ctgaccgcgactggttctgtgtggacgccgaccccaccctgcgccagctggaccacgtcatcaccaccgtggagctgggcaacatcttca 
               
               
                   
               
               
                 aggagcgcggcatcaacctggccgagctgcccgagggcgagtgggacaatccaatgggcgtgggctcgggcgccggcgtgctgttcg 
               
               
                   
               
               
                 gcaccaccggcggtgtcatggaggcggcgctgcgcacggcctatgagctgttcacgggcacgccgctgccgcgcctgagcctgagcga 
               
               
                   
               
               
                 ggtgcgcggcttggacggcatcaaggaggcgtccgtgacgctggtccccgctccgggctccaagttcgccgagctggtggcggagcgc 
               
               
                   
               
               
                 ctggcgcacaaggtcgaggaggcggccgcggctgaggcggcggcggcggtggagggcgccgtgaagccgcccatcgcgtacgacg 
               
               
                   
               
               
                 gcggccagggtttctccacggatgacggcaggggcggcatcacactgcgcgtggccgtggccaacgggctgggcaacgccaagaagc 
               
               
                   
               
               
                 tgatcaccaagatgcaggccggcgaggccaagtacgactttgtggagatcatggcctgccccgcgggctgtgtgggcggcggcggcca 
               
               
                   
               
               
                 gccccgctccaccgacaaggccatcacgcagaagcggcaggeggcgctgtacaacctggacgagaagtccacgctgcgccgcagcca 
               
               
                   
               
               
                 cgagaacccgtccatccgcgagctgtacgacacgtacctcggagagccgctgggccacaaggcgcacgagctgctgcacacccactac 
               
               
                   
               
               
                 gtggccggcggcgtggaggagaaggacgagaagaagtga 
               
             
          
         
       
     
     Hydrogen content in the vapor space versus time is shown in  FIG. 35  for the best performing bacterial mutant (C16A) identified so far, along with the results for the wild type control. Hydrogen produced by the wild type was noted to be linear with time, with the slope varying proportionally to the initial concentration of sodium dithionate (SD). This indicates a first order reaction, with SD being the rate controlling reactant. After four hours of reaction time, under the standard conditions, the wild-type control ceased to hydrogen production. Comparable results were observed with mutants that produced hydrogen at a rate near that of the wild-type control. 
     In contrast, hydrogen production by the C16A mutant increased exponentially with time over the four hours of evaluation indicating a higher order reaction mechanism. C16A continued to produce hydrogen up to 71 hours, albeit at a much slower rate due to ever decreasing amounts of SD in the batch reactor. Other mutants which displayed much higher rates of hydrogen production than the wild type control also yielded hydrogen production correlations versus time which were similar to that observed for C16A. 
     Under standard condition, the maximum hydrogen production at four hours for the wild-type control was 25 counts. The time required for the C16A mutant to reach 25 counts under the same conditions was about 0.01 hours. Hence, the efficiency of C16A, the best performing mutant, is about 400 times that of the wild type. 
     The computational model described herein which uses the ratio of positive to negative electrostatic potential surface areas (EPSA) as the correlating parameter, fit all the results with a Chi squared error of 0.41 ( FIG. 36 ). The estimated maximum bacterial mutant efficiency, mutant/wild-type hydrogen production ratio, is about 540 under standard conditions. Initial results have indicated that nearly double this ratio can be obtained with higher initial concentrations of SD and methyl viologen MV. 
     Over 100 algal mutants were experimentally evaluated for hydrogen production. Among the mutants tested, about 35% have yielded hydrogen, with six mutants yielding hydrogen in excess of that produced by the algal wild-type control. All algal hydrogenases were evaluated using the same operating conditions as those used in the bacterial experimental evaluations 
     Most of the algal hydrogenases tested exhibit EPSA ratios in the range of 1-5. The best algal mutant (C80) thus far identified had a hydrogen production rate of 2.5 times that of the algal control and exhibited a linear hydrogen production with time. Further testing of the algal mutants and identification of the best candidates using the computational model, is expected to identify algal mutants showing an exponential rise in hydrogen production, as was observed with the best bacterial mutant hydrogenases. Isolation and testing of algal hydrogenases with EPSA ratios in the expected optimal range has commenced. 
     One skilled in the art would readily appreciate that the methods described in the present disclosure are well adapted to carry out the objects and obtain the ends and advantages mentioned, as well as those inherent therein. The methods, procedures, treatments, described herein are merely representative and exemplary, and are not intended as limitations on the scope of the invention. It will be readily apparent to one skilled in the art that varying substitutions and modifications may be made to the present disclosure disclosed herein without departing from the scope and spirit of the invention. 
     All patents and publications mentioned in the specification are indicative of the levels of those skilled in the art to which the present disclosure pertains. All patents and publications are herein incorporated by reference to the same extent as if each individual publication was specifically and individually indicated to be incorporated by reference.