On the Interpretation of subtilisin Carlsberg Time-Resolved Fluorescence Anisotropy Decays: Modeling with Classical Simulations.

In this work, we discuss the challenging time-resolved fluorescence anisotropy of subtilisin Carlsberg (SC), which contains a single Trp residue and is a model fluorescence system. Experimental decay rates and quenching data suggest that the fluorophore should be exposed to water, but the Trp is partially buried in a hydrophobic pocket in the crystallographic structure. In order to study this inconsistency, molecular dynamics simulations were performed to predict the anisotropy decay rates and emission wavelengths of the Trp. We confirmed the inconsistency of the crystallographic structure with the experimentally observed fluorescence data and performed free energy calculations to show that the buried Trp conformation is 2 orders of magnitude (∼3 kcal/mol) more stable than the solvent-exposed one. However, molecular dynamics simulations in which the Trp side chain was restricted to solvent-exposed conformations displayed a maximum Trp emission wavelength shifted toward lower energies and decay rates compatible with the experimentally probed rates. Therefore, if the solvent-exposed conformations are the most important emitters, the experimental anisotropy can be compatibilized with the crystallographic structure. The most likely explanation is that the fluorescence of the most probable conformation in solution, observed in the crystal, is quenched, and this is consistent with the low quantum yield of Trp113 of SC. Additionally, some experiments might have probed denatured or lysed SC structures. SC anisotropy provides an interesting target for the study of fluorescence anisotropy using simulations, which can be used to test and exemplify how modeling can aid the interpretation of experimental data in a system where structure and solution experiments appear to be inconsistent.

Parametric models to compute tryptophan fluorescence wavelengths from classical protein simulations.

Fluorescence spectroscopy is an important method to study protein conformational dynamics and solvation structures. Tryptophan (Trp) residues are the most important and practical intrinsic probes for protein fluorescence due to the variability of their fluorescence wavelengths: Trp residues emit in wavelengths ranging from 308 to 360 nm depending on the local molecular environment. Fluorescence involves electronic transitions, thus its computational modeling is a challenging task. We show that it is possible to predict the wavelength of emission of a Trp residue from classical molecular dynamics simulations by computing the solvent-accessible surface area or the electrostatic interaction between the indole group and the rest of the system. Linear parametric models are obtained to predict the maximum emission wavelengths with standard errors of the order 5 nm. In a set of 19 proteins with emission wavelengths ranging from 308 to 352 nm, the best model predicts the maximum wavelength of emission with a standard error of 4.89 nm and a quadratic Pearson correlation coefficient of 0.81. These models can be used for the interpretation of fluorescence spectra of proteins with multiple Trp residues, or for which local Trp environmental variability exists and can be probed by classical molecular dynamics simulations. © 2018 Wiley Periodicals, Inc.