1
|
Pereira AF, Martínez L. Helical Content Correlations and Hydration Structures of the Folding Ensemble of the B Domain of Protein A. J Chem Inf Model 2024; 64:3350-3359. [PMID: 38566451 DOI: 10.1021/acs.jcim.3c01822] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
The B domain of protein A (BdpA), a small three-helix bundle, folds on a time scale of a few microseconds with heterogeneous native and unfolded states. It is widely used as a model for understanding protein folding mechanisms. In this work, we use structure-based models (SBMs) and atomistic simulations to comprehensively investigate how BdpA folding is associated with the formation of its secondary structure. The energy landscape visualization method (ELViM) was used to characterize the pathways that connect the folded and unfolded states of BdpA as well as the sets of structures displaying specific ellipticity patterns. We show that the native state conformational diversity is due mainly to the conformational variability of helix I. Helices I, II, and III occur in a weakly correlated manner, with Spearman's rank correlation coefficients of 0.1539 (I and II), 0.1259 (I and III), and 0.2561 (II and III). These results, therefore, suggest the highest cooperativity between helices II and III. Our results allow the clustering of partially folded structures of folding of the B domain of protein A on the basis of its secondary structure, paving the way to an understanding of environmental factors in the relative stability of the basins of the folding ensemble, which are illustrated by the structural dependency of the protein hydration structures, as computed with minimum-distance distribution functions.
Collapse
Affiliation(s)
- Ander Francisco Pereira
- Institute of Chemistry and Center for Computing in Engineering & Science, Universidade Estadual de Campinas (UNICAMP), 13083-861 Campinas, SP, Brazil
| | - Leandro Martínez
- Institute of Chemistry and Center for Computing in Engineering & Science, Universidade Estadual de Campinas (UNICAMP), 13083-861 Campinas, SP, Brazil
| |
Collapse
|
2
|
Yanaka S, Yagi-Utsumi M, Kato K, Kuwajima K. The B domain of protein A retains residual structures in 6 M guanidinium chloride as revealed by hydrogen/deuterium-exchange NMR spectroscopy. Protein Sci 2023; 32:e4569. [PMID: 36659853 PMCID: PMC9926473 DOI: 10.1002/pro.4569] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2022] [Revised: 01/10/2023] [Accepted: 01/13/2023] [Indexed: 01/21/2023]
Abstract
The characterization of residual structures persistent in unfolded proteins is an important issue in studies of protein folding, because the residual structures present, if any, may form a folding initiation site and guide the subsequent folding reactions. Here, we studied the residual structures of the isolated B domain (BDPA) of staphylococcal protein A in 6 M guanidinium chloride. BDPA is a small three-helix-bundle protein, and until recently its folding/unfolding reaction has been treated as a simple two-state process between the native and the fully unfolded states. We employed a dimethylsulfoxide (DMSO)-quenched hydrogen/deuterium (H/D)-exchange 2D NMR techniques with the use of spin desalting columns, which allowed us to investigate the H/D-exchange behavior of individually identified peptide amide (NH) protons. We obtained H/D-exchange protection factors of the 21 NH protons that form an α-helical hydrogen bond in the native structure, and the majority of these NH protons were significantly protected with a protection factor of 2.0-5.2 in 6 M guanidinium chloride, strongly suggesting that these weakly protected NH protons form much stronger hydrogen bonds under native folding conditions. The results can be used to deduce the structure of an early folding intermediate, when such an intermediate is shown by other methods. Among three native helical regions, the third helix in the C-terminal side was highly protected and stabilized by side-chain salt bridges, probably acting as the folding initiation site of BDPA. The present results are discussed in relation to previous experimental and computational findings on the folding mechanisms of BDPA.
Collapse
Affiliation(s)
- Saeko Yanaka
- Exploratory Research Center on Life and Living Systems (ExCELLS) and Institute for Molecular Science (IMS), National Institutes of Natural Sciences, Myodaiji, Okazaki, Aichi, Japan.,Department of Functional Molecular Science, School of Physical Sciences, SOKENDAI (the Graduate University for Advanced Studies), Myodaiji, Okazaki, Aichi, Japan.,Graduate School of Pharmaceutical Sciences, Nagoya City University, Mizuho-ku, Nagoya, Aichi, Japan
| | - Maho Yagi-Utsumi
- Exploratory Research Center on Life and Living Systems (ExCELLS) and Institute for Molecular Science (IMS), National Institutes of Natural Sciences, Myodaiji, Okazaki, Aichi, Japan.,Department of Functional Molecular Science, School of Physical Sciences, SOKENDAI (the Graduate University for Advanced Studies), Myodaiji, Okazaki, Aichi, Japan.,Graduate School of Pharmaceutical Sciences, Nagoya City University, Mizuho-ku, Nagoya, Aichi, Japan
| | - Koichi Kato
- Exploratory Research Center on Life and Living Systems (ExCELLS) and Institute for Molecular Science (IMS), National Institutes of Natural Sciences, Myodaiji, Okazaki, Aichi, Japan.,Department of Functional Molecular Science, School of Physical Sciences, SOKENDAI (the Graduate University for Advanced Studies), Myodaiji, Okazaki, Aichi, Japan.,Graduate School of Pharmaceutical Sciences, Nagoya City University, Mizuho-ku, Nagoya, Aichi, Japan
| | - Kunihiro Kuwajima
- Department of Physics, School of Science, University of Tokyo, Bunkyo-ku, Tokyo, Japan
| |
Collapse
|
3
|
Leavens MJ, Spang LE, Cherney MM, Bowler BE. Denatured State Conformational Biases in Three-Helix Bundles Containing Divergent Sequences Localize near Turns and Helix Capping Residues. Biochemistry 2021; 60:3071-3085. [PMID: 34606713 PMCID: PMC8751257 DOI: 10.1021/acs.biochem.1c00400] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Rhodopseudomonas palustris cytochrome c', a four-helix bundle, and the second ubiquitin-associated domain, UBA(2), a three-helix bundle from the human homologue of yeast Rad23, HHR23A, deviate from random coil behavior under denaturing conditions in a fold-specific manner. The random coil deviations in each of these folds occur near interhelical turns and loops in their tertiary structures. Here, we examine an additional three-helix bundle with an identical fold to UBA(2), but a highly divergent sequence, the first ubiquitin-associated domain, UBA(1), of HHR23A. We use histidine-heme loop formation methods, employing eight single histidine variants, to probe for denatured state conformational bias of a UBA(1) domain fused to the N-terminus of iso-1-cytochrome c (iso-1-Cytc). Guanidine hydrochloride (GuHCl) denaturation shows that the iso-1-Cytc domain unfolds first, followed by the UBA(1) domain. Denatured state (4 and 6 M GuHCl) histidine-heme loop formation studies show that as the size of the histidine-heme loop increases, loop stability decreases, as expected for the Jacobson-Stockmayer relationship. However, loops formed with His35, His31, and His15, of UBA(1), are 0.6-1.1 kcal/mol more stable than expected from the Jacobson-Stockmayer relationship, confirming the importance of deviations of the denatured state from random coil behavior near interhelical turns of helical domains for facilitating folding to the correct topology. For UBA(1) and UBA(2), hydrophobic clusters on either side of the turns partially explain deviations from random coil behavior; however, helix capping also appears to be important.
Collapse
Affiliation(s)
- Moses J. Leavens
- Department of Chemistry & Biochemistry, University of Montana, Missoula, Montana 59812, United States
- Center for Biomolecular Structure & Dynamics, University of Montana, Missoula, Montana 59812, United States
| | - Lisa E. Spang
- Department of Chemistry & Biochemistry, University of Montana, Missoula, Montana 59812, United States
- Center for Biomolecular Structure & Dynamics, University of Montana, Missoula, Montana 59812, United States
| | - Melisa M. Cherney
- Department of Chemistry & Biochemistry, University of Montana, Missoula, Montana 59812, United States
- Center for Biomolecular Structure & Dynamics, University of Montana, Missoula, Montana 59812, United States
| | - Bruce E. Bowler
- Department of Chemistry & Biochemistry, University of Montana, Missoula, Montana 59812, United States
- Center for Biomolecular Structure & Dynamics, University of Montana, Missoula, Montana 59812, United States
| |
Collapse
|
4
|
Mazigi O, Schofield P, Langley DB, Christ D. Protein A superantigen: structure, engineering and molecular basis of antibody recognition. Protein Eng Des Sel 2019; 32:359-366. [PMID: 31641749 DOI: 10.1093/protein/gzz026] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2019] [Revised: 07/02/2019] [Accepted: 07/11/2019] [Indexed: 12/22/2022] Open
Abstract
Staphylococcus aureus interacts with the human immune system through the production of secreted factors. Key among these is protein A, a B-cell superantigen capable of interacting with both antibody Fc and VH regions. Here, we review structural and molecular features of this important example of naturally occurring bacterial superantigens, as well as engineered variants and their application in biotechnology.
Collapse
Affiliation(s)
- Ohan Mazigi
- Department of Immunology, Garvan Institute of Medical Research, 384 Victoria Street, Darlinghurst, Sydney, NSW 2010, Australia.,Faculty of Medicine, St Vincent's Clinical School, UNSW Sydney, Darlinghurst, Sydney, NSW 2010, Australia
| | - Peter Schofield
- Department of Immunology, Garvan Institute of Medical Research, 384 Victoria Street, Darlinghurst, Sydney, NSW 2010, Australia.,Faculty of Medicine, St Vincent's Clinical School, UNSW Sydney, Darlinghurst, Sydney, NSW 2010, Australia
| | - David B Langley
- Department of Immunology, Garvan Institute of Medical Research, 384 Victoria Street, Darlinghurst, Sydney, NSW 2010, Australia
| | - Daniel Christ
- Department of Immunology, Garvan Institute of Medical Research, 384 Victoria Street, Darlinghurst, Sydney, NSW 2010, Australia.,Faculty of Medicine, St Vincent's Clinical School, UNSW Sydney, Darlinghurst, Sydney, NSW 2010, Australia
| |
Collapse
|
5
|
Martin OA, Vorobjev Y, Scheraga HA, Vila JA. Outline of an experimental design aimed to detect a protein A mirror image in solution. PEERJ PHYSICAL CHEMISTRY 2019; 1. [PMID: 34079958 DOI: 10.7717/peerj-pchem.2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
There is abundant theoretical evidence indicating that a mirror image of Protein A may occur during the protein folding process. However, as to whether such mirror image exists in solution is an unsolved issue. Here we provide outline of an experimental design aimed to detect the mirror image of Protein A in solution. The proposal is based on computational simulations indicating that the use of a mutant of protein A, namely Q10H, could be used to detect the mirror image conformation in solution. Our results indicate that the native conformation of the protein A should have a pKa, for the Q10H mutant, at ≈6.2, while the mirror-image conformation should have a pKa close to ≈7.3. Naturally, if all the population is in the native state for the Q10H mutant, the pKa should be ≈6.2, while, if all are in the mirror-image state, it would be ≈7.3, and, if it is a mixture, the pKa should be largerthan 6.2, presumably in proportion to the mirror population. In addition, evidence is provided indicating the tautomeric distribution of H10 must also change between the native and mirror conformations. Although this may not be completely relevant for the purpose of determining whether the protein A mirror image exists in solution, it could provide valuable information to validate the pKa findings. We hope this proposal will foster experimental work on this problem either by direct application of our proposed experimental design or serving as inspiration and motivation for other experiments.
Collapse
Affiliation(s)
- Osvaldo A Martin
- Instituto de Matemática Aplicada San Luis, UNSL-CONICET, San Luis, Argentina
| | - Yury Vorobjev
- Institute of Chemical Biology and Fundamental Medicine, Siberian Branch of the Russian Academy of Science, Novosibirsk, Russia
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, United States
| | - Jorge A Vila
- Instituto de Matemática Aplicada San Luis, UNSL-CONICET, San Luis, Argentina.,Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, United States
| |
Collapse
|
6
|
Otosu T, Yamaguchi S. Two-Dimensional Fluorescence Lifetime Correlation Spectroscopy: Concepts and Applications. Molecules 2018; 23:E2972. [PMID: 30441830 PMCID: PMC6278346 DOI: 10.3390/molecules23112972] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2018] [Revised: 11/07/2018] [Accepted: 11/13/2018] [Indexed: 01/03/2023] Open
Abstract
We review the basic concepts and recent applications of two-dimensional fluorescence lifetime correlation spectroscopy (2D FLCS), which is the extension of fluorescence correlation spectroscopy (FCS) to analyze the correlation of fluorescence lifetime in addition to fluorescence intensity. Fluorescence lifetime is sensitive to the microenvironment and can be a "molecular ruler" when combined with FRET. Utilization of fluorescence lifetime in 2D FLCS thus enables us to quantify the inhomogeneity of the system and the interconversion dynamics among different species with a higher time resolution than other single-molecule techniques. Recent applications of 2D FLCS to various biological systems demonstrate that 2D FLCS is a unique and promising tool to quantitatively analyze the microsecond conformational dynamics of macromolecules at the single-molecule level.
Collapse
Affiliation(s)
- Takuhiro Otosu
- Department of Applied Chemistry, Graduate School of Science and Engineering, Saitama University, 255 Shimo-Okubo, Sakura, Saitama 338-8570, Japan.
| | - Shoichi Yamaguchi
- Department of Applied Chemistry, Graduate School of Science and Engineering, Saitama University, 255 Shimo-Okubo, Sakura, Saitama 338-8570, Japan.
| |
Collapse
|
7
|
Shao Q, Zhu W. The effects of implicit modeling of nonpolar solvation on protein folding simulations. Phys Chem Chem Phys 2018; 20:18410-18419. [PMID: 29946610 DOI: 10.1039/c8cp03156h] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Implicit solvent models, in which the polar and nonpolar solvation free-energies of solute molecules are treated separately, have been widely adopted for molecular dynamics simulation of protein folding. While the development of the implicit models is mainly focused on the methodological improvement and key parameter optimization for polar solvation, nonpolar solvation has been either ignored or described by a simplistic surface area (SA) model. In this work, we performed the folding simulations of multiple β-hairpin and α-helical proteins with varied surface tension coefficients embedded in the SA model to clearly demonstrate the effects of nonpolar solvation treated by a popular SA model on protein folding. The results indicate that the change in the surface tension coefficient does not alter the ability of implicit solvent simulations to reproduce a protein native structure but indeed controls the components of the equilibrium conformational ensemble and modifies the energetic characterization of the folding transition pathway. The suitably set surface tension coefficient can yield explicit solvent simulations and/or experimentally suggested folding mechanism of protein. In addition, the implicit treatment of both polar and nonpolar components of solvation free-energy contributes to the overestimation of the secondary structure in implicit solvent simulations.
Collapse
Affiliation(s)
- Qiang Shao
- Drug Discovery and Design Center, CAS Key Laboratory of Receptor Research, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203, China.
| | | |
Collapse
|
8
|
Shao Q, Zhu W. How Well Can Implicit Solvent Simulations Explore Folding Pathways? A Quantitative Analysis of α-Helix Bundle Proteins. J Chem Theory Comput 2017; 13:6177-6190. [DOI: 10.1021/acs.jctc.7b00726] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Qiang Shao
- Drug
Discovery and Design Center, CAS Key Laboratory of Receptor Research,
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
- University of
Chinese Academy of Sciences, Beijing 100049, China
| | - Weiliang Zhu
- Drug
Discovery and Design Center, CAS Key Laboratory of Receptor Research,
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai 201203, China
- University of
Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
9
|
Rahamim G, Amir D, Haas E. Simultaneous Determination of Two Subdomain Folding Rates Using the "Transfer-Quench" Method. Biophys J 2017; 112:1786-1796. [PMID: 28494950 DOI: 10.1016/j.bpj.2017.01.037] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2016] [Revised: 12/21/2016] [Accepted: 01/06/2017] [Indexed: 11/29/2022] Open
Abstract
The investigation of the mechanism of protein folding is complicated by the context dependence of the rates of intramolecular contact formation. Methods based on site-specific labeling and ultrafast spectroscopic detection of fluorescence signals were developed for monitoring the rates of individual subdomain folding transitions in situ, in the context of the whole molecule. However, each site-specific labeling modification might affect rates of folding of near-neighbor structural elements, and thus limit the ability to resolve fine differences in rates of folding of these elements. Therefore, it is highly desirable to be able to study the rates of folding of two or more neighboring subdomain structures using a single mutant to facilitate resolution of the order and interdependence of such steps. Here, we report the development of the "Transfer-Quench" method for measuring the rate of formation of two structural elements using a single triple-labeled mutant. This method is based on Förster resonance energy transfer combined with fluorescence quenching. We placed the donor and acceptor at the loop ends, and a quencher at an α-helical element involved in the node forming the loop. The folding of the triple-labeled mutant is monitored by the acceptor emission. The formation of nonlocal contact (loop closure) increases the time-dependent acceptor emission, while the closure of the labeled helix turn reduces this emission. The method was applied in a study of the folding mechanism of the common model protein, the B domain of staphylococcal protein A. Only natural amino acids were used as probes, and thus possible structural perturbations were minimized. Tyr and Trp residues served as donor and acceptor at the ends of a long loop between helices I and II, and a Cys residue as a quencher for the acceptor. We found that the closure of the loop (segment 14-33) occurs with the same rate constant as the nucleation of helix HII (segment 33-29), in line with the nucleation-condensation model.
Collapse
Affiliation(s)
- Gil Rahamim
- The Goodman Faculty of Life Sciences Bar Ilan University, Ramat Gan, Israel
| | - Dan Amir
- The Goodman Faculty of Life Sciences Bar Ilan University, Ramat Gan, Israel
| | - Elisha Haas
- The Goodman Faculty of Life Sciences Bar Ilan University, Ramat Gan, Israel.
| |
Collapse
|
10
|
Otosu T, Ishii K, Oikawa H, Arai M, Takahashi S, Tahara T. Highly Heterogeneous Nature of the Native and Unfolded States of the B Domain of Protein A Revealed by Two-Dimensional Fluorescence Lifetime Correlation Spectroscopy. J Phys Chem B 2017; 121:5463-5473. [DOI: 10.1021/acs.jpcb.7b00546] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
| | - Kunihiko Ishii
- Ultrafast Spectroscopy Research Team, RIKEN Center for Advanced Photonics (RAP), 2-1, Hirosawa, Wako, Saitama 351-0198, Japan
| | - Hiroyuki Oikawa
- Institute
of Multidisciplinary Research for Advanced Materials, Tohoku University, 2-1-1
Katahira, Aoba, Sendai, Miyagi 980-8577, Japan
| | - Munehito Arai
- Department
of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan
| | - Satoshi Takahashi
- Institute
of Multidisciplinary Research for Advanced Materials, Tohoku University, 2-1-1
Katahira, Aoba, Sendai, Miyagi 980-8577, Japan
| | - Tahei Tahara
- Ultrafast Spectroscopy Research Team, RIKEN Center for Advanced Photonics (RAP), 2-1, Hirosawa, Wako, Saitama 351-0198, Japan
| |
Collapse
|
11
|
Harada R, Takano Y, Shigeta Y. Common folding processes of mini-proteins: Partial formations of secondary structures initiate the immediate protein folding. J Comput Chem 2017; 38:790-797. [DOI: 10.1002/jcc.24748] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2016] [Revised: 12/20/2016] [Accepted: 01/12/2017] [Indexed: 11/10/2022]
Affiliation(s)
- Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba; 1-1-1 Tennodai Tsukuba Ibaraki 305-8577 Japan
| | - Yu Takano
- Department of Biomedical Information Sciences; Graduate School of Information Sciences, Hiroshima City University; 3-4-1 Ozuka-Higashi, Asa-Minami-Ku Hiroshima 731-3194 Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba; 1-1-1 Tennodai Tsukuba Ibaraki 305-8577 Japan
| |
Collapse
|
12
|
Pancsa R, Raimondi D, Cilia E, Vranken WF. Early Folding Events, Local Interactions, and Conservation of Protein Backbone Rigidity. Biophys J 2017; 110:572-583. [PMID: 26840723 DOI: 10.1016/j.bpj.2015.12.028] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2015] [Revised: 12/21/2015] [Accepted: 12/29/2015] [Indexed: 01/20/2023] Open
Abstract
Protein folding is in its early stages largely determined by the protein sequence and complex local interactions between amino acids, resulting in lower energy conformations that provide the context for further folding into the native state. We compiled a comprehensive data set of early folding residues based on pulsed labeling hydrogen deuterium exchange experiments. These early folding residues have corresponding higher backbone rigidity as predicted by DynaMine from sequence, an effect also present when accounting for the secondary structures in the folded protein. We then show that the amino acids involved in early folding events are not more conserved than others, but rather, early folding fragments and the secondary structure elements they are part of show a clear trend toward conserving a rigid backbone. We therefore propose that backbone rigidity is a fundamental physical feature conserved by proteins that can provide important insights into their folding mechanisms and stability.
Collapse
Affiliation(s)
- Rita Pancsa
- Structural Biology Brussels, Vrije Universiteit Brussel, Brussels, Belgium
| | - Daniele Raimondi
- Structural Biology Brussels, Vrije Universiteit Brussel, Brussels, Belgium
| | - Elisa Cilia
- Structural Biology Brussels, Vrije Universiteit Brussel, Brussels, Belgium
| | - Wim F Vranken
- Structural Biology Brussels, Vrije Universiteit Brussel, Brussels, Belgium.
| |
Collapse
|
13
|
Dissecting protein architecture with communication blocks and communicating segment pairs. BMC Bioinformatics 2016; 17 Suppl 2:13. [PMID: 26823083 PMCID: PMC4959365 DOI: 10.1186/s12859-015-0855-y] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Proteins adapt to environmental conditions by changing their shape and motions. Characterising protein conformational dynamics is increasingly recognised as necessary to understand how proteins function. Given a conformational ensemble, computational tools are needed to extract in a systematic way pertinent and comprehensive biological information. RESULTS Here, we present a method, Communication Mapping (COMMA), to decipher the dynamical architecture of a protein. The method first extracts residue-based dynamic properties from all-atom molecular dynamics simulations. Then, it integrates them in a graph theoretic framework, where it identifies groups of residues or protein regions that mediate short- and long-range communication. COMMA introduces original concepts to contrast the different roles played by these regions, namely communication blocks and communicating segment pairs, and evaluates the connections and communication strengths between them. We show the utility and capabilities of COMMA by applying it to three archetypal proteins, namely protein A, the tyrosine kinase KIT and the tumour suppressor p53. CONCLUSION Our method permits to compare in a direct way the dynamical behaviour either of proteins with different characteristics or of the same protein in different conditions. It is useful to identify residues playing a key role in protein allosteric regulation and to explain the effects of deleterious mutations in a mechanistic way. COMMA is a fully automated tool with broad applicability. It is freely available to the community at www.lcqb.upmc.fr/COMMA .
Collapse
|
14
|
Oikawa H, Kamagata K, Arai M, Takahashi S. Complexity of the folding transition of the B domain of protein A revealed by the high-speed tracking of single-molecule fluorescence time series. J Phys Chem B 2015; 119:6081-91. [PMID: 25938341 DOI: 10.1021/acs.jpcb.5b00414] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The equilibrium unfolding transition of the B domain of protein A (BdpA) was investigated by using single-molecule fluorescence spectroscopy based on line-confocal detection of fast-flowing samples. The method achieved the time resolution of 120 μs and the observation time of a few milliseconds in the single-molecule time-series measurements of fluorescence resonance energy transfer (FRET). Two samples of BdpA doubly labeled with donor and acceptor fluorophores, the first possessing fluorophores at residues 22 and 55 (sample 1) and the second at residues 5 and 55 (sample 2), were prepared. The equilibrium unfolding transition induced by guanidium chloride (GdmCl) was monitored by bulk measurements and demonstrated that the both samples obey the apparent two-state unfolding. In the absence of GdmCl, the single-molecule FRET measurements for the both samples showed a single peak assignable to the native state (N). The FRET efficiency for N shifts to lower values as the increase of GdmCl concentration, suggesting the swelling of the native state structure. At the higher concentration of GdmCl, the both samples convert to the unfolded state (U). Near the unfolding midpoint for sample 1, the kinetic exchange between N and U causes the averaging of the two states and the higher values of the relative fluctuation. The time series for different molecules in U showed slightly different FRET efficiencies, suggesting the apparent heterogeneity. Thus, the high-speed tracking of fluorescence signals from single molecules revealed a complexity and heterogeneity hidden in the apparent two-state behavior of protein folding.
Collapse
Affiliation(s)
- Hiroyuki Oikawa
- †Institute of Multidisciplinary Research for Advanced Materials, Tohoku University, 2-1-1 Katahira, Aoba, Sendai, Miyagi 980-8577, Japan
| | - Kiyoto Kamagata
- †Institute of Multidisciplinary Research for Advanced Materials, Tohoku University, 2-1-1 Katahira, Aoba, Sendai, Miyagi 980-8577, Japan
| | - Munehito Arai
- ‡Department of Life Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan.,§PRESTO, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Satoshi Takahashi
- †Institute of Multidisciplinary Research for Advanced Materials, Tohoku University, 2-1-1 Katahira, Aoba, Sendai, Miyagi 980-8577, Japan
| |
Collapse
|
15
|
Rimratchada S, McLeish TCB, Radford SE, Paci E. The role of high-dimensional diffusive search, stabilization, and frustration in protein folding. Biophys J 2014; 106:1729-40. [PMID: 24739172 DOI: 10.1016/j.bpj.2014.01.051] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2013] [Revised: 01/17/2014] [Accepted: 01/31/2014] [Indexed: 11/25/2022] Open
Abstract
Proteins are polymeric molecules with many degrees of conformational freedom whose internal energetic interactions are typically screened to small distances. Therefore, in the high-dimensional conformation space of a protein, the energy landscape is locally relatively flat, in contrast to low-dimensional representations, where, because of the induced entropic contribution to the full free energy, it appears funnel-like. Proteins explore the conformation space by searching these flat subspaces to find a narrow energetic alley that we call a hypergutter and then explore the next, lower-dimensional, subspace. Such a framework provides an effective representation of the energy landscape and folding kinetics that does justice to the essential characteristic of high-dimensionality of the search-space. It also illuminates the important role of nonnative interactions in defining folding pathways. This principle is here illustrated using a coarse-grained model of a family of three-helix bundle proteins whose conformations, once secondary structure has formed, can be defined by six rotational degrees of freedom. Two folding mechanisms are possible, one of which involves an intermediate. The stabilization of intermediate subspaces (or states in low-dimensional projection) in protein folding can either speed up or slow down the folding rate depending on the amount of native and nonnative contacts made in those subspaces. The folding rate increases due to reduced-dimension pathways arising from the mere presence of intermediate states, but decreases if the contacts in the intermediate are very stable and introduce sizeable topological or energetic frustration that needs to be overcome. Remarkably, the hypergutter framework, although depending on just a few physically meaningful parameters, can reproduce all the types of experimentally observed curvature in chevron plots for realizations of this fold.
Collapse
Affiliation(s)
| | - Tom C B McLeish
- Department of Physics and Biophysical Sciences Institute, Durham University, Durham, United Kingdom.
| | - Sheena E Radford
- Astbury Centre for Structural Molecular Biology, University of Leeds, Leeds, United Kingdom
| | - Emanuele Paci
- School of Physics and Astronomy, University of Leeds, Leeds, United Kingdom; Department of Physics and Biophysical Sciences Institute, Durham University, Durham, United Kingdom.
| |
Collapse
|
16
|
Accounting for a mirror-image conformation as a subtle effect in protein folding. Proc Natl Acad Sci U S A 2014; 111:8458-63. [PMID: 24912167 DOI: 10.1073/pnas.1407837111] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
By using local (free-energy profiles along the amino acid sequence and (13)C(α) chemical shifts) and global (principal component) analyses to examine the molecular dynamics of protein-folding trajectories, generated with the coarse-grained united-residue force field, for the B domain of staphylococcal protein A, we are able to (i) provide the main reason for formation of the mirror-image conformation of this protein, namely, a slow formation of the second loop and part of the third helix (Asp29-Asn35), caused by the presence of multiple local conformational states in this portion of the protein; (ii) show that formation of the mirror-image topology is a subtle effect resulting from local interactions; (iii) provide a mechanism for how protein A overcomes the barrier between the metastable mirror-image state and the native state; and (iv) offer a plausible reason to explain why protein A does not remain in the metastable mirror-image state even though the mirror-image and native conformations are at least energetically compatible.
Collapse
|
17
|
Krokhotin A, Liwo A, Maisuradze GG, Niemi AJ, Scheraga HA. Kinks, loops, and protein folding, with protein A as an example. J Chem Phys 2014; 140:025101. [PMID: 24437917 PMCID: PMC3899063 DOI: 10.1063/1.4855735] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2013] [Accepted: 12/10/2013] [Indexed: 02/02/2023] Open
Abstract
The dynamics and energetics of formation of loops in the 46-residue N-terminal fragment of the B-domain of staphylococcal protein A has been studied. Numerical simulations have been performed using coarse-grained molecular dynamics with the united-residue (UNRES) force field. The results have been analyzed in terms of a kink (heteroclinic standing wave solution) of a generalized discrete nonlinear Schrödinger (DNLS) equation. In the case of proteins, the DNLS equation arises from a C(α)-trace-based energy function. Three individual kink profiles were identified in the experimental three-α-helix structure of protein A, in the range of the Glu16-Asn29, Leu20-Asn29, and Gln33-Asn44 residues, respectively; these correspond to two loops in the native structure. UNRES simulations were started from the full right-handed α-helix to obtain a clear picture of kink formation, which would otherwise be blurred by helix formation. All three kinks emerged during coarse-grained simulations. It was found that the formation of each is accompanied by a local free energy increase; this is expressed as the change of UNRES energy which has the physical sense of the potential of mean force of a polypeptide chain. The increase is about 7 kcal/mol. This value can thus be considered as the free energy barrier to kink formation in full α-helical segments of polypeptide chains. During the simulations, the kinks emerge, disappear, propagate, and annihilate each other many times. It was found that the formation of a kink is initiated by an abrupt change in the orientation of a pair of consecutive side chains in the loop region. This resembles the formation of a Bloch wall along a spin chain, where the C(α) backbone corresponds to the chain, and the amino acid side chains are interpreted as the spin variables. This observation suggests that nearest-neighbor side chain-side chain interactions are responsible for initiation of loop formation. It was also found that the individual kinks are reflected as clear peaks in the principal modes of the analyzed trajectory of protein A, the shapes of which resemble the directional derivatives of the kinks along the chain. These observations suggest that the kinks of the DNLS equation determine the functionally important motions of proteins.
Collapse
Affiliation(s)
- Andrey Krokhotin
- Department of Physics and Astronomy and Science for Life Laboratory, Uppsala University, P.O. Box 803, S-75108 Uppsala, Sweden
| | - Adam Liwo
- Faculty of Chemistry, University of Gdansk, ul. Sobieskiego 18, 80-952 Gdansk, Poland
| | - Gia G Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, USA
| | - Antti J Niemi
- Department of Physics and Astronomy and Science for Life Laboratory, Uppsala University, P.O. Box 803, S-75108 Uppsala, Sweden
| | - Harold A Scheraga
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, USA
| |
Collapse
|
18
|
Davis CM, Dyer RB. Dynamics of an ultrafast folding subdomain in the context of a larger protein fold. J Am Chem Soc 2013; 135:19260-7. [PMID: 24320936 PMCID: PMC3949483 DOI: 10.1021/ja409608r] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Small fast folding subdomains with low contact order have been postulated to facilitate the folding of larger proteins. We have tested this idea by determining how the fastest folding linear β-hairpin, CLN025, which folds on the nanosecond time scale, folds within the context of a two-hairpin WW domain system, which folds on the microsecond time scale. The folding of the wild type FBP28 WW domain was compared to constructs in which each of the loops was replaced by CLN025. A combination of FTIR spectroscopy and laser-induced temperature-jump coupled with infrared spectroscopy was used to probe changes in the peptide backbone. The relaxation dynamics of the β-sheets and β-turn were measured independently by probing the corresponding bands assigned in the amide I region. The folding rate of the CLN025 β-hairpin is unchanged within the larger protein. Insertion of the β-hairpin into the second loop results in an overall stabilization of the WW domain and a relaxation lifetime five times faster than the parent WW domain. In both mutants, folding is initiated in the turns and the β-sheets form last. These results demonstrate that fast folding subdomains can be used to speed the folding of more complex proteins, and that the folding dynamics of the subdomain is unchanged within the context of the larger protein.
Collapse
Affiliation(s)
- Caitlin M. Davis
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| | - R. Brian Dyer
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
19
|
Glyakina AV, Pereyaslavets LB, Galzitskaya OV. Right- and left-handed three-helix proteins. I. Experimental and simulation analysis of differences in folding and structure. Proteins 2013; 81:1527-41. [DOI: 10.1002/prot.24301] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2012] [Revised: 03/27/2013] [Accepted: 03/28/2013] [Indexed: 11/08/2022]
Affiliation(s)
- Anna V. Glyakina
- Institute of Protein Research; Russian Academy of Sciences; Pushchino, Moscow Region 142290 Russia
- Institute of Mathematical Problems of Biology; Russian Academy of Sciences; Pushchino, Moscow Region 142290 Russia
| | - Leonid B. Pereyaslavets
- Institute of Protein Research; Russian Academy of Sciences; Pushchino, Moscow Region 142290 Russia
| | - Oxana V. Galzitskaya
- Institute of Protein Research; Russian Academy of Sciences; Pushchino, Moscow Region 142290 Russia
| |
Collapse
|
20
|
A Hamiltonian replica exchange molecular dynamics (MD) method for the study of folding, based on the analysis of the stabilization determinants of proteins. Int J Mol Sci 2013; 14:12157-69. [PMID: 23743827 PMCID: PMC3709779 DOI: 10.3390/ijms140612157] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2013] [Revised: 04/29/2013] [Accepted: 05/14/2013] [Indexed: 11/16/2022] Open
Abstract
Herein, we present a novel Hamiltonian replica exchange protocol for classical molecular dynamics simulations of protein folding/unfolding. The scheme starts from the analysis of the energy-networks responsible for the stabilization of the folded conformation, by means of the energy-decomposition approach. In this framework, the compact energetic map of the native state is generated by a preliminary short molecular dynamics (MD) simulation of the protein in explicit solvent. This map is simplified by means of an eigenvalue decomposition. The highest components of the eigenvector associated with the lowest eigenvalue indicate which sites, named "hot spots", are likely to be responsible for the stability and correct folding of the protein. In the Hamiltonian replica exchange protocol, we use modified force-field parameters to treat the interparticle non-bonded potentials of the hot spots within the protein and between protein and solvent atoms, leaving unperturbed those relative to all other residues, as well as solvent-solvent interactions. We show that it is possible to reversibly simulate the folding/unfolding behavior of two test proteins, namely Villin HeadPiece HP35 (35 residues) and Protein A (62 residues), using a limited number of replicas. We next discuss possible implications for the study of folding mechanisms via all atom simulations.
Collapse
|
21
|
Lobanov MY, Suvorina MY, Dovidchenko NV, Sokolovskiy IV, Surin AK, Galzitskaya OV. A novel web server predicts amino acid residue protection against hydrogen–deuterium exchange. Bioinformatics 2013; 29:1375-81. [DOI: 10.1093/bioinformatics/btt168] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
22
|
Shao Q, Zhu W, Gao YQ. Robustness in Protein Folding Revealed by Thermodynamics Calculations. J Phys Chem B 2012; 116:13848-56. [DOI: 10.1021/jp307684h] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Qiang Shao
- Institute of Theoretical and
Computational Chemistry, College of Chemistry and Molecular Engineering,
Beijing National Laboratory of Molecular Sciences, Peking University, Beijing 100871, China
- Drug Discovery and Design Center,
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203,
China
| | - Weiliang Zhu
- Drug Discovery and Design Center,
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203,
China
| | - Yi Qin Gao
- Institute of Theoretical and
Computational Chemistry, College of Chemistry and Molecular Engineering,
Beijing National Laboratory of Molecular Sciences, Peking University, Beijing 100871, China
| |
Collapse
|
23
|
Dib L, Carbone A. Protein fragments: functional and structural roles of their coevolution networks. PLoS One 2012; 7:e48124. [PMID: 23139761 PMCID: PMC3489791 DOI: 10.1371/journal.pone.0048124] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2012] [Accepted: 09/27/2012] [Indexed: 11/19/2022] Open
Abstract
Small protein fragments, and not just residues, can be used as basic building blocks to reconstruct networks of coevolved amino acids in proteins. Fragments often enter in physical contact one with the other and play a major biological role in the protein. The nature of these interactions might be multiple and spans beyond binding specificity, allosteric regulation and folding constraints. Indeed, coevolving fragments are indicators of important information explaining folding intermediates, peptide assembly, key mutations with known roles in genetic diseases, distinguished subfamily-dependent motifs and differentiated evolutionary pressures on protein regions. Coevolution analysis detects networks of fragments interaction and highlights a high order organization of fragments demonstrating the importance of studying at a deeper level this structure. We demonstrate that it can be applied to protein families that are highly conserved or represented by few sequences, enlarging in this manner, the class of proteins where coevolution analysis can be performed and making large-scale coevolution studies a feasible goal.
Collapse
Affiliation(s)
- Linda Dib
- Université Pierre et Marie Curie, UMR 7238, Équipe de Génomique Analytique, Paris, France
- CNRS, UMR 7238, Laboratoire de Génomique des Microorganismes, Paris, France
| | - Alessandra Carbone
- Université Pierre et Marie Curie, UMR 7238, Équipe de Génomique Analytique, Paris, France
- CNRS, UMR 7238, Laboratoire de Génomique des Microorganismes, Paris, France
| |
Collapse
|
24
|
Shao Q, Fan Y, Yang L, Gao YQ. Counterion Effects on the Denaturing Activity of Guanidinium Cation to Protein. J Chem Theory Comput 2012; 8:4364-73. [DOI: 10.1021/ct3002267] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Qiang Shao
- Institute of Theoretical and
Computational Chemistry, College of Chemistry and Molecular Engineering,
Beijing National Laboratory of Molecular Sciences, Peking University, Beijing 100871, China
- Drug Discovery and Design Center,
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, 555 Zuchongzhi Road, Shanghai, 201203,
China
| | - Yubo Fan
- Bioinformatics and
Bioengineering
Program, The Methodist Hospital Research Institute, Weill Cornell Medical College, Houston, Texas 77030, United
States
| | - Lijiang Yang
- Institute of Theoretical and
Computational Chemistry, College of Chemistry and Molecular Engineering,
Beijing National Laboratory of Molecular Sciences, Peking University, Beijing 100871, China
| | - Yi Qin Gao
- Institute of Theoretical and
Computational Chemistry, College of Chemistry and Molecular Engineering,
Beijing National Laboratory of Molecular Sciences, Peking University, Beijing 100871, China
| |
Collapse
|
25
|
Yin Y, Maisuradze GG, Liwo A, Scheraga HA. Hidden protein folding pathways in free-energy landscapes uncovered by network analysis. J Chem Theory Comput 2012; 8:1176-1189. [PMID: 22715321 DOI: 10.1021/ct200806n] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
A network analysis is used to uncover hidden folding pathways in free-energy landscapes usually defined in terms of such arbitrary order parameters as root-mean-square deviation from the native structure, radius of gyration, etc. The analysis has been applied to molecular dynamics (MD) trajectories of the B-domain of staphylococcal protein A, generated with the coarse-grained united-residue (UNRES) force field in a broad range of temperatures (270K ≤ T ≤ 325K). Thousands of folding pathways have been identified at each temperature. Out of these many folding pathways, several most probable ones were selected for investigation of the conformational transitions during protein folding. Unlike other conformational space network (CSN) methods, a node in the CSN variant implemented in this work is defined according to the nativelikeness class of the structure, which defines the similarity of segments of the compared structures in terms of secondary-structure, contact-pattern, and local geometry, as well as the overall geometric similarity of the conformation under consideration to that of the reference (experimental) structure. Our previous findings, regarding the folding model and conformations found at the folding-transition temperature for protein A (Maisuradze et al., J. Am. Chem. Soc. 132, 9444, 2010), were confirmed by the conformational space network analysis. In the methodology and in the analysis of the results, the shortest path identified by using the shortest-path algorithm corresponds to the most probable folding pathway in the conformational space network.
Collapse
Affiliation(s)
- Yanping Yin
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York, 14850-1301
| | | | | | | |
Collapse
|
26
|
Shao Q, Qin Gao Y. The relative helix and hydrogen bond stability in the B domain of protein A as revealed by integrated tempering sampling molecular dynamics simulation. J Chem Phys 2011; 135:135102. [DOI: 10.1063/1.3630127] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
|
27
|
Kmiecik S, Kolinski A. Simulation of chaperonin effect on protein folding: a shift from nucleation-condensation to framework mechanism. J Am Chem Soc 2011; 133:10283-9. [PMID: 21618995 PMCID: PMC3132998 DOI: 10.1021/ja203275f] [Citation(s) in RCA: 37] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The iterative annealing mechanism (IAM) of chaperonin-assisted protein folding is explored in a framework of a well-established coarse-grained protein modeling tool, which enables the study of protein dynamics in a time-scale well beyond classical all-atom molecular mechanics. The chaperonin mechanism of action is simulated for two paradigm systems of protein folding, B domain of protein A (BdpA) and B1 domain of protein G (GB1), and compared to chaperonin-free simulations presented here for BdpA and recently published for GB1. The prediction of the BdpA transition state ensemble (TSE) is in perfect agreement with experimental findings. It is shown that periodic distortion of the polypeptide chains by hydrophobic chaperonin interactions can promote rapid folding and leads to a decrease in folding temperature. It is also demonstrated how chaperonin action prevents kinetically trapped conformations and modulates the observed folding mechanisms from nucleation-condensation to a more framework-like.
Collapse
Affiliation(s)
- Sebastian Kmiecik
- Faculty of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland.
| | | |
Collapse
|
28
|
Maisuradze GG, Liwo A, Ołdziej S, Scheraga HA. Evidence, from simulations, of a single state with residual native structure at the thermal denaturation midpoint of a small globular protein. J Am Chem Soc 2010; 132:9444-52. [PMID: 20568747 DOI: 10.1021/ja1031503] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
The folding of the B-domain of staphylococcal protein A has been studied by coarse-grained canonical and multiplexed replica-exchange molecular dynamics simulations with the UNRES force field in a broad range of temperatures (270 K < or = T < or = 350 K). In canonical simulations, the folding was found to occur either directly to the native state or through kinetic traps, mainly the topological mirror image of the native three-helix bundle. The latter folding scenario was observed more frequently at low temperatures. With increase of temperature, the frequency of the transitions between the folded and misfolded/unfolded states increased and the folded state became more diffuse with conformations exhibiting increased root-mean-square deviations from the experimental structure (from about 4 A at T = 300 K to 8.7 A at T = 325 K). An analysis of the equilibrium conformational ensemble determined from multiplexed replica exchange simulations at the folding-transition temperature (T(f) = 325 K) showed that the conformational ensemble at this temperature is a collection of conformations with residual secondary structures, which possess native or near-native clusters of nonpolar residues in place, and not a 50-50% mixture of fully folded and fully unfolded conformations. These findings contradict the quasi-chemical picture of two- or multistate protein folding, which assumes an equilibrium between the folded, unfolded, and intermediate states, with equilibrium shifting with temperature but with the native conformations remaining essentially unchanged. Our results also suggest that long-range hydrophobic contacts are the essential factor to keep the structure of a protein thermally stable.
Collapse
Affiliation(s)
- Gia G Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, USA
| | | | | | | |
Collapse
|
29
|
Structural basis for high-affinity HER2 receptor binding by an engineered protein. Proc Natl Acad Sci U S A 2010; 107:15039-44. [PMID: 20696930 DOI: 10.1073/pnas.1005025107] [Citation(s) in RCA: 154] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
The human epidermal growth factor receptor 2 (HER2) is specifically overexpressed in tumors of several cancers, including an aggressive form of breast cancer. It is therefore a target for both cancer diagnostics and therapy. The 58 amino acid residue Zher2 affibody molecule was previously engineered as a high-affinity binder of HER2. Here we determined the structure of Zher2 in solution and the crystal structure of Zher2 in complex with the HER2 extracellular domain. Zher2 binds to a conformational epitope on HER2 that is distant from those recognized by the therapeutic antibodies trastuzumab and pertuzumab. Its small size and lack of interference may provide Zher2 with advantages for diagnostic use or even for delivery of therapeutic agents to HER2-expressing tumors when trastuzumab or pertuzumab are already employed. Biophysical characterization shows that Zher2 is thermodynamically stable in the folded state yet undergoing conformational interconversion on a submillisecond time scale. The data suggest that it is the HER2-binding conformation that is formed transiently prior to binding. Still, binding is very strong with a dissociation constant K(D) = 22 pM, and perfect conformational homogeneity is therefore not necessarily required in engineered binding proteins. A comparison of the original Z domain scaffold to free and bound Zher2 structures reveals how high-affinity binding has evolved during selection and affinity maturation and suggests how a compromise between binding surface optimization and stability and dynamics of the unbound state has been reached.
Collapse
|
30
|
Maisuradze GG, Senet P, Czaplewski C, Liwo A, Scheraga HA. Investigation of protein folding by coarse-grained molecular dynamics with the UNRES force field. J Phys Chem A 2010; 114:4471-85. [PMID: 20166738 PMCID: PMC2849147 DOI: 10.1021/jp9117776] [Citation(s) in RCA: 82] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Coarse-grained molecular dynamics simulations offer a dramatic extension of the time-scale of simulations compared to all-atom approaches. In this article, we describe the use of the physics-based united-residue (UNRES) force field, developed in our laboratory, in protein-structure simulations. We demonstrate that this force field offers about a 4000-times extension of the simulation time scale; this feature arises both from averaging out the fast-moving degrees of freedom and reduction of the cost of energy and force calculations compared to all-atom approaches with explicit solvent. With massively parallel computers, microsecond folding simulation times of proteins containing about 1000 residues can be obtained in days. A straightforward application of canonical UNRES/MD simulations, demonstrated with the example of the N-terminal part of the B-domain of staphylococcal protein A (PDB code: 1BDD, a three-alpha-helix bundle), discerns the folding mechanism and determines kinetic parameters by parallel simulations of several hundred or more trajectories. Use of generalized-ensemble techniques, of which the multiplexed replica exchange method proved to be the most effective, enables us to compute thermodynamics of folding and carry out fully physics-based prediction of protein structure, in which the predicted structure is determined as a mean over the most populated ensemble below the folding-transition temperature. By using principal component analysis of the UNRES folding trajectories of the formin-binding protein WW domain (PDB code: 1E0L; a three-stranded antiparallel beta-sheet) and 1BDD, we identified representative structures along the folding pathways and demonstrated that only a few (low-indexed) principal components can capture the main structural features of a protein-folding trajectory; the potentials of mean force calculated along these essential modes exhibit multiple minima, as opposed to those along the remaining modes that are unimodal. In addition, a comparison between the structures that are representative of the minima in the free-energy profile along the essential collective coordinates of protein folding (computed by principal component analysis) and the free-energy profile projected along the virtual-bond dihedral angles gamma of the backbone revealed the key residues involved in the transitions between the different basins of the folding free-energy profile, in agreement with existing experimental data for 1E0L .
Collapse
Affiliation(s)
- Gia G Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, USA
| | | | | | | | | |
Collapse
|
31
|
Feldwisch J, Tolmachev V, Lendel C, Herne N, Sjöberg A, Larsson B, Rosik D, Lindqvist E, Fant G, Höidén-Guthenberg I, Galli J, Jonasson P, Abrahmsén L. Design of an optimized scaffold for affibody molecules. J Mol Biol 2010; 398:232-47. [PMID: 20226194 DOI: 10.1016/j.jmb.2010.03.002] [Citation(s) in RCA: 118] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2009] [Revised: 02/26/2010] [Accepted: 03/03/2010] [Indexed: 12/15/2022]
Abstract
Affibody molecules are non-immunoglobulin-derived affinity proteins based on a three-helical bundle protein domain. Here, we describe the design process of an optimized Affibody molecule scaffold with improved properties and a surface distinctly different from that of the parental scaffold. The improvement was achieved by applying an iterative process of amino acid substitutions in the context of the human epidermal growth factor receptor 2 (HER2)-specific Affibody molecule Z(HER2:342). Replacements in the N-terminal region, loop 1, helix 2 and helix 3 were guided by extensive structural modeling using the available structures of the parent Z domain and Affibody molecules. The effect of several single substitutions was analyzed followed by combination of up to 11 different substitutions. The two amino acid substitutions N23T and S33K accounted for the most dramatic improvements, including increased thermal stability with elevated melting temperatures of up to +12 degrees C. The optimized scaffold contains 11 amino acid substitutions in the nonbinding surface and is characterized by improved thermal and chemical stability, as well as increased hydrophilicity, and enables generation of identical Affibody molecules both by chemical peptide synthesis and by recombinant bacterial expression. A HER2-specific Affibody tracer, [MMA-DOTA-Cys61]-Z(HER2:2891)-Cys (ABY-025), was produced by conjugating MMA-DOTA (maleimide-monoamide-1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid) to the peptide produced either chemically or in Escherichia coli. ABY-025 showed high affinity and specificity for HER2 (equilibrium dissociation constant, K(D), of 76 pM) and detected HER2 in tissue sections of SKOV-3 xenograft and human breast tumors. The HER2-binding capacity was fully retained after three cycles of heating to 90 degrees C followed by cooling to room temperature. Furthermore, the binding surfaces of five Affibody molecules targeting other proteins (tumor necrosis factor alpha, insulin, Taq polymerase, epidermal growth factor receptor or platelet-derived growth factor receptor beta) were grafted onto the optimized scaffold, resulting in molecules with improved thermal stability and a more hydrophilic nonbinding surface.
Collapse
|
32
|
Maisuradze GG, Liwo A, Scheraga HA. Relation between free energy landscapes of proteins and dynamics. J Chem Theory Comput 2010; 6:583-595. [PMID: 23620713 DOI: 10.1021/ct9005745] [Citation(s) in RCA: 113] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
By examining the molecular dynamics (MD) of protein folding trajectories, generated with the coarse-grained UNRES force field, for the B-domain of staphylococcal protein A and the triple β-strand WW domain from the Formin binding protein 28 (FBP), by principal component analysis (PCA), it is demonstrated how different free energy landscapes (FELs) and folding pathways of trajectories can be, even though they appear to be very similar by visual inspection of the time-dependence of the root-mean-square deviation (rmsd). Approaches to determine the minimal dimensionality of FELs for a correct description of protein folding dynamics are discussed. The correlation between the amplitude of the fluctuations of proteins and the dimensionality of the FELs is shown. The advantage of internal coordinate PCA over Cartesian PCA for small proteins is also illustrated.
Collapse
Affiliation(s)
- Gia G Maisuradze
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University Ithaca, New York 14853-1301
| | | | | |
Collapse
|
33
|
Localized thermodynamic coupling between hydrogen bonding and microenvironment polarity substantially stabilizes proteins. Nat Struct Mol Biol 2009; 16:684-90. [PMID: 19525973 PMCID: PMC2754385 DOI: 10.1038/nsmb.1610] [Citation(s) in RCA: 155] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2009] [Accepted: 04/23/2009] [Indexed: 11/08/2022]
Abstract
The energetic contributions of hydrogen bonding to protein folding are still unclear, despite more than 70 years of study. This is due partly to the difficulty of extracting thermodynamic information about specific interactions from protein mutagenesis data and partly to the context dependence of hydrogen bond strengths. Herein, we test the hypothesis that hydrogen bond strengths depend on the polarity of their microenvironment, with stronger hydrogen bonds forming in nonpolar surroundings. Double-mutant cycle analysis using a combination of amide-to-ester backbone mutagenesis and traditional side chain mutagenesis revealed that hydrogen bonds can be stronger by up to 1.2 kcal mol(-1) when they are sequestered in hydrophobic surroundings than when they are solvent exposed. Such large coupling energies between hydrogen bond strengths and local polarity suggest that the context dependence of hydrogen bond strengths must be accounted for in any comprehensive account of the forces responsible for protein folding.
Collapse
|
34
|
Ulutaş B, Haliloglu T, Bozma I. Folding pathways explored with artificial potential functions. Phys Biol 2009; 6:036008. [DOI: 10.1088/1478-3975/6/3/036008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
35
|
Baxa MC, Freed KF, Sosnick TR. Quantifying the structural requirements of the folding transition state of protein A and other systems. J Mol Biol 2008; 381:1362-81. [PMID: 18625237 PMCID: PMC2742318 DOI: 10.1016/j.jmb.2008.06.067] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2008] [Revised: 06/18/2008] [Accepted: 06/23/2008] [Indexed: 11/22/2022]
Abstract
The B-domain of protein A is a small three-helix bundle that has been the subject of considerable experimental and theoretical investigation. Nevertheless, a unified view of the structure of the transition-state ensemble (TSE) is still lacking. To characterize the TSE of this surprisingly challenging protein, we apply a combination of psi analysis (which probes the role of specific side-chain to side-chain contacts) and kinetic H/D amide isotope effects (which measures hydrogen-bond content), building upon previous studies using mutational phi analysis (which probes the energetic influence of side-chain substitutions). The second helix is folded in the TSE, while helix formation appears just at the carboxy and amino termini of the first and third helices, respectively. The experimental data suggest a homogenous yet plastic TS with a native-like topology. This study generalizes our earlier conclusion, based on two larger alpha/beta proteins, that the TSEs of most small proteins achieve approximately 70% of their native state's relative contact order. This high percentage limits the degree of possible TS heterogeneity and requires a reevaluation of the structural content of the TSE of other proteins, especially when they are characterized as small or polarized.
Collapse
Affiliation(s)
- Michael C. Baxa
- Department of Physics, University of Chicago, 929 E. 57th St., Chicago, IL 60637
- Institute for Biophysical Dynamics, University of Chicago, 929 E. 57th St., Chicago, IL 60637
| | - Karl F. Freed
- James Franck Institute and Department of Chemistry, University of Chicago, 929 E. 57th St., Chicago, IL 60637
| | - Tobin R. Sosnick
- Institute for Biophysical Dynamics, University of Chicago, 929 E. 57th St., Chicago, IL 60637
- Department of Biochemistry and Molecular Biology, University of Chicago, 929 E. 57th St., Chicago, IL 60637
| |
Collapse
|
36
|
Lei H, Wu C, Wang ZX, Zhou Y, Duan Y. Folding processes of the B domain of protein A to the native state observed in all-atom ab initio folding simulations. J Chem Phys 2008; 128:235105. [PMID: 18570534 DOI: 10.1063/1.2937135] [Citation(s) in RCA: 32] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
Reaching the native states of small proteins, a necessary step towards a comprehensive understanding of the folding mechanisms, has remained a tremendous challenge to ab initio protein folding simulations despite the extensive effort. In this work, the folding process of the B domain of protein A (BdpA) has been simulated by both conventional and replica exchange molecular dynamics using AMBER FF03 all-atom force field. Started from an extended chain, a total of 40 conventional (each to 1.0 micros) and two sets of replica exchange (each to 200.0 ns per replica) molecular dynamics simulations were performed with different generalized-Born solvation models and temperature control schemes. The improvements in both the force field and solvent model allowed successful simulations of the folding process to the native state as demonstrated by the 0.80 A C(alpha) root mean square deviation (RMSD) of the best folded structure. The most populated conformation was the native folded structure with a high population. This was a significant improvement over the 2.8 A C(alpha) RMSD of the best nativelike structures from previous ab initio folding studies on BdpA. To the best of our knowledge, our results demonstrate, for the first time, that ab initio simulations can reach the native state of BdpA. Consistent with experimental observations, including Phi-value analyses, formation of helix II/III hairpin was a crucial step that provides a template upon which helix I could form and the folding process could complete. Early formation of helix III was observed which is consistent with the experimental results of higher residual helical content of isolated helix III among the three helices. The calculated temperature-dependent profile and the melting temperature were in close agreement with the experimental results. The simulations further revealed that phenylalanine 31 may play critical to achieve the correct packing of the three helices which is consistent with the experimental observation. In addition to the mechanistic studies, an ab initio structure prediction was also conducted based on both the physical energy and a statistical potential. Based on the lowest physical energy, the predicted structure was 2.0 A C(alpha) RMSD away from the experimentally determined structure.
Collapse
Affiliation(s)
- Hongxing Lei
- UC Davis Genome Center and Department of Applied Science, University of California at Davis, One Shields Avenue, Davis, California 95616, USA
| | | | | | | | | |
Collapse
|
37
|
Morra G, Colombo G. Relationship between energy distribution and fold stability: Insights from molecular dynamics simulations of native and mutant proteins. Proteins 2008; 72:660-72. [PMID: 18247351 DOI: 10.1002/prot.21963] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
Most proteins must fold to a well-defined structure with a minimal stability to perform their function. Here we use a simple, molecular dynamics-based, energy decomposition approach to map the principal energetic interactions in a set of proteins representative of different folds. This work involves the all-atom simulation and analysis of the native structures and mutants of five different proteins representative of an all-alpha (yACPB, Protein A), all-beta (SH3), and a mixed alpha/beta fold (Proteins G and L). Given a certain structure, a native sequence and a set of mutants, we show that our model discriminates the ability of a mutation to yield a more or less stable protein, in agreement with experimental data, catching the principal energetic determinants of protein stabilization. Our approach identifies the interaction determinants responsible to define a fold and shows that mutations can either modulate the strength of pair-wise coupling between residues important for folding, or modify the profile of the principal interactions. Furthermore, we address the question of how to evaluate the fitness of a sequence to a given structure by comparing the information contained in the energy map, which recapitulates the chemistry of the sequence, to that contained in the contact map, which recapitulates the fold topology. The results show that the better fit between the energetic properties of the sequence and the fold topology corresponds to a higher stabilization of the protein. We discuss the relevance of these observations to the analysis of protein designability and to the rational evolution of new sequences.
Collapse
Affiliation(s)
- Giulia Morra
- Istituto di Chimica del Riconoscimento Molecolare, CNR, Via Mario Bianco 9, 20131, Milano, Italy
| | | |
Collapse
|
38
|
St-Pierre JF, Mousseau N, Derreumaux P. The complex folding pathways of protein A suggest a multiple-funnelled energy landscape. J Chem Phys 2008; 128:045101. [PMID: 18248008 DOI: 10.1063/1.2812562] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
Folding proteins into their native states requires the formation of both secondary and tertiary structures. Many questions remain, however, as to whether these form into a precise order, and various pictures have been proposed that place the emphasis on the first or the second level of structure in describing folding. One of the favorite test models for studying this question is the B domain of protein A, which has been characterized by numerous experiments and simulations. Using the activation-relaxation technique coupled with a generic energy model (optimized potential for efficient peptide structure prediction), we generate more than 50 folding trajectories for this 60-residue protein. While the folding pathways to the native state are fully consistent with the funnel-like description of the free energy landscape, we find a wide range of mechanisms in which secondary and tertiary structures form in various orders. Our nonbiased simulations also reveal the presence of a significant number of non-native beta and alpha conformations both on and off pathway, including the visit, for a non-negligible fraction of trajectories, of fully ordered structures resembling the native state of nonhomologous proteins.
Collapse
Affiliation(s)
- Jean-Francois St-Pierre
- Département de Physique, Université de Montréal, C.P. 6128, Succursale Centre-Ville, Montréal, Québec H3C 3J7, Canada
| | | | | |
Collapse
|
39
|
Universality and diversity of folding mechanics for three-helix bundle proteins. Proc Natl Acad Sci U S A 2008; 105:895-900. [PMID: 18195374 DOI: 10.1073/pnas.0707284105] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
In this study we evaluate, at full atomic detail, the folding processes of two small helical proteins, the B domain of protein A and the Villin headpiece. Folding kinetics are studied by performing a large number of ab initio Monte Carlo folding simulations using a single transferable all-atom potential. Using these trajectories, we examine the relaxation behavior, secondary structure formation, and transition-state ensembles (TSEs) of the two proteins and compare our results with experimental data and previous computational studies. To obtain a detailed structural information on the folding dynamics viewed as an ensemble process, we perform a clustering analysis procedure based on graph theory. Moreover, rigorous p(fold) analysis is used to obtain representative samples of the TSEs and a good quantitative agreement between experimental and simulated Phi values is obtained for protein A. Phi values for Villin also are obtained and left as predictions to be tested by future experiments. Our analysis shows that the two-helix hairpin is a common partially stable structural motif that gets formed before entering the TSE in the studied proteins. These results together with our earlier study of Engrailed Homeodomain and recent experimental studies provide a comprehensive, atomic-level picture of folding mechanics of three-helix bundle proteins.
Collapse
|
40
|
Gu Z, Rao MK, Forsyth WR, Finke JM, Matthews CR. Structural analysis of kinetic folding intermediates for a TIM barrel protein, indole-3-glycerol phosphate synthase, by hydrogen exchange mass spectrometry and Gō model simulation. J Mol Biol 2007; 374:528-46. [PMID: 17942114 DOI: 10.1016/j.jmb.2007.09.024] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2007] [Revised: 09/06/2007] [Accepted: 09/07/2007] [Indexed: 10/22/2022]
Abstract
The structures of partially folded states appearing during the folding of a (betaalpha)(8) TIM barrel protein, the indole-3-glycerol phosphate synthase from Sulfolobus solfataricus (sIGPS), was assessed by hydrogen exchange mass spectrometry (HX-MS) and Gō model simulations. HX-MS analysis of the peptic peptides derived from the pulse-labeled product of the sub-millisecond folding reaction from the urea-denatured state revealed strong protection in the (betaalpha)(4) region, modest protection in the neighboring (betaalpha)(1-3) and (betaalpha)(5)beta(6) segments and no significant protection in the remaining N and C-terminal segments. These results demonstrate that this species is not a collapsed form of the unfolded state under native-favoring conditions nor is it the native state formed via fast-track folding. However, the striking contrast of these results with the strong protection observed in the (betaalpha)(2-5)beta(6) region after 5 s of folding demonstrates that these species represent kinetically distinct folding intermediates that are not identical as previously thought. A re-examination of the kinetic folding mechanism by chevron analysis of fluorescence data confirmed distinct roles for these two species: the burst-phase intermediate is predicted to be a misfolded, off-pathway intermediate, while the subsequent 5 s intermediate corresponds to an on-pathway equilibrium intermediate. Comparison with the predictions using a C(alpha) Gō model simulation of the kinetic folding reaction for sIGPS shows good agreement with the core of the structure offering protection against exchange in the on-pathway intermediate(s). Because the native-centric Gō model simulations do not explicitly include sequence-specific information, the simulation results support the hypothesis that the topology of TIM barrel proteins is a primary determinant of the folding free energy surface for the productive folding reaction. The early misfolding reaction must involve aspects of non-native structure not detected by the Gō model simulation.
Collapse
Affiliation(s)
- Zhenyu Gu
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts Medical School, Worcester, MA 01605, USA
| | | | | | | | | |
Collapse
|
41
|
Sato S, Fersht AR. Searching for Multiple Folding Pathways of a Nearly Symmetrical Protein: Temperature Dependent Φ-Value Analysis of the B Domain of Protein A. J Mol Biol 2007; 372:254-67. [PMID: 17628591 DOI: 10.1016/j.jmb.2007.06.043] [Citation(s) in RCA: 32] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2007] [Revised: 05/25/2007] [Accepted: 06/15/2007] [Indexed: 11/19/2022]
Abstract
The B domain of protein A (BdpA) is a popular paradigm for simulating protein folding pathways. The discrepancies between so many simulations and subsequent experimental testing may be attributable to the protein being highly symmetrical: changing experimental conditions could perturb the subtle interplay between the effects of symmetry in the native structure and the effects of asymmetry from specific interactions in a given sequence. If the protein folds via multiple pathways, perturbations, such as temperature, denaturant concentration, and mutation, should change the flux of micro pathways, leading to changes in the bulk properties of the transition state. We tested this hypothesis by conducting a Phi-analysis of BdpA as a function of temperature from 25.0 degrees C to 60.0 degrees C. The Phi-values had no significant dependence on temperature and the values at 55.0 degrees C (denaturing conditions) are very similar to those at 25.0 degrees C (folding conditions), indicating the structure of the transition state does not significantly change although the experimental conditions are considerably altered. The results suggest that BdpA folds via a single dominant folding pathway.
Collapse
Affiliation(s)
- Satoshi Sato
- MRC Centre for Protein Engineering, Hills Road, Cambridge CB2 2QH, UK
| | | |
Collapse
|
42
|
Ozkan SB, Wu GA, Chodera JD, Dill KA. Protein folding by zipping and assembly. Proc Natl Acad Sci U S A 2007; 104:11987-92. [PMID: 17620603 PMCID: PMC1924571 DOI: 10.1073/pnas.0703700104] [Citation(s) in RCA: 119] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2006] [Indexed: 11/18/2022] Open
Abstract
How do proteins fold so quickly? Some denatured proteins fold to their native structures in only microseconds, on average, implying that there is a folding "mechanism," i.e., a particular set of events by which the protein short-circuits a broader conformational search. Predicting protein structures using atomically detailed physical models is currently challenging. The most definitive proof of a putative folding mechanism would be whether it speeds up protein structure prediction in physical models. In the zipping and assembly (ZA) mechanism, local structuring happens first at independent sites along the chain, then those structures either grow (zip) or coalescence (assemble) with other structures. Here, we apply the ZA search mechanism to protein native structure prediction by using the AMBER96 force field with a generalized Born/surface area implicit solvent model and sampling by replica exchange molecular dynamics. Starting from open denatured conformations, our algorithm, called the ZA method, converges to an average of 2.2 A from the Protein Data Bank native structures of eight of nine proteins that we tested, which ranged from 25 to 73 aa in length. In addition, experimental Phi values, where available on these proteins, are consistent with the predicted routes. We conclude that ZA is a viable model for how proteins physically fold. The present work also shows that physics-based force fields are quite good and that physics-based protein structure prediction may be practical, at least for some small proteins.
Collapse
Affiliation(s)
| | | | - John D. Chodera
- Graduate Group in Biophysics, University of California, San Francisco, CA 94143
| | | |
Collapse
|
43
|
Chowdhury S, Lei H, Duan Y. Denatured-state ensemble and the early-stage folding of the G29A mutant of the B-domain of protein A. J Phys Chem B 2007; 109:9073-81. [PMID: 16852080 DOI: 10.1021/jp0449814] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The folding mechanism of the G29A mutant of the B-domain of protein A (BdpA) has been studied by all-atom molecular dynamics simulation using AMBER force field (ff03) and generalized Born continuum solvent model. Started from the extended chain conformation, a total of 16 simulations (400 ns each) at 300 K captured some early folding events of the G29A mutant of BdpA. In one of the 16 trajectories, the G29A mutant folded within 2.8 A (root mean square) of the wild-type NMR structure. We observed that the fast burial of hydrophobic residues was the driving force to bring the distant residues into close proximity. The initiation of the helix I and III occurred during the stage of hydrophobic collapse. The initiation and growth of the helix II was slow. Both the secondary structure formation and the development of the native tertiary contacts suggested a multistage folding process. Clustering analysis indicated that two helix species (helices I and III) could be intermediates. Further analysis revealed that the hydrophobic residues of partially folded helix II formed nativelike hydrophobic contacts with helices I and III that stabilized a nativelike state and delayed the completion of folding of the entire protein. The details of the early folding process were compared with other theoretical and experimental studies. It was found that a nativelike hydrophobic cluster was formed by residues including F(30), I(31), L(34), L(44), L(45), and A(48) that prevented further development of the native structures, and breaking the hydrophobic cluster like this one contributed to the rate-limiting step. This was in complete agreement with the recent kinetic measurements in which mutations of these residues to Gly and Ala substantially increased the folding rates by as much as 60 times. Apparently, destabilization of nonnative states dramatically enhanced the folding rates.
Collapse
Affiliation(s)
- Shibasish Chowdhury
- University of California-Davis Genome Center and Bioinformatics Program, Department of Applied Science, University of California, One Shields Avenue, Davis, California 95616, USA
| | | | | |
Collapse
|
44
|
Dyer RB. Ultrafast and downhill protein folding. Curr Opin Struct Biol 2007; 17:38-47. [PMID: 17223539 DOI: 10.1016/j.sbi.2007.01.001] [Citation(s) in RCA: 56] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2006] [Revised: 12/05/2006] [Accepted: 01/03/2007] [Indexed: 11/24/2022]
Abstract
Ultrafast folding proteins have served an important role in benchmarking molecular dynamics simulations and testing protein folding theories. These proteins are simple enough and fold fast enough that realistic simulations are possible, which facilitates the direct comparison of absolute folding rates and folding mechanisms with those observed experimentally. Such comparisons have achieved remarkable success, but have also revealed the shortcomings that remain in experiment, theory and simulation alike. Some ultrafast folding proteins may fold without encountering an activation barrier (downhill folding), allowing the exploration of the molecular timescale of folding and the roughness of the energy landscape. The biological significance of ultrafast folding remains uncertain, but its practical significance is crucial to progress in understanding how proteins fold.
Collapse
Affiliation(s)
- R Brian Dyer
- Chemistry Division, MS J567, Los Alamos National Laboratory, Los Alamos, NM 87545, USA.
| |
Collapse
|
45
|
Jagielska A, Scheraga HA. Influence of temperature, friction, and random forces on folding of the B-domain of staphylococcal protein A: All-atom molecular dynamics in implicit solvent. J Comput Chem 2007; 28:1068-82. [PMID: 17279497 DOI: 10.1002/jcc.20631] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
The influences of temperature, friction, and random forces on the folding of protein A have been analyzed. A series of all-atom molecular dynamics folding simulations with the Amber ff99 potential and Generalized Born solvation, starting from the fully extended chain, were carried out for temperatures from 300 to 500 K, using (a) the Berendsen thermostat (with no explicit friction or random forces) and (b) Langevin dynamics (with friction and stochastic forces explicitly present in the system). The simulation temperature influences the relative time scale of the major events on the folding pathways of protein A. At lower temperatures, helix 2 folds significantly later than helices 1 and 3. However, with increasing temperature, the folding time of helix 2 approaches the folding times of helices 1 and 3. At lower temperatures, the complete formation of secondary and tertiary structure is significantly separated in time whereas, at higher temperatures, they occur simultaneously. These results suggest that some earlier experimental and theoretical observations of folding events, e.g., the order of helix formation, could depend on the temperature used in those studies. Therefore, the differences in temperature used could be one of the reasons for the discrepancies among published experimental and computational studies of the folding of protein A. Friction and random forces do not change the folding pathway that was observed in the simulations with the Berendsen thermostat, but their explicit presence in the system extends the folding time of protein A.
Collapse
Affiliation(s)
- Anna Jagielska
- Baker Laboratory of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853-1301, USA
| | | |
Collapse
|
46
|
Kameda T, Takada S. Secondary structure provides a template for the folding of nearby polypeptides. Proc Natl Acad Sci U S A 2006; 103:17765-70. [PMID: 17101976 PMCID: PMC1693821 DOI: 10.1073/pnas.0602632103] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Although protein structures are primarily encoded by their sequences, they are also critically dependent on environmental factors such as solvents and interactions with other molecules. Here we investigate how the folding-energy landscape of a short peptide is altered by interactions with another peptide, by performing atomistic replica-exchange molecular dynamics simulations of polyalanines in various environments. We analyzed the free-energy landscapes of Ala7 and Ala8 in isolation, near an alpha-helix template, and near a beta-strand template. The isolated Ala7 and Ala8 at 270 K were mainly in polyproline II helix conformations and in equilibrium between the alpha-helix and polyproline II helix, respectively, in harmony with the experiment. Interestingly, we found remarkably strong secondary-structure "templating"; namely, the alpha-helix template enhanced alpha-helix conformation and the beta-strand template induced beta-strand conformation in the simulated Ala8. The alpha-helix template lowered the nearby dielectric constant, which strengthened hydrogen bonds in the simulated Ala8, leading to alpha-helix stabilization. The beta-strand template provided hydrogen bond positions to the simulated Ala8, sharply inducing beta-strand structure. With or without templates, the energy landscape of Ala8 is always funnel-like and centered at the alpha-helix conformation, whereas entropic contribution disfavors the alpha-helix, leading to subtle competition. Secondary-structure templating may play a critical role in protein conformation dynamics in the cellular environment.
Collapse
Affiliation(s)
- Tomoshi Kameda
- *Graduate School of Science and Technology, Kobe University, Rokkodai, Nada, Kobe 657-8501, Japan
- Computational Biology Research Center, Advanced Industrial Science and Technology, 2-43 Aomi, Koto, Tokyo 135-0064, Japan; and
| | - Shoji Takada
- *Graduate School of Science and Technology, Kobe University, Rokkodai, Nada, Kobe 657-8501, Japan
- Core Research for Evolutionary Science and Technology, Japan Science and Technology Corp., Rokkodai, Nada, Kobe 657-8501, Japan
- To whom correspondence should be addressed at:
Department of Chemistry, Faculty of Science, Kobe University, Rokkodai, Nada, Kobe 657-8501, Japan. E-mail:
| |
Collapse
|
47
|
Jayachandran G, Vishal V, García AE, Pande VS. Local structure formation in simulations of two small proteins. J Struct Biol 2006; 157:491-9. [PMID: 17098444 PMCID: PMC1931413 DOI: 10.1016/j.jsb.2006.10.001] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2006] [Revised: 09/28/2006] [Accepted: 10/04/2006] [Indexed: 11/30/2022]
Abstract
Massively parallel all-atom, explicit solvent molecular dynamics simulations were used to explore the formation and existence of local structure in two small alpha-helical proteins, the villin headpiece and the helical fragment B of protein A. We report on the existence of transient helices and combinations of helices in the unfolded ensemble, and on the order of formation of helices, which appears to largely agree with previous experimental results. Transient local structure is observed even in the absence of overall native structure. We also calculate sets of residue-residue pairs that are statistically predictive of the formation of given local structures in our simulations.
Collapse
Affiliation(s)
- Guha Jayachandran
- Computer Science Department, Stanford University, Stanford, CA 94305, USA
| | | | | | | |
Collapse
|
48
|
Linhananta A, Boer J, MacKay I. The equilibrium properties and folding kinetics of an all-atom Go model of the Trp-cage. J Chem Phys 2006; 122:114901. [PMID: 15836251 DOI: 10.1063/1.1874812] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The ultrafast-folding 20-residue Trp-cage protein is quickly becoming a new benchmark for molecular dynamics studies. Already several all-atom simulations have probed its equilibrium and kinetic properties. In this work an all-atom Go model is used to accurately represent the side-chain packing and native atomic contacts of the Trp-cage. The model reproduces the hallmark thermodynamics cooperativity of small proteins. Folding simulations observe that in the fast-folding dominant pathway, partial alpha-helical structure forms before hydrophobic core collapse. In the slow-folding secondary pathway, partial core collapse occurs before helical structure. The slow-folding rate of the secondary pathway is attributed to the loss of side-chain rotational freedom, due to the early core collapse, which impedes the helix formation. A major finding is the observation of a low-temperature kinetic intermediate stabilized by a salt bridge between residues Asp-9 and Arg-16. Similar observations [R. Zhou, Proc. Natl. Acad. Sci. U.S.A. 100, 13280 (2003)] were reported in a recent study using an all-atom model of the Trp-cage in explicit water, in which the salt-bridge stabilized intermediate was hypothesized to be the origin of the ultrafast-folding mechanism. A theoretical mutation that eliminates the Asp-9-Arg-16 salt bridge, but leaves the residues intact, is performed. Folding simulations of the mutant Trp-cage observe a two-state free-energy landscape with no kinetic intermediate and a significant decrease in the folding rate, in support of the hypothesis.
Collapse
|
49
|
Sato S, Religa TL, Fersht AR. Phi-analysis of the folding of the B domain of protein A using multiple optical probes. J Mol Biol 2006; 360:850-64. [PMID: 16782128 DOI: 10.1016/j.jmb.2006.05.051] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2006] [Revised: 03/23/2006] [Accepted: 05/19/2006] [Indexed: 10/24/2022]
Abstract
We examined the co-operativity of ultra-fast folding of a protein and whether the Phi-value analysis of its transition state depended on the location of the optical probe. We incorporated in turn a tryptophan residue into each of the three helices of the B domain of Protein A. Each Trp mutant of the three-helix bundle protein was used as a pseudo-wild-type parent for Phi-analysis in which the intrinsic Trp fluorescence probed the formation of each helix during the transition state. Apart from local effects in the immediate vicinity of the probe, the three separate sets of Phi-values were in excellent agreement, demonstrating the overall co-operativity of folding and the robustness of the Phi-analysis. The transition state of folding of Protein A contains the second helix being well formed with many stabilizing tertiary hydrophobic interactions. In contrast, the first and the third helices are more poorly structured in the transition state. The mechanism of folding thus involves the concurrent formation of secondary and tertiary interactions, and is towards the nucleation-condensation extreme in the nucleation-condensation-framework continuum of mechanism, with helix 2 being the nucleus. We provide an error analysis of Phi-values derived purely from the kinetics of two-state chevron plots.
Collapse
Affiliation(s)
- Satoshi Sato
- MRC Centre for Protein Engineering, Hills Road, Cambridge CB2 2QH, UK
| | | | | |
Collapse
|
50
|
Affiliation(s)
- Yawen Bai
- Laboratory of Biochemistry, National Cancer Institute, Building 37, Room 6114E, National Institutes of Health, Bethesda, Maryland 20892, USA.
| |
Collapse
|