Mary A. Rohrdanz

A long-range-corrected density functional that performs well for both ground-state properties and time-dependent density functional theory excitation energies, including charge-transfer excited states

We introduce a hybrid density functional that asymptotically incorporates full Hartree-Fock exchange, based on the long-range-corrected exchange-hole model of Henderson et al. [J. Chem. Phys. 128, 194105 (2008)]. The performance of this functional, for ground-state properties and for vertical excitation energies within time-dependent density functional theory, is systematically evaluated, and optimal values are determined for the range-separation parameter, omega, and for the fraction of short-range Hartree-Fock exchange. We denote the new functional as LRC-omegaPBEh, since it reduces to the standard PBEh hybrid functional (also known as PBE0 or PBE1PBE) for a certain choice of its two parameters. Upon optimization of these parameters against a set of ground- and excited-state benchmarks, the LRC-omegaPBEh functional fulfills three important requirements: (i) It outperforms the PBEh hybrid functional for ground-state atomization energies and reaction barrier heights; (ii) it yields statistical errors comparable to PBEh for valence excitation energies in both small and medium-sized molecules; and (iii) its performance for charge-transfer excitations is comparable to its performance for valence excitations. LRC-omegaPBEh, with the parameters determined herein, is the first density functional that satisfies all three criteria. Notably, short-range Hartree-Fock exchange appears to be necessary in order to obtain accurate ground-state properties and vertical excitation energies using the same value of omega.

Simultaneous benchmarking of ground- and excited-state properties with long-range-corrected density functional theory

We present benchmark calculations using several long-range-corrected (LRC) density functionals, in which Hartree-Fock exchange is incorporated asymptotically using a range-separated Coulomb operator, while local exchange is attenuated using an ansatz introduced by Iikura et al. [J. Chem. Phys. 115, 3540 (2001)]. We calculate ground-state atomization energies, reaction barriers, ionization energies, and electron affinities, each as a function of the range-separation parameter mu. In addition, we calculate excitation energies of small- and medium-sized molecules, again as a function of mu, by applying the LRC to time-dependent density functional theory. Representative examples of both pure and hybrid density functionals are tested. On the basis of these results, there does not appear to be a single range-separation parameter that is reasonable for both ground-state properties and vertical excitation energies. Reasonable errors in atomization energies and barrier heights are achieved only at the expense of excessively high excitation energies, at least for the medium-sized molecules, whereas values of mu that afford reasonable excitation energies yield some of the largest errors for ground-state atomization energies and barrier heights in small molecules. Notably, this conclusion is obscured if the database of excitation energies includes only small molecules, as has been the case in previous benchmark studies of LRC functionals.

Determination of reaction coordinates via locally scaled diffusion map.

We present a multiscale method for the determination of collective reaction coordinates for macromolecular dynamics based on two recently developed mathematical techniques: diffusion map and the determination of local intrinsic dimensionality of large datasets. Our method accounts for the local variation of molecular configuration space, and the resulting global coordinates are correlated with the time scales of the molecular motion. To illustrate the approach, we present results for two model systems: all-atom alanine dipeptide and coarse-grained src homology 3 protein domain. We provide clear physical interpretation for the emerging coordinates and use them to calculate transition rates. The technique is general enough to be applied to any system for which a Boltzmann-sampled set of molecular configurations is available.

Charge-Transfer Excited States in a π-Stacked Adenine Dimer, As Predicted Using Long-Range-Corrected Time-Dependent Density Functional Theory

The lowest few electronic excitations of a pi-stacked adenine dimer in its B-DNA geometry are investigated, in the gas phase and in a water cluster, using a long-range-corrected version of time-dependent density functional theory (TD-DFT) that asymptotically incorporates Hartree-Fock exchange. Long-range correction is shown to eliminate the catastrophic underestimation of charge-transfer (CT) excitation energies that plagues conventional TD-DFT, at the expense of introducing one adjustable parameter, mu, that determines the length scale on which Hartree-Fock exchange is turned on. This parameter allows us to interpolate smoothly between hybrid density functionals and time-dependent Hartree-Fock theory. Excitation energies for CT states (in which an electron is transferred from one adenine molecule to the other) are found to increase dramatically as a function of mu. Uncorrected hybrid functionals underestimate the CT excitation energies, placing them well below the valence excitations, while time-dependent Hartree-Fock calculations place these states well above the valence states. Values for mu determined from certain benchmark calculations place the CT states well above the valence pipi* and npi* states at the Franck-Condon point.

Discovering Mountain Passes via Torchlight: Methods for the Definition of Reaction Coordinates and Pathways in Complex Macromolecular Reactions

The long-timescale dynamics of macromolecular systems can be oftentimes viewed as a reaction connecting metastable states of the system. In the past decade, various approaches have been developed to discover the collective motions associated with this dynamics. The corresponding collective variables are used in many applications, e.g., to understand the reaction mechanism, to quantify the systems free energy landscape, to enhance the sampling of the reaction path, and to determine the reaction rate. In this review we focus on a number of key developments in this field, providing an overview of several methods along with their relative regimes of applicability.

Rapid exploration of configuration space with diffusion-map-directed molecular dynamics.

The gap between the time scale of interesting behavior in macromolecular systems and that which our computational resources can afford often limits molecular dynamics (MD) from understanding experimental results and predicting what is inaccessible in experiments. In this paper, we introduce a new sampling scheme, named diffusion-map-directed MD (DM-d-MD), to rapidly explore molecular configuration space. The method uses a diffusion map to guide MD on the fly. DM-d-MD can be combined with other methods to reconstruct the equilibrium free energy, and here, we used umbrella sampling as an example. We present results from two systems: alanine dipeptide and alanine-12. In both systems, we gain tremendous speedup with respect to standard MD both in exploring the configuration space and reconstructing the equilibrium distribution. In particular, we obtain 3 orders of magnitude of speedup over standard MD in the exploration of the configurational space of alanine-12 at 300 K with DM-d-MD. The method is reaction coordinate free and minimally dependent on a priori knowledge of the system. We expect wide applications of DM-d-MD to other macromolecular systems in which equilibrium sampling is not affordable by standard MD.

Delineation of Folding Pathways of a β-Sheet Miniprotein

Several methods have been developed in the past few years for the analysis of molecular dynamics simulations of biological (macro)molecules whose complexity is difficult to capture by simple projections of the free-energy surface onto one or two geometric variables. The locally scaled diffusion map (LSDMap) method is a nonlinear dimensionality reduction technique for describing the dynamics of complex systems in terms of a few collective coordinates. Here, we compare LSDMap to two previously developed approaches for the characterization of the configurational landscape associated with the folding dynamics of a three-stranded antiparallel β-sheet peptide, termed Beta3s. The analysis is aided by an improved procedure for extracting pathways from the equilibrium transition network, which enables calculation of pathway-specific cut-based free energy profiles. We find that the results from LSDMap are consistent with analysis based on transition networks and allow a coherent interpretation of metastable states and folding pathways in terms of different time scales of transitions between minima on the free energy projections.

Polymer reversal rate calculated via locally scaled diffusion map

A recent study on the dynamics of polymer reversal inside a nanopore by Huang and Makarov [J. Chem. Phys. 128, 114903 (2008)] demonstrated that the reaction rate cannot be reproduced by projecting the dynamics onto a single empirical reaction coordinate, a result suggesting the dynamics of this system cannot be correctly described by using a single collective coordinate. To further investigate this possibility we have applied our recently developed multiscale framework, locally scaled diffusion map (LSDMap), to obtain collective reaction coordinates for this system. Using a single diffusion coordinate, we obtain a reversal rate via Kramers expression that is in good agreement with the exact rate obtained from the simulations. Our mathematically rigorous approach accounts for the local heterogeneity of molecular configuration space in constructing a diffusion map, from which collective coordinates emerge. We believe this approach can be applied in general to characterize complex macromolecular dynamics by providing an accurate definition of the collective coordinates associated with processes at different time scales.

Molecular recognition of DNA by ligands: Roughness and complexity of the free energy profile

Understanding the molecular mechanism by which probes and chemotherapeutic agents bind to nucleic acids is a fundamental issue in modern drug design. From a computational perspective, valuable insights are gained by the estimation of free energy landscapes as a function of some collective variables (CVs), which are associated with the molecular recognition event. Unfortunately the choice of CVs is highly non-trivial because of DNAs high flexibility and the presence of multiple association-dissociation events at different locations and/or sliding within the grooves. Here we have applied a modified version of Locally-Scaled Diffusion Map (LSDMap), a nonlinear dimensionality reduction technique for decoupling multiple-timescale dynamics in macromolecular systems, to a metadynamics-based free energy landscape calculated using a set of intuitive CVs. We investigated the binding of the organic drug anthramycin to a DNA 14-mer duplex. By performing an extensive set of metadynamics simulations, we observed sliding of anthramycin along the full-length DNA minor groove, as well as several detachments from multiple sites, including the one identified by X-ray crystallography. As in the case of equilibrium processes, the LSDMap analysis is able to extract the most relevant collective motions, which are associated with the slow processes within the system, i.e., ligand diffusion along the minor groove and dissociation from it. Thus, LSDMap in combination with metadynamics (and possibly every equivalent method) emerges as a powerful method to describe the energetics of ligand binding to DNA without resorting to intuitive ad hoc reaction coordinates.

Multiscale Approach to the Determination of the Photoactive Yellow Protein Signaling State Ensemble

The nature of the optical cycle of photoactive yellow protein (PYP) makes its elucidation challenging for both experiment and theory. The long transition times render conventional simulation methods ineffective, and yet the short signaling-state lifetime makes experimental data difficult to obtain and interpret. Here, through an innovative combination of computational methods, a prediction and analysis of the biological signaling state of PYP is presented. Coarse-grained modeling and locally scaled diffusion map are first used to obtain a rough birds-eye view of the free energy landscape of photo-activated PYP. Then all-atom reconstruction, followed by an enhanced sampling scheme; diffusion map-directed-molecular dynamics are used to focus in on the signaling-state region of configuration space and obtain an ensemble of signaling state structures. To the best of our knowledge, this is the first time an all-atom reconstruction from a coarse grained model has been performed in a relatively unexplored region of molecular configuration space. We compare our signaling state prediction with previous computational and more recent experimental results, and the comparison is favorable, which validates the method presented. This approach provides additional insight to understand the PYP photo cycle, and can be applied to other systems for which more direct methods are impractical.