Ultrafast hydration dynamics in melittin folding and aggregation: helix formation and tetramer self-assembly

Melittin, an amphipathic peptide from honeybee venom, consists of 26 amino acid residues and adopts different conformations from a random coil, to an alpha-helix, and to a self-assembled tetramer under certain aqueous environments. We report here our systematic studies of the hydration dynamics in these conformations using single intrinsic tryptophan (W19) as a molecular probe. With femtosecond resolution, we observed the solvation dynamics occurring in 0.62 and 14.7 ps in a random-coiled primary structure. The former represents bulklike water motion, and the latter reflects surface-type hydration dynamics of proteins. As a comparison, a model tripeptide (KWK) was also studied. At a membrane-water interface, melittin folds into a secondary alpha-helical structure, and the interfacial water motion was found to take as long as 114 ps, indicating a well-ordered water structure along the membrane surface. In high-salt aqueous solution, the dielectric screening and ionic solvation promote the hydrophobic core collapse in melittin aggregation and facilitate the tetramer formation. This self-assembled tertiary structure is also stabilized by the strong hydrophilic interactions of charged C-terminal residues and associated ions with water molecules in the two assembled regions. The hydration dynamics was observed to occur in 87 ps, significantly slower than typical water relaxation at protein surfaces but similar to water motion at membrane interfaces. Thus, the observed time scale of approximately 100 ps probably implies appropriate water mobility for mediating the formation of high-order structures of melittin in an alpha-helix and a self-assembled tetramer. These results elucidate the critical role of hydration dynamics in peptide conformational transitions and protein structural stability and integrity.     

Relaxation dynamics of 2,7- and 3,6-distyrylcarbazoles in solutions and in solid films: mechanism for efficient nonradiative deactivation in the 3,6-linked carbazole

We performed time-resolved spectral investigations of two distyrylcarbazole derivatives, 2,7- and 3,6-distyrylcarbazole (2,7-DPVTCz and 3,6-DPVTCz, respectively), in dilute toluene solution and in solid films mixed with poly(methyl methacrylate) (PMMA). The lifetime of 2,7-DPVTCz in its excited state in solution is approximately 100 times as great as that of 3,6-DPVTCz, consistent with their photophysical nature. The former shows intense emission, but the latter is nearly nonfluorescent in a free environment. Moreover, the lifetime of 3,6-DPVTCz in its excited state increased also approximately 100 times when the molecule was encapsulated in a 3,6-DPVTCz/PMMA solid film, indicating that intramolecular motion of the molecule significantly affects the observed relaxation dynamics in a confined environment. Calculations on the excited states indicate that an efficient intersystem crossing is activated upon twisting of the bridged C-C single bond in a free 3,6-linked carbazole; such efficient deactivation is impractical in 2,7-linked carbazole or for 3,6-linked carbazole in a PMMA matrix. Information obtained from experiments on femtosecond fluorescence enables us to distinguish crucial relaxation processes in the excited state for a profound understanding of the details of vibrational and electronic relaxations of 3,6-DPVTCz in solution.     

Nonadiabatic molecular dynamics simulations of correlated electrons in solution. 1. Full configuration interaction (CI) excited-state relaxation dynamics of hydrated dielectrons

The hydrated dielectron is composed of two excess electrons dissolved in liquid water that occupy a single cavity; in both its singlet and triplet spin states there is a significant exchange interaction so the two electrons cannot be considered to be independent. In this paper and the following paper,we present the results of mixed quantum/classical molecular dynamics simulations of the nonadiabatic relaxation dynamics of photoexcited hydrated dielectrons, where we use full configuration interaction (CI) to solve for the two-electron wave function at every simulation time step. To the best of our knowledge, this represents the first systematic treatment of excited-state solvation dynamics where the multiple-electron problem is solved exactly. The simulations show that the effects of exchange and correlation contribute significantly to the relaxation dynamics. For example, spin-singlet dielectrons relax to the ground state on a time scale similar to that of single electrons excited at the same energy, but spin-triplet dielectrons relax much faster. The difference in relaxation dynamics is caused by exchange and correlation: The Pauli exclusion principle imposes very different electronic structure when the electrons' spins are singlet paired than when they are triplet paired, altering the available nonadiabatic relaxation pathways. In addition, we monitor how electronic correlation changes dynamically during nonadiabatic relaxation and show that solvent dynamics cause electron correlation to evolve quite differently for singlet and triplet dielectrons. Despite such differences, our calculations show that both spin states are stable to excited-state dissociation, but that the excited-state stability has different origins for the two spin states. For singlet dielectrons, the stability depends on whether the solvent structure can rearrange to create a second cavity before the ground state is reached. For triplet dielectrons, in contrast, electronic correlation ensures that the two electrons do not dissociate, even if the dielectron is artificially kept from reaching the ground state. In addition, both singlet and triplet dielectrons change shape dramatically during relaxation, so that linear response fails to describe the solvation dynamics for either spin state. In the following paper (Larsen, R. E.; Schwartz, B. J. J. Phys. Chem. B 2006, 110, 9692), we use these simulations to calculate the pump-probe spectroscopic signal expected for photoexcited hydrated dielectrons and to predict an experiment to observe hydrated dielectrons directly.     

Diffusive dynamics on multidimensional rough free energy surfaces

The dynamics of processes relevant to chemistry and biophysics on rough free energy landscapes is investigated using a recently developed algorithm to solve the Smoluchowski equation. Two different processes are considered: ligand rebinding in MbCO and protein folding. For the rebinding dynamics of carbon monoxide (CO) to native myoglobin (Mb) from locations around the active site, the two-dimensional free energy surface (FES) is constructed using extensive molecular dynamics simulations. The surface describes the minima in the A state (bound MbCO), CO in the distal pocket and in the Xe4 pocket, and the transitions between these states and allows to study the diffusion of CO in detail. For the folding dynamics of protein G, a previously determined two-dimensional FES was available. To follow the diffusive dynamics on these rough free energy surfaces, the Smoluchowski equation is solved using the recently developed hierarchical discrete approximation method. From the relaxation of the initial nonequilibrium distribution, experimentally accessible quantities such as the rebinding time for CO or the folding time for protein G can be calculated. It is found that the free energy barrier for CO in the Xe4 pocket and in the distal pocket (B state) closer to the heme iron is approximately 6 kcal/mol which is considerably larger than the inner barrier which separates the bound state and the B state. For the folding of protein G, a barrier of approximately 10 kcal/mol between the unfolded and the folded state is consistent with folding times of the order of milliseconds.     

Laser pulse trains for controlling excited state dynamics of adenine in water

We investigate theoretically the control of the ultrafast excited state dynamics of adenine in water by laser pulse trains, with the aim to extend the excited state lifetime and to suppress nonradiative relaxation processes. For this purpose, we introduce the combination of our field-induced surface hopping method (FISH) with the quantum mechanical-molecular mechanical (QM/MM) technique for simulating the laser-driven dynamics in the condensed phase under explicit inclusion of the solvent environment. Moreover, we employ parametric pulse shaping in the frequency domain in order to design simplified laser pulse trains allowing to establish a direct link between the pulse parameters and the controlled dynamics. We construct pulse trains which achieve a high excitation efficiency and at the same time keep a high excited state population for a significantly extended time period compared to the uncontrolled dynamics. The control mechanism involves a sequential cycling of the population between the lowest and higher excited states, thereby utilizing the properties of the corresponding potential energy surfaces to avoid conical intersections and thus to suppress the nonradiative decay to the ground state. Our findings provide a means to increase the fluorescence yield of molecules with an intrinsically very short excited state lifetime, which can lead to novel applications of shaped laser fields in the context of biosensing.     

Energetics and dynamics in MbCN: CN--vibrational relaxation from molecular dynamics simulations

The dynamics of the cyanide anion bound to sperm-whale myoglobin is investigated using atomistic simulations. With density-functional theory, a 2D potential energy surface for the cyanide-heme complex is calculated. Two deep minima with a stabilization energy of approximately 50 kcal/mol corresponding to two different binding orientations (Fe-CN and Fe-NC) of the ligand are found. The Fe-CN conformation is favored over Fe-NC by several kcal/mol. Mixed quantum mechanics/molecular mechanics calculations show that the binding orientation affects the bond strength of the ligand, with a significantly different bond length and a 25 cm-1 shift in the fundamental CN-frequency. For the molecular dynamics (MD) simulations, a 3-center fluctuating charge model for the Fe-CN unit is developed that captures polarization and ligand-metal charge transfer. Stability arguments based on the energetics around the active site and the CN- frequency shifts suggest that the Fe-CN conformation with epsilon-protonation of His epsilon 64 are most likely, which is in agreement with experiment. Both equilibrium and nonequilibrium MD simulations are carried out to investigate the relaxation time scale and possible relaxation pathways in bound MbCN. The nonequilibrium MD simulations with a vibrationally excited ligand reveal that vibrational relaxation takes place on a time scale of hundreds of picoseconds within the active site. This finding supports the hypothesis that the experimentally observed relaxation rate (3.6 ps) reflects the repopulation of the electronic ground state.     

Molecular modeling and simulation of conjugated polymer oligomers: ground and excited state chain dynamics of PPV in the gas phase

The ground and excited state dynamics of poly(p-phenylenevinylene) (PPV) chains is studied through an implementation of mixed quantum/classical molecular dynamics simulation. The model used in the simulations combines the semiempirical Pariser-Parr-Pople (PPP) Hamiltonian to treat the pi molecular electronic structure with a mechanical force field capturing all other aspects. Nuclear degrees of freedom are treated classically. We first validate the model by simulating PPV chains of various length, and evaluate the absorption spectra. The thermal disorder contribution to the breadth of the first absorption band is estimated to be 0.2 eV at T = 300 K. To investigate the relationship between the emission and chain conformation, we simulate an isolated ten unit chain of PPV in the ground and the lowest excited state. The emission spectrum, red-shifted with respect to absorption of about 0.2 eV as found in experiments, shows a structured line shape that we relate to the photoinduced CC bond distortions. In accord with earlier studies, the exciton self-traps in the middle of the chain. We introduce two collective variables that reflect geometrical distortion, and find these to be effective in describing the contribution of chain conformation to the emission spectrum. The collective variables are also shown to be effective in describing the bond relaxation dynamics after photoexcitation. Such a relaxation is found to occur in approximately 100 fs and is guided by a compensatory release of energy between the double and single bonds in the vinylene junctions and p-phenyl rings. Finally, we find that the chain has a very slight preference for a more planar conformation in the excited state, compared to the ground state. However, the thermal motions induce the chain to explore out-of-plane conformations in both the ground and the excited states with an amplitude significantly greater than this difference.     

Excimer formation dynamics in the isolated tetracene dimer

The understanding of excimer formation and its interplay with the singlet-correlated triplet pair state ¹(TT) is of high significance for the development of efficient organic electronics. Here, we study the photoinduced dynamics of the tetracene dimer in the gas phase by time-resolved photoionisation and photoion imaging experiments as well as nonadiabatic dynamics simulations in order to obtain mechanistic insight into the excimer formation dynamics. The experiments are performed using a picosecond laser system for excitation into the S₂ state and reveal a biexponential time dependence. The time constants, obtained as a function of excess energy, lie in the range between ≈10 ps and 100 ps and are assigned to the relaxation of the excimer on the S₁ surface and to its deactivation to the ground state. Simulations of the quantum-classical photodynamics are carried out in the frame of the semi-empirical CISD and TD-lc-DFTB methods. Both theoretical approaches reveal a dominating relaxation pathway that is characterised by the formation of a perfectly stacked excimer. TD-lc-DFTB simulations have also uncovered a second relaxation channel into a less stable dimer conformation in the S₁ state. Both methods have consistently shown that the electronic and geometric relaxation to the excimer state is completed in less than 10 ps. The inclusion of doubly excited states in the CISD dynamics and their diabatisation further allowed to observe a transient population of the ¹(TT) state, which, however, gets depopulated on a timescale of 8 ps, leading finally to the trapping in the excimer minimum.

The understanding of excimer formation and its interplay with the singlet-correlated triplet pair state ¹(TT) is of high significance for the development of efficient organic electronics.     

Charge-transfer-to-solvent-driven dissolution dynamics of I- (H2O)2-5 upon excitation: excited-state ab initio molecular dynamics simulations

In contrast to the extensive theoretical investigation of the solvation phenomena, the dissolution phenomena have hardly been investigated theoretically. Upon the excitation of hydrated halides, which are important substances in atmospheric chemistry, an excess electron transfers from the anionic precursor (halide anion) to the solvent and is stabilized by the water cluster. This results in the dissociation of hydrated halides into halide radicals and electron-water clusters. Here we demonstrate the charge-transfer-to-solvent (CTTS)-driven femtosecond-scale dissolution dynamics for I-(H2O)n=2-5 clusters using excited state (ES) ab initio molecular dynamics (AIMD) simulations employing the complete-active-space self-consistent-field (CASSCF) method. This study shows that after the iodine radical is released from I-(H2O)n=2-5, a simple population decay is observed for small clusters (2 相似文献   

Energy relaxation dynamics of the hydration complex of hydroxide

We use polarization-resolved mid-infrared pump-probe spectroscopy to study the dynamics of the hydration shells of hydroxide ions (OH(-)). We excite the OH stretch vibrations of H(2)O molecules solvating the OH(-) ion and observe that this excitation decays with a relaxation time constant T(1) of 200 fs. This relaxation is followed by a thermalization process that becomes slower with increasing concentration of OH(-). The prethermalized state is observed to be anisotropic, showing that the energy of the excited OH stretch vibrations is dissipated within the hydration complex. The anisotropy of the prethermalized state decays both as a result of the reorientation of the OH(-) hydration complex and heat diffusion from the excited complexes to unexcited complexes. Modeling the anisotropy data at different concentrations allows for an accurate estimate of the number of water molecules in the hydration shell of OH(-), the reorientation dynamics of the OH(-) hydration complex, and the molecular-scale heat diffusivity.     

Nanosecond folding dynamics of a three-stranded beta-sheet

We studied conformational stability and folding kinetics of a three-stranded beta-sheet containing two rigid turns. Static infrared measurements indicate that this beta-sheet undergoes a broad but cooperative thermal unfolding transition with a midpoint at approximately 53 degrees C. Interestingly, time-resolved infrared experiments show that its relaxation kinetics in response to a temperature-jump (T-jump) occur on the nanosecond time scale (e.g., the relaxation time is approximately 140 ns at 35.0 degrees C), thereby suggesting that the conformational relaxation encounters only a small free energy barrier or even proceeds in a downhill manner. Further Langevin dynamics simulations suggest that the observed T-jump relaxation kinetics could be modeled by a conformational diffusion process along a single-well free energy profile, which allowed us to determine the effective diffusion constant and also the roughness of the folding energy landscape.     

Single and multi-exciton dynamics in aqueous protochlorophyllide aggregates

In plants, the oxidoreductase enzyme POR reduces protochlorophyllide (Pchlide) into chlorophyllide (Chlide), using NADPH as a cofactor. The reduction involves the transfer of two electrons and two protons to the C17═C18 double bond of Pchlide, and the reaction is initiated by the absorption of light by Pchlide itself. In this work we have studied the excited state dynamics of Pchlide dissolved in water, where it forms excitonically coupled aggregates, by ultrafast time-resolved transient absorption and fluorescence experiments performed in the 480-720 nm visible region and in the 1780-1590 cm(-1) mid-IR region. The ground state visible absorption spectrum of aqueous Pchlide red shifts and broadens in comparison to the spectrum of monomeric Pchlide in organic solvents. The population of the one-exciton state occurs at low excitation densities, of <1 photon per aggregate. We characterized the multiexciton manifolds spectra by measuring the absorption difference spectra at increasingly higher photon densities. The multiexciton states are characterized by blue-shifted stimulated emission and red-shifted excited state absorption in comparison to those of the one-exciton manifold. The relaxation dynamics of the multiexciton manifolds into the one-exciton manifold is found to occur in ～10 ps. This surprisingly slow rate we suggest is due to the intrinsic charge transfer character of the PChlide excited state that leads to solvation, stabilizing the CT state, and subsequent charge recombination, which limits the exciton relaxation.     

Elucidating the initial dynamics of electron photodetachment from atoms in liquids using variably-time-delayed resonant multiphoton ionization

We study the photodetachment of electrons from sodium anions in room temperature liquid tetrahydrofuran (THF) using a new type of three-pulse pump-probe spectroscopy. Our experiments use two variably-time-delayed pulses for excitation in what is essentially a resonant 1+1 two-photon ionization: By varying the arrival time of the second excitation pulse, we can directly observe how solvent motions stabilize and trap the excited electron prior to electron detachment. Moreover, by varying the arrival times of the ionization (excitation) and probe pulses, we also can determine the fate of the photoionized electrons and the distance they are ejected from their parent Na atoms. We find that as solvent reorganization proceeds, the second excitation pulse becomes less effective at achieving photoionization, and that the solvent motions that stabilize the excited electron following the first excitation pulse occur over a time of approximately 450 fs. We also find that there is no spectroscopic evidence for significant solvent relaxation after detachment of the electron is complete. In combination with the results of previous experiments and molecular dynamics simulations, the data provide new insight into the role of the solvent in solution-phase electron detachment and charge-transfer-to-solvent reactions.     

Time-resolved relaxation dynamics of Hgn- (11 <or = n <or = 16,n = 18) clusters following intraband excitation at 1.5 eV

Electron-nuclear relaxation dynamics are studied in Hg(n) (-) (11 相似文献   

Hyper-aging dynamics of nanoclay suspension

Aqueous suspension of nanoclay Laponite undergoes structural evolution as a function of time, which enhances its elasticity and relaxation time. In this work, we employ an effective time approach to investigate long-term relaxation dynamics by carrying out creep experiments. Typically, we observe that the monotonic evolution of elastic modulus shifts to lower aging times, while maxima in viscous moduli get progressively broader for experiments carried out on a later date after preparation (idle time) of the nanoclay suspension. Application of effective time theory produces a superposition of all the creep curves irrespective of their initial state. The resulting dependence of the relaxation time on aging time shows very strong hyper-aging dynamics at short idle times, which progressively weakens to demonstrate a linear dependence in the limit of very long idle times. Remarkably, this behavior of nanoclay suspensions is akin to that observed for polymeric glasses. Consideration of aging as a first-order process suggests that continued hyper-aging dynamics causes cessation of aging. The dependence of relaxation time on aging time, therefore, must attenuate eventually producing linear or weaker dependence on time in order to approach a progressively low-energy state in the limit of very long times as observed experimentally. We also develop a simple scaling model based on a concept of aging of an energy well, which qualitatively captures various experimental observations very well, leading to profound insight into the hyper-aging dynamics of nanoclay suspensions.     

Molecular dynamics simulations of electron-alkali cation pairs in bulk water

The structural, dynamic, and thermodynamic properties of an excess electron interacting with an alkali cation (Na+, K+, Li+) in bulk water were investigated by means of a mixed quantum-classical molecular dynamics simulation technique. This study includes a reparametrization of the electron-cation pseudopotentials. The free energy calculations for all three systems show that a contact electron-cation pair can be observed, which is either as stable as the dissociated pair (Li+) or more stable by only a few kT (Na+, K+). Given that the dissociation barrier is also quite small, we suggest that the average cation-electron distance in the experiments at room temperature will not depend on this free energy profile but rather on the minimization of the Coulombic repulsive interaction between like charges in the solvent medium. This enables us to compare the present molecular dynamics simulations with the spectroscopic data obtained for different ionic strengths. The overall trend of the UV-vis hydrated absorption spectra, namely, the shift toward shorter wavelengths at high ionic strengths, is fairly well reproduced. This confirms our hypothesis of statistical distribution of the cations and solvated electrons.     

Ab initio nonadiabatic dynamics study of ultrafast radiationless decay over conical intersections illustrated on the Na3F cluster

We present a theoretical approach for the ultrafast nonadiabatic dynamics based on the ab initio molecular dynamics carried out "on the fly" in the framework of the configuration interaction method combined with Tully's surface hopping algorithm for nonadiabatic transitions. This approach combined with our Wigner distribution approach allows us to perform accurate simulations of femtosecond pump-probe spectra in the systems where radiationless transitions among electronic states take place. In this paper we illustrate this by theoretical simulation of ultrafast processes and nonradiative relaxation in the Na(3)F cluster, involving three excited states and the ground electronic state. Furthermore, we show that our accurate simulation of the photoionization pump-probe spectrum is in full agreement with the experimental signal. Based on the nonadiabatic dynamics at high level of accuracy and taking into account all degrees of freedom, the nonradiative lifetime for the 1 (1)B(1) excited state of Na(3)F has been determined to be approximately 900 fs.     

Excited‐State Deactivation of Branched Phthalocyanine Compounds

The excited‐state relaxation dynamics and chromophore interactions in two phthalocyanine compounds (bis‐ and trisphthalocyanines) are studied by using steady‐state and femtosecond transient absorption spectral measurements, where the excited‐state energy‐transfer mechanism is explored. By exciting phthalocyanine compounds to their second electronically excited states and probing the subsequent relaxation dynamics, a multitude of deactivation pathways are identified. The transient absorption spectra show the relaxation pathway from the exciton state to excimer state and then back to the ground state in bisphthalocyanine (bis‐Pc). In trisphthalocyanine (tris‐Pc), the monomeric and dimeric subunits are excited and the excitation energy transfers from the monomeric vibrationally hot S1 state to the exciton state of a pre‐associated dimer, with subsequent relaxation to the ground state through the excimer state. The theoretical calculations and steady‐state spectra also show a face‐to‐face conformation in bis‐Pc, whereas in tris‐Pc, two of the three phthalocyanine branches form a pre‐associated face‐to‐face dimeric conformation with the third one acting as a monomeric unit; this is consistent with the results of the transient absorption experiments from the perspective of molecular structure. The detailed structure–property relationships in phthalocyanine compounds is useful for exploring the function of molecular aggregates in energy migration of natural photosynthesis systems.