1
|
Cassidy CK, Himes BA, Sun D, Ma J, Zhao G, Parkinson JS, Stansfeld PJ, Luthey-Schulten Z, Zhang P. Structure and dynamics of the E. coli chemotaxis core signaling complex by cryo-electron tomography and molecular simulations. Commun Biol 2020; 3:24. [PMID: 31925330 PMCID: PMC6954272 DOI: 10.1038/s42003-019-0748-0] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2019] [Accepted: 12/19/2019] [Indexed: 01/08/2023] Open
Abstract
To enable the processing of chemical gradients, chemotactic bacteria possess large arrays of transmembrane chemoreceptors, the histidine kinase CheA, and the adaptor protein CheW, organized as coupled core-signaling units (CSU). Despite decades of study, important questions surrounding the molecular mechanisms of sensory signal transduction remain unresolved, owing especially to the lack of a high-resolution CSU structure. Here, we use cryo-electron tomography and sub-tomogram averaging to determine a structure of the Escherichia coli CSU at sub-nanometer resolution. Based on our experimental data, we use molecular simulations to construct an atomistic model of the CSU, enabling a detailed characterization of CheA conformational dynamics in its native structural context. We identify multiple, distinct conformations of the critical P4 domain as well as asymmetries in the localization of the P3 bundle, offering several novel insights into the CheA signaling mechanism.
Collapse
Affiliation(s)
- C Keith Cassidy
- Department of Biochemistry, University of Oxford, Oxford, OX1 3QU, UK.
- Department of Physics and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA.
| | - Benjamin A Himes
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, 15260, USA
| | - Dapeng Sun
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, 15260, USA
| | - Jun Ma
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, 15260, USA
| | - Gongpu Zhao
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, 15260, USA
| | - John S Parkinson
- School of Biological Sciences, University of Utah, Salt Lake City, UT, 84112, USA
| | - Phillip J Stansfeld
- Department of Biochemistry, University of Oxford, Oxford, OX1 3QU, UK
- School of Life Sciences & Department of Chemistry, University of Warwick, Coventry, CV4 7AL, UK
| | - Zaida Luthey-Schulten
- Department of Physics and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA
- Department of Chemistry and Center for the Physics of Living Cells, University of Illinois at Urbana-Champaign, Urbana, IL, 61801, USA
| | - Peijun Zhang
- Department of Structural Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA, 15260, USA.
- Division of Structural Biology, Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, OX3 7BN, UK.
- Electron Bio-Imaging Centre, Diamond Light Sources, Harwell Science and Innovation Campus, Didcot, OX11 0DE, UK.
| |
Collapse
|
2
|
Nanjundaiah S, Mutturi S, Bhatt P. Modeling of caffeine degradation kinetics during cultivation of Fusarium solani using sucrose as co-substrate. Biochem Eng J 2017. [DOI: 10.1016/j.bej.2017.05.018] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
3
|
Yang B, Lin Z. Systematic search of conformations of five tetrapeptides and a divide and conquer strategy for the predictions of peptide structures. COMPUT THEOR CHEM 2017. [DOI: 10.1016/j.comptc.2017.03.014] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
4
|
Protein remains stable at unusually high temperatures when solvated in aqueous mixtures of amino acid based ionic liquids. J Mol Model 2016; 22:258. [DOI: 10.1007/s00894-016-3123-9] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2016] [Accepted: 09/08/2016] [Indexed: 10/20/2022]
|
5
|
Bernardi RC, Melo MCR, Schulten K. Enhanced sampling techniques in molecular dynamics simulations of biological systems. Biochim Biophys Acta Gen Subj 2014; 1850:872-877. [PMID: 25450171 DOI: 10.1016/j.bbagen.2014.10.019] [Citation(s) in RCA: 410] [Impact Index Per Article: 41.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2014] [Revised: 10/17/2014] [Accepted: 10/20/2014] [Indexed: 01/29/2023]
Abstract
BACKGROUND Molecular dynamics has emerged as an important research methodology covering systems to the level of millions of atoms. However, insufficient sampling often limits its application. The limitation is due to rough energy landscapes, with many local minima separated by high-energy barriers, which govern the biomolecular motion. SCOPE OF REVIEW In the past few decades methods have been developed that address the sampling problem, such as replica-exchange molecular dynamics, metadynamics and simulated annealing. Here we present an overview over theses sampling methods in an attempt to shed light on which should be selected depending on the type of system property studied. MAJOR CONCLUSIONS Enhanced sampling methods have been employed for a broad range of biological systems and the choice of a suitable method is connected to biological and physical characteristics of the system, in particular system size. While metadynamics and replica-exchange molecular dynamics are the most adopted sampling methods to study biomolecular dynamics, simulated annealing is well suited to characterize very flexible systems. The use of annealing methods for a long time was restricted to simulation of small proteins; however, a variant of the method, generalized simulated annealing, can be employed at a relatively low computational cost to large macromolecular complexes. GENERAL SIGNIFICANCE Molecular dynamics trajectories frequently do not reach all relevant conformational substates, for example those connected with biological function, a problem that can be addressed by employing enhanced sampling algorithms. This article is part of a Special Issue entitled Recent developments of molecular dynamics.
Collapse
Affiliation(s)
| | - Marcelo C R Melo
- Center for Biophysics and Computational Biology, University of Illinois, Urbana, IL 61801, USA
| | - Klaus Schulten
- Beckman Institute, University of Illinois, Urbana, IL 61801, USA; Center for Biophysics and Computational Biology, University of Illinois, Urbana, IL 61801, USA; Department of Physics, University of Illinois, Urbana, IL 61801, USA.
| |
Collapse
|
6
|
Chaotic multiquenching annealing applied to the protein folding problem. ScientificWorldJournal 2014; 2014:364352. [PMID: 24790563 PMCID: PMC3980881 DOI: 10.1155/2014/364352] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2013] [Accepted: 01/19/2014] [Indexed: 11/17/2022] Open
Abstract
The Chaotic Multiquenching Annealing algorithm (CMQA) is proposed. CMQA is a new algorithm, which is applied to protein folding problem (PFP). This algorithm is divided into three phases: (i) multiquenching phase (MQP), (ii) annealing phase (AP), and (iii) dynamical equilibrium phase (DEP). MQP enforces several stages of quick quenching processes that include chaotic functions. The chaotic functions can increase the exploration potential of solutions space of PFP. AP phase implements a simulated annealing algorithm (SA) with an exponential cooling function. MQP and AP are delimited by different ranges of temperatures; MQP is applied for a range of temperatures which goes from extremely high values to very high values; AP searches for solutions in a range of temperatures from high values to extremely low values. DEP phase finds the equilibrium in a dynamic way by applying least squares method. CMQA is tested with several instances of PFP.
Collapse
|
7
|
Melo MCR, Bernardi RC, Fernandes TVA, Pascutti PG. GSAFold: a new application of GSA to protein structure prediction. Proteins 2012; 80:2305-10. [PMID: 22622959 DOI: 10.1002/prot.24120] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2012] [Revised: 05/08/2012] [Accepted: 05/21/2012] [Indexed: 11/07/2022]
Abstract
The folding process defines three-dimensional protein structures from their amino acid chains. A protein's structure determines its activity and properties; thus knowing such conformation on an atomic level is essential for both basic and applied studies of protein function and dynamics. However, the acquisition of such structures by experimental methods is slow and expensive, and current computational methods mostly depend on previously known structures to determine new ones. Here we present a new software called GSAFold that applies the generalized simulated annealing (GSA) algorithm on ab initio protein structure prediction. The GSA is a stochastic search algorithm employed in energy minimization and used in global optimization problems, especially those that depend on long-range interactions, such as gravity models and conformation optimization of small molecules. This new implementation applies, for the first time in ab initio protein structure prediction, an analytical inverse for the Visitation function of GSA. It also employs the broadly used NAMD Molecular Dynamics package to carry out energy calculations, allowing the user to select different force fields and parameterizations. Moreover, the software also allows the execution of several simulations simultaneously. Applications that depend on protein structures include rational drug design and structure-based protein function prediction. Applying GSAFold in a test peptide, it was possible to predict the structure of mastoparan-X to a root mean square deviation of 3.00 Å.
Collapse
Affiliation(s)
- Marcelo C R Melo
- Laboratório de Biotecnologia/DIPRO, Instituto Nacional de Metrologia, Qualidade e Tecnologia, Rio de Janeiro, Brasil
| | | | | | | |
Collapse
|
8
|
Pandey PR, Roy S. Distinctions in early stage unwinding mechanisms of zwitterionic, capped, and neutral forms of different α-helical homopolymeric peptides. J Phys Chem B 2012; 116:4731-40. [PMID: 22448707 DOI: 10.1021/jp301556x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Molecular dynamics simulations of α-helical polyalanine, polyleucine, polylysine, and poly(glutamic acid) with different forms of terminal groups in water at 300 K showed sharp distinctions in their unwinding mechanisms. Zwitterionic, capped, and neutral forms of polyalanine, polyleucine, and polylysine have been explored to elucidate their unwinding mechanism at very early stage, e.g., initial time window. Role of water in the unwinding mechanisms of the various helices has been envisaged. Also, it is evident from our calculations that the short- and long-range nonbonded interactions among the side chains is an important factor determining the unwinding mechanisms of the various homopolymeric α-helices. These findings can be helpful in constructing predictive models for understanding of the unwinding of α-helical proteins and peptides.
Collapse
Affiliation(s)
- Prithvi Raj Pandey
- Physical Chemistry Division, National Chemical Laboratory , Pune 411008, India
| | | |
Collapse
|
9
|
Balamurugan K, Gopalakrishnan R, Raman SS, Subramanian V. Exploring the Changes in the Structure of α-Helical Peptides Adsorbed onto a Single Walled Carbon Nanotube Using Classical Molecular Dynamics Simulation. J Phys Chem B 2010; 114:14048-58. [DOI: 10.1021/jp106177n] [Citation(s) in RCA: 51] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- K. Balamurugan
- Chemical Laboratory, Central Leather Research Institute, Council of Scientific and Industrial Research, Adyar, Chennai- 600020, India
| | - R. Gopalakrishnan
- Chemical Laboratory, Central Leather Research Institute, Council of Scientific and Industrial Research, Adyar, Chennai- 600020, India
| | - S. Sundar Raman
- Chemical Laboratory, Central Leather Research Institute, Council of Scientific and Industrial Research, Adyar, Chennai- 600020, India
| | - V. Subramanian
- Chemical Laboratory, Central Leather Research Institute, Council of Scientific and Industrial Research, Adyar, Chennai- 600020, India
| |
Collapse
|
10
|
Crokidakis N, Soares-Pinto DO, Reis MS, Souza AM, Sarthour RS, Oliveira IS. Finite-size analysis of a two-dimensional Ising model within a nonextensive approach. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:051101. [PMID: 20364941 DOI: 10.1103/physreve.80.051101] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/07/2009] [Revised: 08/28/2009] [Indexed: 05/29/2023]
Abstract
In this work we present a thorough analysis of the phase transitions that occur in a ferromagnetic two-dimensional Ising model, with only nearest-neighbors interactions, in the framework of the Tsallis nonextensive statistics. We performed Monte Carlo simulations on square lattices with linear sizes L ranging from 32 up to 512. The statistical weight of the Metropolis algorithm was changed according to the nonextensive statistics. Discontinuities in the m(T) curve are observed for q< or =0.5 . However, we have verified only one peak on the energy histograms at the critical temperatures, indicating the occurrence of continuous phase transitions. For the 0.5<q< or =1.0 regime, we have found continuous phase transitions between the ordered and the disordered phases, and determined the critical exponents via finite-size scaling. We verified that the critical exponents alpha , beta , and gamma depend on the entropic index q in the range 0.5<q< or =1.0 in the form alpha(q)=(10q2-33q+23)/20 , beta(q)=(2q-1)/8 , and gamma(q)=(q2-q+7)/4 . On the other hand, the critical exponent nu does not depend on q . This suggests a violation of the scaling relations 2beta+gamma=dnu and alpha+2beta+gamma=2 and a nonuniversality of the critical exponents along the ferro-paramagnetic frontier.
Collapse
Affiliation(s)
- N Crokidakis
- Instituto de Física, Universidade Federal Fluminense, Av. Litorânea s/n, 24210-340 Niterói, RJ, Brazil.
| | | | | | | | | | | |
Collapse
|
11
|
Moret MA, Santana MC, Zebende GF, Pascutti PG. Self-similarity and protein compactness. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2009; 80:041908. [PMID: 19905343 DOI: 10.1103/physreve.80.041908] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/13/2009] [Revised: 08/14/2009] [Indexed: 05/28/2023]
Abstract
The hydrophobic effect is the major factor that drives a protein toward collapse and folding. As a consequence of the folding process a hydrophobic core is shielded by the solvent-accessible surface area of the protein. We analyze the solvent-accessible surface area of 1825 nonhomolog protein chains deposited in the Brookhaven Protein Data Bank. This solvent-accessible surface area presents an intrinsic self-similarity behavior. The comparison between the accessible surface area as function of the number of amino acids and the accessible surface area as function of gyration radius supplies a measure of the scaling exponent close to the one observed by volume as function of radius of gyration or by mass-size exponent. The present finding indicates that the fractal analysis describes the protein compactness as an object packing between random spheres in percolation threshold and crumpled wires.
Collapse
Affiliation(s)
- M A Moret
- Programa de Modelagem Computacional, SENAI CIMATEC, Salvador, BA, Brazil.
| | | | | | | |
Collapse
|
12
|
Abstract
The development of affordable, high-throughput sequencing technology has led to a flood of publicly available bacterial genome-sequence data. The availability of multiple genome sequences presents both an opportunity and a challenge for microbiologists, and new computational approaches are needed to extract the knowledge that is required to address specific biological problems and to analyse genomic data. The field of e-Science is maturing, and Grid-based technologies can help address this challenge.
Collapse
|
13
|
Raman SS, Vijayaraj R, Parthasarathi R, Subramanian V. Helix forming tendency of valine substituted poly-alanine: a molecular dynamics investigation. J Phys Chem B 2008; 112:9100-4. [PMID: 18597521 DOI: 10.1021/jp7119813] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
In this study, classical molecular dynamics simulations have been carried out on the valine (guest) substituted poly alanine (host) using the host-guest peptide approach to understand the role of valine in the formation and stabilization of helix. Valine has been substituted in the host peptide starting from N terminal to C terminal. Various structural parameters have been obtained from the molecular dynamics simulation to understand the tolerance of helical motif to valine. Depending on the position of valine in the host peptide, it stabilizes (or destabilizes) the formation of the helical structure. The substitution of valine in the poly alanine at some positions has no effect on the helix formation (deformation). It is interesting to observe the coexistence of 3 10 and alpha-helix in the peptides due to the dynamical nature of the hydrogen bonding interaction and sterical interactions.
Collapse
Affiliation(s)
- S Sundar Raman
- Chemical Laboratory, Central Leather Research Institute, Adyar, Chennai 600 020 India
| | | | | | | |
Collapse
|
14
|
Moon JH, Sun M, Kim MS. Efficient and reliable calculation of Rice-Ramsperger-Kassel-Marcus unimolecular reaction rate constants for biopolymers: modification of Beyer-Swinehart algorithm for degenerate vibrations. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2007; 18:1063-9. [PMID: 17448674 DOI: 10.1016/j.jasms.2007.03.012] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/20/2007] [Revised: 03/23/2007] [Accepted: 03/26/2007] [Indexed: 05/15/2023]
Abstract
The Beyer-Swinehart (BS) algorithm, which calculates vibrational state density and sum, was modified for simultaneous treatment of degenerate vibrations. The modified algorithm was used in the grouped-frequency mode of the Rice-Ramsperger-Kassel-Marcus (RRKM) unimolecular reaction rate constant calculation for proteins with relative molecular mass as large as 100,000. Compared to the original BS method, reduction in computation time by a factor of around 3000 was achieved. Even though large systematic errors arising from frequency grouping were observed for state densities and sums, they more or less canceled each other, thus enabling reliable rate constant calculation. The present method is thought to be adequate for efficient and reliable RRKM calculations for any macromolecule in the gas phase such as the molecular ions of proteins, nucleic acids, and carbohydrates generated inside a mass spectrometer. The algorithm can also be used to calculate the internal energy distribution of a macromolecule at thermal equilibrium.
Collapse
Affiliation(s)
- Jeong Hee Moon
- Systemic Proteomics Research Center, KRIBB, Daejeon, Korea
| | | | | |
Collapse
|
15
|
Mikhonin AV, Asher SA, Bykov SV, Murza A. UV Raman spatially resolved melting dynamics of isotopically labeled polyalanyl peptide: slow alpha-helix melting follows 3(10)-helices and pi-bulges premelting. J Phys Chem B 2007; 111:3280-92. [PMID: 17388440 DOI: 10.1021/jp0654009] [Citation(s) in RCA: 30] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
We used UV resonance Raman (UVRR) to examine the spatial dependence of the T-jump secondary structure relaxation of an isotopically labeled 21-residue mainly Ala peptide, AdP. The AdP penultimate Ala residues were perdeuterated, leaving the central residues hydrogenated, to allow separate monitoring of melting of the middle versus the end peptide bonds. For 5 to 30 degrees C T-jumps, the central peptide bonds show a approximately 2-fold slower relaxation time (189 +/- 31 ns) than do the exterior peptide bonds (97 +/- 15 ns). In contrast, for a 20 to 40 degrees C T-jump, the central peptide bond relaxation appears to be faster (56 +/- 6 ns) than that of the penultimate peptide bonds (131 +/- 46 ns). We show that, if the data are modeled as a two-state transition, we find that only exterior peptide bonds show anti-Arrhenius folding behavior; the middle peptide bonds show both normal Arrhenius-like folding and unfolding. This anti-Arrhenius behavior results from the involvement of pi-bulges/helices and 3(10)-helix states in the melting. The unusual temperature dependence of the (un)folding rates of the interior and exterior peptide bonds is due to the different relative (un)folding rates of 3(10)-helices, alpha-helices, and pi-bulges/helices. Pure alpha-helix unfolding rates are approximately 12-fold slower (approximately 1 micros) than that of pi-bulges and 3(10)-helices. In addition, we also find that the alpha-helix is most stable at the AdP N-terminus where eight consecutive Ala occur, whereas the three hydrophilic Arg located in the middle and at the C-terminus destabilize the alpha-helix in these regions and induce defects such as pi-bulges and 3(10)-helices.
Collapse
Affiliation(s)
- Aleksandr V Mikhonin
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | | | | |
Collapse
|
16
|
Moret MA, Zebende GF. Amino acid hydrophobicity and accessible surface area. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2007; 75:011920. [PMID: 17358197 DOI: 10.1103/physreve.75.011920] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/02/2005] [Revised: 10/20/2006] [Indexed: 05/14/2023]
Abstract
It is well known that the hydrophobic effect is the major factor that drives a protein toward collapse and folding. We analyze the variation of the solvent-accessible surface area of amino acids in small fragments of protein (3N45) . In this way, we look into 5526 protein chains deposited in the Brookhaven Protein Data Bank. The accessible surface area behaves as a power law for N9 . The comparison between the loss of accessible area and the self-similar behavior gives us a measure of the possibility of an amino acid to have apolar or polar side chain. It is therefore possible to infer about amino acid hydrophobicity, i.e., if one amino acid has a hydrophobic side chain or if it has a hydrophilic one. Furthermore, the present findings indicate that the variation of the accessible surface area describes an alternative hydrophobicity scale.
Collapse
Affiliation(s)
- M A Moret
- Programa de Modelagem Computacional-Fundação Visconde de Cairu, 40226-900 Salvador, Bahia, Brazil.
| | | |
Collapse
|