1
|
Szell PMJ, Zablotny S, Bryce DL. Halogen bonding as a supramolecular dynamics catalyst. Nat Commun 2019; 10:916. [PMID: 30796220 PMCID: PMC6385366 DOI: 10.1038/s41467-019-08878-8] [Citation(s) in RCA: 50] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2018] [Accepted: 02/04/2019] [Indexed: 01/23/2023] Open
Abstract
Dynamic processes have many implications in functional molecules, including catalysts, enzymes, host-guest complexes, and molecular machines. Here, we demonstrate via deuterium NMR relaxation experiments how halogen bonding directly impacts the dynamics in solid 2,3,5,6-tetramethylpyrazine cocrystals, catalyzing the methyl group rotation. On average, we observe a reduction of 56% in the rotational activation energy of the methyl groups in the halogen bonded cocrystals, contrasting the reduction of 36% in the hydrogen bonded cocrystals, with respect to pure 2,3,5,6-tetramethylpyrazine. Density functional theory calculations attribute this superior catalytic ability of the halogen bond to the simultaneous destabilization of the staggered conformation and stabilization of the gauche conformation, overall reducing the rotational energy barrier. Furthermore, the calculations suggest that the catalytic ability of the halogen bond may be tuneable, with stronger halogen bond donors acting as superior dynamics catalysts. Thus, halogen bonding may play a role in both assembly and promoting dynamical processes. The halogen bond is well known for its ability to assemble supramolecules. Here, using NMR experiments, the authors reveal the role of these bonds in dynamic processes, finding that the halogen bond directly catalyzes dynamical rotation in solid cocrystals by reducing the associated energy barrier.
Collapse
Affiliation(s)
- Patrick M J Szell
- Department of Chemistry and Biomolecular Sciences & Centre for Catalysis Research and Innovation, University of Ottawa, Ottawa, ON, K1N 6N5, Canada
| | - Scott Zablotny
- Department of Chemistry and Biomolecular Sciences & Centre for Catalysis Research and Innovation, University of Ottawa, Ottawa, ON, K1N 6N5, Canada
| | - David L Bryce
- Department of Chemistry and Biomolecular Sciences & Centre for Catalysis Research and Innovation, University of Ottawa, Ottawa, ON, K1N 6N5, Canada.
| |
Collapse
|
2
|
Rajeshwar T R, Krishnan M. Direct Determination of Site-Specific Noncovalent Interaction Strengths of Proteins from NMR-Derived Fast Side Chain Motional Parameters. J Phys Chem B 2017; 121:5174-5186. [PMID: 28452484 DOI: 10.1021/acs.jpcb.7b01402] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
A novel approach to accurately determine residue-specific noncovalent interaction strengths (ξ) of proteins from NMR-measured fast side chain motional parameters (Oaxis2) is presented. By probing the environmental sensitivity of side chain conformational energy surfaces of individual residues of a diverse set of proteins, the microscopic connections between ξ, Oaxis2, conformational entropy (Sconf), conformational barriers, and rotamer stabilities established here are found to be universal among proteins. The results reveal that side chain flexibility and conformational entropy of each residue decrease with increasing ξ and that for each residue type there exists a critical range of ξ, determined primarily by the mean side chain conformational barriers, within which flexibility of any residue can be reversibly tuned from highly flexible (with Oaxis2 ∼ 0) to highly restricted (with Oaxis2 ∼ 1) by increasing ξ by ∼3 kcal/mol. Beyond this critical range of ξ, both side chain flexibility and conformational entropy are insensitive to ξ. The interrelationships between conformational dynamics, conformational entropy, and noncovalent interactions of protein side chains established here open up new avenues to probe perturbation-induced (for example, ligand-binding, temperature, pressure) changes in fast side chain dynamics and thermodynamics of proteins by comparing their conformational energy surfaces in the native and perturbed states.
Collapse
Affiliation(s)
- Rajitha Rajeshwar T
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology , Gachibowli, Hyderabad 500 032, India
| | - Marimuthu Krishnan
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology , Gachibowli, Hyderabad 500 032, India
| |
Collapse
|
3
|
Zandarashvili L, Nguyen D, Anderson KM, White MA, Gorenstein DG, Iwahara J. Entropic Enhancement of Protein-DNA Affinity by Oxygen-to-Sulfur Substitution in DNA Phosphate. Biophys J 2016; 109:1026-37. [PMID: 26331260 DOI: 10.1016/j.bpj.2015.07.032] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2015] [Revised: 07/07/2015] [Accepted: 07/15/2015] [Indexed: 10/23/2022] Open
Abstract
Dithioation of DNA phosphate is known to enhance binding affinities, at least for some proteins. We mechanistically characterized this phenomenon for the Antennapedia homeodomain-DNA complex by integrated use of fluorescence, isothermal titration calorimetry, NMR spectroscopy, and x-ray crystallography. By fluorescence and isothermal titration calorimetry, we found that this affinity enhancement is entropy driven. By NMR, we investigated the ionic hydrogen bonds and internal motions of lysine side-chain NH3(+) groups involved in ion pairs with DNA. By x-ray crystallography, we compared the structures of the complexes with and without dithioation of the phosphate. Our NMR and x-ray data show that the lysine side chain in contact with the DNA phosphate becomes more dynamic upon dithioation. Our thermodynamic, structural, and dynamic investigations collectively suggest that the affinity enhancement by the oxygen-to-sulfur substitution in DNA phosphate is largely due to an entropic gain arising from mobilization of the intermolecular ion pair at the protein-DNA interface.
Collapse
Affiliation(s)
- Levani Zandarashvili
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas
| | - Dan Nguyen
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas
| | - Kurtis M Anderson
- Department of NanoMedicine and Biomedical Engineering and Institute of Molecular Medicine, University of Texas Health Science Center at Houston, Houston, Texas
| | - Mark A White
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas
| | - David G Gorenstein
- Department of NanoMedicine and Biomedical Engineering and Institute of Molecular Medicine, University of Texas Health Science Center at Houston, Houston, Texas
| | - Junji Iwahara
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas.
| |
Collapse
|
4
|
Hong L, Petridis L, Smith JC. Biomolecular Structure and Dynamics with Neutrons: The View from Simulation. Isr J Chem 2014. [DOI: 10.1002/ijch.201300137] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
5
|
Jones EM, Balakrishnan G, Squier TC, Spiro TG. Distinguishing unfolding and functional conformational transitions of calmodulin using ultraviolet resonance Raman spectroscopy. Protein Sci 2014; 23:1094-101. [PMID: 24895328 DOI: 10.1002/pro.2495] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2014] [Revised: 05/27/2014] [Accepted: 05/28/2014] [Indexed: 01/09/2023]
Abstract
Calmodulin (CaM) is a ubiquitous moderator protein for calcium signaling in all eukaryotic cells. This small calcium-binding protein exhibits a broad range of structural transitions, including domain opening and folding-unfolding, that allow it to recognize a wide variety of binding partners in vivo. While the static structures of CaM associated with its various binding activities are fairly well-known, it has been challenging to examine the dynamics of transition between these structures in real-time, due to a lack of suitable spectroscopic probes of CaM structure. In this article, we examine the potential of ultraviolet resonance Raman (UVRR) spectroscopy for clarifying the nature of structural transitions in CaM. We find that the UVRR spectral change (with 229 nm excitation) due to thermal unfolding of CaM is qualitatively different from that associated with opening of the C-terminal domain in response to Ca(2+) binding. This spectral difference is entirely due to differences in tertiary contacts at the interdomain tyrosine residue Tyr138, toward which other spectroscopic methods are not sensitive. We conclude that UVRR is ideally suited to identifying the different types of structural transitions in CaM and other proteins with conformation-sensitive tyrosine residues, opening a path to time-resolved studies of CaM dynamics using Raman spectroscopy.
Collapse
Affiliation(s)
- Eric M Jones
- Department of Chemistry, University of Washington, Seattle, Washington, 98195-1700
| | | | | | | |
Collapse
|
6
|
T. RR, Smith JC, Krishnan M. Hidden Regularity and Universal Classification of Fast Side Chain Motions in Proteins. J Am Chem Soc 2014; 136:8590-605. [DOI: 10.1021/ja5024783] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Rajitha Rajeshwar T.
- Center
for Computational Natural Sciences and Bioinformatics, International
Institute of Information Technology, Gachibowli, Hyderabad 500 032, India
| | - Jeremy C. Smith
- UT/ORNL
Center for Molecular Biophysics, Oak Ridge National Laboratory, P.O. Box 2008 Oak Ridge, Tennessee 37831-6309, United States
- Department
of Biochemistry and Molecular and Cellular Biology, University of Tennessee, M407 Walters Life Sciences, 1414 Cumberland Avenue, Knoxville, Tennessee 37996, United States
| | - Marimuthu Krishnan
- Center
for Computational Natural Sciences and Bioinformatics, International
Institute of Information Technology, Gachibowli, Hyderabad 500 032, India
| |
Collapse
|
7
|
Glass DC, Krishnan M, Smith JC, Baudry J. Three entropic classes of side chain in a globular protein. J Phys Chem B 2013; 117:3127-34. [PMID: 23421556 DOI: 10.1021/jp400564q] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The relationship between the NMR methyl group axial order parameter and the side chain conformational entropy is investigated in inhibitor-bound and apo human HIV protease using molecular dynamics simulation. Three distinct entropic classes of methyl-bearing side chains, determined by the topological distance of the methyl group from the protein backbone (i.e., the number of χ-bonds between the Cα and the carbon of the CH3 group), are revealed by atomistic trajectory analyses performed in the local frame of reference of individual methyl probes. The results demonstrate that topologically equivalent methyl groups experience similar nonbonded microenvironments regardless of the type of residues to which they are attached. Similarly, methyl groups that belong to the same side chain but that are not topologically equivalent exhibit different thermodynamic and dynamic properties. The two-parameter classification (based upon entropy and methyl axial order parameter) of side chains described here permits improved estimates of the conformational entropies of proteins from NMR motional parameters.
Collapse
Affiliation(s)
- Dennis C Glass
- UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, P.O. Box 2008 Oak Ridge, Tennessee 37831-6309, United States
| | | | | | | |
Collapse
|
8
|
Anderson KM, Esadze A, Manoharan M, Brüschweiler R, Gorenstein DG, Iwahara J. Direct observation of the ion-pair dynamics at a protein-DNA interface by NMR spectroscopy. J Am Chem Soc 2013; 135:3613-9. [PMID: 23406569 PMCID: PMC3721336 DOI: 10.1021/ja312314b] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Ion pairing is one of the most fundamental chemical interactions and is essential for molecular recognition by biological macromolecules. From an experimental standpoint, very little is known to date about ion-pair dynamics in biological macromolecular systems. Absorption, infrared, and Raman spectroscopic methods were previously used to characterize dynamic properties of ion pairs, but these methods can be applied only to small compounds. Here, using NMR (15)N relaxation and hydrogen-bond scalar (15)N-(31)P J-couplings ((h3)J(NP)), we have investigated the dynamics of the ion pairs between lysine side-chain NH3(+) amino groups and DNA phosphate groups at the molecular interface of the HoxD9 homeodomain-DNA complex. We have determined the order parameters and the correlation times for C-N bond rotation and reorientation of the lysine NH3(+) groups. Our data indicate that the NH3(+) groups in the intermolecular ion pairs are highly dynamic at the protein-DNA interface, which should lower the entropic costs for protein-DNA association. Judging from the C-N bond-rotation correlation times along with experimental and quantum-chemically derived (h3)J(NP) hydrogen-bond scalar couplings, it seems that breakage of hydrogen bonds in the ion pairs occurs on a sub-nanosecond time scale. Interestingly, the oxygen-to-sulfur substitution in a DNA phosphate group was found to enhance the mobility of the NH3(+) group in the intermolecular ion pair. This can partially account for the affinity enhancement of the protein-DNA association by the oxygen-to-sulfur substitution, which is a previously observed but poorly understood phenomenon.
Collapse
Affiliation(s)
- Kurtis M. Anderson
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, TX 77555
- Department of NanoMedicine and Biomedical Engineering, Institute of Molecular Medicine, University of Texas Health Science Center at Houston, Houston, TX 77225
| | - Alexandre Esadze
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, TX 77555
| | - Mariappan Manoharan
- Chemical Sciences Laboratory, Department of Chemistry and Biochemistry, and National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306
| | - Rafael Brüschweiler
- Chemical Sciences Laboratory, Department of Chemistry and Biochemistry, and National High Magnetic Field Laboratory, Florida State University, Tallahassee, Florida 32306
| | - David G. Gorenstein
- Department of NanoMedicine and Biomedical Engineering, Institute of Molecular Medicine, University of Texas Health Science Center at Houston, Houston, TX 77225
| | - Junji Iwahara
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, TX 77555
| |
Collapse
|
9
|
Krishnan M, Smith JC. Reconstruction of protein side-chain conformational free energy surfaces from NMR-derived methyl axis order parameters. J Phys Chem B 2012; 116:4124-33. [PMID: 22401582 DOI: 10.1021/jp2104853] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
An analytical approach is developed for reconstructing site-specific methyl-bearing protein side-chain conformational energy surfaces from NMR methyl axis order parameters (O(axis)(2)). Application of an enhanced sampling algorithm (adaptive biasing force) to molecular dynamics simulation of a protein, calcium-bound calmodulin, reveals a nonlinear correlation between O(axis)(2) and the populations of rotamer states of protein side-chains, permitting the rotamer populations to be extracted directly from O(axis)(2). The analytical approach yields side-chain conformational distributions that are in excellent agreement with those obtained from the enhanced-sampling MD results.
Collapse
Affiliation(s)
- Marimuthu Krishnan
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500 032, India.
| | | |
Collapse
|
10
|
Zhang S, Xu L, Dong J, Cheng P, Zhou Z, Fu J. Collision-induced dissociation of singly and doubly charged CuII–cytidine complexes in the gas phase: an experimental and computational study. RSC Adv 2012. [DOI: 10.1039/c2ra01293f] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
|
11
|
Polyansky AA, Zubac R, Zagrovic B. Estimation of conformational entropy in protein-ligand interactions: a computational perspective. Methods Mol Biol 2012; 819:327-53. [PMID: 22183546 DOI: 10.1007/978-1-61779-465-0_21] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Conformational entropy is an important component of the change in free energy upon binding of a ligand to its target protein. As a consequence, development of computational techniques for reliable estimation of conformational entropies is currently receiving an increased level of attention in the context of computational drug design. Here, we review the most commonly used techniques for conformational entropy estimation from classical molecular dynamics simulations. Although by-and-large still not directly used in practical drug design, these techniques provide a golden standard for developing other, computationally less-demanding methods for such applications, in addition to furthering our understanding of protein-ligand interactions in general. In particular, we focus on the quasi-harmonic approximation and discuss different approaches that can be used to go beyond it, most notably, when it comes to treating anharmonic and/or correlated motions. In addition to reviewing basic theoretical formalisms, we provide a concrete set of steps required to successfully calculate conformational entropy from molecular dynamics simulations, as well as discuss a number of practical issues that may arise in such calculations.
Collapse
Affiliation(s)
- Anton A Polyansky
- Laboratory of Computational Biophysics, Department of Structural and Computational Biology, Max Perutz Laboratories, Vienna, Austria
| | | | | |
Collapse
|
12
|
Hembree WI, Baudry J. Three-dimensional mapping of microenvironmental control of methyl rotational barriers. J Phys Chem B 2011; 115:8575-80. [PMID: 21634366 DOI: 10.1021/jp201887v] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Sterical (van der Waals-induced) rotational barriers of methyl groups are investigated theoretically, using ab initio and empirical force field calculations, for various three-dimensional microenvironmental conditions around the methyl group rotator of a model neopentane molecule. The destabilization (reducing methyl rotational barriers) or stabilization (increasing methyl rotational barriers) of the staggered conformation of the methyl rotator depends on a combination of microenvironmental contributions from (i) the number of atoms around the rotator, (ii) the distance between the rotator and the microenvironmental atoms, and (iii) the dihedral angle between the stator, rotator, and molecular environment around the rotator. These geometrical criteria combine their respective effects in a linearly additive fashion, with no apparent cooperative effects, and their combination in space around a rotator may increase, decrease, or leave the rotator's rotational barrier unmodified. This is exemplified in a geometrical analysis of the alanine dipeptide crystal where microenvironmental effects on methyl rotators' barrier of rotation fit the geometrical mapping described in the neopentane model.
Collapse
Affiliation(s)
- William I Hembree
- UT/ORNL Center for Molecular Biophysics, Oak Ridge, Tennessee 37831, United States
| | | |
Collapse
|
13
|
Mauldin RV, Lee AL. Nuclear magnetic resonance study of the role of M42 in the solution dynamics of Escherichia coli dihydrofolate reductase. Biochemistry 2010; 49:1606-15. [PMID: 20073522 DOI: 10.1021/bi901798g] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
It is widely recognized that key positions throughout a protein's structure contribute unequally to function. In light of recent studies that suggest protein dynamics are required for function, a number of these residues may serve to promote motions required for ligand binding and catalysis. In this nuclear magnetic resonance (NMR) study, the conformational dynamics of the dihydrofolate reductase (DHFR) mutant M42W, in the presence of methotrexate and NADPH, are characterized and compared to those of the wild-type enzyme. M42 is distal to the active site, yet the M42W substitution regulates catalysis and ligand affinity and is therefore analogous to an allosteric modulator of DHFR function. To gain understanding of how this mutation regulates activity, we employ a "pandynamic" strategy by measuring conformational fluctuations of backbone amide and side-chain methyl groups on multiple time scales. Changes in pico- to nanosecond dynamics indicate that the mutational effects are propagated throughout a network of interacting residues within DHFR, consistent with a role for M42 as a dynamic communication hub. On the micro- to millisecond time scale, mutation increases the rate of switching in the catalytic core. Mutation also introduces switching in the adenosine binding subdomain that occurs at a higher frequency than in the catalytic core and which correlates with the rate of product release for M42W-DHFR. Finally, a structurally inferred analysis of side-chain dynamics suggests that the M42W mutation dampens motional contributions from nonlocal sources. These data show that the M42W mutation alters the dynamics of DHFR and are consistent with theoretical analysis that suggests this mutation disrupts motion that promotes catalysis.
Collapse
Affiliation(s)
- Randall V Mauldin
- Department of Biochemistry and Biophysics, School of Medicine, University of North Carolina, Chapel Hill, North Carolina 27599-7568, USA
| | | |
Collapse
|