1
|
Edholm F, Nandy A, Reinhardt CR, Kastner DW, Kulik HJ. Protein3D: Enabling analysis and extraction of metal-containing sites from the Protein Data Bank with molSimplify. J Comput Chem 2024; 45:352-361. [PMID: 37873926 DOI: 10.1002/jcc.27242] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2023] [Revised: 09/27/2023] [Accepted: 10/03/2023] [Indexed: 10/25/2023]
Abstract
Metalloenzymes catalyze a wide range of chemical transformations, with the active site residues playing a key role in modulating chemical reactivity and selectivity. Unlike smaller synthetic catalysts, a metalloenzyme active site is embedded in a larger protein, which makes interrogation of electronic properties and geometric features with quantum mechanical calculations challenging. Here we implement the ability to fetch crystallographic structures from the Protein Data Bank and analyze the metal binding sites in the program molSimplify. We show the usefulness of the newly created protein3D class to extract the local environment around non-heme iron enzymes containing a two histidine motif and prepare 372 structures for quantum mechanical calculations. Our implementation of protein3D serves to expand the range of systems molSimplify can be used to analyze and will enable high-throughput study of metal-containing active sites in proteins.
Collapse
Affiliation(s)
- Freya Edholm
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
| | - Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
| | - Clorice R Reinhardt
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
| | - David W Kastner
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts, USA
| |
Collapse
|
2
|
Zhou S, Liu Y, Wang S, Wang L. Chemical features and machine learning assisted predictions of protein-ligand short hydrogen bonds. Sci Rep 2023; 13:13741. [PMID: 37612311 PMCID: PMC10447522 DOI: 10.1038/s41598-023-40614-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2023] [Accepted: 08/14/2023] [Indexed: 08/25/2023] Open
Abstract
There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 Å closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.
Collapse
Affiliation(s)
| | - Yuanhao Liu
- Department of Statistics, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA
| | - Sijian Wang
- Department of Statistics, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA.
| | - Lu Wang
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA.
| |
Collapse
|
3
|
Gutermuth T, Sieg J, Stohn T, Rarey M. Modeling with Alternate Locations in X-ray Protein Structures. J Chem Inf Model 2023; 63:2573-2585. [PMID: 37018549 DOI: 10.1021/acs.jcim.3c00100] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/07/2023]
Abstract
In many molecular modeling applications, the standard procedure is still to handle proteins as single, rigid structures. While the importance of conformational flexibility is widely known, handling it remains challenging. Even the crystal structure of a protein usually contains variability exemplified in alternate side chain orientations or backbone segments. This conformational variability is encoded in PDB structure files by so-called alternate locations (AltLocs). Most modeling approaches either ignore AltLocs or resolve them with simple heuristics early on during structure import. We analyzed the occurrence and usage of AltLocs in the PDB and developed an algorithm to automatically handle AltLocs in PDB files enabling all structure-based methods using rigid structures to take the alternative protein conformations described by AltLocs into consideration. A respective software tool named AltLocEnumerator can be used as a structure preprocessor to easily exploit AltLocs. While the amount of data makes it difficult to show impact on a statistical level, handling AltLocs has a substantial impact on a case-by-case basis. We believe that the inspection and consideration of AltLocs is a very valuable approach in many modeling scenarios.
Collapse
Affiliation(s)
- Torben Gutermuth
- Universität Hamburg, ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Jochen Sieg
- Universität Hamburg, ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Tim Stohn
- Universität Hamburg, ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| | - Matthias Rarey
- Universität Hamburg, ZBH - Center for Bioinformatics, Bundesstraße 43, 20146 Hamburg, Germany
| |
Collapse
|
4
|
Villot C, Lao KU. Electronic structure theory on modeling short-range noncovalent interactions between amino acids. J Chem Phys 2023; 158:094301. [PMID: 36889981 DOI: 10.1063/5.0138032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023] Open
Abstract
While short-range noncovalent interactions (NCIs) are proving to be of importance in many chemical and biological systems, these atypical bindings happen within the so-called van der Waals envelope and pose an enormous challenge for current computational methods. We introduce SNCIAA, a database of 723 benchmark interaction energies of short-range noncovalent interactions between neutral/charged amino acids originated from protein x-ray crystal structures at the "gold standard" coupled-cluster with singles, doubles, and perturbative triples/complete basis set [CCSD(T)/CBS] level of theory with a mean absolute binding uncertainty less than 0.1 kcal/mol. Subsequently, a systematic assessment of commonly used computational methods, such as the second-order Møller-Plesset theory (MP2), density functional theory (DFT), symmetry-adapted perturbation theory (SAPT), composite electronic-structure methods, semiempirical approaches, and the physical-based potentials with machine learning (IPML) on SNCIAA is carried out. It is shown that the inclusion of dispersion corrections is essential even though these dimers are dominated by electrostatics, such as hydrogen bondings and salt bridges. Overall, MP2, ωB97M-V, and B3LYP+D4 turned out to be the most reliable methods for the description of short-range NCIs even in strongly attractive/repulsive complexes. SAPT is also recommended in describing short-range NCIs only if the δMP2 correction has been included. The good performance of IPML for dimers at close-equilibrium and long-range conditions is not transferable to the short-range. We expect that SNCIAA will assist the development/improvement/validation of computational methods, such as DFT, force-fields, and ML models, in describing NCIs across entire potential energy surfaces (short-, intermediate-, and long-range NCIs) on the same footing.
Collapse
Affiliation(s)
- Corentin Villot
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, USA
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, USA
| |
Collapse
|
5
|
Kastner DW, Nandy A, Mehmood R, Kulik HJ. Mechanistic Insights into Substrate Positioning That Distinguish Non-heme Fe(II)/α-Ketoglutarate-Dependent Halogenases and Hydroxylases. ACS Catal 2023. [DOI: 10.1021/acscatal.2c06241] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- David W. Kastner
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Rimsha Mehmood
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
6
|
Dutta D, Sharma P, Gomila RM, Frontera A, Barcelo-Oliver M, Verma AK, Baishya T, Bhattacharyya MK. Supramolecular assemblies involving unconventional non-covalent contacts in pyrazole-based coordination compounds of Co(II) and Cu(II) pyridinedicarboxylates: Antiproliferative evaluation and theoretical studies. Polyhedron 2022. [DOI: 10.1016/j.poly.2022.116025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
7
|
Nandy A, Adamji H, Kastner DW, Vennelakanti V, Nazemi A, Liu M, Kulik HJ. Using Computational Chemistry To Reveal Nature’s Blueprints for Single-Site Catalysis of C–H Activation. ACS Catal 2022. [DOI: 10.1021/acscatal.2c02096] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Husain Adamji
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - David W. Kastner
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Biological Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Vyshnavi Vennelakanti
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Azadeh Nazemi
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Mingjie Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
8
|
Yang Z, Hajlasz N, Kulik HJ. Computational Modeling of Conformer Stability in Benenodin-1, a Thermally Actuated Lasso Peptide Switch. J Phys Chem B 2022; 126:3398-3406. [PMID: 35481742 DOI: 10.1021/acs.jpcb.2c00762] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Benenodin-1 is a thermally actuated lasso peptide rotaxane switch with two primary translational isomers that differ in the relative position of the residue Gln15. The conversion from one conformer to the other involves substantial enthalpy-entropy compensation: one conformer is energetically favored and the other is entropically favored. Here, we take a multi-scale quantum mechanical (QM) and classical molecular dynamic (MD) approach to reveal residue-specific sources of these differences in stability. QM reveals that the two benenodin-1 conformers involve distinct hydrogen bonding networks, with the enthalpically favored conformer having more intra-peptide hydrogen bonds between the Gln15 side chain and nearby residues. The evaluation of configurational entropy over the MD-sampled geometries reveals that the entropically favored conformer has enhanced conformational flexibility. By computing the by-residue-sum entropies, we identify the role of Gln15 and neighboring Glu14 in mediating the entropic variation during the switching process. These computational insights help explain the effects of Glu14Ala and Gln15Ala mutations on the conformational population of benenodin-1 observed experimentally.
Collapse
Affiliation(s)
- Zhongyue Yang
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Natalia Hajlasz
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
9
|
Demapan D, Kussmann J, Ochsenfeld C, Cui Q. Factors That Determine the Variation of Equilibrium and Kinetic Properties of QM/MM Enzyme Simulations: QM Region, Conformation, and Boundary Condition. J Chem Theory Comput 2022; 18:2530-2542. [PMID: 35226489 PMCID: PMC9652774 DOI: 10.1021/acs.jctc.1c00714] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
To analyze the impact of various technical details on the results of quantum mechanical (QM)/molecular mechanical (MM) enzyme simulations, including the QM region size, catechol-O-methyltransferase (COMT) is studied as a model system using an approximate QM/MM method (DFTB3/CHARMM). The results show that key equilibrium and kinetic properties for methyl transfer in COMT exhibit limited variations with respect to the size of the QM region, which ranges from ∼100 to ∼500 atoms in this study. With extensive sampling, local and global structural characteristics of the enzyme are largely conserved across the studied QM regions, while the nature of the transition state (e.g., secondary kinetic isotope effect) and reaction exergonicity are largely maintained. Deviations in the free energy profile with different QM region sizes are similar in magnitude to those observed with changes in other simulation protocols, such as different initial enzyme conformations and boundary conditions. Electronic structural properties, such as the covariance matrix of residual charge fluctuations, appear to exhibit rather long-range correlations, especially when the peptide backbone is included in the QM region; this observation holds when a range-separated DFT approach is used as the QM region, suggesting that delocalization error is unlikely the origin. Overall, the analyses suggest that multiple simulation details determine the results of QM/MM enzyme simulations with comparable contributions.
Collapse
Affiliation(s)
- Darren Demapan
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 7 (C), D-81377 Munich, Germany.,Department of Chemistry, University of Wisconsin, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Jörg Kussmann
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 7 (C), D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 7 (C), D-81377 Munich, Germany
| | - Qiang Cui
- Departments of Chemistry, Physics and Biomedical Engineering, Boston University, 590 Commonwealth Avenue, Boston, Massachusetts 02215, United States
| |
Collapse
|
10
|
Zhou S, Liu Y, Wang S, Wang L. Effective prediction of short hydrogen bonds in proteins via machine learning method. Sci Rep 2022; 12:469. [PMID: 35013487 PMCID: PMC8748993 DOI: 10.1038/s41598-021-04306-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Accepted: 12/20/2021] [Indexed: 12/11/2022] Open
Abstract
Short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms lie within 2.7 Å, exhibit prominent quantum mechanical characters and are connected to a wide range of essential biomolecular processes. However, exact determination of the geometry and functional roles of SHBs requires a protein to be at atomic resolution. In this work, we analyze 1260 high-resolution peptide and protein structures from the Protein Data Bank and develop a boosting based machine learning model to predict the formation of SHBs between amino acids. This model, which we name as machine learning assisted prediction of short hydrogen bonds (MAPSHB), takes into account 21 structural, chemical and sequence features and their interaction effects and effectively categorizes each hydrogen bond in a protein to a short or normal hydrogen bond. The MAPSHB model reveals that the type of the donor amino acid plays a major role in determining the class of a hydrogen bond and that the side chain Tyr-Asp pair demonstrates a significant probability of forming a SHB. Combining electronic structure calculations and energy decomposition analysis, we elucidate how the interplay of competing intermolecular interactions stabilizes the Tyr-Asp SHBs more than other commonly observed combinations of amino acid side chains. The MAPSHB model, which is freely available on our web server, allows one to accurately and efficiently predict the presence of SHBs given a protein structure with moderate or low resolution and will facilitate the experimental and computational refinement of protein structures.
Collapse
Affiliation(s)
- Shengmin Zhou
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA
| | - Yuanhao Liu
- Department of Statistics, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA
| | - Sijian Wang
- Department of Statistics, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA
| | - Lu Wang
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, NJ, 08854, USA.
| |
Collapse
|
11
|
Mehmood R, Vennelakanti V, Kulik HJ. Spectroscopically Guided Simulations Reveal Distinct Strategies for Positioning Substrates to Achieve Selectivity in Nonheme Fe(II)/α-Ketoglutarate-Dependent Halogenases. ACS Catal 2021. [DOI: 10.1021/acscatal.1c03169] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Rimsha Mehmood
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Vyshnavi Vennelakanti
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
12
|
Boonserm P, Puthong S, Wichai T, Noitang S, Khunrae P, Sooksai S, Komolpis K. Investigation of major amino acid residues of anti-norfloxacin monoclonal antibodies responsible for binding with fluoroquinolones. Sci Rep 2021; 11:17140. [PMID: 34433868 PMCID: PMC8387498 DOI: 10.1038/s41598-021-96466-6] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2021] [Accepted: 08/10/2021] [Indexed: 12/03/2022] Open
Abstract
It is important to understand the amino acid residues that govern the properties of the binding between antibodies and ligands. We studied the binding of two anti-norfloxacins, anti-nor 132 and anti-nor 155, and the fluoroquinolones norfloxacin, enrofloxacin, ciprofloxacin, and ofloxacin. Binding cross-reactivities tested by an indirect competitive enzyme-linked immunosorbent assay indicated that anti-nor 132 (22–100%) had a broader range of cross-reactivity than anti-nor 155 (62–100%). These cross-reactivities correlated with variations in the numbers of interacting amino acid residues and their positions. Molecular docking was employed to investigate the molecular interactions between the fluoroquinolones and the monoclonal antibodies. Homology models of the heavy chain and light chain variable regions of each mAb 3D structure were docked with the fluoroquinolones targeting the crucial part of the complementarity-determining regions. The fluoroquinolone binding site of anti-nor 155 was a region of the HCDR3 and LCDR3 loops in which hydrogen bonds were formed with TYR (H:35), ASN (H:101), LYS (H:106), ASN (L:92), and ASN (L:93). These regions were further away in anti-nor 132 and could not contact the fluoroquinolones. Another binding region consisting of HIS (L:38) and ASP (H:100) was found for norfloxacin, enrofloxacin, and ciprofloxacin, whereas only ASP (H:100) was found for ofloxacin.
Collapse
Affiliation(s)
- Patamalai Boonserm
- Program in Biotechnology, Faculty of Science, Chulalongkorn University, Bangkok, Thailand
| | - Songchan Puthong
- Institute of Biotechnology and Genetic Engineering, Chulalongkorn University, Bangkok, Thailand
| | - Thanaporn Wichai
- Institute of Biotechnology and Genetic Engineering, Chulalongkorn University, Bangkok, Thailand
| | - Sajee Noitang
- Institute of Biotechnology and Genetic Engineering, Chulalongkorn University, Bangkok, Thailand
| | - Pongsak Khunrae
- King Mongkut's University of Technology Thonburi, Bangkok, Thailand
| | - Sarintip Sooksai
- Institute of Biotechnology and Genetic Engineering, Chulalongkorn University, Bangkok, Thailand.
| | - Kittinan Komolpis
- Institute of Biotechnology and Genetic Engineering, Chulalongkorn University, Bangkok, Thailand. .,Food Risk Hub, Research Unit of Chulalongkorn University, Bangkok, Thailand.
| |
Collapse
|
13
|
Harder, better, faster, stronger: Large-scale QM and QM/MM for predictive modeling in enzymes and proteins. Curr Opin Struct Biol 2021; 72:9-17. [PMID: 34388673 DOI: 10.1016/j.sbi.2021.07.004] [Citation(s) in RCA: 36] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2021] [Revised: 06/25/2021] [Accepted: 07/05/2021] [Indexed: 11/23/2022]
Abstract
Computational prediction of enzyme mechanism and protein function requires accurate physics-based models and suitable sampling. We discuss recent advances in large-scale quantum mechanical (QM) modeling of biochemical systems that have reduced the cost of high-accuracy models. Tradeoffs between sampling and accuracy have motivated modeling with molecular mechanics (MM) in a multiscale QM/MM or iterative approach. Limitations to both conventional density-functional theory and classical MM force fields remain for describing noncovalent interactions in comparison to experiment or wavefunction theory. Because predictions of enzyme action (i.e. electrostatics), free energy barriers, and mechanisms are sensitive to the protocol and embedding method in QM/MM, convergence tests and systematic methods for quantifying QM-level interactions are a needed, active area of development.
Collapse
|
14
|
Takaya D, Watanabe C, Nagase S, Kamisaka K, Okiyama Y, Moriwaki H, Yuki H, Sato T, Kurita N, Yagi Y, Takagi T, Kawashita N, Takaba K, Ozawa T, Takimoto-Kamimura M, Tanaka S, Fukuzawa K, Honma T. FMODB: The World's First Database of Quantum Mechanical Calculations for Biomacromolecules Based on the Fragment Molecular Orbital Method. J Chem Inf Model 2021; 61:777-794. [PMID: 33511845 DOI: 10.1021/acs.jcim.0c01062] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
We developed the world's first web-based public database for the storage, management, and sharing of fragment molecular orbital (FMO) calculation data sets describing the complex interactions between biomacromolecules, named FMO Database (https://drugdesign.riken.jp/FMODB/). Each entry in the database contains relevant background information on how the data was compiled as well as the total energy of each molecular system and interfragment interaction energy (IFIE) and pair interaction energy decomposition analysis (PIEDA) values. Currently, the database contains more than 13 600 FMO calculation data sets, and a comprehensive search function implemented at the front-end. The procedure for selecting target proteins, preprocessing the experimental structures, construction of the database, and details of the database front-end were described. Then, we demonstrated a use of the FMODB by comparing IFIE value distributions of hydrogen bond, ion-pair, and XH/π interactions obtained by FMO method to those by molecular mechanics approach. From the comparison, the statistical analysis of the data provided standard reference values for the three types of interactions that will be useful for determining whether each interaction in a given system is relatively strong or weak compared to the interactions contained within the data in the FMODB. In the final part, we demonstrate the use of the database to examine the contribution of halogen atoms to the binding affinity between human cathepsin L and its inhibitors. We found that the electrostatic term derived by PIEDA greatly correlated with the binding affinities of the halogen containing cathepsin L inhibitors, indicating the importance of QM calculation for quantitative analysis of halogen interactions. Thus, the FMO calculation data in FMODB will be useful for conducting statistical analyses to drug discovery, for conducting molecular recognition studies in structural biology, and for other studies involving quantum mechanics-based interactions.
Collapse
Affiliation(s)
- Daisuke Takaya
- RIKEN Center for Biosystems Dynamics Research, 1-7-22 Suehiro-cho Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Chiduru Watanabe
- RIKEN Center for Biosystems Dynamics Research, 1-7-22 Suehiro-cho Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan.,JST PRESTO, 4-1-8, Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Shunpei Nagase
- RIKEN Center for Biosystems Dynamics Research, 1-7-22 Suehiro-cho Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Kikuko Kamisaka
- RIKEN Center for Biosystems Dynamics Research, 1-7-22 Suehiro-cho Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Yoshio Okiyama
- RIKEN Center for Biosystems Dynamics Research, 1-7-22 Suehiro-cho Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan.,Division of Medicinal Safety Science, National Institute of Health Sciences, 3-25-26 Tonomachi, Kawasaki-ku, Kawasaki, Kanagawa 210-9501, Japan
| | - Hirotomo Moriwaki
- RIKEN Center for Biosystems Dynamics Research, 1-7-22 Suehiro-cho Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Hitomi Yuki
- RIKEN Center for Biosystems Dynamics Research, 1-7-22 Suehiro-cho Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Tomohiro Sato
- RIKEN Center for Biosystems Dynamics Research, 1-7-22 Suehiro-cho Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| | - Noriyuki Kurita
- Department of Computer Science and Engineering, Toyohashi University of Technology, 1-1 Hibarigaoka Tempaku-cho, Toyohashi, Aichi 441-8580, Japan
| | - Yoichiro Yagi
- Graduate School of Engineering, Okayama University of Science, Okayama, 1-1 Ridai-cho, Okayama 700-0005, Japan
| | - Tatsuya Takagi
- Graduate School of Pharmaceutical Sciences, Osaka University, 1-6 Yamadaoka, Suita, Osaka 565-0871, Japan
| | - Norihito Kawashita
- Faculty of Science and Engineering, Kindai University, 3-4-1 Kowakae, Higashiosaka, Osaka 577-8502, Japan
| | - Kenichiro Takaba
- Pharmaceutical Research Center, Laboratory for Medicinal Chemistry, Asahi Kasei Pharma Corporation, 632-1 Mifuku, Izunokuni, Shizuoka 410-2321, Japan
| | - Tomonaga Ozawa
- Kissei Pharmaceutical Co., LTD., Frontier Technology Research Lab., Research Div. 4365-1 Hotaka Kashiwabara, Azumino, Nagano 399-8304, Japan
| | - Midori Takimoto-Kamimura
- Teijin Institute for Biomedical Research, Teijin Pharma Ltd., 4-3-2 Asahigaoka, Hino, Tokyo 191-8512, Japan
| | - Shigenori Tanaka
- Graduate School of System Informatics, Department of Computational Science, Kobe University, 1-1 Rokkodai, Kobe, Hyogo 657-8501, Japan
| | - Kaori Fukuzawa
- School of Pharmacy and Pharmaceutical Sciences, Hoshi University, 2-4-41 Ebara, Shinagawa, Tokyo 142-8501, Japan.,Department of Biomolecular Engineering, Graduate School of Engineering, Tohoku University, 6-6-11 Aoba, Aramaki, Sendai, Miyagi 980-8579, Japan
| | - Teruki Honma
- RIKEN Center for Biosystems Dynamics Research, 1-7-22 Suehiro-cho Tsurumi-ku, Yokohama, Kanagawa 230-0045, Japan
| |
Collapse
|
15
|
Vennelakanti V, Qi HW, Mehmood R, Kulik HJ. When are two hydrogen bonds better than one? Accurate first-principles models explain the balance of hydrogen bond donors and acceptors found in proteins. Chem Sci 2021; 12:1147-1162. [PMID: 35382134 PMCID: PMC8908278 DOI: 10.1039/d0sc05084a] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2020] [Accepted: 11/18/2020] [Indexed: 01/02/2023] Open
Abstract
Hydrogen bonds (HBs) play an essential role in the structure and catalytic action of enzymes, but a complete understanding of HBs in proteins challenges the resolution of modern structural (i.e., X-ray diffraction) techniques and mandates computationally demanding electronic structure methods from correlated wavefunction theory for predictive accuracy. Numerous amino acid sidechains contain functional groups (e.g., hydroxyls in Ser/Thr or Tyr and amides in Asn/Gln) that can act as either HB acceptors or donors (HBA/HBD) and even form simultaneous, ambifunctional HB interactions. To understand the relative energetic benefit of each interaction, we characterize the potential energy surfaces of representative model systems with accurate coupled cluster theory calculations. To reveal the relationship of these energetics to the balance of these interactions in proteins, we curate a set of 4000 HBs, of which >500 are ambifunctional HBs, in high-resolution protein structures. We show that our model systems accurately predict the favored HB structural properties. Differences are apparent in HBA/HBD preference for aromatic Tyr versus aliphatic Ser/Thr hydroxyls because Tyr forms significantly stronger O–H⋯O HBs than N–H⋯O HBs in contrast to comparable strengths of the two for Ser/Thr. Despite this residue-specific distinction, all models of residue pairs indicate an energetic benefit for simultaneous HBA and HBD interactions in an ambifunctional HB. Although the stabilization is less than the additive maximum due both to geometric constraints and many-body electronic effects, a wide range of ambifunctional HB geometries are more favorable than any single HB interaction. Correlated wavefunction theory predicts and high-resolution crystal structure analysis confirms the important, stabilizing effect of simultaneous hydrogen bond donor and acceptor interactions in proteins.![]()
Collapse
Affiliation(s)
- Vyshnavi Vennelakanti
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
- Department of Chemistry
| | - Helena W. Qi
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
- Department of Chemistry
| | - Rimsha Mehmood
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
- Department of Chemistry
| | - Heather J. Kulik
- Department of Chemical Engineering
- Massachusetts Institute of Technology
- Cambridge
- USA
| |
Collapse
|
16
|
Auffinger P, Ennifar E, D'Ascenzo L. Deflating the RNA Mg 2+ bubble. Stereochemistry to the rescue! RNA (NEW YORK, N.Y.) 2020; 27:rna.076067.120. [PMID: 33268500 PMCID: PMC7901845 DOI: 10.1261/rna.076067.120] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/25/2020] [Accepted: 11/20/2020] [Indexed: 05/03/2023]
Abstract
Proper evaluation of the ionic structure of biomolecular systems through X ray and cryo-EM techniques remains challenging but is essential for advancing our understanding of the underlying structure/activity/solvent relationships. However, numerous studies overestimate the number of Mg2+ in deposited structures due to assignment errors finding their origin in improper consideration of stereochemical rules. Herein, to tackle such issues, we re-evaluate the PDBid 6QNR and 6SJ6 models of the ribosome ionic structure. We establish that stereochemical principles need to be carefully pondered when evaluating ion binding features, even when K+ anomalous signals are available as it is the case for the 6QNR PDB entry. For ribosomes, assignment errors can result in misleading conceptions of their solvent structure. For instance, present stereochemical analysis result in a significant decrease of the number of assigned Mg2+ in 6QNR, suggesting that K+ and not Mg2+ is the prevalent ion in the ribosome 1st solvation shell. We stress that the use of proper stereochemical guidelines in combination or not with other identification techniques, such as those pertaining to the detection of transition metals, of some anions and of K+ anomalous signals, is critical for deflating the current Mg2+ bubble witnessed in many ribosome and other RNA structures. We also stress that for the identification of lighter ions such as Mg2+, Na+, …, for which no anomalous signals can be detected, stereochemistry coupled with high resolution structures (<2.4 Å) remain the best currently available option.
Collapse
|
17
|
Wu Z, Newstead S, Biggin PC. The KDEL trafficking receptor exploits pH to tune the strength of an unusual short hydrogen bond. Sci Rep 2020; 10:16903. [PMID: 33037300 PMCID: PMC7547670 DOI: 10.1038/s41598-020-73906-3] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2020] [Accepted: 09/23/2020] [Indexed: 12/27/2022] Open
Abstract
The endoplasmic reticulum (ER) is the main site of protein synthesis in eukaryotic cells and requires a high concentration of luminal chaperones to function. During protein synthesis, ER luminal chaperones are swept along the secretory pathway and must be retrieved to maintain cell viability. ER protein retrieval is achieved by the KDEL receptor, which recognises a C-terminal Lys-Asp-Glu-Leu (KDEL) sequence. Recognition of ER proteins by the KDEL receptor is pH dependent, with binding occurring under acidic conditions in the Golgi and release under conditions of higher pH in the ER. Recent crystal structures of the KDEL receptor in the apo and peptide bound state suggested that peptide binding drives the formation of a short-hydrogen bond that locks the KDEL sequence in the receptor and activates the receptor for COPI binding in the cytoplasm. Using quantum mechanical calculations we demonstrate that the strength of this short hydrogen bond is reinforced following protonation of a nearby histidine, providing a conceptual link between receptor protonation and KDEL peptide binding. Protonation also controls the water networks adjacent to the peptide binding site, leading to a conformational change that ultimately allows the receptor-complex to be recognized by the COPI system.
Collapse
Affiliation(s)
- Zhiyi Wu
- Department of Biochemistry, South Parks Road, Oxford, OX1 3QU, UK
| | - Simon Newstead
- Department of Biochemistry, South Parks Road, Oxford, OX1 3QU, UK.
| | - Philip C Biggin
- Department of Biochemistry, South Parks Road, Oxford, OX1 3QU, UK.
| |
Collapse
|
18
|
Zhou S, Wang L. Quantum effects and 1H NMR chemical shifts of a bifurcated short hydrogen bond. J Chem Phys 2020; 153:114301. [PMID: 32962368 DOI: 10.1063/5.0024734] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
The monoprotonated compound N,N',N''-tris(p-tolyl)azacalix[3](2,6)pyridine (TAPH) contains an intramolecular hydrogen bond that is formed from three N atoms in its cavity. Constrained by the macrocyclic molecular structure, the separations between the N atoms in this bifurcated hydrogen bond are about 2.6 Å, considerably shorter than those typically observed for hydrogen bonded systems in the condensed phases. As such, TAPH exhibits significantly elongated N-H lengths in its hydrogen bond and a downfield 1H NMR chemical shift of 22.1 ppm. In this work, we carry out ab initio molecular dynamics and ab initio path integral molecular dynamics simulations of TAPH in the acetonitrile solution to reveal the geometry and proton sharing conditions of the bifurcated short hydrogen bond and uncover how the interplay of electronic and nuclear quantum effects gives rise to its far downfield 1H chemical shift. Taking a linear short hydrogen bond as a reference, we demonstrate the distinct features of competing quantum effects and electronic shielding effects in the bifurcated hydrogen bond of TAPH. We further use the degree of deshielding on the proton as a measure of the hydrogen bonding interactions and evaluate the strength of the bifurcated short hydrogen bond as compared to its linear counterpart.
Collapse
Affiliation(s)
- Shengmin Zhou
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, New Jersey 08854, USA
| | - Lu Wang
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, New Jersey 08854, USA
| |
Collapse
|
19
|
Andersson CD, Mishra BK, Forsgren N, Ekström F, Linusson A. Physical Mechanisms Governing Substituent Effects on Arene-Arene Interactions in a Protein Milieu. J Phys Chem B 2020; 124:6529-6539. [PMID: 32610016 PMCID: PMC7467712 DOI: 10.1021/acs.jpcb.0c03778] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2020] [Revised: 06/26/2020] [Indexed: 11/28/2022]
Abstract
Arene-arene interactions play important roles in protein-ligand complex formation. Here, we investigate the characteristics of arene-arene interactions between small organic molecules and aromatic amino acids in protein interiors. The study is based on X-ray crystallographic data and quantum mechanical calculations using the enzyme acetylcholinesterase and selected inhibitory ligands as a model system. It is shown that the arene substituents of the inhibitors dictate the strength of the interaction and the geometry of the resulting complexes. Importantly, the calculated interaction energies correlate well with the measured inhibitor potency. Non-hydrogen substituents strengthened all interaction types in the protein milieu, in keeping with results for benzene dimer model systems. The interaction energies were dispersion-dominated, but substituents that induced local dipole moments increased the electrostatic contribution and thus yielded more strongly bound complexes. These findings provide fundamental insights into the physical mechanisms governing arene-arene interactions in the protein milieu and thus into molecular recognition between proteins and small molecules.
Collapse
Affiliation(s)
| | - Brijesh Kumar Mishra
- International
Institute of Information Technology, Bangalore, Karnataka 560003, India
| | - Nina Forsgren
- CBRN
Defense and Security, Swedish Defense Research
Agency, SE-90621 Umeå, Sweden
| | - Fredrik Ekström
- CBRN
Defense and Security, Swedish Defense Research
Agency, SE-90621 Umeå, Sweden
| | - Anna Linusson
- Department
of Chemistry, Umeå University, SE-90187 Umeå, Sweden
| |
Collapse
|
20
|
Mehmood R, Kulik HJ. Both Configuration and QM Region Size Matter: Zinc Stability in QM/MM Models of DNA Methyltransferase. J Chem Theory Comput 2020; 16:3121-3134. [PMID: 32243149 DOI: 10.1021/acs.jctc.0c00153] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Quantum-mechanical/molecular-mechanical (QM/MM) methods are essential to the study of metalloproteins, but the relative importance of sampling and degree of QM treatment in achieving quantitative predictions is poorly understood. We study the relative magnitude of configurational and QM-region sensitivity of energetic and electronic properties in a representative Zn2+ metal binding site of a DNA methyltransferase. To quantify property variations, we analyze snapshots extracted from 250 ns of molecular dynamics simulation. To understand the degree of QM-region sensitivity, we perform analysis using QM regions ranging from a minimal 49-atom region consisting only of the Zn2+ metal and its four coordinating Cys residues up to a 628-atom QM region that includes residues within 12 Å of the metal center. Over the configurations sampled, we observe that illustrative properties (e.g., rigid Zn2+ removal energy) exhibit large fluctuations that are well captured with even minimal QM regions. Nevertheless, for both energetic and electronic properties, we observe a slow approach to asymptotic limits with similarly large changes in absolute values that converge only with larger (ca. 300-atom) QM region sizes. For the smaller QM regions, the electronic description of Zn2+ binding is incomplete: the metal binds too tightly and is too stabilized by the strong electrostatic potential of MM point charges, and the Zn-S bond covalency is overestimated. Overall, this work suggests that efficient sampling with QM/MM in small QM regions is an effective method to explore the influence of enzyme structure on target properties. At the same time, accurate descriptions of electronic and energetic properties require a larger QM region than the minimal metal-coordinating residues in order to converge treatment of both metal-local bonding and the overall electrostatic environment.
Collapse
Affiliation(s)
- Rimsha Mehmood
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
21
|
Zhou S, Wang L. Symmetry and 1H NMR chemical shifts of short hydrogen bonds: impact of electronic and nuclear quantum effects. Phys Chem Chem Phys 2020; 22:4884-4895. [DOI: 10.1039/c9cp06840f] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
Electronic and nuclear quantum effects determine the symmetry and highly downfield 1H NMR chemical shifts of short hydrogen bonds.
Collapse
Affiliation(s)
- Shengmin Zhou
- Department of Chemistry and Chemical Biology
- Institute for Quantitative Biomedicine
- Rutgers University
- Piscataway
- USA
| | - Lu Wang
- Department of Chemistry and Chemical Biology
- Institute for Quantitative Biomedicine
- Rutgers University
- Piscataway
- USA
| |
Collapse
|
22
|
Kuhn B, Gilberg E, Taylor R, Cole J, Korb O. How Significant Are Unusual Protein-Ligand Interactions? Insights from Database Mining. J Med Chem 2019; 62:10441-10455. [PMID: 31730345 DOI: 10.1021/acs.jmedchem.9b01545] [Citation(s) in RCA: 32] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
We present a new approach to derive interaction propensities of protein-ligand atom pairs from mining of the Protein Data Bank. To ensure solid statistics, we use a line-of-sight contact filter and normalize the observed frequency of hits by a statistical null model based on exposed surface areas of atom types in the protein-ligand binding site. This allows us to investigate which intermolecular interactions and geometries are found more often than expected by chance in protein-ligand complexes. We focus our study on some of the unusual interactions that were postulated to be favorable, including σ-hole bonding of halogen and sulfur atoms, weak hydrogen bonding with fluorine as acceptor, and different types of dipolar interactions. Our results confirm some and challenge other common assumptions on these interactions and highlight other contact types that are yet underexplored in structure-based drug design.
Collapse
Affiliation(s)
- Bernd Kuhn
- Roche Pharma Research and Early Development, Roche Innovation Center Basel , F. Hoffmann-La Roche Ltd. , Grenzacherstrasse 124 , CH-4070 Basel , Switzerland
| | - Erik Gilberg
- Roche Pharma Research and Early Development, Roche Innovation Center Basel , F. Hoffmann-La Roche Ltd. , Grenzacherstrasse 124 , CH-4070 Basel , Switzerland
| | - Robin Taylor
- Cambridge Crystallographic Data Centre , 12 Union Road , Cambridge CB2 1EZ , U.K
| | - Jason Cole
- Cambridge Crystallographic Data Centre , 12 Union Road , Cambridge CB2 1EZ , U.K
| | - Oliver Korb
- Roche Pharma Research and Early Development, Roche Innovation Center Basel , F. Hoffmann-La Roche Ltd. , Grenzacherstrasse 124 , CH-4070 Basel , Switzerland
| |
Collapse
|
23
|
Qi HW, Kulik HJ. Reply to "Comment on 'Evaluating Unexpectedly Short Non-covalent Distances in X-ray Crystal Structures of Proteins with Electronic Structure Analysis'". J Chem Inf Model 2019; 59:3609-3610. [PMID: 31424928 DOI: 10.1021/acs.jcim.9b00606] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Helena W Qi
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States.,Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Heather J Kulik
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
24
|
Kruse H, Sponer J, Auffinger P. Comment on “Evaluating Unexpectedly Short Non-covalent Distances in X-ray Crystal Structures of Proteins with Electronic Structure Analysis”. J Chem Inf Model 2019; 59:3605-3608. [DOI: 10.1021/acs.jcim.9b00473] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, CZ-61265 Brno, Czech Republic
| | - Jiri Sponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, CZ-61265 Brno, Czech Republic
| | - Pascal Auffinger
- Architecture et Réactivité de l’ARN, Université de Strasbourg, Institut de Biologie Moléculaire et Cellulaire du CNRS, Strasbourg 67084, France
| |
Collapse
|
25
|
Yang Z, Liu F, Steeves AH, Kulik HJ. Quantum Mechanical Description of Electrostatics Provides a Unified Picture of Catalytic Action Across Methyltransferases. J Phys Chem Lett 2019; 10:3779-3787. [PMID: 31244268 DOI: 10.1021/acs.jpclett.9b01555] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Methyl transferases (MTases) are a well-studied class of enzymes for which competing enzymatic enhancement mechanisms have been suggested, ranging from structural methyl group CH···X hydrogen bonds (HBs) to electrostatic- and charge-transfer-driven stabilization of the transition state (TS). We identified all Class I MTases for which reasonable resolution (<2.0 Å) crystal structures could be used to form catalytically competent ternary complexes for multiscale (i.e., quantum-mechanical/molecular-mechanical or QM/MM) simulation of the SN2 methyl transfer reaction coordinate. The four Class I MTases studied have both distinct functions (e.g., protein repair or biosynthesis) and substrate nucleophiles (i.e., C, N, or O). While CH···X HBs stabilize all reactant complexes, no universal TS stabilization role is found for these interactions in MTases. A consistent picture is instead obtained through analysis of charge transfer and electrostatics, wherein much of cofactor-substrate charge separation is maintained in the TS region, and electrostatic potential is correlated with substrate nucleophilicity (i.e., intrinsic reactivity).
Collapse
Affiliation(s)
- Zhongyue Yang
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Fang Liu
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Adam H Steeves
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Heather J Kulik
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
26
|
Zhou S, Wang L. Unraveling the structural and chemical features of biological short hydrogen bonds. Chem Sci 2019; 10:7734-7745. [PMID: 31588321 PMCID: PMC6764281 DOI: 10.1039/c9sc01496a] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2019] [Accepted: 06/30/2019] [Indexed: 02/06/2023] Open
Abstract
Short hydrogen bonds are ubiquitous in biological macromolecules and exhibit distinctive proton potential energy surfaces and proton sharing properties.
The three-dimensional architecture of biomolecules often creates specialized structural elements, notably short hydrogen bonds that have donor–acceptor separations below 2.7 Å. In this work, we statistically analyze 1663 high-resolution biomolecular structures from the Protein Data Bank and demonstrate that short hydrogen bonds are prevalent in proteins, protein–ligand complexes and nucleic acids. From these biological macromolecules, we characterize the preferred location, connectivity and amino acid composition in short hydrogen bonds and hydrogen bond networks, and assess their possible functional importance. Using electronic structure calculations, we further uncover how the interplay of the structural and chemical features determines the proton potential energy surfaces and proton sharing conditions in biological short hydrogen bonds.
Collapse
Affiliation(s)
- Shengmin Zhou
- Department of Chemistry and Chemical Biology , Institute for Quantitative Biomedicine , Rutgers University , Piscataway , NJ 08854 , USA .
| | - Lu Wang
- Department of Chemistry and Chemical Biology , Institute for Quantitative Biomedicine , Rutgers University , Piscataway , NJ 08854 , USA .
| |
Collapse
|
27
|
Cachau RE, Zhu J, Nicklaus MC. The upcoming subatomic resolution revolution. Curr Opin Struct Biol 2019; 58:53-58. [PMID: 31233975 DOI: 10.1016/j.sbi.2019.05.013] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2019] [Revised: 05/12/2019] [Accepted: 05/13/2019] [Indexed: 10/26/2022]
Abstract
Subatomic resolution macromolecular crystallography has been revealing the most fascinating details of macromolecular structures for many years. This most extreme form of macromolecular crystallography is going through rapid changes. A new generation of superbrilliant X-ray sources and detectors is facilitating the rapid acquisition of high-quality datasets. Equally important, a new breed of methods and highly integrated advanced computational tools for structure refinement and analysis is poised to change the way we use subatomic resolution data and reposition high-resolution macromolecular crystallography in medicinal chemistry studies. Subatomic resolution macromolecular crystallography may soon be a routine source of detailed molecular information besides precise geometries, including binding energies and other chemical descriptors, opening new possibilities of application.
Collapse
Affiliation(s)
- Raul E Cachau
- Advanced Biomedical Computational Science, Frederick National Laboratory for Cancer Research, Leidos Biomedical Inc., Frederick, MD 21702, USA.
| | - Jianghai Zhu
- Advanced Biomedical Computational Science, Frederick National Laboratory for Cancer Research, Leidos Biomedical Inc., Frederick, MD 21702, USA
| | - Marc C Nicklaus
- Chemical Biology Laboratory, National Cancer Institute, Frederick, MD 21702, USA
| |
Collapse
|
28
|
Wahab HA, Amaro RE, Cournia Z. A Celebration of Women in Computational Chemistry. J Chem Inf Model 2019; 59:1683-1692. [DOI: 10.1021/acs.jcim.9b00368] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
| | - Rommie E. Amaro
- Department of Chemistry and Biochemistry, University of California, San Diego, 3234 Urey Hall, #0340, 9500 Gilman Drive, La Jolla, California 92093-0340, United States
| | - Zoe Cournia
- Biomedical Research Foundation, Academy of Athens, 11527 Athens, Greece
| |
Collapse
|
29
|
Mehmood R, Qi HW, Steeves AH, Kulik HJ. The Protein’s Role in Substrate Positioning and Reactivity for Biosynthetic Enzyme Complexes: The Case of SyrB2/SyrB1. ACS Catal 2019. [DOI: 10.1021/acscatal.9b00865] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|