1
|
Kamenik AS, Handle PH, Hofer F, Kahler U, Kraml J, Liedl KR. Polarizable and non-polarizable force fields: Protein folding, unfolding, and misfolding. J Chem Phys 2021; 153:185102. [PMID: 33187403 DOI: 10.1063/5.0022135] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Molecular dynamics simulations are an invaluable tool to characterize the dynamic motions of proteins in atomistic detail. However, the accuracy of models derived from simulations inevitably relies on the quality of the underlying force field. Here, we present an evaluation of current non-polarizable and polarizable force fields (AMBER ff14SB, CHARMM 36m, GROMOS 54A7, and Drude 2013) based on the long-standing biophysical challenge of protein folding. We quantify the thermodynamics and kinetics of the β-hairpin formation using Markov state models of the fast-folding mini-protein CLN025. Furthermore, we study the (partial) folding dynamics of two more complex systems, a villin headpiece variant and a WW domain. Surprisingly, the polarizable force field in our set, Drude 2013, consistently leads to destabilization of the native state, regardless of the secondary structure element present. All non-polarizable force fields, on the other hand, stably characterize the native state ensembles in most cases even when starting from a partially unfolded conformation. Focusing on CLN025, we find that the conformational space captured with AMBER ff14SB and CHARMM 36m is comparable, but the ensembles from CHARMM 36m simulations are clearly shifted toward disordered conformations. While the AMBER ff14SB ensemble overstabilizes the native fold, CHARMM 36m and GROMOS 54A7 ensembles both agree remarkably well with experimental state populations. In addition, GROMOS 54A7 also reproduces experimental folding times most accurately. Our results further indicate an over-stabilization of helical structures with AMBER ff14SB. Nevertheless, the presented investigations strongly imply that reliable (un)folding dynamics of small proteins can be captured in feasible computational time with current additive force fields.
Collapse
Affiliation(s)
- Anna S Kamenik
- Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Philip H Handle
- Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Florian Hofer
- Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Ursula Kahler
- Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Johannes Kraml
- Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | - Klaus R Liedl
- Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| |
Collapse
|
2
|
Design and structural characterisation of monomeric water-soluble α-helix and β-hairpin peptides: State-of-the-art. Arch Biochem Biophys 2019; 661:149-167. [DOI: 10.1016/j.abb.2018.11.014] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2018] [Revised: 11/06/2018] [Accepted: 11/14/2018] [Indexed: 02/06/2023]
|
3
|
Anderson JM, Kier BL, Jurban B, Byrne A, Shu I, Eidenschink LA, Shcherbakov AA, Hudson M, Fesinmeyer RM, Andersen NH. Aryl-aryl interactions in designed peptide folds: Spectroscopic characteristics and optimal placement for structure stabilization. Biopolymers 2016; 105:337-356. [PMID: 26850220 PMCID: PMC5638712 DOI: 10.1002/bip.22821] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2015] [Revised: 02/01/2016] [Accepted: 02/03/2016] [Indexed: 01/27/2023]
Abstract
We have extended our studies of Trp/Trp to other Aryl/Aryl through-space interactions that stabilize hairpins and other small polypeptide folds. Herein we detail the NMR and CD spectroscopic features of these types of interactions. NMR data remains the best diagnostic for characterizing the common T-shape orientation. Designated as an edge-to-face (EtF or FtE) interaction, large ring current shifts are produced at the edge aryl ring hydrogens and, in most cases, large exciton couplets appear in the far UV circular dichroic (CD) spectrum. The preference for the face aryl in FtE clusters is W ≫ Y ≥ F (there are some exceptions in the Y/F order); this sequence corresponds to the order of fold stability enhancement and always predicts the amplitude of the lower energy feature of the exciton couplet in the CD spectrum. The CD spectra for FtE W/W, W/Y, Y/W, and Y/Y pairs all include an intense feature at 225-232 nm. An additional couplet feature seen for W/Y, W/F, Y/Y, and F/Y clusters, is a negative feature at 197-200 nm. Tyr/Tyr (as well as F/Y and F/F) interactions produce much smaller exciton couplet amplitudes. The Trp-cage fold was employed to search for the CD effects of other Trp/Trp and Trp/Tyr cluster geometries: several were identified. In this account, we provide additional examples of the application of cross-strand aryl/aryl clusters for the design of stable β-sheet models and a scale of fold stability increments associated with all possible FtE Ar/Ar clusters in several structural contexts. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 337-356, 2016.
Collapse
Affiliation(s)
- Jordan M Anderson
- Department of Chemistry, University of Washington, Seattle, WA, 98195
| | - Brandon L Kier
- Department of Chemistry, University of Washington, Seattle, WA, 98195
| | - Brice Jurban
- Department of Chemistry, University of Washington, Seattle, WA, 98195
| | - Aimee Byrne
- Department of Chemistry, University of Washington, Seattle, WA, 98195
| | - Irene Shu
- Department of Chemistry, University of Washington, Seattle, WA, 98195
| | | | | | - Mike Hudson
- Department of Chemistry, University of Washington, Seattle, WA, 98195
| | - R M Fesinmeyer
- Department of Chemistry, University of Washington, Seattle, WA, 98195
| | - Niels H Andersen
- Department of Chemistry, University of Washington, Seattle, WA, 98195
| |
Collapse
|
4
|
Otvos L, Kovalszky I, Olah J, Coroniti R, Knappe D, Nollmann FI, Hoffmann R, Wade JD, Lovas S, Surmacz E. Optimization of adiponectin-derived peptides for inhibition of cancer cell growth and signaling. Biopolymers 2015; 104:156-66. [DOI: 10.1002/bip.22627] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2015] [Revised: 01/27/2015] [Accepted: 02/09/2015] [Indexed: 12/16/2022]
Affiliation(s)
- Laszlo Otvos
- Department of Biology; Temple University; Philadelphia PA 19122
| | - Ilona Kovalszky
- 1st Department of Pathology and Experimental Cancer Research; Semmelweis University; Budapest 1085 Hungary
| | - Julia Olah
- 1st Department of Pathology and Experimental Cancer Research; Semmelweis University; Budapest 1085 Hungary
| | - Roberta Coroniti
- Sbarro Institute for Cancer Research and Molecular Medicine; Temple University; Philadelphia PA 19122
| | - Daniel Knappe
- Institute of Bioanalytical Chemistry; Leipzig University; Leipzig 04103 Germany
| | | | - Ralf Hoffmann
- Institute of Bioanalytical Chemistry; Leipzig University; Leipzig 04103 Germany
| | - John D. Wade
- Florey Neurosciences Institutes; University of Melbourne; Victoria 3010 Australia
- School of Chemistry; University of Melbourne; Victoria 3010 Australia
| | - Sandor Lovas
- Department of Biomedical Sciences; Creighton University; Omaha NE 68178
| | - Eva Surmacz
- Sbarro Institute for Cancer Research and Molecular Medicine; Temple University; Philadelphia PA 19122
| |
Collapse
|
5
|
Nishio M, Umezawa Y, Fantini J, Weiss MS, Chakrabarti P. CH-π hydrogen bonds in biological macromolecules. Phys Chem Chem Phys 2015; 16:12648-83. [PMID: 24836323 DOI: 10.1039/c4cp00099d] [Citation(s) in RCA: 335] [Impact Index Per Article: 37.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
This is a sequel to the previous Perspective "The CH-π hydrogen bond in chemistry. Conformation, supramolecules, optical resolution and interactions involving carbohydrates", which featured in a PCCP themed issue on "Weak Hydrogen Bonds - Strong Effects?": Phys. Chem. Chem. Phys., 2011, 13, 13873-13900. Evidence that weak hydrogen bonds play an enormously important role in chemistry and biochemistry has now accumulated to an extent that the rigid classical concept of hydrogen bonds formulated by Pauling needs to be seriously revised and extended. The concept of a more generalized hydrogen bond definition is indispensable for understanding the folding mechanisms of proteins. The CH-π hydrogen bond, a weak molecular force occurring between a soft acid CH and a soft base π-electron system, among all is one of the most important and plays a functional role in defining the conformation and stability of 3D structures as well as in many molecular recognition events. This concept is also valuable in structure-based drug design efforts. Despite their frequent occurrence in organic molecules and bio-molecules, the importance of CH-π hydrogen bonds is still largely unknown to many chemists and biochemists. Here we present a review that deals with the evidence, nature, characteristics and consequences of the CH-π hydrogen bond in biological macromolecules (proteins, nucleic acids, lipids and polysaccharides). It is hoped that the present Perspective will show the importance of CH-π hydrogen bonds and stimulate interest in the interactions of biological macromolecules, one of the most fascinating fields in bioorganic chemistry. Implication of this concept is enormous and valuable in the scientific community.
Collapse
Affiliation(s)
- Motohiro Nishio
- The CHPI Institute, 705-6-338, Minamioya, Machida-shi, Tokyo 194-0031, Japan.
| | | | | | | | | |
Collapse
|
6
|
Yasuda S, Hayashi T, Kinoshita M. Physical origins of the high structural stability of CLN025 with only ten residues. J Chem Phys 2014; 141:105103. [DOI: 10.1063/1.4894753] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Satoshi Yasuda
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Tomohiko Hayashi
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Masahiro Kinoshita
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| |
Collapse
|
7
|
Abstract
Since the first report in 1993 (JACS 115, 5887-5888) of a peptide able to form a monomeric β-hairpin structure in aqueous solution, the design of peptides forming either β-hairpins (two-stranded antiparallel β-sheets) or three-stranded antiparallel β-sheets has become a field of growing interest and activity. These studies have yielded great insights into the principles governing the stability and folding of β-hairpins and antiparallel β-sheets. This chapter provides an overview of the reported β-hairpin/β-sheet peptides focussed on the applied design criteria, reviews briefly the factors contributing to β-hairpin/β-sheet stability, and describes a protocol for the de novo design of β-sheet-forming peptides based on them. Guidelines to select appropriate turn and strand residues and to avoid self-association are provided. The methods employed to check the success of new designed peptides are also summarized. Since NMR is the best technique to that end, NOEs and chemical shifts characteristic of β-hairpins and three-stranded antiparallel β-sheets are given.
Collapse
Affiliation(s)
- M Angeles Jiménez
- Consejo Superior de Investigaciones Científicas (CSIC), Instituto de Química Física Rocasolano (IQFR), Serrano 119, 28006, Madrid, Spain,
| |
Collapse
|
8
|
Mandal M, Mukhopadhyay C. Microsecond molecular dynamics simulation of guanidinium chloride induced unfolding of ubiquitin. Phys Chem Chem Phys 2014; 16:21706-16. [DOI: 10.1039/c4cp01657b] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
All atom molecular dynamics simulations have been used to explore the atomic detail mechanism of guanidinium induced unfolding of the protein ubiquitin.
Collapse
Affiliation(s)
- Manoj Mandal
- Department of Chemistry
- University of Calcutta
- Kolkata – 700 009, India
| | | |
Collapse
|
9
|
Leveles I, Németh V, Szabó JE, Harmat V, Nyíri K, Bendes ÁÁ, Papp-Kádár V, Zagyva I, Róna G, Ozohanics O, Vékey K, Tóth J, Vértessy BG. Structure and enzymatic mechanism of a moonlighting dUTPase. ACTA CRYSTALLOGRAPHICA SECTION D: BIOLOGICAL CRYSTALLOGRAPHY 2013; 69:2298-308. [PMID: 24311572 DOI: 10.1107/s0907444913021136] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Received: 03/14/2013] [Accepted: 07/29/2013] [Indexed: 02/08/2023]
Abstract
Genome integrity requires well controlled cellular pools of nucleotides. dUTPases are responsible for regulating cellular dUTP levels and providing dUMP for dTTP biosynthesis. In Staphylococcus, phage dUTPases are also suggested to be involved in a moonlighting function regulating the expression of pathogenicity-island genes. Staphylococcal phage trimeric dUTPase sequences include a specific insertion that is not found in other organisms. Here, a 2.1 Å resolution three-dimensional structure of a ϕ11 phage dUTPase trimer with complete localization of the phage-specific insert, which folds into a small β-pleated mini-domain reaching out from the dUTPase core surface, is presented. The insert mini-domains jointly coordinate a single Mg2+ ion per trimer at the entrance to the threefold inner channel. Structural results provide an explanation for the role of Asp95, which is suggested to have functional significance in the moonlighting activity, as the metal-ion-coordinating moiety potentially involved in correct positioning of the insert. Enzyme-kinetics studies of wild-type and mutant constructs show that the insert has no major role in dUTP binding or cleavage and provide a description of the elementary steps (fast binding of substrate and release of products). In conclusion, the structural and kinetic data allow insights into both the phage-specific characteristics and the generally conserved traits of ϕ11 phage dUTPase.
Collapse
Affiliation(s)
- Ibolya Leveles
- Institute of Enzymology, Research Centre for Natural Sciences, Hungarian Academy of Sciences, 29 Karolina Street, 1113 Budapest, Hungary
| | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|
10
|
Lovas S, Zhang Y, Yu J, Lyubchenko YL. Molecular mechanism of misfolding and aggregation of Aβ(13-23). J Phys Chem B 2013; 117:6175-86. [PMID: 23642026 DOI: 10.1021/jp402938p] [Citation(s) in RCA: 40] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
The misfolding and self-assembly of the amyloid-beta (Aβ) peptide into aggregates is a molecular signature of the development of Alzheimer's disease, but molecular mechanisms of the peptide aggregation remain unknown. Here, we combined Atomic Force Microscopy (AFM) and Molecular Dynamics (MD) simulations to characterize the misfolding process of an Aβ peptide. Dynamic force spectroscopy AFM analysis showed that the peptide forms stable dimers with a lifetime of ∼1 s. During MD simulations, isolated monomers gradually adopt essentially similar nonstructured conformations independent from the initial structure. However, when two monomers approach their structure changes dramatically, and the conformational space for the two monomers become restricted. The arrangement of monomers in antiparallel orientation leads to the cooperative formation of β-sheet conformation. Interactions, including hydrogen bonds, salt bridges, and weakly polar interactions of side chains stabilize the structure of the dimer. Under the applied force, the dimer, as during the AFM experiments, dissociates in a cooperative manner. Thus, misfolding of the Aβ peptide proceeds via the loss of conformational flexibility and formation of stable dimers suggesting their key role in the subsequent Aβ aggregation process.
Collapse
Affiliation(s)
- Sándor Lovas
- Department of Biomedical Sciences, Creighton University, Omaha, Nebraska 68178, United States.
| | | | | | | |
Collapse
|
11
|
Xia Z, Das P, Shakhnovich EI, Zhou R. Collapse of unfolded proteins in a mixture of denaturants. J Am Chem Soc 2012; 134:18266-74. [PMID: 23057830 DOI: 10.1021/ja3031505] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Both urea and guanidinium chloride (GdmCl) are frequently used as protein denaturants. Given that proteins generally adopt extended or unfolded conformations in either aqueous urea or GdmCl, one might expect that the unfolded protein chains will remain or become further extended due to the addition of another denaturant. However, a collapse of denatured proteins is revealed using atomistic molecular dynamics simulations when a mixture of denaturants is used. Both hen egg-white lysozyme and protein L are found to undergo collapse in the denaturant mixture. The collapse of the protein conformational ensembles is accompanied by a decreased solubility and increased non-native self-interactions of hydrophobic residues in the urea/GdmCl mixture. The increase of non-native interactions rather than the native contacts indicates that the proteins experience a simple collapse transition from the fully denatured states. During the protein collapse, the relatively stronger denaturant GdmCl displays a higher tendency to be absorbed onto the protein surface due to their stronger electrostatic interactions with proteins. At the same time, urea molecules also accumulate near the protein surface, resulting in an enhanced "local crowding" for the protein near its first solvation shell. This rearrangement of denaturants near the protein surface and crowded local environment induce the protein collapse, mainly by burying their hydrophobic residues. These findings from molecular simulations are then further explained by a simple analytical model based on statistical mechanics.
Collapse
Affiliation(s)
- Zhen Xia
- Computational Biology Center, IBM Thomas J. Watson Research Center, Yorktown Heights, New York 10598, USA
| | | | | | | |
Collapse
|
12
|
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
|