Theoretical investigation of infrared spectra and pocket dynamics of photodissociated carbonmonoxy myoglobin

Molecular dynamics simulations of the photodissociated state of carbonmonoxy myoglobin (MbCO) are presented using a fluctuating charge model for CO. A new three-point charge model is fitted to high-level ab initio calculations of the dipole and quadrupole moment functions taken from the literature. The infrared spectrum of the CO molecule in the heme pocket is calculated using the dipole moment time autocorrelation function and shows good agreement with experiment. In particular, the new model reproduces the experimentally observed splitting of the CO absorption spectrum. The splitting of 3–7 cm−1 (compared to the experimental value of 10 cm−1) can be directly attributed to the two possible orientations of CO within the docking site at the edge of the distal heme pocket (the B states), as previously suggested on the basis of experimental femtosecond time-resolved infrared studies. Further information on the time evolution of the position and orientation of the CO molecule is obtained and analyzed. The calculated difference in the free energy between the two possible orientations (Fe···CO and Fe···OC) is 0.3 kcal mol−1 and agrees well with the experimentally estimated value of 0.29 kcal mol−1. A comparison of the new fluctuating charge model with an established fixed charge model reveals some differences that may be critical for the correct prediction of the infrared spectrum and energy barriers. The photodissociation of CO from the myoglobin mutant L29F using the new model shows rapid escape of CO from the distal heme pocket, in good agreement with recent experimental data. The effect of the protein environment on the multipole moments of the CO ligand is investigated and taken into account in a refined model. Molecular dynamics simulations with this refined model are in agreement with the calculations based on the gas-phase model. However, it is demonstrated that even small changes in the electrostatics of CO alter the details of the dynamics.

A new interpretation of the bonding properties and UV–vis spectra of [M3(CO)12] clusters (M = Ru, Os): a TD-DFT study

DFT and TD-DFT calculations (ADF program) were performed in order to analyze the electronic structure of the [M-3(CO)(12)] clusters (M = Ru, Os) and interpret their electronic spectra. The highest occupied molecular orbitals are M-M bonding (sigma) involving different M-M bonds, both for Ru and Os. They participate in low-energy excitation processes and their depopulation should weaken M-M bonds in general. While the LUMO is M-NI and M-CO anti-bonding (sigma*), the next, higher-lying empty orbitals have a main contribution from CO (pi*) and either a small (Ru) or an almost negligible one (Os) from the metal atoms. The main difference between the two clusters comes from the different nature of these low-energy unoccupied orbitals that have a larger metal contribution in the case of ruthenium. The photochemical reactivity of the two clusters is reexamined and compared to earlier interpretations.

Infrared absorption spectra, radiative efficiencies, and global warming potentials of perfluorocarbons: Comparison between experiment and theory

Experimentally and theoretically determined infrared spectra are reported for a series of straight-chain perfluorocarbons: C2F6, C3F8, C4F10, C5F12, C6F14, and C8F18. Theoretical spectra were determined using both density functional (DFT) and ab initio methods. Radiative efficiencies (REs) were determined using the method of Pinnock et al. (1995) and combined with atmospheric lifetimes from the literature to determine global warming potentials (GWPs). Theoretically determined absorption cross sections were within 10% of experimentally determined values. Despite being much less computationally expensive, DFT calculations were generally found to perform better than ab initio methods. There is a strong wavenumber dependence of radiative forcing in the region of the fundamental C-F vibration, and small differences in wavelength between band positions determined by theory and experiment have a significant impact on the REs. We apply an empirical correction to the theoretical spectra and then test this correction on a number of branched chain and cyclic perfluoroalkanes. We then compute absorption cross sections, REs, and GWPs for an additional set of perfluoroalkenes.

Tuning the self-assembly of the bioactive dipeptide L-carnosine by incorporation of a bulky aromatic substituent

The dipeptide L-carnosine has a number of important biological properties. Here, we explore the effect of attachment of a bulky hydrophobic aromatic unit, Fmoc [N-(fluorenyl-9-methoxycarbonyl)] on the self-assembly of Fmoc-L-carnosine, i.e., Fmoc-Beta-alanine-histidine (Fmoc-BetaAH). It is shown that Fmoc-BetaAH forms well-defined amyloid fibril containing Beta sheets above a critical aggregation concentration, which is determined from pyrene and ThT fluorescence experiments. Twisted fibrils were imaged by cryogenic transmission electron microscopy. The zinc-binding properties of Fmoc-BetaAH were investigated by FTIR and Raman spectroscopy since the formation of metal ion complexes with the histidine residue in carnosine is well-known, and important to its biological roles. Observed changes in the spectra may reflect differences in the packing of the Fmoc-dipeptides due to electrostatic interactions. Cryo-TEM shows that this leads to changes in the fibril morphology. Hydrogelation is also induced by addition of an appropriate concentration of zinc ions. Our work shows that the Fmoc motif can be employed to drive the self-assembly of carnosine into amyloid fibrils.

Origin of electronic absorption spectra of MLCT-excited and one-electron reduced 2,2′-bipyridine and 1,10-phenanthroline complexes

UV–Vis absorption spectra of one-electron reduction products and 3MLCT excited states of [ReICl(CO)3- (N,N)] (N,N = 2,20-bipyridine, bpy; 1,10-phenanthroline, phen) have been measured by low-temperature spectroelectrochemistry and UV–Vis transient absorption spectroscopy, respectively, and assigned by open-shell TD-DFT calculations. The characters of the electronic transitions are visualized and analyzed using electron density redistribution maps. It follows that reduced and excited states can be approximately formulated as [ReICl(CO)3(N,Nÿ)]ÿ and ⁄[ReIICl(CO)3(N,Nÿ)], respectively. UV–Vis spectra of the reduced complexes are dominated by IL transitions, plus weaker MLCT contributions. Excited-state spectra show an intense band in the UV region of 50% IL origin mixed with LMCT (bpy, 373 nm) or MLCT (phen, 307 nm) excitations. Because of the significant IL contribution, this spectral feature is akin to the principal IL band of the anions. In contrast, the excited-state visible spectral pattern arises from predominantly LMCT transitions, any resemblance with the reduced-state visible spectra being coincidental. The Re complexes studied herein are representatives of a broad class of metal a-diimines, for which similar spectroscopic behavior can be expected.

Structure and spectra of dichloro(hydroxy-methoxy-di(2-pyridyl)methane)copper(II)

The non-electrolyte dichloro(hydroxy-methoxy-di(2-pyridylmethane)copper(II), resulting from the reaction of di(2-pyridyl)ketone and copper(II) chloride in methanol solution, was isolated and characterized and its structure was determined by X-ray diffraction. The pyridyl nitrogens and the chloride anions virtually from a basal plane in which lies the copper atom, while the oxygen of the methoxy group is in an apical position at a distance of 2.497 (3)Å. The nitrogenous base adopts the boat conformation with the pyridyl rings forming a dihedral angle of 108.72 (14)°. The nearest interatomic copper distance of 3.940(3)Å precludes copper-copper interactions, while the proximity of copper to the out-of-plane chlorine atoms [3.109(3)Å] suggests weakly bound chloro-bridged dimers. Spectral changes indicate that protic molecules displace the methoxy group and water affords the corresponding 1,1-diol.

Influence of hydrogen bonding in the structure of mixed-ligand copper(II) chelates. Spectra and structure of [aquo((2,2'-dipyridylamine)-(3-chloro-2,4-pentanedionato)copper(II)] nitrate

The IR, the ligand field spectra and the crystal structure of the mixed-ligand compound [(aquo)2,2P1 , a = 8.718(5), b = 9.407(5), c = 13.484 (7) Å, = 94.17(4)°, = 105.12(5)°, = 119.75(5)°, Z = 2, R = 0.0332, R W = 0.0869).

Synthesis, characterisation, Mössbauer spectra, and structures of some trinuclear iron-tin clusters

Reactions of [Fe3(CO)12] with diaryltin species SnR2(R1= 2,4,6-triisopropylphenyl, R2= 2,6-diethylphenyl, R3= pentamethylphenyl) and with Sn[CH(PPh2)2]2 have been investigated. The tin reagents SnR2(R = R1 or R2) reacted under mild conditions to give in moderate yields the trinuclear species [Fe2(CO)8(µ-SnR12)]1 or [Fe2(CO)8(µ-SnR22)]2, as orange-red crystalline solids, which decompose in air on prolonged exposure. The compound [Fe2(CO)8(µ-SnR42)]3(R4= 2,4,6-triphenylphenyl) can be similarly obtained. Prolonged treatment of the carbonyl with the novel tin reagent SnR32, by contrast, afforded the known compound spiro-[(OC)8Fe2SnFe2(CO)8]4 for which data are briefly reported. Reactions with tin or lead reagents M[CH(PPh2)2]2(M = Sn or Pb) afforded [Fe2(CO)6(µ-CO)(µ-dppm)][dppm = 1,2-bis(diphenylphosphino)methane] rapidly and almost quantitatively. Full crystal and molecular structural data are reported for [Fe2(CO)8(µ-SnR12)] and [Fe2(CO)8(µ-SnR22)]. Mössbauer data are also presented for compounds 1–3, and interpreted in terms of the structural data for these and other systems.

The troposphere-to-stratosphere transition in kinetic energy spectra and nonlinear spectral fluxes as seen in ECMWF analyses

Global horizontal wavenumber kinetic energy spectra and spectral fluxes of rotational kinetic energy and enstrophy are computed for a range of vertical levels using a T799 ECMWF operational analysis. Above 250 hPa, the kinetic energy spectra exhibit a distinct break between steep and shallow spectral ranges, reminiscent of dual power-law spectra seen in aircraft data and high-resolution general circulation models. The break separates a large-scale ‘‘balanced’’ regime in which rotational flow strongly dominates divergent flow and a mesoscale ‘‘unbalanced’’ regime where divergent energy is comparable to or larger than rotational energy. Between 230 and 100 hPa, the spectral break shifts to larger scales (from n 5 60 to n 5 20, where n is spherical harmonic index) as the balanced component of the flow preferentially decays. The location of the break remains fairly stable throughout the stratosphere. The spectral break in the analysis occurs at somewhat larger scales than the break seen in aircraft data. Nonlinear spectral fluxes defined for the rotational component of the flow maximize between about 300 and 200 hPa. Large-scale turbulence thus centers on the extratropical tropopause region, within which there are two distinct mechanisms of upscale energy transfer: eddy–eddy interactions sourcing the transient energy peak in synoptic scales, and zonal mean–eddy interactions forcing the zonal flow. A well-defined downscale enstrophy flux is clearly evident at these altitudes. In the stratosphere, the transient energy peak moves to planetary scales and zonal mean–eddy interactions become dominant.

Comments on some recent measurements of anomalously steep N2O and O3tracer spectra in the stratospheric surf zone

Recent aircraft measurements, primarily in the extratropics, of the horizontal variance of nitrous oxide (N2O) and ozone (O3) in the middle stratosphere indicate that horizontal spectra of the tracer variance scale nearly as k−2, where k is the spatial wavenumber along the aircraft flight track [Strahan and Mahlman, 1994; Bacmeister et al., 1996]. This spectral scaling has been regarded as inconsistent with the accepted picture of stratospheric tracer motion; large-scale quasi-two-dimensional tracer advection typically yields a k−1 scaling (i.e., the classical Batchelor spectrum). In this paper it is argued that the nearly k−2 scaling seen in the measurements is a natural outcome of quasi-two-dimensional filamentation of the polar vortex edge. The accepted picture of stratospheric tracer motion can thus be retained: no additional physical processes are needed to account for deviations from the Batchelor spectrum. Our argument is based on the finite lifetime of tracer filaments and on the “singularity spectrum” associated with a one-dimensional field composed of randomly spaced jumps in concentration.

Differentiation between Neotropical rainforest, dry forest and savannah ecosystems by their modern pollen spectra and implications for the fossil pollen record

Accurate differentiation between tropical forest and savannah ecosystems in the fossil pollen record is hampered by the combination of: i) poor taxonomic resolution in pollen identification, and ii) the high species diversity of many lowland tropical families, i.e. with many different growth forms living in numerous environmental settings. These barriers to interpreting the fossil record hinder our understanding of the past distributions of different Neotropical ecosystems and consequently cloud our knowledge of past climatic, biodiversity and carbon storage patterns. Modern pollen studies facilitate an improved understanding of how ecosystems are represented by the pollen their plants produce and therefore aid interpretation of fossil pollen records. To understand how to differentiate ecosystems palynologically, it is essential that a consistent sampling method is used across ecosystems. However, to date, modern pollen studies from tropical South America have employed a variety of methodologies (e.g. pollen traps, moss polsters, soil samples). In this paper, we present the first modern pollen study from the Neotropics to examine the modern pollen rain from moist evergreen tropical forest (METF), semi-deciduous dry tropical forest (SDTF) and wooded savannah (cerradão) using a consistent sampling methodology (pollen traps). Pollen rain was sampled annually in September for the years 1999–2001 from within permanent vegetation study plots in, or near, the Noel Kempff Mercado National Park (NKMNP), Bolivia. Comparison of the modern pollen rain within these plots with detailed floristic inventories allowed estimates of the relative pollen productivity and dispersal for individual taxa to be made (% pollen/% vegetation or ‘p/v’). The applicability of these data to interpreting fossil records from lake sediments was then explored by comparison with pollen assemblages obtained from five lake surface samples.