Hydration dynamics and IR spectroscopy of 4-fluorophenol.

Halogenated groups are relevant in pharmaceutical applications and potentially useful spectroscopic probes for infrared spectroscopy. In this work, the structural dynamics and infrared spectroscopy of para-fluorophenol (F-PhOH) and phenol (PhOH) is investigated in the gas phase and in water using a combination of experiment and molecular dynamics (MD) simulations. The gas phase and solvent dynamics around F-PhOH and PhOH is characterized from atomistic simulations using empirical energy functions with point charges or multipoles for the electrostatics, Machine Learning (ML) based parametrizations and with full ab initio (QM) and mixed Quantum Mechanical/Molecular Mechanics (QM/MM) simulations with a particular focus on the CF- and OH-stretch region. The CF-stretch band is heavily mixed with other modes whereas the OH-stretch in solution displays a characteristic high-frequency peak around 3600 cm-1 most likely associated with the -OH group of PhOH and F-PhOH together with a characteristic progression below 3000 cm-1 due to coupling with water modes which is also reproduced by several of the simulations. Solvent and radial distribution functions indicate that the CF-site is largely hydrophobic except for simulations using point charges which renders them unsuited for correctly describing hydration and dynamics around fluorinated sites. The hydrophobic character of the CF-group is particularly relevant for applications in pharmaceutical chemistry with a focus on local hydration and interaction with the surrounding protein.

Folding and Unfolding of the Tryptophan Zipper in the Presence of Two Thioamide Substitutions.

We studied the stability and folding and unfolding kinetics of the tryptophan zipper, containing different double thioamide subsitutions. Conformation change was triggered by photoisomerization of an integrated AMPP photoswitch in the turn region of the hairpin, and transient spectra were recorded in the deep UV and the mid-IR, covering the time window of the (un)folding transition from picoseconds to tens of microseconds. Thio-substitution of inward-pointing backbone carbonyls was found to strongly destabilize the ß-hairpin structures, whereas molecules with two outward pointing thio-carbonyls showed similar or enhanced stability with respect to the unsubstituted sequence, which we attribute to stronger interstrand hydrogen bonding. Thiolation of the two Trp residues closest to the turn can even prevent the opening of the hairpin after cis-trans isomerization of the switch. The circular dichroism due to the two thioamide ππ* transitions is spectrally well-separated from the aromatic tryptophan signal. It changes upon photoswitching, reflecting a local change in coupling and geometry.

Broad-Band Ultraviolet CD Spectroscopy of Ultrafast Peptide Backbone Conformational Dynamics.

The far-UV spectral window widely used for the conformational analysis of biomolecules is not easily covered with broad-band lasers. This has made it difficult to use circular dichroism (CD) spectroscopy to directly follow fast structure changes. By combining transient CD spectroscopy in the deep-UV with thioamide substitution, we demonstrate a method to overcome this difficulty. We investigated a dipeptide whose two carbonyl oxygen atoms were replaced by sulfur, red-shifting the strong lowest-lying ππ* transitions into the more accessible 250-370 nm spectral window. Coupling of the two thioamide units cannot be resolved by achiral 2D-UV spectroscopy, but it gives rise to a pronounced bisignate CD spectrum. The transient CD spectra reveal weakening of this coupling in the electronically excited state, where conformational constraints are released. Our results show that direct local probing of fast backbone conformational change via CD spectroscopy is possible in combination with site-selective thio substitution in peptides and proteins.

Kinetic response of a photoperturbed allosteric protein.

By covalently linking an azobenzene photoswitch across the binding groove of a PDZ domain, a conformational transition, similar to the one occurring upon ligand binding to the unmodified domain, can be initiated on a picosecond timescale by a laser pulse. The protein structures have been characterized in the two photoswitch states through NMR spectroscopy and the transition between them through ultrafast IR spectroscopy and molecular dynamics simulations. The binding groove opens on a 100-ns timescale in a highly nonexponential manner, and the molecular dynamics simulations suggest that the process is governed by the rearrangement of the water network on the protein surface. We propose this rearrangement of the water network to be another possible mechanism of allostery.

Photocontrol of reversible amyloid formation with a minimal-design peptide.

Amyloid aggregates are highly ordered fibrillar assemblies of polypeptides involved in a number of neurodegenerative diseases. Very little is known on the pathways of self-assembly of peptides into the final amyloid fibrils, which is due in part to the difficulty of triggering the aggregation process in a controlled manner. Here we present the design and validation of a cross-linked hexapeptide that reversibly aggregates and dissociates under ultraviolet light irradiation control. First molecular dynamics simulations were carried out to identify, among hundreds of possible sequences, those with the highest propensity to form ordered (ß-sheet) oligomers in the trans state of the azobenzene cross-linker, and at the same time with the highest solubility in the cis state. In the simulations, the peptides were observed to spontaneously form ordered oligomers with cross-ß contacts when the cross-linker was in the trans state, whereas in the cis state they self-assemble into amorphous aggregates. For the most promising sequence emerging from the simulations (Ac-Cys-His-Gly-Gln-Cys-Lys-NH(2) cross-linked at the two cysteine residues), the photoisomerization of the azobenzene group was shown to induce reversible aggregation by time-resolved light scattering and fluorescence measurements. The amyloid-like fibrillar topology was confirmed by electron microscopy. Potential applications of minimally designed peptides with photoswitchable amyloidogenic propensity are briefly discussed.

2D-IR study of a photoswitchable isotope-labeled alpha-helix.

A series of photoswitchable, alpha-helical peptides were studied using two-dimensional infrared spectroscopy (2D-IR). Single-isotope labeling with (13)C(18)O at various positions in the sequence was employed to spectrally isolate particular backbone positions. We show that a single (13)C(18)O label can give rise to two bands along the diagonal of the 2D-IR spectrum, one of which is from an amide group that is hydrogen-bonded internally, or to a solvent molecule, and the other from a non-hydrogen-bonded amide group. The photoswitch enabled examination of both the folded and unfolded state of the helix. For most sites, unfolding of the peptide caused a shift of intensity from the hydrogen-bonded peak to the non-hydrogen-bonded peak. The relative intensity of the two diagonal peaks gives an indication of the fraction of molecules hydrogen-bonded at a certain location along the sequence. As this fraction varies quite substantially along the helix, we conclude that the helix is not uniformly folded. Furthermore, the shift in hydrogen bonding is much smaller than the change of helicity measured by CD spectroscopy, indicating that non-native hydrogen-bonded or mis-folded loops are formed in the unfolded ensemble.

Dynamical transition in a small helical peptide and its implication for vibrational energy transport.

The two-dimensional infrared spectrum of an octameric helical peptide in chloroform was measured as a function of temperature. Isotope labeling of the carbonyl group of one of the amino acids was used to obtain information for an isolated vibration. The antidiagonal width of the 2D-IR signal, which is a measure of the homogeneous dephasing time T(2), is constant from 220 to 260 K (within experimental error), and increases steeply above. The homogeneous dephasing time of the carbonyl vibration is attributed to the flexibility of the system and/or its immediate surrounding. The system undergoes a dynamical transition at about 270 K, with similarities to the protein dynamical transition. Furthermore, the temperature dependence of the antidiagonal width strongly resembles that of the efficiency of vibrational energy transport along the helix, which has been studied in a recent paper (J. Phys. Chem. B 2008, 112, 15487). The connection between the two processes, structural flexibility and energy transport mechanism, is discussed.

Synthesis, characterization and applicability of three isotope labeled azobenzene photoswitches.

We describe a short, efficient approach for the synthesis of three novel isotope labeled azobenzene photoswitches. The synthesis is based on commercially available fully isotope labeled precursors. The target molecules have been obtained in good yields, checked for purity, and identified by NMR and IR spectroscopy and a variety of standard analytical methods (UV-vis, mp, ESI-MS, elemental analysis). Using conventional coupling techniques the three isotope labeled photoswitches can be incorporated very easily in biomacromolecules.

Energy transport in peptide helices: a comparison between high- and low-energy excitations.

Energy transport in a short helical peptide in chloroform solution is studied by time-resolved femtosecond spectroscopy and accompanying nonequilibrium molecular dynamics (MD) simulations. In particular, the heat transport after excitation of an azobenzene chromophore attached to one terminus of the helix with 3 eV (UV) photons is compared to the excitation of a peptide C=O oscillator with 0.2 eV (IR) photons. The heat in the helix is detected at various distances from the heat source as a function of time by employing vibrational pump-probe spectroscopy. As a result, the carbonyl oscillators at different positions along the helix act as local thermometers. The experiments show that heat transport through the peptide after excitation with low-energy photons is at least 4 times faster than after UV excitation. On the other hand, the heat transport obtained by nonequilibrium MD simulations is largely insensitive to the kind of excitation. The calculations agree well with the experimental results for the low-frequency case; however, they give a factor of 5 too fast energy transport for the high-energy case. Employing instantaneous normal mode calculations of the MD trajectories, a simple harmonic model of heat transport is adopted, which shows that the heat diffusivity decreases significantly at temperatures initially reached by high-energy excitation. This finding suggests that the photoinduced energy gets trapped, if it is deposited in high amounts. The various competing mechanisms, such as vibrational T(1) relaxation, resonant transfer between excitonic states, cascading down relaxation, and low-frequency mode transfer, are discussed in detail.

Time-resolved infrared spectroscopy of thiopeptide isomerization and hydrogen-bond breaking.

The photoisomerization of the protected tetrathioxopeptide Boc-Ala-Gly(=S)-Ala-Aib-OMe was followed using time-resolved infrared spectroscopy in the amide I region in combination with isotope labeling. In acetonitrile at room temperature, approximately half of the molecules are found in a loop conformation, restrained by an intramolecular hydrogen bond, while the other half adopts more extended conformations. UV-excitation of the thioxopeptide unit immediately weakens the intramolecular hydrogen bond. After the molecules have relaxed to the electronic ground state with a 130 ps time-constant, a delayed re-formation of the intramolecular hydrogen bond is observed for molecules returning to the initial trans conformation of the thioamide bond, while the loop structure is permanently broken when the molecules isomerize to the cis conformation.

Coexistence of hydrogen-bonded loop and extended tetrapeptide conformations.

The conformations of a protected tetrathiopeptide have been analysed by isotope labelling and two-dimensional infrared spectroscopy (2D-IR). It has been found that Boc-Ala-Gly(=S)-Ala-Aib-OMe in acetonitrile, as well as its oxopeptide analogue, can adopt a hydrogen-bonded loop conformation in coexistence with less restricted structures. The two types of conformations interconvert too quickly to be resolved on the nuclear magnetic resonance (NMR) timescale, but give rise to different cross peaks in two-dimensional infrared spectra. The hydrogen bond between the Boc terminal group and the amide proton of Aib can be broken by photoisomerisation of the thioamide bond.

Energy transport in peptide helices.

We investigate energy transport through an alpha-aminoisobutyric acid-based 3(10)-helix dissolved in chloroform in a combined experimental-theoretical approach. Vibrational energy is locally deposited at the N terminus of the helix by ultrafast internal conversion of a covalently attached, electronically excited, azobenzene moiety. Heat flow through the helix is detected with subpicosecond time resolution by employing vibrational probes as local thermo meters at various distances from the heat source. The experiment is supplemented by detailed nonequilibrium molecular dynamics (MD) simulations of the process, revealing good qualitative agreement with experiment: Both theory and experiment exhibit an almost instantaneous temperature jump of the reporter units next to the heater which is attributed to the direct impact of the isomerizing azobenzene moiety. After this impact event, helix and azobenzene moiety appear to be thermally decoupled. The energy deposited in the helix thermalizes on a subpicosecond timescale and propagates along the helix in a diffusive-like process, accompanied by a significant loss into the solvent. However, in terms of quantitative numbers, theory and experiment differ. In particular, the MD simulation seems to overestimate the heat diffusion constant (2 A(2) ps(-1) from the experiment) by a factor of five.

Folding and unfolding of a photoswitchable peptide from picoseconds to microseconds.

Using time-resolved IR spectroscopy, we monitored the kinetics of folding and unfolding processes of a photoswitchable 16-residue alanine-based alpha-helical peptide on a timescale from few picoseconds to almost 40 micros and over a large temperature range (279-318 K). The folding and unfolding processes were triggered by an ultrafast laser pulse that isomerized the cross linker within a few picoseconds. The main folding and unfolding times (700 ns and 150 ns, respectively, at room temperature) are in line with previous T-jump experiments obtained from similar peptides. However, both processes show complex, strongly temperature-dependent spectral kinetics that deviate clearly from a single-exponential behavior. Whereas in the unfolding experiment the ensemble starts from a well defined folded state, the starting ensemble in the folding experiment is more heterogeneous, which leads to distinctly different kinetics of the experiments, because they are sensitive to different regions of the energy surface. A qualitative agreement with the experimental data-set can be obtained by a model where the unfolded states act as a hub connected to several separated "misfolded" states with a distribution of rates. We conclude that a rather large spread of rates (k(1) : k(n) approximately 9) is needed to explain the experimentally observed stretched exponential response with stretching factor beta = 0.8 at 279 K.

A fast photoswitch for minimally perturbed peptides: investigation of the trans-->cis photoisomerization of N-methylthioacetamide.

Thio amino acids can be integrated into the backbone of peptides without significantly perturbing their structure. In this contribution we use ultrafast infrared and visible spectroscopy as well as state-of-the-art ab initio computations to investigate the photoisomerization of the trans form of N-methylthioacetamide (NMTAA) as a model conformational photoswitch. Following the S2 excitation of trans-NMTAA in water, the return of the molecule into the trans ground state and the formation of the cis isomer is observed on a dual time scale, with a fast component of 8-9 ps and a slow time constant of approximately 250 ps. On both time scales the probability of isomerization to the cis form is found to be 30-40%, independently of excitation wavelength. Ab initio CASPT2//CASSCF photochemical reaction path calculations indicate that, in vacuo, the trans-->cis isomerization event takes place on the S1 and/or T1 triplet potential energy surfaces and is controlled by very small energy barriers, in agreement with the experimentally observed picosecond time scale. Furthermore, the calculations identify one S2/S1 and four nearly isoenergetic S1/S0 conical intersection decay channels. In line with the observed isomerization probability, only one of the S1/S0 conical intersections yields the cis conformation upon S1-->S0 decay. A substantially equivalent excited-state relaxation results from four T1/S0 intersystem crossing points.

Occurrence and fate modeling of MTBE and BTEX compounds in a Swiss Lake used as drinking water supply.

Emissions of fuel components from boating use on multiple-use lakes and reservoirs are of high concern with regard to the drinking water supply from such water bodies. We report results of a detailed study on the occurrence, sources and fate of aromatic hydrocarbons and methyl tert-butyl ether (MTBE) in a typical holomictic lake, Lake Zurich, that supplies drinking water for the largest Swiss city. Emphasis of the investigation was on the fuel oxygenate MTBE, which was found in concentrations up to 1.4 microg/L in the epilimnion and up to 0.05microg/L in the hypolimnion of the lake. The concentration difference was due to the stratification of the lake during the boating season with very limited water exchange across the thermocline. MTBE and BTEX nearly completely volatilized before vertical lake mixing occurred in winter. Spatial and temporal variations of MTBE concentrations in the lake were observed and successfully predicted using two complementary box models (MASAS Light and Aquasim). The drinking water supply from holomictic lakes is not at risk for the scenarios studied if water is extracted from well below the thermocline. Since emissions of unburned gasoline into such water bodies are caused predominantly by boating activities, restrictions of highly emitting two-stroke engines could substantially reduce the MTBE and BTEX load of the epilimnion during the boating season.