1
|
Application of conventional molecular dynamics simulation in evaluating the stability of apomyoglobin in urea solution. Sci Rep 2017; 7:44651. [PMID: 28300210 PMCID: PMC5353640 DOI: 10.1038/srep44651] [Citation(s) in RCA: 61] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2016] [Accepted: 02/09/2017] [Indexed: 01/02/2023] Open
Abstract
In this study, we had exploited the advancement in computer technology to determine the stability of four apomyoglobin variants namely wild type, E109A, E109G and G65A/G73A by conducting conventional molecular dynamics simulations in explicit urea solution. Variations in RMSD, native contacts and solvent accessible surface area of the apomyoglobin variants during the simulation were calculated to probe the effect of mutation on the overall conformation of the protein. Subsequently, the mechanism leading to the destabilization of the apoMb variants was studied through the calculation of correlation matrix, principal component analyses, hydrogen bond analyses and RMSF. The results obtained here correlate well with the study conducted by Baldwin and Luo which showed improved stability of apomyoglobin with E109A mutation and contrariwise for E109G and G65A/G73A mutation. These positive observations showcase the feasibility of exploiting MD simulation in determining protein stability prior to protein expression.
Collapse
|
2
|
Structure and dynamics of a free aquaporin (AQP1) by a coarse-grained Monte Carlo simulation. Struct Chem 2016. [DOI: 10.1007/s11224-016-0836-4] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
3
|
Abstract
Motivation: Comparing protein tertiary structures is a fundamental procedure in structural biology and protein bioinformatics. Structure comparison is important particularly for evaluating computational protein structure models. Most of the model structure evaluation methods perform rigid body superimposition of a structure model to its crystal structure and measure the difference of the corresponding residue or atom positions between them. However, these methods neglect intrinsic flexibility of proteins by treating the native structure as a rigid molecule. Because different parts of proteins have different levels of flexibility, for example, exposed loop regions are usually more flexible than the core region of a protein structure, disagreement of a model to the native needs to be evaluated differently depending on the flexibility of residues in a protein. Results: We propose a score named FlexScore for comparing protein structures that consider flexibility of each residue in the native state of proteins. Flexibility information may be extracted from experiments such as NMR or molecular dynamics simulation. FlexScore considers an ensemble of conformations of a protein described as a multivariate Gaussian distribution of atomic displacements and compares a query computational model with the ensemble. We compare FlexScore with other commonly used structure similarity scores over various examples. FlexScore agrees with experts’ intuitive assessment of computational models and provides information of practical usefulness of models. Availability and implementation:https://bitbucket.org/mjamroz/flexscore Contact:dkihara@purdue.edu Supplementary information:Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Michal Jamroz
- Department of Chemistry, University of Warsaw, Warsaw, 02-093, Poland
| | - Andrzej Kolinski
- Department of Chemistry, University of Warsaw, Warsaw, 02-093, Poland
| | - Daisuke Kihara
- Department of Biological Sciences Department of Computer Science, Purdue University, West Lafayette, IN 47907, USA
| |
Collapse
|
4
|
Morgan BR, Deveau LM, Massi F. Probing the structural and dynamical effects of the charged residues of the TZF domain of TIS11d. Biophys J 2015; 108:1503-1515. [PMID: 25809263 DOI: 10.1016/j.bpj.2015.01.039] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2014] [Revised: 01/15/2015] [Accepted: 01/21/2015] [Indexed: 12/25/2022] Open
Abstract
A member of the TTP family of proteins, TIS11d binds RNA with high specificity using a pair of CCCH-type tandem zinc fingers separated by a 18 residue long linker. Our previous work showed that the formation of hydrogen bonds between the C-terminal residue E220 and the residues of the linker region stabilized a compact structure of TIS11d in the absence of RNA. To investigate the role of the C-terminal residues in the structure of unbound TIS11d, the E220A mutant and the truncation mutant lacking the last two residues (D219/E220) were studied using molecular dynamics, NMR spectroscopy, and biochemical methods. This study confirmed the importance of the charged residues D219 and E220 in maintaining structural stability in unbound TIS11d and elucidated the underlying physical mechanisms. We observed a greater structural heterogeneity for the residues of the linker in the molecular dynamics trajectories of both mutant proteins relative to the wild-type. This heterogeneity was more pronounced in the D219/E220 deletion mutant than in the E220A mutant, indicating that a greater reduction of the charge of the C-terminus results in greater flexibility. In agreement with the increased flexibility and the reduced number of negatively charged residues of the D219/E220 deletion mutant, we measured more unfavorable entropic and a more favorable enthalpic contribution to the free energy of RNA binding in the mutant than in the wild-type protein. The relative orientation of the zinc fingers was stabilized by the electrostatic interaction between E220 and positively charged residues of the linker in TIS11d. In the E220A mutant, the relative orientation of the zinc fingers was less constrained, whereas in the D219/E220 deletion mutant, little orientational preference was observed. We posit that favorable electrostatic interactions provide a mechanism to promote preferential orientation of separate domains without imposing structural rigidity.
Collapse
Affiliation(s)
- Brittany R Morgan
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts, Worcester, Massachusetts
| | - Laura M Deveau
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts, Worcester, Massachusetts
| | - Francesca Massi
- Department of Biochemistry and Molecular Pharmacology, University of Massachusetts, Worcester, Massachusetts.
| |
Collapse
|
5
|
Uyar A, Kantarci-Carsibasi N, Haliloglu T, Doruker P. Features of large hinge-bending conformational transitions. Prediction of closed structure from open state. Biophys J 2015; 106:2656-66. [PMID: 24940783 DOI: 10.1016/j.bpj.2014.05.017] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2013] [Revised: 04/22/2014] [Accepted: 05/08/2014] [Indexed: 12/11/2022] Open
Abstract
We performed a detailed analysis of conformational transition pathways for a set of 10 proteins, which undergo large hinge-bending-type motions with 4-12 Å RMSD (root mean-square distance) between open and closed crystal structures. Anisotropic network model-Monte Carlo (ANM-MC) algorithm generates a targeted pathway between two conformations, where the collective modes from the ANM are used for deformation at each iteration and the conformational energy of the deformed structure is minimized via an MC algorithm. The target structure was approached successfully with an RMSD of 0.9-4.1 Å when a relatively low cutoff radius of 10 Å was used in ANM. Even though one predominant mode (first or second) directed the open-to-closed conformational transition, changes in the dominant mode character were observed for most cases along the transition. By imposing radius of gyration constraint during mode selection, it was possible to predict the closed structure for eight out of 10 proteins (with initial 4.1-7.1 Å and final 1.7-2.9 Å RMSD to target). Deforming along a single mode leads to most successful predictions. Based on the previously reported free energy surface of adenylate kinase, deformations along the first mode produced an energetically favorable path, which was interestingly facilitated by a change in mode shape (resembling second and third modes) at key points. Pathway intermediates are provided in our database of conformational transitions (http://safir.prc.boun.edu.tr/anmmc/method/1).
Collapse
Affiliation(s)
- Arzu Uyar
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, Istanbul, Turkey
| | - Nigar Kantarci-Carsibasi
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, Istanbul, Turkey
| | - Turkan Haliloglu
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, Istanbul, Turkey.
| | - Pemra Doruker
- Department of Chemical Engineering and Polymer Research Center, Bogazici University, Istanbul, Turkey.
| |
Collapse
|
6
|
Ha-Duong T. Coarse-grained models of the proteins backbone conformational dynamics. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2014; 805:157-69. [PMID: 24446361 DOI: 10.1007/978-3-319-02970-2_7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Coarse-grained models are more and more frequently used in the studies of the proteins structural and dynamic properties, since the reduced number of degrees of freedom allows to enhance the conformational space exploration. This chapter attempts to provide an overview of the various coarse-grained models that were applied to study the functional conformational changes of the polypeptides main chain around their native state. It will more specifically discuss the methods used to represent the protein backbone flexibility and to account for the physico-chemical interactions that stabilize the secondary structure elements.
Collapse
Affiliation(s)
- Tap Ha-Duong
- BIOCIS - UMR CNRS 8076, Faculté de Pharmacie - Université Paris Sud, 5 rue Jean-Baptiste Clément, 92296, Châtenay-Malabry, France,
| |
Collapse
|
7
|
Moskovitz Y, Srebnik S. Conformational changes of globular proteins upon adsorption on a hydrophobic surface. Phys Chem Chem Phys 2014; 16:11698-707. [DOI: 10.1039/c4cp00354c] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Coarse-grained Monte Carlo simulations are used to study thermal denaturation of small globular proteins adsorbed on a hydrophobic surface. Though helices are more stable than sheets, they are highly deformed in the adsorbed protein.
Collapse
Affiliation(s)
- Yevgeny Moskovitz
- Department of Chemistry
- Scientific Computing Research Unit
- University of Cape Town
- Rondebosch 7701, South Africa
| | - Simcha Srebnik
- Department of Chemical Engineering
- Technion – Israel Institute of Technology
- Haifa 32000, Israel
| |
Collapse
|
8
|
Allosteric control in a metalloprotein dramatically alters function. Proc Natl Acad Sci U S A 2012; 110:948-53. [PMID: 23271805 DOI: 10.1073/pnas.1208286110] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Metalloproteins (MPs) comprise one-third of all known protein structures. This diverse set of proteins contain a plethora of unique inorganic moieties capable of performing chemistry that would otherwise be impossible using only the amino acids found in nature. Most of the well-studied MPs are generally viewed as being very rigid in structure, and it is widely thought that the properties of the metal centers are primarily determined by the small fraction of amino acids that make up the local environment. Here we examine both theoretically and experimentally whether distal regions can influence the metal center in the diabetes drug target mitoNEET. We demonstrate that a loop (L2) 20 Å away from the metal center exerts allosteric control over the cluster binding domain and regulates multiple properties of the metal center. Mutagenesis of L2 results in significant shifts in the redox potential of the [2Fe-2S] cluster and orders of magnitude effects on the rate of [2Fe-2S] cluster transfer to an apo-acceptor protein. These surprising effects occur in the absence of any structural changes. An examination of the native basin dynamics of the protein using all-atom simulations shows that twisting in L2 controls scissoring in the cluster binding domain and results in perturbations to one of the cluster-coordinating histidines. These allosteric effects are in agreement with previous folding simulations that predicted L2 could communicate with residues surrounding the metal center. Our findings suggest that long-range dynamical changes in the protein backbone can have a significant effect on the functional properties of MPs.
Collapse
|
9
|
Moskovitz Y, Srebnik S. Thermal stability limits of proteins in solution and adsorbed on a hydrophobic surface. Phys Chem Chem Phys 2012; 14:8013-22. [DOI: 10.1039/c2cp00005a] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
10
|
Nechushtai R, Lammert H, Michaeli D, Eisenberg-Domovich Y, Zuris JA, Luca MA, Capraro DT, Fish A, Shimshon O, Roy M, Schug A, Whitford PC, Livnah O, Onuchic JN, Jennings PA. Allostery in the ferredoxin protein motif does not involve a conformational switch. Proc Natl Acad Sci U S A 2011; 108:2240-2245. [PMID: 21266547 PMCID: PMC3038707 DOI: 10.1073/pnas.1019502108] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2023] Open
Abstract
Regulation of protein function via cracking, or local unfolding and refolding of substructures, is becoming a widely recognized mechanism of functional control. Oftentimes, cracking events are localized to secondary and tertiary structure interactions between domains that control the optimal position for catalysis and/or the formation of protein complexes. Small changes in free energy associated with ligand binding, phosphorylation, etc., can tip the balance and provide a regulatory functional switch. However, understanding the factors controlling function in single-domain proteins is still a significant challenge to structural biologists. We investigated the functional landscape of a single-domain plant-type ferredoxin protein and the effect of a distal loop on the electron-transfer center. We find the global stability and structure are minimally perturbed with mutation, whereas the functional properties are altered. Specifically, truncating the L1,2 loop does not lead to large-scale changes in the structure, determined via X-ray crystallography. Further, the overall thermal stability of the protein is only marginally perturbed by the mutation. However, even though the mutation is distal to the iron-sulfur cluster (∼20 Å), it leads to a significant change in the redox potential of the iron-sulfur cluster (57 mV). Structure-based all-atom simulations indicate correlated dynamical changes between the surface-exposed loop and the iron-sulfur cluster-binding region. Our results suggest intrinsic communication channels within the ferredoxin fold, composed of many short-range interactions, lead to the propagation of long-range signals. Accordingly, protein interface interactions that involve L1,2 could potentially signal functional changes in distal regions, similar to what is observed in other allosteric systems.
Collapse
Affiliation(s)
- Rachel Nechushtai
- Life Science Institute and The Wolfson Centre for Applied Structural Biology, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat Ram, Jerusalem 91904, Israel
| | - Heiko Lammert
- Center for Theoretical Biological Physics and the Department of Physics, University of California at San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0375
| | - Dorit Michaeli
- Life Science Institute and The Wolfson Centre for Applied Structural Biology, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat Ram, Jerusalem 91904, Israel
| | - Yael Eisenberg-Domovich
- Life Science Institute and The Wolfson Centre for Applied Structural Biology, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat Ram, Jerusalem 91904, Israel
| | - John A. Zuris
- Department of Chemistry and Biochemistry, University of California at San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0375
| | - Maria A. Luca
- Department of Chemistry and Biochemistry, University of California at San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0375
| | - Dominique T. Capraro
- Department of Chemistry and Biochemistry, University of California at San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0375
| | - Alex Fish
- Life Science Institute and The Wolfson Centre for Applied Structural Biology, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat Ram, Jerusalem 91904, Israel
| | - Odelia Shimshon
- Life Science Institute and The Wolfson Centre for Applied Structural Biology, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat Ram, Jerusalem 91904, Israel
| | - Melinda Roy
- Department of Chemistry and Biochemistry, University of California at San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0375
| | - Alexander Schug
- Center for Theoretical Biological Physics and the Department of Physics, University of California at San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0375
- Department of Chemistry, Umeå University, Umeå, Sweden; and
| | - Paul C. Whitford
- Los Alamos National Laboratory, Theoretical Biology and Biophysics, MS K710, Los Alamos, NM 87545
| | - Oded Livnah
- Life Science Institute and The Wolfson Centre for Applied Structural Biology, The Hebrew University of Jerusalem, Edmond J. Safra Campus, Givat Ram, Jerusalem 91904, Israel
| | - José N. Onuchic
- Center for Theoretical Biological Physics and the Department of Physics, University of California at San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0375
| | - Patricia A. Jennings
- Department of Chemistry and Biochemistry, University of California at San Diego, 9500 Gilman Drive, La Jolla, CA 92093-0375
| |
Collapse
|
11
|
Ha-Duong T. Protein Backbone Dynamics Simulations Using Coarse-Grained Bonded Potentials and Simplified Hydrogen Bonds. J Chem Theory Comput 2010; 6:761-73. [DOI: 10.1021/ct900408s] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Affiliation(s)
- Tap Ha-Duong
- Laboratoire Analyse et Modélisation pour la Biologie et l’Environnement Université d’Evry-Val-d’Essonne Rue du Pere André Jarlan, 91025 Evry Cedex, France
| |
Collapse
|
12
|
|
13
|
Yogurtcu ON, Gur M, Erman B. Statistical thermodynamics of residue fluctuations in native proteins. J Chem Phys 2009; 130:095103. [PMID: 19275429 DOI: 10.1063/1.3078517] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Statistical thermodynamics of residue fluctuations of native proteins in a temperature, pressure, and force reservoir is formulated. The general theory is discussed in terms of harmonic and anharmonic fluctuations of residues. The two elastic network models based on the harmonic approximation, the anisotropic network and the Gaussian network models are discussed as the limiting cases of the general theory. The heat capacity and the correlations between the energy fluctuations and residue fluctuations are obtained for the harmonic approximation. The formulation is extended to large fluctuations of residues in order to account for effects of anharmonicity. The fluctuation probability function is constructed for this purpose as a tensorial Hermite series expansion with higher order moments of fluctuations as coefficients. Evaluation of the higher order moments using the proposed statistical thermodynamics model is explained. The formulation is applied to a hexapeptide and the fluctuations of residues obtained by molecular dynamics simulations are characterized in the framework of the model developed. In particular, coupling of two different modes in the nonlinear model is discussed in detail.
Collapse
Affiliation(s)
- Osman N Yogurtcu
- Center for Computational Biology and Bioinformatics, Koc University, Istanbul, Turkey
| | | | | |
Collapse
|
14
|
Conformational transition pathways explored by Monte Carlo simulation integrated with collective modes. Biophys J 2008; 95:5862-73. [PMID: 18676657 DOI: 10.1529/biophysj.107.128447] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
Conformational transitions between open/closed or free/bound states in proteins possess functional importance. We propose a technique in which the collective modes obtained from an anisotropic network model (ANM) are used in conjunction with a Monte Carlo (MC) simulation approach, to investigate conformational transition pathways and pathway intermediates. The ANM-MC technique is applied to adenylate kinase (AK) and hemoglobin. The iterative method, in which normal modes are continuously updated during the simulation, proves successful in accomplishing the transition between open-closed conformations of AK and tense-relaxed forms of hemoglobin (C(alpha)-root mean square deviations between two end structures of 7.13 A and 3.55 A, respectively). Target conformations are reached by root mean-square deviations of 2.27 A and 1.90 A for AK and hemoglobin, respectively. The intermediate conformations overlap with crystal structures from the AK family within a 3.0-A root mean-square deviation. In the case of hemoglobin, the transition of tense-to-relaxed passes through the relaxed state. In both cases, the lowest-frequency modes are effective during transitions. The targeted Monte Carlo approach is used without the application of collective modes. Both the ANM-MC and targeted Monte Carlo techniques can explore sequences of events in transition pathways with an efficient yet realistic conformational search.
Collapse
|
15
|
Kuznetsov IB. Ordered conformational change in the protein backbone: Prediction of conformationally variable positions from sequence and low-resolution structural data. Proteins 2008; 72:74-87. [DOI: 10.1002/prot.21899] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
|
16
|
Curcó D, Nussinov R, Aleman C. Coarse-grained representation of beta-helical protein building blocks. J Phys Chem B 2007; 111:10538-49. [PMID: 17691836 DOI: 10.1021/jp072832q] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A general strategy to develop coarse-grained models of beta-helical protein fragments is presented. The procedure has been applied to a building block formed by a two-turn repeat motif from E. coli galactoside acetyltransferase, which is able to provide a very stable self-assembled tubular nanoconstruct upon stacking of its replicas. For this purpose, first, we have developed a computational scheme to sample very efficiently the configurational space of the building block. This method, which is inspired by a strategy recently designed to study amorphous polymers and by an advanced Monte Carlo algorithm, provides a large ensemble of uncorrelated configurations at a very reasonable computational cost. The atomistic configurations provided by this method have been used to obtain a coarse-grained model that describes the amino acids with fewer particles than those required for full atomistic detail, i.e., two, three, or four depending on the chemical nature of the amino acid. Coarse-grained potentials have been developed considering the following types of interactions: (i) electrostatic and van der Waals interactions between residues i and i + n with n >/= 2; (ii) interactions between residues i and i + 1; and (c) intra-residue interactions. The reliability of the proposed model has been tested by comparing the atomistic and coarse-grained energies calculated for a large number of independent configurations of the beta-helical building block.
Collapse
Affiliation(s)
- David Curcó
- Departament d'Enginyeria Química, Facultat de Química, Universitat de Barcelona, Martí i Franquees 1, Barcelona E-08028, Spain.
| | | | | |
Collapse
|
17
|
Güner U, Arkun Y, Erman B. Optimum folding pathways of proteins: Their determination and properties. J Chem Phys 2006; 124:134911. [PMID: 16613481 DOI: 10.1063/1.2181976] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We develop a dynamic optimization technique for determining optimum folding pathways of proteins starting from different initial configurations. A coarse-grained Go model is used. Forces acting on each bead are (i) the friction force, (ii) forces from bond length constraints, (iii) excluded volume constraints, and (iv) attractive forces between residue pairs that are in contact in the native state. An objective function is defined as the total attractive energy between nonbonded residues, which are neighbors in the native state. The objective function is minimized over all feasible paths, satisfying bond length and excluded volume constraints. The optimization problem is nonconvex and contains a large number of constraints. An augmented Lagrangian method with a penalty barrier function was used to solve the problem. The method is applied to a 36-residue protein, chicken villin headpiece. Sequences of events during folding of the protein are determined for various pathways and analyzed. The relative time scales are compared and scaled according to experimentally measured events. Formation times of the helices, turn, and the loop agree with experimental data. We obtain the overall folding time of the protein in the range of 600 ns-1.2 micros that is smaller than the experimental result of 4-5 micros, showing that the optimal folding times that we obtain may be possible lower bounds. Time dependent variables during folding and energies associated with short- and long-range interactions between secondary structures are analyzed in modal space using Karhunen-Loeve expansion.
Collapse
Affiliation(s)
- Uğur Güner
- College of Engineering, Koc University, Rumeli Feneri Yolu, 34450 Sariyer, Istanbul, Turkey
| | | | | |
Collapse
|
18
|
Shental-Bechor D, Kirca S, Ben-Tal N, Haliloglu T. Monte Carlo studies of folding, dynamics, and stability in alpha-helices. Biophys J 2005; 88:2391-402. [PMID: 15653741 PMCID: PMC1305339 DOI: 10.1529/biophysj.104.050708] [Citation(s) in RCA: 37] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Folding simulations of polyalanine peptides were carried out using an off-lattice Monte Carlo simulation technique. The peptide was represented as a chain of residues, each of which contains two interaction sites: one corresponding to the C(alpha) atom and the other to the side chain. A statistical potential was used to describe the interaction between these sites. The preferred conformations of the peptide chain on the energy surface, starting from several initial conditions, were searched by perturbations on its generalized coordinates with the Metropolis criterion. We observed that, at low temperatures, the effective energy was low and the helix content high. The calculated helix propagation (s) and nucleation (sigma) parameters of the Zimm-Bragg model were in reasonable agreement with the empirical data. Exploration of the energy surface of the alanine-based peptides (AAQAA)(3) and AAAAA(AAARA)(3)A demonstrated that their behavior is similar to that of polyalanine, in regard to their effective energy, helix content, and the temperature-dependence of their helicity. In contrast, stable secondary structures were not observed for (Gly)(20) at similar temperatures, which is consistent with the nonfolder nature of this peptide. The fluctuations in the slowest dynamics mode, which describe the elastic behavior of the chain, showed that as the temperature decreases, the polyalanine peptides become stiffer and retain conformations with higher helix content. Clustering of conformations during the folding phase implied that polyalanine folds into a helix through fewer numbers of intermediate conformations as the temperature decreases.
Collapse
Affiliation(s)
- Dalit Shental-Bechor
- Department of Biochemistry, George S. Wise Faculty of Life Sciences, Tel Aviv University, Ramat Aviv, Israel
| | | | | | | |
Collapse
|
19
|
Kessel A, Shental-Bechor D, Haliloglu T, Ben-Tal N. Interactions of hydrophobic peptides with lipid bilayers: Monte Carlo simulations with M2delta. Biophys J 2004; 85:3431-44. [PMID: 14645040 PMCID: PMC1303652 DOI: 10.1016/s0006-3495(03)74765-1] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
Abstract
We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.
Collapse
Affiliation(s)
- Amit Kessel
- Department of Biochemistry, George S Wise Faculty of Life Sciences, Tel Aviv University, Ramat Aviv, Israel
| | | | | | | |
Collapse
|
20
|
Kolinski A, Klein P, Romiszowski P, Skolnick J. Unfolding of globular proteins: monte carlo dynamics of a realistic reduced model. Biophys J 2004; 85:3271-8. [PMID: 14581227 PMCID: PMC1303603 DOI: 10.1016/s0006-3495(03)74745-6] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022] Open
Abstract
Reduced lattice models of proteins and Monte Carlo dynamics were used to simulate the initial stages of the unfolding of several proteins of various structural types, and the results were compared to experiment. The models semiquantitatively reproduce the approximate order of events of unfolding as well as subtle mutation effects and effects resulting from differences in sequences of similar folds. The short-time mobility of particular residues, observed in simulations, correlates with the crystallographic temperature factor. The main factor controlling unfolding is the native state topology, with sequence playing a less important role. The correlation with various experiments, especially for sequence-specific effects, strongly suggests that properly designed reduced models of proteins can be used for qualitative studies (or prediction) of protein unfolding pathways.
Collapse
Affiliation(s)
- Andrzej Kolinski
- Laboratory of Theory of Biopolymers, Faculty of Chemistry, Warsaw University, Warsaw, Poland.
| | | | | | | |
Collapse
|
21
|
Ozkan S, Dalgýn GS, Haliloglu T. Unfolding events of Chymotrypsin Inhibitor 2 (CI2) revealed by Monte Carlo (MC) simulations and their consistency from structure-based analysis of conformations. POLYMER 2004. [DOI: 10.1016/j.polymer.2003.10.092] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
22
|
Kurt N, Haliloglu T, Schiffer CA. Structure-based prediction of potential binding and nonbinding peptides to HIV-1 protease. Biophys J 2003; 85:853-63. [PMID: 12885633 PMCID: PMC1303207 DOI: 10.1016/s0006-3495(03)74525-1] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022] Open
Abstract
HIV-1 protease is a major drug target against AIDS as it permits viral maturation by processing the gag and pol polyproteins of the virus. The cleavage sites in these polyproteins do not have obvious sequence homology or a binding motif and the specificity of the protease is not easily determined. We used various threading approaches, together with the crystal structures of substrate complexes which served as template structures, to study the substrate specificity of HIV-1 protease with the aim of obtaining a better differentiation between binding and nonbinding sequences. The predictions from threading improved when distance-dependent interaction energy functions were used instead of contact matrices. To rank the peptides and properly account for the peptide's conformation in the total energy, the results from using short-range potentials on multiple template structures were averaged. Finally, a dynamic threading approach is introduced which is potentially useful for cases when there is only one template structure available. The conformational energy of the peptide-especially the term accounting for the side chains-was found to be important in differentiating between binding and nonbinding sequences. Hence, the substrate specificity, and thus the ability of the virus to mature, is affected by the compatibility of the substrate peptide to fit within the limited conformational space of the active site groove.
Collapse
Affiliation(s)
- Nese Kurt
- Polymer Research Center, Bogazici University, Bebek, Istanbul, Turkey
| | | | | |
Collapse
|
23
|
Abstract
6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) belongs to a class of catalytic enzymes involved in phosphoryl transfer and is a new target for the development of novel antimicrobial agents. In the present study, the fundamental consideration is to view the overall structure of HPPK as a network of interacting residues and to extract the most cooperative collective motions that define its global dynamics. A coarse-grained model, harmonically constrained according to HPPK's crystal structure is used. Four crystal structures of HPPK (one apo and three holo forms with different nucleotide and pterin analogs) are studied with the goal of providing insights about the function-dynamic correlation and ligand induced conformational changes. The dynamic differences are examined between HPPK's apo- and holo-forms, because they are involved in the catalytic reaction steps. Our results indicate that the palm-like structure of HPPK is nearly rigid, whereas the two flexible loops: L2 (residues 43-53) and L3 (residues 82-92) exhibit the most concerted motions for ligand recognition and presumably, catalysis. These two flexible loops are involved in the recognition of HPPKs nucleotide and pterin ligands, whereas the rigid palm region is associated with binding of these cognate ligands. Six domains of collective motions are identified, comprised of structurally close but not necessarily sequential residues. Two of these domains correspond to the flexible loops (L2 and L3), whereas the remaining domains correspond to the rigid part of the molecule.
Collapse
Affiliation(s)
- O Keskin
- Computational Technologies Laboratory, Screening Technologies Branch, Developmental Therapeutics Program, National Cancer Institute-Frederick, NIH, Frederick, Maryland, USA
| | | | | | | |
Collapse
|
24
|
Kurt N, Haliloglu T. Distribution of cooperative interactions in barnase at different time windows by coarse-grained simulations. POLYMER 2002. [DOI: 10.1016/s0032-3861(01)00431-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
25
|
Massi F, Straub JE. Probing the origins of increased activity of the E22Q "Dutch" mutant Alzheimer's beta-amyloid peptide. Biophys J 2001; 81:697-709. [PMID: 11463618 PMCID: PMC1301546 DOI: 10.1016/s0006-3495(01)75734-7] [Citation(s) in RCA: 38] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022] Open
Abstract
The amyloid peptide congener A beta(10--35)-NH(2) is simulated in an aqueous environment in both the wild type (WT) and E22Q "Dutch" mutant forms. The origin of the noted increase in deposition activity resulting from the Dutch mutation is investigated. Multiple nanosecond time scale molecular dynamics trajectories were performed and analyzed using a variety of measures of the peptide's average structure, hydration, conformational fluctuations, and dynamics. The results of the study support the conclusions that 1) the E22Q mutant and WT peptide are both stable in "collapsed coil" conformations consistent with the WT structure of, J. Struct. Biol. 130:130--141); 2) the E22Q peptide is more flexible in solution, supporting early claims that its equilibrium structural fluctuations are larger than those of the WT peptide; and 3) the local E22Q mutation leads to a change in the first solvation layer in the region of the peptide's "hydrophobic patch," resulting in a more hydrophobic solvation of the mutant peptide. The simulation results support the view that the noted increase in activity due to the Dutch mutation results from an enhancement of the desolvation process that is an essential step in the aggregation of the peptide.
Collapse
Affiliation(s)
- F Massi
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | | |
Collapse
|
26
|
Kurt N, Haliloğlu T. Conformational dynamics of subtilisin-chymotrypsin inhibitor 2 complex by coarse-grained simulations. J Biomol Struct Dyn 2001; 18:713-31. [PMID: 11334109 DOI: 10.1080/07391102.2001.10506702] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Abstract
An off-lattice dynamic Monte Carlo (MC) method is used to investigate the conformational dynamics of chymotrypsin inhibitor 2 (CI2) and subtilisin in both free and complex forms over two time windows, referring to short and long time scales. The conformational dynamics of backbone bonds analysed from several independent trajectories reveal that: Both the inhibitor and the enzyme are restricted in their bond rotations, excluding a few bonds, upon binding; the effect being greatest for the loop regions, and for the inhibitor. A cooperativity in the near-neighbor bond rotations are observed on both time scales, whereas the cooperative rotations of the bonds far along the sequence appear only in the long time window, and the latter time window is where most of the interactions between the inhibitor and the enzyme are observed. Upon binding, the cooperatively rotating parts of the inhibitor and the enzyme are readjusted compared to their free forms, and new correlations appear. The binding loop, although it is the closest contact region, is not the only part of the inhibitor involved in the interactions with the enzyme. Loops 3 and 8 and the helices F and G in bound enzyme and the binding loop of the inhibitor contribute at the most to the collective motions of whole structure on the slow time scale and are apparently important for enzyme-inhibitor interactions and function. The results in general provide evidence for the contribution of the loops with cooperative motions to the extensive communication network of the complex.
Collapse
Affiliation(s)
- N Kurt
- Polymer Research Center and Chemical Engineering Department, Boğaziçi University, Istanbul, Turkey
| | | |
Collapse
|
27
|
Higo J, Sugimoto Y, Wakabayashi K, Nakamura H. Collective motions of myosin head derived from backbone molecular dynamics and combination with X-ray solution scattering data. J Comput Chem 2001. [DOI: 10.1002/jcc.1147] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
28
|
|
29
|
Haliloglu T, Bahar I. Structure-based analysis of protein dynamics: comparison of theoretical results for hen lysozyme with X-ray diffraction and NMR relaxation data. Proteins 1999; 37:654-67. [PMID: 10651280 DOI: 10.1002/(sici)1097-0134(19991201)37:4<654::aid-prot15>3.0.co;2-j] [Citation(s) in RCA: 82] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
An analytical approach based on Gaussian network model (GNM) is proposed for predicting the rotational dynamics of proteins. The method, previously shown to successfully reproduce X-ray crystallographic temperature factors for a series of proteins is extended here to predict bond torsional mobilities and reorientation of main chain amide groups probed by 15N-H nuclear magnetic resonance (NMR) relaxation. The dynamics of hen egg-white lysozyme (HEWL) in the folded state is investigated using the proposed approach. Excellent agreement is observed between theoretical results and experimental (X-ray diffraction and NMR relaxation) data. The analysis reveals the important role of coupled rotations, or cross-correlations between dihedral angle librations, in defining the relaxation mechanism on a local scale. The crystal and solution structures exhibit some differences in their local motions, but their global motions are identical. Hinge residues mediating the cooperative movements of the alpha- and beta-domains are identified, which comprise residues in helix C, Glu35 and Ser36 on the loop succeeding helix B, Ile55 and Leu56 at the turn between strands II and III. The central part of the beta-domain long loop and the turn between strands I and II display an enhanced mobility. Finally, kinetically hot residues and key interactions are identified, which point at helix B and beta-strand III as the structural elements underlying the stability of the tertiary structure.
Collapse
Affiliation(s)
- T Haliloglu
- Polymer Research Center and Chemical Engineering Department, Bogazici University, Istanbul, Turkey
| | | |
Collapse
|
30
|
Abstract
A coarse-grained dynamic Monte Carlo (MC) simulation method is used to investigate the conformational dynamics of chymotrypsin inhibitor 2 (CI2). Each residue is represented therein by two interaction sites, one at the alpha-carbon and the other on the amino acid side-chain. The energy and geometry parameters extracted from databank structures are used. The calculated rms fluctuations of alpha-carbon atoms are in good agreement with crystallographic temperature factors. The two regions of the protein that pack against each other to form the main hydrophobic core exhibit negatively correlated fluctuations. The conformational dynamics could efficiently be probed by the time-delayed orientational and conformational correlation functions of the virtual bonds: the active site loop, excluding the active site bond, the turn region, and the N-terminal of the alpha-helix are relatively more mobile regions of the structure. A correlation is observed between the hydrogen/deuterium (H/D) exchange behavior and the long-time orientational and conformational autocorrelation function values for CI2. A cooperativity in the rotations of the bonds near in sequence is observed at all time windows, whereas the cooperative rotations of the bonds far along the sequence appear at long time windows; these correlations contribute to the stability of the secondary structures and the tertiary structure, respectively.
Collapse
Affiliation(s)
- N Kurt
- Chemical Engineering Department, Boğaziçi University, Bebek, Istanbul, Turkey
| | | |
Collapse
|
31
|
Lecomte JT, Sukits SF, Bhattacharya S, Falzone CJ. Conformational properties of native sperm whale apomyoglobin in solution. Protein Sci 1999; 8:1484-91. [PMID: 10422837 PMCID: PMC2144374 DOI: 10.1110/ps.8.7.1484] [Citation(s) in RCA: 46] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Abstract
Apomyoglobin from sperm whale is often used for studies of ligand binding, protein folding, and protein stability. In an effort to describe its conformational properties in solution, homonuclear and heteronuclear (13C and 15N) NMR methods were applied to the protein in its native state. Assignments were confirmed for nuclear Overhauser effects (NOEs) involving side chain and backbone protons in the folded regions of the structure. These NOEs were used to derive distance restraints. The shifts induced by the hydrophobic dye 8-anilino-1-naphthalenesulfonic acid (ANS) were inspected in the regions remote from its binding site and served as an indicator of conformational flexibility. 3JalphaH-NH values were obtained to assess dihedral angle averaging and to provide additional restraints. A family of structures was calculated with X-PLOR and an ab initio simulated annealing protocol using holomyoglobin as a template. Where the structure appeared well defined by chemical shift, line width, ANS perturbation, and density of NOEs, the low resolution model of apomyoglobin provides a valid approximation for the structure. The new model offers an improved representation of the folded regions of the protein, which encompass the A, B, E, helices as well as parts of the G and H helices. Regions that are less well defined at this stage of calculations include the CD corner and the end of the H-helix. The EF-F-FG segment remains uncharacterized.
Collapse
Affiliation(s)
- J T Lecomte
- Department of Chemistry and the Center for Biomolecular Structure and Function, the Pennsylvania State University, University Park 16802, USA.
| | | | | | | |
Collapse
|
32
|
Kikuchi T. Study of protein fluctuation with an effective inter-C? atomic potential derived from average distances between amino acids in proteins. J Comput Chem 1999; 20:713-719. [DOI: 10.1002/(sici)1096-987x(199905)20:7<713::aid-jcc6>3.0.co;2-s] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/1998] [Accepted: 01/08/1999] [Indexed: 11/09/2022]
|
33
|
|