Design of stable circular permutants of the GroEL chaperone apical domain.

Enhancing protein stability holds paramount significance in biotechnology, therapeutics, and the food industry. Circular permutations offer a distinctive avenue for manipulating protein stability while keeping intra-protein interactions intact. Amidst the creation of circular permutants, determining the optimal placement of the new N- and C-termini stands as a pivotal, albeit largely unexplored, endeavor. In this study, we employed PONDR-FIT's predictions of disorder propensity to guide the design of circular permutants for the GroEL apical domain (residues 191-345). Our underlying hypothesis posited that a higher predicted disorder value would correspond to reduced stability in the circular permutants, owing to the increased likelihood of fluctuations in the novel N- and C-termini. To substantiate this hypothesis, we engineered six circular permutants, positioning glycines within the loops as locations for the new N- and C-termini. We demonstrated the validity of our hypothesis along the set of the designed circular permutants, as supported by measurements of melting temperatures by circular dichroism and differential scanning microcalorimetry. Consequently, we propose a novel computational methodology that rationalizes the design of circular permutants with projected stability. Video Abstract.

In Silico Simulations Reveal Molecular Mechanism of Uranyl Ion Toxicity towards DNA-Binding Domain of PARP-1 Protein.

The molecular toxicity of the uranyl ion (UO22+) in living cells is primarily determined by its high affinity to both native and potential metal-binding sites that commonly occur in the structure of biomolecules. Recent advances in computational and experimental research have shed light on the structural properties and functional impacts of uranyl binding to proteins, organic ligands, nucleic acids, and their complexes. In the present work, we report the results of the computational investigation of the uranyl-mediated loss of DNA-binding activity of PARP-1, a eukaryotic enzyme that participates in DNA repair, cell differentiation, and the induction of inflammation. The latest experimental studies have shown that the uranyl ion directly interacts with its DNA-binding subdomains, zinc fingers Zn1 and Zn2, and alters their tertiary structure. Here, we propose an atomistic mechanism underlying this process and compute the free energy change along the suggested pathway. Our Quantum Mechanics/Molecular Mechanics (QM/MM) simulations of the Zn2-UO22+ complex indicate that the uranyl ion replaces zinc in its native binding site. However, the resulting state is destroyed due to the spontaneous internal hydrolysis of the U-Cys162 coordination bond. Despite the enthalpy of hydrolysis being +2.8 kcal/mol, the overall reaction free energy change is -0.6 kcal/mol, which is attributed to the loss of domain's native tertiary structure originally maintained by a zinc ion. The subsequent reorganization of the binding site includes the association of the uranyl ion with the Glu190/Asp191 acidic cluster and significant perturbations in the domain's tertiary structure driven by a further decrease in the free energy by 6.8 kcal/mol. The disruption of the DNA-binding interface revealed in our study is consistent with previous experimental findings and explains the loss of PARP-like zinc fingers' affinity for nucleic acids.

In-Frame Deletion of Dystrophin Exons 8-50 Results in DMD Phenotype.

Mutations that prevent the production of proteins in the DMD gene cause Duchenne muscular dystrophy. Most frequently, these are deletions leading to reading-frame shift. The "reading-frame rule" states that deletions that preserve ORF result in a milder Becker muscular dystrophy. By removing several exons, new genome editing tools enable reading-frame restoration in DMD with the production of BMD-like dystrophins. However, not every truncated dystrophin with a significant internal loss functions properly. To determine the effectiveness of potential genome editing, each variant should be carefully studied in vitro or in vivo. In this study, we focused on the deletion of exons 8-50 as a potential reading-frame restoration option. Using the CRISPR-Cas9 tool, we created the novel mouse model DMDdel8-50, which has an in-frame deletion in the DMD gene. We compared DMDdel8-50 mice to C57Bl6/CBA background control mice and previously generated DMDdel8-34 KO mice. We discovered that the shortened protein was expressed and correctly localized on the sarcolemma. The truncated protein, on the other hand, was unable to function like a full-length dystrophin and prevent disease progression. On the basis of protein expression, histological examination, and physical assessment of the mice, we concluded that the deletion of exons 8-50 is an exception to the reading-frame rule.

Using AlphaFold to predict the impact of single mutations on protein stability and function.

AlphaFold changed the field of structural biology by achieving three-dimensional (3D) structure prediction from protein sequence at experimental quality. The astounding success even led to claims that the protein folding problem is "solved". However, protein folding problem is more than just structure prediction from sequence. Presently, it is unknown if the AlphaFold-triggered revolution could help to solve other problems related to protein folding. Here we assay the ability of AlphaFold to predict the impact of single mutations on protein stability (ΔΔG) and function. To study the question we extracted the pLDDT and metrics from AlphaFold predictions before and after single mutation in a protein and correlated the predicted change with the experimentally known ΔΔG values. Additionally, we correlated the same AlphaFold pLDDT metrics with the impact of a single mutation on structure using a large scale dataset of single mutations in GFP with the experimentally assayed levels of fluorescence. We found a very weak or no correlation between AlphaFold output metrics and change of protein stability or fluorescence. Our results imply that AlphaFold may not be immediately applied to other problems or applications in protein folding.

Best templates outperform homology models in predicting the impact of mutations on protein stability.

MOTIVATION: Prediction of protein stability change upon mutation (ΔΔG) is crucial for facilitating protein engineering and understanding of protein folding principles. Robust prediction of protein folding free energy change requires the knowledge of protein three-dimensional (3D) structure. In case, protein 3D structure is not available, one can predict the structure from protein sequence; however, the perspectives of ΔΔG predictions for predicted protein structures are unknown. The accuracy of using 3D structures of the best templates for the ΔΔG prediction is also unclear. RESULTS: To investigate these questions, we used a representative set of seven diverse and accurate publicly available tools (FoldX, Eris, Rosetta, DDGun, ACDC-NN, ThermoNet and DynaMut) for stability change prediction combined with AlphaFold or I-Tasser for protein 3D structure prediction. We found that best templates perform consistently better than (or similar to) homology models for all ΔΔG predictors. Our findings imply using the best template structure for the prediction of protein stability change upon mutation if the protein 3D structure is not available. AVAILABILITY AND IMPLEMENTATION: The data are available at https://github.com/ivankovlab/template-vs-model. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

Inhibition of hyaluronan secretion by novel coumarin compounds and chitin synthesis inhibitors.

Elevated plasma levels of hyaluronic acid (HA) is a disease marker in liver pathology and other inflammatory disorders. Inhibition of HA synthesis with coumarin 4-methylumbelliferone (4MU) has a beneficial effect in animal models of fibrosis, inflammation, cancer and metabolic syndrome. 4MU is an active compound of approved choleretic drug hymecromone with low bioavailability and a broad spectrum of action. New, more specific and efficient inhibitors of hyaluronan synthases (HAS) are required. We have tested several newly synthesized coumarin compounds and commercial chitin synthesis inhibitors to inhibit HA production in cell culture assay. Coumarin derivative compound VII (10'-methyl-6'-phenyl-3'H-spiro[piperidine-4,2'-pyrano[3,2-g]chromene]-4',8'-dione) demonstrated inhibition of HA secretion by NIH3T3 cells with the half-maximal inhibitory concentration (IC50) = 1.69 ± 0.75 µΜ superior to 4MU (IC50 = 8.68 ± 1.6 µΜ). Inhibitors of chitin synthesis, etoxazole, buprofezin, triflumuron, reduced HA deposition with IC50 of 4.21 ± 3.82 µΜ, 1.24 ± 0.87 µΜ and 1.48 ± 1.44 µΜ, respectively. Etoxazole reduced HA production and prevented collagen fibre formation in the CCl4 liver fibrosis model in mice similar to 4MU. Bioinformatics analysis revealed homology between chitin synthases and HAS enzymes, particularly in the pore-forming domain, containing the proposed site for etoxazole binding.

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.

HypercubeME: two hundred million combinatorially complete datasets from a single experiment.

MOTIVATION: Epistasis, the context-dependence of the contribution of an amino acid substitution to fitness, is common in evolution. To detect epistasis, fitness must be measured for at least four genotypes: the reference genotype, two different single mutants and a double mutant with both of the single mutations. For higher-order epistasis of the order n, fitness has to be measured for all 2n genotypes of an n-dimensional hypercube in genotype space forming a "combinatorially complete dataset". So far, only a handful of such datasets have been produced by manual curation. Concurrently, random mutagenesis experiments have produced measurements of fitness and other phenotypes in a high-throughput manner, potentially containing a number of combinatorially complete datasets. RESULTS: We present an effective recursive algorithm for finding all hypercube structures in random mutagenesis experimental data. To test the algorithm, we applied it to the data from a recent HIS3 protein dataset and found all 199,847,053 unique combinatorially complete genotype combinations of dimensionality ranging from two to twelve. The algorithm may be useful for researchers looking for higher-order epistasis in their high-throughput experimental data. AVAILABILITY: https://github.com/ivankovlab/HypercubeME.git. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.