The aqueous solvent profoundly influences protein folding, yet its effects are relatively poorly understood. In this study, we investigate the impact of solvation on the folding of lattice proteins by using Monte Carlo simulations. The proteins are modelled as self-avoiding 27-mer chains on a cubic lattice, with compact native states and structure-based Gō potentials. Each residue that makes no contacts with other residues in a given protein conformation is assigned a solvation energy ε_s , representing its full exposure to the solvent. We find that a negative ε_s , indicating a favorable solvation, increases the cooperativity of the folding transition by lowering the free energy of the unfolded state, increasing the folding free energy barrier, and narrowing the folding routes. This favorable solvation also significantly improves the correlation between folding rates and the native topology, measured by the relative contact order. Our results suggest that Gō model may overestimate the importance of native interactions and a solvation potential countering the native bias can play a significant role. The solvation energy in our model can be related to the polar interaction between w
Folding can transform mundane objects such as napkins into stunning works of art. However, finding new folding transformations for sheet materials is a challenging problem that requires expertise and real-world experimentation. In this paper, we present Modal Folding -- an automated approach for discovering energetically optimal folding transformations, i.e., large deformations that require little mechanical work. For small deformations, minimizing internal energy for fixed displacement magnitudes leads to the well-known elastic eigenmodes. While linear modes provide promising directions for bending, they cannot capture the rotational motion required for folding. To overcome this limitation, we introduce strain-space modes -- nonlinear analogues of elastic eigenmodes that operate on per-element curvatures instead of vertices. Using strain-space modes to determine target curvatures for bending elements, we can generate complex nonlinear folding motions by simply minimizing the sheet's internal energy. Our modal folding approach offers a systematic and automated way to create complex designs. We demonstrate the effectiveness of our method with simulation results for a range of shapes
Origami is the art of paper folding, and it borrows its name from two Japanese words \emph{ori} and \emph{kami}. In Japanese, {ori} means folding, and the paper is called {kami}. While origami is just a hobby to most, there is a lot more to it. If you fold a square sheet of paper into any of the traditional origami model (for example the flapping bird) and unfold it, you can see crease patterns. These crease patterns tell us that there is a lot of geometry hidden behind the folds. In this article, we investigate the symbiotic relationship between mathematics and origami. The first part of this article explores the utility of origami in education. We will see how origami could become an effective way of teaching methods of geometry, mainly because of its experiential nature. Complex origami patterns cannot be created out of thin air. They usually involve understanding deep mathematical theories and the ability to apply them to paper folding. In the second part of the article, we attempt to provide a glimpse of this beautiful connection between origami and mathematics.
Folding is emerging as a promising manufacturing process to transform flat materials into functional structures, offering efficiency by reducing the need for welding, gluing, and molding, while minimizing waste and enabling automation. Designing target shapes requires not only to determine cuts and folds, but also folding pathways. Simple combinatorics is impractical as the possibilities grow factorially with the number of folds. To address this, we present a graph-based algorithm for polyhedral shapes. By representing the target shape as a graph, where nodes correspond to faces and edges represent adjacency, the algorithm identifies all possible fold sequences and maps the configuration space into a discrete set of intermediate configurations. This systematic mapping is critical for the design of optimized processes, the simplifying of folding operations, the reduction of failures, and the improvement of manufacturing reliability.
The rates of protein folding with photon absorption or emission and the cross section of photon -protein inelastic scattering are calculated from the quantum folding theory by use of standard field-theoretical method. All these protein photo-folding processes are compared with common protein folding without interaction of photons (nonradiative folding). It is demonstrated that there exists a common factor (thermo-averaged overlap integral of vibration wave function, TAOI) for protein folding and protein photo-folding. Based on this finding it is predicted that: 1) the stimulated photo-folding rates show the same temperature dependence as protein folding; 2) the spectral line of electronic transition is broadened to a band which includes abundant vibration spectrum without and with conformational transition and the width of the vibration spectral line is largely reduced; 3) the resonance fluorescence cross section changes with temperature obeying the same law (Luo-Lu's law). The particular form of the folding rate - temperature relation and the abundant spectral structure imply the existence of a set of quantum oscillators in the transition process and these oscillators are mainly o
It is well known that the set of origami constructible numbers is larger than the classical straight-edge and compass constructible numbers. However, the Huzita-Justin-Hatori origami constructible numbers remain algebraic so that the transcendental number $π$ can only be approximated using a finite number of straight line folds. Using these methods we give a convergent sequence for folding $π$ as well as other methods to approximate $π$. Folding along curved creases, however, allows for the construction of transcendental numbers. We here give a method to construct $π$ exactly by folding along a parabola, and we discuss generalizations for folding other transcendental numbers such as $Γ(1/4)$.
We adopt the point of view that analysis of the stability of the protein folding process is central to understanding the underlying physics of folding. Stability of the folding process means that many perturbations do not disrupt the progress from the random coil to the native state. In this paper we explore the stability of folding using established methods from physics and mathematics. Our result is a preliminary theory of the physics of folding. We suggest some tests of these ideas using folding simulations. We begin by supposing that folding events are related in some way to mechanical waves on the molecule. We adopt an analytical approach to the physics which was pioneered by M.V. Berry, (in another context), based upon mathematics developed mainly by R. Thom and V.I. Arnold. We find that the stability of the folding process can be understood in terms of structures known as caustics, which occur in many kinds of wave phenomena. The picture that emerges is that natural selection has given us a set of protein molecules which have mechanical waves that propagate according to several mathematically specific restrictions. Successful simulations of folding can be used to test and co
We develop an operator algebraic framework for infinite games with a continuum of agents and prove that regret based learning dynamics governed by a noncommutative continuity equation converge to a unique quantal response equilibrium under mild regularity assumptions. The framework unifies functional analysis, coarse geometry and game theory by assigning to every game a von Neumann algebra that represents collective strategy evolution. A reflective regret operator within this algebra drives the flow of strategy distributions and its fixed point characterises equilibrium. We introduce the ordinal folding index, a computable ordinal valued metric that measures the self referential depth of the dynamics, and show that it bounds the transfinite time needed for convergence, collapsing to zero on coarsely amenable networks. The theory yields new invariant subalgebra rigidity results, establishes existence and uniqueness of envy free and maximin share allocations in continuum economies, and links analytic properties of regret flows with empirical stability phenomena in large language models. These contributions supply a rigorous mathematical foundation for large scale multi agent systems
The Ordinal Folding Index (OFI) is a new, fully computable yard-stick that measures how many rounds of self-reference a statement, protocol or position must unfold before its truth or outcome stabilises. By turning this abstract 'fold-back' depth into a single ordinal number, OFI forges a direct link between areas that are usually studied in isolation: the closure stages of fixed-point logics, the time-to-win values of infinite parity games, and the ordinal progressions that calibrate the strength of formal theories. We prove that OFI refines all classical game-theoretic and logical metrics while remaining algorithmically enumerable, supply a polynomial-time approximation scheme on finite arenas, and show how the index coincides exactly with the length of the shortest winning strategy in the associated evaluation game. Alongside the theory we outline five open problems from the completeness of the computable-ordinal spectrum to the possibility of 'compressing' deep self-reference that chart a research programme at the intersection of computer-aided logic, algorithmic game theory and ordinal analysis. OFI thus invites game theorists and logicians alike to view infinite play, transfi
We prove several hardness results on folding origami crease patterns. Flat-folding finite crease patterns is fixed-parameter tractable in the ply of the folded pattern (how many layers overlap at any point) and the treewidth of an associated cell adjacency graph. Under the exponential time hypothesis, the singly-exponential dependence of our algorithm on treewidth is necessary, even for bounded ply. Improving the dependence on ply would require progress on the unsolved map folding problem. Finding the shape of a polyhedron folded from a net with triangular faces and integer edge lengths is not possible in algebraic computation tree models of computation that at each tree node allow either the computation of arbitrary integer roots of real numbers, or the extraction of roots of polynomials with bounded degree and integer coefficients. For a model of reconfigurable origami with origami squares are attached at one edge by a hinge to a rigid surface, moving from one flat-folded state to another by changing the position of one square at a time is PSPACE-complete, and counting flat-folded states is #P-complete. For self-similar square crease patterns with infinitely many folds, testing f
Repeat proteins are made with tandem copies of similar amino acid stretches that fold into elongated architectures. Due to their symmetry, these proteins constitute excellent model systems to investigate how evolution relates to structure, folding and function. Here, we propose a scheme to map evolutionary information at the sequence level to a coarse-grained model for repeat-protein folding and use it to investigate the folding of thousands of repeat-proteins. We model the energetics by a combination of an inverse Potts model scheme with an explicit mechanistic model of duplications and deletions of repeats to calculate the evolutionary parameters of the system at single residue level. This is used to inform an Ising-like model that allows for the generation of folding curves, apparent domain emergence and occupation of intermediate states that are highly compatible with experimental data in specific case studies. We analyzed the folding of thousands of natural Ankyrin-repeat proteins and found that a multiplicity of folding mechanisms are possible. Fully cooperative all-or-none transition are obtained for arrays with enough sequence-similar elements and strong interactions betwee
RNA co-transcriptional folding has long been suspected to play an active role in helping proper native folding of ribozymes and structured regulatory motifs in mRNA untranslated regions. Yet, the underlying mechanisms and coding requirements for efficient co-transcriptional folding remain unclear. Traditional approaches have intrinsic limitations to dissect RNA folding paths, as they rely on sequence mutations or circular permutations that typically perturb both RNA folding paths and equilibrium structures. Here, we show that exploiting sequence symmetries instead of mutations can circumvent this problem by essentially decoupling folding paths from equilibrium structures of designed RNA sequences. Using bistable RNA switches with symmetrical helices conserved under sequence reversal, we demonstrate experimentally that native and transiently formed helices can guide efficient co-transcriptional folding into either long-lived structure of these RNA switches. Their folding path is controlled by the order of helix nucleations and subsequent exchanges during transcription, and may also be redirected by transient antisense interactions. Hence, transient intra- and intermolecular base pai
Extending Stallings' foldings of trees, we show in this article that every parallel-preserving map between median graphs factors as an isometric embedding through a sequence of elementary transformations which we call foldings and swellings. This new construction proposes a unified point of view on Beeker and Lazarovich's work on folding pocsets and on Ben-Zvi, Kropholler, and Lyman's work on folding nonpositively curved cube complexes.
A general theoretical framework is developed using free energy functional methods to understand the effects of heterogeneity in the folding of a well-designed protein. Native energetic heterogeneity arising from non-uniformity in native stability, as well as entropic heterogeneity intrinsic to the topology of the native structure are both investigated as to their impact on the folding free energy landscape and resulting folding mechanism. Given a minimally frustrated protein, both structural and energetic heterogeneity lower the thermodynamic barrier to folding, and designing in sufficient heterogeneity can eliminate the barrier at the folding transition temperature. Sequences with different distributions of stability throughout the protein and correspondingly different folding mechanisms may still be good folders to the same structure. This theoretical framework allows for a systematic study of the coupled effects of energetics and topology in protein folding, and provides interpretations and predictions for future experiments which may investigate these effects.
An interior penalty discontinuous Galerkin method is devised to approximate minimizers of a linear folding model by discontinuous isoparametric finite element functions that account for an approximation of a folding arc. The numerical analysis of the discrete model includes an a priori error estimate in case of an accurate representation of the folding curve by the isoparametric mesh. Additional estimates show that geometric consistency errors may be controlled separately if the folding arc is approximated by piecewise polynomial curves. Various numerical experiments are carried out to validate the a priori error estimate for the folding model.
The foldings of a connected graph $G$ are defined as follows. First, $G$ is a folding of itself. Let $G'$ be a graph obtained from $G$ by identifying two vertices at distance 2 in $G$. Then every folding of $G'$ is a folding of $G$. The folding number of $G$ is the minimum order of a complete folding of $G$. Theorem: The folding number of every graph equals its chromatic number.
A simple statistical mechanical model proposed by Wako and Sait$\hat{\rm o}$ has explained the aspects of protein folding surprisingly well. This model was systematically applied to multiple proteins by Mu$\tilde{\rm n}$oz and Eaton and has since been referred to as the Wako-Sait$\hat{\rm o}$-Mu$\tilde{\rm n}$oz-Eaton (WSME) model. The success of the WSME model in explaining the folding of many proteins has verified the hypothesis that the folding is dominated by native interactions, which makes the energy landscape globally biased toward native conformation. Using the WSME and other related models, Sait$\hat{\rm o}$ emphasized the importance of the hierarchical pathway in protein folding; folding starts with the creation of contiguous segments having a native-like configuration and proceeds as growth and coalescence of these segments. The $φ$-values calculated for barnase with the WSME model suggested that segments contributing to the folding nucleus are similar to the structural modules defined by the pattern of native atomic contacts. The WSME model was extended to explain folding of multi-domain proteins having a complex topology, which opened the way to comprehensively underst
Why is it difficult to refold a previously folded sheet of paper? We show that even crease patterns with only one designed folding motion inevitably contain an exponential number of `distractor' folding branches accessible from a bifurcation at the flat state. Consequently, refolding a sheet requires finding the ground state in a glassy energy landscape with an exponential number of other attractors of higher energy, much like in models of protein folding (Levinthal's paradox) and other NP-hard satisfiability (SAT) problems. As in these problems, we find that refolding a sheet requires actuation at multiple carefully chosen creases. We show that seeding successful folding in this way can be understood in terms of sub-patterns that fold when cut out (`folding islands'). Besides providing guidelines for the placement of active hinges in origami applications, our results point to fundamental limits on the programmability of energy landscapes in sheets.
How proteins fold remains a central unsolved problem in biology. While the idea of a folding code embedded in the amino acid sequence was introduced more than 6 decades ago, this code remains undefined. While we now have powerful predictive tools to predict the final native structure of proteins, we still lack a predictive framework for how sequences dictate folding pathways. Two main conceptual models dominate as explanations of folding mechanism: the funnel model, in which folding proceeds through many alternative routes on a rugged, hyperdimensional energy landscape; and the foldon model, which proposes a hierarchical sequence of discrete intermediates. Recent advances on two fronts are now enabling folding studies in unprecedented ways. Powerful experimental approaches; in particular, single-molecule force spectroscopy and hydrogen (deuterium exchange assays) allow time-resolved tracking of the folding process at high resolution. At the same time, computational breakthroughs culminating in algorithms such as AlphaFold have revolutionized static structure prediction, opening opportunities to extend machine learning toward dynamics. Together, these developments mark a turning poi
One of the most puzzling and unsolved challenges in molecular biology is understanding how proteins fold. Despite having advanced predictive tools that can accurately estimate the native structures of proteins, we still lack a comprehensive model that explains how amino acid sequences dictate folding pathways and trajectories. This manuscript takes a fresh approach to this problem by resorting to the principle of least action. This approach enables us to explore an intriguing question: how does a protein achieve its native state at a constant folding rate and within a time frame that is biologically plausible? A response to this inquiry will help us understand why proteins must fold along specific pathways and identify the boundary conditions that restrict their availability. It will also clarify why different folding pathways could be characterized by a common effective folding trajectory. Finally, it will provide a clear explanation for Levinthal's paradox. Our results are expected to pave the way for a more profound understanding of how proteins fold, shedding light on how the amino acid sequence and its surrounding environment encode the protein's folding pathways and, conseque