1
|
Chaturvedi SS, Vargas S, Ajmera P, Alexandrova AN. Directed Evolution of Protoglobin Optimizes the Enzyme Electric Field. J Am Chem Soc 2024. [PMID: 38848547 DOI: 10.1021/jacs.4c03914] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2024]
Abstract
To unravel why computational design fails in creating viable enzymes, while directed evolution (DE) succeeds, our research delves into the laboratory evolution of protoglobin. DE has adapted this protein to efficiently catalyze carbene transfer reactions. We show that the previously proposed enhanced substrate access and binding alone cannot account for increased yields during DE. The 3D electric field in the entire active site is tracked through protein dynamics, clustered using the affinity propagation algorithm, and subjected to principal component analysis. This analysis reveals notable changes in the electric field with DE, where distinct field topologies influence transition state energetics and mechanism. A chemically meaningful field component emerges and takes the lead during DE and facilitates crossing the barrier to carbene transfer. Our findings underscore intrinsic electric field dynamic's influence on enzyme function, the ability of the field to switch mechanisms within the same protein, and the crucial role of the field in enzyme design.
Collapse
Affiliation(s)
- Shobhit S Chaturvedi
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| | - Santiago Vargas
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| | - Pujan Ajmera
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| | - Anastassia N Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| |
Collapse
|
2
|
Xu J, Hao J, Bu C, Meng Y, Xiao H, Zhang M, Li C. XMECP: Reaching State-of-the-Art MECP Optimization in Multiscale Complex Systems. J Chem Theory Comput 2024; 20:3590-3600. [PMID: 38651739 DOI: 10.1021/acs.jctc.4c00033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
The Python-based program, XMECP, is developed for realizing robust, efficient, and state-of-the-art minimum energy crossing point (MECP) optimization in multiscale complex systems. This article introduces the basic capabilities of the XMECP program by theoretically investigating the MECP mechanism of several example systems including (1) the photosensitization mechanism of benzophenone, (2) photoinduced proton-coupled electron transfer in the cytosine-guanine base pair in DNA, (3) the spin-flip process in oxygen activation catalyzed by an iron-containing 2-oxoglutarate-dependent oxygenase (Fe/2OGX), and (4) the photochemical pathway of flavoprotein adjusted by the intensity of an external electric field. MECPs related to multistate reaction and multistate reactivity in large-scale complex biochemical systems can be well-treated by workflows suggested by the XMECP program. The branching plane updating the MECP optimization algorithm is strongly recommended as it provides derivative coupling vector (DCV) with explicit calculation and can equivalently evaluate contributions from non-QM residues to DCV, which can be nonadiabatic coupling or spin-orbit coupling in different cases. In the discussed QM/MM examples, we also found that the influence on the QM region by DCV can occur through noncovalent interactions and decay with distance. In the example of DNA base pairs, the nonadiabatic coupling occurs across the π-π stacking structure formed in the double-helix system. In contrast to general intuition, in the example of Fe/2OGX, the central ferrous and oxygen part contribute little to the spin-orbit coupling; however, a nearby arginine residue, which is treated by molecular mechanics in the QM/MM method, contributes significantly via two hydrogen bonds formed with α-ketoglutarate (α-KG). This indicates that the arginine residue plays a significant role in oxygen activation, driving the initial triplet state toward the productive quintet state, which is more than the previous knowledge that the arginine residue can bind α-KG at the reaction site by hydrogen bonds.
Collapse
Affiliation(s)
- Jiawei Xu
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
- School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Jian Hao
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
- School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Caijie Bu
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
- College of Chemistry and Materials Science, Fujian Normal University, Fuzhou 350117, Fujian, P. R. China
| | - Yajie Meng
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
- School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Han Xiao
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
| | - Minyi Zhang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
| | - Chunsen Li
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, Fuzhou 350002, Fujian, P. R. China
- School of Chemical Sciences, University of Chinese Academy of Sciences, Beijing 100049, P. R. China
- Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Xiamen University, Xiamen 361005, Fujian, P. R. China
| |
Collapse
|
3
|
Wang H, Abe I. Recent developments in the enzymatic modifications of steroid scaffolds. Org Biomol Chem 2024; 22:3559-3583. [PMID: 38639195 DOI: 10.1039/d4ob00327f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/20/2024]
Abstract
Steroids are an important family of bioactive compounds. Steroid drugs are renowned for their multifaceted pharmacological activities and are the second-largest category in the global pharmaceutical market. Recent developments in biocatalysis and biosynthesis have led to the increased use of enzymes to enhance the selectivity, efficiency, and sustainability for diverse modifications of steroids. This review discusses the advancements achieved over the past five years in the enzymatic modifications of steroid scaffolds, focusing on enzymatic hydroxylation, reduction, dehydrogenation, cascade reactions, and other modifications for future research on the synthesis of novel steroid compounds and related drugs, and new therapeutic possibilities.
Collapse
Affiliation(s)
- Huibin Wang
- Graduate School of Pharmaceutical Sciences, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan.
| | - Ikuro Abe
- Graduate School of Pharmaceutical Sciences, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan.
- Collaborative Research Institute for Innovative Microbiology, The University of Tokyo, Yayoi 1-1-1, Bunkyo-ku, Tokyo 113-8657, Japan
| |
Collapse
|
4
|
Zheng C, Ji Z, Mathews II, Boxer SG. Enhanced active-site electric field accelerates enzyme catalysis. Nat Chem 2023; 15:1715-1721. [PMID: 37563323 PMCID: PMC10906027 DOI: 10.1038/s41557-023-01287-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2023] [Accepted: 06/29/2023] [Indexed: 08/12/2023]
Abstract
The design and improvement of enzymes based on physical principles remain challenging. Here we demonstrate that the principle of electrostatic catalysis can be leveraged to substantially improve a natural enzyme's activity. We enhanced the active-site electric field in horse liver alcohol dehydrogenase by replacing the serine hydrogen-bond donor with threonine and replacing the catalytic Zn2+ with Co2+. Based on the electric field enhancement, we make a quantitative prediction of rate acceleration-50-fold faster than the wild-type enzyme-which was in close agreement with experimental measurements. The effects of the hydrogen bonding and metal coordination, two distinct chemical forces, are described by a unified physical quantity-electric field, which is quantitative, and shown here to be additive and predictive. These results suggest a new design paradigm for both biological and non-biological catalysts.
Collapse
Affiliation(s)
- Chu Zheng
- Department of Chemistry, Stanford University, Stanford, CA, USA
| | - Zhe Ji
- Department of Chemistry, Stanford University, Stanford, CA, USA
| | | | - Steven G Boxer
- Department of Chemistry, Stanford University, Stanford, CA, USA.
| |
Collapse
|
5
|
Pan X, Van R, Pu J, Nam K, Mao Y, Shao Y. Free Energy Profile Decomposition Analysis for QM/MM Simulations of Enzymatic Reactions. J Chem Theory Comput 2023; 19:8234-8244. [PMID: 37943896 PMCID: PMC10835707 DOI: 10.1021/acs.jctc.3c00973] [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/12/2023]
Abstract
In enzyme mechanistic studies and mutant design, it is highly desirable to know the individual residue contributions to the reaction free energy and barrier. In this work, we show that such free energy contributions from each residue can be readily obtained by postprocessing ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulation trajectories. Specifically, through a mean force integration along the minimum free energy pathway, one can obtain the electrostatic, polarization, and van der Waals contributions from each residue to the free energy barrier. Separately, a similar analysis procedure allows us to assess the contribution from different collective variables along the reaction coordinate. The chorismate mutase reaction is used to demonstrate the utilization of these two trajectory analysis tools.
Collapse
Affiliation(s)
- Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Richard Van
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
- Laboratory of Computational Biology, National, Heart, Lung and Blood Institute, National Institutes of Health, Bethesda, Maryland 20824, United States
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis, Indianapolis, Indiana 46202, United States
| | - Kwangho Nam
- Department of Chemistry and Biochemistry, University of Texas at Arlington, Arlington, Texas 76019, United States
| | - Yuezhi Mao
- Department of Chemistry and Biochemistry, San Diego State University, San Diego, California 92182, United States
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| |
Collapse
|
6
|
Zhu C, Pham LN, Yuan X, Ouyang H, Coote ML, Zhang X. High Electric Fields on Water Microdroplets Catalyze Spontaneous and Fast Reactions in Halogen-Bond Complexes. J Am Chem Soc 2023; 145:21207-21212. [PMID: 37724917 DOI: 10.1021/jacs.3c08818] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/21/2023]
Abstract
The use of external electric fields as green and efficient catalysts in synthetic chemistry has recently received significant attention for their ability to deliver remarkable control of reaction selectivity and acceleration of reaction rates. Technically, methods of generating high electric fields in the range of 1-10 V/nm are limited, as in-vacuo techniques have obvious scalability issues. The spontaneous high fields at various interfaces promise to solve this problem. In this study, we take advantage of the spontaneous high electric field at the air-water interface of sprayed water microdroplets in the reactions of several halogen bond systems: Nu:--X-X, where Nu: is pyridine or quinuclidine and X is bromine or iodine. The field facilitates ultrafast electron transfer from Nu:, yielding a Nu-X covalent bond and causing the X-X bond to cleave. This reaction occurs in microseconds in microdroplets but takes days to weeks in bulk solution. Density functional theory calculations predict that the reaction becomes barrier-free in the presence of oriented external electric fields, supporting the notion that the electric fields in the water droplets are responsible for the catalysis. We anticipate that microdroplet chemistry will be an avenue rich in opportunities in the reactions facilitated by high electric fields and provides an alternative way to tackle the scalability problem.
Collapse
Affiliation(s)
- Chenghui Zhu
- College of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education), Renewable Energy Conversion and Storage Centre, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Frontiers Science Centre for New Organic Matter, Nankai University, Tianjin, 300071, China
- Beijing National Laboratory for Molecular Sciences, Beijing, 100190, China
| | - Le Nhan Pham
- Institute for Nanoscale Science & Technology, Flinders University, Adelaide, South Australia 5042, Australia
| | - Xu Yuan
- College of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education), Renewable Energy Conversion and Storage Centre, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Frontiers Science Centre for New Organic Matter, Nankai University, Tianjin, 300071, China
- Beijing National Laboratory for Molecular Sciences, Beijing, 100190, China
| | - Haoran Ouyang
- College of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education), Renewable Energy Conversion and Storage Centre, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Frontiers Science Centre for New Organic Matter, Nankai University, Tianjin, 300071, China
- Beijing National Laboratory for Molecular Sciences, Beijing, 100190, China
| | - Michelle L Coote
- Institute for Nanoscale Science & Technology, Flinders University, Adelaide, South Australia 5042, Australia
| | - Xinxing Zhang
- College of Chemistry, Key Laboratory of Advanced Energy Materials Chemistry (Ministry of Education), Renewable Energy Conversion and Storage Centre, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Frontiers Science Centre for New Organic Matter, Nankai University, Tianjin, 300071, China
- Beijing National Laboratory for Molecular Sciences, Beijing, 100190, China
| |
Collapse
|
7
|
Demchenko AP. Proton transfer reactions: from photochemistry to biochemistry and bioenergetics. BBA ADVANCES 2023. [DOI: 10.1016/j.bbadva.2023.100085] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/12/2023] Open
|
8
|
Eberhart ME, Wilson TR, Johnston NW, Alexandrova AN. Geometry of Charge Density as a Reporter on the Role of the Protein Scaffold in Enzymatic Catalysis: Electrostatic Preorganization and Beyond. J Chem Theory Comput 2023; 19:694-704. [PMID: 36562645 DOI: 10.1021/acs.jctc.2c01060] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Enzymes host active sites inside protein macromolecules, which have diverse, often incredibly complex, and atom-expensive structures. It is an outstanding question what the role of these expensive scaffolds might be in enzymatic catalysis. Answering this question is essential to both enzymology and the design of artificial enzymes with proficiencies that will match those of the best natural enzymes. Protein rigidifying the active site, contrasted with the dynamics and vibrational motion promoting the reaction, as well as long-range electrostatics (also known as electrostatic preorganization) were all proposed as central contributions of the scaffold to the catalysis. Here, we show that all these effects inevitably produce changes in the quantum mechanical electron density in the active site, which in turn defines the reactivity. The phenomena are therefore fundamentally inseparable. The geometry of the electron density-a scalar field characterized by a number of mathematical features such as critical points-is a rigorous and convenient descriptor of enzymatic catalysis and a reporter on the role of the protein. We show how this geometry can be analyzed, linked to the reaction barriers, and report in particular on intramolecular electric fields in enzymes. We illustrate these tools on the studies of electrostatic preorganization in several representative enzyme classes, both natural and artificial. We highlight the forward-looking aspects of the approach.
Collapse
Affiliation(s)
- Mark E Eberhart
- Department of Chemistry, Colorado School of Mines, 1500 Illinois Street, Golden, Colorado 80401, United States
| | - Timothy R Wilson
- Department of Chemistry, Colorado School of Mines, 1500 Illinois Street, Golden, Colorado 80401, United States
| | - Nathaniel W Johnston
- Department of Chemistry, and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095, United States
| | - Anastassia N Alexandrova
- Department of Chemistry, and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095, United States
| |
Collapse
|
9
|
Brunner syndrome caused by point mutation explained by multiscale simulation of enzyme reaction. Sci Rep 2022; 12:21889. [PMID: 36536002 PMCID: PMC9763434 DOI: 10.1038/s41598-022-26296-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2022] [Accepted: 12/13/2022] [Indexed: 12/23/2022] Open
Abstract
Brunner syndrome is a disorder characterized by intellectual disability and impulsive, aggressive behavior associated with deficient function of the monoamine oxidase A (MAO-A) enzyme. These symptoms (along with particularly high serotonin levels) have been reported in patients with two missense variants in MAO-A (p.R45W and p.E446K). Herein, we report molecular simulations of the rate-limiting step of MAO-A-catalyzed serotonin degradation for these variants. We found that the R45W mutation causes a 6000-fold slowdown of enzymatic function, whereas the E446K mutation causes a 450-fold reduction of serotonin degradation rate, both of which are practically equivalent to a gene knockout. In addition, we thoroughly compared the influence of enzyme electrostatics on the catalytic function of both the wild type MAO-A and the p.R45W variant relative to the wild type enzyme, revealing that the mutation represents a significant electrostatic perturbation that contributes to the barrier increase. Understanding genetic disorders is closely linked to understanding the associated chemical mechanisms, and our research represents a novel attempt to bridge the gap between clinical genetics and the underlying chemical physics.
Collapse
|
10
|
Vaissier Welborn V, Archer WR, Schulz MD. Characterizing Ion-Polymer Interactions in Aqueous Environment with Electric Fields. J Chem Inf Model 2022; 63:2030-2036. [PMID: 36533730 DOI: 10.1021/acs.jcim.2c01048] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Polymers make the basis of highly tunable materials that could be designed and optimized for metal recovery from aqueous environments. While experimental studies show that this approach has potential, it suffers from a limited knowledge of the detailed molecular interaction between polymers and target metal ions. Here, we propose to calculate intrinsic electric fields from polarizable force field molecular dynamics simulations to characterize the driving force behind Eu3+ motion in the presence of poly(ethylenimine methylenephosphonate), a specifically designed metal chelating polymer. Focusing on the metal chelation initiation step (i.e., before binding), we can rationalize the role of each molecule on ion dynamics by projecting these electric fields along the direction of ion motion. We find that the polymer functional groups act indirectly, and the polymer-metal ion interaction is actually mediated by water. This result is consistent with the experimental observation that metal sequestration by these polymers is entropically driven. This study suggests that electric field calculations can help the design of metal chelating polymers, for example, by seeking to optimize polymer-solvent interactions rather than polymer-ion interactions.
Collapse
Affiliation(s)
- Valerie Vaissier Welborn
- Department of Chemistry, Macromolecules Innovation Institute (MII), Virginia Tech, Blacksburg, Virginia24060, United States
| | - William R. Archer
- Department of Chemistry, Macromolecules Innovation Institute (MII), Virginia Tech, Blacksburg, Virginia24060, United States
| | - Michael D. Schulz
- Department of Chemistry, Macromolecules Innovation Institute (MII), Virginia Tech, Blacksburg, Virginia24060, United States
| |
Collapse
|
11
|
Ji Z, Boxer SG. β-Lactamases Evolve against Antibiotics by Acquiring Large Active-Site Electric Fields. J Am Chem Soc 2022; 144:22289-22294. [PMID: 36399691 PMCID: PMC10075085 DOI: 10.1021/jacs.2c10791] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
A compound bound covalently to an enzyme active site can act either as a substrate if the covalent linkage is readily broken up by the enzyme or as an inhibitor if the bond dissociates slowly. We tracked the reactivity of such bonds associated with the rise of the resistance to penicillin G (PenG) in protein evolution from penicillin-binding proteins (PBPs) to TEM β-lactamases and with the development of avibactam (Avb) to overcome the resistance. We found that the ester linkage in PBP-PenG is resistant to hydrolysis mainly due to the small electric fields present in the protein active site. Conversely, the same linkage in the descendant TEM-PenG experiences large electric fields that stabilize the more charge-separated transition state and thus lower the free energy barrier to hydrolysis. Specifically, the electric fields were improved from -59 to -140 MV/cm in an ancient evolution dating back billions of years, contributing 5 orders of magnitude rate acceleration. This trend continues today in the nullification of newly developed antibiotic drugs. The fast linkage hydrolysis acquired from evolution is counteracted by the upgrade of PenG to Avb whose linkage escapes from the hydrolysis by returning to a low-field environment. Using the framework of electrostatic catalysis, the electric field, an observable from vibrational spectroscopy, provides a unifying physical metric to understand protein evolution and to guide the design of covalent drugs.
Collapse
Affiliation(s)
- Zhe Ji
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Steven G Boxer
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
12
|
Wilson TR, Morgenstern A, Alexandrova AN, Eberhart ME. Bond Bundle Analysis of Ketosteroid Isomerase. J Phys Chem B 2022; 126:9443-9456. [PMID: 36383139 DOI: 10.1021/acs.jpcb.2c03638] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
Abstract
Bond bundle analysis is used to investigate enzymatic catalysis in the ketosteroid isomerase (KSI) active site. We identify the unique bonding regions in five KSI systems, including those exposed to applied oriented electric fields and those with amino acid mutations, and calculate the precise redistribution of electron density and other regional properties that accompanies either enhancement or inhibition of KSI catalytic activity. We find that catalytic enhancement results from promoting both inter- and intra-molecular electron density redistribution, between bond bundles and bond wedges within the KSI-docked substrate molecule, in the forward direction of the catalyzed reaction. Though the redistribution applies to both types of perturbed systems and is thus suggestive of a general catalytic role, we observe that bond properties (e.g., volume vs energy vs electron count) can respond independently and disproportionately depending on the type of perturbation. We conclude that the resulting catalytic enhancement/inhibition proceeds via different mechanisms, where some bond properties are utilized more by one type of perturbation than the other. Additionally, we find that the correlations between bond wedge properties and catalyzed reaction barrier energies are additive to predict those of bond bundles and atomic basins, providing a rigorous grounding for connecting changes in local charge density to resulting shifts in reaction barrier energy.
Collapse
Affiliation(s)
- Timothy R Wilson
- Department of Chemistry, Colorado School of Mines, 1500 Illinois Street, Golden, Colorado 80004, United States
| | - Amanda Morgenstern
- Department of Chemistry & Biochemistry, UCCS, 1420 Austin Bluffs Pkwy, Colorado Springs, Colorado 80918, United States
| | - Anastassia N Alexandrova
- Department of Chemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095, United States
| | - M E Eberhart
- Department of Chemistry, Colorado School of Mines, 1500 Illinois Street, Golden, Colorado 80004, United States
| |
Collapse
|
13
|
Diao W, Yan S, Farrell JD, Wang B, Ye F, Wang Z. Preorganized Internal Electric Field Powers Catalysis in the Active Site of Uracil-DNA Glycosylase. ACS Catal 2022. [DOI: 10.1021/acscatal.2c02886] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Wenwen Diao
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China
| | - Shengheng Yan
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - James D. Farrell
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Binju Wang
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Fangfu Ye
- Wenzhou Institute, University of Chinese Academy of Sciences, Wenzhou, Zhejiang 325000, China
- Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China
- Oujiang Laboratory (Zhejiang Lab for Regenerative Medicine, Vision and Brain Health), Wenzhou, Zhejiang 325001, China
| | - Zhanfeng Wang
- Center for Advanced Materials Research, Beijing Normal University, Zhuhai 519087, China
| |
Collapse
|
14
|
Wright D, Sangtarash S, Mueller NS, Lin Q, Sadeghi H, Baumberg JJ. Vibrational Stark Effects: Ionic Influence on Local Fields. J Phys Chem Lett 2022; 13:4905-4911. [PMID: 35623089 PMCID: PMC9189927 DOI: 10.1021/acs.jpclett.2c01048] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/11/2023]
Abstract
Molecules containing vibrational Stark shift reporters provide a useful tool for measuring DC electric fields in situ. To quantify this effect theoretically, density functional theory (DFT) calculations are usually utilized in a uniform electric field. However, using a combined theoretical and experimental study, we demonstrate here that uniform field DFT cannot simultaneously model the behavior of the three strongest vibrational modes in molecules forming a monolayer on an electrode. We show, by directly modeling ionic movement, that the measured Stark shifts are explained by partial electrical double-layer penetration into the molecular layer. This effect is sensitive to the local environment, and the Stark shifts can be fully suppressed experimentally by introducing a mixed molecular layer that prevents ionic double-layer penetration.
Collapse
Affiliation(s)
- Demelza Wright
- NanoPhotonics
Centre, Department of Physics, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, U.K.
| | - Sara Sangtarash
- Device
Modelling Group, School of Engineering, University of Warwick, Coventry CV4 7AL, U.K.
| | - Niclas S. Mueller
- NanoPhotonics
Centre, Department of Physics, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, U.K.
| | - Qianqi Lin
- NanoPhotonics
Centre, Department of Physics, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, U.K.
| | - Hatef Sadeghi
- Device
Modelling Group, School of Engineering, University of Warwick, Coventry CV4 7AL, U.K.
| | - Jeremy J. Baumberg
- NanoPhotonics
Centre, Department of Physics, Cavendish Laboratory, University of Cambridge, Cambridge CB3 0HE, U.K.
| |
Collapse
|
15
|
Lawal MM, Vaissier Welborn V. Structural dynamics support electrostatic interactions in the active site of Adenylate Kinase. Chembiochem 2022; 23:e202200097. [PMID: 35303385 DOI: 10.1002/cbic.202200097] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2022] [Revised: 03/17/2022] [Indexed: 11/12/2022]
Abstract
Electrostatic preorganization as well as structural and dynamic heterogeneity are often used to rationalize the remarkable catalytic efficiency of enzymes. However, they are often presented as incompatible because the generation of permanent electrostatic effects implies that the protein structure remains rigid. Here, we use a metric, electric fields, that can treat electrostatic contributions and dynamics effects on equal footing, for a unique perspective on enzymatic catalysis. We find that the residues that contribute the most to electrostatic interactions with the substrate in the active site of Adenylate Kinase (our working example) are also the most flexible residues. Further, entropy-tuning mutations raise flexibility at the picosecond timescale where more conformations can be visited on short time periods, thereby softening the sharp heterogeneity normally visible at the microsecond timescale.
Collapse
Affiliation(s)
| | - Valerie Vaissier Welborn
- Virginia Polytechnic Institute and State University, Chemistry, Davidson 421A, 1040 Drillfield Drive, 24073, Blacksburg, UNITED STATES
| |
Collapse
|
16
|
Why Monoamine Oxidase B Preferably Metabolizes N-Methylhistamine over Histamine: Evidence from the Multiscale Simulation of the Rate-Limiting Step. Int J Mol Sci 2022; 23:ijms23031910. [PMID: 35163835 PMCID: PMC8836602 DOI: 10.3390/ijms23031910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2022] [Revised: 02/04/2022] [Accepted: 02/05/2022] [Indexed: 11/18/2022] Open
Abstract
Histamine levels in the human brain are controlled by rather peculiar metabolic pathways. In the first step, histamine is enzymatically methylated at its imidazole Nτ atom, and the produced N-methylhistamine undergoes an oxidative deamination catalyzed by monoamine oxidase B (MAO-B), as is common with other monoaminergic neurotransmitters and neuromodulators of the central nervous system. The fact that histamine requires such a conversion prior to oxidative deamination is intriguing since MAO-B is known to be relatively promiscuous towards monoaminergic substrates; its in-vitro oxidation of N-methylhistamine is about 10 times faster than that for histamine, yet this rather subtle difference appears to be governing the decomposition pathway. This work clarifies the MAO-B selectivity toward histamine and N-methylhistamine by multiscale simulations of the rate-limiting hydride abstraction step for both compounds in the gas phase, in aqueous solution, and in the enzyme, using the established empirical valence bond methodology, assisted by gas-phase density functional theory (DFT) calculations. The computed barriers are in very good agreement with experimental kinetic data, especially for relative trends among systems, thereby reproducing the observed MAO-B selectivity. Simulations clearly demonstrate that solvation effects govern the reactivity, both in aqueous solution as well as in the enzyme although with an opposing effect on the free energy barrier. In the aqueous solution, the transition-state structure involving histamine is better solvated than its methylated analog, leading to a lower barrier for histamine oxidation. In the enzyme, the higher hydrophobicity of N-methylhistamine results in a decreased number of water molecules at the active side, leading to decreased dielectric shielding of the preorganized catalytic electrostatic environment provided by the enzyme. This renders the catalytic environment more efficient for N-methylhistamine, giving rise to a lower barrier relative to histamine. In addition, the transition state involving N-methylhistamine appears to be stabilized by the surrounding nonpolar residues to a larger extent than with unsubstituted histamine, contributing to a lower barrier with the former.
Collapse
|
17
|
Hao H, Leven I, Head-Gordon T. Can electric fields drive chemistry for an aqueous microdroplet? Nat Commun 2022; 13:280. [PMID: 35022410 PMCID: PMC8755715 DOI: 10.1038/s41467-021-27941-x] [Citation(s) in RCA: 87] [Impact Index Per Article: 43.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2021] [Accepted: 12/14/2021] [Indexed: 11/20/2022] Open
Abstract
Reaction rates of common organic reactions have been reported to increase by one to six orders of magnitude in aqueous microdroplets compared to bulk solution, but the reasons for the rate acceleration are poorly understood. Using a coarse-grained electron model that describes structural organization and electron densities for water droplets without the expense of ab initio methods, we investigate the electric field distributions at the air-water interface to understand the origin of surface reactivity. We find that electric field alignments along free O-H bonds at the surface are ~16 MV/cm larger on average than that found for O-H bonds in the interior of the water droplet. Furthermore, electric field distributions can be an order of magnitude larger than the average due to non-linear coupling of intramolecular solvent polarization with intermolecular solvent modes which may contribute to even greater surface reactivity for weakening or breaking chemical bonds at the droplet surface.
Collapse
Affiliation(s)
- Hongxia Hao
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, University of California, Berkeley, CA, 94720, USA
- Pitzer Center for Theoretical Chemistry, University of California, Berkeley, CA, 94720, USA
- Departments of Chemistry, University of California, Berkeley, CA, 94720, USA
| | - Itai Leven
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, University of California, Berkeley, CA, 94720, USA
- Pitzer Center for Theoretical Chemistry, University of California, Berkeley, CA, 94720, USA
- Departments of Chemistry, University of California, Berkeley, CA, 94720, USA
| | - Teresa Head-Gordon
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, University of California, Berkeley, CA, 94720, USA.
- Pitzer Center for Theoretical Chemistry, University of California, Berkeley, CA, 94720, USA.
- Departments of Chemistry, University of California, Berkeley, CA, 94720, USA.
- Departments of Bioengineering, University of California, Berkeley, CA, 94720, USA.
- Departments of Chemical and Biomolecular Engineering, University of California, Berkeley, CA, 94720, USA.
| |
Collapse
|
18
|
Wilson TR, Alexandrova AN, Eberhart ME. Electron Density Geometry and the Quantum Theory of Atoms in Molecules. J Phys Chem A 2021; 125:10622-10631. [PMID: 34905923 DOI: 10.1021/acs.jpca.1c09359] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A novel form of charge density analysis, that of isosurface curvature redistribution, is formulated and applied to the toy problem of carbonyl oxygen activation in formaldehyde. The isosurface representation of the electron charge density allows us to incorporate the rigorous geometric constraints of closed surfaces toward the analysis and chemical interpretation of the charge density response to perturbations. Visual inspection of 2D isosurface motion resulting from applied external electric fields reveals how the isosurface curvature flows within and between atoms and that a molecule can be uniquely and completely partitioned into chemically significant regions of positive and negative curvatures. These concepts reveal that carbonyl oxygen activation proceeds primarily through curvature and charge redistribution within rather than between Bader atoms. Using gradient bundle analysis─the partitioning of formaldehyde into infinitesimal volume elements bounded by QTAIM zero-flux surfaces─the observations from visual isosurface inspection are verified. The results of the formaldehyde carbonyl analysis are then shown to be transferable to the substrate carbonyl in the ketosteroid isomerase enzyme, laying the groundwork for extending this approach to the problems of enzymatic catalysis.
Collapse
Affiliation(s)
- Timothy R Wilson
- Department of Chemistry, Colorado School of Mines, Golden, Colorado 80401, United States
| | | | - M E Eberhart
- Department of Chemistry, Colorado School of Mines, Golden, Colorado 80401, United States
| |
Collapse
|
19
|
Prah A, Mavri J, Stare J. An electrostatic duel: subtle differences in the catalytic performance of monoamine oxidase A and B isoenzymes elucidated at the residue level using quantum computations. Phys Chem Chem Phys 2021; 23:26459-26467. [PMID: 34806105 DOI: 10.1039/d1cp03993h] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The origin of the immense catalytic power of enzymes remains one of the biggest unresolved questions in biochemistry, with electrostatics being one of the main contenders. Herein, we report results that not only confirm that electrostatics is the driving force behind enzyme catalysis, but also that it is capable of tuning subtle differences in the catalytic performance between structurally similar enzymes, as demonstrated using the example of isoenzymes, monoamine oxidases A and B. Using our own computationally efficient multiscale model [A. Prah, et al., ACS Catal., 2019, 9, 1231] we analyzed the rate-limiting step of the reaction between phenylethylamine and both isoenzymes and deduced that the electrostatic environment provided by isoenzyme B has a perceivably higher catalytic influence on all the considered parameters of the reaction (energy barrier, charge transfer, dipole moment, and HOMO-LUMO gap). This is in full agreement with the available experimental kinetic data and with our own simulations of the reaction in question. In-depth analysis of individual amino acid contributions of both isoenzymes to the barrier (based on the interaction between the electric field provided by the enzyme and the dipole moment of the reacting moiety) shows that the majority of the difference between the isoenzymes can be attributed to a small number of sizable differences between the aligned amino acid pairs, whereas in most of the pairs the difference in contribution to the barrier is vanishingly small. These results suggest that electrostatics largely controls the substrate selectivity of enzymes and validates our approach as being capable of discerning fine nuances in the selectivity of structurally related isoenzymes.
Collapse
Affiliation(s)
- Alja Prah
- Theory Department, National Institute of Chemistry, Slovenia. .,University of Ljubljana, Faculty of Pharmacy, Slovenia
| | - Janez Mavri
- Theory Department, National Institute of Chemistry, Slovenia.
| | - Jernej Stare
- Theory Department, National Institute of Chemistry, Slovenia.
| |
Collapse
|
20
|
Kędzierski P, Moskal M, Sokalski WA. Catalytic Fields as a Tool to Analyze Enzyme Reaction Mechanism Variants and Reaction Steps. J Phys Chem B 2021; 125:11606-11616. [PMID: 34648705 PMCID: PMC8558854 DOI: 10.1021/acs.jpcb.1c05256] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2021] [Revised: 09/19/2021] [Indexed: 11/29/2022]
Abstract
Catalytic fields representing the topology of the optimal molecular environment charge distribution that reduces the activation barrier have been used to examine alternative reaction variants and to determine the role of conserved catalytic residues for two consecutive reactions catalyzed by the same enzyme. Until now, most experimental and conventional top-down theoretical studies employing QM/MM or ONIOM methods have focused on the role of enzyme electric fields acting on broken bonds of reactants. In contrast, our bottom-up approach dealing with a small reactant and transition-state model allows the analysis of the opposite effects: how the catalytic field resulting from the charge redistribution during the enzyme reaction acts on conserved amino acid residues and contributes to the reduction of the activation barrier. This approach has been applied to the family of histidyl tRNA synthetases involved in the translation of the genetic code into the protein amino acid sequence. Activation energy changes related to conserved charged amino acid residues for 12 histidyl tRNA synthetases from different biological species allowed to compare on equal footing the catalytic residues involved in ATP aminoacylation and tRNA charging reactions and to analyze different reaction mechanisms proposed in the literature. A scan of the library of atomic multipoles for amino acid side-chain rotamers within the catalytic field pointed out the change in the Glu83 conformation as the critical catalytic effect, providing, at low computational cost, insight into the electrostatic preorganization of the enzyme catalytic site at a level of detail that has not yet been accessible in conventional experimental or theoretical methods. This opens the way for rational reverse biocatalyst design at a very limited computational cost without resorting to empirical methods.
Collapse
Affiliation(s)
- Paweł Kędzierski
- Department of Chemistry, Wrocław
University of Science and Technology, Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - Martyna Moskal
- Department of Chemistry, Wrocław
University of Science and Technology, Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland
| | - W. Andrzej Sokalski
- Department of Chemistry, Wrocław
University of Science and Technology, Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland
| |
Collapse
|
21
|
Advances in optimizing enzyme electrostatic preorganization. Curr Opin Struct Biol 2021; 72:1-8. [PMID: 34280872 DOI: 10.1016/j.sbi.2021.06.006] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Revised: 05/28/2021] [Accepted: 06/05/2021] [Indexed: 12/19/2022]
Abstract
Utilizing electric fields to catalyze chemical reactions is not a new idea, but in enzymology it undergoes a renaissance, inspired by Warhsel's concept of electrostatic preorganization. According to this concept, the source of the immense catalytic efficiency of enzymes is the intramolecular electric field that permanently favors the reaction transition state over the reactants. Within enzyme design, computational efforts have fallen short in designing enzymes with natural-like efficacy. The outcome could improve if long-range electrostatics (often omitted in current protocols) would be optimized. Here, we highlight the major developments in methods for analyzing and designing electric fields generated by the protein scaffolds, in order to both better understand how natural enzymes function, and aid artificial enzyme design.
Collapse
|
22
|
Bím D, Alexandrova AN. Local Electric Fields as a Natural Switch of Heme-Iron Protein Reactivity. ACS Catal 2021; 11:6534-6546. [PMID: 34413991 DOI: 10.1021/acscatal.1c00687] [Citation(s) in RCA: 32] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Abstract
Heme-iron oxidoreductases operating through the high-valent FeIVO intermediates perform crucial and complicated transformations, such as oxidations of unreactive saturated hydrocarbons. These enzymes share the same Fe coordination, only differing by the axial ligation, e.g., Cys in P450 oxygenases, Tyr in catalases, and His in peroxidases. By examining ~200 heme-iron proteins, we show that the protein hosts exert highly specific intramolecular electric fields on the active sites, and there is a strong correlation between the direction and magnitude of this field and the protein function. In all heme proteins, the field is preferentially aligned with the Fe-O bond ( Fz ). The Cys-ligated P450 oxygenases have the highest average Fz of 28.5 MV cm-1, i.e., most enhancing the oxyl-radical character of the oxo group, and consistent with the ability of these proteins to activate strong C-H bonds. In contrast, in Tyr-ligated proteins, the average Fz is only 3.0 MV cm-1, apparently suppressing single-electron off-pathway oxidations, and in His-ligated proteins, Fz is -8.7 MV cm-1. The operational field range is given by the trade-off between the low reactivity of the FeIVO Compound I at the more negative Fz , and the low selectivity at the more positive Fz . Consequently, a heme-iron site placed in the field characteristic of another heme-iron protein class loses its canonical function, and gains an adverse one. Thus, electric fields produced by the protein scaffolds, together with the nature of the axial ligand, control all heme-iron chemistry.
Collapse
Affiliation(s)
- Daniel Bím
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095-1569, United States
| | - Anastassia N. Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095-1569, United States
- California NanoSystems Institute, University of California, Los Angeles, 570 Westwood Plaza, Los Angeles, California 90095-1569, United States
| |
Collapse
|
23
|
Hennefarth MR, Alexandrova AN. Heterogeneous Intramolecular Electric Field as a Descriptor of Diels–Alder Reactivity. J Phys Chem A 2021; 125:1289-1298. [DOI: 10.1021/acs.jpca.1c00181] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
Affiliation(s)
- Matthew R. Hennefarth
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095-1569, United States
| | - Anastassia N. Alexandrova
- Department of Chemistry and Biochemistry, University of California, Los Angeles, 607 Charles E. Young Drive East, Los Angeles, California 90095-1569, United States
- California NanoSystems Institute, University of California, Los Angeles, 570 Westwood Plaza, Los Angeles, California 90095-1569, United States
| |
Collapse
|
24
|
Li WL, Head-Gordon T. Catalytic Principles from Natural Enzymes and Translational Design Strategies for Synthetic Catalysts. ACS CENTRAL SCIENCE 2021; 7:72-80. [PMID: 33532570 PMCID: PMC7844850 DOI: 10.1021/acscentsci.0c01556] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/18/2020] [Indexed: 05/19/2023]
Abstract
As biocatalysts, enzymes are characterized by their high catalytic efficiency and strong specificity but are relatively fragile by requiring narrow and specific reactive conditions for activity. Synthetic catalysts offer an opportunity for more chemical versatility operating over a wider range of conditions but currently do not reach the remarkable performance of natural enzymes. Here we consider some new design strategies based on the contributions of nonlocal electric fields and thermodynamic fluctuations to both improve the catalytic step and turnover for rate acceleration in arbitrary synthetic catalysts through bioinspired studies of natural enzymes. With a focus on the enzyme as a whole catalytic construct, we illustrate the translational impact of natural enzyme principles to synthetic enzymes, supramolecular capsules, and electrocatalytic surfaces.
Collapse
Affiliation(s)
- Wan-Lu Li
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Kenneth
S. Pitzer Center for Theoretical Chemistry, University of California Berkeley, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Chemical
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- Kenneth
S. Pitzer Center for Theoretical Chemistry, University of California Berkeley, Berkeley, California 94720, United States
- Department of Chemistry, Department of Chemical and Biomolecular Engineering, and Department of
Bioengineering, University of California
Berkeley, Berkeley, California 94720, United States
| |
Collapse
|