Search

Topological methods are indispensable in theoretical studies of particle physics, condensed matter physics, and gravity. These powerful techniques have also been applied to biological physics. For example, knowledge of DNA topology is pivotal to the understanding as to how living cells function. Here, the biophysical repertoire of topological methods is extended, with the aim to understand and characterize the global structure of a folded protein. For this, the elementary concept of winding number of a vector field on a plane is utilized to introduce a topological quantity called the folding index of a crystallographic protein. It is observed that in the case of high resolution protein crystals, the folding index, when evaluated over the entire length of the crystallized protein backbone, has a very clear and strong propensity towards integer values. The observation proposes that the way how a protein folds into its biologically active conformation is a structural self-organization process with a topological facet that relates to the concept of solitons. It is proposed that the folding index has a potential to become a useful tool for the global, topological characterization of the folding pathways.

In all-atom molecular simulation studies of proteins, each atom in the protein is represented by a point mass and interactions are defined in terms of the atomic positions. In recent years, various simplified approaches have been proposed. These approaches aim to improve computational efficiency and to provide a better physical insight. The simplified models can differ widely in their description of the geometry and the interactions inside the protein. This study explores the most fundamental choice in the simplified protein models: the choice of a coordinate set defining the protein structure. A simplified model can use fewer point masses than the all-atom model and/or eliminate some of the internal coordinates of the molecule by setting them to an average or ideal value. We look at the implications of such choices for the overall protein structure. We find that care must be taken for angular coordinates, where even very small variations can lead to significant changes in the positions of far away atoms. In particular, we show that the phi/psi torsion angles are not a sufficient coordinate set, whereas another coordinate set with two degrees of freedom per residue, virtual C-alpha backbone bond, and torsion angles performs satisfactorily.

The theory of string-like continuous curves and discrete chains have numerous important physical applications. Here we develop a general geometrical approach, to systematically derive Hamiltonian energy functions for these objects. In the case of continuous curves, we demand that the energy function must be invariant under local frame rotations, and it should also transform covariantly under reparametrizations of the curve. This leads us to consider energy functions that are constructed from the conserved quantities in the hierarchy of the integrable nonlinear Schrödinger equation. We point out the existence of a Weyl transformation that we utilize to introduce a dual hierarchy to the standard nonlinear Schrödinger equation hierarchy. We propose that the dual hierarchy is also integrable, and we confirm this to the first nontrivial order. In the discrete case the requirement of reparametrization invariance is void. But the demand of invariance under local frame rotations prevails, and we utilize it to introduce a discrete variant of the Zakharov-Shabat recursion relation. We use this relation to derive frame-independent quantities that we propose are the essentially unique and as such natural candidates for constructing energy functions for piecewise linear polygonal chains. We also investigate the discrete version of the Weyl duality transformation. We confirm that in the continuum limit the discrete energy functions go over to their continuum counterparts, including the perfect derivative contributions.

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 Schrodinger 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 Theta-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.

Protein collapse from a random chain to the native state involves a dynamical phase transition. During the process, new scales and collective variables become excited while old ones recede and fade away. The presence of different phases and many scales causes formidable computational bottle-necks in approaches that are based on full atomic scale scrutiny. Here we propose a way to describe the folding and unfolding processes effectively, using only the biologically relevant time and distance scales. We merge a coarse grained Landau theory that models the static collapsed protein in the low-temperature limit with a Glauber protocol that describes finite-temperature relaxation dynamics in a statistical system which is out of thermal equilibrium. As an example we inspect the collapse of a HP35 chicken villin headpiece subdomain, a paradigm specimen in protein folding studies. We simulate the folding and unfolding process by repeated heating and cooling cycles between a given low-temperature, i.e. bad solvent, environment where the protein is collapsed and various different high-temperature, i.e. good solvent, environments. We find that as long as the high temperature value stays below a value in the range that separates the random walk phase from the self-avoiding walk phase, we consistently recover the native state upon cooling. But, when heated to sufficiently high temperatures, the native state practically never recurs. Our result confirms Anfinsen's thermodynamical hypothesis and estimates a temperature range for its validity, in the case of villin.

The enterobacteria lambda phage is a paradigm temperate bacteriophage. Its lysogenic and lytic life cycles echo competition between the DNA binding lambda-repressor (CI) and CRO proteins. Here we scrutinize the structure, stability, and folding pathways of the lambda-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 Schrodinger 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 lambda-repressor protein. We introduce a coarse-grained energy function to model the backbone in terms of the C-alpha atoms and the side chains in terms of the relative orientation of the C-beta 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.

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.

We investigate the fine structure of the sp3 hybridized covalent bond geometry that governs the tetrahedral architecture around the central C-alpha 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-beta carbon orientations. The energy function can model the side-chain geometry with a subatomic precision. As an example we construct the C-alpha-C-beta structure of HP35 chicken villin headpiece. We obtain a configuration that deviates less than 0.4 angstrom in root-mean-square distance from the experimental x-ray structure.

We introduce a novel generalization of the discrete nonlinear Schrodinger equation. It supports solitons that we utilize to model chiral polymers in the collapsed phase and, in particular, proteins in their native state. As an example we consider the villin headpiece HP35, an archetypal protein for testing both experimental and theoretical approaches to protein folding. We use its backbone as a template to explicitly construct a two-soliton configuration. Each of the two solitons describe well over 7.000 supersecondary structures of folded proteins in the Protein Data Bank with sub-angstrom accuracy suggesting that these solitons are common in nature.

Recently it has been pointed out that the so-called Faddeev-Niemi equations that describe the Yang-Mills equations of motion in terms of a decomposed gauge field, can have solutions that obey the standard Yang-Mills equations with a source term. Here we argue that the source term is covariantly constant. Furthermore, we find that there are solutions of the Yang-Mills equation with a covariantly constant source term that are not solutions to the Faddeev-Niemi equations. We also present a general class of gauge field configurations that obey the Faddeev-Niemi equation but do not solve the Yang-Mills equation. We propose that these configurations might have physical relevance in a strongly coupled phase, where spin-charge separation takes place and the Yang-Mills theory cannot be described in terms of a Landau liquid of asymptotically free gluons.

We present a numerical Monte Carlo analysis of a continuos spin Ising chain that can describe the statistical proterties of folded proteins. We find that depending on the value of the Metropolis temperature, the model displays the three known nontrivial phases of polymers: At low temperatures the model is in a collapsed phase, at medium temperatures it is in a random walk phase, and at high temperatures it enters the self-avoiding random walk phase. By investigating the temperature dependence of the specific energy we confirm that the transition between the collapsed phase and the random walk phase is a phase transition, while the random walk phase and self-avoiding random walk phase are separated from each other by a cross-over transition. We also compare the predictions of the model to a phenomenological elastic energy formula, proposed by Huang and Lei to describe 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-beta 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-beta 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 role in determining the loop structure of folded proteins.

We present a numerical Monte Carlo analysis of the phase structure in a continuous spin Ising chain that describes chiral homopolymers. We find that depending on the value of the Metropolis temperature, the model displays the three known nontrivial phases of polymers: At low temperatures the model is in a collapsed phase, at medium temperatures it is in a random walk phase, and at high temperatures it enters the self-avoiding random walk phase. By investigating the temperature dependence of the specific energy we confirm that the transition between the collapsed phase and the random walk phase is a phase transition, while the random walk phase and self-avoiding random walk phase are separated from each other by a crossover transition. We propose that the model can be applied to characterize the statistical properties of protein folding. For this we compare the predictions of the model to a phenomenological elastic energy formula, proposed by J. Lei and K. Huang [e-print arXiv:1002.5013; Europhys. Lett. 88, 68004 (2009)] to describe folded proteins.

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.

Previously, it has been shown that in a spin-charge separated SU(2) Yang– Mills theory, (Euclidean) spacetime rotation invariance can be broken by an infinitesimal 1-cocycle that appears in the S O(4) boosts. Here we study in detail the structure of this 1-cocycle. In particular, we show that its non- triviality relates to the presence of a (Dirac) magnetic monopole bundle. We also compute the finite version of the cocycle.

We introduce a new supercurrent in the electroweak sector of the standard model. Its interaction with the hypergauge field influences the mass of the Z boson but has no effect on the W-+/- boson masses. In the leptonic sector it affects the numerical value of the vector and axial coupling constants between neutral currents and the Z boson, and a comparison with experimental values yields an upper bound to the strength of the coupling between the supercurrent and the hypergauge field. In the baryonic sector the supercurrent gives a new contribution to the anomaly equation for baryon number current. As a consequence it may have an effect on baryogenesis.

A closed framed curve can be characterized by the value of its self-linking number. However, in general the self-linking number is not uniquely determined. In particular, it depends on the way how the curve has been framed. Moreover, the self-linking number can also change, when a bifurcation called perestroika takes place. Here we devise a simple Hamiltonian energy function to study perestroikas during the time evolution of a closed, piecewise linear polygonal chain. We analyze several examples to follow the progress of the discrete Frenet framing during the time evolution of the chain, and we observe how self-linking number changing perestroikas occur whenever there is an inflection point along the chain.