1
|
Zheng X, Pei Q, Tan J, Bai S, Luo Y, Ye S. Local electric field in nanocavities dictates the vibrational relaxation dynamics of interfacial molecules. Chem Sci 2024; 15:11507-11514. [PMID: 39055024 PMCID: PMC11268483 DOI: 10.1039/d4sc02463j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2024] [Accepted: 06/16/2024] [Indexed: 07/27/2024] Open
Abstract
Plasmonic nanocavities enable the generation of strong light-matter coupling and exhibit great potential in plasmon-mediated chemical reactions (PMCRs). Although an electric field generated by nanocavities (E n) has recently been reported, its effect on the vibrational energy relaxation (VER) of the molecules in the nanocavities has not been explored. In this study, we reveal the impact of an electric field sensed by molecules (para-substituted thiophenol derivatives) in a nanocavity (E f) on VER processes by employing advanced time-resolved femtosecond sum frequency generation vibrational spectroscopy (SFG-VS) supplemented by electrochemical measurements. The magnitude of E n is almost identical (1.0 ± 0.2 V nm-1) beyond the experimental deviation while E f varies from 0.3 V nm-1 to 1.7 V nm-1 depending on the substituent. An exponential correlation between E f and the complete recovery time of the ground vibrational C[double bond, length as m-dash]C state (T 2) of the phenyl ring is observed. Substances with a smaller T 2 are strongly correlated with the reported macroscopic chemical reactivity. This finding may aid in enriching the current understanding of PMCRs and highlights the possibility of regulating vibrational energy flow into desired reaction coordinates by using a local electric field.
Collapse
Affiliation(s)
- Xiaoxuan Zheng
- Hefei National Research Center for Physical Sciences at the Microscale, Department of Chemical Physics, University of Science and Technology of China Hefei Anhui 230026 China
| | - Quanbing Pei
- Hefei National Research Center for Physical Sciences at the Microscale, Department of Chemical Physics, University of Science and Technology of China Hefei Anhui 230026 China
| | - Junjun Tan
- Hefei National Laboratory, University of Science and Technology of China Hefei Anhui 230088 China
| | - Shiyu Bai
- Hefei National Research Center for Physical Sciences at the Microscale, Department of Chemical Physics, University of Science and Technology of China Hefei Anhui 230026 China
| | - Yi Luo
- Hefei National Research Center for Physical Sciences at the Microscale, Department of Chemical Physics, University of Science and Technology of China Hefei Anhui 230026 China
- Hefei National Laboratory, University of Science and Technology of China Hefei Anhui 230088 China
| | - Shuji Ye
- Hefei National Research Center for Physical Sciences at the Microscale, Department of Chemical Physics, University of Science and Technology of China Hefei Anhui 230026 China
- Hefei National Laboratory, University of Science and Technology of China Hefei Anhui 230088 China
| |
Collapse
|
2
|
Cruz R, Ataka K, Heberle J, Kozuch J. Evaluating aliphatic CF, CF2, and CF3 groups as vibrational Stark effect reporters. J Chem Phys 2024; 160:204308. [PMID: 38814010 DOI: 10.1063/5.0198303] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Accepted: 05/09/2024] [Indexed: 05/31/2024] Open
Abstract
Given the extensive use of fluorination in molecular design, it is imperative to understand the solvation properties of fluorinated compounds and the impact of the C-F bond on electrostatic interactions. Vibrational spectroscopy can provide direct insights into these interactions by using the C-F bond stretching [v(C-F)] as an electric field probe through the vibrational Stark effect (VSE). In this work, we explore the VSE of the three basic patterns of aliphatic fluorination, i.e., mono-, di-, and trifluorination in CF, CF2, and CF3 groups, respectively, and compare their response to the well-studied aromatic v(C-F). Magnitudes (i.e., Stark tuning rates) and orientations of the difference dipole vectors of the v(C-F)-containing normal modes were determined using density functional theory and a molecular dynamics (MD)-assisted solvatochromic analysis of model compounds in solvents of varying polarity. We obtain Stark tuning rates of 0.2-0.8 cm-1/(MV/cm), with smallest and largest electric field sensitivities for CFaliphatic and CF3,aliphatic, respectively. While average electric fields of solvation were oriented along the main symmetry axis of the CFn, and thus along its static dipole, the Stark tuning rate vectors were tilted by up to 87° potentially enabling to map electrostatics in multiple dimensions. We discuss the influence of conformational heterogeneity on spectral shifts and point out the importance of multipolar and/or polarizable MD force fields to describe the electrostatics of fluorinated molecules. The implications of this work are of direct relevance for studies of fluorinated molecules as found in pharmaceuticals, fluorinated peptides, and proteins.
Collapse
Affiliation(s)
- R Cruz
- Fachbereich Physik, Freie Universität Berlin, Berlin 14195, Germany
| | - K Ataka
- Fachbereich Physik, Freie Universität Berlin, Berlin 14195, Germany
| | - J Heberle
- Fachbereich Physik, Freie Universität Berlin, Berlin 14195, Germany
- Forschungsbau SupraFAB, Freie Universität Berlin, Berlin 14195, Germany
| | - J Kozuch
- Fachbereich Physik, Freie Universität Berlin, Berlin 14195, Germany
- Forschungsbau SupraFAB, Freie Universität Berlin, Berlin 14195, Germany
| |
Collapse
|
3
|
Antila HS, Dixit S, Kav B, Madsen JJ, Miettinen MS, Ollila OHS. Evaluating Polarizable Biomembrane Simulations against Experiments. J Chem Theory Comput 2024; 20:4325-4337. [PMID: 38718349 PMCID: PMC11137822 DOI: 10.1021/acs.jctc.3c01333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2023] [Revised: 04/10/2024] [Accepted: 04/10/2024] [Indexed: 05/29/2024]
Abstract
Owing to the increase of available computational capabilities and the potential for providing a more accurate description, polarizable molecular dynamics force fields are gaining popularity in modeling biomolecular systems. It is, however, crucial to evaluate how much precision is truly gained with increasing cost and complexity of the simulation. Here, we leverage the NMRlipids open collaboration and Databank to assess the performance of available polarizable lipid models─the CHARMM-Drude and the AMOEBA-based parameters─against high-fidelity experimental data and compare them to the top-performing nonpolarizable models. While some improvement in the description of ion binding to membranes is observed in the most recent CHARMM-Drude parameters, and the conformational dynamics of AMOEBA-based parameters are excellent, the best nonpolarizable models tend to outperform their polarizable counterparts for each property we explored. The identified shortcomings range from inaccuracies in describing the conformational space of lipids to excessively slow conformational dynamics. Our results provide valuable insights for the further refinement of polarizable lipid force fields and for selecting the best simulation parameters for specific applications.
Collapse
Affiliation(s)
- Hanne S. Antila
- Department
of Theory and Bio-Systems, Max Planck Institute
of Colloids and Interfaces, Potsdam 14476, Germany
- Department
of Biomedicine, University of Bergen, Bergen 5020, Norway
- Computational
Biology Unit, Department of Informatics, University of Bergen, Bergen 5008, Norway
| | - Sneha Dixit
- Department
of Theory and Bio-Systems, Max Planck Institute
of Colloids and Interfaces, Potsdam 14476, Germany
| | - Batuhan Kav
- Institute
of Biological Information Processing: Structural Biochemistry (IBI-7), Forschungszentrum Jülich, Jïulich 52428, Germany
| | - Jesper J. Madsen
- Department
of Molecular Medicine, Morsani College of Medicine, University of South Florida, Tampa, Florida 33612, United States
- Center
for Global Health and Infectious Diseases Research, Global and Planetary
Health, College of Public Health, University
of South Florida, Tampa, Florida 33612, United States of America
| | - Markus S. Miettinen
- Department
of Theory and Bio-Systems, Max Planck Institute
of Colloids and Interfaces, Potsdam 14476, Germany
- Computational
Biology Unit, Department of Informatics, University of Bergen, Bergen 5008, Norway
- Department
of Chemistry, University of Bergen, Bergen 5007, Norway
| | - O. H. Samuli Ollila
- VTT Technical
Research Centre of Finland, Espoo 02044, Finland
- Institute
of Biotechnology, University of Helsinki, Helsinki 00014, Finland
| |
Collapse
|
4
|
Sakpal S, Chakrabarty S, Reddy KD, Deshmukh SH, Biswas R, Bagchi S, Ghosh A. Perturbation of Fermi Resonance on Hydrogen-Bonded > C═O: 2D IR Studies of Small Ester Probes. J Phys Chem B 2024. [PMID: 38686937 DOI: 10.1021/acs.jpcb.3c06698] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/02/2024]
Abstract
We utilized linear and 2D infrared spectroscopy to analyze the carbonyl stretching modes of small esters in different solvents. Particularly noteworthy were the distinct carbonyl spectral line shapes in aqueous solutions, prompting our investigation of the underlying factors responsible for these differences. Through our experimental and theoretical calculations, we identified the presence of the hydrogen-bond-induced Fermi resonance as the primary contributor to the varied line shapes of small esters in aqueous solutions. Furthermore, our findings revealed that the skeletal deformation mode plays a crucial role in the Fermi resonance for all small esters. Specifically, the first overtone band of the skeletal deformation mode intensifies when hydrogen bonds form with the carbonyl group of esters, whereas such coupling is rare in aprotic organic solvents. These spectral insights carry significant implications for the utilization of esters as infrared probes in both biological and chemical systems.
Collapse
Affiliation(s)
- Sushil Sakpal
- Physical and Materials Chemistry Division, CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| | - Suranjana Chakrabarty
- Department of Condensed Matter Physics and Materials Sciences, S. N. Bose National Centre for Basic Sciences, Kolkata 700106, India
| | - Kambham Devendra Reddy
- Department of Chemistry, Indian Institute of Technology Tirupati, Tirupati, Andhra Pradesh, India 517619
| | - Samadhan H Deshmukh
- Physical and Materials Chemistry Division, CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| | - Rajib Biswas
- Department of Chemistry, Indian Institute of Technology Tirupati, Tirupati, Andhra Pradesh, India 517619
| | - Sayan Bagchi
- Physical and Materials Chemistry Division, CSIR-National Chemical Laboratory, Dr. Homi Bhabha Road, Pune 411008, India
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
| | - Anup Ghosh
- Department of Condensed Matter Physics and Materials Sciences, S. N. Bose National Centre for Basic Sciences, Kolkata 700106, India
- Department of Chemical Science, Bose Institute, Kolkata 700091, India
| |
Collapse
|
5
|
Chaturvedi SS, Bím D, Christov CZ, Alexandrova AN. From random to rational: improving enzyme design through electric fields, second coordination sphere interactions, and conformational dynamics. Chem Sci 2023; 14:10997-11011. [PMID: 37860658 PMCID: PMC10583697 DOI: 10.1039/d3sc02982d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2023] [Accepted: 09/11/2023] [Indexed: 10/21/2023] Open
Abstract
Enzymes are versatile and efficient biological catalysts that drive numerous cellular processes, motivating the development of enzyme design approaches to tailor catalysts for diverse applications. In this perspective, we investigate the unique properties of natural, evolved, and designed enzymes, recognizing their strengths and shortcomings. We highlight the challenges and limitations of current enzyme design protocols, with a particular focus on their limited consideration of long-range electrostatic and dynamic effects. We then delve deeper into the impact of the protein environment on enzyme catalysis and explore the roles of preorganized electric fields, second coordination sphere interactions, and protein dynamics for enzyme function. Furthermore, we present several case studies illustrating successful enzyme-design efforts incorporating enzyme strategies mentioned above to achieve improved catalytic properties. Finally, we envision the future of enzyme design research, spotlighting the challenges yet to be overcome and the synergy of intrinsic electric fields, second coordination sphere interactions, and conformational dynamics to push the state-of-the-art boundaries.
Collapse
Affiliation(s)
- Shobhit S Chaturvedi
- Department of Chemistry and Biochemistry, University of California, Los Angeles California 90095 USA
| | - Daniel Bím
- Department of Chemistry and Biochemistry, University of California, Los Angeles California 90095 USA
| | - Christo Z Christov
- Department of Chemistry, Michigan Technological University Houghton Michigan 49931 USA
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles California 90095 USA
| |
Collapse
|
6
|
Lin YC, Ren P, Webb LJ. AMOEBA Force Field Predicts Accurate Hydrogen Bond Counts of Nitriles in SNase by Revealing Water-Protein Interaction in Vibrational Absorption Frequencies. J Phys Chem B 2023; 127:5609-5619. [PMID: 37339399 PMCID: PMC10851345 DOI: 10.1021/acs.jpcb.3c02060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/22/2023]
Abstract
Precisely quantifying the magnitude and direction of electric fields in proteins has long been an outstanding challenge in understanding biological functions. Nitrile vibrational Stark effect probes have been shown to be minimally disruptive to the protein structure and can be better direct reporters of local electrostatic field in the native state of a protein than other measures such as pKa shifts of titratable residues. However, interpretations of the connection between measured vibrational energy and electric field rely on the accurate molecular understanding of interactions of the nitrile group and its environment, particularly from hydrogen bonding. In this work, we compared the extent of hydrogen bonding calculated in two common force fields, the fixed charge force field Amber03 and polarizable force field AMOEBA, at 10 locations of cyanocysteine (CNC) in staphylococcal nuclease (SNase) against the experimental nitrile absorption frequency in terms of full width at half-maximum (FWHM) and frequency temperature line slope (FTLS). We observed that the number of hydrogen bonds correlated well in AMOEBA trajectories with respect to both the FWHM (r = 0.88) and the FTLS (r = -0.85), whereas the correlation of Amber03 trajectories was less reliable because the Amber03 force field predicted more hydrogen bonds in some mutants. Moreover, we demonstrated that contributions from the interactions between CNC and nearby water molecules were significant in AMOEBA trajectories but were not predicted by Amber03. We conclude that although the nitrile absorption peak shape could be qualitatively predicted by the fixed charge Amber03 force field, the detailed electrostatic environment measured by the nitrile probe in terms of the extent of hydrogen bonding could only be accurately observed in the AMOEBA trajectories, where the permanent dipole, quadrupole, and dipole-induced-dipole polarizable interactions were all taken into account. The significance of this finding to the goal of accurately predicting electric fields in complex biomolecular environments is discussed.
Collapse
Affiliation(s)
- Yu-Chun Lin
- Department of Chemistry, Texas Materials Institute, and Interdisciplinary Life Sciences Program, The University of Texas at Austin, 105 E 24th St. STOP A5300, Austin, TX, 78712, USA
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, TX, 78712, USA
| | - Lauren J. Webb
- Department of Chemistry, Texas Materials Institute, and Interdisciplinary Life Sciences Program, The University of Texas at Austin, 105 E 24th St. STOP A5300, Austin, TX, 78712, USA
| |
Collapse
|
7
|
Datar A, Wright C, Matthews DA. Theoretical Investigation of the X-ray Stark Effect in Small Molecules. J Phys Chem A 2023; 127:1576-1587. [PMID: 36787229 DOI: 10.1021/acs.jpca.2c08311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023]
Abstract
We have studied the Stark effect in the soft x-ray region for various small molecules by calculating the field-dependent x-ray absorption spectra. This effect is explained in terms of the response of molecular orbitals (core and valence), the molecular dipole moment, and the molecular geometry to the applied electric field. A number of consistent trends are observed linking the computed shifts in absorption energies and intensities with specific features of the molecular electronic structure. We find that both the virtual molecular orbitals (valence and/or Rydberg) as well as the core orbitals contribute to observed trends in a complementary fashion. This initial study highlights the potential impact of x-ray Stark spectroscopy as a tool to study electronic structure and environmental perturbations at a submolecular scale.
Collapse
Affiliation(s)
- Avdhoot Datar
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Catherine Wright
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Devin A Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| |
Collapse
|
8
|
Yankova R, Tankov I, Tsaneva T. Crystal structure, intermolecular interactions and NLO properties for imidazolium hydrogen sulfate ionic liquid. J Mol Struct 2023. [DOI: 10.1016/j.molstruc.2022.134307] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|
9
|
Fica-Contreras SM, Charnay AP, Pan J, Fayer MD. Rethinking Vibrational Stark Spectroscopy: Peak Shifts, Line Widths, and the Role of Non-Stark Solvent Coupling. J Phys Chem B 2023; 127:717-731. [PMID: 36629314 DOI: 10.1021/acs.jpcb.2c06071] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
A vibration's transition frequency is partly determined by the first-order Stark effect, which accounts for the electric field experienced by the mode. Using ultrafast infrared pump-probe and FT-IR spectroscopies, we characterized both the 0 → 1 and 1 → 2 vibrational transitions' field-dependent peak positions and line widths of the CN stretching mode of benzonitrile (BZN) and phenyl selenocyanate (PhSeCN) in ten solvents. We present a theoretical model that decomposes the observed line width into a field-dependent Stark contribution and a field-independent non-Stark solvent coupling contribution (NSC). The model demonstrates that the field-dependent peak position is independent of the line width, even when the NSC dominates the latter. Experiments show that when the Stark tuning rate is large compared to the NSC (PhSeCN), the line width has a field dependence, albeit with major NSC-induced excursions from linearity. When the Stark tuning rate is small relative to the NSC (BZN), the line width is field-independent. BZN's line widths are substantially larger for the 1 → 2 transition, indicating a 1 → 2 transition enhancement of the NSC. Additionally, we examine, theoretically and experimentally, the difference in the 0 → 1 and 1 → 2 transitions' Stark tuning rates. Second-order perturbation theory combined with density functional theory explain the difference and show that the 1 → 2 transition's Stark tuning rate is ∼10% larger. The Stark tuning rate of PhSeCN is larger than BZN's for both transitions, consistent with the theoretical calculations. This study provides new insights into vibrational line shape components and a more general understanding of the vibrational response to external electric fields.
Collapse
Affiliation(s)
| | - Aaron P Charnay
- Department of Chemistry, Stanford University, Stanford, California94305, United States
| | - Junkun Pan
- Department of Chemistry, Stanford University, Stanford, California94305, United States
| | - Michael D Fayer
- Department of Chemistry, Stanford University, Stanford, California94305, United States
| |
Collapse
|
10
|
Martins-Costa MTC, Ruiz-López MF. Electrostatics and Chemical Reactivity at the Air-Water Interface. J Am Chem Soc 2023; 145:1400-1406. [PMID: 36622259 DOI: 10.1021/jacs.2c12089] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
It has been recently discovered that chemical reactions at aqueous interfaces can be orders of magnitude faster compared to conventional bulk phase reactions, but despite its wide-ranging implications, which extend from atmospheric to synthetic chemistry or technological applications, the phenomenon is still incompletely understood. The role of strong electric fields due to space asymmetry and the accumulation of ions at the interface has been claimed as a possible cause from some experiments, but the reorganization of the solvent around the reactive system should provide even greater additional electrostatic contributions that have not yet been analyzed. In this study, with the help of first-principles molecular dynamics simulations, we go deeper into this issue by a careful assessment of solvation electrostatics at the air-water interface. Our simulations confirm that electrostatic forces can indeed be a key factor in rate acceleration compared to bulk solution. Remarkably, the study reveals that the effect cannot simply be attributed to the magnitude of the local electric field and that the fluctuations of the full electrostatic potential resulting from unique dynamical behavior of the solvation shells at the interface must be accounted for. This finding paves the way for future applications of the phenomenon in organic synthesis, especially for charge transfer or redox reactions in thin films and microdroplets.
Collapse
Affiliation(s)
- Marilia T C Martins-Costa
- Laboratoire de Physique et Chimie Théoriques, UMR CNRS 7019, University of Lorraine, CNRS, BP 70239, 54506 Vandoeuvre-lès-Nancy, France
| | - Manuel F Ruiz-López
- Laboratoire de Physique et Chimie Théoriques, UMR CNRS 7019, University of Lorraine, CNRS, BP 70239, 54506 Vandoeuvre-lès-Nancy, France
| |
Collapse
|
11
|
Song C. State averaged CASSCF in AMOEBA polarizable water model for simulating nonadiabatic molecular dynamics with nonequilibrium solvation effects. J Chem Phys 2023; 158:014101. [PMID: 36610973 DOI: 10.1063/5.0131689] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
This paper presents a state-averaged complete active space self-consistent field (SA-CASSCF) in the atomic multipole optimized energetics for biomolecular application (AMOEBA) polarizable water model, which enables rigorous simulation of non-adiabatic molecular dynamics with nonequilibrium solvation effects. The molecular orbital and configuration interaction coefficients of the solute wavefunction, and the induced dipoles on solvent atoms, are solved by minimizing the state averaged energy variationally. In particular, by formulating AMOEBA water models and the polarizable continuum model (PCM) in a unified way, the algorithms developed for computing SA-CASSCF/PCM energies, analytical gradients, and non-adiabatic couplings in our previous work can be generalized to SA-CASSCF/AMOEBA by properly substituting a specific list of variables. Implementation of this method will be discussed with the emphasis on how the calculations of different terms are partitioned between the quantum chemistry and molecular mechanics codes. We will present and discuss results that demonstrate the accuracy and performance of the implementation. Next, we will discuss results that compare three solvent models that work with SA-CASSCF, i.e., PCM, fixed-charge force fields, and the newly implemented AMOEBA. Finally, the new SA-CASSCF/AMOEBA method has been interfaced with the ab initio multiple spawning method to carry out non-adiabatic molecular dynamics simulations. This method is demonstrated by simulating the photodynamics of the model retinal protonated Schiff base molecule in water.
Collapse
Affiliation(s)
- Chenchen Song
- Department of Chemistry, University of California Davis, Davis, California 95616, USA
| |
Collapse
|
12
|
Lin YC, Ren P, Webb LJ. AMOEBA Force Field Trajectories Improve Predictions of Accurate p Ka Values of the GFP Fluorophore: The Importance of Polarizability and Water Interactions. J Phys Chem B 2022; 126:7806-7817. [PMID: 36194474 PMCID: PMC10851343 DOI: 10.1021/acs.jpcb.2c03642] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Precisely quantifying the magnitude, direction, and biological functions of electric fields in proteins has long been an outstanding challenge in the field. The most widely implemented experimental method to measure such electric fields at a particular residue in a protein has been through changes in pKa of titratable residues. While many computational strategies exist to predict these values, it has been difficult to do this accurately or connect predicted results to key structural or mechanistic features of the molecule. Here, we used experimentally determined pKa values of the fluorophore in superfolder green fluorescent protein (GFP) with amino acid mutations made at position Thr 203 to evaluate the pKa prediction ability of molecular dynamics (MD) simulations using a polarizable force field, AMOEBA. Structure ensembles from AMOEBA were used to calculate pKa values of the GFP fluorophore. The calculated pKa values were then compared to trajectories using a conventional fixed charge force field (Amber03 ff). We found that the position of water molecules included in the pKa calculation had opposite effects on the pKa values between the trajectories from AMOEBA and Amber03 force fields. In AMOEBA trajectories, the inclusion of water molecules within 35 Å of the fluorophore decreased the difference between the predicted and experimental values, resulting in calculated pKa values that were within an average of 0.8 pKa unit from the experimental results. On the other hand, in Amber03 trajectories, including water molecules that were more than 5 Å from the fluorophore increased the differences between the calculated and experimental pKa values. The inaccuracy of pKa predictions determined from Amber03 trajectories was caused by a significant stabilization of the deprotonated chromophore's free energy compared to the result in AMOEBA. We rationalize the cutoffs for explicit water molecules when calculating pKa to better predict the electrostatic environment surrounding the fluorophore buried in GFP. We discuss how the results from this work will assist the prospective prediction of pKa values or other electrostatic effects in a wide variety of folded proteins.
Collapse
Affiliation(s)
- Yu-Chun Lin
- Department of Chemistry, Texas Materials Institute, and Interdisciplinary Life Sciences Program, The University of Texas at Austin, 105 E 24th St. STOP A5300, Austin, TX 78712-1224
| | - Pengyu Ren
- Department of Chemistry, Texas Materials Institute, and Interdisciplinary Life Sciences Program, The University of Texas at Austin, 105 E 24th St. STOP A5300, Austin, TX 78712-1224
| | - Lauren J. Webb
- Department of Chemistry, Texas Materials Institute, and Interdisciplinary Life Sciences Program, The University of Texas at Austin, 105 E 24th St. STOP A5300, Austin, TX 78712-1224
| |
Collapse
|
13
|
Fried SDE, Zheng C, Mao Y, Markland TE, Boxer SG. Solvent Organization and Electrostatics Tuned by Solute Electronic Structure: Amide versus Non-Amide Carbonyls. J Phys Chem B 2022; 126:5876-5886. [PMID: 35901512 PMCID: PMC10081530 DOI: 10.1021/acs.jpcb.2c03095] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The ability to exploit carbonyl groups to measure electric fields in enzymes and other complex reactive environments by using the vibrational Stark effect has inspired growing interest in how these fields can be measured, tuned, and ultimately designed. Previous studies have concentrated on the role of the solvent in tuning the fields exerted on the solute. Here, we explore instead the role of the solute electronic structure in modifying the local solvent organization and electric field exerted on the solute. By measuring the infrared absorption spectra of amide-containing molecules, as prototypical peptides, and contrasting them with non-amide carbonyls in a wide range of solvents, we show that these solutes experience notable differences in their frequency shifts in polar solvents. Using vibrational Stark spectroscopy and molecular dynamics simulations, we demonstrate that while some of these differences can be rationalized by using the distinct intrinsic Stark tuning rates of the solutes, the larger frequency shifts for amides and dimethylurea primarily result from the larger solvent electric fields experienced by their carbonyl groups. These larger fields arise due to their stronger p-π conjugation, which results in larger C═O bond dipole moments that further induce substantial solvent organization. Using electronic structure calculations, we decompose the electric fields into contributions from solvent molecules that are in the first solvation shell and those from the bulk and show that both of these contributions are significant and become larger with enhanced conjugation in solutes. These results show that structural modifications of a solute can be used to tune both the solvent organization and electrostatic environment, indicating the importance of a solute-centric paradigm in modulating and designing the electrostatic environment in condensed-phase chemical processes.
Collapse
Affiliation(s)
- Steven D E Fried
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Chu Zheng
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Steven G Boxer
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
14
|
Zheng C, Mao Y, Kozuch J, Atsango AO, Ji Z, Markland TE, Boxer SG. A two-directional vibrational probe reveals different electric field orientations in solution and an enzyme active site. Nat Chem 2022; 14:891-897. [PMID: 35513508 PMCID: PMC10082611 DOI: 10.1038/s41557-022-00937-w] [Citation(s) in RCA: 34] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Accepted: 03/25/2022] [Indexed: 12/26/2022]
Abstract
The catalytic power of an electric field depends on its magnitude and orientation with respect to the reactive chemical species. Understanding and designing new catalysts for electrostatic catalysis thus requires methods to measure the electric field orientation and magnitude at the molecular scale. We demonstrate that electric field orientations can be extracted using a two-directional vibrational probe by exploiting the vibrational Stark effect of both the C=O and C-D stretches of a deuterated aldehyde. Combining spectroscopy with molecular dynamics and electronic structure partitioning methods, we demonstrate that, despite distinct polarities, solvents act similarly in their preference for electrostatically stabilizing large bond dipoles at the expense of destabilizing small ones. In contrast, we find that for an active-site aldehyde inhibitor of liver alcohol dehydrogenase, the electric field orientation deviates markedly from that found in solvents, which provides direct evidence for the fundamental difference between the electrostatic environment of solvents and that of a preorganized enzyme active site.
Collapse
Affiliation(s)
- Chu Zheng
- Department of Chemistry, Stanford University, Stanford, CA, USA
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, CA, USA
| | - Jacek Kozuch
- Experimental Molecular Biophysics, Department of Physics, Freie Univeresität Berlin, Berlin, Germany
| | | | - Zhe Ji
- Department of Chemistry, Stanford University, Stanford, CA, USA
| | | | - Steven G Boxer
- Department of Chemistry, Stanford University, Stanford, CA, USA.
| |
Collapse
|
15
|
Hu K, Matsuura H, Shirakashi R. Stochastic Analysis of Molecular Dynamics Reveals the Rotation Dynamics Distribution of Water around Lysozyme. J Phys Chem B 2022; 126:4520-4530. [PMID: 35675630 DOI: 10.1021/acs.jpcb.2c00970] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Water dynamics is essential to biochemical processes by mediating all such reactions, including biomolecular degeneration in solutions. To disentangle the molecular-scale distribution of water dynamics around a solute biomolecule, we investigated here the rotational dynamics of water around lysozyme by combining molecular dynamics (MD) simulations and broadband dielectric spectroscopy (BDS). A statistical analysis using the relaxation times and trajectories of every single water molecule was proposed, and the two-dimensional probability distribution of water at a distance from the lysozyme surface with a rotational relaxation time was given. For the observed lysozyme solutions of 34-284 mg/mL, we discovered that the dielectric relaxation time obtained from this distribution agrees well with the measured γ relaxation time, which suggests that rotational self-correlation of water molecules underlies the gigahertz domain of the dielectric spectra. Regardless of protein concentration, water rotational relaxation time versus the distance from the lysozyme surface revealed that the water rotation is severely retarded within 3 Å from the lysozyme surface and is nearly comparable to pure water when farther than 10 Å. The dimension of the first hydration layer was subsequently identified in terms of the relationship between the acceleration of water rotation and the distance from the protein surface.
Collapse
Affiliation(s)
- Kang Hu
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro City, Tokyo 153-8505, Japan.,Department of Mechanical Engineering, The University of Tokyo, 7-3-1, Bunkyo-ku, Tokyo 113-0033, Japan
| | - Hiroaki Matsuura
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro City, Tokyo 153-8505, Japan.,Research Fellow of the Japan Society for the Promotion of Science, Chiyoda-ku, Tokyo 102-0083, Japan
| | - Ryo Shirakashi
- Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro City, Tokyo 153-8505, Japan
| |
Collapse
|
16
|
Tankov I, Yankova R. Crystal structure, chemical bonds nature and thermodynamic functions for the new ionic liquid pyridinium bis(dihydrogen phosphate). J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2021.131363] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
17
|
Eggert A, Etling C, Millbrodt L, Schulz G, Kalesse M. Oxazaborolidinone-Mediated Asymmetric Bisvinylogous Mukaiyama Aldol Reaction. Org Lett 2021; 23:8722-8726. [PMID: 34738464 DOI: 10.1021/acs.orglett.1c03165] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A bisvinylogous Mukaiyama aldol reaction using oxazaborolidinones as a source of chirality was developed. This methodology allows the fast assembly of conjugated dienols by expanding the vinylogy principle by two additional carbons, and can be conducted using a readily available Lewis acid at reasonable reaction times. A broad range of aromatic and aliphatic aldehydes can be used providing access to complex building blocks for polyketide synthesis.
Collapse
Affiliation(s)
- Alina Eggert
- Institute for Organic Chemistry, Leibniz University Hannover, Schneiderberg 1b, 30167 Hannover, Germany.,Centre of Biomolecular Drug Research (BMWZ), Schneiderberg 38, 30167 Hannover, Germany
| | - Christoph Etling
- Institute for Organic Chemistry, Leibniz University Hannover, Schneiderberg 1b, 30167 Hannover, Germany.,Centre of Biomolecular Drug Research (BMWZ), Schneiderberg 38, 30167 Hannover, Germany
| | - Lucas Millbrodt
- Institute for Organic Chemistry, Leibniz University Hannover, Schneiderberg 1b, 30167 Hannover, Germany.,Centre of Biomolecular Drug Research (BMWZ), Schneiderberg 38, 30167 Hannover, Germany
| | - Göran Schulz
- Institute for Organic Chemistry, Leibniz University Hannover, Schneiderberg 1b, 30167 Hannover, Germany
| | - Markus Kalesse
- Institute for Organic Chemistry, Leibniz University Hannover, Schneiderberg 1b, 30167 Hannover, Germany.,Centre of Biomolecular Drug Research (BMWZ), Schneiderberg 38, 30167 Hannover, Germany.,Helmholtz Centre for Infection Research GmbH (HZI), Inhoffenstraße 7, 38124 Braunschweig, Germany
| |
Collapse
|
18
|
Kathmann SM. Electric fields and potentials in condensed phases. Phys Chem Chem Phys 2021; 23:23836-23849. [PMID: 34647950 DOI: 10.1039/d1cp03571a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The electric fields and potentials inside and at the interface of matter are relevant to many branches of physics, chemistry, and biology. Accurate quantification of these fields and/or potentials is essential to control and exploit chemical and physical transformations. Before we understand the response of matter to external fields, it is first important to understand the intrinsic interior and interfacial fields and potentials, both classically and quantum mechanically, as well as how they are probed experimentally. Here we compare and contrast, beginning with the hydrogen atom in vacuum and ending with concentrated aqueous NaCl electrolyte, both classical and quantum mechanical electric potentials and fields. We make contact with experimental vibrational Stark, electrochemical, X-ray, and electron spectroscopic probes of these potentials and fields, outline relevant conceptual difficulties, and underscore the advantage of electron holography as a basis to better understand electrostatics in matter.
Collapse
Affiliation(s)
- Shawn M Kathmann
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, WA, 99354, USA.
| |
Collapse
|
19
|
Drexler CI, Cracchiolo OM, Myers RL, Okur HI, Serrano AL, Corcelli SA, Cremer PS. Local Electric Fields in Aqueous Electrolytes. J Phys Chem B 2021; 125:8484-8493. [PMID: 34313130 DOI: 10.1021/acs.jpcb.1c03257] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Vibrational Stark shifts were explored in aqueous solutions of organic molecules with carbonyl- and nitrile-containing constituents. In many cases, the vibrational resonances from these moieties shifted toward lower frequency as salt was introduced into solution. This is in contrast to the blue-shift that would be expected based upon Onsager's reaction field theory. Salts containing well-hydrated cations like Mg2+ or Li+ led to the most pronounced Stark shift for the carbonyl group, while poorly hydrated cations like Cs+ had the greatest impact on nitriles. Moreover, salts containing I- gave rise to larger Stark shifts than those containing Cl-. Molecular dynamics simulations indicated that cations and anions both accumulate around the probe in an ion- and probe-dependent manner. An electric field was generated by the ion pair, which pointed from the cation to the anion through the vibrational chromophore. This resulted from solvent-shared binding of the ions to the probes, consistent with their positions in the Hofmeister series. The "anti-Onsager" Stark shifts occur in both vibrational spectroscopy and fluorescence measurements.
Collapse
Affiliation(s)
| | - Olivia M Cracchiolo
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | | | - Halil I Okur
- Department of Chemistry and National Nanotechnology Research Center (UNAM), Bilkent University, 06800 Ankara, Turkey
| | - Arnaldo L Serrano
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Steven A Corcelli
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | | |
Collapse
|
20
|
Bhattacharyya D, Videla PE, Cattaneo M, Batista VS, Lian T, Kubiak CP. Vibrational Stark shift spectroscopy of catalysts under the influence of electric fields at electrode-solution interfaces. Chem Sci 2021; 12:10131-10149. [PMID: 34377403 PMCID: PMC8336477 DOI: 10.1039/d1sc01876k] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2021] [Accepted: 07/13/2021] [Indexed: 11/24/2022] Open
Abstract
External control of chemical processes is a subject of widespread interest in chemical research, including control of electrocatalytic processes with significant promise in energy research. The electrochemical double-layer is the nanoscale region next to the electrode/electrolyte interface where chemical reactions typically occur. Understanding the effects of electric fields within the electrochemical double layer requires a combination of synthesis, electrochemistry, spectroscopy, and theory. In particular, vibrational sum frequency generation (VSFG) spectroscopy is a powerful technique to probe the response of molecular catalysts at the electrode interface under bias. Fundamental understanding can be obtained via synthetic tuning of the adsorbed molecular catalysts on the electrode surface and by combining experimental VSFG data with theoretical modelling of the Stark shift response. The resulting insights at the molecular level are particularly valuable for the development of new methodologies to control and characterize catalysts confined to electrode surfaces. This Perspective article is focused on how systematic modifications of molecules anchored to surfaces report information concerning the geometric, energetic, and electronic parameters of catalysts under bias attached to electrode surfaces. Heterogeneous electrocatalysis: characterization of interfacial electric field within the electrochemical double layer.![]()
Collapse
Affiliation(s)
- Dhritiman Bhattacharyya
- Department of Chemistry, Emory University 1515 Dickey Drive Northeast Atlanta Georgia 30322 USA
| | - Pablo E Videla
- Department of Chemistry and Energy Sciences Institute, Yale University 225 Prospect Street New Haven Connecticut 06520 USA
| | - Mauricio Cattaneo
- INQUINOA-UNT-CONICET, Facultad de Bioquímica, Química y Farmacia, Instituto de Química Física, Universidad Nacional de Tucumán Ayacucho 471 (4000) San Miguel de Tucumán Argentina
| | - Victor S Batista
- Department of Chemistry and Energy Sciences Institute, Yale University 225 Prospect Street New Haven Connecticut 06520 USA
| | - Tianquan Lian
- Department of Chemistry, Emory University 1515 Dickey Drive Northeast Atlanta Georgia 30322 USA
| | - Clifford P Kubiak
- Department of Chemistry and Biochemistry, University of California, San Diego 9500 Gilman Drive, MC 0358 La Jolla California 92093 USA
| |
Collapse
|
21
|
Kozuch J, Schneider SH, Zheng C, Ji Z, Bradshaw RT, Boxer SG. Testing the Limitations of MD-Based Local Electric Fields Using the Vibrational Stark Effect in Solution: Penicillin G as a Test Case. J Phys Chem B 2021; 125:4415-4427. [PMID: 33900769 PMCID: PMC8522303 DOI: 10.1021/acs.jpcb.1c00578] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Noncovalent interactions underlie nearly all molecular processes in the condensed phase from solvation to catalysis. Their quantification within a physically consistent framework remains challenging. Experimental vibrational Stark effect (VSE)-based solvatochromism can be combined with molecular dynamics (MD) simulations to quantify the electrostatic forces in solute-solvent interactions for small rigid molecules and, by extension, when these solutes bind in enzyme active sites. While generalizing this approach toward more complex (bio)molecules, such as the conformationally flexible and charged penicillin G (PenG), we were surprised to observe inconsistencies in MD-based electric fields. Combining synthesis, VSE spectroscopy, and computational methods, we provide an intimate view on the origins of these discrepancies. We observe that the electric fields are correlated to conformation-dependent effects of the flexible PenG side chain, including both the local solvation structure and solute conformational sampling in MD. Additionally, we identified that MD-based electric fields are consistently overestimated in three-point water models in the vicinity of charged groups; this cannot be entirely ameliorated using polarizable force fields (AMOEBA) or advanced water models. This work demonstrates the value of the VSE as a direct method for experiment-guided refinements of MD force fields and establishes a general reductionist approach to calibrating vibrational probes for complex (bio)molecules.
Collapse
Affiliation(s)
- Jacek Kozuch
- Department of Chemistry, Stanford University, Stanford, California 94305-5012, United States
| | - Samuel H Schneider
- Department of Chemistry, Stanford University, Stanford, California 94305-5012, United States
| | - Chu Zheng
- Department of Chemistry, Stanford University, Stanford, California 94305-5012, United States
| | - Zhe Ji
- Department of Chemistry, Stanford University, Stanford, California 94305-5012, United States
| | - Richard T Bradshaw
- Department of Chemistry, King's College London, Britannia House, 7 Trinity Street, London SE1 1DB, U.K
| | - Steven G Boxer
- Department of Chemistry, Stanford University, Stanford, California 94305-5012, United States
| |
Collapse
|
22
|
Yu S, Vermeeren P, Hamlin TA, Bickelhaupt FM. How Oriented External Electric Fields Modulate Reactivity. Chemistry 2021; 27:5683-5693. [PMID: 33289179 PMCID: PMC8049047 DOI: 10.1002/chem.202004906] [Citation(s) in RCA: 32] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Revised: 12/04/2020] [Indexed: 01/27/2023]
Abstract
A judiciously oriented external electric field (OEEF) can catalyze a wide range of reactions and can even induce endo/exo stereoselectivity of cycloaddition reactions. The Diels-Alder reaction between cyclopentadiene and maleic anhydride is studied by using quantitative activation strain and Kohn-Sham molecular orbital theory to pinpoint the origin of these catalytic and stereoselective effects. Our quantitative model reveals that an OEEF along the reaction axis induces an enhanced electrostatic and orbital interaction between the reactants, which in turn lowers the reaction barrier. The stronger electrostatic interaction originates from an increased electron density difference between the reactants at the reactive center, and the enhanced orbital interaction arises from the promoted normal electron demand donor-acceptor interaction driven by the OEEF. An OEEF perpendicular to the plane of the reaction axis solely stabilizes the exo pathway of this reaction, whereas the endo pathway remains unaltered and efficiently steers the endo/exo stereoselectivity. The influence of the OEEF on the inverse electron demand Diels-Alder reaction is also investigated; unexpectedly, it inhibits the reaction, as the electric field now suppresses the critical inverse electron demand donor-acceptor interaction.
Collapse
Affiliation(s)
- Song Yu
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
| | - Pascal Vermeeren
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
| | - Trevor A. Hamlin
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
| | - F. Matthias Bickelhaupt
- Department of Theoretical ChemistryAmsterdam Institute of Molecular and Life Sciences (AIMMS)Amsterdam Center for Multiscale Modeling (ACMM)Vrije Universiteit AmsterdamDe Boelelaan 10831081HVAmsterdamThe Netherlands
- Institute for Molecules and Materials (IMM)Radboud UniversityHeyendaalseweg 1356525 AJNijmegenThe Netherlands
| |
Collapse
|
23
|
Drobizhev M, Molina RS, Callis PR, Scott JN, Lambert GG, Salih A, Shaner NC, Hughes TE. Local Electric Field Controls Fluorescence Quantum Yield of Red and Far-Red Fluorescent Proteins. Front Mol Biosci 2021; 8:633217. [PMID: 33763453 PMCID: PMC7983054 DOI: 10.3389/fmolb.2021.633217] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 01/06/2021] [Indexed: 12/17/2022] Open
Abstract
Genetically encoded probes with red-shifted absorption and fluorescence are highly desirable for imaging applications because they can report from deeper tissue layers with lower background and because they provide additional colors for multicolor imaging. Unfortunately, red and especially far-red fluorescent proteins have very low quantum yields, which undermines their other advantages. Elucidating the mechanism of nonradiative relaxation in red fluorescent proteins (RFPs) could help developing ones with higher quantum yields. Here we consider two possible mechanisms of fast nonradiative relaxation of electronic excitation in RFPs. The first, known as the energy gap law, predicts a steep exponential drop of fluorescence quantum yield with a systematic red shift of fluorescence frequency. In this case the relaxation of excitation occurs in the chromophore without any significant changes of its geometry. The second mechanism is related to a twisted intramolecular charge transfer in the excited state, followed by an ultrafast internal conversion. The chromophore twisting can strongly depend on the local electric field because the field can affect the activation energy. We present a spectroscopic method of evaluating local electric fields experienced by the chromophore in the protein environment. The method is based on linear and two-photon absorption spectroscopy, as well as on quantum-mechanically calculated parameters of the isolated chromophore. Using this method, which is substantiated by our molecular dynamics simulations, we obtain the components of electric field in the chromophore plane for seven different RFPs with the same chromophore structure. We find that in five of these RFPs, the nonradiative relaxation rate increases with the strength of the field along the chromophore axis directed from the center of imidazolinone ring to the center of phenolate ring. Furthermore, this rate depends on the corresponding electrostatic energy change (calculated from the known fields and charge displacements), in quantitative agreement with the Marcus theory of charge transfer. This result supports the dominant role of the twisted intramolecular charge transfer mechanism over the energy gap law for most of the studied RFPs. It provides important guidelines of how to shift the absorption wavelength of an RFP to the red, while keeping its brightness reasonably high.
Collapse
Affiliation(s)
- Mikhail Drobizhev
- Department of Cell Biology and Neuroscience, Montana State University, Bozeman, MT, United States
| | - Rosana S Molina
- Department of Cell Biology and Neuroscience, Montana State University, Bozeman, MT, United States
| | - Patrik R Callis
- Department of Chemistry and Biochemistry, Montana State University, Bozeman, MT, United States
| | | | - Gerard G Lambert
- Department of Neurosciences, UC San Diego, San Diego, CA, United States
| | - Anya Salih
- Antares & Fluoresci Research, Dangar Island, NSW, Australia
| | - Nathan C Shaner
- Department of Neurosciences, UC San Diego, San Diego, CA, United States
| | - Thomas E Hughes
- Department of Cell Biology and Neuroscience, Montana State University, Bozeman, MT, United States
| |
Collapse
|
24
|
Raiteri P, Kraus P, Gale JD. Molecular dynamics simulations of liquid-liquid interfaces in an electric field: The water-1,2-dichloroethane interface. J Chem Phys 2020; 153:164714. [PMID: 33138425 DOI: 10.1063/5.0027876] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023] Open
Abstract
The polarized interface between two immiscible liquids plays a central role in many technological processes. In particular, for electroanalytical and ion extraction applications, an external electric field is typically used to selectively induce the transfer of ionic species across the interfaces. Given that it is experimentally challenging to obtain an atomistic insight into the ion transfer process and the structure of liquid-liquid interfaces, atomistic simulations have often been used to fill this knowledge gap. However, due to the long-range nature of the electrostatic interactions and the use of 3D periodic boundary conditions, the use of external electric fields in molecular dynamics simulations requires special care. Here, we show how the simulation setup affects the dielectric response of the materials and demonstrate how by a careful design of the system it is possible to obtain the correct electric field on both sides of a liquid-liquid interface when using standard 3D Ewald summation methods. In order to prove the robustness of our approach, we ran extensive molecular dynamics simulations with a rigid-ion and polarizable force field of the water/1,2-dichloroethane interface in the presence of weak external electric fields.
Collapse
Affiliation(s)
- Paolo Raiteri
- Curtin Institute for Computation, School of Molecular and Life Sciences, Curtin University, P.O. Box U1987, Perth, WA 6845, Australia
| | - Peter Kraus
- Curtin Institute for Computation, School of Molecular and Life Sciences, Curtin University, P.O. Box U1987, Perth, WA 6845, Australia
| | - Julian D Gale
- Curtin Institute for Computation, School of Molecular and Life Sciences, Curtin University, P.O. Box U1987, Perth, WA 6845, Australia
| |
Collapse
|
25
|
Ahsan M, Pindi C, Senapati S. Electrostatics Plays a Crucial Role in HIV-1 Protease Substrate Binding, Drugs Fail to Take Advantage. Biochemistry 2020; 59:3316-3331. [PMID: 32822154 DOI: 10.1021/acs.biochem.0c00341] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
HIV-1 protease (HIVPR) is an important drug target for combating AIDS. This enzyme is an aspartyl protease that is functionally active in its dimeric form. Nuclear magnetic resonance reports have convincingly shown that a pseudosymmetry exists at the HIVPR active site, where only one of the two aspartates remains protonated over the pH range of 2.5-7.0. To date, all HIVPR-targeted drug design strategies focused on maximizing the size-shape complementarity and van der Waals interactions of the small molecule drugs with the deprotonated, symmetric active site envelope of crystallized HIVPR. However, these strategies were ineffective with the emergence of drug resistant protease variants, primarily due to the steric clashes at the active site. In this study, we traced a specificity in the substrate binding motif that emerges primarily from the asymmetrical electrostatic potential present in the protease active site due to the uneven protonation. Our detailed results from atomistic molecular dynamics simulations show that while such a specific mode of substrate binding involves significant electrostatic interactions, none of the existing drugs or inhibitors could utilize this electrostatic hot spot. As the electrostatic is long-range interaction, it can provide sufficient binding strength without the necessity of increasing the bulkiness of the inhibitors. We propose that introducing the electrostatic component along with optimal fitting at the binding pocket could pave the way for promising designs that might be more effective against both wild type and HIVPR resistant variants.
Collapse
Affiliation(s)
- Mohd Ahsan
- Department of Biotechnology and BJM School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Chinmai Pindi
- Department of Biotechnology and BJM School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| | - Sanjib Senapati
- Department of Biotechnology and BJM School of Biosciences, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
26
|
Kempfer-Robertson EM, Thompson LM. Effect of Oriented External Electric Fields on the Photo and Thermal Isomerization of Azobenzene. J Phys Chem A 2020; 124:3520-3529. [PMID: 32286821 DOI: 10.1021/acs.jpca.0c00492] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Azobenzene is a prototype molecule with potential applications in molecular switches, solar thermal batteries, sensors, photoresponsive membranes, molecular electronics, data storage, and nonlinear optics. Photo and thermal isomerization pathways exhibit different charge-transfer character and dipole moments, implying that the use of electric fields can be used to modulate the reactivity of azobenzene. This article examines the differential effect of orientated electric fields on the rotation and inversion thermal and photoisomerization pathways of azobenzene to explore the feasibility of using electric fields in the design of azobenzene-based molecular devices. Our findings demonstrate that the application of orientated electric fields modifies the accessibility of the S0/S1 seam of electronic degeneracy, as well as changes the energetically favored relaxation pathway in the branching space to yield different photoproducts. In addition, we observed strong-field dipole-inversion effects that cause a topographical change in the response of the potential energy surface to the applied field and can result in geometric minima that do not exist under field-free conditions. On the S0 surface, transition barriers can be modified on the order of ±10 kcal mol-1, enabling control of thermal isomerization rates.
Collapse
Affiliation(s)
| | - Lee M Thompson
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40205, United States
| |
Collapse
|
27
|
Verma N, Tao Y, Zou W, Chen X, Chen X, Freindorf M, Kraka E. A Critical Evaluation of Vibrational Stark Effect (VSE) Probes with the Local Vibrational Mode Theory. SENSORS 2020; 20:s20082358. [PMID: 32326248 PMCID: PMC7219233 DOI: 10.3390/s20082358] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/20/2020] [Revised: 04/14/2020] [Accepted: 04/15/2020] [Indexed: 02/06/2023]
Abstract
Over the past two decades, the vibrational Stark effect has become an important tool to measure and analyze the in situ electric field strength in various chemical environments with infrared spectroscopy. The underlying assumption of this effect is that the normal stretching mode of a target bond such as CO or CN of a reporter molecule (termed vibrational Stark effect probe) is localized and free from mass-coupling from other internal coordinates, so that its frequency shift directly reflects the influence of the vicinal electric field. However, the validity of this essential assumption has never been assessed. Given the fact that normal modes are generally delocalized because of mass-coupling, this analysis was overdue. Therefore, we carried out a comprehensive evaluation of 68 vibrational Stark effect probes and candidates to quantify the degree to which their target normal vibration of probe bond stretching is decoupled from local vibrations driven by other internal coordinates. The unique tool we used is the local mode analysis originally introduced by Konkoli and Cremer, in particular the decomposition of normal modes into local mode contributions. Based on our results, we recommend 31 polyatomic molecules with localized target bonds as ideal vibrational Stark effect probe candidates.
Collapse
Affiliation(s)
- Niraj Verma
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX 75275-0314, USA; (N.V.); (Y.T.); (M.F.)
| | - Yunwen Tao
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX 75275-0314, USA; (N.V.); (Y.T.); (M.F.)
| | - Wenli Zou
- Institute of Modern Physics, Northwest University, Xi’an 710127, China;
| | - Xia Chen
- Hubei Key Laboratory of Natural Medicinal Chemistry and Resource Evaluation, School of Pharmacy, Tongji Medical College, Huazhong University of Science and Technology, Wuhan 430030, China;
| | - Xin Chen
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, China;
| | - Marek Freindorf
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX 75275-0314, USA; (N.V.); (Y.T.); (M.F.)
| | - Elfi Kraka
- Department of Chemistry, Southern Methodist University, 3215 Daniel Avenue, Dallas, TX 75275-0314, USA; (N.V.); (Y.T.); (M.F.)
- Correspondence:
| |
Collapse
|
28
|
Bradshaw RT, Dziedzic J, Skylaris CK, Essex JW. The Role of Electrostatics in Enzymes: Do Biomolecular Force Fields Reflect Protein Electric Fields? J Chem Inf Model 2020; 60:3131-3144. [DOI: 10.1021/acs.jcim.0c00217] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Richard T. Bradshaw
- School of Chemistry, University of Southampton, Highfield Campus, Southampton, SO17 1BJ, United Kingdom
| | - Jacek Dziedzic
- School of Chemistry, University of Southampton, Highfield Campus, Southampton, SO17 1BJ, United Kingdom
- Faculty of Applied Physics and Mathematics, Gdańsk University of Technology, 80-233 Gdańsk, Poland
| | - Chris-Kriton Skylaris
- School of Chemistry, University of Southampton, Highfield Campus, Southampton, SO17 1BJ, United Kingdom
| | - Jonathan W. Essex
- School of Chemistry, University of Southampton, Highfield Campus, Southampton, SO17 1BJ, United Kingdom
| |
Collapse
|
29
|
Abstract
Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.
Collapse
|
30
|
Acosta-Silva C, Bertran J, Branchadell V, Oliva A. Kemp Elimination Reaction Catalyzed by Electric Fields. Chemphyschem 2020; 21:295-306. [PMID: 31840917 DOI: 10.1002/cphc.201901155] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2019] [Revised: 12/12/2019] [Indexed: 12/20/2022]
Abstract
The Kemp elimination reaction is the most widely used in the de novo design of new enzymes. The effect of two different kinds of electric fields in the reactions of acetate as a base with benzisoxazole and 5-nitrobenzisoxazole as substrates have been theoretically studied. The effect of the solvent reaction field has been calculated using the SMD continuum model for several solvents; we have shown that solvents inhibit both reactions, the decrease of the reaction rate being larger as far as the dielectric constant is increased. The diminution of the reaction rate is especially remarkable between aprotic organic solvents and protic solvents as water, the electrostatic term of the hydrogen bonds being the main factor for the large inhibitory effect of water. The presence of an external electric field oriented in the direction of the charge transfer (z axis) increases it and, so, the reaction rate. In the reaction of the nitro compound, if the electric field is oriented in an orthogonal direction (x axis) the charge transfer to the NO2 group is favored and there is a subsequent increase of the reaction rate. However, this increase is smaller than the one produced by the field in the z axis. It is worthwhile mentioning that one of the main effects of external electric fields of intermediate intensity is the reorientation of the reactants. Finally, the implications of our results in the de novo design of enzymes are discussed.
Collapse
Affiliation(s)
- Carles Acosta-Silva
- Departament de Química, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Joan Bertran
- Departament de Química, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Vicenç Branchadell
- Departament de Química, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Antoni Oliva
- Departament de Química, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| |
Collapse
|
31
|
Kelly BD, Smith WR. Alchemical Hydration Free-Energy Calculations Using Molecular Dynamics with Explicit Polarization and Induced Polarity Decoupling: An On–the–Fly Polarization Approach. J Chem Theory Comput 2020; 16:1146-1161. [DOI: 10.1021/acs.jctc.9b01139] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
- Braden D. Kelly
- Department of Mathematics and Statistics, University of Guelph, Guelph, ON N1G 2W1, Canada
| | - William R. Smith
- Department of Mathematics and Statistics, University of Guelph, Guelph, ON N1G 2W1, Canada
- Department of Chemical Engineering, University of Waterloo, Waterloo, ON N2L 3G1, Canada
- Faculty of Science, Ontario Tech University, Oshawa, ON L1H 7K4, Canada
| |
Collapse
|
32
|
Milovanović B, Stanković IM, Petković M, Etinski M. Elucidating Solvent Effects on Strong Intramolecular Hydrogen Bond: DFT-MD Study of Dibenzoylmethane in Methanol Solution. Chemphyschem 2019; 20:2852-2859. [PMID: 31544323 DOI: 10.1002/cphc.201900704] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2019] [Revised: 08/31/2019] [Indexed: 11/10/2022]
Abstract
The dynamic aspect of solvation plays a crucial role in determining properties of strong intramolecular hydrogen bonds since solvent fluctuations modify instantaneous hydrogen-bonded proton transfer barriers. Previous studies pointed out that solvent-solute interactions in the first solvation shell govern the position of the proton but the ability of the electric field due to other solvent molecules to localize the proton remains an important issue. In this work, we examine the structure of the O-H⋅⋅⋅O intramolecular hydrogen bond of dibenzoylmethane in methanol solution by employing density functional theory-based molecular dynamics and quantum chemical calculations. Our computations showed that homogeneous electric fields with intensities corresponding to those found in polar solvents are able to considerably alter the proton transfer barrier height in the gas phase. In methanol solution, the proton position is correlated with the difference in electrostatic potentials on the oxygen atoms of dibenzoylmethane even when dibenzoylmethane-methanol hydrogen bonding is lacking. On a timescale of our simulation, the hydrogen bonding and solvent electrostatics tend to localize the proton on different oxygen atoms. These findings provide an insight into the importance of the solvent electric field on the structure of a strong intramolecular hydrogen bond.
Collapse
Affiliation(s)
- Branislav Milovanović
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11000, Belgrade, Serbia
| | | | - Milena Petković
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11000, Belgrade, Serbia
| | - Mihajlo Etinski
- Faculty of Physical Chemistry, University of Belgrade, Studentski trg 12-16, 11000, Belgrade, Serbia
| |
Collapse
|
33
|
Welborn VV, Head-Gordon T. Fluctuations of Electric Fields in the Active Site of the Enzyme Ketosteroid Isomerase. J Am Chem Soc 2019; 141:12487-12492. [DOI: 10.1021/jacs.9b05323] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
- Valerie Vaissier Welborn
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
34
|
Kundu A, Dahms F, Fingerhut BP, Nibbering ETJ, Pines E, Elsaesser T. Hydrated Excess Protons in Acetonitrile/Water Mixtures: Solvation Species and Ultrafast Proton Motions. J Phys Chem Lett 2019; 10:2287-2294. [PMID: 30999753 DOI: 10.1021/acs.jpclett.9b00756] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The solvation structure of protons in aqueous media is highly relevant to electric properties and to proton transport in liquids and membranes. At ambient temperature, polar liquids display structural fluctuations on femto- to picosecond time scales with a direct impact on proton solvation. We use two-dimensional infrared (2D-IR) spectroscopy to follow proton dynamics in acetonitrile/water mixtures with the Zundel cation H5O2+ prepared in neat acetonitrile as a benchmark. The 2D-IR spectra of the proton transfer mode of H5O2+ demonstrate stochastic large-amplitude motions in the double-minimum proton potential, driven by fluctuating electric fields. In all cases, the excess proton is embedded in a water dimer, forming an H5O2+ complex as the major solvation species. This observation is rationalized by quantum mechanics/molecular mechanics molecular dynamics simulations including up to four water molecules embedded in acetonitrile. The Zundel motif interacts with its closest water neighbor in an H7O3+ unit without persistent proton localization.
Collapse
Affiliation(s)
- Achintya Kundu
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| | - Fabian Dahms
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| | - Benjamin P Fingerhut
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| | - Erik T J Nibbering
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| | - Ehud Pines
- Department of Chemistry , Ben Gurion University of the Negev , Beer-Sheva 84105 , Israel
| | - Thomas Elsaesser
- Max-Born-Institut für Nichtlineare Optik und Kurzzeitspektroskopie , Berlin 12489 , Germany
| |
Collapse
|
35
|
Affiliation(s)
- He Yin
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, P. R. China
| | - Hui Li
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, P. R. China
| | - Adam Grofe
- Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, P. R. China
| | - Jiali Gao
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
- Lab of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
- Shenzhen Bay Laboratory, Shenzhen 518055, China
| |
Collapse
|
36
|
Albaugh A, Tuckerman ME, Head-Gordon T. Combining Iteration-Free Polarization with Large Time Step Stochastic-Isokinetic Integration. J Chem Theory Comput 2019; 15:2195-2205. [PMID: 30830768 DOI: 10.1021/acs.jctc.9b00072] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In order to accelerate molecular dynamics simulations using polarizable force fields, we combine a new extended Lagrangian approach that eliminates the self-consistent field step (iEL/0-SCF) with a stochastic integration scheme that allows for a long time step using a multiple time scale algorithm (SIN(R)). We consider different algorithms for the combined scheme that places different components of the nonbonded forces into different time scales, as well as splitting individual nonbonded forces across time scales, to demonstrate that the combined method works well for bulk water as well as for a concentrated salt solution, aqueous peptide, and solvated protein. Depending on system and desired accuracy, the iEL/0-SCF and SIN(R) combination yields lower bound computational speed-ups of ∼6-8 relative to a molecular dynamics Verlet integration using a standard SCF solver implemented in the reference program TINKER 8.1. The combined approach embodies a significant advance for equilibrium simulations in the canonical ensemble of many-body potential energy surfaces for condensed phase systems with speed-ups that exceed what is possible by either method alone.
Collapse
Affiliation(s)
| | - Mark E Tuckerman
- NYU-ECNU , Center for Computational Chemistry at NYU, Shanghai , Shanghai 200062 , China
| | | |
Collapse
|
37
|
Shi MW, Thomas SP, Hathwar VR, Edwards AJ, Piltz RO, Jayatilaka D, Koutsantonis GA, Overgaard J, Nishibori E, Iversen BB, Spackman MA. Measurement of Electric Fields Experienced by Urea Guest Molecules in the 18-Crown-6/Urea (1:5) Host-Guest Complex: An Experimental Reference Point for Electric-Field-Assisted Catalysis. J Am Chem Soc 2019; 141:3965-3976. [PMID: 30761898 DOI: 10.1021/jacs.8b12927] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
High-resolution synchrotron and neutron single-crystal diffraction data of 18-crown-6/(pentakis)urea measured at 30 K are combined, with the aim of better appreciating the electrostatics associated with intermolecular interactions in condensed matter. With two 18-crown-6 molecules and five different urea molecules in the crystal, this represents the most ambitious combined X-ray/synchrotron and neutron experimental charge density analysis to date on a cocrystal or host-guest system incorporating such a large number of unique molecules. The dipole moments of the five urea guest molecules in the crystal are enhanced considerably compared to values determined for isolated molecules, and 2D maps of the electrostatic potential and electric field show clearly how the urea molecules are oriented with dipole moments aligned along the electric field exerted by their molecular neighbors. Experimental electric fields in the range of 10-19 GV m-1, obtained for the five different urea environments, corroborate independent measurements of electric fields in the active sites of enzymes and provide an important experimental reference point for recent discussions focused on electric-field-assisted catalysis.
Collapse
Affiliation(s)
- Ming W Shi
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| | - Sajesh P Thomas
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia.,Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark
| | - Venkatesha R Hathwar
- Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark.,Division of Physics, Faculty of Pure and Applied Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8571 , Japan
| | - Alison J Edwards
- Australian Nuclear Science and Technology Organization , Australian Centre for Neutron Scattering , New Illawarra Road , Lucas Heights , New South Wales 2234 , Australia
| | - Ross O Piltz
- Australian Nuclear Science and Technology Organization , Australian Centre for Neutron Scattering , New Illawarra Road , Lucas Heights , New South Wales 2234 , Australia
| | - Dylan Jayatilaka
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| | - George A Koutsantonis
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| | - Jacob Overgaard
- Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark
| | - Eiji Nishibori
- Division of Physics, Faculty of Pure and Applied Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8571 , Japan
| | - Bo B Iversen
- Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark
| | - Mark A Spackman
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| |
Collapse
|
38
|
Ramos S, Horness RE, Collins JA, Haak D, Thielges MC. Site-specific 2D IR spectroscopy: a general approach for the characterization of protein dynamics with high spatial and temporal resolution. Phys Chem Chem Phys 2019; 21:780-788. [PMID: 30548035 PMCID: PMC6360950 DOI: 10.1039/c8cp06146g] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The conformational heterogeneity and dynamics of protein side chains contribute to function, but investigating exactly how is hindered by experimental challenges arising from the fast timescales involved and the spatial heterogeneity of protein structures. The potential of two-dimensional infrared (2D IR) spectroscopy for measuring conformational heterogeneity and dynamics with unprecedented spatial and temporal resolution has motivated extensive effort to develop amino acids with functional groups that have frequency-resolved absorptions to serve as probes of their protein microenvironments. We demonstrate the full advantage of the approach by selective incorporation of the probe p-cyanophenylalanine at six distinct sites in a Src homology 3 domain and the application of 2D IR spectroscopy to site-specifically characterize heterogeneity and dynamics and their contribution to cognate ligand binding. The approach revealed a wide range of microenvironments and distinct responses to ligand binding, including at the three adjacent, conserved aromatic residues that form the recognition surface of the protein. Molecular dynamics simulations performed for all the labeled proteins provide insight into the underlying heterogeneity and dynamics. Similar application of 2D IR spectroscopy and site-selective probe incorporation will allow for the characterization of heterogeneity and dynamics of other proteins, how heterogeneity and dynamics are affected by solvation and local structure, and how they might contribute to biological function.
Collapse
Affiliation(s)
- Sashary Ramos
- Department of Chemistry, Indiana University, Bloomington, IN 47405, USA.
- Indiana University School of Medicine, Indianapolis, IN 46202, USA
| | - Rachel E. Horness
- Department of Chemistry, Indiana University, Bloomington, IN 47405, USA.
- Indiana University School of Medicine, Indianapolis, IN 46202, USA
| | - Jessica A. Collins
- Department of Chemistry, Indiana University, Bloomington, IN 47405, USA.
- Indiana University School of Medicine, Indianapolis, IN 46202, USA
| | - David Haak
- Indiana University School of Medicine, Indianapolis, IN 46202, USA
| | - Megan C. Thielges
- Department of Chemistry, Indiana University, Bloomington, IN 47405, USA.
- Indiana University School of Medicine, Indianapolis, IN 46202, USA
| |
Collapse
|
39
|
Ghosh A, Cohn B, Prasad AK, Chuntonov L. Quantifying conformations of ester vibrational probes with hydrogen-bond-induced Fermi resonances. J Chem Phys 2018; 149:184501. [PMID: 30441918 DOI: 10.1063/1.5055041] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
Solvatochromic shifts of local vibrational probes report on the strength of the surrounding electric fields and the probe's hydrogen bonding status. Stretching vibrational mode of the ester carbonyl group is a popular solvatochromic reporter used in the studies of peptides and proteins. Small molecules, used to calibrate the response of the vibrational probes, sometimes involve Fermi resonances (FRs) induced by inter-molecular interactions. In the present work, we focus on the scenario where FR does not appear in the infrared spectrum of the ester carbonyl stretching mode in aprotic solvents; however, it is intensified when a hydrogen bond with the reporter is established. When two molecules form hydrogen bonds to the same carbonyl oxygen atom, FR leads to strong hybridization of the involved modes and splitting of the absorption peak. Spectral overlap between the Fermi doublets associated with singly and doubly hydrogen-bonded carbonyl groups significantly complicates quantifying different hydrogen-bonded conformations. We employed a combination of linear and third-order (2DIR) infrared spectroscopy with chemometrics analysis to reveal the individual line shapes and to estimate the occupations of the hydrogen-bonded conformations in methyl acetate, a model small molecule. We identified a hydrogen-bond-induced FR in complexes of methyl acetate with alcohols and water and found that FR is lifted in larger molecules used for control experiments-cholesteryl stearate and methyl cyanoacetate. Applying this methodology to analyze acetonitrile-water solutions revealed that when dissolved in neat water, methyl acetate occupies a single hydrogen-bonding conformation, which is in contrast to the conclusions of previous studies. Our approach can be generally used when FRs prevent direct quantification of the hydrogen bonding status of the vibrational probe.
Collapse
Affiliation(s)
- Anup Ghosh
- Schulich Faculty of Chemistry and Solid State Institute, Technion-Israel Institute of Technology, Haifa 3200003, Israel
| | - Bar Cohn
- Schulich Faculty of Chemistry and Solid State Institute, Technion-Israel Institute of Technology, Haifa 3200003, Israel
| | - Amit K Prasad
- Schulich Faculty of Chemistry and Solid State Institute, Technion-Israel Institute of Technology, Haifa 3200003, Israel
| | - Lev Chuntonov
- Schulich Faculty of Chemistry and Solid State Institute, Technion-Israel Institute of Technology, Haifa 3200003, Israel
| |
Collapse
|
40
|
Błasiak B. One-particle density matrix polarization susceptibility tensors. J Chem Phys 2018; 149:164115. [PMID: 30384720 DOI: 10.1063/1.5051995] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The electric field-induced change in the one-electron density has been expressed as a series of the one-particle density matrix susceptibilities interacting with the spatial distribution of the electric field. The analytic approximate expressions are derived at the Hartree-Fock theory, which serves as a basis for the construction of the generalized model that is designed for an arbitrary form of wavefunction and any type of one-particle density matrix. It is shown that it is possible to accurately predict the changes in the one-electron ground-state density of water molecule in a spatially uniform electric field, as well as in spatially non-uniform electric field distribution generated by point charges. When both linear and quadratic terms with respect to the electric field are accounted for, the electric field-induced polarization energies, dipole moments, and quadrupole moments are quantitatively described by the present theory in electric fields ranging from weak to very strong (0.001-0.07 a.u.). It is believed that the proposed model could open new routes in quantum chemistry for fast and efficient calculations of molecular properties in condensed phases.
Collapse
Affiliation(s)
- Bartosz Błasiak
- Department of Physical and Quantum Chemistry, Faculty of Chemistry, Wrocław University of Science and Technology, Wybrzeże Wyspiańskiego 27, Wrocław 50-370, Poland
| |
Collapse
|
41
|
Biava H, Schreiber T, Katz S, Völler JS, Stolarski M, Schulz C, Michael N, Budisa N, Kozuch J, Utesch T, Hildebrandt P. Long-Range Modulations of Electric Fields in Proteins. J Phys Chem B 2018; 122:8330-8342. [PMID: 30109934 DOI: 10.1021/acs.jpcb.8b03870] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Electrostatic interactions are essential for controlling the protein structure and function. Whereas so far experimental and theoretical efforts focused on the effect of local electrostatics, this work aims at elucidating the long-range modulation of electric fields in proteins upon binding to charged surfaces. The study is based on cytochrome c (Cytc) variants carrying nitrile reporters for the vibrational Stark effect that are incorporated into the protein via genetic engineering and chemical modification. The Cytc variants were thoroughly characterized with respect to possible structural perturbations due to labeling. For the proteins in solution, the relative hydrogen bond occupancy and the calculated electric fields, both obtained from molecular dynamics (MD) simulations, and the experimental nitrile stretching frequencies were used to develop a relationship for separating hydrogen-bonding and non-hydrogen-bonding electric field effects. This relationship provides an excellent description for the stable Cytc variants in solution. For the proteins bound to Au electrodes coated with charged self-assembled monolayers (SAMs), the underlying MD simulations can only account for the electric field changes Δ Eads due to the formation of the electrostatic SAM-Cytc complexes but not for the additional contribution, Δ Eint, representing the consequences of the potential drops over the electrode/SAM/protein interfaces. Both Δ Eads and Δ Eint, determined at distances between 20 and 30 Å with respect to the SAM surface, are comparable in magnitude to the non-hydrogen-bonding electric field in the unbound protein. This long-range modulation of the internal electric field may be of functional relevance for proteins in complexes with partner proteins (Δ Eads) and attached to membranes (Δ Eads + Δ Eint).
Collapse
Affiliation(s)
- Hernan Biava
- Institut für Chemie , Technische Universität Berlin , Sekr. L1, Müller-Breslau-Straße 10 , D-10623 Berlin , Germany
| | - Toni Schreiber
- Institut für Chemie , Technische Universität Berlin , Sekr. PC14, Straße des 17. Juni 135 , D-10623 Berlin , Germany
| | - Sagie Katz
- Institut für Chemie , Technische Universität Berlin , Sekr. PC14, Straße des 17. Juni 135 , D-10623 Berlin , Germany
| | - Jan-Stefan Völler
- Institut für Chemie , Technische Universität Berlin , Sekr. L1, Müller-Breslau-Straße 10 , D-10623 Berlin , Germany
| | - Michael Stolarski
- Institut für Chemie , Technische Universität Berlin , Sekr. PC14, Straße des 17. Juni 135 , D-10623 Berlin , Germany
| | - Claudia Schulz
- Institut für Chemie , Technische Universität Berlin , Sekr. PC14, Straße des 17. Juni 135 , D-10623 Berlin , Germany
| | - Norbert Michael
- Institut für Chemie , Technische Universität Berlin , Sekr. PC14, Straße des 17. Juni 135 , D-10623 Berlin , Germany
| | - Nediljko Budisa
- Institut für Chemie , Technische Universität Berlin , Sekr. L1, Müller-Breslau-Straße 10 , D-10623 Berlin , Germany
| | - Jacek Kozuch
- Institut für Chemie , Technische Universität Berlin , Sekr. PC14, Straße des 17. Juni 135 , D-10623 Berlin , Germany
| | - Tillmann Utesch
- Institut für Chemie , Technische Universität Berlin , Sekr. PC14, Straße des 17. Juni 135 , D-10623 Berlin , Germany
| | - Peter Hildebrandt
- Institut für Chemie , Technische Universität Berlin , Sekr. PC14, Straße des 17. Juni 135 , D-10623 Berlin , Germany
| |
Collapse
|
42
|
Grafton AB, Cheatum CM. Two-dimensional infrared study of the C D and C O stretching vibrations in strongly hydrogen-bonded complexes. Chem Phys 2018. [DOI: 10.1016/j.chemphys.2018.05.025] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
43
|
Boateng HA. Mesh-free hierarchical clustering methods for fast evaluation of electrostatic interactions of point multipoles. J Chem Phys 2018; 147:164104. [PMID: 29096477 DOI: 10.1063/1.4990552] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Electrostatic interactions involving point multipoles are being increasingly implemented to achieve higher accuracy in molecular simulations. A major drawback of multipolar electrostatics is the increased computational cost. Here we develop and compare two Cartesian tree algorithms which employ Taylor approximations and hierarchical clustering to speed up the evaluation of point multipole interactions. We present results from applying the algorithms to compute the free space Coulomb potential and forces of different sets of interacting point multipoles with different densities. The methods achieve high accuracy and speedup of more than an order of magnitude over direct sum calculations and scale well in parallel.
Collapse
Affiliation(s)
- H A Boateng
- Department of Mathematics, Bates College, 2 Andrews Rd., Lewiston, Maine 04240, USA
| |
Collapse
|
44
|
Dalton SR, Vienneau AR, Burstein SR, Xu RJ, Linse S, Londergan CH. Cyanylated Cysteine Reports Site-Specific Changes at Protein-Protein-Binding Interfaces Without Perturbation. Biochemistry 2018; 57:3702-3712. [PMID: 29787228 PMCID: PMC6034165 DOI: 10.1021/acs.biochem.8b00283] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
![]()
To investigate the
cyanylated cysteine vibrational probe group’s
ability to report on binding-induced changes along a protein–protein
interface, the probe group was incorporated at several sites in a
peptide of the calmodulin (CaM)-binding domain of skeletal muscle
myosin light chain kinase. Isothermal titration calorimetry was used
to determine the binding thermodynamics between calmodulin and each
peptide. For all probe positions, the binding affinity was nearly
identical to that of the unlabeled peptide. The CN stretching infrared
band was collected for each peptide free in solution and bound to
calmodulin. Binding-induced shifts in the IR spectral frequencies
were correlated with estimated solvent accessibility based on molecular
dynamics simulations. This work generally suggests (1) that site-specific
incorporation of this vibrational probe group does not cause major
perturbations to its local structural environment and (2) that this
small probe group might be used quite broadly to map dynamic protein-binding
interfaces. However, site-specific perturbations due to artificial
labeling groups can be somewhat unpredictable and should be evaluated
on a site-by-site basis through complementary measurements. A fully
quantitative, simulation-based interpretation of the rich probe IR
spectra is still needed but appears to be possible given recent advances
in simulation techniques.
Collapse
Affiliation(s)
- Shannon R Dalton
- Department of Chemistry , Haverford College , 370 Lancaster Ave , Haverford , Pennsylvania 19041-1392 , United States
| | - Alice R Vienneau
- Department of Chemistry , Haverford College , 370 Lancaster Ave , Haverford , Pennsylvania 19041-1392 , United States
| | - Shana R Burstein
- Department of Chemistry , Haverford College , 370 Lancaster Ave , Haverford , Pennsylvania 19041-1392 , United States
| | - Rosalind J Xu
- Department of Chemistry , Haverford College , 370 Lancaster Ave , Haverford , Pennsylvania 19041-1392 , United States
| | - Sara Linse
- Department of Chemistry and Biochemistry , Lund University , Kemicentrum, Box 118 , 221 00 Lund , Sweden
| | - Casey H Londergan
- Department of Chemistry , Haverford College , 370 Lancaster Ave , Haverford , Pennsylvania 19041-1392 , United States
| |
Collapse
|
45
|
Che F, Gray JT, Ha S, Kruse N, Scott SL, McEwen JS. Elucidating the Roles of Electric Fields in Catalysis: A Perspective. ACS Catal 2018. [DOI: 10.1021/acscatal.7b02899] [Citation(s) in RCA: 142] [Impact Index Per Article: 23.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Fanglin Che
- The Gene and Linda Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, Washington 99164, United States
| | - Jake T. Gray
- The Gene and Linda Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, Washington 99164, United States
| | - Su Ha
- The Gene and Linda Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, Washington 99164, United States
| | - Norbert Kruse
- The Gene and Linda Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, Washington 99164, United States
- Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Susannah L. Scott
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
- Department of Chemical Engineering, University of California, Santa Barbara, California 93106, United States
| | - Jean-Sabin McEwen
- The Gene and Linda Voiland School of Chemical Engineering and Bioengineering, Washington State University, Pullman, Washington 99164, United States
- Institute for Integrated Catalysis, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
- Department of Physics and Astronomy, Washington State University, Pullman, Washington 99164, United States
- Department of Chemistry, Washington State University, Pullman, Washington 99164, United States
- Department of Biological Systems Engineering, Washington State University, Pullman, Washington 99164, United States
| |
Collapse
|
46
|
Mahakrishnan S, Chakraborty S, Vijay A. Anomalies in the equilibrium and nonequilibrium properties of correlated ions in complex molecular environments. Phys Rev E 2018; 96:052133. [PMID: 29347760 DOI: 10.1103/physreve.96.052133] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2017] [Indexed: 11/07/2022]
Abstract
Emergent statistical attributes, and therefore the equations of state, of an assembly of interacting charge carriers embedded within a complex molecular environment frequently exhibit a variety of anomalies, particularly in the high-density (equivalently, the concentration) regime, which are not well understood, because they do not fall under the low-concentration phenomenologies of Debye-Hückel-Onsager and Poisson-Nernst-Planck, including their variants. To go beyond, we here use physical concepts and mathematical tools from quantum scattering theory, transport theory with the Stosszahlansatz of Boltzmann, and classical electrodynamics (Lorentz gauge) and obtain analytical expressions both for the average and the frequency-wave vector-dependent longitudinal and transverse current densities, diffusion coefficient, and the charge density, and therefore the analytical expressions for (a) the chemical potential, activity coefficient, and the equivalent conductivity for strong electrolytes and (b) the current-voltage characteristics for ion-transport processes in complex molecular environments. Using a method analogous to the notion of Debye length and thence the electrical double layer, we here identify a pair of characteristic length scales (longitudinal and the transverse), which, being wave vector and frequency dependent, manifestly exhibit nontrivial fluctuations in space-time. As a unifying theme, we advance a quantity (inverse length dimension), g_{scat}^{(a)}, which embodies all dynamical interactions, through various quantum scattering lengths, relevant to molecular species a, and the analytical behavior which helps us to rationalize the properties of strong electrolytes, including anomalies, in all concentration regimes. As an example, the behavior of g_{scat}^{(a)} in the high-concentration regime explains the anomalous increase of the Debye length with concentration, as seen in a recent experiment on electrolyte solutions. We also put forth an extension of the standard diffusion equation, which manifestly incorporates the effects arising from the underlying microscopic collisions among constituent molecular species. Furthermore, we show a nontrivial connection between the current-voltage characteristics of electrolyte solutions and the Landauer's approach to electrical conduction in mesoscopic solids and thereby establish a definite conceptual bridge between the two disjoint subjects. For numerical insight, we present results on the aqueous solution of KCl as an example of strong electrolyte, and the transport (conduction as well as diffusion) of K^{+} ions in water, as an example of ion transport across the voltage-gated channels in biological cells.
Collapse
Affiliation(s)
- Sathiya Mahakrishnan
- Department of Chemistry, Indian Institute of Technology Madras, Chennai 600036, India
| | - Subrata Chakraborty
- Department of Chemistry, Indian Institute of Technology Madras, Chennai 600036, India
| | - Amrendra Vijay
- Department of Chemistry, Indian Institute of Technology Madras, Chennai 600036, India
| |
Collapse
|
47
|
Warde U, Sekar N. Fluorescent Benzocoumarin-π-Extended Styryl Hybrids: Solvatochromism, Excess Dipole Moment, NLO Properties and DFT Study. J Fluoresc 2017; 28:293-309. [DOI: 10.1007/s10895-017-2192-1] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2017] [Accepted: 10/27/2017] [Indexed: 11/29/2022]
|
48
|
Abstract
The structure and function of biomolecules are strongly influenced by their hydration shells. Structural fluctuations and molecular excitations of hydrating water molecules cover a broad range in space and time, from individual water molecules to larger pools and from femtosecond to microsecond time scales. Recent progress in theory and molecular dynamics simulations as well as in ultrafast vibrational spectroscopy has led to new and detailed insight into fluctuations of water structure, elementary water motions, electric fields at hydrated biointerfaces, and processes of vibrational relaxation and energy dissipation. Here, we review recent advances in both theory and experiment, focusing on hydrated DNA, proteins, and phospholipids, and compare dynamics in the hydration shells to bulk water.
Collapse
Affiliation(s)
- Damien Laage
- École
Normale Supérieure, PSL Research University, UPMC Univ Paris
06, CNRS, Département de Chimie,
PASTEUR, 24 rue Lhomond, 75005 Paris, France
- Sorbonne
Universités, UPMC Univ Paris 06, ENS, CNRS, PASTEUR, 75005 Paris, France
| | - Thomas Elsaesser
- Max-Born-Institut
für Nichtlineare Optik und Kurzzeitspektroskopie, D-12489 Berlin, Germany
| | - James T. Hynes
- École
Normale Supérieure, PSL Research University, UPMC Univ Paris
06, CNRS, Département de Chimie,
PASTEUR, 24 rue Lhomond, 75005 Paris, France
- Sorbonne
Universités, UPMC Univ Paris 06, ENS, CNRS, PASTEUR, 75005 Paris, France
- Department
of Chemistry and Biochemistry, University
of Colorado, Boulder, Colorado 80309, United
States
| |
Collapse
|
49
|
Sun LZ, Zhang JX, Chen SJ. MCTBI: a web server for predicting metal ion effects in RNA structures. RNA (NEW YORK, N.Y.) 2017; 23:1155-1165. [PMID: 28450533 PMCID: PMC5513060 DOI: 10.1261/rna.060947.117] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/29/2017] [Accepted: 04/16/2017] [Indexed: 05/27/2023]
Abstract
Metal ions play critical roles in RNA structure and function. However, web servers and software packages for predicting ion effects in RNA structures are notably scarce. Furthermore, the existing web servers and software packages mainly neglect ion correlation and fluctuation effects, which are potentially important for RNAs. We here report a new web server, the MCTBI server (http://rna.physics.missouri.edu/MCTBI), for the prediction of ion effects for RNA structures. This server is based on the recently developed MCTBI, a model that can account for ion correlation and fluctuation effects for nucleic acid structures and can provide improved predictions for the effects of metal ions, especially for multivalent ions such as Mg2+ effects, as shown by extensive theory-experiment test results. The MCTBI web server predicts metal ion binding fractions, the most probable bound ion distribution, the electrostatic free energy of the system, and the free energy components. The results provide mechanistic insights into the role of metal ions in RNA structure formation and folding stability, which is important for understanding RNA functions and the rational design of RNA structures.
Collapse
Affiliation(s)
- Li-Zhen Sun
- Department of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri 65211, USA
- Department of Applied Physics, Zhejiang University of Technology, Hangzhou 310023, China
| | - Jing-Xiang Zhang
- School of Science and Technology, Zhejiang International Studies University, Hangzhou 310012, China
| | - Shi-Jie Chen
- Department of Physics, Department of Biochemistry, and Informatics Institute, University of Missouri, Columbia, Missouri 65211, USA
| |
Collapse
|
50
|
Albaugh A, Niklasson AMN, Head-Gordon T. Accurate Classical Polarization Solution with No Self-Consistent Field Iterations. J Phys Chem Lett 2017; 8:1714-1723. [PMID: 28350167 DOI: 10.1021/acs.jpclett.7b00450] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
We present a new solution for classical polarization that does not require any self-consistent field iterations, the aspect of classical polarization that makes it computationally expensive. The new method builds upon our iEL/SCF Lagrangian scheme that defines a set of auxiliary induced dipoles whose original purpose was to serve as a time-reversible initial guess to the SCF solution of the set of real induced dipoles. In the new iEL/0-SCF approach the auxiliary dipoles now drive the time evolution of the real induced dipoles such that they stay close to the Born-Oppenheimer surface in order to achieve a truly SCF-less method. We show that the iEL/0-SCF exhibits no loss of simulation accuracy when analyzed across bulk water, low to high concentration salt solutions, and small solutes to large proteins in water. In addition, iEL/0-SCF offers significant computational savings over more expensive SCF calculations based on traditional 1 fs time step integration using symplectic integrators and is as fast as reversible reference system propagator algorithms with an outer 2 fs time step.
Collapse
Affiliation(s)
| | - Anders M N Niklasson
- Theoretical Division, Los Alamos National Laboratory , Los Alamos, New Mexico 87545, United States
| | - Teresa Head-Gordon
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, University of California , Berkeley, California 94720, United States
| |
Collapse
|