Towards quantitative classification of folded proteins in terms of elementary functions

Exploring Structural Flexibility and Stability of α-Synuclein by the Landau-Ginzburg-Wilson Approach

α-Synuclein (αS) is the principal protein component of the Lewy body and Lewy neurite deposits that are found in the brains of the victims of one of the most prevalent neurodegenerative disorders, Parkinson's disease. αS can be qualified as a chameleon protein because of the large number of different conformations that it is able to adopt: it is disordered under physiological conditions in solution, in equilibrium with a minor α-helical tetrameric form in the cytoplasm, and is α-helical when bound to a cell membrane. Also, in vitro, αS forms polymorphic amyloid fibrils with unique arrangements of cross-β-sheet motifs. Therefore, it is of interest to elucidate the origins of the structural flexibility of αS and what makes αS stable in different conformations. We address these questions here by analyzing the experimental structures of the micelle-bound, tetrameric, and fibrillar αS in terms of a kink (heteroclinic standing wave solution) of a generalized discrete nonlinear Schrödinger equation. It is illustrated that without molecular dynamics simulations the kinks are capable of identifying the key residues causing structural flexibility of αS. Also, the stability of the experimental structures of αS is investigated by simulating heating/cooling trajectories using the Glauber algorithm. The findings are consistent with experiments.

ESCASA: Analytical estimation of atomic coordinates from coarse-grained geometry for nuclear-magnetic-resonance-assisted protein structure modeling. I. Backbone and Hβ protons

A method for the estimation of coordinates of atoms in proteins from coarse-grained geometry by simple analytical formulas (ESCASA), for use in nuclear-magnetic-resonance (NMR) data-assisted coarse-grained simulations of proteins is proposed. In this paper, the formulas for the backbone H^α and amide (H^N ) protons, and the side-chain H^β protons, given the C^α -trace, have been derived and parameterized, by using the interproton distances calculated from a set of 140 high-resolution non-homologous protein structures. The mean standard deviation over all types of proton pairs in the set was 0.44 Å after fitting. Validation against a set of 41 proteins with NMR-determined structures, which were not considered in parameterization, resulted in average standard deviation from average proton-proton distances of the NMR-determined structures of 0.25 Å, compared to 0.21 Å obtained with the PULCHRA all-atom-chain reconstruction algorithm and to the 0.12 Å standard deviation of the average-structure proton-proton distance of NMR-determined ensembles. The formulas provide analytical forces and can, therefore, be used in coarse-grained molecular dynamics.

Investigation of Phosphorylation-Induced Folding of an Intrinsically Disordered Protein by Coarse-Grained Molecular Dynamics

Apart from being the most common mechanism of regulating protein function and transmitting signals throughout the cell, phosphorylation has an ability to induce disorder-to-order transition in an intrinsically disordered protein. In particular, it was shown that folding of the intrinsically disordered protein, eIF4E-binding protein isoform 2 (4E-BP2), can be induced by multisite phosphorylation. Here, the principles that govern the folding of phosphorylated 4E-BP2 (pT37pT46 4E-BP2_18-62) are investigated by analyzing canonical and replica exchange molecular dynamics trajectories, generated with the coarse-grained united-residue force field, in terms of local and global motions and the time dependence of formation of contacts between C^αs of selected pairs of residues. The key residues involved in the folding of the pT37pT46 4E-BP2_18-62 are elucidated by this analysis. The correlations between local and global motions are identified. Moreover, for a better understanding of the physics of the formation of the folded state, the experimental structure of the pT37pT46 4E-BP2_18-62 is analyzed in terms of a kink (heteroclinic standing wave solution) of a generalized discrete nonlinear Schrödinger equation. It is shown that without molecular dynamics simulations the kinks are able to identify not only the phosphorylated sites of protein, the key players in folding, but also the reasons for the weak stability of the pT37pT46 4E-BP2_18-62.

New Insights into Folding, Misfolding, and Nonfolding Dynamics of a WW Domain

Intermediate states in protein folding are associated with formation of amyloid fibrils, which are responsible for a number of neurodegenerative diseases. Therefore, prevention of the aggregation of folding intermediates is one of the most important problems to overcome. Recently, we studied the origins and prevention of formation of intermediate states with the example of the Formin binding protein 28 (FBP28) WW domain. We demonstrated that the replacement of Leu26 by Asp26 or Trp26 (in ∼15% of the folding trajectories) can alter the folding scenario from three-state folding, a major folding scenario for the FBP28 WW domain (WT) and its mutants, toward two-state or downhill folding at temperatures below the melting point. Here, for a better understanding of the physics of the formation/elimination of intermediates, (i) the dynamics and energetics of formation of β-strands in folding, misfolding, and nonfolding trajectories of these mutants (L26D and L26W) is investigated; (ii) the experimental structures of WT, L26D, and L26W are analyzed in terms of a kink (heteroclinic standing wave solution) of a generalized discrete nonlinear Schrödinger equation. We show that the formation of each β-strand in folding trajectories is accompanied by the emergence of kinks in internal coordinate space as well as a decrease in local free energy. In particular, the decrease in downhill folding trajectory is ∼7 kcal/mol, while it varies between 31 and 48 kcal/mol for the three-state folding trajectory. The kink analyses of the experimental structures give new insights into formation of intermediates, which may become a useful tool for preventing aggregation.

Topological Indices of Proteins

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 folded states of proteins, for example, using the self-linking number of the corresponding framed curves. In this paper we extend this classification to the discrete version, taking advantage of the “randomness” of such curves. 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. However, 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. This works well for the proteins in our analysis, for which 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.

Solution x-ray scattering and structure formation in protein dynamics

We propose a computationally effective approach that builds on Landau mean-field theory in combination with modern nonequilibrium statistical mechanics 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 analyzing experimental data in the case of Engrailed homeodomain protein as an example. We demonstrate how to interpret S/WAXS data qualitatively with a good precision and over an extended temperature range. We explain experimental observations in terms of protein phase structure, and we make predictions for future experiments and for how to analyze data at different ambient temperature values. We conclude that the approach we propose has the potential to become a highly accurate, computationally effective, and predictive tool for analyzing S/WAXS data. For this, we compare our results with those obtained previously in an all-atom molecular dynamics simulation.

Bloch spin waves and emergent structure in protein folding with HIV envelope glycoprotein as an example

We inquire how structure emerges during the process of protein folding. For this we scrutinize collective 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 Schrödinger equation. The example we consider is an α-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 supersecondary 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 precision of tenths to several angstroms in terms of a soliton solution to a discrete nonlinear Schrödinger equation.

Clustering and percolation in protein loop structures

Background

High precision protein loop modelling remains a challenge, both in template based and template independent approaches to protein structure prediction.

Method

We introduce the concepts of protein loop clustering and percolation, to develop a quantitative approach to systematically classify the modular building blocks of loops in crystallographic folded proteins. These fragments are all different parameterisations of a unique kink solution to a generalised discrete nonlinear Schrödinger (DNLS) equation. Accordingly, the fragments are also local energy minima of the ensuing energy function.

Results

We show how the loop fragments cover practically all ultrahigh resolution crystallographic protein structures in Protein Data Bank (PDB), with a 0.2 Ångström root-mean-square (RMS) precision. We find that no more than 12 different loop fragments are needed, to describe around 38 % of ultrahigh resolution loops in PDB. But there is also a large number of loop fragments that are either unique, or very rare, and examples of unique fragments are found even in the structure of a myoglobin.

Conclusions

Protein loops are built in a modular fashion. The loops are composed of fragments that can be modelled by the kink of the DNLS equation. The majority of loop fragments are also common, which are shared by many proteins. These common fragments are probably important for supporting the overall protein conformation. But there are also several fragments that are either unique to a given protein, or very rare. Such fragments are probably related to the function of the protein. Furthermore, we have found that the amino acid sequence does not determine the structure in a unique fashion. There are many examples of loop fragments with an identical amino acid sequence, but with a very different structure.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-015-0049-x) contains supplementary material, which is available to authorized users.

Aspects of structural landscape of human islet amyloid polypeptide

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 nuclear magnetic resonance 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 α-helical conformations that might cause the transition to amyloidosis, by nucleation.

Kinks, loops, and protein folding, with protein A as an example

The dynamics and energetics of formation of loops in the 46-residue N-terminal fragment of the B-domain of staphylococcal protein A has been studied. Numerical simulations have been performed using coarse-grained molecular dynamics with the united-residue (UNRES) force field. The results have been analyzed in terms of a kink (heteroclinic standing wave solution) of a generalized discrete nonlinear Schrödinger (DNLS) equation. In the case of proteins, the DNLS equation arises from a C(α)-trace-based energy function. Three individual kink profiles were identified in the experimental three-α-helix structure of protein A, in the range of the Glu16-Asn29, Leu20-Asn29, and Gln33-Asn44 residues, respectively; these correspond to two loops in the native structure. UNRES simulations were started from the full right-handed α-helix to obtain a clear picture of kink formation, which would otherwise be blurred by helix formation. All three kinks emerged during coarse-grained simulations. It was found that the formation of each is accompanied by a local free energy increase; this is expressed as the change of UNRES energy which has the physical sense of the potential of mean force of a polypeptide chain. The increase is about 7 kcal/mol. This value can thus be considered as the free energy barrier to kink formation in full α-helical segments of polypeptide chains. During the simulations, the kinks emerge, disappear, propagate, and annihilate each other many times. It was found that the formation of a kink is initiated by an abrupt change in the orientation of a pair of consecutive side chains in the loop region. This resembles the formation of a Bloch wall along a spin chain, where the C(α) backbone corresponds to the chain, and the amino acid side chains are interpreted as the spin variables. This observation suggests that nearest-neighbor side chain-side chain interactions are responsible for initiation of loop formation. It was also found that the individual kinks are reflected as clear peaks in the principal modes of the analyzed trajectory of protein A, the shapes of which resemble the directional derivatives of the kinks along the chain. These observations suggest that the kinks of the DNLS equation determine the functionally important motions of proteins.

Gauge fields, strings, solitons, anomalies, and the speed of life

Joel Cohen proposed 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 that were 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 stringlike configuration in the far infrared. We observe that the energy supports topological solitons that relate to an anomaly similarly to 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 multisolitons with experimental precision and investigate the nonequilibrium dynamics of proteins under temperature variations. We simulate the folding process of a protein at in vivo speed and with close to picoscale accuracy using a standard laptop computer. With picobiology as next pursuit of mathematical physics, things can only get better.

On the role of thermal backbone fluctuations in myoglobin ligand gate dynamics

We construct an energy function that describes the crystallographic structure of sperm whale 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 nonlinear Schrödinger equation. Likewise, ours supports topological solitons as local minimum energy configurations. We describe the 1ABS backbone in terms of topological 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 topological multi-soliton solution to extensive numerical heating and cooling experiments, over a very wide range of temperatures. We concentrate in particular to temperatures above 300 K and below the Θ-point unfolding temperature, which is around 348 K. We confirm that the behavior of the topological multi-soliton is fully consistent with Anfinsen's thermodynamic 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. 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 to thermally induced conformational changes. Our approach provides a novel perspective to the important problem of ligand entry and exit.

Coexistence of phases in a protein heterodimer

A heterodimer consisting of two or more different kinds of proteins can display an enormous number of distinct molecular architectures. The conformational entropy is an essential ingredient in the Helmholtz free energy and, consequently, these heterodimers can have a very complex phase structure. Here, it is proposed that there is a state of proteins, in which the different components of a heterodimer exist in different phases. For this purpose, the structures in the protein data bank (PDB) have been analyzed, with radius of gyration as the order parameter. Two major classes of heterodimers with their protein components coexisting in different phases have been identified. An example is the PDB structure 3DXC. This is a transcriptionally active dimer. One of the components is an isoform of the intra-cellular domain of the Alzheimer-disease related amyloid precursor protein (AICD), and the other is a nuclear multidomain adaptor protein in the Fe65 family. It is concluded from the radius of gyration that neither of the two components in this dimer is in its own collapsed phase, corresponding to a biologically active protein. The UNRES energy function has been utilized to confirm that, if the two components are separated from each other, each of them collapses. The results presented in this work show that heterodimers whose protein components coexist in different phases, can have intriguing physical properties with potentially important biological consequences.

Correlation between protein secondary structure, backbone bond angles, and side-chain orientations

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 the 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 subatomic 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.

Solitons and collapse in the λ-repressor protein

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, which controls the transition from the lysogenic to the lytic state. We first investigate the supersecondary 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.

Soliton concepts and protein structure

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 that the modular building blocks are made of the dark soliton solution of a generalized discrete nonlinear Schrödinger equation. We find 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 compute their statistical distribution in the Protein Data Bank (PDB). We explicitly construct a collection of 200 sets of parameters, each determining a soliton profile that describes a different short loop. The ensuing profiles cover practically all those proteins in PDB that have a resolution which is better than 2.0 Å, with a precision such that the average root-mean-square distance between the loop and its soliton is less than the experimental B-factor fluctuation distance. We also present two examples that describe how the loop library can be employed both to model and to analyze folded proteins.

Discrete Frenet frame, inflection point solitons, and curve visualization with applications to folded proteins

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. This 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 reproduces the generalized Frenet equation. In particular, we draw attention to the conceptual similarity between inflection points where the curvature vanishes and topologically stable solitons. 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 (a.k.a. soliton centers) can be located in the loops and clarify their distinctive rôle in determining the loop structure of folded proteins.