Modeling the bacterial photosynthetic reaction center. 4. The structural, electrochemical, and hydrogen-bonding properties of 22 mutants of Rhodobacter sphaeroides

Detangling Solvatochromic Effects by the Effective Fragment Potential Method

Understanding molecular interactions in complex systems opens avenues for the efficient design of new materials with target properties. Energy decomposition methods provide a means to obtain a detailed picture of intermolecular interactions. This work introduces a molecular modeling approach for decomposing the solvatochromic shifts of the electronic excited states into the contributions of the individual molecular fragments of the environment surrounding the chromophore. The developed approach is implemented for the QM/EFP (quantum mechanics/effective fragment potential) model that provides a rigorous first-principles-based description of the electronic states of the chromophores in complex polarizable environments. On the example of two model systems, water pentamer and hydrated uracil, we show how the decomposition of the solvatochromic shifts into the contributions of individual solvent water molecules provides a detailed picture of the intermolecular interactions in the ground and excited states of these systems. The analysis also demonstrates the nonadditivity of solute-solvent interactions and the significant contribution of solute polarization to the total values of solvatochromic shifts.

Photosynthetic reaction center variants made via genetic code expansion show Tyr at M210 tunes the initial electron transfer mechanism

Photosynthetic reaction centers (RCs) from Rhodobacter sphaeroides were engineered to vary the electronic properties of a key tyrosine (M210) close to an essential electron transfer component via its replacement with site-specific, genetically encoded noncanonical amino acid tyrosine analogs. High fidelity of noncanonical amino acid incorporation was verified with mass spectrometry and X-ray crystallography and demonstrated that RC variants exhibit no significant structural alterations relative to wild type (WT). Ultrafast transient absorption spectroscopy indicates the excited primary electron donor, P*, decays via a ∼4-ps and a ∼20-ps population to produce the charge-separated state P⁺H_A ^- in all variants. Global analysis indicates that in the ∼4-ps population, P⁺H_A ^- forms through a two-step process, P*→ P⁺B_A ^-→ P⁺H_A ^-, while in the ∼20-ps population, it forms via a one-step P* → P⁺H_A ^- superexchange mechanism. The percentage of the P* population that decays via the superexchange route varies from ∼25 to ∼45% among variants, while in WT, this percentage is ∼15%. Increases in the P* population that decays via superexchange correlate with increases in the free energy of the P⁺B_A ^- intermediate caused by a given M210 tyrosine analog. This was experimentally estimated through resonance Stark spectroscopy, redox titrations, and near-infrared absorption measurements. As the most energetically perturbative variant, 3-nitrotyrosine at M210 creates an ∼110-meV increase in the free energy of P⁺B_A ^- along with a dramatic diminution of the 1,030-nm transient absorption band indicative of P⁺B_A ^- formation. Collectively, this work indicates the tyrosine at M210 tunes the mechanism of primary electron transfer in the RC.

Studying hydrogen bonding and dynamics of the acetylate groups of the Special Pair of Rhodobacter sphaeroides WT

Although the cofactors in the bacterial reaction centre of Rhodobacter sphaeroides wild type (WT) are arranged almost symmetrically in two branches, the light-induced electron transfer occurs selectively in one branch. As origin of this functional symmetry break, a hydrogen bond between the acetyl group of P_L in the primary donor and His-L168 has been discussed. In this study, we investigate the existence and rigidity of this hydrogen bond with solid-state photo-CIDNP MAS NMR methods offering information on the local electronic structure due to highly sensitive and selective NMR experiments. On the time scale of the experiment, the hydrogen bond between P_L and His-L168 appears to be stable and not to be affected by illumination confirming a structural asymmetry within the Special Pair.

Local and Global Electric Field Asymmetry in Photosynthetic Reaction Centers

The origin of unidirectional electron transfer in photosynthetic reaction centers (RCs) has been widely discussed. Despite the high level of structural similarity between the two branches of pigments that participate in the initial electron transfer steps of photosynthesis, electron transfer only occurs along one branch. One possible explanation for this functional asymmetry is the differences in the electrostatic environment between the active and the inactive branches arising from the charges and dipoles of the organized protein structure. We present an analysis of electric fields in the RC of the purple bacterium Rhodobacter sphaeroides using the intrinsic carbonyl groups of the pigments as vibrational reporters whose vibrational frequency shifts can be converted into electric fields based on the vibrational Stark effect and also provide Stark effect data for plant pigments that can be used in future studies. The carbonyl stretches of the isolated pigments show pronounced Stark effects. We use these data, solvatochromism, molecular dynamics simulations, and data in the literature from IR and Raman spectra to evaluate differences in fields at symmetry-related positions, in particular at the 9-keto and 2-acetyl positions of the pigments involved in primary charge separation.

Putting David Craig’s Legacy to Work in Nanotechnology and Biotechnology

David Craig (1919–2015) left us with a lasting legacy concerning basic understanding of chemical spectroscopy and bonding. This is expressed in terms of some of the recent achievements of my own research career, with a focus on integration of Craig’s theories with those of Noel Hush to solve fundamental problems in photosynthesis, molecular electronics (particularly in regard to the molecules synthesized by Maxwell Crossley), and self-assembled monolayer structure and function. Reviewed in particular is the relation of Craig’s legacy to: the 50-year struggle to assign the visible absorption spectrum of arguably the world’s most significant chromophore, chlorophyll; general theories for chemical bonding and structure extending Hush’s adiabatic theory of electron-transfer processes; inelastic electron-tunnelling spectroscopy (IETS); chemical quantum entanglement and the Penrose–Hameroff model for quantum consciousness; synthetic design strategies for NMR quantum computing; Gibbs free-energy measurements and calculations for formation and polymorphism of organic self-assembled monolayers on graphite surfaces from organic solution; and understanding the basic chemical processes involved in the formation of gold surfaces and nanoparticles protected by sulfur-bound ligands, ligands whose form is that of Au0-thiyl rather than its commonly believed AuI-thiolate tautomer.

Putative hydrogen bond to tyrosine M208 in photosynthetic reaction centers from Rhodobacter capsulatus significantly slows primary charge separation

Slow, ∼50 ps, P* → P⁺H_A^– electron transfer is observed in Rhodobacter capsulatus reaction centers (RCs) bearing the native Tyr residue at M208 and the single amino acid change of isoleucine at M204 to glutamic acid. The P* decay kinetics are unusually homogeneous (single exponential) at room temperature. Comparative solid-state NMR of [4′-¹³C]Tyr labeled wild-type and M204E RCs show that the chemical shift of Tyr M208 is significantly altered in the M204E mutant and in a manner consistent with formation of a hydrogen bond to the Tyr M208 hydroxyl group. Models based on RC crystal structure coordinates indicate that if such a hydrogen bond is formed between the Glu at M204 and the M208 Tyr hydroxyl group, the −OH would be oriented in a fashion expected (based on the calculations by Alden et al., J. Phys. Chem.1996, 100, 16761–16770) to destabilize P⁺B_A^– in free energy. Alteration of the environment of Tyr M208 and B_A by Glu M204 via this putative hydrogen bond has a powerful influence on primary charge separation.

QM and QM/MM simulations of proteins

Molecular dynamics simulations of biomolecules have matured into powerful tools of structural biology. In addition to the commonly used empirical force field potentials, quantum mechanical descriptions are gaining popularity for structure optimization and dynamic simulations of peptides and proteins. In this chapter, we introduce methodological developments such as the QM/MM framework and linear-scaling QM that make efficient calculations on large biomolecules possible. We identify the most common scenarios in which quantum descriptions of peptides and proteins are employed, such as structural refinement, force field development, treatment of unusual residues, and predicting spectroscopic and exited state properties. The benefits and shortcomings of QM potentials, in comparison to classical force fields, are discussed, with special emphasis on the sampling problems of protein conformational space. Finally, recent examples of QM/MM calculations in light-sensitive membrane proteins illustrate typical applications of the reviewed methods.

Inter-porphyrin coupling: how strong should it be for molecular electronics applications?

Porphyrins and phthalocyanines have now been assembled in a multitude of different architectures, each of which may be identified with a different scenario of the coupling acting between the porphyrins. The synthetic flexibility of these compounds makes possible the design of particular molecules for specific applications in molecular electronics, both in naturally occurring and synthetic devices. Here, we form an overview of these features and focus on the coupling strength, considering what values are appropriate for different molecular electronics applications. In particular, we focus on model compounds that have been prepared as mimics of naturally occurring photosynthetic functional units, oligoporphyrins molecular wires, and stacked systems in which small changes in geometry can affect significant changes in the inter-porphyrin coupling and hence produce dramatic changes in device properties.

The electronic structure of the primary electron donor of reaction centers of purple bacteria at atomic resolution as observed by photo-CIDNP 13C NMR

Composed of the two bacteriochlorophyll cofactors, P(L) and P(M), the special pair functions as the primary electron donor in bacterial reaction centers of purple bacteria of Rhodobacter sphaeroides. Under light absorption, an electron is transferred to a bacteriopheophytin and a radical pair is produced. The occurrence of the radical pair is linked to the production of enhanced nuclear polarization called photochemically induced dynamic nuclear polarization (photo-CIDNP). This effect can be used to study the electronic structure of the special pair at atomic resolution by detection of the strongly enhanced nuclear polarization with laser-flash photo-CIDNP magic-angle spinning NMR on the carotenoid-less mutant R26. In the electronic ground state, P(L) is strongly disturbed, carrying a slightly negative charge. In the radical cation state, the ratio of total electron spin densities between P(L) and P(M) is 2:1, although it is 2.5:1 for the pyrrole carbons, 2.2:1 for all porphyrinic carbons, and 4:1 for the pyrrole nitrogen. It is shown that the symmetry break between the electronic structures in the electronic ground state and in the radical cation state is an intrinsic property of the special pair supermolecule, which is particularly attributable to a modification of the structure of P(L). The significant difference in electron density distribution between the ground and radical cation states is explained by an electric polarization effect of the nearby histidine.

Stark spectroscopy of mixed-valence systems

Many mixed-valence systems involve two or more states with different electric dipole moments whose magnitudes depend upon the charge transfer distance and the degree of delocalization; these systems can be interconverted by excitation of an intervalence charge transfer transition. Stark spectroscopy involves the interaction between the change in dipole moment of a transition and an electric field, so the Stark spectra of mixed-valence systems are expected to provide quantitative information on the degree of delocalization. In limiting cases, a classical Stark analysis can be used, but in intermediate cases the analysis is much more complex because the field affects not only the band position but also the intrinsic bandshape. Such non-classical Stark effects lead to widely different bandshapes. Several examples of both classes are discussed. Because electric fields are applied to immobilized samples, complications arise from inhomogeneous broadening, along with other effects that limit our ability to extract unique parameters in some cases. In the case of the radical cation of the special pair in photosynthetic reaction centres, where the mixed-valence system is in a very complex but structurally well-defined environment, a detailed analysis can be performed.

Assignment of the Qy absorption spectrum of photosystem-I from Thermosynechococcus elongatus based on CAM-B3LYP calculations at the PW91-optimized protein structure

The Qy absorption spectrum of Photosystem-I from Thermosynecochoccus elongatus (formerly Synecochoccus elongatus) is calculated using the CAM-B3LYP density functional and INDO schemes based on a quantum-mechanically refined structure for the entire photosystem obtained using the PW91 density functional. These methods present a priori predictions of the absorption and linear dichroism spectra and include protein electrostatic effects, short range inductive effects, long-range and short-range exciton couplings, and superexchange effects involving aromatic residues and carotenes. CAM-B3LYP is used as it is the only known density functional that correctly describes the Q bands of chlorophylls, all other methods contaminating them with erroneous charge-transfer excitations. A critical feature is found to be the use of fully optimized heavy-atom coordinates, with those obtained from just X-ray crystallography providing a poor description of the electronic properties of the chromophores. The result is a realistic first-principles prediction of the observed absorption band that identifies the nature of the red-shifted chlorophylls as well as the energies of the reaction-center chlorophylls and the exciton couplings acting between them. The "special pair" appears more like a dimer of dimers than a self-contained functional unit, with the exciton couplings between its members and the accessory chlorophylls exceeding the internal coupling.

Charge Delocalization in the Special-Pair Radical Cation of Mutant Reaction Centers of Rhodobacter sphaeroides from Stark Spectra and Nonadiabatic Spectral Simulations

Stark and absorption spectra for the hole-transfer band of the bacteriochlorophyll special pair in the wild-type and L131LH, M160LH, and L131LH/M160LH mutants of the bacterial reaction center of Rhodobacter sphaeroides are presented, along with extensive analyses based on nonadiabatic spectral simulations. Dramatic changes in the Stark spectra are induced by the mutations, changes that are readily interpreted in terms of the redox-energy asymmetry and degree of charge localization in the special-pair radical cation. The effect of mutagenesis on key properties such as the electronic coupling within the special pair and the reorganization energy associated with intervalence hole transfer are determined for the first time. Results for the L131LH and M160LH/L131LH mutants indicate that these species can be considered as influencing the special pair primarily through modulation of the redox asymmetry, as is usually conceptualized, but M160LH is shown to develop a wide range of effects that can be interpreted in terms of significant mutation-induced structural changes in and around the special pair. The nonadiabatic spectra simulations are performed using both a simple two-state 1-mode and an extensive four-state 70-mode model, which includes the descriptions of additional electronic states and explicitly treats the major vibrational modes involved. Excellent agreement between the two simulation approaches is obtained. The simple model is shown to reproduce key features of the Stark effect of the main intervalence transition, while the extensive model quantitatively reproduces most features of the observed spectra for both the electronic and the phase-phonon regions, thus giving a more comprehensive description of the effect of the mutations on the properties of the special-pair radical cation. These results for a series of closely related mixed-valence complexes show that the Stark spectra provide a sensitive indicator for the properties of the mixed-valence complexes and should serve as an instructive example on the application of nonadiabatic simulations to the study of mixed-valence complexes in general as well as other chemical systems akin to the photosynthetic special pair.

Are AM1 ligand-protein binding enthalpies good enough for use in the rational design of new drugs?

We have examined the performance of semiempirical quantum mechanical methods in solving the problem of accurately predicting protein-ligand binding energies and geometries. Firstly, AM1 and PM3 geometries and binding enthalpies between small molecules that simulate typical ligand-protein interactions were compared with high level quantum mechanical techniques that include electronic correlation (e.g., MP2 or B3LYP). Species studied include alkanes, aromatic systems, molecules including groups with hypervalent sulfur or with donor or acceptor hydrogen bonding capability, as well as ammonium or carboxylate ions. B3LYP/6-311+G(2d,p) binding energies correlated very well with the BSSE corrected MP2/6-31G(d) values. AM1 binding enthalpies also showed good correlation with MP2 values, and their systematic deviation is acceptable when enthalpies are used for the comparison of interaction energies between ligands and a target. PM3 otherwise gave erratic energy differences in comparison to the B3LYP or MP2 approaches. As one would expect, the geometries of the binding complexes showed the known limitations of the semiempirical and DFT methods. AM1 calculations were subsequently applied to a test set consisting of "real" protein active site-ligand complexes. Preliminary results indicate that AM1 could be a valuable tool for the design of new drugs using proteins as templates. This approach also has a reasonable computational cost. The ligand-protein X-ray structures were reasonably reproduced by AM1 calculations and the corresponding AM1 binding enthalpies are in agreement with the results from the "small molecules" test set.

A Unified Description of the Electrochemical, Charge Distribution, and Spectroscopic Properties of the Special-Pair Radical Cation in Bacterial Photosynthesis

We apply our four-state 70-vibration vibronic-coupling model for the properties of the photosynthetic special-pair radical cation to: (1) interpret the observed correlations between the midpoint potential and the distribution of spin density between the two bacteriochlorophylls for 30 mutants of Rhodobacter sphaeroides, (2) interpret the observed average intervalence hole-transfer absorption energies as a function of spin density for six mutants, and (3) simulate the recently obtained intervalence electroabsorption Stark spectrum of the wild-type reaction center. While three new parameters describing the location of the sites of mutation with respect to the special pair are required to describe the midpoint-potential data, a priori predictions are made for the transition energies and the Stark spectrum. In general, excellent predictions are made of the observed quantities, with deviations being typically of the order of twice the experimental uncertainties. A unified description of many chemical and spectroscopic properties of the bacterial reaction center is thus provided. Central to the analysis is the assumption that the perturbations made to the reaction center, either via mutations of protein residues or by application of an external electric field, act only to independently modify the oxidation potentials of the two halves of the special pair and hence the redox asymmetry E0. While this appears to be a good approximation, clear evidence is presented that effects of mutation can be more extensive than what is allowed for. A thorough set of analytical equations describing the observed properties is obtained using the Born-Oppenheimer adiabatic approximation. These equations are generally appropriate for intervalence charge-transfer problems and include, for the first time, full treatment of both symmetric and antisymmetric vibrational motions. The limits of validity of the adiabatic approach to the full nonadiabatic problem are obtained.

Structural, dynamic, and energetic aspects of long-range electron transfer in photosynthetic reaction centers

Intramolecular electron transfer within proteins plays an essential role in biological energy transduction. Electron donor and acceptor cofactors are bound in the protein matrix at specific locations, and protein-cofactor interactions as well as protein conformational changes can markedly influence the electron transfer rates. To assess these effects, we have investigated charge recombination from the primary quinone acceptor to the special pair bacteriochlorophyll dimer in wild-type reaction centers of Rhodobacter sphaeroides and four mutants with widely modified free energy gaps. After light-induced charge separation, the recombination kinetics were measured in the light- and dark-adapted forms of the protein from 10 to 300 K. The data were analyzed by using the spin-boson model, which allowed us to self-consistently determine the electronic coupling energy, the distribution of energy gaps, the spectral density of phonons, and the reorganization energy. The analysis revealed slow changes of the energy gap after charge separation. Interesting correlations of the control parameters governing electron transfer were found and related to structural and dynamic properties of the protein.

Identification of intramembrane hydrogen bonding between 131 keto group of bacteriochlorophyll and serine residue α27 in the LH2 light-harvesting complex

Intramembrane hydrogen bonding and its effect on the structural integrity of purple bacterial light-harvesting complex 2, LH2, have been assessed in the native membrane environment. A novel hydrogen bond has been identified by Raman resonance spectroscopy between a serine residue of the membrane-spanning region of LH2 alpha-subunit, and the C-13(1) keto carbonyl of bacteriochlorophyll (BChl) B850 bound to the beta-subunit. Replacement of the serine by alanine disrupts this strong hydrogen bond, but this neither alters the strongly red-shifted absorption nor the structural arrangement of the BChls, as judged from circular dichroism. It also decreases only slightly the thermal stability of the mutated LH2 in the native membrane environment. The possibility is discussed that weak H-bonding between the C-13(1) keto carbonyl and a methyl hydrogen of the alanine replacing serine(-4) or the imidazole group of the nearby histidine maintains structural integrity in this very stable bacterial light-harvesting complex. A more widespread occurrence of H-bonding to C-13(1) not only in BChl, but also in chlorophyll proteins, is indicated by a theoretical analysis of chlorophyll/polypeptide contacts at <3.5 A in the high-resolution structure of Photosystem I. Nearly half of the 96 chlorophylls have aa residues suitable as hydrogen bond donors to their keto groups.

Combined quantum and molecular mechanics calculations on metalloproteins

The combination of quantum mechanics and molecular mechanics (QM/MM) methods is one of the most promising approaches to study the structure, function and properties of proteins. The number of QM/MM applications on metalloproteins is steadily increasing, especially studies with density functional methods on redox-active metal centres. Recent developments include new parameterised methods to treat covalent bonds between the quantum and classical systems, methods to obtain free energy from QM/MM results, and the combination of quantum chemistry and protein crystallography.