Deciphering crystal packing effects in the spin crossover of six [FeII(2-pic)3]Cl2 solvatomorphs.

In this work we investigate how the crystal packing affects the spin crossover (SCO) of [FeII(2-pic)3]Cl2. Four alcohol solvatomorphs of this compound have been reported to present a transition from a singlet low spin (LS, S = 0) to a quintet high spin (HS, S = 5) state at different temperatures, whereas two other solvatomorphs remain in the HS state along the entire range of temperatures explored. Given that all these solvatomorphs are isostructural, this behavior hints at crystal packing effects playing an important role in the spin transition. With the aim of deciphering their origin, we performed periodic DFT+U+D2 computations on the crystal lattices of the six solvatomorphs to quantify the importance of all possible molecular and intermolecular contributions to their spin-state energetics. We demonstrate that the spin crossover of the alcohol solvatomorphs of [FeII(2-pic)3]Cl2 (1) is strongly influenced by the change of intermolecular interactions. We prove that isolated molecules of 1 would not undergo SCO, whereas intermolecular interactions are crucial to allow the spin transition in the solid state. The key interactions are those between (i) different SCO molecules (stabilizing the LS) and (ii) between SCO molecules and counterions (stabilizing the HS). Moreover, we show that these contributions are proportional to the crystal expansion associated with the phase change, in a way that larger volume expansions imply larger changes in these contributions, with the overall effect of stabilizing the LS state of the material. Finally, we discuss how our results challenge the common perspective that associates the cooperativity of a spin transition with the number- and strength- of intermolecular contacts in the crystal.

Cooperativity in Spin Crossover Systems. An Atomistic Perspective on the Devil's Staircase.

Cooperativity is key in defining the shape (i.e., gradual, abrupt, or hysteretic) of thermally driven spin transitions in magnetic switches. Despite its importance, there is very little information on its atomistic origin, which hinders the rational design of materials displaying a bistability region (i.e., hysteresis). In this paper, we investigate the spin transition of two solvatomorphs of [Fe(2-pic)3]Cl2, an Fe(II)-complex displaying thermal spin crossover (SCO) from a low-spin (LS) to a high-spin (HS) state with either gradual or abrupt two-step character. To do it, we apply a novel computational protocol to study the cooperativity of SCO compounds from DFT calculations, which does not rely on a priori assumptions on the studied system. The distinct shape of the spin transition is successfully captured, and the atomistic origin of cooperativity is traced back to geometrical distortions of the Fe-N6 core in case of the solvatomorph exhibiting an abrupt transition. According to our calculations, HS and LS molecules contribute differently to cooperativity, which results in a complex energetic evolution of the spin transition that cannot be described by the common Slichter-Drickamer model. The present work opens new avenues for the study of cooperativity of SCO systems having a chemically oriented perspective and demonstrates that quantum chemistry calculations can discriminate the shape of a spin transition, while providing insight into the atomistic underlying factors that contribute to its cooperative behavior.

Single molecule FRET investigation of pressure-driven unfolding of cold shock protein A.

We demonstrate that fused silica capillaries are suitable for single molecule fluorescence resonance energy transfer (smFRET) measurements at high pressure with an optical quality comparable to the measurement on microscope coverslips. Therefore, we optimized the imaging conditions in a standard square fused silica capillary with an adapted arrangement and evaluated the performance by imaging the focal volume, fluorescence correlation spectroscopy benchmarks, and FRET measurements. We demonstrate single molecule FRET measurements of cold shock protein A unfolding at a pressure up to 2000 bars and show that the unfolded state exhibits an expansion almost independent of pressure.

Structure-Based Insights into the Dynamics and Function of Two-Domain SlpA from Escherichia coli.

SlpA (SlyD-like protein A) comprises two domains, a FK506 binding domain (FKBP fold) of moderate prolyl cis/trans-isomerase activity and an inserted in flap (IF) domain that hosts its chaperone activity. Here we present the nuclear magnetic resonance (NMR) solution structure of apo Escherichia coli SlpA determined by NMR that mirrors the structural properties seen for various SlyD homologues. Crucial structural differences in side-chain orientation arise for F37, which points directly into the hydrophobic core of the active site. It forms a prominent aromatic stacking with F15, one of the key residues for PPIase activity, thus giving a possible explanation for the inherently low PPIase activity of SlpA. The IF domain reveals the highest stability within the FKBP-IF protein family, most likely arising from an aromatic cluster formed by four phenylalanine residues. Both the thermodynamic stability and the PPIase and chaperone activity let us speculate that SlpA is a backup system for homologous bacterial systems under unfavorable conditions.

Cooperativity of Spin Crossover Complexes: Combining Periodic Density Functional Calculations and Monte Carlo Simulation.

The total enthalpies of the 16 different spin configurations that can be realized in the unit cell of the archetype spin crossover complex [Fe(phen)2(NCS)2] (phen = 1,2-phenanthroline) were calculated, applying periodic density functional theory combined with the Hubbard model and the Grimme-D2 dispersion correction (DFT+U+D2). The obtained enthalpy differences between the individual spin configurations were used to determine spin couplings of an Ising-like model, and subsequent Monte Carlo simulations for this model allowed the estimation of the phenomenological interaction parameter Γ of the Slichter-Drickamer model, which is commonly used to describe the cooperativity of the spin transition. The calculation procedure described here-which led to an estimate of about 3 kJ·mol-1 for Γ, in good agreement with experiment-may be used to predict from first principles how modifications of spin crossover complexes can change the character of the spin transition from gradual to abrupt and vice versa.

The cysteine-rich region of type VII collagen is a cystine knot with a new topology.

Collagens are a group of extracellular matrix proteins with essential functions for skin integrity. Anchoring fibrils are made of type VII collagen (Col7) and link different skin layers together: the basal lamina and the underlying connective tissue. Col7 has a central collagenous domain and two noncollagenous domains located at the N and C terminus (NC1 and NC2), respectively. A cysteine-rich region of hitherto unknown function is located at the transition of the NC1 domain to the collagenous domain. A synthetic model peptide of this region was investigated by CD and NMR spectroscopy. The peptide folds into a collagen triple helix, and the cysteine residues form disulfide bridges between the different strands. The eight cystine knot topologies that are characterized by exclusively intermolecular disulfide bridges have been analyzed by molecular modeling. Two cystine knots are energetically preferred; however, all eight disulfide bridge arrangements are essentially possible. This novel cystine knot is present in type IX collagen, too. The conserved motif of the cystine knot is CX3CP. The cystine knot is N-terminal to the collagen triple helix in both collagens and therefore probably impedes unfolding of the collagen triple helix from the N terminus.

The spectroscopic ruler revisited at 77 K.

We report on single-molecule FRET measurements on the classical molecular ruler polyproline immersed in amorphous ice at 77 K. Confocal scanning microscopy allows for observation of the fluorescence of the dyes Alexa Fluor 488 as the donor and Alexa Fluor 594 as the acceptor in the FRET experiment. The dipole orientations are fixed in amorphous ice leading to a broad distribution of transfer efficiencies. The measured transfer efficiency distributions are mainly in accordance with simulated distributions, with a significant deviation at low transfer efficiencies, which are absent in the experimental distribution, in stark contrast to the theory.

Nuclear inelastic scattering and Mössbauer spectroscopy as local probes for ligand binding modes and electronic properties in proteins: vibrational behavior of a ferriheme center inside a ß-barrel protein.

In this work, we present a study of the influence of the protein matrix on its ability to tune the binding of small ligands such as NO, cyanide (CN(-)), and histamine to the ferric heme iron center in the NO-storage and -transport protein Nitrophorin 2 (NP2) from the salivary glands of the blood-sucking insect Rhodnius prolixus. Conventional Mössbauer spectroscopy shows a diamagnetic ground state of the NP2-NO complex and Type I and II electronic ground states of the NP2-CN(-) and NP2-histamine complex, respectively. The change in the vibrational signature of the protein upon ligand binding has been monitored by Nuclear Inelastic Scattering (NIS), also called Nuclear Resonant Vibrational Spectroscopy (NRVS). The NIS data thus obtained have also been calculated by quantum mechanical (QM) density functional theory (DFT) coupled with molecular mechanics (MM) methods. The calculations presented here show that the heme ruffling in NP2 is a consequence of the interaction with the protein matrix. Structure optimizations of the heme and its ligands with DFT retain the characteristic saddling and ruffling only if the protein matrix is taken into account. Furthermore, simulations of the NIS data by QM/MM calculations suggest that the pH dependence of the binding of NO, but not of CN(-) and histamine, might be a consequence of the protonation state of the heme carboxyls.

Exonucleolytic degradation of high-density labeled DNA studied by fluorescence correlation spectroscopy.

The exonucleolytic degradation of high-density labeled DNA by exonuclease III was monitored using two-color fluorescence correlation spectroscopy (FCS). One strand of the double stranded template DNA was labeled on either one or two base types and additionally at one end via a 5' Cy5 tagged primer. Exonucleolytic degradation was followed via the diffusion time, the brightness of the remaining DNA as well as the concentration of released labeled bases. We found a hydrolyzation rate of about 11 to 17 nucleotides per minute per enzyme (nt/min/enzyme) for high-density labeled DNA, which is by a factor of about 4 slower than for unlabeled DNA. The exonucleolytic degradation of a 488 base pair long double stranded DNA resulted in a short double stranded DNA segment of 112 ± 40 base pairs (bp) length with two single-stranded tails.

Interpretation of nuclear resonant vibrational spectra of rubredoxin using a combined quantum mechanics and molecular mechanics approach.

Nuclear resonant vibrational spectra of the reduced and oxidized form of a mutant of rubredoxin from Pyrococcus abyssii were measured and are compared with simulated spectra that were calculated by a combined quantum mechanics (QM) and molecular mechanics (MM) method. Density functional theory was used for the QM level. Calculations were performed for different models of rubredoxin. Realistic spectra were simulated with reduced models that include at least the iron center, the four cysteins coordinating it, and the residues connected to the cysteins together with a QM layer that comprises the first two coordination shells of the iron center. Larger QM layers did not lead to significant changes of the simulated spectra.

Fe(II) complex with the octadentate btpa ligand: a DFT study on a spin-crossover system that reveals two distinct high-spin states.

Density functional theory calculations (DFT) were performed for the spin-crossover system [Fe(btpa)]2+ (btpa = N,N,N',N'-tetrakis(2-pyridylmethyl)-6,6'-bis(aminomethyl)-2,2'-bipyridine), and for the predominantly low-spin [Fe(b(bdpa))]2+ complex (in the solid state) (b(bdpa) = N,N'-bis(benzyl)-N,N'-bis(2-pyridylmethyl)-6,6'-bis(aminomethyl)-2,2'-bipyridine). The calculations confirmed that the former complex exhibits two high-spin isomers of the complexes, i.e. with C1 quasi hepta-coordinated (long-lived isomer) and C2 hexa-coordinated (short-lived isomer) structures that have been suggested previously based on time-resolved Raman and flash photolysis experiments. Application of B3LYP and B3LYP* functionals together with the CEP-31G basis yielded reasonable estimates of electronic energies (E(el) = E(el)(HS) - E(el)(LS)) for both isomers (calculated E(el) of ca. 24 and 31 kJ mol(-1) for long- and short-lived HS isomers, respectively, vs. the experimentally determined value of 27.5 kJ mol(-1)). Further calculations yielded the electronic structure of the low-spin isomer together with lowest lying singlet and triplet excited states of the [Fe(btpa)]2+ as well as the energy profile of the C2 <--> C1 isomerisation pathway for the high-spin [Fe(btpa)]2+ within the framework of the QST (quadratic synchronous transit) approach. The data obtained are discussed in relation to the observed ultrafast intersystem crossing in Fe(II) polypyridine complexes. The importance of ligand strain in relation to the destabilisation of the low-spin isomers is also discussed. In that context, calculations for a further 15 Fe(II) spin-crossover complexes of hexa-coordinating nitrogen-donor ligands have shown that the LS-HS conversion is associated with a release of ligand stress of 95 +/-16 kJ mol(-1), on average. On the basis of the calculations presented in this paper we propose that octahedral high-spin d-6 isomers are far more elastic regarding the angular distortions (equatorial and meridional strain, i.e. the declination of cis- and trans-L-M-L angle from the regular values of 90 and 180 degrees, respectively) than their low-spin counterparts.

Sub-picosecond time resolved infrared spectroscopy of high-spin state formation in Fe(ii) spin crossover complexes.

The photoinduced low-spin (S = 0) to high-spin (S = 2) transition of the iron(ii) spin-crossover systems [Fe(btpa)](PF(6))(2) and [Fe(b(bdpa))](PF(6))(2) in solution have been studied for the first time by means of ultrafast transient infrared spectroscopy at room temperature. Negative and positive infrared difference bands between 1000 and 1065 cm(-1) that appear within the instrumental system response time of 350 fs after excitation at 387 nm display the formation of the vibrationally unrelaxed and hot high-spin (5)T(2) state. Vibrational relaxation is observed and characterized by the time constants 9.4 +/- 0.7 ps for [Fe(btpa)](PF(6))(2)/acetone and 12.7 +/- 0.7 ps for both [Fe(btpa)](PF(6))(2)/acetonitrile and [Fe(b(bdpa)](PF(6))(2)/acetonitrile. Vibrational analysis has been performed via DFT calculations of the low-spin and high-spin state normal modes of both compounds as well as their respective infrared absorption cross sections. The simulated infrared difference spectra are dominated by an increase of the absorption cross section upon high-spin state formation in accordance with the experimental infrared spectra.

Estimate of the vibrational contribution to the entropy change associated with the spin transition in the d4 systems [MnIII(pyrol)3tren] and [CrII(depe)2I2].

The vibrational contribution to DeltaS of the low-spin ((3)T(1)) to high-spin ((5)E) spin transition in two 3d(4) octahedral systems [Mn(III)(pyrol)(3)tren] and [Cr(depe)(2)I(2)] have been estimated by means of DFT calculations (B3LYP/CEP-31G) of the vibrational normal-modes frequencies. The obtained value at the transition temperature for the Mn(iii) complex is DeltaS(vib)(44 K) = 6.3 J K(-1) mol(-1), which is comparable with the proposed Jahn-Teller contribution of R ln3 = 9.1 J K(-1) mol(-1) and which is approximately half of the experimentally determined 13.8 J K(-1) mol(-1). The corresponding value for the Cr(ii) complex is DeltaS(vib)(171.45 K) = 46.5 J K(-1) mol(-1), as compared to the experimental value of 39.45 J K(-1) mol(-1). The analysis of the vibrational normal modes reveals that for the d(4) systems under study, contrary to Fe(ii) d(6) systems, not all metal-ligand stretching vibrations make a contribution. For the Mn(iii) complex, the only vibration that contributes to DeltaS(vib) involve the nitrogens occupying the Jahn-Teller axis, while in the case of Cr(ii) the contributing vibrations involve the Cr-I bonds. Low-frequency modes due to ring vibrations, metal-ligand bending and movement of the molecule as a whole also contribute to the vibrational entropy associated with the spin transition.

Vibrational spectrum of the spin crossover complex [Fe(phen)(2)(NCS)(2)] studied by IR and Raman spectroscopy, nuclear inelastic scattering and DFT calculations.

The vibrational modes of the low-spin and high-spin isomers of the spin crossover complex [Fe(phen)(2)(NCS)(2)] (phen = 1,10-phenanthroline) have been measured by IR and Raman spectroscopy and by nuclear inelastic scattering. The vibrational frequencies and normal modes and the IR and Raman intensities have been calculated by density functional methods. The vibrational entropy difference between the two isomers, DeltaS(vib), which is--together with the electronic entropy difference DeltaS(el)--the driving force for the spin-transition, has been determined from the measured and from the calculated frequencies. The calculated difference (DeltaS(vib) = 57-70 J mol(-1) K(-1), depending on the method) is in qualitative agreement with experimental values (20-36 J mol(-1) K(-1)). Only the low energy vibrational modes (20% of the 147 modes of the free molecule) contribute to the entropy difference and about three quarters of the vibrational entropy difference are due to the 15 modes of the central FeN(6) octahedron.

A range of spin-crossover temperature T1/2>300 K results from out-of-sphere anion exchange in a series of ferrous materials based on the 4-(4-imidazolylmethyl)-2-(2-imidazolylmethyl)imidazole (trim) ligand, [Fe(trim)2]X2 (X=F, Cl, Br, I): comparison of experimental results with those derived from density functional theory calculations.

The synthesis and characterization of [FeII(trim)2]Cl2 (2), [FeII(trim)2]Br2MeOH (3), and [FeII(trim)2]I2MeOH (4), including the X-ray crystal structure determinations of 2 (50 and 293 K) and 4 (293 K), have been performed and their properties have been examined. In agreement with the magnetic susceptibility results, the Mössbauer data show the presence of high-spin (HS) to low-spin (LS) crossover with a range of T1/2 larger than 300 K (from approximately 20 K for [FeII(trim)2]F2 (1) to approximately 380 K for 4). All complexes in this series include the same [Fe(trim)2]2+ complex cation: the ligand field comprises a constant contribution from the trim ligands and a variable one originating from the out-of-sphere anions, which is transmitted to the metal center by the connecting imidazole rings and hydrogen bonds. The impressive variation in the intrinsic characteristics of the spin-crossover (SCO) phenomenon in this series is then interpreted as an inductive effect of the anions transmitted to the nitrogen donors through the hydrogen bonds. Based on this qualitative analysis, an increased inductive effect of the out-of-sphere anion corresponds to a decreased SCO temperature T1/2, in agreement with the experimental results. Electronic structure calculations with periodic boundary conditions have been performed that show the importance of intermolecular effects in tuning the ligand field, and thus in determining the transition temperature. Starting with the geometries obtained from the X-ray studies, the [FeII(trim)2]X2 complex molecules 1-4 have been investigated both for the single molecules and the crystal lattices with the local density approximation of density functional theory. The bulk geometries of the complex cations deduced from the X-ray studies and those calculated are in fair agreement for both approaches. However, the trend observed for the transition temperatures of 1-4 disagrees with the trend for the spin-state splittings ES (difference EHS-ELS between the energy of the HS and LS isomers) calculated for the isolated molecules, whereas it agrees with the trend for ES calculated with periodic boundary conditions. The latter calculations predict the strongest stabilization of the HS state for the fluoride complex, which actually is essentially HS above T=50 K, while the most pronounced stabilization of the LS state is predicted for 4, in line with the experimental results.

Models of the membrane-bound cytochromes: mössbauer spectra of crystalline low-spin ferriheme complexes having axial ligand plane dihedral angles ranging from 0 degree to 90 degrees.

Crystalline samples of four low-spin Fe(III) octaalkyltetraphenylporphyrinate and two low-spin Fe(III) tetramesitylporphyrinate complexes, all of which are models of the bis-histidine-coordinated cytochromes of mitochondrial complexes II, III, and IV and chloroplast complex b(6)f, and whose molecular structures and EPR spectra have been reported previously, have been investigated in detail by Mössbauer spectroscopy. The six complexes and the dihedral angles between axial ligand planes of each are [(TMP)Fe(1-MeIm)(2)]ClO(4) (0 degree), paral-[(OMTPP)Fe(1-MeIm)(2)]Cl (19.5 degrees), paral-[(TMP)Fe(5-MeHIm)(2)]ClO(4) (26 degrees, 30 degrees for two molecules in the unit cell whose EPR spectra overlap), [(OETPP)Fe(4-Me(2)NPy)(2)]Cl (70 degrees), perp-[(OETPP)Fe(1-MeIm)(2)]Cl (73 degrees), and perp-[(OMTPP)Fe(1-MeIm)(2)]Cl (90 degrees). Of these, the first three have been shown to exhibit normal rhombic EPR spectra, each with three clearly resolved g-values, while the last three have been shown to exhibit "large g(max)" EPR spectra at 4.2 K. It is found that the hyperfine coupling constants of the complexes are consistent with those reported previously for low-spin ferriheme systems, with the largest-magnitude hyperfine coupling constant, A(zz), being considerably smaller for the "parallel" complexes (400-540 kG) than for the strictly perpendicular complex (902 kG), A(xx) being negative for all six complexes, and A(zz) and A(xx) being of similar magnitude for the "parallel" complexes (for example, for [(TMP)Fe(1-MeIm)(2)]Cl, A(zz) = 400 kG, A(xx) = -400 kG). In all cases, A(yy) is small but difficult to estimate with accuracy. With results for six structurally characterized model systems, we find for the first time qualitative correlations of g(zz), A(zz), and DeltaE(Q) with axial ligand plane dihedral angle Deltavarphi.

High-spin to low-spin relaxation kinetics in the [Fe(TRIM)2]Cl2 complex.

A quasi-quantitative photo-induced low-spin (LS)-->high-spin (HS) conversion of FeII ions has been observed in the [Fe(TRIM)2]Cl2 complex by irradiating the sample with blue light (488 nm) at 10 K. The time dependence of the HS-->LS relaxation has been studied between 10 K and 44 K by means of magnetic susceptibility measurements. These relaxation curves could be satisfactorily fitted by mono-exponential decays including tunnelling effect except for temperatures below 30 K. The introduction of a distribution of vibrational frequencies into this model improved significantly the fits in the low-temperature range and gave a good agreement with the experimental data in the whole temperature range suggesting a multi-rate relaxation process in this complex.

Strain-induced substitutional lability in a Ru(II) complex of a hypodentate polypyridine ligand.

The ruthenium(II) complex of heptadentate N,N,N',N'-tetrakis(2-pyridylmethyl)-2,6-bis(aminomethyl)pyridine (tpap) was isolated as the hexafluorophosphate complex Ru(tpap)(PF6)2. The crystal structure has been determined for the triflate salt Ru(tpap)(CF3SO3)2.2H2O, which crystallizes in the monoclinic space group P2(1)/n. The structure was refined to a final R value of 0.0549 for 5894 observed reflections. The heptadentate ligand coordinates with six nitrogens, i.e. with two tertiary nitrogens and four pyridine nitrogens, one of the pyridines remaining un-coordinated. The resulting structure is significantly distorted from octahedral geometry with an equatorial Nsp3-Ru-Npyridine angle of 120 degrees. The consequence of the above steric strain is a labilization of the system and fluxional behaviour involving exchange between equatorially coordinated and non-coordinated pyridines has been observed by 1H NMR for Ru(tpap)(PF6)2 in d6-acetone solution. The activation parameters of DeltaG(not equal to 298) = 53 kJ mol(-1), DeltaH(not equal) = 56 +/- 1 kJ mol(-1) and DeltaS(not equal) = -10 +/- 3 J mol(-1) K(-1) were determined on the basis of NMR experiments. In addition electronic structure calculations applying density functional theory (DFT) have been performed in order to identify a transition state and to estimate the activation barrier. On the basis of NMR and DFT results the mechanism of isoexchange involving a hepta-coordinated intermediate has been proposed.

Metastable isonitrosyl structure of the nitroprusside anion confirmed by nuclear inelastic scattering.

Nuclear inelastic scattering (NIS) measurements were performed on a guanidium nitroprusside ((CN(3)H(6))(2)[Fe(CN)(5)NO], GNP) monocrystal at 77 K after the sample was illuminated with blue light (450 nm) at 50 K to populate the two metastable states, MS(1) and MS(2), of the nitroprusside anion. A second measurement was performed at 77 K after warming up the illuminated crystal to 250 K where the metastable states decay to the groundstate. The measured spectra were compared with simulated NIS spectra that were calculated by using density functional methods. Comparison of measured and simulated spectra provides strong evidence for the isonitrosyl structure of the metastable MS(1) state proposed by Carducci et al. (Carducci, M. D.; Pressprich, M. R.; Coppens, P. J. Am. Chem. Soc. 1997, 119, 2669-2678).