Capturing protein dynamics with time-resolved luminescence spectroscopy

The presented doctoral research utilizes time-resolved spectroscopy to characterize protein dynamics and folding mechanisms. We resolve millisecond-timescale folding by coupling time-resolved fluorescence energy transfer (trFRET) to a continuous flow microfluidic mixer to obtain intramolecular distance distributions throughout the folding process. We have elucidated the folding mechanisms of two cytochromes---one that exhibits two-state folding (cytochrome cb₅₆₂) and one that has both a kinetic refolding intermediate ensemble and a distinct equilibrium unfolding intermediate (cytochrome c₅₅₂). Our data reveal that the distinct structural features of cytochrome c₅₅₂ contribute to its thermostability.

We have also investigated intrachain contact dynamics in unfolded cytochrome cb₅₆₂ by monitoring electron transfer, which occurs as the heme collides with a ruthenium photosensitizer, covalently bound to residues along the polypeptide. Intrachain diffusion for chemically denatured proteins proceeds on the microsecond timescale with an upper limit of 0.1 microseconds. The power-law dependence (slope = -1.5) of the rate constants on the number of peptide bonds between the heme and Ru complex indicate that cytochrome cb₅₆₂ is minimally frustrated.

In addition, we have explored the pathway dependence of electron tunneling rates between metal sites in proteins. Our research group has converted cytochrome b₅₆₂ to a c-type cytochrome with the porphyrin covalently bound to cysteine sidechains. We have investigated the effects of the changes to the protein structure (i.e., increased rigidity and potential new equatorial tunneling pathways) on the electron transfer rates, measured by transient absorption, in a series of ruthenium photosensitizer-modified proteins.

Crystal Coulomb energies: I. Organic donor-acceptor complexes--a review. II. Classical Ewald calculations of the Coulomb binding energy of four donor-acceptor crystals and Wurster's blue perchlorate. III. A rapid-convergence quantum-mechanical formalism for crystal electronic energies

Chapter I

Theories for organic donor-acceptor (DA) complexes in solution and in the solid state are reviewed, and compared with the available experimental data. As shown by McConnell et al. (Proc. Natl. Acad. Sci. U.S., 53, 46-50 (1965)), the DA crystals fall into two classes, the holoionic class with a fully or almost fully ionic ground state, and the nonionic class with little or no ionic character. If the total lattice binding energy 2ε₁ (per DA pair) gained in ionizing a DA lattice exceeds the cost 2ε_o of ionizing each DA pair, ε₁ + ε_o less than 0, then the lattice is holoionic. The charge-transfer (CT) band in crystals and in solution can be explained, following Mulliken, by a second-order mixing of states, or by any theory that makes the CT transition strongly allowed, and yet due to a small change in the ground state of the non-interacting components D and A (or D⁺ and A^-). The magnetic properties of the DA crystals are discussed.

Chapter II

A computer program, EWALD, was written to calculate by the Ewald fast-convergence method the crystal Coulomb binding energy E_C due to classical monopole-monopole interactions for crystals of any symmetry. The precision of E_C values obtained is high: the uncertainties, estimated by the effect on E_C of changing the Ewald convergence parameter η, ranged from ± 0.00002 eV to ± 0.01 eV in the worst case. The charge distribution for organic ions was idealized as fractional point charges localized at the crystallographic atomic positions: these charges were chosen from available theoretical and experimental estimates. The uncertainty in E_C due to different charge distribution models is typically ± 0.1 eV (± 3%): thus, even the simple Hückel model can give decent results.

E_C for Wurster's Blue Perchl orate is -4.1 eV/molecule: the crystal is stable under the binding provided by direct Coulomb interactions. E_C for N-Methylphenazinium Tetracyanoquino- dimethanide is 0.1 eV: exchange Coulomb interactions, which cannot be estimated classically, must provide the necessary binding.

EWALD was also used to test the McConnell classification of DA crystals. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine: 7,7,8,8-Tetracyanoquinodimethan) E_C = -4.0 eV while 2ε_o = 4.6₅ eV: clearly, exchange forces must provide the balance. For the holoionic (1:1)-(N,N,N',N'-Tetramethyl-para- phenylenediamine:para-Chloranil) E_C = -4.4 eV, while 2ε_o = 5.0 eV: again EC falls short of 2ε₁. As a Gedankenexperiment, two nonionic crystals were assumed to be ionized: for (1:1)-(Hexamethyl- benzene:para-Chloranil) E_C = -4.5 eV, 2ε_o = 6.6 eV; for (1:1)- (Napthalene:Tetracyanoethylene) E_C = -4.3 eV, 2ε_o = 6.5 eV. Thus, exchange energies in these nonionic crystals must not exceed 1 eV.

Chapter III

A rapid-convergence quantum-mechanical formalism is derived to calculate the electronic energy of an arbitrary molecular (or molecular-ion) crystal: this provides estimates of crystal binding energies which include the exchange Coulomb inter- actions. Previously obtained LCAO-MO wavefunctions for the isolated molecule(s) ("unit cell spin-orbitals") provide the starting-point. Bloch's theorem is used to construct "crystal spin-orbitals". Overlap between the unit cell orbitals localized in different unit cells is neglected, or is eliminated by Löwdin orthogonalization. Then simple formulas for the total kinetic energy Q^(XT)_λ, nuclear attraction [λ/λ]^XT, direct Coulomb [λλ/λ'λ']^XT and exchange Coulomb [λλ'/λ'λ]^XT integrals are obtained, and direct-space brute-force expansions in atomic wavefunctions are given. Fourier series are obtained for [λ/λ]^XT, [λλ/λ'λ']^XT, and [λλ/λ'λ]^XT with the help of the convolution theorem; the Fourier coefficients require the evaluation of Silverstone's two-center Fourier transform integrals. If the short-range interactions are calculated by brute-force integrations in direct space, and the long-range effects are summed in Fourier space, then rapid convergence is possible for [λ/λ]^XT, [λλ/λ'λ']^XT and [λλ'/λ'λ]^XT. This is achieved, as in the Ewald method, by modifying each atomic wavefunction by a "Gaussian convergence acceleration factor", and evaluating separately in direct and in Fourier space appropriate portions of [λ/λ]^XT, etc., where some of the portions contain the Gaussian factor.

Physical characterization of the rack effect and hydrogen bond networks in blue copper proteins

A summary of previous research is presented that indicates that the purpose of a blue copper protein's fold and hydrogen bond network, aka, the rack effect, enforce a copper(II) geometry around the copper(I) ion in the metal site. In several blue copper proteins, the C-terminal histidine ligand becomes protonated and detaches from the copper in the reduced forms. Mutants of amicyanin from Paracoccus denitrificans were made to alter the hydrogen bond network and quantify the rack effect by pK_a shifts.

The pK_a's of mutant amicyanins have been measured by pH-dependent electrochemistry. P94F and P94A mutations loosen the Northern loop, allowing the reduced copper to adopt a relaxed conformation: the ability to relax drives the reduction potentials up. The measured potentials are 265 (wild type), 380 (P94A), and 415 (P94F) mV vs. NHE. The measured pK_a's are 7.0 (wild type), 6.3 (P94A), and 5.0 (P94F). The additional hydrogen bond to the thiolate in the mutants is indicated by a red-shift in the blue copper absorption and an increase in the parallel hyperfine splitting in the EPR spectrum. This hydrogen bond is invoked as the cause for the increased stability of the C-terminal imidazole.

Melting curves give a measure of the thermal stability of the protein. A thermodynamic intermediate with pH-dependent reversibility is revealed. Comparisons with the electrochemistry and apoamicyanin suggest that the intermediate involves the region of the protein near the metal site. This region is destabilized in the P94F mutant; coupled with the evidence that the imidazole is stabilized under the same conditions confirms an original concept of the rack effect: a high energy configuration is stabilized at a cost to the rest of the protein.

Luminescence properties of electron-hole-droplets in pure and doped germanium

This work contains 4 topics dealing with the properties of the luminescence from Ge.

The temperature, pump-power and time dependences of the photoluminescence spectra of Li-, As-, Ga-, and Sb-doped Ge crystals were studied. For impurity concentrations less than about 10¹⁵cm^-3, emissions due to electron-hole droplets can clearly be identified. For impurity concentrations on the order of 10¹⁶cm^-3, the broad lines in the spectra, which have previously been attributed to the emission from the electron-hole-droplet, were found to possess pump-power and time dependent line shape. These properties show that these broad lines cannot be due to emission of electron-hole-droplets alone. We interpret these lines to be due to a combination of emissions from (1) electron-hole- droplets, (2) broadened multiexciton complexes, (3) broadened bound-exciton, and (4) plasma of electrons and holes. The properties of the electron-hole-droplet in As-doped Ge were shown to agree with theoretical predictions.

The time dependences of the luminescence intensities of the electron-hole-droplet in pure and doped Ge were investigated at 2 and 4.2°K. The decay of the electron-hole-droplet in pure Ge at 4.2°K was found to be pump-power dependent and too slow to be explained by the widely accepted model due to Pokrovskii and Hensel et al. Detailed study of the decay of the electron-hole-droplets in doped Ge were carried out for the first time, and we find no evidence of evaporation of excitons by electron-hole-droplets at 4.2°K. This doped Ge result is unexplained by the model of Pokrovskii and Hensel et al. It is shown that a model based on a cloud of electron-hole-droplets generated in the crystal and incorporating (1) exciton flow among electron-hole-droplets in the cloud and (2) exciton diffusion away from the cloud is capable of explaining the observed results.

It is shown that impurities, introduced during device fabrication, can lead to the previously reported differences of the spectra of laser-excited high-purity Ge and electrically excited Ge double injection devices. By properly choosing the device geometry so as to minimize this Li contamination, it is shown that the Li concentration in double injection devices may be reduced to less than about 10¹⁵cm^-3 and electrically excited luminescence spectra similar to the photoluminescence spectra of pure Ge may be produced. This proves conclusively that electron-hole-droplets may be created in double injection devices by electrical excitation.

The ratio of the LA- to TO-phonon-assisted luminescence intensities of the electron-hole-droplet is demonstrated to be equal to the high temperature limit of the same ratio of the exciton for Ge. This result gives one confidence to determine similar ratios for the electron-hole-droplet from the corresponding exciton ratio in semiconductors in which the ratio for the electron-hole-droplet cannot be determined (e.g., Si and GaP). Knowing the value of this ratio for the electron-hole-droplet, one can obtain accurate values of many parameters of the electron-hole-droplet in these semiconductors spectroscopically.

I. Phosphorescence and the true lifetime of triplet states in fluid solutions. II. Why is condensed oxygen blue?

I. PHOSPHORESCENCE AND THE TRUE LIFETIME OF TRIPLET STATES IN FLUID SOLUTIONS

Phosphorescence has been observed in a highly purified fluid solution of naphthalene in 3-methylpentane (3-MP). The phosphorescence lifetime of C₁₀H₈ in 3-MP at -45 °C was found to be 0.49 ± 0.07 sec, while that of C₁₀D₈ under identical conditions is 0.64 ± 0.07 sec. At this temperature 3-MP has the same viscosity (0.65 centipoise) as that of benzene at room temperature. It is believed that even these long lifetimes are dominated by impurity quenching mechanisms. Therefore it seems that the radiationless decay times of the lowest triplet states of simple aromatic hydrocarbons in liquid solutions are sensibly the same as those in the solid phase. A slight dependence of the phosphorescence lifetime on solvent viscosity was observed in the temperature region, -60° to -18°C. This has been attributed to the diffusion-controlled quenching of the triplet state by residual impurity, perhaps oxygen. Bimolecular depopulation of the triplet state was found to be of major importance over a large part of the triplet decay.

The lifetime of triplet C₁₀H₈ at room temperature was also measured in highly purified benzene by means of both phosphorescence and triplet-triplet absorption. The lifetime was estimated to be at least ten times shorter than that in 3-MP. This is believed to be due not only to residual impurities in the solvent but also to small amounts of impurities produced through unavoidable irradiation by the excitation source. In agreement with this idea, lifetime shortening caused by intense flashes of light is readily observed. This latter result suggests that experiments employing flash lamp techniques are not suitable for these kinds of studies.

The theory of radiationless transitions, based on Robinson's theory, is briefly outlined. A simple theoretical model which is derived from Fano's autoionization gives identical result.

Il. WHY IS CONDENSED OXYGEN BLUE?

The blue color of oxygen is mostly derived from double transitions. This paper presents a theoretical calculation of the intensity of the double transition (a ¹Δ_g) (a ¹Δ_g)←(X ³Σ_g^-) (X ³Σ_g^-), using a model based on a pair of oxygen molecules at a fixed separation of 3.81 Å. The intensity enhancement is assumed to be derived from the mixing (a ¹Δ_g) (a ¹Δ_g) _~~~ (X ³Σ_g^-) (X ³Σ_u^-) and (a ¹Δ_g) (¹Δ_u) _~~~ (X ³Σ_g^-) (X ³Σ_g^-). Matrix elements for these interactions are calculated using a π-electron approximation for the pair system. Good molecular wavefunctions are used for all but the perturbing (B ³Σ_u^-) state, which is approximated in terms of ground state orbitals. The largest contribution to the matrix elements arises from large intramolecular terms multiplied by intermolecular overlap integrals. The strength of interaction depends not only on the intermolecular separation of the two oxygen molecules, but also as expected on the relative orientation. Matrix elements are calculated for different orientations, and the angular dependence is fit to an analytical expression. The theory therefore not only predicts an intensity dependence on density but also one on phase at constant density. Agreement between theory and available experimental results is satisfactory considering the nature of the approximation, and indicates the essential validity of the overall approach to this interesting intensity enhancement problem.