1
|
Zhang J, Kriebel CN, Wan Z, Shi M, Glaubitz C, He X. Automated Fragmentation Quantum Mechanical Calculation of 15N and 13C Chemical Shifts in a Membrane Protein. J Chem Theory Comput 2023; 19:7405-7422. [PMID: 37788419 DOI: 10.1021/acs.jctc.3c00621] [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: 10/05/2023]
Abstract
In this work, we developed an accurate and cost-effective automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) method to calculate the chemical shifts of 15N and 13C of membrane proteins. The convergence of the AF-QM/MM method was tested using Krokinobacter eikastus rhodopsin 2 as a test case. When the distance threshold of the QM region is equal to or larger than 4.0 Å, the results of the AF-QM/MM calculations are close to convergence. In addition, the effects of selected density functionals, basis sets, and local chemical environment of target atoms on the chemical shift calculations were systematically investigated. Our results demonstrate that the predicted chemical shifts are more accurate when important environmental factors including cross-protomer interactions, lipid molecules, and solvent molecules are taken into consideration, especially for the 15N chemical shift prediction. Furthermore, with the presence of sodium ions in the environment, the chemical shift of residues, retinal, and retinal Schiff base are affected, which is consistent with the results of the solid-state nuclear magnetic resonance (NMR) experiment. Upon comparing the performance of various density functionals (namely, B3LYP, B3PW91, M06-2X, M06-L, mPW1PW91, OB95, and OPBE), the results show that mPW1PW91 is a suitable functional for the 15N and 13C chemical shift prediction of the membrane proteins. Meanwhile, we find that the improved accuracy of the 13Cβ chemical shift calculations can be achieved by the employment of the triple-ζ basis set. However, the employment of the triple-ζ basis set does not improve the accuracy of the 15N and 13Cα chemical shift calculations nor does the addition of a diffuse function improve the overall prediction accuracy of the chemical shifts. Our study also underscores that the AF-QM/MM method has significant advantages in predicting the chemical shifts of key ligands and nonstandard residues in membrane proteins than most widely used empirical models; therefore, it could be an accurate computational tool for chemical shift calculations on various types of biological systems.
Collapse
Affiliation(s)
- Jinhuan Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Clara Nassrin Kriebel
- Institute of Biophysical Chemistry and Center for Biomolecular Magnetic Resonance, Goethe University Frankfurt, 60438 Frankfurt am Main, Germany
| | - Zheng Wan
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Man Shi
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
| | - Clemens Glaubitz
- Institute of Biophysical Chemistry and Center for Biomolecular Magnetic Resonance, Goethe University Frankfurt, 60438 Frankfurt am Main, Germany
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, China
| |
Collapse
|
2
|
Kovács D, Bodor A. The influence of random-coil chemical shifts on the assessment of structural propensities in folded proteins and IDPs. RSC Adv 2023; 13:10182-10203. [PMID: 37006359 PMCID: PMC10065145 DOI: 10.1039/d3ra00977g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Accepted: 03/15/2023] [Indexed: 04/03/2023] Open
Abstract
In studying secondary structural propensities of proteins by nuclear magnetic resonance (NMR) spectroscopy, secondary chemical shifts (SCSs) serve as the primary atomic scale observables. For SCS calculation, the selection of an appropriate random coil chemical shift (RCCS) dataset is a crucial step, especially when investigating intrinsically disordered proteins (IDPs). The scientific literature is abundant in such datasets, however, the effect of choosing one over all the others in a concrete application has not yet been studied thoroughly and systematically. Hereby, we review the available RCCS prediction methods and to compare them, we conduct statistical inference by means of the nonparametric sum of ranking differences and comparison of ranks to random numbers (SRD-CRRN) method. We try to find the RCCS predictors best representing the general consensus regarding secondary structural propensities. The existence and the magnitude of resulting differences on secondary structure determination under varying sample conditions (temperature, pH) are demonstrated and discussed for globular proteins and especially IDPs.
Collapse
Affiliation(s)
- Dániel Kovács
- ELTE, Eötvös Loránd University, Institute of Chemistry, Analytical and BioNMR Laboratory Pázmány Péter sétány 1/A Budapest 1117 Hungary
- Eötvös Loránd University, Hevesy György PhD School of Chemistry Pázmány Péter sétány 1/A Budapest 1117 Hungary
| | - Andrea Bodor
- ELTE, Eötvös Loránd University, Institute of Chemistry, Analytical and BioNMR Laboratory Pázmány Péter sétány 1/A Budapest 1117 Hungary
| |
Collapse
|
3
|
Efficient construction of a diverse conformational library for amyloid-β as an intrinsically disordered protein. J Mol Graph Model 2019; 88:183-193. [DOI: 10.1016/j.jmgm.2019.01.009] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2018] [Revised: 01/11/2019] [Accepted: 01/16/2019] [Indexed: 12/16/2022]
|
4
|
Chan-Yao-Chong M, Durand D, Ha-Duong T. Molecular Dynamics Simulations Combined with Nuclear Magnetic Resonance and/or Small-Angle X-ray Scattering Data for Characterizing Intrinsically Disordered Protein Conformational Ensembles. J Chem Inf Model 2019; 59:1743-1758. [PMID: 30840442 DOI: 10.1021/acs.jcim.8b00928] [Citation(s) in RCA: 34] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The concept of intrinsically disordered proteins (IDPs) has emerged relatively slowly, but over the past 20 years, it has become an intense research area in structural biology. Indeed, because of their considerable flexibility and structural heterogeneity, the determination of IDP conformational ensemble is particularly challenging and often requires a combination of experimental measurements and computational approaches. With the improved accuracy of all-atom force fields and the increasing computing performances, molecular dynamics (MD) simulations have become more and more reliable to generate realistic conformational ensembles. And the combination of MD simulations with experimental approaches, such as nuclear magnetic resonance (NMR) and/or small-angle X-ray scattering (SAXS) allows one to converge toward a more accurate and exhaustive description of IDP structures. In this Review, we discuss the state of the art of MD simulations of IDP conformational ensembles, with a special focus on studies that back-calculated and directly compared theoretical and experimental NMR or SAXS observables, such as chemical shifts (CS), 3J-couplings (3Jc), residual dipolar couplings (RDC), or SAXS intensities. We organize the review in three parts. In the first section, we discuss the studies which used NMR and/or SAXS data to test and validate the development of force fields or enhanced sampling techniques. In the second part, we explore different methods for the refinement of MD-derived structural ensembles, such as NMR or SAXS data-restrained MD simulations or ensemble reweighting to better fit experiments. Finally, we survey some recent studies combining MD simulations with NMR and/or SAXS measurements to investigate the relationship between IDP conformational ensemble and biological activity, as well as their implication in human diseases. From this review, we noticed that quite a few studies compared MD-generated conformational ensembles with both NMR and SAXS measurements to validate IDP structures at both local and global levels. Yet, beside the IDP propensity to form local secondary structures, their dynamic extension or compactness also appears important for their activity. Thus, we believe that a close synergy between MD simulations, NMR, and SAXS experiments would be greatly appropriate to address the challenges of characterizing the disordered structures of proteins and their complexes, relative to their biological functions.
Collapse
Affiliation(s)
- Maud Chan-Yao-Chong
- BioCIS, Université Paris-Sud, CNRS , Université Paris-Saclay , 92290 Châtenay-Malabry , France.,Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Univ. Paris-Sud , Université Paris-Saclay , 91198 , Gif-sur-Yvette cedex, France
| | - Dominique Durand
- Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Univ. Paris-Sud , Université Paris-Saclay , 91198 , Gif-sur-Yvette cedex, France
| | - Tâp Ha-Duong
- BioCIS, Université Paris-Sud, CNRS , Université Paris-Saclay , 92290 Châtenay-Malabry , France
| |
Collapse
|
5
|
Vacogne CD, Wei C, Tauer K, Schlaad H. Self-Assembly of α-Helical Polypeptides into Microscopic and Enantiomorphic Spirals. J Am Chem Soc 2018; 140:11387-11394. [DOI: 10.1021/jacs.8b06503] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Charlotte D. Vacogne
- University of Potsdam, Institute of Chemistry, Karl-Liebknecht-Straße 24-25, 14476 Potsdam, Germany
- Max Planck Institute of Colloids and Interfaces, Research Campus Golm, 14424 Potsdam, Germany
| | - Chunxiang Wei
- Max Planck Institute of Colloids and Interfaces, Research Campus Golm, 14424 Potsdam, Germany
| | - Klaus Tauer
- Max Planck Institute of Colloids and Interfaces, Research Campus Golm, 14424 Potsdam, Germany
| | - Helmut Schlaad
- University of Potsdam, Institute of Chemistry, Karl-Liebknecht-Straße 24-25, 14476 Potsdam, Germany
| |
Collapse
|
6
|
Tran L, Basdevant N, Prévost C, Ha-Duong T. Structure of ring-shaped Aβ₄₂ oligomers determined by conformational selection. Sci Rep 2016; 6:21429. [PMID: 26868929 PMCID: PMC4751476 DOI: 10.1038/srep21429] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2015] [Accepted: 01/22/2016] [Indexed: 11/28/2022] Open
Abstract
The oligomerization of amyloid beta (Aβ) peptides into soluble non-fibrillar species plays a critical role in the pathogenesis of Alzheimer’s disease. However, it has been challenging to characterize the tertiary and quaternary structures of Aβ peptides due to their disordered nature and high aggregation propensity. In this work, replica exchange molecular dynamics simulations were used to explore the conformational space of Aβ42 monomer. Among the most populated transient states, we identified a particular conformation which was able to generate ring-shaped pentamers and hexamers, when docked onto itself. The structures of these aggregates were stable during microsecond all-atom MD simulations in explicit solvent. In addition to high resolution models of these oligomers, this study provides support for the conformational selection mechanism of Aβ peptide self-assembly.
Collapse
Affiliation(s)
- Linh Tran
- BIOCIS, Univ. Paris-Sud, CNRS, Université Paris-Saclay, Châtenay-Malabry, 92290, France
| | | | - Chantal Prévost
- LBT, Univ. Paris-Diderot, CNRS, Sorbonne Paris-Cité, Paris, 75005 France
| | - Tâp Ha-Duong
- BIOCIS, Univ. Paris-Sud, CNRS, Université Paris-Saclay, Châtenay-Malabry, 92290, France
| |
Collapse
|
7
|
Robustelli P, Trbovic N, Friesner RA, Palmer AG. Conformational Dynamics of the Partially Disordered Yeast Transcription Factor GCN4. J Chem Theory Comput 2013; 9. [PMID: 24409105 DOI: 10.1021/ct400654r] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
Molecular dynamics (MD) simulations have been employed to study the conformational dynamics of the partially disordered DNA binding basic leucine zipper domain of the yeast transcription factor GCN4. We demonstrate that back-calculated NMR chemical shifts and spin-relaxation data provide complementary probes of the structure and dynamics of disordered protein states and enable comparisons of the accuracy of multiple MD trajectories. In particular, back-calculated chemical shifts provide a sensitive probe of the populations of residual secondary structure elements and helix capping interactions, while spin-relaxation calculations are sensitive to a combination of dynamic and structural factors. Back calculated chemical shift and spin-relaxation data can be used to evaluate the populations of specific interactions in disordered states and identify regions of the phase space that are inconsistent with experimental measurements. The structural interactions that favor and disfavor helical conformations in the disordered basic region of the GCN4 bZip domain were analyzed in order to assess the implications of the structure and dynamics of the apo form for the DNA binding mechanism. The structural couplings observed in these experimentally validated simulations are consistent with a mechanism where the binding of a preformed helical interface would induce folding in the remainder of the protein, supporting a hybrid conformational selection / induced folding binding mechanism.
Collapse
Affiliation(s)
- Paul Robustelli
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York NY 10032, USA
| | - Nikola Trbovic
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York NY 10032, USA
| | | | - Arthur G Palmer
- Department of Biochemistry and Molecular Biophysics, Columbia University, New York NY 10032, USA
| |
Collapse
|
8
|
Chong SH, Ham S. Assessing the influence of solvation models on structural characteristics of intrinsically disordered protein. COMPUT THEOR CHEM 2013. [DOI: 10.1016/j.comptc.2013.05.029] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
9
|
Chong SH, Yim J, Ham S. Structural heterogeneity in familial Alzheimer's disease mutants of amyloid-beta peptides. MOLECULAR BIOSYSTEMS 2013; 9:997-1003. [DOI: 10.1039/c2mb25457c] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
10
|
Lin YS, Bowman GR, Beauchamp KA, Pande VS. Investigating how peptide length and a pathogenic mutation modify the structural ensemble of amyloid beta monomer. Biophys J 2012; 102:315-24. [PMID: 22339868 DOI: 10.1016/j.bpj.2011.12.002] [Citation(s) in RCA: 100] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2011] [Revised: 10/12/2011] [Accepted: 12/01/2011] [Indexed: 12/21/2022] Open
Abstract
The aggregation of amyloid beta (Aβ) peptides plays an important role in the development of Alzheimer's disease. Despite extensive effort, it has been difficult to characterize the secondary and tertiary structure of the Aβ monomer, the starting point for aggregation, due to its hydrophobicity and high aggregation propensity. Here, we employ extensive molecular dynamics simulations with atomistic protein and water models to determine structural ensembles for Aβ(42), Aβ(40), and Aβ(42)-E22K (the Italian mutant) monomers in solution. Sampling of a total of >700 microseconds in all-atom detail with explicit solvent enables us to observe the effects of peptide length and a pathogenic mutation on the disordered Aβ monomer structural ensemble. Aβ(42) and Aβ(40) have crudely similar characteristics but reducing the peptide length from 42 to 40 residues reduces β-hairpin formation near the C-terminus. The pathogenic Italian E22K mutation induces helix formation in the region of residues 20-24. This structural alteration may increase helix-helix interactions between monomers, resulting in altered mechanism and kinetics of Aβ oligomerization.
Collapse
Affiliation(s)
- Yu-Shan Lin
- Department of Chemistry, Stanford University, Stanford, California, USA
| | | | | | | |
Collapse
|
11
|
Olubiyi OO, Strodel B. Structures of the amyloid β-peptides Aβ1-40 and Aβ1-42 as influenced by pH and a D-peptide. J Phys Chem B 2012; 116:3280-91. [PMID: 22300010 DOI: 10.1021/jp2076337] [Citation(s) in RCA: 85] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
In this simulation study, we present a comparison of the secondary structure of the two major alloforms of the Alzheimer's peptide (Aβ(1-40) and Aβ(1-42)) on the basis of molecular dynamics (MD) simulations on thea microsecond time scale using the two GROMOS96 force fields ffG43a2 and ffG53a6. We observe peptide and force-field related differences in the sampled conformations of Aβ(1-40) and Aβ(1-42), which we characterize in terms of NMR chemical shifts calculated from the MD trajectories and validate against the corresponding experimental NMR results. From this analysis, we can conclude that ffG53a6 is better able to model the structural propensities of Aβ(1-40) and Aβ(1-42) than ffG43a2. Furthermore, we provide a description of the influences of pH and binding of D3, a 12-residue D-enantiomeric peptide with demonstrated antiamyloid effects, on the structure of Aβ(1-42). We demonstrate that, under slightly acidic conditions, protonation of the three histidine residues in Aβ(1-42) promotes the formation of β-sheets via a reduction in electrostatic repulsion between the two terminal regions. Our studies further reveal that the binding between D3 and Aβ(1-42) is driven by electrostatic interactions between negatively charged Aβ(1-42) residues and the five positively charged arginine residues of D3. The binding of D3 was found to induce large conformational changes in the amyloid peptide, with a reduction in β-sheet units being the most significant effect recorded, possibly explaining the observed amyloid-inhibiting properties of the D-peptide.
Collapse
Affiliation(s)
- Olujide O Olubiyi
- Institute of Complex Systems: Structural Biochemistry, Research Centre Jülich, Germany
| | | |
Collapse
|
12
|
Zhu T, He X, Zhang JZH. Fragment density functional theory calculation of NMR chemical shifts for proteins with implicit solvation. Phys Chem Chem Phys 2012; 14:7837-45. [PMID: 22314755 DOI: 10.1039/c2cp23746f] [Citation(s) in RCA: 71] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/07/2022]
Abstract
Fragment density functional theory (DFT) calculation of NMR chemical shifts for several proteins (Trp-cage, Pin1 WW domain, the third IgG-binding domain of Protein G (GB3) and human ubiquitin) has been carried out. The present study is based on a recently developed automatic fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach but the solvent effects are included by using the PB (Poisson-Boltzmann) model. Our calculated chemical shifts of (1)H and (13)C for these four proteins are in excellent agreement with experimentally measured values and represent clear improvement over that from the gas phase calculation. However, although the inclusion of the solvent effect also improves the computed chemical shifts of (15)N, the results do not agree with experimental values as well as (1)H and (13)C. Our study also demonstrates that AF-QM/MM calculated results accurately reproduce the separation of α-helical and β-sheet chemical shifts for (13)C(α) atoms in proteins, and using the (1)H chemical shift to discriminate the native structure of proteins from decoys is quite remarkable.
Collapse
Affiliation(s)
- Tong Zhu
- State Key Laboratory of Precision Spectroscopy and Department of Physics, Institute of Theoretical and Computational Science, East China Normal University, Shanghai, China 200062
| | | | | |
Collapse
|
13
|
Chong SH, Park M, Ham S. Structural and Thermodynamic Characteristics That Seed Aggregation of Amyloid-β Protein in Water. J Chem Theory Comput 2012; 8:724-34. [DOI: 10.1021/ct200757a] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023]
Affiliation(s)
- Song-Ho Chong
- Department
of Chemistry, Sookmyung Women’s University,
Hyochangwon-gil 52, Yongsan-gu, Seoul, 140-742, Korea
| | - Mirae Park
- Department
of Chemistry, Sookmyung Women’s University,
Hyochangwon-gil 52, Yongsan-gu, Seoul, 140-742, Korea
| | - Sihyun Ham
- Department
of Chemistry, Sookmyung Women’s University,
Hyochangwon-gil 52, Yongsan-gu, Seoul, 140-742, Korea
| |
Collapse
|