Correction: Finkelstein et al. How Can Ice Emerge at 0 °C? Biomolecules 2022, 12, 981.

We regret to state that our article "How Can Ice Emerge at 0 °C?" [...].

Protein structure prediction using the evolutionary algorithm USPEX.

Protein structure prediction is one of major problems of modern biophysics: current attempts to predict the tertiary protein structure from amino acid sequence are successful mostly when the use of big data and machine learning allows one to reduce the "prediction problem" to the "problem of recognition". Compared with recent successes of deep learning, classical predictive methods lag behind in their accuracy for the prediction of stable conformations. Therefore, in this work we extended the evolutionary algorithm USPEX to predict protein structure based on global optimization starting with the amino acid sequence. Moreover, we compared frequently used force fields for the task of protein structure prediction. Protein structure relaxation and energy calculations were performed using Tinker (with several different force fields) and Rosetta (with REF2015 force field) codes. To create new protein structure models in the USPEX algorithm, we developed novel variation operators. The test of the new method on seven proteins having (for simplicity) no cis-proline (with ω ≈ 0°) residues, and a length of up to 100 residues, revealed that our algorithm predicts tertiary structures of proteins with high accuracy. The comparison of the final potential energies of the predicted protein structures obtained using the USPEX and the Rosetta Abinitio approach showed that in most cases the developed algorithm found structures with close or even lower energy (Amber/Charmm/Oplsaal) and scoring function (REF2015). While USPEX has clearly demonstrated its ability to find very deep energy minima, our study showed that the existing force fields are not sufficiently accurate for accurate blind prediction of protein structures without further experimental verification.

Physics of Ice Nucleation and Antinucleation: Action of Ice-Binding Proteins.

Ice-binding proteins are crucial for the adaptation of various organisms to low temperatures. Some of these, called antifreeze proteins, are usually thought to inhibit growth and/or recrystallization of ice crystals. However, prior to these events, ice must somehow appear in the organism, either coming from outside or forming inside it through the nucleation process. Unlike most other works, our paper is focused on ice nucleation and not on the behavior of the already-nucleated ice, its growth, etc. The nucleation kinetics is studied both theoretically and experimentally. In the theoretical section, special attention is paid to surfaces that bind ice stronger than water and thus can be "ice nucleators", potent or relatively weak; but without them, ice cannot be nucleated in any way in calm water at temperatures above -30 °C. For experimental studies, we used: (i) the ice-binding protein mIBP83, which is a previously constructed mutant of a spruce budworm Choristoneura fumiferana antifreeze protein, and (ii) a hyperactive ice-binding antifreeze protein, RmAFP1, from a longhorn beetle Rhagium mordax. We have shown that RmAFP1 (but not mIBP83) definitely decreased the ice nucleation temperature of water in test tubes (where ice originates at much higher temperatures than in bulk water and thus the process is affected by some ice-nucleating surfaces) and, most importantly, that both of the studied ice-binding proteins significantly decreased the ice nucleation temperature that had been significantly raised in the presence of potent ice nucleators (CuO powder and ice-nucleating bacteria Pseudomonas syringae). Additional experiments on human cells have shown that mIBP83 is concentrated in some cell regions of the cooled cells. Thus, the ice-binding protein interacts not only with ice, but also with other sites that act or potentially may act as ice nucleators. Such ice-preventing interaction may be the crucial biological task of ice-binding proteins.

The Molten Globule State of a Globular Protein in a Cell Is More or Less Frequent Case Rather than an Exception.

Quite a long time ago, Oleg B. Ptitsyn put forward a hypothesis about the possible functional significance of the molten globule (MG) state for the functioning of proteins. MG is an intermediate between the unfolded and the native state of a protein. Its experimental detection and investigation in a cell are extremely difficult. In the last decades, intensive studies have demonstrated that the MG-like state of some globular proteins arises from either their modifications or interactions with protein partners or other cell components. This review summarizes such reports. In many cases, MG was evidenced to be functionally important. Thus, the MG state is quite common for functional cellular proteins. This supports Ptitsyn's hypothesis that some globular proteins may switch between two active states, rigid (N) and soft (MG), to work in solution or interact with partners.

How Can Ice Emerge at 0 °C?

The classical nucleation theory shows that bulk water freezing does not occur at temperatures above ≈ -30 °C, and that at higher temperatures ice nucleation requires the presence of some ice-binding surfaces. The temperature and rate of ice nucleation depend on the size and level of complementarity between the atomic structure of these surfaces and various H-bond-rich/depleted crystal planes. In our experiments, the ice nucleation temperature was within a range from -8 °C to -15 °C for buffer and water in plastic test tubes. Upon the addition of ice-initiating substances (i.e., conventional AgI or CuO investigated here), ice appeared in a range from -3 °C to -7 °C, and in the presence of the ice-nucleating bacterium Pseudomonas syringae from -1 °C to -2 °C. The addition of an antifreeze protein inhibited the action of the tested ice-initiating agents.

Calculation of Crystal-Solution Dissociation Constants.

The calculation of dissociation constants is an important problem in molecular biophysics. For such a calculation, it is important to correctly calculate both terms of the binding free energy; that is, the enthalpy and entropy of binding. Both these terms can be computed using molecular dynamics simulations, but this approach is very computationally expensive, and entropy calculations are especially slow. We develop an alternative very fast method of calculating the binding entropy and dissociation constants. The main part of our approach is based on the evaluation of movement ranges of molecules in the bound state. Then, the range of molecular movements in the bound state (here, in molecular crystals) is used for the calculation of the binding entropies and, then (using, in addition, the experimentally measured sublimation enthalpies), the crystal-to-vapor dissociation constants. Previously, we considered the process of the reversible sublimation of small organic molecules from crystals to vapor. In this work, we extend our approach by considering the dissolution of molecules, in addition to their sublimation. Similar to the sublimation case, our method shows a good correlation with experimentally measured dissociation constants at the dissolution of crystals.

Protein folding problem: enigma, paradox, solution.

The ability of protein chains to spontaneously form their three-dimensional structures is a long-standing mystery in molecular biology. The most conceptual aspect of this mystery is how the protein chain can find its native, "working" spatial structure (which, for not too big protein chains, corresponds to the global free energy minimum) in a biologically reasonable time, without exhaustive enumeration of all possible conformations, which would take billions of years. This is the so-called "Levinthal's paradox." In this review, we discuss the key ideas and discoveries leading to the current understanding of protein folding kinetics, including folding landscapes and funnels, free energy barriers at the folding/unfolding pathways, and the solution of Levinthal's paradox. A special role here is played by the "all-or-none" phase transition occurring at protein folding and unfolding and by the point of thermodynamic (and kinetic) equilibrium between the "native" and the "unfolded" phases of the protein chain (where the theory obtains the simplest form). The modern theory provides an understanding of key features of protein folding and, in good agreement with experiments, it (i) outlines the chain length-dependent range of protein folding times, (ii) predicts the observed maximal size of "foldable" proteins and domains. Besides, it predicts the maximal size of proteins and domains that fold under solely thermodynamic (rather than kinetic) control. Complementarily, a theoretical analysis of the number of possible protein folding patterns, performed at the level of formation and assembly of secondary structures, correctly outlines the upper limit of protein folding times.

Solution of Levinthal's Paradox and a Physical Theory of Protein Folding Times.

"How do proteins fold?" Researchers have been studying different aspects of this question for more than 50 years. The most conceptual aspect of the problem is how protein can find the global free energy minimum in a biologically reasonable time, without exhaustive enumeration of all possible conformations, the so-called "Levinthal's paradox." Less conceptual but still critical are aspects about factors defining folding times of particular proteins and about perspectives of machine learning for their prediction. We will discuss in this review the key ideas and discoveries leading to the current understanding of folding kinetics, including the solution of Levinthal's paradox, as well as the current state of the art in the prediction of protein folding times.

Life in Phases: Intra- and Inter- Molecular Phase Transitions in Protein Solutions.

Proteins, these evolutionarily-edited biological polymers, are able to undergo intramolecular and intermolecular phase transitions. Spontaneous intramolecular phase transitions define the folding of globular proteins, whereas binding-induced, intra- and inter- molecular phase transitions play a crucial role in the functionality of many intrinsically-disordered proteins. On the other hand, intermolecular phase transitions are the behind-the-scenes players in a diverse set of macrosystemic phenomena taking place in protein solutions, such as new phase nucleation in bulk, on the interface, and on the impurities, protein crystallization, protein aggregation, the formation of amyloid fibrils, and intermolecular liquid-liquid or liquid-gel phase transitions associated with the biogenesis of membraneless organelles in the cells. This review is dedicated to the systematic analysis of the phase behavior of protein molecules and their ensembles, and provides a description of the major physical principles governing intramolecular and intermolecular phase transitions in protein solutions.

The Reverse Side of a Coin: "Factor-Free" Ribosomal Protein Synthesis In Vitro is a Consequence of the In Vivo Proofreading Mechanism.

This paper elucidates a close connection between two well-known facts that until now have seemed independent: (i) the quality control ("proofreading") of the emerging amino acid sequence, occurring during the normal, elongation-factor-dependent ribosomal biosynthesis, which is performed by removing those Aa-tRNAs (aminoacyl tRNAs) whose anticodons are not complementary to the exhibited mRNA codons, and (ii) the in vitro discovered existence of the factor-free ribosomal synthesis of polypeptides. It is shown that a biological role of proofreading is played by a process that is exactly opposite to the step of factor-free binding of Aa-tRNA to the ribosome-exposed mRNA: a factor-free removal of that Aa-tRNA whose anticodon is not complementary to the ribosome-exhibited mRNA codon.

What is Responsible for Atypical Dependence of the Rate of Amyloid Formation on Protein Concentration: Fibril-Catalyzed Initiation of New Fibrils or Competition with Oligomers?

An abnormal dependence of the rate of amyloid formation on protein concentration has been recently observed by Meisl et al. for Aß40 peptides associated with Alzheimer's disease. To explain this effect, Meisl et al. proposed a novel mechanism of fibril growth: the fibril-catalyzed initiation of fibril formation. In this paper we offer an alternative explanation of the observed anomalous kinetics: formation of metastable oligomers competing with fibril formation by decreasing the concentration of the fibril-forming free monomers. Here we show that the oligomer sizes resulting from the anomalous dependence of the fibril growth rate on protein concentration are close to the sizes of oligomers observed by electron microscopy.

Sublimation Entropy and Dissociation Constants Prediction by Quantitative Evaluation of Molecular Mobility in Crystals.

Prediction of binding free energies (or dissociation constants) is a crucial challenge for computational biochemistry. One of the main problems here consists in fast and accurate evaluation of binding entropy, which is far more time-consuming than evaluation of binding enthalpy. Here, we offer a fast and rather accurate approach to evaluate the sublimation entropy (i.e., entropy of binding of a vapor molecule to a crystal, taken with the opposite sign) from the average range of molecular movements in the solid state. To estimate this range (and the corresponding amplitude), we considered equilibrium sublimation of small organic molecules from molecular crystals. The calculations were based on experimental data on the sublimation enthalpy, pressure of saturated vapor, and structural characteristics of the molecule in question. The resulting average amplitude (0.17 ± 0.01 Å) of molecular movements in crystals was used to predict sublimation entropies and dissociation constants for sublimation of 28 molecular crystals. The results of these predictions are in close agreement with experimental values.

There and back again: Two views on the protein folding puzzle.

The ability of protein chains to spontaneously form their spatial structures is a long-standing puzzle in molecular biology. Experimentally measured folding times of single-domain globular proteins range from microseconds to hours: the difference (10-11 orders of magnitude) is the same as that between the life span of a mosquito and the age of the universe. This review describes physical theories of rates of overcoming the free-energy barrier separating the natively folded (N) and unfolded (U) states of protein chains in both directions: "U-to-N" and "N-to-U". In the theory of protein folding rates a special role is played by the point of thermodynamic (and kinetic) equilibrium between the native and unfolded state of the chain; here, the theory obtains the simplest form. Paradoxically, a theoretical estimate of the folding time is easier to get from consideration of protein unfolding (the "N-to-U" transition) rather than folding, because it is easier to outline a good unfolding pathway of any structure than a good folding pathway that leads to the stable fold, which is yet unknown to the folding protein chain. And since the rates of direct and reverse reactions are equal at the equilibrium point (as follows from the physical "detailed balance" principle), the estimated folding time can be derived from the estimated unfolding time. Theoretical analysis of the "N-to-U" transition outlines the range of protein folding rates in a good agreement with experiment. Theoretical analysis of folding (the "U-to-N" transition), performed at the level of formation and assembly of protein secondary structures, outlines the upper limit of protein folding times (i.e., of the time of search for the most stable fold). Both theories come to essentially the same results; this is not a surprise, because they describe overcoming one and the same free-energy barrier, although the way to the top of this barrier from the side of the unfolded state is very different from the way from the side of the native state; and both theories agree with experiment. In addition, they predict the maximal size of protein domains that fold under solely thermodynamic (rather than kinetic) control and explain the observed maximal size of the "foldable" protein domains.

Reduction of the Search Space for the Folding of Proteins at the Level of Formation and Assembly of Secondary Structures: A New View on the Solution of Levinthal's Paradox.

The complete volume of the protein conformation space is, by many orders of magnitude, smaller at the level of secondary structure elements than that at the level of amino acid residues; the latter, according to Levinthal's estimate, scales approximately as 10(2 L), with L being the number of residues in the chain, whereas the former, as demonstrated in this paper, scales no faster than ∼L(N), with N being the number of the secondary structure elements, which is approximately equal to L/15. This drastic decrease in the exponent (L/15 instead of 2 L) explains why sampling of the conformation space does not contradict the ability of the protein chain to find its most stable fold.

Time to overcome the high, long, and bumpy free energy barrier in a multi-stage process: the generalized steady-state approach.

This paper presents a simple but general analytical estimate of the characteristic time of crossing a long, high, and arbitrary bumpy free energy barrier in the course of a sequential (one-dimensional) chemical, biochemical, or physical reaction. This estimate is based on a generalized steady-state approach developed in this paper.

A structural perspective of compensatory evolution.

The study of molecular evolution is important because it reveals how protein functions emerge and evolve. Recently, several types of studies indicated that substitutions in molecular evolution occur in a compensatory manner, whereby the occurrence of a substitution depends on the amino acid residues at other sites. However, a molecular or structural basis behind the compensation often remains obscure. Here, we review studies on the interface of structural biology and molecular evolution that revealed novel aspects of compensatory evolution. In many cases structural studies benefit from evolutionary data while structural data often add a functional dimension to the study of molecular evolution.