-
From Feynman's ratchet to timecrystalline molecular motors
Authors:
Jianmei Wang,
Jin Dai,
Antti J. Niemi,
Xubiao Peng
Abstract:
Cats use the connection governing parallel transport in the space of shapes to land safely on their feet. Here we argue that this connection also explains the impressive performance of molecular motors by enabling molecules to evade conclusions of Feynman's ratchet-and-pawl analysis. We first demonstrate, using simple molecular models, how directed rotational motion can emerge from shape changes e…
▽ More
Cats use the connection governing parallel transport in the space of shapes to land safely on their feet. Here we argue that this connection also explains the impressive performance of molecular motors by enabling molecules to evade conclusions of Feynman's ratchet-and-pawl analysis. We first demonstrate, using simple molecular models, how directed rotational motion can emerge from shape changes even without angular momentum. We then computationally design knotted polyalanine molecules and show how their shape space connection organizes individual atom thermal vibrations into collective rotational motion, independently of angular momentum. Our simulations show that rotational motion arises effortlessly even in ambient water, making the molecule an effective theory time crystal. Our findings have potential for practical molecular motor design and engineering and can be verified through high-precision nuclear magnetic resonance measurements.
△ Less
Submitted 24 April, 2023;
originally announced April 2023.
-
Topological Indices of Proteins
Authors:
Dmitry Melnikov,
Antti J Niemi,
Ara Sedrakyan
Abstract:
Protein molecules can be approximated by discrete polygonal chains of amino acids. Standard topological tools can be applied to the smoothening of the polygons to introduce a topological classification of proteins, for example, using the self-linking number of the corresponding framed curves. In this paper we add new details to the standard classification. Known definitions of the self-linking num…
▽ More
Protein molecules can be approximated by discrete polygonal chains of amino acids. Standard topological tools can be applied to the smoothening of the polygons to introduce a topological classification of proteins, for example, using the self-linking number of the corresponding framed curves. In this paper we add new details to the standard classification. Known definitions of the self-linking number apply to non-singular framings: for example, the Frenet framing cannot be used if the curve has inflection points. Meanwhile in the discrete proteins the special points are naturally resolved. Consequently, a separate integer topological characteristics can be introduced, which takes into account the intrinsic features of the special points. For large number of proteins we compute integer topological indices associated with the singularities of the Frenet framing. We show how a version of the Calugareanu's theorem is satisfied for the associated self-linking number of a discrete curve. Since the singularities of the Frenet framing correspond to the structural motifs of proteins, we propose topological indices as a technical tool for the description of the folding dynamics of proteins.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
Can all-atom protein dynamics be reconstructed from the knowledge of C-alpha time evolution?
Authors:
Jiaojiao Liu,
Jin Dai,
Jianfeng He,
Xubiao Peng,
Antti J. Niemi
Abstract:
We inquire to what extent protein peptide plane and side chain dynamics can be reconstructed from knowledge of C-alpha dynamics. Due to lack of experimental data we analyze all atom molecular dynamics trajectories from Anton supercomputer, and for clarity we limit our attention to the peptide plane O atoms and side chain C-beta atoms. We try and reconstruct their dynamics using four different appr…
▽ More
We inquire to what extent protein peptide plane and side chain dynamics can be reconstructed from knowledge of C-alpha dynamics. Due to lack of experimental data we analyze all atom molecular dynamics trajectories from Anton supercomputer, and for clarity we limit our attention to the peptide plane O atoms and side chain C-beta atoms. We try and reconstruct their dynamics using four different approaches. Three of these are the publicly available reconstruction programs Pulchra, Remo Scwrl4. The fourth, Statistical Method, builds entirely on statistical analysis of Protein Data Bank (PDB) structures. All four methods place the O and C-beta atoms accurately along the Anton trajectories. However, the Statistical Method performs best. The results suggest that under physiological conditions, the all atom dynamics is slaved to that of C-alpha atoms. The results can help improve all atom force fields, and advance reconstruction and refinement methods for reduced protein structures. The results provide impetus for development of effective coarse grained force fields in terms of reduced coordinates.
△ Less
Submitted 21 January, 2019;
originally announced January 2019.
-
Virtual reality analysis of intrinsic protein geometry with applications to cis peptide planes
Authors:
Yanzhen Hou,
Jin Dai,
Nevena Ilieva,
Antti J. Niemi,
Xubiao Peng,
Jianfeng He
Abstract:
A protein is traditionally visualised as a piecewise linear discrete curve, and its geometry is conventionally characterised by the extrinsically determined Ramachandran angles. However, a protein backbone has also two independent intrinsic geometric structures, due to the peptide planes and the side chains. Here we adapt and develop modern 3D virtual reality techniques to scrutinize the atomic ge…
▽ More
A protein is traditionally visualised as a piecewise linear discrete curve, and its geometry is conventionally characterised by the extrinsically determined Ramachandran angles. However, a protein backbone has also two independent intrinsic geometric structures, due to the peptide planes and the side chains. Here we adapt and develop modern 3D virtual reality techniques to scrutinize the atomic geometry along a protein backbone, in the vicinity of a peptide plane. For this we compare backbone geometry-based (extrinsic) and structure-based (intrinsic) coordinate systems, and as an example we inspect the trans and cis peptide planes. We reveal systematics in the way how a cis peptide plane deforms the neighbouring atomic geometry, and we develop a virtual reality based visual methodology that can identify the presence of a cis peptide plane from the arrangement of atoms in its vicinity. Our approach can easily detect exceptionally placed atoms in crystallographic structures. Thus it can be employed as a powerful visual refinement tool which is applicable also in the case when resolution of the protein structure is limited and whenever refinement is needed. As concrete examples we identify a number of crystallographic protein structures in Protein Data Bank (PDB) that display exceptional atomic positions around their cis peptide planes.
△ Less
Submitted 6 June, 2017; v1 submitted 5 June, 2017;
originally announced June 2017.
-
Multiple scales and phases in discrete chains with application to folded proteins
Authors:
Anna Sinelnikova,
Antti J. Niemi,
Johan Nilsson,
Maksim Ulybyshev
Abstract:
Chiral heteropolymers such as larger globular proteins can simultaneously support multiple length scales. The interplay between different scales brings about conformational diversity, and governs the structure of the energy landscape. Multiple scales produces also complex dynamics, which in the case of proteins sustains live matter. However, thus far no clear understanding exist, how to distinguis…
▽ More
Chiral heteropolymers such as larger globular proteins can simultaneously support multiple length scales. The interplay between different scales brings about conformational diversity, and governs the structure of the energy landscape. Multiple scales produces also complex dynamics, which in the case of proteins sustains live matter. However, thus far no clear understanding exist, how to distinguish the various scales that determine the structure and dynamics of a complex protein. Here we propose a systematic method to identify the scales in chiral heteropolymers such as a protein. For this we introduce a novel order parameter, that not only reveals the scales but also probes the phase structure. In particular, we argue that a chiral heteropolymer can simultaneously display traits of several different phases, contingent on the length scale at which it is scrutinized. Our approach builds on a variant of Kadanoff's block-spin transformation that we employ to coarse grain piecewise linear chains such as the C$α$ backbone of a protein. We derive analytically and then verify numerically a number of properties that the order parameter can display. We demonstrate how, in the case of crystallographic protein structures in Protein Data Bank, the order parameter reveals the presence of different length scales, and we propose that a relation must exist between the scales, phases, and the complexity of folding pathways.
△ Less
Submitted 26 May, 2017;
originally announced May 2017.
-
Towards multistage modelling of protein dynamics with monomeric Myc oncoprotein as an example
Authors:
Jiaojiao Liu,
Jin Dai,
Jianfeng He,
Antti J. Niemi,
Nevena Ilieva
Abstract:
We propose to combine a mean field approach with all atom molecular dynamics into a multistage algorithm that can model protein folding and dynamics over very long time periods yet with atomic level precision. As an example we investigate an isolated monomeric Myc oncoprotein that has been implicated in carcinomas including those in colon, breast and lungs. Under physiological conditions a monomer…
▽ More
We propose to combine a mean field approach with all atom molecular dynamics into a multistage algorithm that can model protein folding and dynamics over very long time periods yet with atomic level precision. As an example we investigate an isolated monomeric Myc oncoprotein that has been implicated in carcinomas including those in colon, breast and lungs. Under physiological conditions a monomeric Myc is presumed to be an example of intrinsically disordered proteins, that pose a serious challenge to existing modelling techniques. We argue that a room temperature monomeric Myc is in a dynamical state, it oscillates between different conformations that we identify. For this we adopt the C-alpha backbone of Myc in a crystallographic heteromer as an initial Ansatz for the monomeric structure. We construct a multisoliton of the pertinent Landau free energy, to describe the C-alpha profile with ultra high precision. We use Glauber dynamics to resolve how the multisoliton responds to repeated increases and decreases in ambient temperature. We confirm that the initial structure is unstable in isolation. We reveal a highly degenerate ground state landscape, an attractive set towards which Glauber dynamics converges in the limit of vanishing ambient temperature. We analyse the thermal stability of this Glauber attractor using room temperature molecular dynamics. We identify and scrutinise a particularly stable subset in which the two helical segments of the original multisoliton align in parallel, next to each other. During the MD time evolution of a representative structure from this subset, we observe intermittent quasiparticle oscillations along the C-terminal alpha-helix, some of which resemble a translating Davydov's Amide-I soliton. We propose that the presence of oscillatory motion is in line with the expected intrinsically disordered character of Myc.
△ Less
Submitted 5 December, 2016;
originally announced December 2016.
-
Solution X-ray scattering (S/WAXS) and structure formation in protein dynamics
Authors:
Alexandr Nasedkin,
Jan Davidsson,
Antti J. Niemi,
Xubiao Peng
Abstract:
We propose to develop mean field theory in combination with Glauber algorithm, to model and interpret protein dynamics and structure formation in small to wide angle x-ray scattering (S/WAXS) experiments. We develop the methodology by analysing the Engrailed homeodomain protein as an example. We demonstrate how to interpret S/WAXS data with a good precision and over an extended temperature range.…
▽ More
We propose to develop mean field theory in combination with Glauber algorithm, to model and interpret protein dynamics and structure formation in small to wide angle x-ray scattering (S/WAXS) experiments. We develop the methodology by analysing the Engrailed homeodomain protein as an example. We demonstrate how to interpret S/WAXS data with a good precision and over an extended temperature range. We explain experimentally observed phenomena in terms of protein phase structure, and we make predictions for future experiments how the scattering data behaves at different ambient temperature values. We conclude that a combination of mean field theory with Glauber algorithm has the potential to develop into a highly accurate, computationally effective and predictive tool for analysing S/WAXS data. Finally, we compare our results with those obtained previously in an all-atom molecular dynamics simulation.
△ Less
Submitted 5 June, 2017; v1 submitted 24 November, 2016;
originally announced November 2016.
-
Thermal unfolding of myoglobin in the Landau-Ginzburg-Wilson approach
Authors:
Xubiao Peng,
Adam Sieradzan,
Antti J. Niemi
Abstract:
The Landau-Ginzburg-Wilson paradigm is applied to model the low-temperature crystallographic C$α$ backbone structure of sperm whale myoglobin. The Glauber protocol is employed to simulate its response to an increase in ambient temperature. The myoglobin is found to unfold from its native state by a succession of $α$-helical intermediates, fully in line with the observed folding and unfolding patte…
▽ More
The Landau-Ginzburg-Wilson paradigm is applied to model the low-temperature crystallographic C$α$ backbone structure of sperm whale myoglobin. The Glauber protocol is employed to simulate its response to an increase in ambient temperature. The myoglobin is found to unfold from its native state by a succession of $α$-helical intermediates, fully in line with the observed folding and unfolding patterns in denaturation experiments. In particular, a molten globule intermediate is identified with experimentally correct attributes. A detailed, experimentally testable contact map is constructed to characterise the specifics of the unfolding pathway, including the formation of long range interactions. The results reveal how the unfolding process of a protein is driven by the interplay between, and a successive melting of, its modular secondary structure components.
△ Less
Submitted 7 February, 2016;
originally announced February 2016.
-
Bloch spin waves and emergent structure in protein folding with HIV envelope glycoprotein as an example
Authors:
Jin Dai,
Antti J. Niemi,
Jianfeng He,
Adam Sieradzan,
Nevena Ilieva
Abstract:
We inquire how structure emerges during the process of protein folding. For this we scrutinise col- lective many-atom motions during all-atom molecular dynamics simulations. We introduce, develop and employ various topological techniques, in combination with analytic tools that we deduce from the concept of integrable models and structure of discrete nonlinear Schroedinger equation. The example we…
▽ More
We inquire how structure emerges during the process of protein folding. For this we scrutinise col- lective many-atom motions during all-atom molecular dynamics simulations. We introduce, develop and employ various topological techniques, in combination with analytic tools that we deduce from the concept of integrable models and structure of discrete nonlinear Schroedinger equation. The example we consider is an alpha-helical subunit of the HIV envelope glycoprotein gp41. The helical structure is stable when the subunit is part of the biological oligomer. But in isolation the helix becomes unstable, and the monomer starts deforming. We follow the process computationally. We interpret the evolving structure both in terms of a backbone based Heisenberg spin chain and in terms of a side chain based XY spin chain. We find that in both cases the formation of protein super-secondary structure is akin the formation of a topological Bloch domain wall along a spin chain. During the process we identify three individual Bloch walls and we show that each of them can be modelled with a very high precision in terms of a soliton solution to a discrete nonlinear Schroedinger equation.
△ Less
Submitted 23 November, 2015;
originally announced November 2015.
-
On phase diagram and the pseudogap state in a linear chiral homopolymer model
Authors:
A. Sinelnikova,
A. J. Niemi,
M. Ulybyshev
Abstract:
The phase structure of a homopolymer chain is investigated in terms of a universal theoretical model, designed to describe the infrared limit of slow spatial variations. The effects of chirality are studied and compared with the influence of a short-range attractive interaction between monomers, at various ambient temperature values. In the high temperature limit the homopolymer chain is in the se…
▽ More
The phase structure of a homopolymer chain is investigated in terms of a universal theoretical model, designed to describe the infrared limit of slow spatial variations. The effects of chirality are studied and compared with the influence of a short-range attractive interaction between monomers, at various ambient temperature values. In the high temperature limit the homopolymer chain is in the self-avoiding random walk phase. But at low temperatures, two different phases are possible: When short-range attractive interactions dominate over chirality, the chain collapses into a space- filling conformation. But when the attractive interactions become weaker, there is a low temperature unfolding transition and the chain becomes like a straight rod. Between the high temperature and low temperature limits, several intermediate states are observed. For sufficiently high values of short-range attraction, the conventional θ-regime is observed between the self-avoiding random walk phase and the space filling collapsed phase. But when chirality increases, there is a trasition from the θ-regime to a pseudogap state. Moreover, a regime akin the θ-regime is identified between the pseudogap state, and the low temperature phase where the chain is like a straight rod. Applications to polymers and proteins, in particular collagen, are suggested.
△ Less
Submitted 23 December, 2015; v1 submitted 21 January, 2015;
originally announced January 2015.
-
WHAT IS LIFE - Sub-cellular Physics of Live Matter
Authors:
Antti J. Niemi
Abstract:
This is a set of lectures that I presented at the Les Houches 2014 Summer School "Topological Aspects in Condensed Matter Physics". The lectures are an introduction to physics of proteins. To physicists, and by a physicist. My lectures at les Houches were also celebration of the anniversary of Schroedinger's 1944 lectures, and for that reason I decided to share my title with his book.
This is a set of lectures that I presented at the Les Houches 2014 Summer School "Topological Aspects in Condensed Matter Physics". The lectures are an introduction to physics of proteins. To physicists, and by a physicist. My lectures at les Houches were also celebration of the anniversary of Schroedinger's 1944 lectures, and for that reason I decided to share my title with his book.
△ Less
Submitted 29 December, 2014;
originally announced December 2014.
-
Virtual reality based approach to protein heavy-atom structure reconstruction
Authors:
Xubiao Peng,
Alireza Chenani,
Shuangwei Hu,
Yifan Zhou,
Antti J. Niemi
Abstract:
A commonly recurring problem in structural protein studies, is the determination of all heavy atom positions from the knowledge of the central alpha-carbon coordinates. We employ advances in virtual reality to address the problem. The outcome is a 3D visualisation based technique where all the heavy backbone and side chain atoms are treated on equal footing, in terms of the C-alpha coordinates. Ea…
▽ More
A commonly recurring problem in structural protein studies, is the determination of all heavy atom positions from the knowledge of the central alpha-carbon coordinates. We employ advances in virtual reality to address the problem. The outcome is a 3D visualisation based technique where all the heavy backbone and side chain atoms are treated on equal footing, in terms of the C-alpha coordinates. Each heavy atom can be visualised on the surfaces of the different two-spheres, that are centered at the other heavy backbone and side chain atoms. In particular, the rotamers are visible as clusters which display strong dependence on the underlying backbone secondary structure. Our method easily detects those atoms in a crystallographic protein structure which have been been likely misplaced. Our approach forms a basis for the development of a new generation, visualisation based side chain construction, validation and refinement tools. The heavy atom positions are identified in a manner which accounts for the secondary structure environment, leading to improved accuracy over existing methods.
△ Less
Submitted 26 December, 2014;
originally announced December 2014.
-
Aspects of structural landscape of human islet amyloid polypeptide
Authors:
Jianfeng He,
Jin Dai,
Jing Li,
Xubiao Peng,
Antti J. Niemi
Abstract:
The human islet amyloid polypeptide (hIAPP) co-operates with insulin to maintain glycemic balance. It also constitutes the amyloid plaques that aggregate in the pancreas of type-II diabetic patients. We have performed extensive in silico investigations to analyse the structural landscape of monomeric hIAPP, which is presumed to be intrinsically disordered. For this we construct from first principl…
▽ More
The human islet amyloid polypeptide (hIAPP) co-operates with insulin to maintain glycemic balance. It also constitutes the amyloid plaques that aggregate in the pancreas of type-II diabetic patients. We have performed extensive in silico investigations to analyse the structural landscape of monomeric hIAPP, which is presumed to be intrinsically disordered. For this we construct from first principles a highly predictive energy function that describes a monomeric hIAPP observed in a NMR experiment, as a local energy minimum. We subject our theoretical model of hIAPP to repeated heating and cooling simulations, back and forth between a high temperature regime where the conformation resembles a random walker and a low temperature limit where no thermal motions prevail. We find that the final low temperature conformations display a high level of degeneracy, in a manner which is fully in line with the presumed intrinsically disordered character of hIAPP. In particular, we identify an isolated family of alpha-helical conformations that might cause the transition to amyloidosis, by nucleation.
△ Less
Submitted 26 December, 2014;
originally announced December 2014.
-
Gauge field, strings, solitons, anomalies and the speed of life
Authors:
Antti J. Niemi
Abstract:
It's been said that "mathematics is biology's next microscope, only better; biology is mathematics' next physics, only better". Here we aim for something even better. We try to combine mathematical physics and biology into a picoscope of life. For this we merge techniques which have been introduced and developed in modern mathematical physics, largely by Ludvig Faddeev to describe objects such as…
▽ More
It's been said that "mathematics is biology's next microscope, only better; biology is mathematics' next physics, only better". Here we aim for something even better. We try to combine mathematical physics and biology into a picoscope of life. For this we merge techniques which have been introduced and developed in modern mathematical physics, largely by Ludvig Faddeev to describe objects such as solitons and Higgs and to explain phenomena such as anomalies in gauge fields. We propose a synthesis that can help to resolve the protein folding problem, one of the most important conundrums in all of science. We apply the concept of gauge invariance to scrutinize the extrinsic geometry of strings in three dimensional space. We evoke general principles of symmetry in combination with Wilsonian universality and derive an essentially unique Landau-Ginzburg energy that describes the dynamics of a generic string-like configuration in the far infrared. We observe that the energy supports topological solitons, that pertain to an anomaly in the manner how a string is framed around its inflection points. We explain how the solitons operate as modular building blocks from which folded proteins are composed. We describe crystallographic protein structures by multi-solitons with experimental precision, and investigate the non-equilibrium dynamics of proteins under varying temperature. We simulate the folding process of a protein at in vivo speed and with close to pico-scale accuracy using a standard laptop computer: With pico-biology as mathematical physics' next pursuit, things can only get better.
△ Less
Submitted 5 July, 2014; v1 submitted 29 June, 2014;
originally announced June 2014.
-
Long range correlations and folding angle in polymers with applications to α-helical proteins
Authors:
Andrey Krokhotin,
Stam Nicolis,
Antti J. Niemi
Abstract:
The conformational complexity of linear polymers far exceeds that of point-like atoms and molecules. Polymers can bend, twist, even become knotted. Thus they may also display a much richer phase structure than point particles. But it is not very easy to characterize the phase of a polymer. Essentially, the only attribute is the radius of gyration. The way how it changes when the degree of polymeri…
▽ More
The conformational complexity of linear polymers far exceeds that of point-like atoms and molecules. Polymers can bend, twist, even become knotted. Thus they may also display a much richer phase structure than point particles. But it is not very easy to characterize the phase of a polymer. Essentially, the only attribute is the radius of gyration. The way how it changes when the degree of polymerization becomes different, and how it evolves when the ambient temperature and solvent properties change, discloses the phase of the polymer. Moreover, in any finite length chain there are corrections to scaling, that complicate the detailed analysis of the phase structure. Here we introduce a quantity that we call the folding angle, a novel tool to identify and scrutinize the phases of polymers. We argue for a mean-field relationship between its values and those of the scaling exponent in the radius of gyration. But unlike in the case of the radius of gyration, the value of the folding angle can be evaluated from a single structure. As an example we estimate the value of the folding angle in the case of crystallographic α-helical protein structures in the Protein Data Bank (PDB). We also show how the value can be numerically computed using a theoretical model of α-helical chiral homopolymers.
△ Less
Submitted 22 June, 2013;
originally announced June 2013.
-
On the role of thermal backbone fluctuations in myoglobin ligand gate dynamics
Authors:
Andrey Krokhotin,
Antti J. Niemi,
Xubiao Peng
Abstract:
We construct an energy function that describes the crystallographic structure of spermwhale myoglobin backbone. As a model in our construction, we use the Protein Data Bank entry 1ABS that has been measured at liquid helium temperature. Consequently the thermal B-factor fluctuations are very small, which is an advantage in our construction. The energy function that we utilize resembles that of the…
▽ More
We construct an energy function that describes the crystallographic structure of spermwhale myoglobin backbone. As a model in our construction, we use the Protein Data Bank entry 1ABS that has been measured at liquid helium temperature. Consequently the thermal B-factor fluctuations are very small, which is an advantage in our construction. The energy function that we utilize resembles that of the discrete non-linear Schrodinger equation. Likewise, ours supports solitons as local minimum energy configurations. We describe the 1ABS backbone in terms of solitons with a precision that deviates from 1ABS by an average root-mean-square distance, which is less than the experimentally observed Debye-Waller B-factor fluctuation distance. We then subject the multisoliton solution to extensive numerical heating and cooling experiments, over a very wide range of temperatures. We concentrate in particular to temperatures above 300K and below the theta-point unfolding temperature, which is around 348K. We confirm that the behavior of the multisoliton is fully consistent with Anfinsen's principle, up to very high temperatures. We observe that the structure responds to an increase of temperature consistently in a very similar manner. This enables us to characterize the onset of thermally induced conformational changes in terms of three distinct backbone ligand gates. One of the gates is made of the helix F and the helix E. This is a pathway that is presumed to have a major role in ligand migration between the heme and the exterior. The two other gates are chosen similarly, when open they provide a direct access route for a ligand to reach the heme. We find that out of the three gates we investigate, the one which is formed by helices B and G is the most sensitive one to thermally induced conformational changes. Our approach provides a novel perspective to the important problem of ligand migration.
△ Less
Submitted 11 October, 2012;
originally announced October 2012.
-
Solitons and Collapse in the lambda-repressor protein
Authors:
Andrey Krokhotin,
Martin Lundgren,
Antti J. Niemi
Abstract:
The enterobacteria lambda phage is a paradigm temperate bacteriophage. Its lysogenic and lytic life cycles echo competition between the DNA binding $λ$-repressor (CI) and CRO proteins. Here we scrutinize the structure, stability and folding pathways of the $λ$-repressor protein, that controls the transition from the lysogenic to the lytic state. We first investigate the super-secondary helix-loop-…
▽ More
The enterobacteria lambda phage is a paradigm temperate bacteriophage. Its lysogenic and lytic life cycles echo competition between the DNA binding $λ$-repressor (CI) and CRO proteins. Here we scrutinize the structure, stability and folding pathways of the $λ$-repressor protein, that controls the transition from the lysogenic to the lytic state. We first investigate the super-secondary helix-loop-helix composition of its backbone. We use a discrete Frenet framing to resolve the backbone spectrum in terms of bond and torsion angles. Instead of four, there appears to be seven individual loops. We model the putative loops using an explicit soliton Ansatz. It is based on the standard soliton profile of the continuum nonlinear Schrödinger equation. The accuracy of the Ansatz far exceeds the B-factor fluctuation distance accuracy of the experimentally determined protein configuration. We then investigate the folding pathways and dynamics of the $λ$-repressor protein. We introduce a coarse-grained energy function to model the backbone in terms of the C$_α$ atoms and the side-chains in terms of the relative orientation of the C$_β$ atoms. We describe the folding dynamics in terms of relaxation dynamics, and find that the folded configuration can be reached from a very generic initial configuration. We conclude that folding is dominated by the temporal ordering of soliton formation. In particular, the third soliton should appear before the first and second. Otherwise, the DNA binding turn does not acquire its correct structure. We confirm the stability of the folded configuration by repeated heating and cooling simulations.
△ Less
Submitted 2 September, 2012;
originally announced September 2012.
-
Protein loops, solitons and side-chain visualization with applications to the left-handed helix region
Authors:
Martin Lundgren,
Antti J. Niemi,
Fan Sha
Abstract:
Folded proteins have a modular assembly. They are constructed from regular secondary structures like alpha-helices and beta-strands that are joined together by loops. Here we develop a visualization technique that is adapted to describe this modular structure. In complement to the widely employed Ramachandran plot that is based on toroidal geometry, our approach utilizes the geometry of a two-sphe…
▽ More
Folded proteins have a modular assembly. They are constructed from regular secondary structures like alpha-helices and beta-strands that are joined together by loops. Here we develop a visualization technique that is adapted to describe this modular structure. In complement to the widely employed Ramachandran plot that is based on toroidal geometry, our approach utilizes the geometry of a two-sphere. Unlike the more conventional approaches that only describea given peptide unit, ours is capable of describing the entire backbone environment including the neighboring peptide units. It maps the positions of each atom to the surface of the two-sphere exactly how these atoms are seen by an observer who is located at the position of the central C-alpha atom. At each level of side-chain atoms we observe a strong correlation between the positioning of the atom and the underlying local secondary structure with very little if any variation between the different amino acids. As a concrete example we analyze the left-handed helix region of non-glycyl amino acids. This region corresponds to an isolated and highly localized residue independent sector in the direction of the C-beta carbons on the two-sphere. We show that the residue independent localization extends to C-gamma and C-delta carbons, and to side-chain oxygen and nitrogen atoms in the case of asparagine and aspartic acid. When we extend the analysis to the side-chain atoms of the neighboring residues, we observe that left-handed beta-turns display a regular and largely amino acid independent structure that can extend to seven consecutive residues. This collective pattern is duu to the presence of a backbone soliton. We show how one can use our visualization techniques to analyze and classify the different solitons in terms of selection rules that we describe in detail.
△ Less
Submitted 2 September, 2012;
originally announced September 2012.
-
On correlation between protein secondary structure, backbone bond angles, and side-chain orientations
Authors:
Martin Lundgren,
Antti J. Niemi
Abstract:
We investigate the fine structure of the sp3 hybridized covalent bond geometry that governs the tetrahedral architecture around the central C$_α$ carbon of a protein backbone, and for this we develop new visualization techniques to analyze high resolution X-ray structures in Protein Data Bank. We observe that there is a correlation between the deformations of the ideal tetrahedral symmetry and the…
▽ More
We investigate the fine structure of the sp3 hybridized covalent bond geometry that governs the tetrahedral architecture around the central C$_α$ carbon of a protein backbone, and for this we develop new visualization techniques to analyze high resolution X-ray structures in Protein Data Bank. We observe that there is a correlation between the deformations of the ideal tetrahedral symmetry and the local secondary structure of the protein. We propose a universal coarse grained energy function to describe the ensuing side-chain geometry in terms of the C$_β$ carbon orientations. The energy function can model the side-chain geometry with a sub-atomic precision. As an example we construct the C$_α$-C$_β$ structure of HP35 chicken villin headpiece. We obtain a configuration that deviates less than 0.4 Ȧ in root-mean-square distance from the experimental X-ray structure.
△ Less
Submitted 2 September, 2012;
originally announced September 2012.
-
Soliton driven relaxation dynamics and universality in protein collapse
Authors:
Andrey Krokhotin,
Martin Lundgren,
Antti J. Niemi
Abstract:
Protein collapse can be viewed as a dynamical phase transition, during which new scales and collective variables become excited while the old ones recede and fade away. This causes formidable computational bottle-necks in approaches that are based on atomic scale scrutiny. Here we consider an effective dynamical Landau theory to model the folding process at biologically relevant time and distance…
▽ More
Protein collapse can be viewed as a dynamical phase transition, during which new scales and collective variables become excited while the old ones recede and fade away. This causes formidable computational bottle-necks in approaches that are based on atomic scale scrutiny. Here we consider an effective dynamical Landau theory to model the folding process at biologically relevant time and distance scales. We reach both a substantial decrease in the execution time and improvement in the accuracy of the final configuration, in comparison to more conventional approaches. As an example we inspect the collapse of HP35 chicken villin headpiece subdomain, where there are detailed molecular dynamics simulations to compare with. We start from a structureless, unbend and untwisted initial configuration. In less than one second of wall-clock time on a single processor personal computer we consistently reach the native state with 0.5 Angstrom root mean square distance (RMSD) precision. We confirm that our folding pathways are indeed akin those obtained in recent atomic level molecular dynamics simulations. We conclude that our approach appears to have the potential for a computationally economical method to accurately understand theoretical aspects of protein collapse.
△ Less
Submitted 8 November, 2011;
originally announced November 2011.
-
The protein dynamical transition is a pseudogap changeover
Authors:
Andrey Krokhotin,
Antti J. Niemi
Abstract:
The emergence of biochemical activities in a protein seem to commence with the onset of atomic mean-square displacements along the protein lattice. The ensuing protein dynamical transition has been discussed extensively in the literature, and often with conflicting conclusions. Here we clarify the phenomenon by establishing a deep connection between the dynamical transition and the pseudogap state…
▽ More
The emergence of biochemical activities in a protein seem to commence with the onset of atomic mean-square displacements along the protein lattice. The ensuing protein dynamical transition has been discussed extensively in the literature, and often with conflicting conclusions. Here we clarify the phenomenon by establishing a deep connection between the dynamical transition and the pseudogap state where high-temperature superconductivity comes to its end.
For this we first show how to endow proteins with an order parameter akin the quasiparticle wave function in superconductors. We then present universality arguments to claim that the protein dynamical transition takes place in tandem with a pseudogap transmutation. We confirm that available experimental data fully supports our proposal.
△ Less
Submitted 21 September, 2011;
originally announced September 2011.
-
Soliton concepts and the protein structure
Authors:
Andrei Krokhotin,
Antti J. Niemi,
Xubiao Peng
Abstract:
Structural classification shows that the number of different protein folds is surprisingly small. It also appears that proteins are built in a modular fashion, from a relatively small number of components. Here we propose to identify the modular building blocks of proteins with the dark soliton solution of a generalized discrete nonlinear Schrodinger equation. For this we show that practically all…
▽ More
Structural classification shows that the number of different protein folds is surprisingly small. It also appears that proteins are built in a modular fashion, from a relatively small number of components. Here we propose to identify the modular building blocks of proteins with the dark soliton solution of a generalized discrete nonlinear Schrodinger equation. For this we show that practically all protein loops can be obtained simply by scaling the size and by joining together a number of copies of the soliton, one after another. The soliton has only two loop specific parameters and we identify their possible values in Protein Data Bank. We show that with a collection of 200 sets of parameters, each determining a soliton profile that describes a different short loop, we cover over 90% of all proteins with experimental accuracy. We also present two examples that describe how the loop library can be employed both to model and to analyze the structure of folded proteins.
△ Less
Submitted 18 September, 2011;
originally announced September 2011.
-
Protein Regge Trajectories, Phase Coexistence and Physics of Alzheimer's Disease
Authors:
Andrei Krokhotin,
Antti J. Niemi
Abstract:
Alzheimer's disease causes severe neurodegeneration in the brain that leads to a certain death. The defining factor is the formation of extracellular senile amyloid plaques in the brain. However, therapeutic approaches to remove them have not been effective in humans, and so our understanding of the cause of Alzheimer's disease remains incomplete. Here we investigate physical processes that might…
▽ More
Alzheimer's disease causes severe neurodegeneration in the brain that leads to a certain death. The defining factor is the formation of extracellular senile amyloid plaques in the brain. However, therapeutic approaches to remove them have not been effective in humans, and so our understanding of the cause of Alzheimer's disease remains incomplete. Here we investigate physical processes that might relate to its onset. Instead of the extracellular amyloid, we scrutinize the intracellular domain of its precursor protein. We argue for a phenomenon that has never before been discussed in the context of polymer physics: Like ice and water together, the intracellular domain of the amyloid precursor protein forms a state of phase coexistence with another protein. This leads to an inherent instability that could well be among the missing pieces in the puzzle of Alzheimer's disease.
△ Less
Submitted 6 September, 2011;
originally announced September 2011.
-
Covalent bond symmetry breaking and protein secondary structure
Authors:
Martin Lundgren,
Antti J. Niemi
Abstract:
Both symmetry and organized breaking of symmetry have a pivotal rôle in our understanding of structure and pattern formation in physical systems, including the origin of mass in the Universe and the chiral structure of biological macromolecules. Here we report on a new symmetry breaking phenomenon that takes place in all biologically active proteins, thus this symmetry breaking relates to the ince…
▽ More
Both symmetry and organized breaking of symmetry have a pivotal rôle in our understanding of structure and pattern formation in physical systems, including the origin of mass in the Universe and the chiral structure of biological macromolecules. Here we report on a new symmetry breaking phenomenon that takes place in all biologically active proteins, thus this symmetry breaking relates to the inception of life. The unbroken symmetry determines the covalent bond geometry of a sp3 hybridized carbon atom. It dictates the tetrahedral architecture of atoms around the central carbon of an amino acid. Here we show that in a biologically active protein this symmetry becomes broken. Moreover, we show that the pattern of symmetry breaking is in a direct correspondence with the local secondary structure of the folded protein.
△ Less
Submitted 2 September, 2011;
originally announced September 2011.
-
On universal aspects of the left-handed helix region
Authors:
Martin Lundgren,
Antti J. Niemi,
Fan Sha
Abstract:
We inspect the geometry of proteins by identifying their backbones as framed polygons. We find that the left-handed helix region of the Ramachandran map for non-glycyl residues corresponds to an isolated and highly localized sector in the orientation of the $C_β$ carbons, when viewed in a Frenet frame that is centered at the corresponding $C_α$ carbons. We show that this localization in the orient…
▽ More
We inspect the geometry of proteins by identifying their backbones as framed polygons. We find that the left-handed helix region of the Ramachandran map for non-glycyl residues corresponds to an isolated and highly localized sector in the orientation of the $C_β$ carbons, when viewed in a Frenet frame that is centered at the corresponding $C_α$ carbons. We show that this localization in the orientation persists to $C_γ$ and $C_δ$ carbons. Furthermore, when we extend our analysis to the neighboring residues we conclude that the left-handed helix region reflects a very regular and apparently residue independent collective interplay of at least seven consecutive amino acids.
△ Less
Submitted 12 April, 2011;
originally announced April 2011.
-
Solitons and Physics of the Lysogenic to Lytic Transition in Enterobacteria Lambda Phage
Authors:
Andrei Krokhotine,
Antti J. Niemi
Abstract:
The lambda phage is a paradigm temperate bacteriophage. Its lysogenic and lytic life cycles echo competition between the DNA binding CI and CRO proteins. Here we address the Physics of this transition in terms of an energy function that portrays the backbone as a multi-soliton configuration. The precision of the individual solitons far exceeds the B-factor accuracy of the experimentally determined…
▽ More
The lambda phage is a paradigm temperate bacteriophage. Its lysogenic and lytic life cycles echo competition between the DNA binding CI and CRO proteins. Here we address the Physics of this transition in terms of an energy function that portrays the backbone as a multi-soliton configuration. The precision of the individual solitons far exceeds the B-factor accuracy of the experimentally determined protein conformations giving us confidence to conclude that three of the four loops are each composites of two closely located solitons. The only exception is the repressive DNA binding turn, it is the sole single soliton configuration of the backbone. When we compare the solitons with the Protein Data Bank we find that the one preceding the DNA recognition helix is unique to the CI protein, prompting us to conclude that the lysogenic to lytic transition is due to a saddle-node bifurcation involving a soliton-antisoliton annihilation that removes the first loop.
△ Less
Submitted 12 April, 2011;
originally announced April 2011.
-
The Discrete Frenet Frame, Inflection Point Solitons And Curve Visualization with Applications to Folded Proteins
Authors:
Shuangwei Hu,
Martin Lundgren,
Antti J. Niemi
Abstract:
We develop a transfer matrix formalism to visualize the framing of discrete piecewise linear curves in three dimensional space. Our approach is based on the concept of an intrinsically discrete curve, which enables us to more effectively describe curves that in the limit where the length of line segments vanishes approach fractal structures in lieu of continuous curves. We verify that in the case…
▽ More
We develop a transfer matrix formalism to visualize the framing of discrete piecewise linear curves in three dimensional space. Our approach is based on the concept of an intrinsically discrete curve, which enables us to more effectively describe curves that in the limit where the length of line segments vanishes approach fractal structures in lieu of continuous curves. We verify that in the case of differentiable curves the continuum limit of our discrete equation does reproduce the generalized Frenet equation. As an application we consider folded proteins, their Hausdorff dimension is known to be fractal. We explain how to employ the orientation of $C_β$ carbons of amino acids along a protein backbone to introduce a preferred framing along the backbone. By analyzing the experimentally resolved fold geometries in the Protein Data Bank we observe that this $C_β$ framing relates intimately to the discrete Frenet framing. We also explain how inflection points can be located in the loops, and clarify their distinctive rôle in determining the loop structure of foldel proteins.
△ Less
Submitted 28 February, 2011;
originally announced February 2011.
-
Towards Quantitative Classification of Folded Proteins in Terms of Elementary Functions
Authors:
Shuangwei Hu,
Andrei Krokhotin,
Antti J. Niemi,
Xubiao Peng
Abstract:
A comparative classification scheme provides a good basis for several approaches to understand proteins, including prediction of relations between their structure and biological function. But it remains a challenge to combine a classification scheme that describes a protein starting from its well organized secondary structures and often involves direct human involvement, with an atomary level Phys…
▽ More
A comparative classification scheme provides a good basis for several approaches to understand proteins, including prediction of relations between their structure and biological function. But it remains a challenge to combine a classification scheme that describes a protein starting from its well organized secondary structures and often involves direct human involvement, with an atomary level Physics based approach where a protein is fundamentally nothing more than an ensemble of mutually interacting carbon, hydrogen, oxygen and nitrogen atoms. In order to bridge these two complementary approaches to proteins, conceptually novel tools need to be introduced. Here we explain how the geometrical shape of entire folded proteins can be described analytically in terms of a single explicit elementary function that is familiar from nonlinear physical systems where it is known as the kink-soliton. Our approach enables the conversion of hierarchical structural information into a quantitative form that allows for a folded protein to be characterized in terms of a small number of global parameters that are in principle computable from atomary level considerations. As an example we describe in detail how the native fold of the myoglobin 1M6C emerges from a combination of kink-solitons with a very high atomary level accuracy. We also verify that our approach describes longer loops and loops connecting $α$-helices with $β$-strands, with same overall accuracy.
△ Less
Submitted 2 December, 2010; v1 submitted 13 November, 2010;
originally announced November 2010.
-
Knots and Swelling in Protein Folding
Authors:
Martin Lundgren,
Antti J. Niemi
Abstract:
Proteins can sometimes be knotted, and for many reasons the study of knotted proteins is rapidly becoming very important. For example, it has been proposed that a knot increases the stability of a protein. Knots may also alter enzymatic activities and enhance binding. Moreover, knotted proteins may even have some substantial biomedical significance in relation to illnesses such as Parkinson's di…
▽ More
Proteins can sometimes be knotted, and for many reasons the study of knotted proteins is rapidly becoming very important. For example, it has been proposed that a knot increases the stability of a protein. Knots may also alter enzymatic activities and enhance binding. Moreover, knotted proteins may even have some substantial biomedical significance in relation to illnesses such as Parkinson's disease. But to a large extent the biological role of knots remains a conundrum. In particular, there is no explanation why knotted proteins are so scarce. Here we argue that knots are relatively rare because they tend to cause swelling in proteins that are too short, and presently short proteins are over-represented in the Protein Data Bank (PDB). Using Monte Carlo simulations we predict that the figure-8 knot leads to the most compact protein configuration when the number of amino acids is in the range of 200-600. For the existence of the simplest knot, the trefoil, we estimate a theoretical upper bound of 300-400 amino acids, in line with the available PDB data.
△ Less
Submitted 26 June, 2009;
originally announced June 2009.
-
A Gauge Field Theory of Chirally Folded Homopolymers with Applications to Folded Proteins
Authors:
Ulf H. Danielsson,
Martin Lundgren,
Antti J. Niemi
Abstract:
We combine the principle of gauge invariance with extrinsic string geometry to develop a lattice model that can be employed to theoretically describe properties of chiral, unbranched homopolymers. We find that in its low temperature phase the model is in the same universality class with proteins that are deposited in the Protein Data Bank, in the sense of the compactness index. We apply the model…
▽ More
We combine the principle of gauge invariance with extrinsic string geometry to develop a lattice model that can be employed to theoretically describe properties of chiral, unbranched homopolymers. We find that in its low temperature phase the model is in the same universality class with proteins that are deposited in the Protein Data Bank, in the sense of the compactness index. We apply the model to analyze various statistical aspects of folded proteins. Curiously we find that it can produce results that are a very good good match to the data in the Protein Data Bank.
△ Less
Submitted 26 August, 2010; v1 submitted 17 February, 2009;
originally announced February 2009.