skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Polarizable embedding for simulating redox potentials of biomolecules
Redox reactions play a key role in various biological processes, including photosynthesis and respiration. Quantitative and predictive computational characterization of redox events is therefore highly desirable for enriching our knowledge on mechanistic features of biological redox-active macromolecules. Here, we present a computational protocol exploiting polarizable embedding hybrid quantum-classical approach and resulting in accurate estimates of redox potentials of biological macromolecules. A special attention is paid to fundamental aspects of the theoretical description such as the effects of environment polarization and of the long-range electrostatic interactions on the computed energetic parameters. Environment (protein and the solvent) polarization is shown to be crucial for accurate estimates of the redox potential: hybrid quantum-classical results with and without account for environment polarization differ by 1.4 V. Long-range electrostatic interactions are shown to contribute significantly to the computed redox potential value even at the distances far beyond the protein outer surface. The approach is tested on simulating reduction potential of cryptochrome 1 protein from Arabidopsis thaliana . The theoretical estimate (0.07 V) of the midpoint reduction potential is in good agreement with available experimental data (−0.15 V).  more » « less
Award ID(s):
1800505
PAR ID:
10099479
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
21
Issue:
22
ISSN:
1463-9076
Page Range / eLocation ID:
11642 to 11650
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, Proceedings of the National Academy of Sciences)
    null (Ed.)
    Electrostatic potentials computed from three-dimensional structures of biomolecules by solving the Poisson–Boltzmann equation are widely used in molecular biophysics, structural biology, and medicinal chemistry. Despite the approximate nature of the Poisson–Boltzmann theory, validation of the computed electrostatic potentials around biological macromolecules is rare and methodologically limited. Here, we present a unique and powerful NMR method that allows for straightforward and extensive comparison with electrostatic models for biomolecules and their complexes. This method utilizes paramagnetic relaxation enhancement arising from analogous cationic and anionic cosolutes whose spatial distributions around biological macromolecules reflect electrostatic potentials. We demonstrate that this NMR method enables de novo determination of near-surface electrostatic potentials for individual protein residues without using any structural information. We applied the method to ubiquitin and the Antp homeodomain–DNA complex. The experimental data agreed well with predictions from the Poisson–Boltzmann theory. Thus, our experimental results clearly support the validity of the theory for these systems. However, our experimental study also illuminates certain weaknesses of the Poisson–Boltzmann theory. For example, we found that the theory predicts stronger dependence of near-surface electrostatic potentials on ionic strength than observed in the experiments. Our data also suggest that conformational flexibility or structural uncertainties may cause large errors in theoretical predictions of electrostatic potentials, particularly for highly charged systems. This NMR-based method permits extensive assessment of near-surface electrostatic potentials for various regions around biological macromolecules and thereby may facilitate improvement of the computational approaches for electrostatic potentials. 
    more » « less
  2. ; (, Proceedings of the National Academy of Sciences)
    With a goal of determining an absolute free energy scale for ion hydration, quasi-chemical theory and ab initio quantum mechanical simulations are employed to obtain an accurate value for the bulk hydration free energy of the Na+ion. The free energy is partitioned into three parts: 1) the inner-shell or chemical contribution that includes direct interactions of the ion with nearby waters, 2) the packing free energy that is the work to produce a cavity of size λ in water, and 3) the long-range contribution that involves all interactions outside the inner shell. The interfacial potential contribution to the free energy resides in the long-range term. By averaging cation and anion data for that contribution, cumulant terms of all odd orders in the electrostatic potential are removed. The computed total is then the bulk hydration free energy. Comparison with the experimentally derived real hydration free energy produces an effective surface potential of water in the range −0.4 to −0.5 V. The result is consistent with a variety of experiments concerning acid–base chemistry, ion distributions near hydrophobic interfaces, and electric fields near the surface of water droplets. 
    more » « less
  3. ; (, Journal of Chemical Theory and Computation)
    The hydration shells of proteins mediate interactions, such as small molecule binding, that are vital to their biological function or in some cases their dysfunction. However, even when the structure of a protein is known, the properties of its hydration environment cannot be easily predicted due to the complex interplay between protein surface heterogeneity and the collective structure of water’s hydrogen bonding network. This manuscript presents a theoretical study of the influence of surface charge heterogeneity on the polarization response of the liquid water interface. We focus our attention on classical point charge models of water, where the polarization response is limited to molecular reorientation. We introduce a new computational method for analyzing simulation data that is capable of quantifying water’s collective polarization response and determining the effective surface charge distribution of hydrated surfaces over atomistic length scales. To illustrate the utility of this method, we present the results of molecular dynamics simulations of liquid water in contact with a heterogeneous model surface and the CheY protein. 
    more » « less
  4. ; ; (, Physical Chemistry Chemical Physics)
    The redox potential is a powerful thermodynamic and kinetic tool used to predict numerous chemical and biochemical mechanisms. However, despite the improving predictive power of density functional theory (DFT), chemically accurate theoretical redox potentials are often difficult to achieve with DFT. For example, calculated redox potentials are sensitive to density functional choice and often fall short of the desired accuracy. Thus, ranges of errors for computed redox potentials between different density functionals can become quite large. The current study presents a cost-effective protocol that utilizes effective error cancellation schemes in order to accurately predict the redox potentials of a wide range of organic molecules. This computational protocol, called CBH-Redox, is an extension of the connectivity-based hierarchy (CBH) method, and produces thermochemical data with near-G4 accuracy. Herein, we test the CBH-Redox protocol against both experimental and G4 reference values and compare these results to DFT alone. Considering 46 C, O, N, F, Cl, and S atom-containing molecules, when using the CBH-Redox correction scheme, the MAEs for all eight density functionals tested are within the 0.09 V target accuracy versus both experiment and G4. Moreover, CBH-Redox achieves an impressive accuracy, with a MAE of 0.05 V or below when compared to G4 for six of the eight density functionals tested. In addition, when the CBH correction is applied, the error range across all functionals tested decreases from 0.12 V to about 0.05 V versus G4, and from 0.13 V to 0.04 V versus experiment. 
    more » « less
  5. ; (, The Journal of Chemical Physics)
    PyDFT-QMMM is a Python-based package for performing hybrid quantum mechanics/molecular mechanics (QM/MM) simulations at the density functional level of theory. The program is designed to treat short-range and long-range interactions through user-specified combinations of electrostatic and mechanical embedding procedures within periodic simulation domains, providing necessary interfaces to external quantum chemistry and molecular dynamics software. To enable direct embedding of long-range electrostatics in periodic systems, we have derived and implemented force terms for our previously described QM/MM/PME approach [Pederson and McDaniel, J. Chem. Phys. 156, 174105 (2022)]. Communication with external software packages Psi4 and OpenMM is facilitated through Python application programming interfaces (APIs). The core library contains basic utilities for running QM/MM molecular dynamics simulations, and plug-in entry-points are provided for users to implement custom energy/force calculation and integration routines, within an extensible architecture. The user interacts with PyDFT-QMMM primarily through its Python API, allowing for complex workflow development with Python scripting, for example, interfacing with PLUMED for free energy simulations. We provide benchmarks of forces and energy conservation for the QM/MM/PME and alternative QM/MM electrostatic embedding approaches. We further demonstrate a simple example use case for water solute in a water solvent system, for which radial distribution functions are computed from 100 ps QM/MM simulations; in this example, we highlight how the solvation structure is sensitive to different basis-set choices due to under- or over-polarization of the QM water molecule’s electron density. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email