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

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.

Energy metabolism, enzymatic flux capacities, and metabolic flux rates in flying honeybees.

Honeybees rely primarily on the oxidation of hexose sugars to provide the energy required for flight. Measurement of VCO2 (equal to VO2, because VCO2/VO2 = 1.0 during carbohydrate oxidation) during flight allowed estimation of steady-state flux rates through pathways of flight muscle energy metabolism. Comparison of Vmax values for flight muscle hexokinase, phosphofructokinase, citrate synthase, and cytochrome c oxidase with rates of carbon and O2 flux during flight reveal that these enzymes operate closer to Vmax in the flight muscles of flying honeybees than in other muscles previously studied. Possible mechanistic and evolutionary implications of these findings are discussed.

Increased tricarboxylic acid cycle flux in rat brain during forepaw stimulation detected with 1H[13C]NMR.

NMR spectroscopy was used to test recent proposals that the additional energy required for brain activation is provided through nonoxidative glycolysis. Using localized NMR spectroscopic methods, the rate of C4-glutamate isotopic turnover from infused [1-(13)C]glucose was measured in the somatosensory cortex of rat brain both at rest and during forepaw stimulation. Analysis of the glutamate turnover data using a mathematical model of cerebral glucose metabolism showed that the tricarboxylic acid cycle flux [(V(TCA)] increased from 0.49 +/- 0.03 at rest to 1.48 +/- 0.82 micromol/g/min during stimulation (P < 0.01). The minimum fraction of C4-glutamate derived from C1-glucose was approximately 75%, and this fraction was found in both the resting and stimulated rats. Hence, the percentage increase in oxidative cerebral metabolic rate of glucose use (CMRglc) equals the percentage increases in V(TCA) and cerebral metabolic rate of oxygen consumption (CMRO2). Comparison with previous work for the same rat model, which measured total CMRglc [Ueki, M., Linn, F. & Hossman, K. A. (1988) J. Cereb. Blood Flow Metab. 8, 486-4941, indicates that oxidative CMRglc supplies the majority of energy during sustained brain activation.

Extraordinarily high spider densities on islands: flow of energy from the marine to terrestrial food webs and the absence of predation.

Some islands in the Gulf of California support very high densities of spiders. Spider density is negatively correlated with island size; many small islands support 50-200 spiders per m3 of cactus. Energy for these spiders comes primarily from the ocean and not from in situ productivity by land plants. We explicitly connect the marine and terrestrial systems to show that insular food webs represent one endpoint of the marine web. We describe two conduits for marine energy entering these islands: shore drift and seabird colonies. Both conduits are related to island area, having a much stronger effect on smaller islands. This asymmetric effect helps to explain the exceptionally high spider densities on small islands. Although productivity sets the maximal potential densities, predation (by scorpions) limits realized spider abundance. Thus, prey availability and predation act in concert to set insular spider abundance.

Atmospheric energy for subsurface life on Mars?

The location and density of biologically useful energy sources on Mars will limit the biomass, spatial distribution, and organism size of any biota. Subsurface Martian organisms could be supplied with a large energy flux from the oxidation of photochemically produced atmospheric H2 and CO diffusing into the regolith. However, surface abundance measurements of these gases demonstrate that no more than a few percent of this available flux is actually being consumed, suggesting that biological activity driven by atmospheric H2 and CO is limited in the top few hundred meters of the subsurface. This is significant because the available but unused energy is extremely large: for organisms at 30-m depth, it is 2,000 times previous estimates of hydrothermal and chemical weathering energy and far exceeds the energy derivable from other atmospheric gases. This also implies that the apparent scarcity of life on Mars is not attributable to lack of energy. Instead, the availability of liquid water may be a more important factor limiting biological activity because the photochemical energy flux can only penetrate to 100- to 1,000-m depth, where most H2O is probably frozen. Because both atmospheric and Viking lander soil data provide little evidence for biological activity, the detection of short-lived trace gases will probably be a better indicator of any extant Martian life.

Fluctuations and stimulus-induced changes in blood flow observed in individual capillaries in layers 2 through 4 of rat neocortex

Cortical blood flow at the level of individual capillaries and the coupling of neuronal activity to flow in capillaries are fundamental aspects of homeostasis in the normal and the diseased brain. To probe the dynamics of blood flow at this level, we used two-photon laser scanning microscopy to image the motion of red blood cells (RBCs) in individual capillaries that lie as far as 600 μm below the pia mater of primary somatosensory cortex in rat; this depth encompassed the cortical layers with the highest density of neurons and capillaries. We observed that the flow was quite variable and exhibited temporal fluctuations around 0.1 Hz, as well as prolonged stalls and occasional reversals of direction. On average, the speed and flux (cells per unit time) of RBCs covaried linearly at low values of flux, with a linear density of ≈70 cells per mm, followed by a tendency for the speed to plateau at high values of flux. Thus, both the average velocity and density of RBCs are greater at high values of flux than at low values. Time-locked changes in flow, localized to the appropriate anatomical region of somatosensory cortex, were observed in response to stimulation of either multiple vibrissae or the hindlimb. Although we were able to detect stimulus-induced changes in the flux and speed of RBCs in some single trials, the amplitude of the stimulus-evoked changes in flow were largely masked by basal fluctuations. On average, the flux and the speed of RBCs increased transiently on stimulation, although the linear density of RBCs decreased slightly. These findings are consistent with a stimulus-induced decrease in capillary resistance to flow.

Energy Sources for HCO3− and CO2 Transport in Air-Grown Cells of Synechococcus UTEX 6251

Light-dependent inorganic C (Ci) transport and accumulation in air-grown cells of Synechococcus UTEX 625 were examined with a mass spectrometer in the presence of inhibitors or artificial electron acceptors of photosynthesis in an attempt to drive CO2 or HCO3− uptake separately by the cyclic or linear electron transport chains. In the presence of 3-(3,4-dichlorophenyl)-1,1-dimethylurea, the cells were able to accumulate an intracellular Ci pool of 20 mm, even though CO2 fixation was completely inhibited, indicating that cyclic electron flow was involved in the Ci-concentrating mechanism. When 200 μm N,N-dimethyl-p-nitrosoaniline was used to drain electrons from ferredoxin, a similar Ci accumulation was observed, suggesting that linear electron flow could support the transport of Ci. When carbonic anhydrase was not present, initial CO2 uptake was greatly reduced and the extracellular [CO2] eventually increased to a level higher than equilibrium, strongly suggesting that CO2 transport was inhibited and that Ci accumulation was the result of active HCO3− transport. With 3-(3,4-dichlorophenyl)-1,1-dimethylurea-treated cells, Ci transport and accumulation were inhibited by inhibitors of CO2 transport, such as COS and Na2S, whereas Li+, an HCO3−-transport inhibitor, had little effect. In the presence of N,N-dimethyl-p-nitrosoaniline, Ci transport and accumulation were not inhibited by COS and Na2S but were inhibited by Li+. These results suggest that CO2 transport is supported by cyclic electron transport and that HCO3− transport is supported by linear electron transport.

Creating a dynamic picture of the sliding clamp during T4 DNA polymerase holoenzyme assembly by using fluorescence resonance energy transfer

The coordinated assembly of the DNA polymerase (gp43), the sliding clamp (gp45), and the clamp loader (gp44/62) to form the bacteriophage T4 DNA polymerase holoenzyme is a multistep process. A partially opened toroid-shaped gp45 is loaded around DNA by gp44/62 in an ATP-dependent manner. Gp43 binds to this complex to generate the holoenzyme in which gp45 acts to topologically link gp43 to DNA, effectively increasing the processivity of DNA replication. Stopped-flow fluorescence resonance energy transfer was used to investigate the opening and closing of the gp45 ring during holoenzyme assembly. By using two site-specific mutants of gp45 along with a previously characterized gp45 mutant, we tracked changes in distances across the gp45 subunit interface through seven conformational changes associated with holoenzyme assembly. Initially, gp45 is partially open within the plane of the ring at one of the three subunit interfaces. On addition of gp44/62 and ATP, this interface of gp45 opens further in-plane through the hydrolysis of ATP. Addition of DNA and hydrolysis of ATP close gp45 in an out-of-plane conformation. The final holoenzyme is formed by the addition of gp43, which causes gp45 to close further in plane, leaving the subunit interface open slightly. This open interface of gp45 in the final holoenzyme state is proposed to interact with the C-terminal tail of gp43, providing a point of contact between gp45 and gp43. This study further defines the dynamic process of bacteriophage T4 polymerase holoenzyme assembly.

Centripetal cholesterol flux from extrahepatic organs to the liver is independent of the concentration of high density lipoprotein-cholesterol in plasma.

High density lipoproteins (HDLs) play a role in two processes that include the amelioration of atheroma formation and the centripetal flow of cholesterol from the extrahepatic organs to the liver. This study tests the hypothesis that the flow of sterol from the peripheral organs to the liver is dependent upon circulating HDL concentrations. Transgenic C57BL/6 mice were used that expressed variable amounts of simian cholesteryl ester-transfer protein (CETP). The rate of centripetal cholesterol flux was quantitated as the sum of the rates of cholesterol synthesis and low density lipoprotein-cholesterol uptake in the extrahepatic tissues. Steady-state concentrations of cholesterol carried in HDL (HDL-C) varied from 59 to 15 mg/dl and those of apolipoprotein AI from 138 to 65 mg/dl between the control mice (CETPc) and those maximally expressing the transfer protein (CETP+). There was no difference in the size of the extrahepatic cholesterol pools in the CETPc and CETP+ animals. Similarly, the rates of cholesterol synthesis (83 and 80 mg/day per kg, respectively) and cholesterol carried in low density lipoprotein uptake (4 and 3 mg/day per kg, respectively) were virtually identical in the two groups. Thus, under circumstances where the steady-state concentration of HDL-C varied 4-fold, the centripetal flux of cholesterol from the peripheral organs to the liver was essentially constant at approximately 87 mg/day per kg. These studies demonstrate that neither the concentration of HDL-C or apolipoprotein AI nor the level of CETP activity dictates the magnitude of centripetal cholesterol flux from the extrahepatic organs to the liver, at least in the mouse.

Least activation path for protein folding: investigation of staphylococcal nuclease folding by stopped-flow circular dichroism.

Is the pathway of protein folding determined by the relative stability of folding intermediates, or by the relative height of the activation barriers leading to these intermediates? This is a fundamental question for resolving the Levinthal paradox, which stated that protein folding by a random search mechanism would require a time too long to be plausible. To answer this question, we have studied the guanidinium chloride (GdmCl)-induced folding/unfolding of staphylococcal nuclease [(SNase, formerly EC 3.1.4.7; now called microbial nuclease or endonuclease, EC 3.1.31.1] by stopped-flow circular dichroism (CD) and differential scanning microcalorimetry (DSC). The data show that while the equilibrium transition is a quasi-two-state process, kinetics in the 2-ms to 500-s time range are triphasic. Data support the sequential mechanism for SNase folding: U3 <--> U2 <--> U1 <--> N0, where U1, U2, and U3 are substates of the unfolded protein and N0 is the native state. Analysis of the relative population of the U1, U2, and U3 species in 2.0 M GdmCl gives delta-G values for the U3 --> U2 reaction of +0.1 kcal/mol and for the U2 --> U1 reaction of -0.49 kcal/mol. The delta-G value for the U1 --> N0 reaction is calculated to be -4.5 kcal/mol from DSC data. The activation energy, enthalpy, and entropy for each kinetic step are also determined. These results allow us to make the following four conclusions. (i) Although the U1, U2, and U3 states are nearly isoenergetic, no random walk occurs among them during the folding. The pathway of folding is unique and sequential. In other words, the relative stability of the folding intermediates does not dictate the folding pathway. Instead, the folding is a descent toward the global free-energy minimum of the native state via the least activation path in the vast energy landscape. Barrier avoidance leads the way, and barrier height limits the rate. Thus, the Levinthal paradox is not applicable to the protein-folding problem. (ii) The main folding reaction (U1 --> N0), in which the peptide chain acquires most of its free energy (via van der Waals' contacts, hydrogen bonding, and electrostatic interactions), is a highly concerted process. These energy-acquiring events take place in a single kinetic phase. (iii) U1 appears to be a compact unfolded species; the rate of conversion of U2 to U1 depends on the viscosity of solution. (iv) All four relaxation times reported here depend on GdmCl concentrations: it is likely that none involve the cis/trans isomerization of prolines. Finally, a mechanism is presented in which formation of sheet-like chain conformations and a hydrophobic condensation event precede the main-chain folding reaction.

The acceleration and collimation of jets.

I will discuss several issues related to the acceleration, collimation, and propagation of jets from active galactic nuclei. Hydromagnetic stresses provide the best bet for both accelerating relativistic flows and providing a certain amount of initial collimation. However, there are limits to how much "self-collimation" can be achieved without the help of an external pressurized medium. Moreover, existing models, which postulate highly organized poloidal flux near the base of the flow, are probably unrealistic. Instead, a large fraction of the magnetic energy may reside in highly disorganized "chaotic" fields. Such a field can also accelerate the flow to relativistic speeds, in some cases with greater efficiency than highly organized fields, but at the expense of self-collimation. The observational interpretation of jet physics is still hampered by a dearth of unambiguous diagnostics. Propagating disturbances in flows, such as the oblique shocks that may constitute the kiloparsec-scale "knots" in the M87 jet, may provide a wide range of untapped diagnostics for jet properties.