A polarizable QM/MM approach to the molecular dynamics of amide groups solvated in water.

The infrared (IR) spectra of polypeptides are dominated by the so-called amide bands. Because they originate from the strongly polar and polarizable amide groups (AGs) making up the backbone, their spectral positions sensitively depend on the local electric fields. Aiming at accurate computations of these IR spectra by molecular dynamics (MD) simulations, which derive atomic forces from a hybrid quantum and molecular mechanics (QM/MM) Hamiltonian, here we consider the effects of solvation in bulk liquid water on the amide bands of the AG model compound N-methyl-acetamide (NMA). As QM approach to NMA we choose grid-based density functional theory (DFT). For the surrounding MM water, we develop, largely based on computations, a polarizable molecular mechanics (PMM) model potential called GP6P, which features six Gaussian electrostatic sources (one induced dipole, five static partial charge distributions) and, therefore, avoids spurious distortions of the DFT electron density in hybrid DFT/PMM simulations. Bulk liquid GP6P is shown to have favorable properties at the thermodynamic conditions of the parameterization and beyond. Lennard-Jones (LJ) parameters of the DFT fragment NMA are optimized by comparing radial distribution functions in the surrounding GP6P liquid with reference data obtained from a "first-principles" DFT-MD simulation. Finally, IR spectra of NMA in GP6P water are calculated from extended DFT/PMM-MD trajectories, in which the NMA is treated by three different DFT functionals (BP, BLYP, B3LYP). Method-specific frequency scaling factors are derived from DFT-MD simulations of isolated NMA. The DFT/PMM-MD simulations with GP6P and with the optimized LJ parameters then excellently predict the effects of aqueous solvation and deuteration observed in the IR spectra of NMA. As a result, the methods required to accurately compute such spectra by DFT/PMM-MD also for larger peptides in aqueous solution are now at hand.

Linearly scaling and almost Hamiltonian dielectric continuum molecular dynamics simulations through fast multipole expansions.

Hamiltonian Dielectric Solvent (HADES) is a recent method [S. Bauer et al., J. Chem. Phys. 140, 104103 (2014)] which enables atomistic Hamiltonian molecular dynamics (MD) simulations of peptides and proteins in dielectric solvent continua. Such simulations become rapidly impractical for large proteins, because the computational effort of HADES scales quadratically with the number N of atoms. If one tries to achieve linear scaling by applying a fast multipole method (FMM) to the computation of the HADES electrostatics, the Hamiltonian character (conservation of total energy, linear, and angular momenta) may get lost. Here, we show that the Hamiltonian character of HADES can be almost completely preserved, if the structure-adapted fast multipole method (SAMM) as recently redesigned by Lorenzen et al. [J. Chem. Theory Comput. 10, 3244-3259 (2014)] is suitably extended and is chosen as the FMM module. By this extension, the HADES/SAMM forces become exact gradients of the HADES/SAMM energy. Their translational and rotational invariance then guarantees (within the limits of numerical accuracy) the exact conservation of the linear and angular momenta. Also, the total energy is essentially conserved-up to residual algorithmic noise, which is caused by the periodically repeated SAMM interaction list updates. These updates entail very small temporal discontinuities of the force description, because the employed SAMM approximations represent deliberately balanced compromises between accuracy and efficiency. The energy-gradient corrected version of SAMM can also be applied, of course, to MD simulations of all-atom solvent-solute systems enclosed by periodic boundary conditions. However, as we demonstrate in passing, this choice does not offer any serious advantages.

Utilizing fast multipole expansions for efficient and accurate quantum-classical molecular dynamics simulations.

Recently, a novel approach to hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations has been suggested [Schwörer et al., J. Chem. Phys. 138, 244103 (2013)]. Here, the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 10(3)-10(5) molecules as negative gradients of a DFT/PMM hybrid Hamiltonian. The electrostatic interactions are efficiently described by a hierarchical fast multipole method (FMM). Adopting recent progress of this FMM technique [Lorenzen et al., J. Chem. Theory Comput. 10, 3244 (2014)], which particularly entails a strictly linear scaling of the computational effort with the system size, and adapting this revised FMM approach to the computation of the interactions between the DFT and PMM fragments of a simulation system, here, we show how one can further enhance the efficiency and accuracy of such DFT/PMM-MD simulations. The resulting gain of total performance, as measured for alanine dipeptide (DFT) embedded in water (PMM) by the product of the gains in efficiency and accuracy, amounts to about one order of magnitude. We also demonstrate that the jointly parallelized implementation of the DFT and PMM-MD parts of the computation enables the efficient use of high-performance computing systems. The associated software is available online.

Electrostatics of proteins in dielectric solvent continua. I. An accurate and efficient reaction field description.

We present a reaction field (RF) method which accurately solves the Poisson equation for proteins embedded in dielectric solvent continua at a computational effort comparable to that of an electrostatics calculation with polarizable molecular mechanics (MM) force fields. The method combines an approach originally suggested by Egwolf and Tavan [J. Chem. Phys. 118, 2039 (2003)] with concepts generalizing the Born solution [Z. Phys. 1, 45 (1920)] for a solvated ion. First, we derive an exact representation according to which the sources of the RF potential and energy are inducible atomic anti-polarization densities and atomic shielding charge distributions. Modeling these atomic densities by Gaussians leads to an approximate representation. Here, the strengths of the Gaussian shielding charge distributions are directly given in terms of the static partial charges as defined, e.g., by standard MM force fields for the various atom types, whereas the strengths of the Gaussian anti-polarization densities are calculated by a self-consistency iteration. The atomic volumes are also described by Gaussians. To account for covalently overlapping atoms, their effective volumes are calculated by another self-consistency procedure, which guarantees that the dielectric function ε(r) is close to one everywhere inside the protein. The Gaussian widths σ(i) of the atoms i are parameters of the RF approximation. The remarkable accuracy of the method is demonstrated by comparison with Kirkwood's analytical solution for a spherical protein [J. Chem. Phys. 2, 351 (1934)] and with computationally expensive grid-based numerical solutions for simple model systems in dielectric continua including a di-peptide (Ac-Ala-NHMe) as modeled by a standard MM force field. The latter example shows how weakly the RF conformational free energy landscape depends on the parameters σ(i). A summarizing discussion highlights the achievements of the new theory and of its approximate solution particularly by comparison with so-called generalized Born methods. A follow-up paper describes how the method enables Hamiltonian, efficient, and accurate MM molecular dynamics simulations of proteins in dielectric solvent continua.

Electrostatics of proteins in dielectric solvent continua. II. Hamiltonian reaction field dynamics.

In Paper I of this work [S. Bauer, G. Mathias, and P. Tavan, J. Chem. Phys. 140, 104102 (2014)] we have presented a reaction field (RF) method, which accurately solves the Poisson equation for proteins embedded in dielectric solvent continua at a computational effort comparable to that of polarizable molecular mechanics (MM) force fields. Building upon these results, here we suggest a method for linearly scaling Hamiltonian RF/MM molecular dynamics (MD) simulations, which we call "Hamiltonian dielectric solvent" (HADES). First, we derive analytical expressions for the RF forces acting on the solute atoms. These forces properly account for all those conditions, which have to be self-consistently fulfilled by RF quantities introduced in Paper I. Next we provide details on the implementation, i.e., we show how our RF approach is combined with a fast multipole method and how the self-consistency iterations are accelerated by the use of the so-called direct inversion in the iterative subspace. Finally we demonstrate that the method and its implementation enable Hamiltonian, i.e., energy and momentum conserving HADES-MD, and compare in a sample application on Ac-Ala-NHMe the HADES-MD free energy landscape at 300 K with that obtained in Paper I by scanning of configurations and with one obtained from an explicit solvent simulation.

Piperazine derivatives inhibit PrP/PrP(res) propagation in vitro and in vivo.

Prion diseases are fatal neurodegenerative disorders, which are not curable and no effective treatment exists so far. The major neuropathological change in diseased brains is the conversion of the normal cellular form of the prion protein PrPc(C) into a disease-associated isoform PrP(Sc). PrP(Sc) accumulates into multimeres and fibrillar aggregates, which leads to the formation of amyloid plaques. Increasing evidence indicates a fundamental role of PrP(Sc) species and its aggregation in the pathogenesis of prion diseases, which initiates the pathological cascade and leads to neurodegeneration accompanied by spongiform changes. In search of compounds that have the potential to interfere with PrP(Sc) formation and propagation, we used a cell based assay for the screening of potential aggregation inhibitors. The assay deals with a permanently prion infected cell line that was adapted for a high-throughput screening of a compound library composed of 10,000 compounds (DIVERset 2, ChemBridge). We could detect six different classes of highly potent inhibitors of PrP(Sc) propagation in vitro and identified piperazine derivatives as a new inhibitory lead structure, which increased incubation time of scrapie infected mice.

Polarizable six-point water models from computational and empirical optimization.

Tröster et al. (J. Phys. Chem B 2013, 117, 9486-9500) recently suggested a mixed computational and empirical approach to the optimization of polarizable molecular mechanics (PMM) water models. In the empirical part the parameters of Buckingham potentials are optimized by PMM molecular dynamics (MD) simulations. The computational part applies hybrid calculations, which combine the quantum mechanical description of a H2O molecule by density functional theory (DFT) with a PMM model of its liquid phase environment generated by MD. While the static dipole moments and polarizabilities of the PMM water models are fixed at the experimental gas phase values, the DFT/PMM calculations are employed to optimize the remaining electrostatic properties. These properties cover the width of a Gaussian inducible dipole positioned at the oxygen and the locations of massless negative charge points within the molecule (the positive charges are attached to the hydrogens). The authors considered the cases of one and two negative charges rendering the PMM four- and five-point models TL4P and TL5P. Here we extend their approach to three negative charges, thus suggesting the PMM six-point model TL6P. As compared to the predecessors and to other PMM models, which also exhibit partial charges at fixed positions, TL6P turned out to predict all studied properties of liquid water at p0 = 1 bar and T0 = 300 K with a remarkable accuracy. These properties cover, for instance, the diffusion constant, viscosity, isobaric heat capacity, isothermal compressibility, dielectric constant, density, and the isobaric thermal expansion coefficient. This success concurrently provides a microscopic physical explanation of corresponding shortcomings of previous models. It uniquely assigns the failures of previous models to substantial inaccuracies in the description of the higher electrostatic multipole moments of liquid phase water molecules. Resulting favorable properties concerning the transferability to other temperatures and conditions like the melting of ice are also discussed.

The Microscopic Physical Cause for the Density Maximum of Liquid Water.

The existence of a density maximum at 277 K is probably the most prominent anomaly among the many very special thermodynamic properties of liquid water. While usually attributed to so-called hydrogen bonding, the microscopic physical cause of this prominent anomaly is still elusive. Here we show that the density anomaly is caused by those short-range electrostatic forces, which are generated by the quadrupole and higher moments of the charge distributions present in liquid-phase water molecules. This conclusion derives from 20 ns replica exchange molecular-dynamics simulations with closely related polarizable four-, five-, and six-point water models. As soon as the model complexity suffices to represent the higher electrostatic moments with sufficient accuracy, the density temperature profile n(T) calculated for T ∈ [250,320] K at the standard pressure 1 bar locks in to the experimental observation. The corresponding six-point model is, therefore, the most simple available cartoon for liquid-phase water molecules.

Coupling density functional theory to polarizable force fields for efficient and accurate Hamiltonian molecular dynamics simulations.

Hybrid molecular dynamics (MD) simulations, in which the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 10(3)-10(5) molecules, pose a challenge. A corresponding computational approach should guarantee energy conservation, exclude artificial distortions of the electron density at the interface between the DFT and PMM fragments, and should treat the long-range electrostatic interactions within the hybrid simulation system in a linearly scaling fashion. Here we describe a corresponding Hamiltonian DFT/(P)MM implementation, which accounts for inducible atomic dipoles of a PMM environment in a joint DFT/PMM self-consistency iteration. The long-range parts of the electrostatics are treated by hierarchically nested fast multipole expansions up to a maximum distance dictated by the minimum image convention of toroidal boundary conditions and, beyond that distance, by a reaction field approach such that the computation scales linearly with the number of PMM atoms. Short-range over-polarization artifacts are excluded by using Gaussian inducible dipoles throughout the system and Gaussian partial charges in the PMM region close to the DFT fragment. The Hamiltonian character, the stability, and efficiency of the implementation are investigated by hybrid DFT/PMM-MD simulations treating one molecule of the water dimer and of bulk water by DFT and the respective remainder by PMM.

Polarizable water models from mixed computational and empirical optimization.

Here we suggest a mixed computational and empirical approach serving to optimize the parameters of complex and polarizable molecular mechanics (PMM) models for complicated liquids. The computational part of the parameter optimization relies on hybrid calculations combining density functional theory (DFT) for a solute molecule with a PMM treatment of its solvent environment at well-defined thermodynamic conditions. As an application we have developed PMM models for water featuring ν = 3, 4, and 5 points of force action, a Gaussian inducible dipole and a Buckingham potential at the oxygen, the experimental liquid phase geometry, the experimental gas phase polarizability α(exp)(g) = 1.47 ų, and, for ν = 4 and 5, the gas phase value µ(exp)(g) = 1.855 D for the static dipole moment. The widths of the Gaussian dipoles and, for ν = 4 and 5, also the electrostatic geometries of these so-called TLνP models are derived from self-consistent DFT/PMM calculations, and the parameters of the Buckingham potentials (and the static TL3P dipole moment) are estimated from molecular dynamics (MD) simulations. The high quality of the resulting models is demonstrated for the observables targeted during optimization (potential energy per molecule, pressure, radial distribution functions) and a series of predicted properties (quadrupole moments, density at constant pressure, dielectric constant, diffusivity, viscosity, compressibility, heat capacity) at certain standard conditions. Remaining deficiencies and possible ways for their removal are discussed.

Anle138b: a novel oligomer modulator for disease-modifying therapy of neurodegenerative diseases such as prion and Parkinson's disease.

In neurodegenerative diseases such as Alzheimer's disease (AD), Parkinson's disease (PD) and prion diseases, deposits of aggregated disease-specific proteins are found. Oligomeric aggregates are presumed to be the key neurotoxic agent. Here we describe the novel oligomer modulator anle138b [3-(1,3-benzodioxol-5-yl)-5-(3-bromophenyl)-1H-pyrazole], an aggregation inhibitor we developed based on a systematic high-throughput screening campaign combined with medicinal chemistry optimization. In vitro, anle138b blocked the formation of pathological aggregates of prion protein (PrP(Sc)) and of α-synuclein (α-syn), which is deposited in PD and other synucleinopathies such as dementia with Lewy bodies (DLB) and multiple system atrophy (MSA). Notably, anle138b strongly inhibited all prion strains tested including BSE-derived and human prions. Anle138b showed structure-dependent binding to pathological aggregates and strongly inhibited formation of pathological oligomers in vitro and in vivo both for prion protein and α-synuclein. Both in mouse models of prion disease and in three different PD mouse models, anle138b strongly inhibited oligomer accumulation, neuronal degeneration, and disease progression in vivo. Anle138b had no detectable toxicity at therapeutic doses and an excellent oral bioavailability and blood-brain-barrier penetration. Our findings indicate that oligomer modulators provide a new approach for disease-modifying therapy in these diseases, for which only symptomatic treatment is available so far. Moreover, our findings suggest that pathological oligomers in neurodegenerative diseases share structural features, although the main protein component is disease-specific, indicating that compounds such as anle138b that modulate oligomer formation by targeting structure-dependent epitopes can have a broad spectrum of activity in the treatment of different protein aggregation diseases.

Electronic excitations in long polyenes revisited.

We apply the valence shell model OM2 [W. Weber and W. Thiel, Theor. Chem. Acc. 103, 495, (2000)] combined with multireference configuration interaction (MRCI) to compute the vertical excitation energies and transition dipole moments of the low-energy singlet excitations in the polyenes with 4 ≤ N ≤ 22π-electrons. We find that the OM2/MRCI descriptions closely resemble those of Pariser-Parr-Pople (PPP) π-electron models [P. Tavan and K. Schulten, Phys. Rev. B 36, 4337, (1987)], if equivalent MRCI procedures and regularly alternating model geometries are used. OM2/MRCI optimized geometries are shown to entail improved descriptions particularly for smaller polyenes (N ≤ 12), for which sizeable deviations from the regular model geometries are found. With configuration interaction active spaces covering also the σ- in addition to the π-electrons, OM2/MRCI excitation energies turn out to become smaller by at most 0.35 eV for the ionic and 0.15 eV for the covalent excitations. The particle-hole (ph) symmetry, which in Pariser-Parr-Pople models arises from the zero-differential overlap approximation, is demonstrated to be only weakly broken in OM2 such that the oscillator strengths of the covalent 1B(u)(-) states, which artificially vanish in ph-symmetric models, are predicted to be very small. According to OM2/MRCI and experimental data the 1B(u)(-) state is the third excited singlet state for N < 12 and becomes the second for N ≥ 14. By comparisons with results of other theoretical approaches and experimental evidence we argue that deficiencies of the particular MRCI method employed by us, which show up in a poor size consistency of the covalent excitations for N > 12, are caused by its restriction to at most doubly excited references.

Vibrational spectra of phosphate ions in aqueous solution probed by first-principles molecular dynamics.

We have carried out "first-principles" Born-Oppenheimer molecular dynamics (BOMD) simulations of the phosphate ions H2PO4⁻ and HPO4²â» in liquid water and have calculated their IR spectra by Fourier transform techniques from the trajectories. IR bands were assigned by a so-called "generalized normal coordinate analysis". The effects of including Hartree-Fock (HF) exchange into the density functional theory (DFT) computation of forces were studied by comparing results obtained with the well-known BP, BLYP, and B3LYP functionals. The neglect of dispersion in the functionals was empirically corrected. The inclusion of HF exchange turned out to yield dramatically improved and, thus, quite accurate descriptions of the IR spectra observed for H2PO4⁻ and HPO4²â» in aqueous solution. An analysis of earlier computational results (Klähn, M. et al. J. Phys. Chem. A 2004, 108, 6186-6194) on these vibrational spectra, which had been obtained in a hybrid setting combining a BP description of the respective phosphate with a simple molecular mechanics (MM) model of its aqueous environment, revealed three different sources of error, (i) the BP force field of the phosphates is much too soft and would have required a substantial scaling of frequencies, (ii) the oversimplified water force field entailed incorrect solvation structures and, thus, qualitatively wrong patterns of solvatochromic band shifts, and (iii) quantitative frequency computations additionally required the inclusion of HF exchange. Thus, the results of the B3LYP BOMD simulations do not only characterize physical properties like the IR spectra or the solvation structures of the phosphate systems but also provide clues for the future design of simplified but nevertheless reasonably accurate DFT/MM methods applicable to phosphates.

Optimizing the Accuracy and Efficiency of Fast Hierarchical Multipole Expansions for MD Simulations.

Based on p'th order Cartesian Taylor expansions of Coulomb interactions and on hierarchical decompositions of macromolecular simulation systems into hierarchies of nested, structure-adapted, and adaptively formed clusters of increasing size, fast multipole methods are constructed for rapid and accurate calculations of electrostatic interactions. These so-called SAMMp algorithms are formulated through totally symmetric and traceless tensors describing the multipole moments and the coefficients of local Taylor expansions. Simple recursions for the efficient evaluation and shifting of multipole moments are given. The required tensors are explicitly given up to order p = 4. The SAMMp algorithms are shown to guarantee the reaction principle. For systems with periodic boundaries, a reaction field (RF) correction is applied, which introduces at distances beyond the "minimum image convention" boundary a dielectric continuum surrounding each cluster at the top level of coarse graining. The correctness of the present SAMMp implementation is demonstrated by analyzing the scaling of the residuals and by checking the numerical accuracy of the reaction principle for a pair of distant molecular ions in vacuum. Molecular dynamics simulations of pure water and aqueous solutions containing artificial ions, which are enclosed by periodic boundaries, demonstrate the stability and low-noise behavior of SAMMp/RF.

DFT/MM description of flavin IR spectra in BLUF domains.

A class of photoreceptors occurring in various organisms consists of domains that are blue light sensing using flavin (BLUF). The vibrational spectra of the flavin chromophore are spectroscopically well characterized for the dark-adapted resting states and for the light-adapted signaling states of BLUF domains in solution. Here we present a theoretical analysis of such spectra by applying density functional theory (DFT) to the flavin embedded in molecular mechanics (MM) models of its protein and solvent environment. By DFT/MM we calculate flavin spectra for seven different X-ray and NMR structures of BLUF domains occurring in the transcriptional antirepressor AppA and in the blue light receptor B (BlrB) of the purple bacterium Rb. Sphaeroides as well as in the phototaxis photoreceptor Slr1694 of the cyanobacterium Synechocystis. By considering the dynamical stabilities of associated all-atom simulation models and by comparing calculated with observed vibrational spectra, we show that two of the considered structures (both AppA) are obviously erroneous and that specific features of two further crystal structures (BlrB and Slr1694) cannot represent the states of the respective BLUF domains in solution. Thereby, the conformational transitions elicited by solvation are identified. In this context we demonstrate how hydrogen bonds of varying strengths can tune in BLUF domains the C═O stretching frequencies of the flavin chromophore. Furthermore we show that the DFT/MM spectra of the flavin calculated for two different AppA BLUF conformations, which are called Trp(in) and Met(in), fit very well to the spectroscopic data observed for the dark and light states, respectively, if (i) polarized MM force fields, which are calculated by an iterative DFT/MM procedure, are employed for the flavin binding pockets and (ii) the calculated frequencies are properly scaled. Although the associated analysis indicates that the Trp(in) conformation belongs to the dark state, no clear light vs dark distinction emerges for the Met(in) conformation. In this connection, a number of methodological issues relevant for such complex computations are thoroughly discussed showing, in particular, why our current descriptions could not decide the light vs dark question for Met(in).

From high-throughput cell culture screening to mouse model: identification of new inhibitor classes against prion disease.

Transmissible spongiform encephalopathies (TSE) or prion diseases belong to a category of fatal and so far untreatable neurodegenerative conditions. All prion diseases are characterized by both degeneration in the central nervous system (CNS) in humans and animals and the deposition and accumulation of Proteinase K-resistant prion protein (PrP(res)). Until now, no pharmaceutical product has been available to cure these diseases or to alleviate their associated symptoms. Here, a cell-culture screening system is described that allows for the large-scale analysis of the PrP(res) inhibitory potential of a library of compounds and the identification of structural motifs leading potent compounds able to cause PrP(res) clearance at the cellular level. Based on different scrapie-infected cell lines, 10,000 substances were tested, out of which 530 potential inhibitors were identified. After re-screening and validation using a series of dilutions, 14 compounds were identified as the most effective. These 14 compounds were then used for therapeutic studies in a mouse bioassay to test and verify their in vivo potency. Two compounds exhibited therapeutic potential in the mouse model by significantly extending the survival time of intracerebrally infected mice, when treated 90 days after infection with scrapie.

Diphenylpyrazole-derived compounds increase survival time of mice after prion infection.

Transmissible spongiform encephalopathies (TSEs) represent a group of fatal neurodegenerative disorders that can be transmitted by natural infection or inoculation. TSEs include scrapie in sheep, bovine spongiform encephalopathy (BSE) in cattle, and Creutzfeldt-Jakob disease (CJD) in humans. The emergence of a variant form of CJD (vCJD), which has been associated with BSE, produced strong pressure to search for effective treatments with new drugs. Up to now, however, TSEs have proved incurable, although many efforts have been made both in vitro and in vivo to search for potent therapeutic and prophylactic compounds. For this purpose, we analyzed a compound library consisting of 10,000 compounds with a cell-based high-throughput screening assay dealing with scrapie-infected scrapie mouse brain and ScN(2)A cells and identified a new class of inhibitors consisting of 3,5-diphenylpyrazole (DPP) derivatives. The most effective DPP derivative showed half-maximal inhibition of PrP(Sc) formation at concentrations (IC(50)) of 0.6 and 1.2 µM, respectively. This compound was subsequently subjected to a number of animal experiments using scrapie-infected wild-type C57BL/6 and transgenic Tga20 mice. The DPP derivative induced a significant increase of incubation time both in therapeutic and prophylactic experiments. The onset of the prion disease was delayed by 37 days after intraperitoneal and 42 days after oral application, respectively. In summary, we demonstrate a high in vitro efficiency of DPP derivatives against prion infections that was substantiated in vivo for one of these compounds. These results indicate that the novel class of DPP compounds should comprise excellent candidates for future therapeutic studies.

IR spectra of flavins in solution: DFT/MM description of redox effects.

The functional reactions in blue light photoreceptors generally involve transiently reduced flavins exhibiting characteristic infrared (IR) spectra. To approach a theoretical understanding, here we apply density functional theory (DFT) to flavin radicals embedded in a molecular mechanics (MM) model of an aqueous solution. Combining a DFT/MM approach with instantaneous normal-mode analyses (INMA), we compute the IR solution spectra of anionic and neutral flavin radicals. For a set of mid-IR marker bands, we identify those changes of spectral locations, intensities, and widths, which are caused by sequentially adding an electron and a proton to the oxidized flavin. Comparisons with experimental IR solution spectra of flavin radicals show the accuracy of our DFT/MM-INMA approach and allow us to assign the observed bands. The room temperature ensembles of solvent cages required for the INMA calculations of the IR spectra are generated in an MM setting from molecular dynamics (MD) simulations. For the solvated flavin radicals, these MD simulations employ MM force fields derived from DFT/MM calculations.

Density functional theory combined with molecular mechanics: the infrared spectra of flavin in solution.

The photophysics and photochemistry of flavin dyes determine the functional dynamics of a series of blue light photoreceptors that include the so-called BLUF (blue light sensors using flavin) domains. To enable molecular dynamics (MD) simulation studies of such signaling processes, we derived molecular mechanics (MM) models of flavin chromophores from density functional theory (DFT). Two 300 K ensembles of lumiflavin (LF) in aqueous solution were generated by extended MM-MD simulations using different MM potentials for the water. In a DFT/MM hybrid setting, in which LF was treated by DFT and the polarizing environment at atomistic resolution by MM, we applied instantaneous normal mode analyses (INMA) to these ensembles. From these data we determined the inhomogeneously broadened solution spectra as mixtures of Gaussian bands using a novel automated procedure for mode classification. Comparisons with vibrational spectra available in the literature on native and isotopically labeled flavins in aqueous solution serve us to determine suitable frequency scaling factors and to analyze the accuracy of our scaled DFT/MM-INMA approach. We show that our approach not only agrees with established computational descriptions but also extends such methods substantially by giving access to inhomogeneous line widths and band shapes.