Solvation Dynamics in Protein Environments Studied by Photon Echo Spectroscopy

Wide-Dynamic-Range Control of Quantum-Electrodynamic Electron Transfer Reactions in the Weak Coupling Regime

Catalyzing reactions effectively by vacuum fluctuations of electromagnetic fields is a significant challenge within the realm of chemistry. As opposed to most studies based on vibrational strong coupling, we introduce an innovative catalytic mechanism driven by weakly coupled polaritonic fields. Through the amalgamation of macroscopic quantum electrodynamics (QED) principles with Marcus electron transfer (ET) theory, we predict that ET reaction rates can be precisely modulated across a wide dynamic range by controlling the size and structure of nanocavities. Compared to QED-driven radiative ET rates in free space, plasmonic cavities induce substantial rate enhancements spanning the range from 103- to 10-fold. By contrast, Fabry-Perot cavities engender rate suppression spanning the range from 10-2- to 10-1-fold. This work overcomes the necessity of using strong light-matter interactions in QED chemistry, opening up a new era of manipulating QED-based chemical reactions in a wide dynamic range.

Directed Evolution of a Bright Variant of mCherry: Suppression of Nonradiative Decay by Fluorescence Lifetime Selections

The approximately linear scaling of fluorescence quantum yield (ϕ) with fluorescence lifetime (τ) in fluorescent proteins (FPs) has inspired engineering of brighter fluorophores based on screening for increased lifetimes. Several recently developed FPs such as mTurquoise2, mScarlet, and FusionRed-MQV which have become useful for live cell imaging are products of lifetime selection strategies. However, the underlying photophysical basis of the improved brightness has not been scrutinized. In this study, we focused on understanding the outcome of lifetime-based directed evolution of mCherry, which is a popular red-FP (RFP). We identified four positions (W143, I161, Q163, and I197) near the FP chromophore that can be mutated to create mCherry-XL (eXtended Lifetime: ϕ = 0.70; τ = 3.9 ns). The 3-fold higher quantum yield of mCherry-XL is on par with that of the brightest RFP to date, mScarlet. We examined selected variants within the evolution trajectory and found a near-linear scaling of lifetime with quantum yield and consistent blue-shifts of the absorption and emission spectra. We find that the improvement in brightness is primarily due to a decrease in the nonradiative decay of the excited state. In addition, our analysis revealed the decrease in nonradiative rate is not limited to the blue-shift of the energy gap and changes in the excited state reorganization energy. Our findings suggest that nonradiative mechanisms beyond the scope of energy-gap models such the Englman-Jortner model are suppressed in this lifetime evolution trajectory.

Destructive Photon Echo Formation in Six-Wave Mixing Signals of a MoSe2 Monolayer

Monolayers of transition metal dichalcogenides display a strong excitonic optical response. Additionally encapsulating the monolayer with hexagonal boron nitride allows to reach the limit of a purely homogeneously broadened exciton system. On such a MoSe2 -based system, ultrafast six-wave mixing spectroscopy is performed and a novel destructive photon echo effect is found. This process manifests as a characteristic depression of the nonlinear signal dynamics when scanning the delay between the applied laser pulses. By theoretically describing the process within a local field model, an excellent agreement with the experiment is reached. An effective Bloch vector representation is developed and thereby it is demonstrated that the destructive photon echo stems from a destructive interference of successive repetitions of the heterodyning experiment.

Exciton Delocalization Counteracts the Energy Gap: A New Pathway toward NIR-Emissive Dyes

Exciton coupling between the transition dipole moments of ordered dyes in supramolecular assemblies, so-called J/H-aggregates, leads to shifted electronic transitions. This can lower the excited state energy, allowing for emission well into the near-infrared regime. However, as we show here, it is not only the excited state energy modifications that J-aggregates can provide. A bay-alkylated quaterrylene was synthesized, which was found to form J-aggregates in 1,1,2,2-tetrachloroethane. A combination of superradiance and a decreased nonradiative relaxation rate made the J-aggregate four times more emissive than the monomeric counterpart. A reduced nonradiative relaxation rate is a nonintuitive consequence following the 180 nm (3300 cm^-1) red-shift of the J-aggregate in comparison to the monomeric absorption. However, the energy gap law, which is commonly invoked to rationalize increased nonradiative relaxation rates with increasing emission wavelength, also contains a reorganization energy term. The reorganization energy is highly suppressed in J-aggregates due to exciton delocalization, and the framework of the energy gap law could therefore reproduce our experimental observations. J-Aggregates can thus circumvent the common belief that lowering the excited state energies results in large nonradiative relaxation rates and are thus a pathway toward highly emissive organic dyes in the NIR regime.

Structural Stability of Insulin Oligomers and Protein Association-Dissociation Processes: Free Energy Landscape and Universal Role of Water

Association and dissociation of proteins are important biochemical events. In this Feature Article, we analyze the available studies of these processes for insulin oligomers in aqueous solution. We focus on the solvation of the insulin monomer in water, stability and dissociation of its dimer, and structural integrity of the hexamer. The intricate role of water in solvation of the dimer- and hexamer-forming surfaces, in long-range interactions between the monomers and the stability of the oligomers, is discussed. Ten water molecules inside the central cavity stabilize the structure of the insulin hexamer. We discuss how different order parameters can be used to understand the dissociation of the insulin dimer. The calculation of the rate using a recently computed multidimensional free energy provides considerable insight into the interplay between protein and water dynamics.

Real-time observation of photoionization-induced water migration dynamics in 4-methylformanilide-water by picosecond time-resolved infrared spectroscopy and ab initio molecular dynamics simulations

A novel time-resolved pump-probe spectroscopic approach that enables to keep high resolution in both the time and energy domain, nanosecond excitation-picosecond ionization-picosecond infrared probe (ns-ps-ps TRIR) spectroscopy, has been applied to the trans-4-methylformanilide-water (4MetFA-W) cluster. Water migration dynamics from the CO to the NH binding site in a peptide linkage triggered by photoionization of 4MetFA-W is directly monitored by the ps time evolution of IR spectra, and the presence of an intermediate state is revealed. The time evolution is analyzed by rate equations based on a four-state model of the migration dynamics. Time constants for the initial to the intermediate and hot product and to the final product are obtained. The acceleration of the dynamics by methyl substitution and the strong contribution of intracluster vibrational energy redistribution in the termination of the solvation dynamics is suggested. This picture is well confirmed by the ab initio on-the-fly molecular dynamics simulations. Vibrational assignments of 4MetFA and 4MetFA-W in the neutral (S₀ and S₁) and ionic (D₀) electronic states measured by ns IR dip and electron-impact IR photodissociation spectroscopy are also discussed prior to the results of time-resolved spectroscopy.

Hydrogen-Bond Dynamics and Energetics of Biological Water

Water molecules in the immediate vicinity of biomacromolecules and biomimetic organized assemblies often exhibit a markedly distinct behavior from that of their bulk counterparts. The overall sluggish behavior of biological water substantially affects the stability and integrity of biomolecules, as well as the successful execution of various crucial water-mediated biochemical phenomena. In this Minireview, insights are provided into the features of truncated hydrogen-bond networks that grant biological water its unique characteristics. In particular, experimental results and theoretical investigations, based on chemical kinetics, are presented that have shed light on the dynamics and energetics governing such characteristics. It is emphasized how such details help us to understand the energetics of biological water, an aspect relatively less explored than its dynamics. For instance, when biological water at hydrophilic or charged protein surfaces was explored, the free energy of H-bond breakage was found to be of the order of 0.4 kcal mol^-1 higher than that of bulk water.

Why are Vibrational Lines Narrow in Proteins?

The vibrational Stark effect in proteins yields line shifts indicative of strong internal electric fields up to a few volts per angstrom. These values are supported by numerical simulations of proteins. The simulations also show a significant breadth of field fluctuations translating to inhomogeneous broadening of vibrational lines. According to fluctuation-dissipation arguments, strong internal fields should lead to broad lines. Experimentally reported vibrational lines in proteins are, however, very narrow. This disconnect is explained here in terms of the insufficient (nonergodic) sampling of the protein's configurations on the lifetime of the vibrational probe. The slow component of the electric field fluctuations in proteins relaxes on the time scale of tens of nanoseconds and is dynamically frozen on the vibrational lifetime.

A simple method based on confocal microscopy and thick sections recognizes seven subphases in growth plate chondrocytes

This manuscript reports a novel procedure to imaging growth plate chondrocytes by using confocal microscopy. The method is based on fixed undecalcified bone samples, in-block staining with eosin, epoxy resin embedding and grinding to obtain thick sections. It is simple, inexpensive and provides three-dimensional images of entire chondrocytes inside their native lacunae. Quantitative analysis of volume, shape and cytoplasm density of chondrocytes at different strata of the growth plate allowed to objectively grade chondrocytes of the growth plate in seven different clusters. These seven categories of chondrocytes were subsequently evaluated by immunohistochemistry of some well-defined molecular landmarks of chondrocyte differentiation. Furthermore, immunohistochemical analysis of proteins responsible for ionic changes and water transport allowing chondrocyte swelling during hypertrophy was also performed. Results obtained indicate that four subphases can be defined in the pre-hypertrophic zone and three subphases in the hypertrophic zone, a fact that raises that chondrocytes of the growth plate are less homogeneous than usually considered when different zones are defined according to subjective cell morphological criteria. Results in the present study provide a technological innovation and gives new insights into the complexity of the process of chondrocyte differentiation in the growth plate.

Ultrafast Dynamics of Water-Protein Coupled Motions around the Surface of Eye Crystallin

Water dynamics on the protein surface mediate both protein structure and function. However, many questions remain about the role of the protein hydration layers in protein fluctuations and how the dynamics of these layers relate to specific protein properties. The fish eye lens protein γM7-crystallin (γM7) is found in vivo at extremely high concentrations nearing the packing limit, corresponding to only a few water layers between adjacent proteins. In this study, we conducted a site-specific probing of hydration water motions and side-chain dynamics at nine selected sites around the surface of γM7 using a tryptophan scan with femtosecond spectroscopy and NMR nuclear spin relaxation (NSR). We observed correlated fluctuations between hydration water and protein side chains on the time scales of a few picoseconds and hundreds of picoseconds, corresponding to local reorientations and network restructuring, respectively. These motions are heterogeneous over the protein surface and relate to the various steric and chemical properties of the local protein environment. Overall, we found that γM7 has relatively slower water dynamics within the hydration shell than a similar β-sheet protein, which may contribute to the high packing limit of this unique protein.

Structural Tuning of Quantum Decoherence and Coherent Energy Transfer in Photosynthetic Light Harvesting

Photosynthetic organisms capture energy from solar photons by constructing light-harvesting proteins containing arrays of electronic chromophores. Collective excitations (excitons) arise when energy transfer between chromophores is coherent, or wavelike, in character. Here we demonstrate experimentally that coherent energy transfer to the lowest-energy excitons is principally controlled in a light-harvesting protein by the temporal persistence of quantum coherence rather than by the strength of vibronic coupling. In the peridinin-chlorophyll protein from marine dinoflagellates, broad-band two-dimensional electronic spectroscopy reveals that replacing the native chlorophyll a acceptor chromophores with chlorophyll b slows energy transfer from the carotenoid peridinin to chlorophyll despite narrowing the donor-acceptor energy gap. The formyl substituent on the chlorophyll b macrocycle hastens decoherence by sensing the surrounding electrostatic noise. These findings demonstrate how quantum coherence enhances the efficiency of energy transfer despite being very short lived in light-harvesting proteins at physiological temperatures.

Mechanistic Regimes of Vibronic Transport in a Heterodimer and the Design Principle of Incoherent Vibronic Transport in Phycobiliproteins

Following the observation of coherent oscillations in nonlinear spectra of photosynthetic pigment protein complexes, in particular, phycobilliproteins such as PC645, coherent vibronic transport has been suggested as a design principle for novel light-harvesting materials. Vibronic transport between energetically remote pigments is coherent when the presence of a vibration resonant with the electronic energy gap supports transient delocalization between the electronic excited states. We establish the mechanism of vibronic transport for a model heterodimer across a wide range of molecular parameter values. The resulting mechanistic map demonstrates that the molecular parameters of phycobiliproteins in fact support incoherent vibronic transport. This result points to an important design principle: Incoherent vibronic transport is more efficient than a coherent mechanism when energetic disorder exceeds the coupling between the donor and vibrationally excited acceptor states. Finally, our results suggest that the role of coherent vibronic transport in pigment protein complexes should be reevaluated.

Local protein solvation drives direct down-conversion in phycobiliprotein PC645 via incoherent vibronic transport

The mechanisms controlling excitation energy transport (EET) in light-harvesting complexes remain controversial. Following the observation of long-lived beats in 2D electronic spectroscopy of PC645, vibronic coherence, the delocalization of excited states between pigments supported by a resonant vibration, has been proposed to enable direct excitation transport from the highest-energy to the lowest-energy pigments, bypassing a collection of intermediate states. Here, we instead show that for phycobiliprotein PC645 an incoherent vibronic transport mechanism is at play. We quantify the solvation dynamics of individual pigments using ab initio quantum mechanics/molecular mechanics (QM/MM) nuclear dynamics. Our atomistic spectral densities reproduce experimental observations ranging from absorption and fluorescence spectra to the timescales and selectivity of down-conversion observed in transient absorption measurements. We construct a general model for vibronic dimers and establish the parameter regimes of coherent and incoherent vibronic transport. We demonstrate that direct down-conversion in PC645 proceeds incoherently, enhanced by large reorganization energies and a broad collection of high-frequency vibrations. We suggest that a similar incoherent mechanism is appropriate across phycobiliproteins and represents a potential design principle for nanoscale control of EET.

Real-time observation of the photoionization-induced water rearrangement dynamics in the 5-hydroxyindole-water cluster by time-resolved IR spectroscopy

Solvation plays an essential role in controlling the mechanism and dynamics of chemical reactions in solution. The present study reveals that changes in the local solute-solvent interaction have a great impact on the timescale of solvent rearrangement dynamics. Time-resolved IR spectroscopy has been applied to a hydration rearrangement reaction in the monohydrated 5-hydroxyindole-water cluster induced by photoionization of the solute molecule. The water molecule changes the stable hydration site from the indolic NH site to the substituent OH site, both of which provide a strongly attractive potential for hydration. The rearrangement time constant amounts to 8 ± 2 ns, and is further slowed down by a factor of more than five at lower excess energy. These rearrangement times are slower by about three orders of magnitude than those reported for related systems where the water molecule is repelled from a repulsive part of the interaction potential toward an attractive well. The excess energy dependence of the time constant is well reproduced by RRKM theory. Differences in the reaction mechanism are discussed on the basis of energy relaxation dynamics.

Anomalous water dynamics at surfaces and interfaces: synergistic effects of confinement and surface interactions

In nature, water is often found in contact with surfaces that are extended on the scale of molecule size but small on a macroscopic scale. Examples include lipid bilayers and reverse micelles as well as biomolecules like proteins, DNA and zeolites, to name a few. While the presence of surfaces and interfaces interrupts the continuous hydrogen bond network of liquid water, confinement on a mesoscopic scale introduces new features. Even when extended on a molecular scale, natural and biological surfaces often have features (like charge, hydrophobicity) that vary on the scale of the molecular diameter of water. As a result, many new and exotic features, which are not seen in the bulk, appear in the dynamics of water close to the surface. These different behaviors bear the signature of both water-surface interactions and of confinement. In other words, the altered properties are the result of the synergistic effects of surface-water interactions and confinement. Ultrafast spectroscopy, theoretical modeling and computer simulations together form powerful synergistic approaches towards an understanding of the properties of confined water in such systems as nanocavities, reverse micelles (RMs), water inside and outside biomolecules like proteins and DNA, and also between two hydrophobic walls. We shall review the experimental results and place them in the context of theory and simulations. For water confined within RMs, we discuss the possible interference effects propagating from opposite surfaces. Similar interference is found to give rise to an effective attractive force between two hydrophobic surfaces immersed and kept fixed at a separation of d, with the force showing an exponential dependence on this distance. For protein and DNA hydration, we shall examine a multitude of timescales that arise from frustration effects due to the inherent heterogeneity of these surfaces. We pay particular attention to the role of orientational correlations and modification of the same due to interaction with the surfaces.

Complex Dynamics of Water in Protein Confinement

This paper studies single-molecule and collective dynamics of water confined in protein powders by means of molecular dynamics simulations. The single-particle dynamics show a modest retardation compared to the bulk but become highly stretched in the powder, with the stretching exponent of ≃0.2. The collective dynamics of the total water dipole are affected by intermolecular correlations inside water and by cross-correlations between the water and the protein. The dielectric spectrum of water in the powder has two nearly equal-amplitude peaks: a Debye peak with ≃16 ps relaxation time and a highly stretched peak with the relaxation time of ≃13 ns and a stretching exponent of ≃0.12. The slower relaxation component is not seen in the single-molecule correlation functions and can be assigned to elastic protein motions displacing water in the powder. The loss spectrum of the intermediate scattering function reported by neutron-scattering experiments is also highly stretched, with the high-frequency wing scaling according to a power law. Translational dynamics can become much slower in the powder than in the bulk but are overshadowed by the rotational loss in the overall loss spectrum of neutron scattering.

Ergodicity breaking of iron displacement in heme proteins

We present a model of the dynamical transition of atomic displacements in proteins. Increased mean-square displacement at higher temperatures is caused by the softening of the force constant for atomic/molecular displacements by electrostatic and van der Waals forces from the protein-water thermal bath. Displacement softening passes through a nonergodic dynamical transition when the relaxation time of the force-force correlation function enters, with increasing temperature, the instrumental observation window. Two crossover temperatures are identified. The lower crossover, presently connected to the glass transition, is related to the dynamical unfreezing of rotations of water molecules within nanodomains polarized by charged surface residues of the protein. The higher crossover temperature, usually assigned to the dynamical transition, marks the onset of water translations. All crossovers are ergodicity breaking transitions depending on the corresponding observation windows. Allowing stretched exponential relaxation of the protein-water thermal bath significantly improves the theory-experiment agreement when applied to solid protein samples studied by Mössbauer spectroscopy.

Water Dynamics in the Hydration Shells of Biomolecules

The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.

Theory and Electrochemistry of Cytochrome c

Extensive simulations of cytochrome c in solution are performed to address the apparent contradiction between large reorganization energies of protein electron transfer typically reported by atomistic simulations and much smaller values produced by protein electrochemistry. The two sets of data are reconciled by deriving the activation barrier for electrochemical reaction in terms of an effective reorganization energy composed of half the Stokes shift (characterizing the medium polarization in response to electron transfer) and the variance reorganization energy (characterizing the breadth of electrostatic fluctuations). This effective reorganization energy is much smaller than each of the two components contributing to it and is fully consistent with electrochemical measurements. Calculations in the range of temperatures between 280 and 360 K combine long, classical molecular dynamics simulations with quantum calculations of the protein active site. The results agree with the Arrhenius plots for the reaction rates and with cyclic voltammetry of cytochrome c immobilized on self-assembled monolayers. Small effective reorganization energy, and the resulting small activation barrier, is a general phenomenology of protein electron transfer allowing fast electron transport within biological energy chains.

Conformational Homogeneity in the Pr Isomer of Phytochrome Cph1

Numerous time-resolved studies of the P_r to P_fr photoisomerization in phytochrome Cph1 have revealed multiphasic excited-state decay kinetics. It remains unclear whether these kinetics arise from multiple ground-state conformational subpopulations or from a single ground-state conformation that undergoes an excited-state photoisomerization process-either branching on the excited state or relaxing through multiple sequential intermediates. Many studies have attempted to resolve this debate by fitting the measured dynamics to proposed kinetic models, arriving at different conclusions. Here we probe spectral signatures of ground-state heterogeneity of P_r. Two-dimensional electronic spectra display negligible inhomogeneous line broadening, and vibrational coherence spectra extracted from transient absorption measurements do not contain nodes and phase shifts at the fluorescence maximum. These spectroscopic results support the homogeneous model, in which the primary photochemical transformation of P_r to Lumi-R occurs adiabatically on the excited-state potential energy surface.

On the coupling between the dynamics of protein and water

The solvation entropy of flexible protein regions is higher than that of rigid regions and contributes differently to the overall thermodynamic stability.

Broad-Band Pump-Probe Spectroscopy Quantifies Ultrafast Solvation Dynamics of Proteins and Molecules

In this work, we demonstrate the use of broad-band pump-probe spectroscopy to measure femtosecond solvation dynamics. We report studies of a rhodamine dye in methanol and cryptophyte algae light-harvesting proteins in aqueous suspension. Broad-band impulsive excitation generates a vibrational wavepacket that oscillates on the excited-state potential energy surface, destructively interfering with itself at the minimum of the surface. This destructive interference gives rise to a node at a certain probe wavelength that varies with time. This reveals the Gibbs free-energy changes of the excited-state potential energy surface, which equates to the solvation time correlation function. This method captures the inertial solvent response of water (∼40 fs) and the bimodal inertial response of methanol (∼40 and ∼150 fs) and reveals how protein-buried chromophores are sensitive to the solvent dynamics inside and outside of the protein environment.

Bacteriophytochrome Photoisomerization Proceeds Homogeneously Despite Heterogeneity in Ground State

Phytochromes are red/far-red photoreceptors that are widely distributed in plants and prokaryotes. Ultrafast photoisomerization of a double bond in a biliverdin cofactor or other linear tetrapyrrole drives their photoactivity, but their photodynamics are only partially understood. Multiexponential dynamics were observed in previous ultrafast spectroscopic studies and were attributed to heterogeneous populations of the pigment-protein complex. In this work, two-dimensional photon echo spectroscopy was applied to study dynamics of the bacteriophytochromes RpBphP2 and PaBphP. Two-dimensional photon echo spectroscopy can simultaneously resolve inhomogeneity in ensembles and fast dynamics by correlating pump wavelength with the emitted signal wavelength. The distribution of absorption and emission energies within the same state indicates an ensemble of heterogeneous protein environments that are spectroscopically distinct. However, the lifetimes of the dynamics are uniform across the ensemble, suggesting a homogeneous model involving sequential intermediates for the initial photodynamics of isomerization.

Dynamics and mechanism of ultrafast water-protein interactions

Protein hydration is essential to its structure, dynamics, and function, but water-protein interactions have not been directly observed in real time at physiological temperature to our awareness. By using a tryptophan scan with femtosecond spectroscopy, we simultaneously measured the hydration water dynamics and protein side-chain motions with temperature dependence. We observed the heterogeneous hydration dynamics around the global protein surface with two types of coupled motions, collective water/side-chain reorientation in a few picoseconds and cooperative water/side-chain restructuring in tens of picoseconds. The ultrafast dynamics in hundreds of femtoseconds is from the outer-layer, bulk-type mobile water molecules in the hydration shell. We also found that the hydration water dynamics are always faster than protein side-chain relaxations but with the same energy barriers, indicating hydration shell fluctuations driving protein side-chain motions on the picosecond time scales and thus elucidating their ultimate relationship.

Towards an ab initio description of the optical spectra of light-harvesting antennae: application to the CP29 complex of photosystem II

Light-harvesting pigment-protein complexes (PPC) represent the fundamental units through which the photosynthetic organisms absorb sunlight and funnel the energy to the reaction centre for carrying out the primary energy conversion reactions of photosynthesis. Here we apply a multiscale computational strategy to a specific PPC present in the photosystem II of plants and algae (CP29) to investigate in what detail should the environment effects due to protein and membrane/solvent be included for an accurate description of optical spectra. We find that a refinement of the crystal structure is needed before any meaningful quantum chemical calculations of pigment transition energies can be performed. For this purpose we apply classical molecular dynamics simulations of the PPC within its natural environment and we perform ab initio computations of the exciton Hamiltonian of the complex, including the environment either implicitly by the polarizable continuum model (PCM) or explicitly using the polarizable QM/MM methodology (MMPol). However, PCM essentially leads to an unspecific redshift of all transition energies, and MMPol is able to reveal site-specific changes in the optical properties of the pigments. Based on the latter and the excitonic couplings obtained within a polarizable QM/MM methodology, optical spectra are calculated, which are in good qualitative agreement with experimental data. A weakness of the approach is however found in the overestimation of the fluctuations of the excitonic parameters of the pigments along the MD trajectory. An explanation for such a finding in terms of the limits of the force fields commonly used for protein cofactors is presented and discussed.

Nanoprecipitation and Spectroscopic Characterization of Curcumin-Encapsulated Polyester Nanoparticles

Curcumin-encapsulated polyester nanoparticles (Cur-polyester NPs) of approximately 100 nm diameter with a negatively charged surface were prepared using a one-step nanoprecipitation method. The Cur-polyester NPs were prepared using polylactic acid, poly(D,L-lactic-co-glycolic acid) and poly(ϵ-caprolactone) without any emulsifier or surfactant. The encapsulation of curcumin in these polyester NPs greatly suppresses curcumin degradation in the aqueous environment due to its segregation from water. In addition, the fluorescence of curcumin in polyester NPs has a quantum yield of 4 to 5%, which is higher than that of curcumin in micellar systems and comparable to those in organic solvents, further supporting the idea that the polyester NPs are capable of excluding water from curcumin. Furthermore, the results from femtosecond fluorescence upconversion spectroscopy reveal that there is a decrease in the signal amplitude corresponding to solvent reorganization of excited state curcumin in the polyester NPs compared with curcumin in micellar systems. The Cur-polyester NPs also show a lack of deuterium isotope effect in the fluorescence lifetime. These results indicate that the interaction between curcumin and water in the polyester NPs is significantly weaker than that in micelles. Therefore, the aqueous stability of curcumin is greatly improved due to highly effective segregation from water. The overall outcome suggests that the polyester NPs prepared using the method reported herein are an attractive system for encapsulating and stabilizing curcumin in the aqueous environment.

Solvation Dynamics in Different Phases of the Lyotropic Liquid Crystalline System

Reverse hexagonal (HII) liquid crystalline material based on glycerol monooleate (GMO) is considered as a potential carrier for drugs and other important biomolecules due to its thermotropic phase change and excellent morphology. In this work, the dynamics of encapsulated water, which plays important role in stabilization and formation of reverse hexagonal mesophase, has been investigated by time dependent Stokes shift method using Coumarin-343 as a solvation probe. The formation of the reverse hexagonal mesophase (HII) and transformation to the L2 phase have been monitored using small-angle X-ray scattering and polarized light microscopy experiments. REES studies suggest the existence of different polar regions in both HII and L2 systems. The solvation dynamics study inside the reverse hexagonal (HII) phase reveals the existence of two different types of water molecules exhibiting dynamics on a 120-900 ps time scale. The estimated diffusion coefficients of both types of water molecules obtained from the observed dynamics are in good agreement with the measured diffusion coefficient collected from the NMR study. The calculated activation energy is found to be 2.05 kcal/mol, which is associated with coupled rotational-translational water relaxation dynamics upon the transition from "bound" to "quasi-free" state. The observed ∼2 ns faster dynamics of the L2 phase compared to the HII phase may be associated with both the phase transformation as well as thermotropic effect on the relaxation process. Microviscosities calculated from time-resolved anisotropy studies infer that the interface is almost ∼22 times higher viscous than the central part of the cylinder. Overall, our results reveal the unique dynamical features of water inside the cylinder of reverse hexagonal and inverse micellar phases.

Dipolar Nanodomains in Protein Hydration Shells

The network of hydrogen bonds characteristic of bulk water is significantly disturbed at the protein-water interface, where local fields induce mutually frustrated dipolar domains with potentially novel structure and dynamics. Here the dipolar susceptibility of hydration shells of lysozyme is studied by molecular dynamics simulations in a broad range of temperatures, 140-300 K. The real part of the susceptibility passes through a broad maximum as a function of temperature. The maximum shifts to higher temperatures with increasing frequency of the dielectric experiment. This phenomenology is consistent with that reported for bulk relaxor ferroelectrics, where it is related to the formation of dipolar nanodomains. Nanodomains in the hydration shell extend 12-15 Å from the protein surface into the bulk. Their dynamics are significantly slower than the dynamics of bulk water. The domains dynamically freeze into a ferroelectric glass below 160 K, at which point the Arrhenius plot of the dipolar relaxation time becomes significantly steeper.

The Fenna-Matthews-Olson protein revisited: a fully polarizable (TD)DFT/MM description

We report a combined molecular dynamics and quantum mechanics (QM)/molecular mechanics (MM) analysis of the excitonic properties of the Fenna-Matthews-Olson (FMO) protein by using a polarizable MM model combined with a time-dependent density functional theory description. Overall, our results indicate that structural fluctuations, electrostatic interactions, and short-range quantum effects can significantly modulate the model Hamiltonian parameters (site energies and couplings). We find that the specific interactions with the axial ligand and the hydrogen-bonded residues are responsible for the energy ladder, with their effects being mainly due to electrostatic interactions, but with short-range quantum contributions that are not negligible. In addition, a striking modulation of the screening effects experienced by the BChl pairs, due to the heterogeneous polarizability of the FMO and solvent environment, is observed. Finally, we find that the exciton model gives a reliable description of the delocalized excited states in the complex.

Implementing electrostatic polarization cannot fill the gap between experimental and theoretical measurements for the ultrafast fluorescence decay of myoglobin

Over the past few years, time-dependent ultrafast fluorescence spectroscopy method has been applied to the study of protein dynamics. However, observations from these experiments are in a controversy with other experimental studies. Participating of theoretical methods in this debate has not reconciled the contradiction, because the predicted initial relaxation from computer simulations is one-order faster than the ultrafast fluorescence spectroscopy experiment. In those simulations, pairwise force fields are employed, which have been shown to underestimate the roughness of the free energy landscape. Therefore, the relaxation rate of protein and water molecules under pairwise force fields is falsely exaggerated. In this work, we compared the relaxations of tryptophan/environment interaction under linear response approximation employing pairwise, polarized, and polarizable force fields. Results show that although the relaxation can be slowed down to a certain extent, the large gap between experiment and theory still cannot be filled.

Quantum Calculations in Solution for Large to Very Large Molecules: A New Linear Scaling QM/Continuum Approach

We present a new implementation of continuum solvation models for semiempirical Hamiltonians that allows the description of environmental effects on very large molecular systems. In this approach based on a domain decomposition strategy of the COSMO model (ddCOSMO), the solution to the COSMO equations is no longer the computational bottleneck but becomes a negligible part of the overall computation time. In this Letter, we analyze the computational impact of COSMO on the solution of the SCF equations for large to very large molecules, using semiempirical Hamiltonians, for both the new ddCOSMO implementation and the most recent, linear scaling one, based on the fast multipole method. A further analysis is on the simulation of the UV/visible spectrum of a light-harvesting pigment-protein complex. All of the results show how the new ddCOSMO algorithm paves the way to routine computations for large molecular systems in the condensed phase.

Charge invariant protein-water relaxation in GB1 via ultrafast tryptophan fluorescence

The protein–water interface is a critical determinant of protein structure and function, yet the precise nature of dynamics in this complex system remains elusive. Tryptophan fluorescence has become the probe of choice for such dynamics on the picosecond time scale (especially via fluorescence “upconversion”). In the absence of ultrafast (“quasi-static”) quenching from nearby groups, the TDFSS (time-dependent fluorescence Stokes shift) for exposed Trp directly reports on dipolar relaxation near the interface (both water and polypeptide). The small protein GB1 contains a single Trp (W43) of this type, and its structure is refractory to pH above 3. Thus, it can be used to examine the dependence of dipolar relaxation upon charge reconfiguration with titration. Somewhat surprisingly, the dipolar dynamics in the 100 fs to 100 ps range were unchanged with pH, although nanosecond yield, rates, and access all changed. These results were rationalized with the help of molecular dynamics (including QM-MM) simulations that reveal a balancing, sometimes even countervailing influence of protein and water dipoles. Interestingly, these simulations also showed the dominant influence of water molecules which are associated with the protein interface for up to 30 ps yet free to rotate at approximately “bulk” water rates.