Comparing photosynthetic and photovoltaic efficiencies is not a simple issue. Although both processes harvest the energy in sunlight, they operate in distinctly different ways and produce different types of products: biomass or chemical fuels in the case of natural photosynthesis and nonstored electrical current in the case of photovoltaics. In order to find common ground for evaluating energy-conversion efficiency, we compare natural photosynthesis with present technologies for photovoltaic-driven electrolysis of water to produce hydrogen. Photovoltaic-driven electrolysis is the more efficient process when measured on an annual basis, yet short-term yields for photosynthetic conversion under optimal conditions come within a factor of 2 or 3 of the photovoltaic benchmark. We consider opportunities in which the frontiers of synthetic biology might be used to enhance natural photosynthesis for improved solar energy conversion efficiency.

Protein stability and function relies on residues being in their appropriate ionization states at physiological pH. In situ residue pK a s also provides a sensitive measure of the local protein environment. Multiconformation continuum electrostatics (MCCE) combines continuum electrostatics and molecular mechanics force fields in Monte Carlo sampling to simultaneously calculate side chain ionization and conformation. The response of protein to charges is incorporated both in the protein dielectric constant (prot) of four and by explicit conformational changes. The pK a of 166 residues in 12 proteins was determined. The root mean square error is 0.83 pH units, and 90% have errors of 1 pH units whereas only 3% have errors 2 pH units. Similar results are found with crystal and solution structures, showing that the method's explicit conformational sampling reduces sensitivity to the initial structure. The outcome also changes little with protein dielectric constant (prot 4-20). Multiconformation continuum electrostatics titrations show coupling of conformational flexibility and changes in ionization state. Examples are provided where ionizable side chain position (protein G), Asn orientation (lysozyme), His tautomer distribution (RNase A), and phosphate ion binding (RNase A and H) change with pH. Disallowing these motions changes the calculated pK a .

This paper describes a general method to calculate the pKas of ionizable groups in proteins. Electrostatic calculations are carried out using the finite difference Poisson-Boltzmann (FDPB) method. A formal treatment of the calculation of pKas within the framework of the FDPB method is presented. The major change with respect to previous work is the specific incorporation of the complete charge distribution of both the neutral and charged forms of each ionizable group into the formalism. This is extremely important for the treatment of salt bridges. A hybrid statistical mechanical/Tanford-Roxby method, which is found to be significantly faster than previous treatments, is also introduced. This simplifies the problem of summing over the large number of possible ionization states for a complex polyion. Applications to BPTI and serine proteases suggest that the calculations can be quite reliable. However, the necessity of including bound waters in the treatment of the Asp-70... His-31 salt bridge in T4 lysozyme and experience with other proteins suggest that additional factors ultimately need to be considered in a comprehensive treatment of pKas in proteins.

A method for combining calculations of residue pKa's with changes in the position of polar hydrogens has been developed. The Boltzmann distributions of proton positions in hydroxyls and neutral titratable residues are found in the same Monte Carlo sampling procedure that determines the amino acid ionization states at each pH. Electrostatic, Lennard-Jones potentials, and torsion angle energies are considered at each proton position. Many acidic and basic residues are found to have significant electrostatic interactions with either a water- or hydroxyl-containing side chain. Protonation state changes are coupled to reorientation of the neighboring hydroxyl dipoles, resulting in smaller free energy differences between neutral and ionized residues than when the protein is held rigid. Multiconformation pH titration gives better agreement with the experimental pKa's for triclinic hen egg lysozyme than conventional rigid protein calculations. The hydroxyl motion significantly increases the protein dielectric response, making it sensitive to the composition of the local protein structure. More than one conformer per residue is often found at a given pH, providing information about the distribution of low-energy lysozyme structures.

Reaction centers from Rhodobacter sphaeroides were subjected to Monte Carlo sampling to determine the Boltzmann distribution of side-chain ionization states and positions and buried water orientation and site occupancy. Changing the oxidation states of the bacteriochlorophyll dimer electron donor (P) and primary (QA) and secondary (QB) quinone electron acceptors allows preparation of the ground (all neutral), P+QA-, P+QB-, P0QA-, and P0QB- states. The calculated proton binding going from ground to other oxidation states and the free energy of electron transfer from QA-QB to form QAQB- (DeltaGAB) compare well with experiment from pH 5 to pH 11. At pH 7 DeltaGAB is measured as -65 meV and calculated to be -80 meV. With fixed protein positions as in standard electrostatic calculations, DeltaGAB is +170 meV. At pH 7 approximately 0.2 H+/protein is bound on QA reduction. On electron transfer to QB there is little additional proton uptake, but shifts in side chain protonation and position occur throughout the protein. Waters in channels leading from QB to the surface change site occupancy and orientation. A cluster of acids (GluL212, AspL210, and L213) and SerL223 near QB play important roles. A simplified view shows this cluster with a single negative charge (on AspL213 with a hydrogen bond to SerL233) in the ground state. In the QB- state the cluster still has one negative charge, now on the more distant AspL210. AspL213 and SerL223 move so SerL223 can hydrogen bond to QB-. These rearrangements plus other changes throughout the protein make the reaction energetically favorable.