Damped Response Theory in Combination with Polarizable Environments: The Polarizable Embedding Complex Polarization Propagator Method.

We present a combination of the polarizable embedding (PE) scheme with the complex polarization propagator (CPP) method with the aim of calculating response properties including relaxation for large and complex systems. This new approach, termed PE-CPP, will benefit from the highly advanced description of the environmental electrostatic potential and polarization in the PE method as well as the treatment of near-resonant effects in the CPP approach. The PE-CPP model has been implemented in a Kohn-Sham density functional theory approach, and we present pilot calculations exemplifying the implementation for the UV/vis and carbon K-edge X-ray absorption spectra of the protein plastocyanin. Furthermore, technical details associated with a PE-CPP calculation are discussed.

Photoabsorption of acridine yellow and proflavin bound to human serum albumin studied by means of quantum mechanics/molecular dynamics.

Attempting to unravel mechanisms in optical probing of proteins, we have performed pilot calculations of two cationic chromophores-acridine yellow and proflavin-located at different binding sites within human serum albumin, including the two primary drug binding sites as well as a heme binding site. The computational scheme adopted involves classical molecular dynamics simulations of the ligands bound to the protein and subsequent linear response polarizable embedding density functional theory calculations of the excitation energies. A polarizable embedding potential consisting of point charges fitted to reproduce the electrostatic potential and isotropic atomic polarizabilities computed individually for every residue of the protein was used in the linear response calculations. Comparing the calculated aqueous solution-to-protein shifts of maximum absorption energies to available experimental data, we concluded that the cationic proflavin chromophore is likely not to bind albumin at its drug binding site 1 nor at its heme binding site. Although agreement with experimental data could only be obtained in qualitative terms, our results clearly indicate that the difference in optical response of the two probes is due to deprotonation, and not, as earlier suggested, to different binding sites. The ramifications of this finding for design of molecular probes targeting albumin or other proteins is briefly discussed.

Molecular-Level Insight into the Spectral Tuning Mechanism of the DsRed Chromophore.

We present a detailed study of the protein environmental effects on the one- and two-photon absorption (1PA and 2PA, respectively) properties of the S0-S1 transition in the DsRed protein using the polarizable embedding density functional theory formalism. We find that steric factors and chromophore-protein interactions act in concert to enhance the 2PA activity inside the protein while adversely blue-shifting the 1PA maximum. A two-state model reveals that the 2PA intensity gain is primarily governed by the increased change in the permanent dipole moment between the ground and the excited states acquired inside the protein. Our results indicate that this mainly is attributable to counter-directional contributions stemming from Lys163 and the conserved Arg95 with the former additionally identified as a key residue in the color tuning mechanism. The results provide new insight into the tuning mechanism of DsRed and suggest a possible strategy for simultaneous improvement of its 1PA and 2PA properties.

Computational protocols for prediction of solute NMR relative chemical shifts. a case study of L-tryptophan in aqueous solution.

In this study, we have applied two different spanning protocols for obtaining the molecular conformations of L-tryptophan in aqueous solution, namely a molecular dynamics simulation and a molecular mechanics conformational search with subsequent geometry re-optimization of the stable conformers using a quantum mechanically based method. These spanning protocols represent standard ways of obtaining a set of conformations on which NMR calculations may be performed. The results stemming from the solute-solvent configurations extracted from the MD simulation at 300 K are found to be inferior to the results stemming from the conformations extracted from the MM conformational search in terms of replicating an experimental reference as well as in achieving the correct sequence of the NMR relative chemical shifts of L-tryptophan in aqueous solution. We find this to be due to missing conformations visited during the molecular dynamics run as well as inaccuracies in geometrical parameters generated from the classical molecular dynamics simulations.