Using molecular dynamics and quantum mechanics calculations to model fluorescence observables.

We provide a critical examination of two different methods for generating a donor-acceptor electronic coupling trajectory from a molecular dynamics (MD) trajectory and three methods for sampling that coupling trajectory, allowing the modeling of experimental observables directly from the MD simulation. In the first coupling method we perform a single quantum-mechanical (QM) calculation to characterize the excited state behavior, specifically the transition dipole moment, of the fluorescent probe, which is then mapped onto the configuration space sampled by MD. We then utilize these transition dipoles within the ideal dipole approximation (IDA) to determine the electronic coupling between the probes that mediates the transfer of energy. In the second method we perform a QM calculation on each snapshot and use the complete transition densities to calculate the electronic coupling without need for the IDA. The resulting coupling trajectories are then sampled using three methods ranging from an independent sampling of each trajectory point (the independent snapshot method) to a Markov chain treatment that accounts for the dynamics of the coupling in determining effective rates. The results show that the IDA significantly overestimates the energy transfer rate (by a factor of 2.6) during the portions of the trajectory in which the probes are close to each other. Comparison of the sampling methods shows that the Markov chain approach yields more realistic observables at both high and low FRET efficiencies. Differences between the three sampling methods are discussed in terms of the different mechanisms for averaging over structural dynamics in the system. Convergence of the Markov chain method is carefully examined. Together, the methods for estimating coupling and for sampling the coupling provide a mechanism for directly connecting the structural dynamics modeled by MD with fluorescence observables determined through FRET experiments.

Fretting about FRET: correlation between kappa and R.

Molecular dynamics simulations were used to examine the structural dynamics of two fluorescent probes attached to a typical protein, hen egg-white lysozyme (HEWL). The donor probe (D) was attached via a succinimide group, consistent with the commonly-used maleimide conjugation chemistry, and the acceptor probe (A) was bound into the protein as occurs naturally for HEWL and the dye Eosin Y. The is found to deviate significantly from the theoretical value and high correlation between the orientation factor kappa and the distance R is observed. The correlation is quantified using several possible fixed A orientations and correlation as high as 0.80 is found between kappa and R and as high as 0.68 between kappa(2) and R. The presence of this correlation highlights the fact that essentially all fluorescence-detected resonance energy transfer studies have assumed that kappa and R are independent--an assumption that is clearly not justified in the system studied here. The correlation results in the quantities and < R(-)(6)> differing by a factor of 1.6. The observed correlation between kappa and R is caused by the succinimide linkage between the D and HEWL, which is found to be relatively inflexible.