Vibrational population relaxation of carbon monoxide in the heme pocket of photolyzed carbonmonoxy myoglobin: Comparison of time-resolved mid-IR absorbance experiments and molecular dynamics simulations

The vibrational energy relaxation of carbon monoxide in the heme pocket of sperm whale myoglobin was studied by using molecular dynamics simulation and normal mode analysis methods. Molecular dynamics trajectories of solvated myoglobin were run at 300 K for both the δ- and ɛ-tautomers of the distal His-64. Vibrational population relaxation times of 335 ± 115 ps for the δ-tautomer and 640 ± 185 ps for the ɛ-tautomer were estimated by using the Landau–Teller model. Normal mode analysis was used to identify those protein residues that act as the primary “doorway” modes in the vibrational relaxation of the oscillator. Although the CO relaxation rates in both the ɛ- and δ-tautomers are similar in magnitude, the simulations predict that the vibrational relaxation of the CO is faster in the δ-tautomer with the distal His playing an important role in the energy relaxation mechanism. Time-resolved mid-IR absorbance measurements were performed on photolyzed carbonmonoxy hemoglobin (Hb13CO). From these measurements, a T1 time of 600 ± 150 ps was determined. The simulation and experimental estimates are compared and discussed.

Femtosecond dynamics of the forbidden carotenoid S1 state in light-harvesting complexes of purple bacteria observed after two-photon excitation

Time-resolved excited-state absorption intensities after direct two-photon excitation of the carotenoid S1 state are reported for light-harvesting complexes of purple bacteria. Direct excitation of the carotenoid S1 state enables the measurement of subsequent dynamics on a fs time scale without interference from higher excited states, such as the optically allowed S2 state or the recently discovered dark state situated between S1 and S2. The lifetimes of the carotenoid S1 states in the B800-B850 complex and B800-B820 complex of Rhodopseudomonas acidophila are 7 ± 0.5 ps and 6 ± 0.5 ps, respectively, and in the light-harvesting complex 2 of Rhodobacter sphaeroides ≈1.9 ± 0.5 ps. These results explain the differences in the carotenoid to bacteriochlorophyll energy transfer efficiency after S2 excitation. In Rps. acidophila the carotenoid S1 to bacteriochlorophyll energy transfer is found to be quite inefficient (φET1 <28%) whereas in Rb. sphaeroides this energy transfer is very efficient (φET1 ≈80%). The results are rationalized by calculations of the ensemble averaged time constants. We find that the Car S1 → B800 electronic energy transfer (EET) pathway (≈85%) dominates over Car S1 → B850 EET (≈15%) in Rb. sphaeroides, whereas in Rps. acidophila the Car S1 → B850 EET (≈60%) is more efficient than the Car S1 → B800 EET (≈40%). The individual electronic couplings for the Car S1 → BChl energy transfer are estimated to be approximately 5–26 cm−1. A major contribution to the difference between the energy transfer efficiencies can be explained by different Car S1 energy gaps in the two species.

The vibrational energy flow transition in organic molecules: Theory meets experiment

Most large dynamical systems are thought to have ergodic dynamics, whereas small systems may not have free interchange of energy between degrees of freedom. This assumption is made in many areas of chemistry and physics, ranging from nuclei to reacting molecules and on to quantum dots. We examine the transition to facile vibrational energy flow in a large set of organic molecules as molecular size is increased. Both analytical and computational results based on local random matrix models describe the transition to unrestricted vibrational energy flow in these molecules. In particular, the models connect the number of states participating in intramolecular energy flow to simple molecular properties such as the molecular size and the distribution of vibrational frequencies. The transition itself is governed by a local anharmonic coupling strength and a local state density. The theoretical results for the transition characteristics compare well with those implied by experimental measurements using IR fluorescence spectroscopy of dilution factors reported by Stewart and McDonald [Stewart, G. M. & McDonald, J. D. (1983) J. Chem. Phys. 78, 3907–3915].

Spatiotemporal dynamics of guanosine 3′,5′-cyclic monophosphate revealed by a genetically encoded, fluorescent indicator

To investigate the dynamics of guanosine 3′,5′-cyclic monophosphate (cGMP) in single living cells, we constructed genetically encoded, fluorescent cGMP indicators by bracketing cGMP-dependent protein kinase (cGPK), minus residues 1–77, between cyan and yellow mutants of green fluorescent protein. cGMP decreased fluorescence resonance energy transfer (FRET) and increased the ratio of cyan to yellow emissions by up to 1.5-fold with apparent dissociation constants of ≈2 μM and >100:1 selectivity for cGMP over cAMP. To eliminate constitutive kinase activity, Thr516 of cGPK was mutated to Ala. Emission ratio imaging of the indicators transfected into rat fetal lung fibroblast (RFL)-6 showed cGMP transients resulting from activation of soluble and particulate guanylyl cyclase, respectively, by nitric oxide (NO) and C-type natriuretic peptide (CNP). Whereas all naive cells tested responded to CNP, only 68% responded to NO. Both sets of signals showed large and variable (0.5–4 min) latencies. The phosphodiesterase (PDE) inhibitor 3-isobutyl-1-methylxanthine (IBMX) did not elevate cGMP on its own but consistently amplified responses to NO or CNP, suggesting that basal activity of guanylate cyclase is very low and emphasizing the importance of PDEs in cGMP recycling. A fraction of RFL cells showed slowly propagating tides of cGMP spreading across the cell in response to delocalized application of NO. Biolistically transfected Purkinje neurons showed cGMP responses to parallel fiber activity and NO donors, confirming that single-cell increases in cGMP occur under conditions appropriate to cause synaptic plasticity.

Dynamics of Fos-Jun-NFAT1 complexes

Transcription initiation in eukaryotes is controlled by nucleoprotein complexes formed through cooperative interactions among multiple transcription regulatory proteins. These complexes may be assembled via stochastic collisions or defined pathways. We investigated the dynamics of Fos-Jun-NFAT1 complexes by using a multicolor fluorescence resonance energy transfer assay. Fos-Jun heterodimers can bind to AP-1 sites in two opposite orientations, only one of which is populated in mature Fos-Jun-NFAT1 complexes. We studied the reversal of Fos-Jun binding orientation in response to NFAT1 by measuring the efficiencies of energy transfer from donor fluorophores linked to opposite ends of an oligonucleotide to an acceptor fluorophore linked to one subunit of the heterodimer. The reorientation of Fos-Jun by NFAT1 was not inhibited by competitor oligonucleotides or heterodimers. The rate of Fos-Jun reorientation was faster than the rate of heterodimer dissociation at some binding sites. The facilitated reorientation of Fos-Jun heterodimers therefore can enhance the efficiency of Fos-Jun-NFAT1 complex formation. We also examined the influence of the preferred orientation of Fos-Jun binding on the stability and transcriptional activity of Fos-Jun-NFAT1 complexes. Complexes formed at sites where Fos-Jun favored the same binding orientation in the presence and absence of NFAT1 exhibited an 8-fold slower dissociation rate than complexes formed at sites where Fos-Jun favored the opposite binding orientation. Fos-Jun-NFAT1 complexes also exhibited greater transcription activation at promoter elements that favored the same orientation of Fos-Jun binding in the presence and absence of NFAT1. Thus, the orientation of heterodimer binding can influence both the dynamics and promoter selectivity of multiprotein transcription regulatory complexes.

Ultrafast diffraction and structural dynamics: The nature of complex molecules far from equilibrium

Studies of molecular structures at or near their equilibrium configurations have long provided information on their geometry in terms of bond distances and angles. Far-from-equilibrium structures are relatively unknown—especially for complex systems—and generally, neither their dynamics nor their average geometries can be extrapolated from equilibrium values. For such nonequilibrium structures, vibrational amplitudes and bond distances play a central role in phenomena such as energy redistribution and chemical reactivity. Ultrafast electron diffraction, which was developed to study transient molecular structures, provides a direct method for probing the nature of complex molecules far from equilibrium. Here we present our ultrafast electron diffraction observations of transient structures for two cyclic hydrocarbons. At high internal energies of ≈4 eV, these molecules display markedly different behavior. For 1,3,5-cycloheptatriene, excitation results in the formation of hot ground-state structures with bond distances similar to those of the initial structure, but with nearly three times the average vibrational amplitude. Energy is redistributed within 5 ps, but with a negative temperature characterizing the nonequilibrium population. In contrast, the ring-opening reaction of 1,3-cyclohexadiene is shown to result in hot structures with a C—C bond distance of over 1.7 Å, which is 0.2 Å away from any expected equilibrium value. Even up to 400 ps, energy remains trapped in large-amplitude motions comprised of torsion and asymmetric stretching. These studies promise a new direction for studying structural dynamics in nonequilibrium complex systems.

Configurational diffusion down a folding funnel describes the dynamics of DNA hairpins

Elucidating the mechanism of folding of polynucleotides depends on accurate estimates of free energy surfaces and a quantitative description of the kinetics of structure formation. Here, the kinetics of hairpin formation in single-stranded DNA are measured after a laser temperature jump. The kinetics are modeled as configurational diffusion on a free energy surface obtained from a statistical mechanical description of equilibrium melting profiles. The effective diffusion coefficient is found to be strongly temperature-dependent in the nucleation step as a result of formation of misfolded loops that do not lead to subsequent zipping. This simple system exhibits many of the features predicted from theoretical studies of protein folding, including a funnel-like energy surface with many folding pathways, trapping in misfolded conformations, and non-Arrhenius folding rates.

Excited-state dynamics in photosystem II: Insights from the x-ray crystal structure

The heart of oxygenic photosynthesis is photosystem II (PSII), a multisubunit protein complex that uses solar energy to drive the splitting of water and production of molecular oxygen. The effectiveness of the photochemical reaction center of PSII depends on the efficient transfer of excitation energy from the surrounding antenna chlorophylls. A kinetic model for PSII, based on the x-ray crystal structure coordinates of 37 antenna and reaction center pigment molecules, allows us to map the major energy transfer routes from the antenna chlorophylls to the reaction center chromophores. The model shows that energy transfer to the reaction center is slow compared with the rate of primary electron transport and depends on a few bridging chlorophyll molecules. This unexpected energetic isolation of the reaction center in PSII is similar to that found in the bacterial photosystem, conflicts with the established view of the photophysics of PSII, and may be a functional requirement for primary photochemistry in photosynthesis. In addition, the model predicts a value for the intrinsic photochemical rate constant that is 4 times that found in bacterial reaction centers.

Protonation dynamics of the extracellular and cytoplasmic surface of bacteriorhodopsin in the purple membrane.

The dynamics of proton binding to the extracellular and the cytoplasmic surfaces of the purple membrane were measured by laser-induced proton pulses. Purple membranes, selectively labeled by fluorescein at Lys-129 of bacteriorhodopsin, were pulsed by protons released in the aqueous bulk from excited pyranine (8-hydroxy-1,3,6-pyrenetrisulfonate) and the reaction of protons with the indicators was measured. Kinetic analysis of the data imply that the two faces of the membrane differ in their buffer capacities and in their rates of interaction with bulk protons. The extracellular surface of the purple membrane contains one anionic proton binding site per protein molecule with pK = 5.1. This site is within a Coulomb cage radius (approximately 15 A) from Lys-129. The cytoplasmic surface of the purple membrane bears 4-5 protonable moieties (pK = 5.1) that, due to close proximity, function as a common proton binding site. The reaction of the proton with this cluster is at a very fast rate (3.10(10) M-1.s-1). The proximity between the elements is sufficiently high that even in 100 mM NaCl they still function as a cluster. Extraction of the chromophore retinal from the protein has a marked effect on the carboxylates of the cytoplasmic surface, and two to three of them assume positions that almost bar their reaction with bulk protons. The protonation dynamics determined at the surface of the purple membrane is of relevance both for the vectorial proton transport mechanism of bacteriorhodopsin and for energy coupling, not only in halobacteria, but also in complex chemiosmotic systems such as mitochondrial and thylakoid membranes.

Ultra-fast excited state dynamics in green fluorescent protein: multiple states and proton transfer.

The green fluorescent protein (GFP) of the jellyfish Aequorea Victoria has attracted widespread interest since the discovery that its chromophore is generated by the autocatalytic, posttranslational cyclization and oxidation of a hexapeptide unit. This permits fusion of the DNA sequence of GFP with that of any protein whose expression or transport can then be readily monitored by sensitive fluorescence methods without the need to add exogenous fluorescent dyes. The excited state dynamics of GFP were studied following photo-excitation of each of its two strong absorption bands in the visible using fluorescence upconversion spectroscopy (about 100 fs time resolution). It is shown that excitation of the higher energy feature leads very rapidly to a form of the lower energy species, and that the excited state interconversion rate can be markedly slowed by replacing exchangeable protons with deuterons. This observation and others lead to a model in which the two visible absorption bands correspond to GFP in two ground-state conformations. These conformations can be slowly interconverted in the ground state, but the process is much faster in the excited state. The observed isotope effect suggests that the initial excited state process involves a proton transfer reaction that is followed by additional structural changes. These observations may help to rationalize and motivate mutations that alter the absorption properties and improve the photo stability of GFP.

Probing the interaction between two single molecules: fluorescence resonance energy transfer between a single donor and a single acceptor.

We extend the sensitivity of fluorescence resonance energy transfer (FRET) to the single molecule level by measuring energy transfer between a single donor fluorophore and a single acceptor fluorophore. Near-field scanning optical microscopy (NSOM) is used to obtain simultaneous dual color images and emission spectra from donor and acceptor fluorophores linked by a short DNA molecule. Photodestruction dynamics of the donor or acceptor are used to determine the presence and efficiency of energy transfer. The classical equations used to measure energy transfer on ensembles of fluorophores are modified for single-molecule measurements. In contrast to ensemble measurements, dynamic events on a molecular scale are observable in single pair FRET measurements because they are not canceled out by random averaging. Monitoring conformational changes, such as rotations and distance changes on a nanometer scale, within single biological macromolecules, may be possible with single pair FRET.

The electron-furfural scattering dynamics for 63 energetically open electronic states

We report on integral-, momentum transfer-and differential cross sections for elastic and electronically inelastic electron collisions with furfural (C5H4O2). The calculations were performed with two different theoretical methodologies, the Schwinger multichannel method with pseudopotentials (SMCPP) and the independent atom method with screening corrected additivity rule (IAM-SCAR) that now incorporates a further interference (I) term. The SMCPP with N energetically open electronic states (N-open) at either the static-exchange (N-open ch-SE) or the static-exchange-plus-polarisation (N-open ch-SEP) approximation was employed to calculate the scattering amplitudes at impact energies lying between 5 eV and 50 eV, using a channel coupling scheme that ranges from the 1ch-SEP up to the 63ch-SE level of approximation depending on the energy considered. For elastic scattering, we found very good overall agreement at higher energies among our SMCPP cross sections, our IAM-SCAR+I cross sections and the experimental data for furan (a molecule that differs from furfural only by the substitution of a hydrogen atom in furan with an aldehyde functional group). This is a good indication that our elastic cross sections are converged with respect to the multichannel coupling effect for most of the investigated intermediate energies. However, although the present application represents the most sophisticated calculation performed with the SMCPP method thus far, the inelastic cross sections, even for the low lying energy states, are still not completely converged for intermediate and higher energies. We discuss possible reasons leading to this discrepancy and point out what further steps need to be undertaken in order to improve the agreement between the calculated and measured cross sections. (C) 2016 AIP Publishing LLC.

Subband energy in two-band δ-doped semiconductors

We study electron dynamics in a two-band δ-doped semiconductor within the envelope-function approximation. Using a simple parametrization of the confining potential arising from the ionized donors in the δ -doping layer, we are able to find exact solutions of the Dirac-type equation describing the coupling of host bands. As an application we then consider Si δ -doped GaAs. In particular we find that the ground subband energy scales as a power law of the Si concentration per unit area in a wide range of doping levels. In addition, the coupling of host bands leads to a depression of the subband energy due to nonparabolicity effects.

Influence of a uniform current on collective magnetization dynamics in a ferromagnetic metal

We discuss the influence of a uniform current j⃗ on the magnetization dynamics of a ferromagnetic metal. We find that the magnon energy ε(q⃗) has a current-induced contribution proportional to q⃗⋅J→, where J→ is the spin current, and predict that collective dynamics will be more strongly damped at finite j⃗. We obtain similar results for models with and without local moment participation in the magnetic order. For transition metal ferromagnets, we estimate that the uniform magnetic state will be destabilized for j≳109A cm-2. We discuss the relationship of this effect to the spin-torque effects that alter magnetization dynamics in inhomogeneous magnetic systems.

Spin dynamics of current-driven single magnetic adatoms and molecules

A scanning tunneling microscope can probe the inelastic spin excitations of a single magnetic atom in a surface via spin-flip assisted tunneling in which transport electrons exchange spin and energy with the atomic spin. If the inelastic transport time, defined as the average time elapsed between two inelastic spin flip events, is shorter than the atom spin-relaxation time, the scanning tunnel microscope (STM) current can drive the spin out of equilibrium. Here we model this process using rate equations and a model Hamiltonian that describes successfully spin-flip-assisted tunneling experiments, including a single Mn atom, a Mn dimer, and Fe Phthalocyanine molecules. When the STM current is not spin polarized, the nonequilibrium spin dynamics of the magnetic atom results in nonmonotonic dI/dV curves. In the case of spin-polarized STM current, the spin orientation of the magnetic atom can be controlled parallel or antiparallel to the magnetic moment of the tip. Thus, spin-polarized STM tips can be used both to probe and to control the magnetic moment of a single atom.