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].

Protein fluctuations are sensed by stimulated infrared echoes of the vibrations of carbon monoxide and azide probes

The correlation functions of the fluctuations of vibrational frequencies of azide ions and carbon monoxide in proteins are determined directly from stimulated photon echoes generated with femtosecond infrared pulses. The asymmetric stretching vibration of azide bound to carbonic anhydrase II exhibits a pronounced evolution of its vibrational frequency distribution on the time scale of a few picoseconds, which is attributed to modifications of the ligand structure through interactions with the nearby Thr-199. When azide is bound in hemoglobin, a more complex evolution of the protein structure is required to interchange the different ligand configurations, as evidenced by the much slower relaxation of the frequency distribution in this case. The time evolution of the distribution of frequencies of carbon monoxide bound in hemoglobin occurs on the ≈10-ps time scale and is very nonexponential. The correlation functions of the frequency fluctuations determine the evolution of the protein structure local to the probe and the extent to which the probe can navigate those parts of the energy landscape where the structural configurations are able to modify the local potential energy function of the probe.

Low frequency vibrational modes in proteins: Changes induced by point-mutations in the protein-cofactor matrix of bacterial reaction centers

As a step toward understanding their functional role, the low frequency vibrational motions (<300 cm−1) that are coupled to optical excitation of the primary donor bacteriochlorophyll cofactors in the reaction center from Rhodobacter sphaeroides were investigated. The pattern of hydrogen-bonding interaction between these bacteriochlorophylls and the surrounding protein was altered in several ways by mutation of single amino acids. The spectrum of low frequency vibrational modes identified by femtosecond coherence spectroscopy varied strongly between the different reaction center complexes, including between different mutants where the pattern of hydrogen bonds was the same. It is argued that these variations are primarily due to changes in the nature of the individual modes, rather than to changes in the charge distribution in the electronic states involved in the optical excitation. Pronounced effects of point mutations on the low frequency vibrational modes active in a protein-cofactor system have not been reported previously. The changes in frequency observed indicate a strong involvement of the protein in these nuclear motions and demonstrate that the protein matrix can increase or decrease the fluctuations of the cofactor along specific directions.

Factors controlling the competition among rotational and vibrational energy transfer channels in glyoxal

The state-to-state transfer of rotational and vibrational energy has been studied for S1 glyoxal (CHOCHO) in collisions with D2, N2, CO and C2H4 using crossed molecular beams. A laser is used to pump glyoxal seeded in He to its S1 zero point level with zero angular momentum about its top axis (K′ = 0). The inelastic scattering to each of at least 26 S1 glyoxal rotational and rovibrational levels is monitored by dispersed S1–S0 fluorescence. Various collision partners are chosen to investigate the relative influences of reduced mass and the collision pair interaction potential on the competition among the energy transfer channels. When the data are combined with that obtained previously from other collision partners whose masses range from 2 to 84 amu, it is seen that the channel competition is controlled primarily by the kinematics of the collisional interaction. Variations in the intermolecular potential play strictly a secondary role.

Variations in frequencies of drug resistance in Plasmodium falciparum

Continual exposure of malarial parasite populations to different drugs may have selected not only for resistance to individual drugs but also for genetic traits that favor initiation of resistance to novel unrelated antimalarials. To test this hypothesis, different Plasmodium falciparum clones having varying numbers of preexisting resistance mechanisms were treated with two new antimalarial agents: 5-fluoroorotate and atovaquone. All parasite populations were equally susceptible in small numbers. However, when large populations of these clones were challenged with either of the two compounds, significant variations in frequencies of resistance became apparent. On one extreme, clone D6 from West Africa, which was sensitive to all traditional antimalarial agents, failed to develop resistance under simple nonmutagenic conditions in vitro. In sharp contrast, the Indochina clone W2, which was known to be resistant to all traditional antimalarial drugs, independently acquired resistance to both new compounds as much as a 1,000 times more frequently than D6. Additional clones that were resistant to some (but not all) traditional antimalarial agents acquired resistance to atovaquone at high frequency, but not to 5-fluoroorotate. These findings were unexpected and surprising based on current views of the evolution of drug resistance in P. falciparum populations. Such new phenotypes, named accelerated resistance to multiple drugs (ARMD), raise important questions about the genetic and biochemical mechanisms related to the initiation of drug resistance in malarial parasites. Some potential mechanisms underlying ARMD phenotypes have public health implications that are ominous.

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.

Two-dimensional infrared spectroscopy of peptides by phase-controlled femtosecond vibrational photon echoes

Two-dimensional infrared spectra of peptides are introduced that are the direct analogues of two- and three-pulse multiple quantum NMR. Phase matching and heterodyning are used to isolate the phase and amplitudes of the electric fields of vibrational photon echoes as a function of multiple pulse delays. Structural information is made available on the time scale of a few picoseconds. Line narrowed spectra of acyl-proline-NH2 and cross peaks implying the coupling between its amide-I modes are obtained, as are the phases of the various contributions to the signals. Solvent-sensitive structural differences are seen for the dipeptide. The methods show great promise to measure structure changes in biology on a wide range of time scales.

Site-specific conformational determination in thermal unfolding studies of helical peptides using vibrational circular dichroism with isotopic substitution

Understanding the detailed mechanism of protein folding requires dynamic, site-specific stereochemical information. The short time response of vibrational spectroscopies allows evaluation of the distribution of populations in rapid equilibrium as the peptide unfolds. Spectral shifts associated with isotopic labels along with local stereochemical sensitivity of vibrational circular dichroism (VCD) allow determination of the segment sequence of unfolding. For a series of alanine-rich peptides that form α-helices in aqueous solution, we used isotopic labeling and VCD to demonstrate that the α-helix noncooperatively unwinds from the ends with increasing temperature. For these blocked peptides, the C-terminal is frayed at 5°C. Ab initio level theoretical simulations of the IR and VCD band shapes are used to analyze the spectra and to confirm the conformation of the labeled components. The VCD signals associated with the labeled residues are amplified by coupling to the nonlabeled parts of the molecule. Thus small labeled segments are detectable and stereochemically defined in moderately large peptides in this report of site-specific peptide VCD conformational analysis.

Intramolecular vibrational dephasing obeys a power law at intermediate times

Experimental intramolecular vibrational dephasing transients for several large organic molecules are reanalyzed. Fits to the experimental data, as well as full numerical quantum calculations with a factorized potential surface for all active degrees of freedom of fluorene indicate that power law decays, not exponentials, occur at intermediate times. The results support a proposal that power law decays describe vibrational dephasing dynamics in large molecules at intermediate times because of the local nature of energy flow.

The nature of the excited state of the reaction center of photosystem II of green plants: A high-resolution fluorescence spectroscopy study

We studied the electronically excited state of the isolated reaction center of photosystem II with high-resolution fluorescence spectroscopy at 5 K and compared the obtained spectral features with those obtained earlier for the primary electron donor. The results show that there is a striking resemblance between the emitting and charge-separating states in the photosystem II reaction center, such as a very similar shape of the phonon wing with characteristic features at 19 and 80 cm−1, almost identical frequencies of a number of vibrational modes, a very similar double-Gaussian shape of the inhomogeneous distribution function, and relatively strong electron-phonon coupling for both states. We suggest that the emission at 5 K originates either from an exciton state delocalized over the inactive branch of the photosystem or from a fraction of the primary electron donor that is long-lived at 5 K. The latter possibility can be explained by a distribution of the free energy difference of the primary charge separation reaction around zero. Both possibilities are in line with the idea that the state that drives primary charge separation in the reaction center of photosystem II is a collective state, with contributions from all chlorophyll molecules in the central part of the complex.

Kinetics of formation of hypoxanthine containing base pairs by HIV-RT: RNA template effects on the base substitution frequencies

Hypoxanthine (H), the deamination product of adenine, has been implicated in the high frequency of A to G transitions observed in retroviral and other RNA genomes. Although H·C base pairs are thermodynamically more stable than other H·N pairs, polymerase selection may be determined in part by kinetic factors. Therefore, the hypoxanthine induced substitution pattern resulting from replication by viral polymerases may be more complex than that predicted from thermodynamics. We have examined the steady-state kinetics of formation of base pairs opposite template H in RNA by HIV-RT, and for the incorporation of dITP during first- and second-strand synthesis. Hypoxanthine in an RNA template enhances the k2app for pairing with standard dNTPs by factors of 10–1000 relative to adenine at the same sequence position. The order of base pairing preferences for H in RNA was observed to be H·C >> H·T > H·A > H·G. Steady-state kinetics of insertion for all possible mispairs formed with dITP were examined on RNA and DNA templates of identical sequence. Insertion of dITP opposite all bases occurs 2–20 times more frequently on RNA templates. This bias for higher insertion frequencies on RNA relative to DNA templates is also observed for formation of mispairs at template A. This kinetic advantage afforded by RNA templates for mismatches and pairing involving H suggests a higher induction of mutations at adenines during first-strand synthesis by HIV-RT.

Bathorhodopsin structure in the room-temperature rhodopsin photosequence: picosecond time-resolved coherent anti-Stokes Raman scattering.

Structural changes in the retinal chromophore during the formation of the bathorhodopsin intermediate (bathoRT) in the room-temperature rhodopsin (RhRT) photosequence (i.e., vision) are examined using picosecond time-resolved coherent anti-Stokes Raman scattering. Specifically, the retinal structure assignable to bathoRT following 8-ps excitation of RhRT is measured via vibrational Raman spectroscopy at a 200-ps time delay where the only intermediate present is bathoRT. Significant differences are observed between the C=C stretching frequencies of the retinal chromophore at low temperature where bathorhodopsin is stabilized and at room temperature where bathorhodopsin is a transient species in the RhRT photosequence. These vibrational data are discussed in terms of the formation of bathoRT, an important step in the energy storage/transduction mechanism of RhRT.

High frequency of sex and equal frequencies of mating types in natural populations of the ciliate Tetrahymena thermophila.

In ciliate protists, sex involves the temporary joining of two cells of compatible mating type, followed by meiosis and exchange of gametic nuclei between conjugants. Reproduction is by asexual binary fission following conjugation. For the many ciliates with fixed multiple mating types, frequency-dependent sex-ratio theory predicts equal frequencies of mating types, if sex is common in nature. Here, we report that in natural populations of Tetrahymena thermophila sexually immature cells, indicative of recent conjugation, are found from spring through fall. In addition, the seven mating types occur in approximately equal frequencies, and these frequencies appear to be maintained by interaction between complex, multiple mat alleles and environmental conditions during conjugation. Such genotype-environment interaction determining mating type frequency is rare among ciliates.